We present the results from a large near-infrared spectroscopic survey with Subaru/FMOS (\textit{FastSound}) consisting of $\sim$ 4,000 galaxies at $z\sim1.4$ with significant H$\alpha$ detection. We measure the gas-phase metallicity from the [N~{\sc ii}]$\lambda$6583/H$\alpha$ emission line ratio of the composite spectra in various stellar mass and star-formation rate bins. The resulting mass-metallicity relation generally agrees with previous studies obtained in a similar redshift range to that of our sample. No clear dependence of the mass-metallicity relation with star-formation rate is found. Our result at $z\sim1.4$ is roughly in agreement with the fundamental metallicity relation at $z\sim0.1$ with fiber aperture corrected star-formation rate. We detect significant [S~{\sc ii}]$\lambda\lambda$6716,6731 emission lines from the composite spectra. The electron density estimated from the [S~{\sc ii}]$\lambda\lambda$6716,6731 line ratio ranges from 10 -- 500 cm$^{-3}$, which generally agrees with that of local galaxies. On the other hand, the distribution of our sample on [N~{\sc ii}]$\lambda$6583/H$\alpha$ vs. [S~{\sc ii}]$\lambda\lambda$6716,6731/H$\alpha$ is different from that found locally. We estimate the nitrogen-to-oxygen abundance ratio (N/O) from the N2S2 index, and find that the N/O in galaxies at $z\sim1.4$ is significantly higher than the local values at a fixed metallicity and stellar mass. The metallicity at $z\sim1.4$ recalculated with this N/O enhancement taken into account decreases by 0.1 -- 0.2 dex. The resulting metallicity is lower than the local fundamental metallicity relation.
We analyze three Chandra observations, with a combined exposure time of 99 ks, of the Galactic supernova remnant RCW 103, a young supernova remnant, previously with no clear detection of metal-rich ejecta. Based on our imaging and spectral analyses of these deep Chandra data, we find evidence for metal-rich ejecta emission scattered throughout the remnant. X-ray emission from the shocked ejecta is generally weak, and the shocked circumstellar medium (CSM) is a largely dominant component across the entire remnant. The CSM component shows abundances of ~0.5 solar, while Ne, Mg, Si, S, and Fe abundances of the ejecta are up to a few times solar. Comparison of these ejecta abundances with yields from supernova nucleosynthesis models suggests, together with the existence of a central neutron star, a progenitor mass of ~18-20 M$_\odot$, though the Fe/Si ratios are larger than predicted. The shocked CSM emission suggests a progenitor with high mass-loss rate and subsolar metallicity.
We examine 11 XMM-Newton observations of the giant spiral galaxy NGC 1961, with a total integration time of 289 ks ($\sim 100$ ks after flaring corrections). These deep X-ray data allow us to study the hot gaseous halo of a spiral galaxy in unprecedented detail. We perform both a spatial and a spectral analysis; with the former, the hot halo is detected to at least 80 kpc and with the latter the halo properties can be measured in detail up to 42 kpc. In the region of overlap, there is good agreement between the two methods. We measure the temperature profile of the hot halo, finding a negative gradient as is common for elliptical galaxies. We also measure a rough metallicity profile, which is consistent with being flat at a sub-Solar value ($Z \sim 0.2 Z_{\odot}$). Converting to this metallicity, our deprojected density profile is consistent with previous parametric fits, with no evidence for a break or flattening within the inner 42 kpc (about 10% of the virial radius). We infer pressure and entropy profiles for the hot halo, and use the former to estimate the mass profile of the galaxy assuming hydrostatic equilibrium. Extrapolating these profiles to the virial radius, we infer a hot gaseous halo mass comparable to the stellar mass of the galaxy, and a total baryon fraction from the stars and hot gas of around 30%. We show that the cooling time of the hot gas is orders of magnitude longer than the dynamical time, making the hot halo stable against cooling instabilities, and argue that an extended stream of neutral Hydrogen seen to the NW of this galaxy is likely due to accretion from the intergalactic medium. The low metallicity of the hot halo suggests it too was likely accreted. We compare the hot halo of NGC 1961 to hot halos around isolated elliptical galaxies, and show that the total mass better determines the hot halo properties than the stellar mass.
Supermassive black hole (BH) mergers produce powerful gravitational wave (GW) emission. Asymmetry in this emission imparts a recoil kick to the merged BH, which can eject the BH from its host galaxy altogether. Recoiling BHs could be observed as offset active galactic nuclei (AGN). Several candidates have been identified, but systematic searches have been hampered by large uncertainties regarding their observability. By extracting merging BHs and host galaxy properties from the Illustris cosmological simulations, we have developed a comprehensive model for recoiling AGN. Here, for the first time, we model the effects of BH spin alignment and recoil dynamics based on the gas-richness of host galaxies. For comparable assumptions, we find much higher rates of recoiling AGN than Volonteri & Madau (2008), indicating systematic differences between BH populations in semi-analytic models and cosmological simulations. We predict that if BH spins are not highly aligned, seeing-limited observations could resolve offset AGN, making them promising targets for all-sky surveys. For randomly-oriented spins, less than about 10 spatially-offset AGN may be detectable in HST-COSMOS, and > 10^3 could be found with Pan-STARRS, LSST, Euclid, and WFIRST. Nearly a thousand velocity-offset AGN are predicted within the SDSS footprint; the rarity of large broad-line offsets among SDSS quasars is likely due in part to selection effects but suggests that spin alignment plays a role in suppressing recoils. Nonetheless, in our most physically motivated model where alignment occurs only in gas-rich mergers, hundreds of offset AGN should be found in all-sky surveys. Our findings strongly motivate a dedicated search for recoiling AGN.
Observations at low redshifts thus far fail to account for all of the baryons expected in the Universe according to cosmological constraints. A large fraction of the baryons presumably resides in a thin and warm-hot medium between the galaxies, where they are difficult to observe due to their low densities and high temperatures. Cosmological simulations of structure formation can be used to verify this picture and provide quantitative predictions for the distribution of mass in different large-scale structure components. Here we study the distribution of baryons and dark matter at different epochs using data from the Illustris Simulation. We identify regions of different dark matter density with the primary constituents of large-scale structure, allowing us to measure mass and volume of haloes, filaments and voids. At redshift zero, we find that 49 % of the dark matter and 23 % of the baryons are within haloes. The filaments of the cosmic web host a further 45 % of the dark matter and 46 % of the baryons. The feedback models used in Illustris have a significant impact on the baryon distribution at large-scales, leading to 31 % of the baryons residing in dark matter voids. Categorizing the baryons according to their density and temperature, we find that 17.8 % of them are in a condensed state, 21.6 % are present as cold, diffuse gas, and 53.9 % are found in the state of a warm-hot intergalactic medium.
We show that the speed of the fastest coronal mass ejections (CMEs) that an active region (AR) can produce can be predicted from a vector magnetogram of the AR. This is shown by logarithmic plots of CME speed (from the SOHO LASCO CME catalog) versus each of ten AR-integrated magnetic parameters (AR magnetic flux, three different AR magnetic-twist parameters, and six AR free-magnetic-energy proxies) measured from the vertical and horizontal field components of vector magnetograms (from the {\it Solar Dynamics Observatory's Helioseismic and Magnetic Imager}) of the source ARs of 189 CMEs. These plots show: (1) the speed of the fastest CMEs that an AR can produce increases with each of these whole-AR magnetic parameters, and (2) that one of the AR magnetic-twist parameters and the corresponding free-magnetic-energy proxy each determine the CME-speed upper-limit line somewhat better than any of the other eight whole-AR magnetic parameters.
Aims: Understanding the fragmentation and collapse properties of the dense gas during the onset of high-mass star formation. Methods: We observed the massive (~800M_sun) starless gas clump IRDC18310-4 with the Plateau de Bure Interferometer (PdBI) at sub-arcsecond resolution in the 1.07mm continuum andN2H+(3-2) line emission. Results: Zooming from a single-dish low-resolution map to previous 3mm PdBI data, and now the new 1.07mm continuum observations, the sub-structures hierarchically fragment on the increasingly smaller spatial scales. While the fragment separations may still be roughly consistent with pure thermal Jeans fragmentation, the derived core masses are almost two orders of magnitude larger than the typical Jeans mass at the given densities and temperatures. However, the data can be reconciled with models using non-homogeneous initial density structures, turbulence and/or magnetic fields. While most sub-cores remain (far-)infrared dark even at 70mum, we identify weak 70mum emission toward one core with a comparably low luminosity of ~16L_sun, re-enforcing the general youth of the region. The spectral line data always exhibit multiple spectral components toward each core with comparably small line widths for the individual components (in the 0.3 to 1.0km/s regime). Based on single-dish C18O(2-1) data we estimate a low virial-to-gas-mass ratio <=0.25. We discuss that the likely origin of these spectral properties may be the global collapse of the original gas clump that results in multiple spectral components along each line of sight. Even within this dynamic picture the individual collapsing gas cores appear to have very low levels of internal turbulence.
We present deep Hubble Space Telescope Advanced Camera for Surveys observations of the stellar populations in two fields lying at 20 and 23 kpc from the centre of M31 along the south-west semi-major axis. These data enable the construction of colour-magnitude diagrams reaching the oldest main-sequence turn-offs (~13 Gyr) which, when combined with another field at 25 kpc from our previous work, we use to derive the first precision constraints on the spatially-resolved star formation history of the M31 disc. The star formation rates exhibit temporal as well as field-to-field variations, but are generally always within a factor of two of their time average. There is no evidence of inside-out growth over the radial range probed. We find a median age of ~7.5 Gyr, indicating that roughly half of the stellar mass in the M31 outer disc was formed before z ~ 1. We also find that the age-metallicity relations (AMRs) are smoothly increasing from [Fe/H]~-0.4 to solar metallicity between 10 and 3 Gyr ago, contrary to the flat AMR of the Milky Way disc at a similar number of scale lengths. Our findings provide insight on the roles of stellar feedback and radial migration in the formation and evolution of large disc galaxies.
This article presents a calculation of the mean electromotive force arising from general small-scale magnetohydrodynamical turbulence, within the framework of the second-order correlation approximation. With the goal of improving understanding of the accretion disk dynamo, effects arising through small-scale magnetic fluctuations, velocity gradients, density and turbulence stratification, and rotation, are included. The primary result, which supplements numerical findings, is that an off-diagonal turbulent resistivity due to magnetic fluctuations can produce large-scale dynamo action -- the magnetic analogue of the "shear-current" effect. In addition, consideration of $\alpha$ effects in the stratified regions of disks gives the puzzling result that there is no strong prediction for a sign of $\alpha$, since the effects due to kinetic and magnetic fluctuations, as well as those due to shear and rotation, are each of opposing signs and tend to cancel each other.
We carry out N-body simulations to examine the effects of dynamical interactions on planetary systems in young open star clusters. We explore how the planetary populations in these star clusters evolve, and how this evolution depends on the initial amount of substructure, the virial ratio, the cluster mass and density, and the initial semi-major axis of the planetary systems. The fraction of planetary systems that remains intact as a cluster member, fbps, is generally well-described by the functional form fbps=f0(1+[a/a0]^c)^-1, where (1-f0) is the fraction of stars that escapes from the cluster, a0 the critical semi-major axis for survival, and c a measure for the width of the transition region. The effect of the initial amount of substructure over time can be quantified as fbps=A(t)+B(D), where A(t) decreases nearly linearly with time, and B(D) decreases when the clusters are initially more substructured. Provided that the orbital separation of planetary systems is smaller than the critical value a0, those in clusters with a higher initial stellar density (but identical mass) have a larger probability of escaping the cluster intact. These results help us to obtain a better understanding of the difference between the observed fractions of exoplanets-hosting stars in star clusters and in the Galactic field. It also allows us to make predictions about the free-floating planet population over time in different stellar environments.
We describe digital tracking, a method for asteroid searches that greatly increases the sensitivity of a telescope to faint unknown asteroids. It has been previously used to detect faint Kuiper Belt objects using the Hubble Space Telescope and large ground-based instruments, and to find a small, fast-moving asteroid during a close approach to Earth. We complement this earlier work by developing digital tracking methodology for detecting asteroids using large-format CCD imagers. We demonstrate that the technique enables the ground-based detection of large numbers of new faint asteroids. Our methodology resolves or circumvents all major obstacles to the large-scale application of digital tracking for finding main belt and near-Earth asteroids. We find that for both asteroid populations, digital tracking can deliver a factor of ten improvement over conventional searches. Digital tracking has long been standard practice for deep Kuiper Belt surveys, but even there our methodology enables deeper integrations than have yet been attempted. Our search for main belt asteroids using a one-degree imager on the 0.9m WIYN telescope on Kitt Peak validates our methodology, delivers sensitivity to asteroids in a regime previously probed only with 4-meter and larger instruments, and leads to the detection of 156 previously unknown asteroids and 59 known objects in a single field. Digital tracking has the potential to revolutionize searches for faint moving objects ranging from the Kuiper Belt through main belt and near-Earth asteroids, and perhaps even anthropogenic space debris in low Earth orbit.
The space weather is extremely mild during solar cycle 24: the number of major geomagnetic storms and high-energy solar energetic particle events are at the lowest since the dawn of the space age. Solar wind measurements at 1 AU using Wind and ACE instruments have shown that there is a significant drop in the density, magnetic field, total pressure, and Alfven speed in the inner heliosphere as a result of the low solar activity. The drop in large space weather events is disproportionately high because the number of energetic coronal mass ejections that cause these events has not decreased significantly. For example, the rate of halo CMEs, which is a good indicator of energetic CMEs, is similar to that in cycle 23, even though the sunspot number has declined by about 40%. The mild space weather seems to be a consequence of the anomalous expansion of CMEs due to the low ambient pressure in the heliosphere. The anomalous expansion results in the dilution of the magnetic contents of CMEs, so the geomagnetic storms are generally weak. CME driven shocks propagating through the weak heliospheric field are less efficient in accelerating energetic particles, so the particles do not attain high energies. Finally, we would like to point out that extreme events such as the 2012 July 23 CMEs that occurred on the backside of the Sun and did not affect Earth except for a small proton event.
We investigate a relativistic fluid jet driven by radiation from a shocked accretion disc around a non-rotating black hole approximated by Paczy\'nski-Wiita potential. The sub-Keplerian and Keplerian accretion rates control the shock location and therefore, the radiation field around the accretion disc. We compute the radiative moments with full special relativistic transformation. The effect of a fraction of radiation absorbed by the black hole has been approximated, over and above the special relativistic transformations. We show that the radiative moments around a super massive black hole are different compared to that around a stellar mass black hole. We show that the terminal speed of jets increases with the mass accretion rates,synchrotron emission of the accretion disc and reduction of proton fraction of the flow composition. To obtain relativistic terminal velocities of jets, both thermal and radiative driving are important. We show for very high accretion rates and pair dominated flow, jets around super massive black holes are truly ultra-relativistic, while for jets around stellar mass black holes, terminal Lorentz factor of about $10$ is achievable.
The relation between the Star Formation Rate (SFR) and stellar mass (${\rm M}_{\star}$) of galaxies represents a fundamental constraint on galaxy formation. However, the observed amplitude of the star formation rate - stellar mass relation has not been successfully reproduced in simulations, indicating either that the halo accretion history and baryonic physics are poorly understood or that observations contain biases. In this paper, we examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z\sim 1-4 $ galaxies and display the inconsistency between observed relations that are obtained using different techniques. We employ cosmological hydrodynamic simulations from various groups and compare these with a range of observations. The comparison suggests that using Spectral Energy Distributions (SEDs) to estimate star formation rates, dust corrections and stellar masses produces the most reliable SFR$-{\rm M}_{\star}$ relations. On the contrary, the combination of IR and UV luminosities (UV+IR) overpredicts the SFR and dust corrections at a fixed stellar mass almost by a factor of 5 for $z \sim 1.5-4$. For $z < 1.5$, the SED fitting technique and IR+UV conversion agree well. We find remarkable agreement between the numerical results from various authors who have employed different cosmological codes and run simulations with different resolutions. This is interesting for two reasons. A) simulations can produce realistic populations of galaxies within representative cosmological volumes even at relatively modest resolutions. B) It is likely that current numerical codes that rely on similar subgrid multiphase ISM models and are tuned to reproduce statistical properties of galaxies, produce similar results for the SFR$-{\rm M}_{\star}$ relation by construction, regardless of resolution, box size and, to some extent, the adopted feedback prescriptions.
As both simulations and observations reach the resolution of the star-forming molecular clouds, it becomes important to clarify if these two techniques are discussing the same objects in galaxies. We compare clouds formed in a high resolution galaxy simulation identified as continuous structures within a contour, in the simulator's position-position-position (PPP) co-ordinate space and the observer's position-position-velocity space (PPV). Results indicate that the properties of the cloud populations are similar in both methods and up to 70% of clouds have a single counterpart in the opposite data structure. Comparing individual clouds in a one-to-one match reveals a scatter in properties mostly within a factor of two. However, the small variations in mass, radius and velocity dispersion produce significant differences in derived quantities such as the virial parameter. This makes it difficult to determine if a structure is truely gravitationally bound. The three cloud types originally found in the simulation in Fujimoto et al. (2014) are identified in both data sets, with around 80% of the clouds retaining their type between identification methods. We also compared our results when using a peak decomposition method to identify clouds in both PPP and PPV space. The number of clouds increased with this technique, but the overall cloud properties remained similar. However, the more crowded environment lowered the ability to match clouds between techniques to 40%. The three cloud types also became harder to separate, especially in the PPV data set. The method used for cloud identification therefore plays a critical role in determining cloud properties, but both PPP and PPV can potentially identify the same structures.
Aims: IKT 16 is an X-ray and radio-faint supernova remnant (SNR) in the Small
Magellanic Cloud (SMC). A detailed X-ray study of this SNR with XMM-Newton
confirmed the presence of a hard X-ray source near its centre, indicating the
detection of the first composite SNR in the SMC. With a dedicated Chandra
observation we aim to resolve the point source and confirm its nature. We also
acquire new ATCA observations of the source at 2.1 GHz with improved flux
density estimates and resolution.
Methods: We perform detailed spatial and spectral analysis of the source.
With the highest resolution X-ray and radio image of the centre of the SNR
available today, we resolve the source and confirm its pulsar wind nebula (PWN)
nature. Further, we constrain the geometrical parameters of the PWN and perform
spectral analysis for the point source and the PWN separately. We also test for
the radial variations of the PWN spectrum and its possible east west asymmetry.
Results: The X-ray source at the centre of IKT 16 can be resolved into a
symmetrical elongated feature centering a point source, the putative pulsar.
Spatial modeling indicates an extent of 5.2 arcsec of the feature with its axis
inclined at 82 degree east from north, aligned with a larger radio feature
consisting of two lobes almost symmetrical about the X-ray source. The picture
is consistent with a PWN which has not yet collided with the reverse shock. The
point source is about three times brighter than the PWN and has a hard spectrum
of spectral index 1.1 compared to a value 2.2 for the PWN. This points to the
presence of a pulsar dominated by non-thermal emission. The expected E_{dot} is
~ 10^37 erg s^-1 and spin period < 100 ms. However, the presence of a compact
nebula unresolved by Chandra at the distance of the SMC cannot completely be
ruled out.
We use a spectroscopic sample of 286 star-forming galaxies (SFGs) at 1<z<3 from the GMASS survey to study different star formation rate (SFR) estimators. Infrared (IR) data are used to derive empirical calibrations to correct ultraviolet (UV) and [OII]{\lambda}3727 luminosities for dust extinction and dust-corrected estimates of SFR. In the selection procedure we fully exploit the available spectroscopic information. On the basis of three continuum indices, we are able to identify and exclude from the sample galaxies in which old stellar populations might bring a non-negligible contribution to IR luminosity (LIR) and continuum reddening. Using Spitzer-MIPS and Herschel-PACS data we derive LIR for two-thirds of our sample. The LIR/LUV ratio is used as a probe of effective attenuation (AIRX) to search for correlations with continuum and spectroscopic features. The relation between AIRX and UV continuum slope ({\beta}) was tested for our sample and found to be broadly consistent with the literature results at the same redshift, though with a larger dispersion with respect to UV-selected samples. We find a correlation between the rest-frame equivalent width (EW) of the [OII]{\lambda}3727 line and {\beta}, which is the main result of this work. We therefore propose the [OII]{\lambda}3727 line EW as a dust attenuation probe and calibrate it through AIRX, though the assumption of a reddening curve is still needed to derive the actual attenuation towards the [OII]{\lambda}3727 line. We tested the issue of differential attenuation towards stellar continuum and nebular emission: our results are in line with the traditional prescription of extra attenuation towards nebular lines. A set of relations is provided that allows the recovery of the total unattenuated SFR from UV and [OII]{\lambda}3727 luminosities. (Abridged)
We present an estimate of the third integral of motion for axisymmetric
three-dimensional potentials. This estimate is based on a Staeckel
approximation and is explicitly written as a function of the potential. We
tested this scheme for the Besancon Galactic model and two other disc-halo
models and find that orbits of disc stars have an accurately conserved third
quasi integral.
The accuracy ranges from of 0.1% to 1% for heights varying from z = 0~kpc to
z= 6 kpc and Galactocentric radii R from 5 to 15kpc.
We also tested the usefulness of this quasi integral in analytic distribution
functions of disc stellar populations: we show that the distribution function
remains approximately stationary and that it allows to recover the potential
and forces by applying Jeans equations to its moments.
We improve previous calculations of the CMB spectral distortions due to the decay of primordial magnetic fields. We focus our studies on causally generated magnetic fields at the electroweak and QCD phase transitions. We also consider the decay of helical magnetic fields. We show that the decay of non-helical magnetic fields generated at either the electroweak or QCD scale produce $\mu$ and $y$-type distortions below $10^{-8}$ which are probably not detectable by a future PIXIE-like experiment. We show that magnetic fields generated at the electroweak scale must have a helicity fraction $f_*>10^{-4}$ in order to produce detectable $\mu$-type distortions. Hence a positive detection coming from the decay of magnetic fields would rule out non-helical primordial magnetic fields and provide a lower bound on the magnetic helicity.
Six asteroids including two NEAs, one of which is PHA, accessible for observation in September 2012 were investigated using a low-resolution spectrophotometry in the range 0.35-0.90 um with the aim to study features of their reflectance spectra. For the first time we discovered likely spectral signs (as a maximum at 0.4-0.6 um in reflectance spectra) of simultaneous sublimation activity and presence of a temporal coma on three primitive-type main-belt asteroids, Adeona, Interamnia, and Nina, being at perihelion distances or approaching to it. We suggest that such a cometary-like activity may be a common phenomenon at the highest subsolar surface temperatures for C and close type asteroids including considerable amounts of H2O and CO2 ices beneath the surface. However, excavation of fresh ice at recent impact event(s) could be an alternative explanation of the phenomenon. Similar absorption bands centered at 0.38, 0.44 and 0.67-0.71 um registered in the reflectance spectra of Adeona, Interamnia, and Nina clearly point to predominantly silicate surface matter. To specify its content, we performed laboratory investigations of ground samples of known carbonaceous chondrites (Orguel, Mighei, Murchison, and Boriskino) and seven samples of low-iron Mg serpentines as possible analogs of the asteroids. In particular, we found that the equivalent width of the band centered at 0.44 um in reflectance spectra of the low-Fe serpentine samples has a high correlation with content of Fe3+ (octahedral and tetrahedral). It means that the absorption feature can be used as an indicator of ferric iron in oxidized and hydrated low-Fe silicate compounds on asteroids (abridged).
Infrared Dark Clouds (IRDCs) harbor the earliest phases of massive star
formation, and many of the compact cores in IRDCs, traced by millimeter
continuum or by molecular emission in high critical density lines, host massive
young stellar objects (YSOs). We used the Robert C. Byrd Green Bank Telescope
(GBT) and the Karl G. Jansky Very Large Array (VLA) to map NH${}_{3}$ and CCS
in nine IRDCs to reveal the temperature, density, and velocity structures and
explore chemical evolution in the dense ($>10^{22}$ cm${}^{-2}$) gas. Ammonia
is an excellent molecular tracer for these cold, dense environments. The
internal structure and kinematics of the IRDCs include velocity gradients,
filaments, and possibly colliding clumps that elucidate the formation process
of these structures and their YSOs. We find a wide variety of substructure
including filaments and globules at distinct velocities, sometimes overlapping
at sites of ongoing star formation. It appears that these IRDCs are still being
assembled from molecular gas clumps even as star formation has already begun,
and at least three of them appear consistent with the morphology of
"hub-filament structures" discussed in the literature. Furthermore, we find
that these clumps are typically near equipartition between gravitational and
kinetic energies, so these structures may survive for multiple free-fall times.
Keywords: molecular data -- ISM: clouds -- (ISM:) dust, extinction -- ISM:
molecules -- Stars: formation -- radio lines: ISM
This paper is the second in a series describing the southern Galactic Disk Survey (GDS) performed at the Universit\"atssternwarte Bochum near Cerro Armazones in Chile. Haas et al. (2012, Paper I) presented the survey design and a the characteristics of the observations and data. They identified ~2200 variable stars in an area of 50 square degrees with more than 50 observations in 2011. Here we present the first complete version of the GDS covering all 268 fields with 1323 square degrees along the Galactic disk including revised data from Paper I. The individual fields were observed up to 272 times and comprise a maximum time span between September 2010 and May 2015. We detect a total of 64151 variable sources, which are presented in a catalog including some of their properties and their light curves. A comparison with the International Variable Star Index (VSX) and All Sky Automated Survey (ASAS) indicates that 56794 of these sources are previously unknown variables. Furthermore, we present UBVr'i'z' photometry for all sources within the GDS, resulting in a new multi-color catalog of nearly 16x10^6 sources detected in at least one filter. Both the GDS and the near-infrared VISTA Variables in the Via Lactea survey (VVV) complement each other in the overlap area of about 300 square degrees enabling future comparison studies.
New B and V band monitoring in 2014/2015 reveals that the Seyfert-1 galaxy 3C120 has brightened by 1.4 magnitudes compared to our campaign in 2009/2010. This allows us to check for the debated luminosity and time dependent color variations, claimed for SDSS quasars by Sun et al. (2014) based on an analysis in magnitude units. For our 3C120 data we find that the B/V flux ratio of the variable component in the bright epoch is indistinguishable from the faint one. We do not find any color variability on different timescales ranging from about 1 to 1800 days. We suggest that the luminosity and time dependent color variability by Sun et al. is an artifact caused by analyzing the data in magnitudes instead of fluxes. The flux variation gradients of both epochs yield consistent estimates of the host galaxy contribution to our 7.5" aperture. These results corroborate that the optical flux variation gradient method works well for Seyfert Galaxies.
Quasi-periodic oscillations (QPOs) discovered in soft-gamma repeaters (SGRs) are expected to help us to study the properties of matter in neutron stars. In earlier investigations, we identified the QPOs of frequencies below $\sim100$ Hz observed in giant flares of SGR 1806$-$20 and SGR 1900+14 as the crustal torsional oscillations. For this purpose, we calculated the frequencies of the fundamental torsional oscillations with various angular indices $\ell$, by changing the stellar mass and radius. In this work, we try to explain the additional QPO frequencies recently reported by Huppenkothen et al. (2014a,b) within the same framework as before except that we newly take into account the effect of electron screening, which acts to decrease the frequencies by a small amount. Those QPOs were discovered in two different SGRs, i.e., SGR 1806$-$20 and SGR J1550$-$5418. Then, we find that the newly observed QPO frequency in SGR 1806$-$20 can be still identified as one of the frequencies of the fundamental torsional oscillations, while those in SGR J1550$-$5418 can also be explained in terms of the torsional oscillations although the relevant angular indices are difficult to identify.
We present near-IR K-band intermediate-dispersion spatially-resolved spectroscopic observations of a limited sample of bipolar planetary nebulae (PNe). The spectra have been used to determine the excitation mechanism of the H2 molecule using standard line ratios diagnostics. The H2 molecule is predominantly shock-excited in bipolar PNe with broad equatorial rings, whereas bipolar PNe with narrow equatorial waists present either UV excitation at their cores (e.g., Hb 12) or shock-excitation at their bipolar lobes (e.g., M1-92). The shock-excitation among bipolar PNe with ring is found to be correlated with emission in the H2 1-0 S(1) line brighter than Br{\gamma}. We have extended this investigation to other PNe with available near-IR spectroscopic observations. This confirms that bipolar PNe with equatorial rings are in average brighter in H2 than in Br{\gamma} and show dominant shock excitation.
We combine optical, near-infrared and mid-infrared spectra and photometry to construct expanded spectral energy distributions (SEDs) for 145 field age (\textgreater 500 Myr) and 53 young (lower age estimate \textless 500 Myr) ultracool dwarfs (M6-T9). This range of spectral types includes very low mass stars, brown dwarfs, and planetary mass objects, providing fundamental parameters across both the hydrogen and deuterium burning minimum masses for the largest sample assembled to date. A subsample of 29 objects have well constrained ages as probable members of a nearby young moving group (NYMG). We use 182 parallaxes and 16 kinematic distances to determine precise bolometric luminosities ($L_\text{bol}$) and radius estimates from evolutionary models give semi-empirical effective temperatures ($T_\text{eff}$) for the full range of young and field age late-M, L and T dwarfs. We construct age-sensitive relationships of luminosity, temperature and absolute magnitude as functions of spectral type and absolute magnitude to disentangle the effects of degenerate physical parameters such as $T_\text{eff}$, surface gravity, and clouds on spectral morphology. We report bolometric corrections in $J$ for both field age and young objects and find differences of up to a magnitude for late-L dwarfs. Our correction in $Ks$ shows a larger dispersion but not necessarily a different relationship for young and field age sequences. We also characterize the NIR-MIR reddening of low gravity L dwarfs and identify a systematically cooler $T_\text{eff}$ of up to 300K from field age objects of the same spectral type and 400K cooler from field age objects of the same $M_H$ magnitude.
We consider Horava gravity within the framework of the EFT of dark energy and modified gravity. We work out a complete mapping of the theory into the EFT language for an action including all the operators which are relevant for linear perturbations with up to sixth order spatial derivatives. We then employ an updated version of the EFTCAMB/EFTCosmoMC package to study the cosmology of the low-energy limit of Horava gravity and place constraints on its parameters using several cosmological data sets. In particular we consider two cases: the first in which the three parameters of the low-energy theory are all varied and a second case that is tuned to evade PPN constraints, reducing the number of free parameters to two. We employ data sets which include the CMB TT and lensing power spectra by Planck 2013, WMAP low-l polarization spectra, the WiggleZ galaxy power spectrum, the local Hubble measurements, Supernovae data from SNLS, SDSS and HST and the BAO measurements from BOSS, SDSS and 6dFGS. For both cases we estimate the deviation of the cosmological gravitational constant from the local Newtonian one, getting improved upper bounds with respect to the previous ones from BBN data. At the level of the background, we find a relevant rescaling of the Hubble rate at all epoch, which has a strong impact on the cosmological observables; at the level of perturbations, we discuss all the relevant effects that the modifications of gravity induce, ranging from modifications of the late time ISW effect, the growth of matter perturbations, gravitational lensing and differences in the B-modes of the CMB. In general the quasi-static approximation is not safe to describe the evolution of perturbations in Horava gravity. Overall we find that the effects of the modifications induced by the low-energy Horava gravity action are quite dramatic and current data place tight bounds on the theory parameters.
We consider light propagation in an inhomogeneous irrotational dust universe with vanishing cosmological constant, with initial conditions as in standard linear perturbation theory. A non-perturbative approach to the dynamics of such a universe is combined with a distance formula based on the Sachs optical equations. Then a numerical study implies a redshift-distance relation that roughly agrees with observations. Interpreted in the standard homogeneous setup, our results would appear to imply the currently accepted values for the Hubble rate and the deceleration parameter; furthermore there is consistency with density perturbations at last scattering. The determination of these three quantities relies only on a single parameter related to a cutoff scale. Discrepancies with the existing literature are mainly due to effects beyond second order in perturbation theory.
A recent paper argued that it is not possible to infer the energy scale of inflation from the amplitude of tensor fluctuations in the Cosmic Microwave Background, because the usual connection is substantially altered if there are a large number of universally coupled fields present during inflation, with mass less than the inflationary Hubble scale. We give a simple argument demonstrating that this is incorrect.
We compute the one loop photon contribution to the graviton self-energy on de Sitter background and use it to solve the linearized Einstein equation for a point mass. Our results show that a co-moving observer sees a logarithmic spatial running Newton's constant. Equivalently a static observer reports a secular suppression of the Newtonian potential.
Links to: arXiv, form interface, find, astro-ph, recent, 1508, contact, help (Access key information)
We report the discovery of Antlia B, a faint dwarf galaxy at a projected distance of $\sim$72 kpc from NGC 3109 ($M_{V}$$\sim$$-$15 mag), the primary galaxy of the NGC 3109 dwarf association at the edge of the Local Group. The tip of the red giant branch distance to Antlia B is $D$=1.29$\pm$0.10 Mpc, which is consistent with the distance to NGC 3109. A qualitative analysis indicates the new dwarf's stellar population has both an old, metal poor red giant branch ($\gtrsim$10 Gyr, [Fe/H]$\sim$$-$2), and a younger blue population with an age of $\sim$200-400 Myr, analogous to the original Antlia dwarf, another likely satellite of NGC 3109. Antlia B has \ion{H}{1} gas at a velocity of $v_{helio,HI}$=376 km s$^{-1}$, confirming the association with NGC 3109 ($v_{helio}$=403 km s$^{-1}$). The HI gas mass (M$_{HI}$=2.8$\pm$0.2$\times$10$^{5}$ M$_{\odot}$), stellar luminosity ($M_{V}$=$-$9.7$\pm$0.6 mag) and half light radius ($r_{h}$=273$\pm$29 pc) are all consistent with the properties of dwarf irregular and dwarf spheroidal galaxies in the Local Volume, and is most similar to the Leo P dwarf galaxy. The discovery of Antlia B is the initial result from a Dark Energy Camera survey for halo substructure and faint dwarf companions to NGC 3109 with the goal of comparing observed substructure with expectations from the $\Lambda$+Cold Dark Matter model in the sub-Milky Way regime.
We present the nonlinear 2D galaxy power spectrum, $P(k,\mu)$, in redshift space, measured from the Dark Sky simulations, using galaxy catalogs constructed with both halo occupation distribution and subhalo abundance matching methods, chosen to represent an intermediate redshift sample of luminous red galaxies. We find that the information content in individual $\mu$ (cosine of the angle to the line of sight) bins is substantially richer then multipole moments, and show that this can be used to isolate the impact of nonlinear growth and redshift space distortion (RSD) effects. Using the $\mu<0.2$ simulation data, which we show is not impacted by RSD effects, we can successfully measure the nonlinear bias to an accuracy of $\sim 5$% at $k<0.6 h$Mpc$^{-1}$. This use of individual $\mu $ bins to extract the nonlinear bias successfully removes a large parameter degeneracy when constraining the linear growth rate of structure. We carry out a joint parameter estimation, using the low $\mu$ simulation data to constrain the nonlinear bias, and $\mu\ge0.2$ to constrain the growth rate and show that $f$ can be constrained to $\sim 26\, (22)$% to a $k_{\rm max}< 0.4\, (0.6) h$Mpc$^{-1}$ from clustering alone using a simple dispersion model, for a range of galaxy models. Our analysis of individual $\mu $ bins also reveals interesting physical effects which arise simply from different methods of populating halos with galaxies. We find a prominent turnaround scale, at which RSD damping effects are greater then the nonlinear growth, which differs not only for each $\mu$ bin but also for each galaxy model. These features may provide unique signatures which could be used to shed light on the galaxy-dark matter connection.
Numerous observational studies have revealed the ubiquitous presence of multiple stellar populations in globular clusters and cast many hard challenges for the study of the formation and dynamical history of these stellar systems. In this Letter we present the results of a study of the kinematic properties of multiple populations in NGC 2808 based on high-precision Hubble Space Telescope proper-motion measurements. In a recent study, Milone et al. have identified five distinct populations (A, B, C, D, and E) in NGC 2808. Populations D and E coincide with the helium-enhanced populations in the middle and the blue main sequences (mMS and bMS) previously discovered by Piotto et al.; populations A, B, and C correspond to the redder main sequence (rMS) that in the Piotto et al. was associated with the primordial stellar population. Our analysis shows that, in the outermost regions probed (between about 1.5 and 2 times the cluster half-light radius), the velocity distribution of populations D and E is radially anisotropic (the deviation from an isotropic distribution is significant at the ~3.5-sigma level). Stars of populations D and E have a smaller tangential velocity dispersion than those of populations A, B, and C, while no significant differences are found in the radial-velocity dispersion. We present the results of a numerical simulation showing that the observed differences between the kinematics of these stellar populations are consistent with the expected kinematic fingerprint of the diffusion towards the cluster outer regions of stellar populations initially more centrally concentrated.
The results of simulations of the extragalactic propagation of ultra-high energy cosmic rays (UHECRs) have intrinsic uncertainties due to poorly known physical quantities and approximations used in the codes. We quantify the uncertainties in the simulated UHECR spectrum and composition due to different models for the extragalactic background light (EBL), different photodisintegration setups, approximations concerning photopion production and the use of different simulation codes. We discuss the results for several representative source scenarios with proton, nitrogen or iron at injection. For this purpose we used SimProp and CRPropa, two publicly available codes for Monte Carlo simulations of UHECR propagation. CRPropa is a detailed and extensive simulation code, while SimProp aims to achieve acceptable results using a simpler code. We show that especially the choices for the EBL model and the photodisintegration setup can have a considerable impact on the simulated UHECR spectrum and composition.
Long-exposure spectro-polarimetry in the near-infrared is a preferred method to measure the magnetic field and other physical properties of solar prominences. In the past, it has been very difficult to observe prominences in this way with sufficient spatial resolution to fully understand their dynamical properties. Solar prominences contain highly transient structures, visible only at small spatial scales; hence they must be observed at sub-arcsecond resolution, with a high temporal cadence. An adaptive optics (AO) system capable of directly locking-on to prominence structure away from the solar limb has the potential to allow for diffraction-limited spectro-polarimetry of solar prominences. In this paper, the performance of the off-limb AO system and its expected performance, at the desired science wavelength {\CaII} 8542A, are shown.
We investigate the time evolution of non-parametric morphological quantities and their relationship to major mergers between $4\geq z \geq 2$ in high-resolution cosmological zoom simulations of disk galaxies that implement kinetic wind feedback, $H_2$-based star formation, and minimal ISM pressurisation. We show that the resulting galaxies broadly match basic observed physical properties of $z\sim 2$ objects. We measure the galaxies' concentrations ($C$), asymmetries ($A$), and $Gini$ ($G$) and $M_{20}$ coefficients, and correlate these with major merger events identified from the mass growth history. We find that high values of asymmetry provide the best indicator for identifying major mergers of $>1:4$ mass ratio within our sample, with $Gini$-$M_{20}\,$ merger classification only as effective for face-on systems and much less effective for edge-on or randomly-oriented galaxies. The canonical asymmetry cut of $A\geq0.35$, however, is only able to correctly identify major mergers $\sim 10\%$ of the time, while a higher cut of $A\geq 0.8$ more efficiently picks out mergers at this epoch. We further examine the temporal correlation between morphological statistics and mergers, and show that for randomly-oriented galaxies, half the galaxies with $A\geq0.8$ undergo a merger within $\pm0.2\,{\rm Gyr}$, whereas $Gini$-$M_{20}\,$ identification only identifies about a third correctly. The fraction improves further using $A\geq 1.5$, but about the half the mergers are missed by this stringent cut.
Astronomical observations are affected by several kinds of noise, each with its own causal source; there is photon noise, stochastic source variability, and residuals coming from imperfect calibration of the detector or telescope. The precision of NASA Kepler photometry for exoplanet science---the most precise photometric measurements of stars ever made---appears to be limited by unknown or untracked variations in spacecraft pointing and temperature, and unmodeled stellar variability. Here we present the Causal Pixel Model (CPM) for Kepler data, a data-driven model intended to capture variability but preserve transit signals. The CPM works at the pixel level so that it can capture very fine-grained information about the variation of the spacecraft. The CPM predicts each target pixel value from a large number of pixels of other stars sharing the instrument variabilities while not containing any information on possible transits in the target star. In addition, we use the target star's future and past (auto-regression). By appropriately separating, for each data point, the data into training and test sets, we ensure that information about any transit will be perfectly isolated from the model. The method has four hyper-parameters (the number of predictor stars, the auto-regressive window size, and two L2-regularization amplitudes for model components), which we set by cross-validation. We determine a generic set of hyper-parameters that works well for most of the stars and apply the method to a corresponding set of target stars. We find that we can consistently outperform (for the purposes of exoplanet detection) the Kepler Pre-search Data Conditioning (PDC) method for exoplanet discovery.
A comparative study of the $\Lambda$ hyperon equations of state of Banik, Hempel and Banyopadhyay (BHB) \citep{bhb} and \citet{shen11} (denoted as HShen $\Lambda$) for core collapse supernova (CCSN) simulations is carried out in this work. The dynamical evolution of a protoneutron star (PNS) into a black hole is investigated in core collapse supernova simulations in the general relativistic one dimensional code using the BHB$\Lambda \phi$ and HShen $\Lambda$ equation of state (EoS) tables and different progenitor models from Woosley and Heger \citep{Woos}. Radial profiles of the mass fractions of baryons, the density as well as the temperature in the PNS at different moments in time, are compared for both EoS tables. The behaviour of the central density of the PNS with time is demonstrated for those two $\Lambda$ hyperon EoS tables and compared with their corresponding nuclear EoS tables. It is observed that the black hole formation time is higher in the BHB$\Lambda \phi$ case than in the HShen $\Lambda$ EoS for the entire set of progenitor models adopted here, because the repulsive $\Lambda$-$\Lambda$ interaction makes the BHB$\Lambda \phi$ EoS stiffer. Neutrino emission with the $\Lambda$ hyperon EoS ceases earlier than that of its nuclear counterpart. The long duration evolution of the shock radius and gravitational mass of the PNS after a successful supernova explosion with enhanced neutrino heating are studied with the BHB$\Lambda \phi$ EoS and $s$20WH07 progenitor model. The PNS is found to remain stable for 4 s and might evolve into a cold neutron star.
We calculate the distribution function of astronomical objects (like galaxies and/or smooth halos of different kinds) gravitational fields due to their tidal in- teraction. For that we apply the statistical method of Chandrasekhar (1943), used there to calculate famous Holtzmark distribution. We show that in our approach the distribution function is never Gaussian, its form being dictated by the potential of interaction between objects. This calculation permits us to perform a theoretical analysis of the relation between angular momentum and mass (richness) of the galaxy clusters. To do so, we follow the idea of Catelan & Theuns (1996) and Heavens & Peacock (1988). The main difference is that here we reduce the problem to discrete many-body case, where all physical properties of the system are determined by the interaction potential V(r_ij). The essence of reduction is that we use the multipole (up to quadrupole here) expansion of Newtonian potential so that all hydrodynamic, "extended" characteristics of an object like its density mass are "integrated out" giving its "point-like" charac- teristics like mass and quadrupole moment. In that sense we make no difference between galaxies and smooth components like halos. We compare our theoretical results with observational data.
We investigate the star formation properties of ~800 sources detected in one of the deepest radio surveys at 1.4 GHz. Our sample spans a wide redshift range (~0.1 - 4) and about four orders of magnitude in star formation rate (SFR). It includes both star forming galaxies (SFGs) and active galactic nuclei (AGNs), further divided into radio-quiet and radio-loud objects. We compare the SFR derived from the far infrared luminosity, as traced by Herschel, with the SFR computed from their radio emission. We find that the radio power is a good SFR tracer not only for pure SFGs but also in the host galaxies of RQ AGNs, with no significant deviation with redshift or specific SFR. Moreover, we quantify the contribution of the starburst activity in the SFGs population and the occurrence of AGNs in sources with different level of star formation. Finally we discuss the possibility of using deep radio survey as a tool to study the cosmic star formation history.
We report on a total of 106 nights of optical interferometric observations of the $\epsilon$ Aurigae system taken during the last 14 years by four beam combiners at three different interferometric facilities. This long sequence of data provides an ideal assessment of the system prior to, during, and after the recent 2009-2011 eclipse. We have reconstructed model-independent images from the 10 in-eclipse epochs which show that a disk-like object is indeed responsible for the eclipse. Using new 3D, time-dependent modeling software, we derive the properties of the F-star (diameter, limb darkening), determine previously unknown orbital elements ($\Omega$, $i$), and access the global structures of the optically thick portion of the eclipsing disk using both geometric models and approximations of astrophysically relevant density distributions. These models may be useful in future hydrodynamical modeling of the system. Lastly, we address several outstanding research questions including mid-eclipse brightening, possible shrinking of the F-type primary, and any warps or sub-features within the disk.
High-energy neutrino astronomy will probe the working of the most violent
phenomena in the Universe. The Giant Radio Array for Neutrino Detection (GRAND)
project consists of an array of $\sim10^5$ radio antennas deployed over
$\sim$200000km$^2$ in a mountainous site. It aims at detecting high-energy
neutrinos via the measurement of air showers induced by the decay in the
atmosphere of $\tau$ leptons produced by the interaction of the cosmic
neutrinos under the Earth surface. Our objective with GRAND is to reach a
neutrino sensitivity of
$3\times10^{-11}E^{-2}$GeV$^{-1}$cm$^{-2}$s$^{-1}$sr$^{-1}$ above $3
\times10^{16}$eV. This sensitivity ensures the detection of cosmogenic
neutrinos in the most pessimistic source models, and about 100 events per year
are expected for the standard models. GRAND would also probe the neutrino
signals produced at the potential sources of UHECRs.
We show how our preliminary design should enable us to reach our sensitivity
goals, and present the experimental characteristics. We assess the possibility
to adapt GRAND to other astrophysical radio measurements. We discuss in this
token the technological options for the detector and the steps to be taken to
achieve the GRAND project.
We study the occurrence of delayed SNe~Ia in the single degenerate (SD) scenario. We assume that a massive carbon-oxygen (CO) white dwarf (WD) accretes matter coming from a companion star, making it to spin at the critical rate. We assume uniform rotation due to magnetic field coupling. The carbon ignition mass for non-rotating WDs is M_{ig}^{NR} \approx 1.38 M_{\odot}; while for the case of uniformly rotating WDs it is a few percent larger (M_{ig}^{R} \approx 1.43 M_{\odot}). When accretion rate decreases, the WD begins to lose angular momentum, shrinks, and spins up; however, it does not overflow its critical rotation rate, avoiding mass shedding. Thus, angular momentum losses can lead the CO WD interior to compression and carbon ignition, which would induce an SN~Ia. The delay, largely due to the angular momentum losses timescale, may be large enough to allow the companion star to evolve to a He WD, becoming undetectable at the moment of explosion. This scenario supports the occurrence of delayed SNe~Ia if the final CO WD mass is 1.38 M_{\odot} < M < 1.43 M_{\odot}. We also find that if the delay is longer than ~3 Gyr, the WD would become too cold to explode, rather undergoing collapse.
Multiple or extended turnoffs in young clusters in the Magellanic Clouds have recently received large attention. A number of studies have shown that they may be interpreted as the result of a significant age spread (several 10^8yr in clusters aged 1--2 Gyr), while others attribute them to a spread in stellar rotation. We focus on the cluster NGC 1856, showing a splitting in the upper part of the main sequence, well visible in the color m_{F336W}-m_{F555W}$, and a very wide turnoff region. Using population synthesis available from the Geneva stellar models, we show that the cluster data can be interpreted as superposition of two main populations having the same age (~350Myr), composed for 2/3 of very rapidly rotating stars, defining the upper turnoff region and the redder main sequence, and for 1/3 of slowly/non-rotating stars. Since rapid rotation is a common property of the B-A type stars, the main question raised by this model concerns the origin of the slowly/non-rotating component. Binary synchronization is a possible process behind the slowly/non-rotating population; in this case, many slowly/non-rotating stars should still be part of binary systems with orbital periods in the range from 4 to 500 days. Such periods imply that Roche lobe overflow occurs, during the evolution of the primary off the main sequence, so most primaries may not be able to ignite core helium burning, consistently why the lack of a red clump progeny of the slowly rotating population.
Here we aim to study the physical and kinematical characteristics of the unstudied old planetary nebula (PN) PN G342.0-01.7, which shows evidence of interaction with its surrounding interstellar medium. We used Integral Field Spectra from the Wide Field Spectrograph on the ANU 2.3 m telescope to provide spectroscopy across the whole object covering the spectral range 3400-7000 {\AA}. We formed narrow-band images to investigate the excitation structure. The spectral analysis shows that the object is a distant Peimbert Type I PN of low excitation, formally of excitation class of 0.5. The low electron density, high dynamical age, and low surface brightness of the object confirm that it is observed fairly late in its evolution. It shows clear evidence for dredge-up of CN-processed material characteristic of its class. In addition, the low peculiar velocity of 7 km s$^{-1}$ shows it to be a member of the young disk component of our Galaxy. We built a self-consistent photoionisation model for the PNe matching the observed spectrum, the H$\beta$ luminosity, and the diameter. On the basis of this we derive an effective temperature $\log T_{\rm eff} \sim 5.05$ and luminosity $1.85 < \log L < 2.25$. The temperature is much higher than might have been expected using the excitation class, proving that this can be misleading in classifying evolved PNe. PN G342.0-01.7 is in interaction with its surrounding interstellar medium through which the object is moving in the south-west direction. This interaction drives a slow shock into the outer PN ejecta. A shock model suggests that it only accounts for about 10\% of the total luminosity, but has an important effect on the global spectrum of the PN.
Homologous flares are flares that occur repetitively in the same active region, with similar structure and morphology. A series of at least eight homologous flares occurred in active region NOAA 11237 over 16 - 17 June 2011. A nearby prominence/filament was rooted in the active region, and situated near the bottom of a coronal cavity. The active region was on the southeast solar limb as seen from SDO/AIA, and on the disk as viewed from STEREO/EUVI-B. The dual perspective allows us to study in detail behavior of the prominence/filament material entrained in the magnetic field of the repeatedly-erupting system. Each of the eruptions was mainly confined, but expelled hot material into the prominence/filament cavity system (PFCS). The field carrying and containing the ejected hot material interacted with the PFCS and caused it to inflate, resulting in a step-wise rise of the PFCS approximately in step with the homologous eruptions. The eighth eruption triggered the PFCS to move outward slowly, accompanied by a weak coronal dimming. As this slow PFCS eruption was underway, a final ejective flare occurred in the core of the active region, resulting in strong dimming in the EUVI-B images and expulsion of a coronal mass ejection (CME). A plausible scenario is that the repeated homologous flares could have gradually destabilized the PFCS, and its subsequent eruption removed field above the acitive region and in turn led to the ejective flare, strong dimming, and CME.
Spectra of 34 H II regions in the late-type galaxies NGC1087, NGC2967, NGC3023, NGC4030, NGC4123, and NGC4517A were observed with the South African Large Telescope (SALT). In all 34 H II regions, oxygen abundances were determined through the "counterpart" method (C method). Additionally, in two H II regions in which the auroral lines were detected oxygen abundances were measured through the classic Te method. We also estimated the abundances in our H II regions using the O3N2 and N2 calibrations and compared those with the C-based abundances. With these data we examined the radial abundance distributions in the disks of our target galaxies. We derived surface-brightness profiles and other characteristics of the disks (the surface brightness at the disk center and the disk scale length) in three photometric bands for each galaxy using publicly available photometric imaging data. The radial distributions of the oxygen abundances predicted by the relation between abundance and disk surface brightness in the W1 band obtained for spiral galaxies in our previous study are close to the radial distributions of the oxygen abundances determined from the analysis of the emission line spectra for four galaxies where this relation is applicable. Hence, when the surface-brightness profile of a late-type galaxy is known, this parametric relation can be used to estimate the likely present-day oxygen abundance in its disk.
We present a study of the kinematics and the physical properties of the central region of the Hickson Compact Group 31 (HCG 31), focusing on the HCG 31A+C system, using integral field spectroscopy data taken with the Gemini-south telescope. The main players in the merging event (galaxies A and C) are two dwarf galaxies that in the past have already had one close encounter, given the observed tidal tails, and may now be in their second approach, possibly about to merge. We present new velocity fields and H{\alpha} emission, stellar continuum, velocity dispersion, electron density, H{\alpha} equivalent width and age maps. Considering the high spatial resolution of the IFU data, we were able to measure various components and estimate their physical parameters, spatially resolving the different structures in this region. Our main findings are the following: (1) we report for the first time the presence of a super stellar cluster next to the burst associated to the HCG 31C central blob, related to the high values of velocity dispersion observed in this region as well as to the highest value of stellar continuum emission. This may suggest that this system is cleaning its environment through strong stellar winds that may then trigger a strong star formation event in its neighborhood, (2) among other physical parameters, we estimate an L(H{\alpha})~14x10^{41} erg/s and a SFR~11 Msol/yr for the central merging region of HCG 31 A+C. These values indicate a high star formation density, suggesting that the system is part of a merging object, supporting previous scenarios proposed for this system.
The high-velocity features (HVFs) in optical spectra of type Ia supernovae (SNe Ia) are examined with a large sample including very early-time spectra (e.g., t < -7 days). Multiple Gaussian fits are applied to examine the HVFs and their evolutions, using constraints on expansion velocities for the same species (i.e., SiII 5972 and SiII 6355). We find that strong HVFs tend to appear in SNe Ia with smaller decline rates (e.g., dm15(B)<1.4 mag), clarifying that the finding by Childress et al. (2014) for the Ca-HVFs in near-maximum-light spectra applies both to the Si-HVFs and Ca-HVFs in the earlier phase. The Si-HVFs seem to be more common in fast-expanding SNe Ia, which is different from the earlier result that the Ca-HVFs are associated with SNe Ia having slower SiII 6355 velocities at maximum light (i.e., Vsi). This difference can be due to that the HVFs in fast-expanding SNe Ia usually disappear more rapidly and are easily blended with the photospheric components when approaching the maximum light. Moreover, SNe Ia with both stronger HVFs at early phases and larger Vsi are found to have noticeably redder B-V colors and occur preferentially in the inner regions of their host galaxies, while those with stronger HVFs but smaller Vsi show opposite tendencies, suggesting that these two subclasses have different explosion environments and their HVFs may have different origins. We further examine the relationships between the absorption features of SiII 6355 and CaII IR lines, and find that their photospheric components are well correlated in velocity and strength but the corresponding HVFs show larger scatter. These results cannot be explained with ionization and/or thermal processes alone, and different mechanisms are required for the creation of HVF-forming region in SNe Ia.
We study the radio/X-ray correlation in Cyg X-3. It has been known that the soft and hard X-ray fluxes in the hard spectral state are correlated positively and negatively, respectively, with the radio flux. We show that this implies that the observed $\sim$2--100 keV flux (which is a fair approximation to the bolometric flux) is completely uncorrelated with the radio flux. We can recover a positive correlation (seen in other sources and expected theoretically) if the soft X-rays are strongly absorbed by a local medium. Then, however, the intrinsic X-ray spectrum of Cyg X-3 in its hard state is relatively soft, similar to that of an intermediate spectral state of black-hole binaries, but not to their true hard state. We also find the radio spectra in the hard state of Cyg X-3 to be hard on average, and the flux distributions of the radio emission and soft X-rays to follow a sum of two log-normal functions.
This is a follow-up sensitivity study on r-mode gravitational wave signals from newborn neutron stars illustrating the applicability of machine learning algorithms for the detection of long-lived gravitational-wave transients. In this sensitivity study we examine three machine learning algorithms (MLAs): artificial neural networks (ANNs), support vector machines (SVMs) and constrained subspace classifiers (CSCs). The objective of this study is to compare the detection efficiency that MLAs can achieve with the efficiency of conventional detection algorithms discussed in an earlier paper. Comparisons are made using 2 distinct r-mode waveforms. For the training of the MLAs we assumed that some information about the distance to the source is given so that the training was performed over distance ranges not wider than half an order of magnitude. The results of this study suggest that machine learning algorithms are suitable for the detection of long-lived gravitational-wave transients and that when assuming knowledge of the distance to the source, MLAs are at least as efficient as conventional methods.
We present the GAMA Panchromatic Data Release (PDR) constituting over
230deg$^2$ of imaging with photometry in 21 bands extending from the far-UV to
the far-IR. These data complement our spectroscopic campaign of over 300k
galaxies, and are compiled from observations with a variety of facilities
including: GALEX, SDSS, VISTA, WISE, and Herschel, with the GAMA regions
currently being surveyed by VST and scheduled for observations by ASKAP. These
data are processed to a common astrometric solution, from which photometry is
derived for 221,373 galaxies with r<19.8 mag. Online tools are provided to
access and download data cutouts, or the full mosaics of the GAMA regions in
each band.
We focus, in particular, on the reduction and analysis of the VISTA VIKING
data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining
the data and examining its integrity. Having derived the 21-band photometric
catalogue we proceed to fit the data using the energy balance code MAGPHYS.
These measurements are then used to obtain the first fully empirical
measurement of the 0.1-500$\mu$m energy output of the Universe. Exploring the
Cosmic Spectral Energy Distribution (CSED) across three time-intervals
(0.3-1.1Gyr, 1.1-1.8~Gyr and 1.8---2.4~Gyr), we find that the Universe is
currently generating $(1.5 \pm 0.3) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$, down
from $(2.5 \pm 0.2) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$ 2.3~Gyr ago. More
importantly, we identify significant and smooth evolution in the integrated
photon escape fraction at all wavelengths, with the UV escape fraction
increasing from 27(18)% at z=0.18 in NUV(FUV) to 34(23)% at z=0.06.
The GAMA PDR will allow for detailed studies of the energy production and
outputs of individual systems, sub-populations, and representative galaxy
samples at $z<0.5$. The GAMA PDR can be found at: this http URL
Ongoing and future wide-field galaxy surveys can be used to locate a number of clusters of galaxies with cosmic shear measurement alone. We study constraints on cosmological models using statistics of weak lensing selected galaxy clusters. We extend our previous theoretical framework to model the statistical properties of clusters in variants of cosmological models as well as in the standard LCDM model. Weak lensing selection of clusters does not rely on the conventional assumption such as the relation between luminosity and mass and/or hydrostatic equilibrium, but a number of observational effects compromise robust identification. We use a large set of realistic mock weak-lensing catalogs as well as analytic models to perform a Fisher analysis and make forecast for constraining two competing cosmological models, wCDM model and f(R) model proposed by Hu & Sawicki, with our lensing statistics. We show that weak lensing selected clusters are excellent probe of cosmology when combined with cosmic shear power spectrum even in presence of galaxy shape noise and masked regions. With the information of weak lensing selected clusters, the precision of cosmological parameter estimate can be improved by a factor of ~1.3 and ~10 for wCDM model and f(R) model, respectively. Hyper Suprime-Cam survey with sky coverage of 1250 squared degrees can constrain the equation of state of dark energy w_0 with a level of Delta w_0 ~0.1. It can also constrain the additional scalar degree of freedom in f(R) model with a level of |f_R0| ~6x10^{-6}, when constraints from cosmic microwave background measurements are incorporated. Future weak lensing surveys with sky coverage of 20,000 squared degrees will place tighter constraints on w_0 and |f_R0| even without cosmic microwave background measurements.
We propose a new electromagnetic-emission mechanism in magnetized, force-free plasma, which is driven by the evolution of the underlying dynamical spacetime. With this radiation-generation process, gravitational energy is converted into electromagnetic energy, which is then carried away by both fast-magnetosonic and Alfven waves of the plasma. As an immediate demonstration, we consider compact binary mergers occurring within magnetized plasma, which have been shown by previous numerical studies to produce copious amounts of electromagnetic radiation. The emission power and angular distribution of the two classes of waves are separately determined. When the new process is combined with previously understood mechanisms such as the Blandford-Znajek process and kinetic-motion-driven radiation, one can classify different components of electromagnetic emissions seen in the inspiral stage of compact-binary coalescence.
We present a family of self-consistent, spherical, lowered isothermal models, consisting of one or more mass components, with parameterised prescriptions for the energy truncation and for the amount of radially biased pressure anisotropy. The models are particularly suited to describe the phase-space density of stars in tidally limited, mass-segregated star clusters in all stages of their life-cycle. The models extend a family of isotropic, single-mass models by Gomez-Leyton and Velazquez, of which the well-known Woolley, King and Wilson (in the non-rotating and isotropic limit) models are members. We derive analytic expressions for the density and velocity dispersion components in terms of potential and radius, and introduce a fast model solver in PYTHON (LIMEPY), that can be used for data fitting or for generating discrete samples.
One of the most energetic gamma-ray bursts GRB 110731A, was observed from optical to GeV energy range by Fermi and Swift Observatories, and by the MOA and GROND optical telescopes. The multiwavelength observations over different epochs (from trigger time to more than 800 s) showed that the spectral energy distribution was better fitted by a wind afterglow model. We present a leptonic model based on an early afterglow that evolves in a stellar wind to describe the multiwavelength light curves observations. In particular, the origin of the LAT emission is explained through the superposition of synchrotron radiation from the forward shock and synchrotron self-Compton emission from the reverse shock. The bulk Lorentz factor required in this model is $\Gamma\simeq520$ and the result suggests that the ejecta must be magnetized.
Astronomical data is often uncertain with errors that are heteroscedastic (different for each data point) and covariant between different dimensions. Assuming that a set of D-dimensional data points can be described by a (D - 1)-dimensional plane with intrinsic scatter, we derive the general likelihood function to be maximised to recover the best fitting model. Alongside the mathematical description, we also release the hyper-fit package for the R statistical language (github.com/asgr/hyper.fit) and a user-friendly web interface for online fitting (hyperfit.icrar.org). The hyper-fit package offers access to a large number of fitting routines, includes visualisation tools, and is fully documented in an extensive user manual. Most of the hyper-fit functionality is accessible via the web interface. In this paper we include applications to toy examples and to real astronomical data from the literature: the mass-size, Tully-Fisher, Fundamental Plane, and mass-spin-morphology relations. In most cases the hyper-fit solutions are in good agreement with published values, but uncover more information regarding the fitted model.
Dark matter constitutes the great majority of the matter content in the Universe, but its microscopic nature remains an intriguing mystery, with profound implications for particle physics, astrophysics and cosmology. Here we shed light on the longstanding issue of whether the dark matter is warm or cold by combining the measurements of the galaxy luminosity functions out to high redshifts z~10 from the Hubble Space Telescope with the recent cosmological data on the reionization history of the Universe from the Planck mission. We derive robust and tight bounds on the mass of warm dark matter particle, finding that the current data require it to be in the narrow range between 2 and 3 keV. In addition, we show that a mass not exceeding 3 keV is also concurrently indicated by astrophysical constraints related to the local number of satellites in Milky Way-sized galaxies, though it is in marginal tension with analysis of the Lyman-alpha forest. For warm dark matter masses above 3 keV as well as for cold dark matter, to satisfy the Planck constraints on the optical depth and not to run into the satellite problem would require invoking astrophysical processes that inhibit galaxy formation in halos with mass M_H< a few 10^8 M_sun, corresponding to a limiting UV magnitude M_UV~-11. Anyway, we predict a downturn of the galaxy luminosity function at z~8 faintward of M_UV~-12, and stress that its detailed shape is extremely informative both on particle physics and on the astrophysics of galaxy formation in small halos. These expectations will be tested via the Hubble Frontier Fields and with the advent of the James Webb Space Telescope, which will enable probing the very faint end of the galaxy luminosity function out to z~8-10.
In previous work we identified six Sun-like stars observed by Kepler with exceptionally clear asteroseismic signatures of rotation. Here, we show that five of these stars exhibit surface variability suitable for measuring rotation. In order to further constrain differential rotation, we compare the rotation periods obtained from light-curve variability with those from asteroseismology. The two rotation measurement methods are found to agree within uncertainties, suggesting that radial differential rotation is weak, as is the case for the Sun. Furthermore, we find significant discrepancies between ages from asteroseismology and from three different gyrochronology relations, implying that stellar age estimation is problematic even for Sun-like stars.
The power of micro-arcsecond ($\mu$as) astrometry is about to be unleashed. ESA's Gaia mission, now headed towards the end of the first year of routine science operations, will soon fulfil its promise for revolutionary science in countless aspects of Galactic astronomy and astrophysics. The potential of Gaia position measurements for important contributions to the astrophysics of planetary systems is huge. We focus here on the expectations for detection and improved characterization of 'young' planetary systems in the neighborhood of the Sun using a combination of Gaia $\mu$as astrometry and direct imaging techniques.
We studied the unbiased optical brightness distribution which was calculated from the survival analysis of host galaxies and its relationship with the Swift GRB data of the host galaxies observed by the Keck telescopes. Based on the sample obtained from merging the Swift GRB table and the Keck optical data we also studied the dependence of this distribution on the data of the GRBs. Finally, we compared the HGs distribution with standard galaxies distribution which is in the DEEP2 galaxies catalog.
Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $d\log\rho/d\log r \approx -1$, and some dispersion, which is generally parametrized by the varying index $\alpha$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp - which has only one shape parameter, $\phi_0$ - fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp $\phi_0$, Einasto $\alpha$, or the central logarithmic slope in the Dehnen-Tremaine analytical $\gamma$-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos.
Supermassive black hole dynamics during galaxy mergers is crucial in determining the rate of black hole mergers and cosmic black hole growth. As simulations achieve higher resolution, it becomes important to assess whether the black hole dynamics is influenced by the treatment of the interstellar medium in different simulation codes. We here compare simulations of black hole growth in galaxy mergers with two codes: the Smoothed Particle Hydrodynamics code Gasoline, and the Adaptive Mesh Refinement code Ramses. We seek to identify predictions of these models that are robust despite differences in hydrodynamic methods and implementations of sub-grid physics. We find that the general behavior is consistent between codes. Black hole accretion is minimal while the galaxies are well-separated (and even as they "fly-by" within 10 kpc at first pericenter). At late stages, when the galaxies pass within a few kpc, tidal torques drive nuclear gas inflow that triggers bursts of black hole accretion accompanied by star formation. We also note quantitative discrepancies that are model-dependent: our Ramses simulations show less star formation and black hole growth, and a smoother gas distribution with larger clumps and filaments, than our Gasoline simulations. We attribute these differences primarily to the sub-grid models for black hole fueling and feedback and gas thermodynamics. The main conclusion is that differences exist quantitatively between codes, and this should be kept in mind when making comparisons with observations. However, reassuringly, both codes capture the same dynamical behaviors in terms of triggering of black hole accretion, star formation, and black hole dynamics.
We calculate the cosmic microwave background temperature bispectrum from cosmic strings, for the first time including the contributions from the last scattering surface, using a well-established Gaussian model for the string energy-momentum correlation functions, and a simplified model for the cosmic fluid. We check our approximation for the integrated Sachs-Wolfe (ISW) contribution against the bispectrum obtained from the full sky map of the cosmic string ISW signal used by the Planck team, obtaining good agreement. We validate our model for the last scattering surface contribution by comparing the predicted temperature power spectrum with that obtained from a full Boltzmann code treatment applied to the Unconnected Segment Model of a string network. We find that including the last scattering contribution has only a small impact on the upper limit on the string tension resulting from the bispectrum at Planck resolutions, and argue that the bispectrum is unlikely to be competitive with the power spectrum at any resolution.
The present work is based on a description for the angular mometum loss rate due to magnetic braking for main-sequence stars on the relationship between stellar rotation and age. In general, this loss rate denoted by $\mathrm dJ/\mathrm dt$ depends on angular velocity $\Omega$ in the form $\mathrm dJ/\mathrm dt\propto\Omega^{q}$, where $q$ is a parameter extracted from nonextensive statistical mechanics. Already, in context of stellar rotation, this parameter is greater than unity and it is directly related to the braking index. For $q$ equal to unity, the scenario of saturation of the magnetic field is recovered, otherwise $q$ indicates an unsaturated field. This new approach have been proposed and investigated by de Freitas \& De Medeiros for unsaturated field stars. In present work, we propose a nonextensive approach for the stellar rotational evolution based on the Reiners \& Mohanthy model. In this sense, we developed a nonextensive version of Reiners \& Mohanthy torque and also compare this generalized version with the model proposed in de Freitas \& De Medeiros based on the spin-down Kawaler torque for the main-sequence stars with F and G spectral types. We use the same sample of $\sim16 000$ field stars with rotational velocity $v \sin i$ limited in age and mass. As a result, we show that the Kawaler and Reiners \& Mohanthy models exhibit strong discrepancies, mainly in relation to the domain of validity of the entropic index $q$. These discrepancies are mainly due to sensitivity on the stellar radius. Finally, our results showed that modified Kawaler prescription is compatible with a wider mass range, while the Reiners \& Mohanty model is restricted to masses less than G6 stars.
The results obtained using the temperature monitoring systems of the 6-m BTA telescope primary mirror, dome space, and external environment are reported. We consider the factors that affect the development of microturbulence in the near-mirror air layer and inside the dome space, variation of the telescope focal length with the temperature of its structures, variation of seeing due to temperature gradients inside the primary mirror of the 6-m telescope. The methods used in various observatories for reducing microturbulence are analyzed. We formulate suggestions concerning the improvement of the temperature monitoring system currently in operation and the system of automatic adjustment of the telescope focal length to compensate the thermal drift of the focus during observations.
The emission of supernova remnants reflects the properties of both the progenitor supernovae and the surrounding environment. The complex morphology of the remnants, however, hampers the disentanglement of the two contributions. Here we aim at identifying the imprint of SN 1987A on the X-ray emission of its remnant and at constraining the structure of the environment surrounding the supernova. We performed high-resolution hydrodynamic simulations describing SN 1987A soon after the core-collapse and the following three-dimensional expansion of its remnant between days 1 and 15000 after the supernova. We demonstrated that the physical model reproducing the main observables of SN 1987A during the first 250 days of evolution reproduces also the X-ray emission of the subsequent expanding remnant, thus bridging the gap between supernovae and supernova remnants. By comparing model results with observations, we constrained the explosion energy in the range $1.2-1.4\times 10^{51}$~erg and the envelope mass in the range $15-17 M_{\odot}$. We found that the shape of X-ray lightcurves and spectra at early epochs (<15 years) reflects the structure of outer ejecta: our model reproduces the observations if the outermost ejecta have a post-explosion radial profile of density approximated by a power law with index $\alpha = -8$. At later epochs, the shapes of X-ray lightcurves and spectra reflect the density structure of the nebula around SN 1987A. This enabled us to ascertain the origin of the multi-thermal X-ray emission, to disentangle the imprint of the supernova on the remnant emission from the effects of the remnant interaction with the environment, and to constrain the pre-supernova structure of the nebula.
MIRI (the Mid-Infrared Instrument for the James Webb Space Telescope (JWST)) operates from 5 to 28.5 microns and combines over this range: 1.) unprecedented sensitivity levels; 2.) sub-arcsec angular resolution; 3.) freedom from atmospheric interference; 4.) the inherent stability of observing in space; and 5.) a suite of versatile capabilities including imaging, low and medium resolution spectroscopy (with an integral field unit), and coronagraphy. We illustrate the potential uses of this unique combination of capabilities with various science examples: 1.) imaging exoplanets; 2.) transit and eclipse spectroscopy of exoplanets; 3.) probing the first stages of star and planet formation, including identifying bioactive molecules; 4.) determining star formation rates and mass growth as galaxies are assembled; and 5.) characterizing the youngest massive galaxies. This paper is the introduction to a series of ten covering all aspects of the instrument.
We present the first time-simultaneous high angular resolution spectral energy distribution (SED) of the core of M87 at a scale of 0.4 arcsecs across the electromagnetic spectrum. Two activity periods of the core of M87 are sampled: a quiescent mode, representative of the most common state of M87, and an active one, represented by the outburst occurring in 2005. The main difference between both SEDs is a shift in flux in the active SED by a factor of about two, their shapes remaining similar across the entire spectrum. The shape of the compiled SEDs is remarkably different from those of active galactic nuclei (AGN). It lacks three major AGN features: the IR bump, the inflection point at about 1 micron and the blue bump. The SEDs also differ from the spectrum of a radiatively inefficient accretion flow. Down to the scales of ~12 pc from the centre, we find that the emission from a jet gives an excellent representation of the spectrum over ten orders of magnitude in frequency for both the active and the quiescent phases of M87. The inferred total jet power is one to two orders of magnitude lower than the jet mechanical energy inferred from various methods in the literature. This discrepancy cannot easily be ascribed to variability. Yet, our measurements regard the inner few parsecs which might provide a genuine account of the jet power at the base. We derive a strict upper limit to the accretion rate of 6 x 10E-5 Mo / yr, assuming 10% efficiency. The inferred accretion power can account for M87 radiative luminosity at the jet-frame assuming boosting factors larger than 10, it is however two orders of magnitude below that required to account for M87 jet kinetic power. We thus propose that energy tapped from the black hole spin may be a complementary source to power the jet of M87, a large supply of accreting gas becoming thus unnecessary.
We investigate the transition to Self Organized Criticality in a two-dimensional model of a flux tube with a background flow. The magnetic induction equation, represented by a partial differential equation with a stochastic source term, is discretized and implemented on a two dimensional cellular automaton. The energy released by the automaton during one relaxation event is the magnetic energy. As a result of the simulations we obtain the time evolution of the energy release, of the system control parameter, of the event lifetime distribution and of the event size distribution, respectively, and we establish that a Self Organized Critical state is indeed reached by the system. Moreover, energetic initial impulses in the magnetohydrodynamic flow can lead to one dimensional signatures in the magnetic two dimensional system, once the Self Organized Critical regime is established. The applications of the model for the study of Gamma Ray Bursts is briefly considered, and it is shown that some astrophysical parameters of the bursts, like the light curves, the maximum released energy, and the number of peaks in the light curve can be reproduced and explained, at least on a qualitative level, by working in a framework in which the systems settles in a Self Organized Critical state via magnetic reconnection processes in the magnetized Gamma Ray Burst fireball.
We propose a new approach to the missing baryons problem. Building on the common assumption that the missing baryons are in the form of the Warm Hot Intergalactic Medium (WHIM), we further assumed here that the galaxy luminosity density can be used as a tracer of the WHIM. The latter assumption is supported by our finding of a significant correlation between the WHIM density and the galaxy luminosity density in the hydrodynamical simulations of Cui et al. (2012). We further found that the fraction of the gas mass in the WHIM phase is substantially (by a factor of $\sim$1.6) higher within the large scale galactic filaments, i.e. $\sim$70\%, compared to the average in the full simulation volume of $\sim$0.1\,Gpc$^3$. The relation between the WHIM overdensity and the galaxy luminosity overdensity within the galactic filaments is consistent with linear: $\delta_{\rm whim}\,=\,0.7\,\pm\,0.1\,\times\,\delta_\mathrm{LD}^{0.9 \pm 0.2}$. We applied our procedure to the line of sight to the blazar H2356-309 and found evidence for the WHIM in correspondence of the Sculptor Wall (z $\sim$0.03 and $\log{N_H}$ = $19.9^{+0.1}_{-0.3}$) and Pisces-Cetus superclusters (z $\sim$0.06 and $\log{N_H}$ = $19.7^{+0.2}_{-0.3}$), in agreement with the redshifts and column densities of the X-ray absorbers identified and studied by Fang et al. (2010) and Zappacosta et al. (2010). This agreement indicates that the galaxy luminosity density and galactic filaments are reliable signposts for the WHIM and that our method is robust in estimating the WHIM density. The signal that we detected cannot originate from the halos of the nearby galaxies since they cannot account for the large WHIM column densities that our method and X-ray analysis consistently find in the Sculptor Wall and Pisces-Cetus superclusters.
The three-body problem, which describes three masses interacting through Newtonian gravity without any restrictions imposed on the initial positions and velocities of these masses, has attracted the attention of many scientists for more than 300 years. In this paper, we present a review of the three-body problem in the context of both historical and modern developments. We describe the general and restricted (circular and elliptic) three-body problems, different analytical and numerical methods of finding solutions, methods for performing stability analysis, search for periodic orbits and resonances, and application of the results to some interesting astronomical and space dynamical settings. We also provide a brief presentation of the general and restricted relativistic three-body problems, and discuss their astronomical applications.
SN2014J is the closest supernova of type Ia that occured in the last 40 years. This provides an opportunity for unprecedented observational detail and coverage in many astronomical bands, which will help to better understand the still unknown astrophysics of these supernovae. For the first time, such an event occurs sufficiently nearby so that also gamma rays are able to contribute to such investigations. This is important, as the primary source of the supernova light is the radioactive energy from about 0.5 M$_\odot$ of $^{56}$Ni produced in the explosion, and the gamma rays associated with this decay make the supernova shine for months. The INTEGRAL gamma-ray observatory of ESA has followed the supernova emission for almost 5 months. The characteristic gamma ray lines from the $^{56}$Ni decay chain through $^{56}$Co to $^{56}$Fe have been measured. We discuss these observations, and the implications of the measured gamma-ray line characteristics as they evolve.
We study the stability of cosmic string wakes against the disruption by the dominant Gaussian fluctuations which are present in cosmological models. We find that for a string tension given by $G \mu = 10^{-7}$ wakes remain locally stable until a redshift of $z = 6$, and for a value of $G \mu = 10^{-14}$ they are stable beyond a redshift of $z = 20$. We study a global stability criterion which shows that wakes created by strings at times after $t_{eq}$ are identifiable up to the present time, independent of the value of $G \mu$. Taking into account our criteria it is possible to develop strategies to search for the distinctive position space signals in cosmological maps which are induced by wakes.
The Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST) provides measurements over the wavelength range 5 to 28.5 microns. MIRI has, within a single 'package', four key scientific functions: photometric imaging, coronagraphy, single-source low-spectral resolving power (R ~ 100) spectroscopy, and medium-resolving power (R ~ 1500 to 3500) integral field spectroscopy. An associated cooler system maintains MIRI at its operating temperature of < 6.7 K. This paper describes the driving principles behind the design of MIRI, the primary design parameters, and their realization in terms of the 'as-built' instrument. It also describes the test program that led to delivery of the tested and calibrated Flight Model to NASA in 2012, and the confirmation after delivery of the key interface requirements.
Context. Halo occupation distribution (HOD) is a powerful statistic that allows the study of several aspects of the matter distribution in the Universe, such as evaluating semi-analytic models of galaxy formation or imposing constraints on cosmological models. Consequently, it is important to have a reliable method for estimating this statistic, taking full advantage of the available information on current and future galaxy surveys. Aims. The main goal of this project is to combine photometric and spectroscopic information using a discount method of background galaxies in order to extend the range of absolute magnitudes and to increase the upper limit of masses in which the HOD is estimated. We also evaluate the proposed method and apply it to estimating the HOD on the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) galaxy survey. Methods. We propose the background subtraction technique to mel information provided by spectroscopic galaxy groups and photometric survey of galaxies. To evaluate the feasibility of the method, we implement the proposed technique on a mock catalogue built from a semi-analytic model of galaxy formation. Furthermore, we apply the method to the SDSS DR7 using a galaxy group catalogue taken from spectroscopic version and the corresponding photometric galaxy survey. Results. We demonstrated the validity of the method using the mock catalogue.We applied this technique to obtain the SDSS DR7 HOD in absolute magnitudes ranging from $M=-21.5$ to $M=-16.0$ and masses up to $\simeq 10^{15} M_{\odot}$ throughout this range. On the brighter extreme, we found that our results are in excellent agreement with those obtained in previous works.
The imaging channel on the Mid-Infrared Instrument (MIRI) is equipped with four coronagraphs that provide high contrast imaging capabilities for studying faint point sources and extended emission that would otherwise be overwhelmed by a bright point-source in its vicinity. Such bright sources might include stars that are orbited by exoplanets and circumstellar material, mass-loss envelopes around post-main-sequence stars, the near-nuclear environments in active galaxies, and the host galaxies of distant quasars. This paper describes the coronagraphic observing modes of MIRI, as well as performance estimates based on measurements of the MIRI flight model during cryo-vacuum testing. A brief outline of coronagraphic operations is also provided. Finally, simulated MIRI coronagraphic observations of a few astronomical targets are presented for illustration.
We argue that the stellar velocity dispersion observed in an elliptical galaxy is a good proxy for the halo velocity dispersion. As dark matter halos are almost completely characterized by a single scale parameter, the stellar velocity dispersion tells us the virial radius of the halo and the mass contained within. This permits non-dimensionalizing of the stellar mass and effective radius axes of the stellar mass fundamental plane by the virial radius and halo mass, respectively.
The MIRI Si:As IBC detector arrays extend the heritage technology from the Spitzer IRAC arrays to a 1024 x 1024 pixel format. We provide a short discussion of the principles of operation, design, and performance of the individual MIRI detectors, in support of a description of their operation in arrays provided in an accompanying paper (Ressler et al. (2015)). We then describe modeling of their response. We find that electron diffusion is an important component of their performance, although it was omitted in previous models. Our new model will let us optimize the bias voltage while avoiding avalanche gain. It also predicts the fraction of the IR-active layer that is depleted (and thus contributes to the quantum efficiency) as signal is accumulated on the array amplifier. Another set of models accurately predicts the nonlinearity of the detector-amplifier unit and has guided determination of the corrections for nonlinearity. Finally, we discuss how diffraction at the interpixel gaps and total internal reflection can produce the extended cross-like artifacts around images with these arrays at short wavelengths, ~ 5 microns. The modeling of the behavior of these devices is helping optimize how we operate them and also providing inputs to the development of the data pipeline.
Detailed observational characterization of transiting exoplanet systems has revealed that the spin-axes of massive (M > ~1.2 solar masses) stars often exhibit substantial misalignments with respect to the orbits of the planets they host. Conversely, lower-mass stars tend to only have limited obliquities. A similar trend has recently emerged within the observational dataset of young stars' magnetic field strengths: massive T-Tauri stars tend to have dipole fields that are ~10 times weaker than their less-massive counterparts. Here we show that the associated dependence of magnetic star-disk torques upon stellar mass naturally explains the observed spin-orbit misalignment trend, provided that misalignments are obtained within the disk-hosting phase. Magnetic torques act to realign the stellar spin-axes of lower-mass stars with the disk plane on a timescale significantly shorter than the typical disk lifetime, whereas the same effect operates on a much longer timescale for massive stars. Cumulatively, our results point to a primordial excitation of extrasolar spin-orbit misalignment, signalling consistency with disk-driven migration as the dominant transport mechanism for short-period planets. Furthermore, we predict that spin-orbit misalignments in systems where close-in planets show signatures of dynamical, post-nebular emplacement will not follow the observed correlation with stellar mass.
We present a 30 - 50 GHz survey of Sagittarius B2(N) conducted with the
Australia Telescope Compact Array (ATCA) with 5 - 10 arcsec resolution. This
work releases the survey data and demonstrates the utility of scripts that
perform automated spectral line fitting on broadband line data. We describe the
line-fitting procedure, evaluate the performance of the method, and provide
access to all data and scripts. The scripts are used to characterize the
spectra at the positions of three HII regions, each with recombination line
emission and molecular line absorption. Towards the most line-dense of the
three regions characterised in this work, we detect ~500 spectral line
components of which ~90 per cent are confidently assigned to H and He
recombination lines and to 53 molecular species and their isotopologues.
The data reveal extremely subthermally excited molecular gas absorbing
against the continuum background at two primary velocity components. Based on
the line radiation over the full spectra, the molecular abundances and line
excitation in the absorbing components appear to vary substantially towards the
different positions, possibly indicating that the two gas clouds are located
proximate to the star forming cores instead of within the envelope of Sgr B2.
Furthermore, the spatial distributions of species including CS, OCS, SiO, and
HNCO indicate that the absorbing gas components likely have high UV-flux.
Finally, the data contain line-of-sight absorption by $\sim$15 molecules
observed in translucent gas in the Galactic Center, bar, and intervening spiral
arm clouds, revealing the complex chemistry and clumpy structure of this gas.
Formamide (NH$_2$CHO) is detected for the first time in a translucent cloud.
We recently improved the famous result of Parikh and Wilczek, who found a probability of emission of Hawking radiation which is compatible with a non-strictly thermal spectrum, showing that such a probability of emission is really associated to two non-strictly thermal distributions for boson and fermions. Here we finalize the model by finding the correct value of the pre-factor of the Parikh and Wilczek probability of emission. In fact, that expression has the "of order" sign instead of the equality. In general, in this kind of leading order tunnelling calculations, the exponent arises indeed from the classical action and the pre-factor is an order Planck constant correction. But in the case of emissions of Hawking quanta, the variation of the Bekenstein-Hawking entropy is order 1 for an emitted particle having energy of order the Hawking temperature. As a consequence, the exponent in the Parikh and Wilczek probability of emission is order unity and one asks what is the real significance of that scaling if the pre-factor is unknown. Here we solve the problem assuming the unitarity of the black hole (BH) quantum evaporation and considering the natural correspondence between Hawking radiation and quasi-normal modes (QNMs) of excited BHs, in a "Bohr-like model" that we recently discussed in a series of papers. In that papers, QNMs are interpreted as natural BH quantum levels (the "electron states" in the "Bohr-like model"). Here we find the intriguing result that, although in general it is well approximated by 1, the pre-factor of the Parikh and Wilczek probability of emission depends on the BH quantum level n. We also write down an elegant expression of the probability of emission in terms of the BH quantum levels.
We show that dark energy and dark matter can be described simultaneously by ordinary Einstein gravity interacting with a single scalar field provided the scalar field Lagrangian couples in a symmetric fashion to two different spacetime volume-forms (covariant integration measure densities) on the spacetime manifold - one standard Riemannian given by the square-root of the determinant of the pertinent Riemannian metric and another non-Riemannian volume-form independent of the Riemannian metric, defined in terms of an auxiliary antisymmetric tensor gauge field of maximal rank. Integration of the equations of motion of the latter auxiliary gauge field produce an a priori arbitrary integration constant that plays the role of a dynamically generated cosmological constant or dark energy. Moreover, the above modified scalar field action turns out to possess a hidden Noether symmetry whose associated conserved current describes a pressureless "dust" fluid which we can identify with the dark matter completely decoupled from the dark energy. The form of both the dark energy and dark matter that results from above class of models is insensitive to the specific form of the scalar field Lagrangian. By adding appropriate perturbation, which breaks the above hidden symmetry and along with this it couples dark matter and dark energy, we also suggest a way to obtain growing dark energy in the present universe's epoch without evolution pathologies.
Gravitational waves from neutron star binary inspirals contain information about the equation of state of supranuclear matter. In the absence of definitive experimental evidence that determines the correct equation of state, a number of diverse models that give the pressure in a neutron star as function of its density have been proposed. These models differ not only in the approximations and techniques they use to solve the many-body Schr\"odinger equation, but also in the neutron star composition they assume. We study whether gravitational wave observations of neutron star binaries in quasicircular inspirals will allow us to distinguish between equations of state of differing internal composition, thereby providing important information about the properties of extremely high density matter. We carry out a Bayesian model selection analysis, and find that second generation gravitational wave detectors can heavily constrain equations of state that contain only quark matter, but hybrid stars containing both normal and quark matter are harder to distinguish from normal matter stars. A gravitational wave detection with a signal-to-noise ratio of 30 and masses around $1.4M_{\odot}$ could either detect or rule out strange quark stars with a 20 to 1 confidence. The presence of kaon condensates or hyperons in neutron star inner cores cannot be easily confirmed. For example, for the equations of state studied in this paper, even a gravitational wave signal with a signal-to-noise ratio as high as 60 would not allow us to claim a detection of kaon condensates or hyperons with confidence greater than 5 to 1. On the other hand, if kaon condensates and hyperons do not form in neutron stars, a gravitational wave signal with similar signal-to-noise ratio would be able to constrain their existence with an 11 to 1 confidence for high-mass systems.
Based on the q-kinetic theory in nonextensive statistics, we study the nonextensive q-parameter for the astrophysical systems in an external rotating field. We exactly obtained the equation of the q-parameter for the rotating self-gravitating system. We show that the q-parameter is not only related to the temperature gradient and the gravitational acceleration of the system, but also depends on the inertial centrifugal acceleration and the angular velocity of the rotation, and so the rotation introduces the nonextensivity. This equation of the q-parameter is also applicable to the rotating space plasmas, where an exact expression is presented. We take the Sun, Jupiter and Saturn as examples to illustrate the nonextensive effect introduced by the rotation.
The second generation of ground-based gravitational-wave detectors will begin taking data in September 2015. Sensitive and computationally-efficient data analysis methods will be required to maximize what we learn from their observations. We describe improvements made to the offline analysis pipeline searching for gravitational waves from stellar-mass compact binary coalescences, and assess how these improvements affect search sensitivity. Starting with the two-stage ihope pipeline used in S5, S6 and VSR1-3 and using two weeks of S6/VSR3 data as test periods, we first demonstrate a pipeline with a simpler workflow. This single-stage pipeline performs matched filtering and coincidence testing only once. This simplification allows us to reach much lower false-alarm rates for loud candidate events. We then describe an optimized chi-squared test which minimizes computational cost. Next, we compare methods of generating template banks, demonstrating that a fixed bank may be used for extended stretches of time. Fixing the bank reduces the cost and complexity, compared to the previous method of regenerating a template bank every 2048 s of analyzed data. Creating a fixed bank shared by all detectors also allows us to apply a more stringent coincidence test, whose performance we quantify. With these improvements, we find a 10% increase in sensitive volume with a negligible change in computational cost.
The sensitivity of all-sky searches for gravitational-wave pulsars is primarily limited by the finite availability of computing resources. Semicoherent searches are a widely-used method of maximizing sensitivity to gravitational-wave pulsars at fixed computing cost: the data from a gravitational-wave detector are partitioned into a number of segments, each segment is coherently analyzed, and the analysis results from each segment are summed together. The generation of template banks for the coherent analysis of each segment, and for the summation, requires knowledge of the metrics associated with the coherent and semicoherent parameter spaces respectively. We present a useful approximation to the semicoherent parameter-space metric, analogous to that presented in Wette and Prix [Phys. Rev. D 88, 123005 (2013)] for the coherent metric. The new semicoherent metric is compared to previous work in Pletsch [Phys. Rev. D 82, 042002 (2010)], and Brady and Creighton [Phys. Rev. D 61, 082001 (2000)]. We find that semicoherent all-sky searches require orders of magnitude more templates than previously predicted.
Links to: arXiv, form interface, find, astro-ph, recent, 1508, contact, help (Access key information)
The Dark Energy Survey (DES) is a 5000 sq. degree survey in the southern hemisphere, which is rapidly reducing the existing north-south asymmetry in the census of MW satellites and other stellar substructure. We use the first-year DES data down to previously unprobed photometric depths to search for stellar systems in the Galactic halo, therefore complementing the previous analysis of the same data carried out by our group earlier this year. Our search is based on a matched filter algorithm that produces stellar density maps consistent with stellar population models of various ages, metallicities, and distances over the survey area. The most conspicuous density peaks in these maps have been identified automatically and ranked according to their significance and recurrence for different input models. We report the discovery of one additional stellar system besides those previously found by several authors using the same first-year DES data. The object is compact, and consistent with being dominated by an old and metal-poor population. DES J0034-4902 is found at high significance and appears in the DES images as a compact concentration of faint blue point sources at ~ 87 {kpc}. Its half-light radius of r_h = 9.88 +/- 4.31 {pc} and total luminosity of M_V ~ -3.05_{-0.42}^{+0.69} are consistent with it being a low mass halo cluster. It is also found to have a very elongated shape. In addition, our deeper probe of DES 1st year data confirms the recently reported satellite galaxy candidate Horologium II as a significant stellar overdensity. We also infer its structural properties and compare them to those reported in the literature.
Models for tidal disruption events (TDEs) in which a supermassive black hole disrupts a star commonly assume that the highly eccentric streams of bound stellar debris promptly form a circular accretion disk at the pericenter scale. However, the bolometric peak luminosity of most TDE candidates, $\sim10^{44}\,\rm{erg\,s^{-1}}$, implies that we observe only $\sim1\%$ of the energy expected from accretion. Moreover, recent numerical simulations (Shiokawa et al. 2015) have shown that dissipation via hydrodynamical shocks is insufficient to circularize debris orbits on the pericenter scale, and the debris flow retains its initial semi-major axis scale throughout the first $\sim10$ orbits of the event. Motivated by these numerical results, Piran et al. (2015) suggested that the observed optical TDE emission is powered by shocks at the apocenter between freshly infalling material and earlier-arriving matter. This model explains the small radiated energy, the low temperature, and the large radius implied by the observations as well as the $t^{-5/3}$ light curve. However the question of the system's low efficiency remains unanswered. We suggest that the high orbital energy and low angular momentum of the flow's highly eccentric orbits make it possible for magnetic stresses to reduce the matter's already small angular momentum to the point at which it can pass within the ISCO before these stresses remove enough energy to circularize the orbit. As a result, the typical efficiency is only $\sim3\%$ of a standard accretion disk's efficiency. Thus, the intrinsically high eccentricity of the tidal debris naturally explains why most TDE candidates are fainter than expected.
We investigate the distribution of particle acceleration sites during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type current layers in pair plasma with no guide magnetic field, negligible radiative losses, no initial perturbation, and using periodic boundary conditions. We find that the plasmoids develop a robust internal structure, with colder dense cores and hotter outer shells, that is recovered after each plasmoid merger on a dynamical time scale. We use spacetime diagrams of the reconnection layers to probe the evolution of plasmoids, and in this context we investigate the individual particle histories for a representative sample of energetic electrons. We distinguish three classes of particle acceleration sites associated with (1) magnetic X-points, (2) regions between merging plasmoids, and (3) the trailing edges of accelerating plasmoids. We evaluate the contribution of each class of acceleration sites to the final energy distribution of energetic electrons -- magnetic X-points dominate at moderate energies, and the regions between merging plasmoids dominate at higher energies. We also identify the dominant acceleration scenarios, in order of decreasing importance -- (1) single acceleration between merging plasmoids, (2) single acceleration at a magnetic X-point, and (3) acceleration at a magnetic X-point followed by acceleration in a plasmoid. Particle acceleration is absent only in the vicinity of stationary plasmoids, and it can hardly be associated with magnetic mirrors due to the absence of plasmoid contraction after the initial stage of the simulation.
We show that the new precise measurements of Cosmic Microwave Background (CMB) temperature and polarization anisotropies made by the Planck satellite significantly improves previous constraints on the cosmic gravitational waves background (CGWB) at frequencies $f>10^{-15}$ Hz. On scales smaller than the comoving horizon at the time of decoupling, primordial gravitational waves contribute to the total radiation content of the Universe. Considering adiabatic perturbations, CGWB affects temperature and polarization CMB power spectra and matter power spectrum in a manner identical to relativistic particles. Considering the latest Planck results we constrain the CGWB energy density to $\Omega_{\rm gw} h^2 <2.5\times 10^{-6} $ at 95\% c.l. Combining CMB power spectra with lensing, BAO and primordial Deuterium abundance observations, we obtain $\Omega_{\rm gw} h^2 <1.7\times 10^{-6} $ at 95\% c.l., improving previous cosmological bounds by a factor 5 and the recent direct upper limit from the LIGO and VIRGO experiments by 50\%. A combined analysis of future satellite missions as CORE and EUCLID could improve current bound by a factor $\sim20$.
We report on a multiwavelength observational campaign of the black hole X-ray binary Swift J1753.5-0127 that consists of an ESO/X-shooter spectrum supported by contemporaneous Swift/XRT+UVOT and ATCA data. ISM absorption lines in the X-shooter spectrum allows us to determine E(B-V)=0.45+/-0.02 along the line-of-sight to the source. We also report detection of emission signatures of He II at 4686 angstrom, H alpha, and, for the first time, H I at 10906 angstrom and Paschen Beta. The double-peaked morphology of these four lines is typical of the chromosphere of a rotating accretion disk. Nonetheless, the paucity of disk features points towards a low level of irradiation in the system. This is confirmed through spectral energy distribution modeling and we find that the UVOT+X-shooter continuum mostly stems from the thermal emission of a viscous disk. We speculate that the absence of reprocessing is due to the compactness of an illumination-induced envelope that fails to reflect enough incoming hard X-ray photons back to the outer regions. The disk also marginally contributes to the Compton-dominated X-ray emission and is strongly truncated, with an inner radius about a thousand times larger than the black hole's gravitational radius. A near-infrared excess is present, and we associate it with synchrotron radiation from a compact jet. However, the measured X-ray flux is significantly higher than what can be explained by the optically thin synchrotron jet component. We discuss these findings in the framework of the radio quiet versus X-ray bright hypothesis, favoring the presence of a residual disk, predicted by evaporation models, that contributes to the X-ray emission without enhancing the radio flux.
We present results from a Hubble Space Telescope (HST) program characterizing the atmospheres of the outer two planets, in the HR8799 system. The images were taken over 15 orbits in three near-infrared medium-band filters - F098M, F127M and F139M - using the Wide Field Camera 3. One of the three filters is sensitive to water absorption band inaccessible from ground-based observations, providing a unique probe of the thermal emission from the atmospheres of these young giant planets. The observations were taken at 30 different spacecraft rolls to enable angular differential imaging, and the full data set was analyzed with the Karhunen-Loeve Image Projection (KLIP) routine, an advanced image processing algorithm adapted to work with HST data. To achieve the required high contrast at sub arcsecond resolution, we utilized the pointing accuracy of HST in combination with an improved pipeline designed to combine the dithered, angular differential imaging data with an algorithm designed to both improve the image resolution and accurately measure the photometry. The results include F127M ($J$) detections of the outer planets, HR8799 b and c and the first detection of HR8799 b in the water-band (F139M) filter. The F127M photometry for HR8799 c agrees well with fitted atmospheric models resolving a long standing difficulty to model the near-IR flux for the planet consistently
We have analyzed new and archival time series spectra taken six years apart during transits of the hot Jupiter WASP-33 b, and spectroscopically resolved the line profile perturbation caused by the Rossiter-McLaughlin effect. The motion of this line profile perturbation is determined by the path of the planet across the stellar disk, which we show to have changed between the two epochs due to nodal precession of the planetary orbit. We measured rates of change of the impact parameter and the sky-projected spin-orbit misalignment of $db/dt=-0.0228_{-0.0018}^{+0.0050}$ yr$^{-1}$ and $d\lambda/dt=-0.487_{-0.076}^{+0.089}$ $^{\circ}$ yr$^{-1}$, respectively, corresponding to a rate of nodal precession of $d\Omega/dt=0.117_{-0.029}^{+0.012}$ $^{\circ}$ yr$^{-1}$. This is only the second measurement of nodal precession for a confirmed exoplanet transiting a single star. Finally, we used the rate of precession to set limits on the stellar gravitational quadrupole moment of $0.0017\leq J_2\leq0.011$.
In this paper, we model the observable signatures of tidal disruptions of white dwarf (WD) stars by massive black holes (MBHs) of moderate mass, $\approx 10^3 - 10^5 M_\odot$. When the WD passes deep enough within the MBH's tidal field, these signatures include thermonuclear transients from burning during maximum compression. We combine a hydrodynamic simulation that includes nuclear burning of the disruption of a $0.6 M_\odot$ C/O WD with a Monte Carlo radiative transfer calculation to synthesize the properties of a representative transient. The transient's emission emerges in the optical, with lightcurves and spectra reminiscent of type I SNe. The properties are strongly viewing-angle dependent, and key spectral signatures are $\approx 10,000$ km s$^{-1}$ Doppler shifts due to the orbital motion of the unbound ejecta. Disruptions of He WDs likely produce large quantities of intermediate-mass elements, offering a possible production mechanism for Ca-rich transients. Accompanying multiwavelength transients are fueled by accretion and arise from the nascent accretion disk and relativistic jet. If MBHs of moderate mass exist with number densities similar to those of supermassive MBHs, both high energy wide-field monitors and upcoming optical surveys should detect tens to hundreds of WD tidal disruptions per year. The current best strategy for their detection may therefore be deep optical follow up of high-energy transients of unusually-long duration. The detection rate or the non-detection of these transients by current and upcoming surveys can thus be used to place meaningful constraints on the extrapolation of the MBH mass function to moderate masses.
Two explanations have been put forward to explain the observed conformity between the colours and specific star formation rates (SFR/$M_*$) of galaxies on large scales: 1) the formation times of their surrounding dark matter halos are correlated (commonly referred to as "assembly bias"), 2) gas is heated over large scales at early times, leading to coherent modulation of cooling and star formation between well-separated galaxies (commonly referred to as "pre-heating") . To distinguish between the pre-heating and assembly bias scenarios, we search for relics of energetic feedback events in the neighbourhood of central galaxies with different specific star formation rates. We find a significant excess of very high mass ($\log M_* > 11.3$) galaxies out to a distance of 2.5 Mpc around low SFR/$M_*$ central galaxies compared to control samples of higher SFR/$M_*$ central galaxies with the same stellar mass and redshift. We also find that very massive galaxies in the neighbourhood of low SFR/$M_*$ galaxies have much higher probability of hosting radio loud active galactic nuclei. The radio-loud AGN fraction in neighbours with $\log M_* > 11.3$ is four times higher around passive, non star-forming centrals at projected distances of 1 Mpc and two times higher at projected distances of 4 Mpc. Finally, we carry out an investigation of conformity effects in the recently publicly-released Illustris cosmological hydrodynamical simulation, which includes energetic input both from quasars and from radio mode accretion onto black holes. We do not find conformity effects of comparable amplitude on large scales in the simulations and we propose that gas needs to be pushed out of dark matter halos more efficiently at high redshifts.
The merger of two white dwarfs (WDs) has for many years not been considered as the favoured model for the progenitor system of type Ia supernovae (SNe Ia). But recent years have seen a change of opinion as a number of studies, both observational and theoretical, have concluded that they should contribute significantly to the observed type Ia supernova rate. In this paper, we study the ignition and propagation of detonation through post-merger remnants and we follow the resulting nucleosynthesis up to the point where a homologous expansion is reached. In our study we cover the entire range of WD masses and compositions. For the emergence of a detonation we study several setups, guided by both merger remnants from our own simulations and by results taken from the literature. We carefully compare the nucleosynthetic yields of successful explosions with SN Ia observations. Only three of our models are consistent with all the imposed constraints and potentially lead to a standard type Ia event. The first one, a $0.45\,M_\odot\ {\rm helium} + 0.9 \,M_\odot$ carbon-oxygen WD system produces a sub-luminous, SN 1991bg-like event while the other two, a $0.45\,M_\odot\ {\rm helium} + 1.1\,M_\odot$ oxygen-neon WD system and a $1.05+1.05\,M_\odot$ system with two carbon-oxygen WDs, are good candidates for common SNe Ia.
If advanced extraterrestrial civilizations choose to construct vast numbers of Dyson spheres to harvest radiation energy, this could affect the characteristics of their host galaxies. Potential signatures of such astroengineering projects include reduced optical luminosity, boosted infrared luminosity and morphological anomalies. Here, we apply a technique pioneered by Annis (1999) to search for Kardashev type III civilizations in disk galaxies, based on the predicted offset of these galaxies from the optical Tully-Fisher relation. By analyzing a sample of 1359 disk galaxies, we are able to set a conservative upper limit at 3% on the fraction of local disks subject to Dysonian astroengineering on galaxy-wide scales. However, the available data suggests that a small subset of disk galaxies actually may be underluminous with respect to the Tully-Fisher relation in the way expected for Kardashev type III objects. Based on the optical morphologies and infrared-to-optical luminosity ratios of such galaxies in our sample, we conclude that none of them stand out as strong Kardashev type III candidates and that their inferred properties likely have mundane explanations. This allows us to set a tentative upper limit at 0.3% on the fraction of Karashev type III disk galaxies in the local Universe.
Using new data from the K2 mission, we show that WASP-47, a previously known hot Jupiter host, also hosts two additional transiting planets: a Neptune-sized outer planet and a super-Earth inner companion. We measure planetary properties from the K2 light curve and detect transit timing variations, confirming the planetary nature of the outer planet. We performed a large number of numerical simulations to study the dynamical stability of the system and to find the theoretically expected transit timing variations (TTVs). The theoretically predicted TTVs are in good agreement with those observed, and we use the TTVs to determine the masses of two planets, and place a limit on the third. The WASP-47 planetary system is important because companion planets can both be inferred by TTVs and are also detected directly through transit observations. The depth of the hot Jupiter's transits make ground-based TTV measurements possible, and the brightness of the host star makes it amenable for precise radial velocity measurements. The system serves as a Rosetta Stone for understanding TTVs as a planet detection technique.
We describe the layout and unique features of the focal plane system for MIRI. We begin with the detector array and its readout integrated circuit (combining the amplifier unit cells and the multiplexer), the electronics, and the steps by which the data collection is controlled and the output signals are digitized and delivered to the JWST spacecraft electronics system. We then discuss the operation of this MIRI data system, including detector readout patterns, operation of subarrays, and data formats. Finally, we summarize the performance of the system, including remaining anomalies that need to be corrected in the data pipeline.
The Quadrantid meteor shower is among the strongest annual meteor showers,
and has drawn the attention of scientists for several decades. The stream is
unusual, among others, for several reasons: its very short duration around
maximum activity (~12 - 14 hours) as detected by visual, photographic and radar
observations, its recent onset (around 1835 AD) and because it had been the
only major stream without an obvious parent body until 2003. Ever since, there
have been debates as to the age of the stream and the nature of its proposed
parent body, asteroid 2003 EH1.
In this work, we present results on the most probable age and formation
mechanism of the narrow portion of the Quadrantid meteoroid stream. For the
first time we use data on eight high precision photographic Quadrantids,
equivalent to gram - kilogram size, to constrain the most likely age of the
core of the stream. Out of eight high-precision photographic Quadrantids, five
pertain directly to the narrow portion of the stream. In addition, we also use
data on five high-precision radar Quadrantids, observed within the peak of the
shower.
We performed backwards numerical integrations of the equations of motion of a
large number of 'clones' of both, the eight high-precision photographic and
five radar Quadrantid meteors, along with the proposed parent body, 2003 EH1.
According to our results, from the backward integrations, the most likely age
of the narrow structure of the Quadrantids is between 200 - 300 years. These
presumed ejection epochs, corresponding to 1700 - 1800 AD, are then used for
forward integrations of large numbers of hypothetical meteoroids, ejected from
the parent 2003 EH$_1$, until the present epoch. The aim is to constrain
whether the core of the Quadrantid meteoroid stream is consistent with a
previously proposed relatively young age (~ 200 years).}
We present an estimate of the performance that will be achieved during on orbit operations of the JWST Mid Infrared Instrument, MIRI. The efficiency of the main imager and spectrometer systems in detecting photons from an astronomical target are presented, based on measurements at sub-system and instrument level testing, with the end-to-end transmission budget discussed in some detail. The brightest target fluxes that can be measured without saturating the detectors are provided. The sensitivity for long duration observations of faint sources is presented in terms of the target flux required to achieve a signal to noise ratio of 10 after a 10,000 second observation. The algorithms used in the sensitivity model are presented, including the understanding gained during testing of the MIRI Flight Model and flight-like detectors.
The radiative efficiency of super-Eddington accreting black holes (BHs) is explored for magnetically-arrested disks (MADs), where magnetic flux builds-up to saturation near the BH. Our three-dimensional general relativistic radiation magnetohydrodynamic (GRRMHD) simulation of a spinning BH (spin $a/M=0.8$) accreting at $\sim 50$ times Eddington shows a total efficiency $\sim 50\%$ when time-averaged and total efficiency $\gtrsim 100\%$ in moments. Magnetic compression by the magnetic flux near the rotating BH leads to a thin disk, whose radiation escapes via advection by a magnetized wind and via transport through a low-density channel created by a Blandford-Znajek (BZ) jet. The BZ efficiency is sub-optimal due to inertial loading of field lines by optically thick radiation, leading to BZ efficiency $\sim 40\%$ on the horizon and BZ efficiency $\sim 5\%$ by $r\sim 400r_g$ (gravitational radii) via absorption by the wind. Importantly, radiation escapes at $r\sim 400r_g$ with efficiency $\eta\approx 15\%$ (luminosity $L\sim 50L_{\rm Edd}$), similar to $\eta\approx 12\%$ for a Novikov-Thorne thin disk and beyond $\eta\lesssim 1\%$ seen in prior GRRMHD simulations or slim disk theory. Our simulations show how BH spin, magnetic field, and jet mass-loading affect the radiative and jet efficiencies of super-Eddington accretion.
Context. Solar eruptions and high flare activity often accompany the rapid
rotation of sunspots. The study of sunspot rotation and the mechanisms driving
this motion are therefore key to our understanding of how the solar atmosphere
attains the conditions necessary for large energy release.
Aims. We aim to demonstrate and investigate the rotation of sunspots in a 3D
numerical experiment of the emergence of a magnetic flux tube as it rises
through the solar interior and emerges into the atmosphere. Furthermore, we
seek to show that the sub-photospheric twist stored in the interior is injected
into the solar atmosphere by means of a definitive rotation of the sunspots.
Methods. A numerical experiment is performed to solve the 3D resistive
magnetohydrodynamic (MHD) equations using a Lagrangian-Remap code. We track the
emergence of a toroidal flux tube as it rises through the solar interior and
emerges into the atmosphere investigating various quantities related to both
the magnetic field and plasma.
Results. Through detailed analysis of the numerical experiment, we find clear
evidence that the photospheric footprints or sunspots of the flux tube undergo
a rotation. Significant vertical vortical motions are found to develop within
the two polarity sources after the field emerges. These rotational motions are
found to leave the interior portion of the field untwisted and twist up the
atmospheric portion of the field. This is shown by our analysis of the relative
magnetic helicity as a significant portion of the interior helicity is
transported to the atmosphere. In addition, there is a substantial transport of
magnetic energy to the atmosphere. Rotation angles are also calculated by
tracing selected fieldlines; the fieldlines threading through the sunspot are
found to rotate through angles of up to 353 degrees over the course of the
experiment.
We describe the operations concept and data reduction plan for the Mid- Infrared Instrument (MIRI) for the James Webb Space Telescope (JWST). The overall JWST operations concept is to use Observation Templates (OTs) to provide a straightforward and intuitive way for users to specify observations. MIRI has four OTs that correspond to the four observing modes: 1.) Imaging, 2.) Coronagraphy, 3.) Low Resolution Spectroscopy, and 4.) Medium Resolution Spectroscopy. We outline the user choices and expansion of these choices into detailed instrument operations. The data reduction plans for MIRI are split into three stages, where the specificity of the reduction steps to the observation type increases with stage. The reduction starts with integration ramps: stage 1 yields uncalibrated slope images; stage 2 calibrates the slope images; and then stage 3 combines multiple calibrated slope images into high level data products (e.g. mosaics, spectral cubes, and extracted source information). Finally, we give examples of the data and data products that will be derived from each of the four different OTs.
We study the relation of active galactic nuclei (AGNs) to star formation in their host galaxies. Our sample includes 205 Type-1 and 85 Type-2 AGNs, 162 detected with Herschel, from fields surrounding 30 galaxy clusters in the Local Cluster Substructure Survey (LoCuSS). The sample is identified by optical line widths and ratios after selection to be brighter than 1 mJy at 24 microns. We show that Type-2 AGN [OIII]5007 line fluxes at high z can be contaminated by their host galaxies with typical spectrograph entrance apertures (but our sample is not compromised in this way). We use spectral energy distribution (SED) templates to decompose the galaxy SEDs and estimate star formation rates, AGN luminosities, and host galaxy stellar masses (described in an accompanying paper). The AGNs arise from massive black holes (~ 3 X 10^8 Msun) accreting at ~ 10% of the Eddington rate and residing in galaxies with stellar mass > 3 X 10^{10} Msun; those detected with Herschel have IR luminosity from star formation in the range of 10^{10} -- 10^{12} Lsun. We find that: 1.) the specific star formation rates in the host galaxies are generally consistent with those of normal star-forming (main sequence) galaxies; 2.) there is a strong correlation between the luminosities from star formation and the AGN; and 3.) however, the correlation may not result from a causal connection, but could arise because the black hole mass (and hence AGN Eddington luminosity) and star formation are both correlated with the galaxy mass.
We present a sample of 290 24-micron-selected active galactic nuclei (AGNs) mostly at z ~ 0.3 -- 2.5, within 5.2 square degrees distributed as 25' X 25' fields around each of 30 galaxy clusters in the Local Cluster Substructure Survey (LoCuSS). The sample is nearly complete to 1 mJy at 24 microns, and has a rich multi-wavelength set of ancillary data; 162 are detected by Herschel. We use spectral templates for AGNs, stellar populations, and infrared emission by star forming galaxies to decompose the spectral energy distributions (SEDs) of these AGNs and their host galaxies, and estimate their star formation rates (SFRs), AGN luminosities, and host galaxy stellar masses. The set of templates is relatively simple: a standard Type-1 quasar template; another for the photospheric output of the stellar population; and a far infrared star-forming template. For the Type-2 AGN SEDs, we substitute templates including internal obscuration, and some Type-1 objects require a warm component (T > 50 K). The individually Herschel- detected Type-1 AGNs and a subset of 17 Type-2 ones typically have luminosities > 10^{45} ergs/s, and supermassive black holes of ~ 3 X 10^8 Msun emitting at ~ 10% of the Eddington rate. We find them in about twice the numbers of AGN identified in SDSS data in the same fields, i.e., they represent typical high luminosity AGN, not an infrared-selected minority. These AGNs and their host galaxies are studied further in an accompanying paper.
The spatial curvature ($K$ or $\Omega_K$) is one of the most fundamental parameters of isotropic and homogeneous universe and has a close link to the physics of early universe. Combining the radial and angular diameter distances measured via the baryon acoustic oscillation (BAO) experiments allows us to unambiguously constrain the curvature. The method is primarily based on the metric theory, but not much on the theory of structure formation other than the existence of BAO scale and is free of any model of dark energy. In this paper, we estimate a best-achievable accuracy of constraining the curvature with the BAO experiments. We show that an all-sky, cosmic-variance-limited galaxy survey covering the universe up to $z>4$ enables a precise determination of the curvature to an accuracy of $\sigma(\Omega_K)\simeq 10^{-3}$. When we assume a model of dark energy, either the cosmological constraint or the $(w_0,w_a)$-model, it can achieve a precision of $\sigma(\Omega_K)\simeq \mbox{a few}\times 10^{-4}$. These forecasts require a high sampling density of galaxies, and are degraded by up to a factor of a few for a survey with a finite number density of galaxies.
Using a grid of $\sim 2$ million elements ($\Delta z = 0.005$) adapted from COSMOS photometric redshift (photo-z) searches, we investigate the general properties of template-based photo-z likelihood surfaces. We find these surfaces are filled with numerous local minima and large degeneracies that generally confound rapid but "greedy" optimization schemes, even with additional stochastic sampling methods. In order to robustly and efficiently explore these surfaces, we develop BAD-Z [Brisk Annealing-Driven Redshifts (Z)], which combines ensemble Markov Chain Monte Carlo (MCMC) sampling with simulated annealing to sample arbitrarily large, pre-generated grids in approximately constant time. Using a mock catalog of 384,662 objects, we show BAD-Z samples $\sim 40$ times more efficiently compared to a brute-force counterpart while maintaining similar levels of accuracy. Our results represent first steps toward designing template-fitting photo-z approaches limited mainly by memory constraints rather than computation time.
In this article, we describe the MIRI Imager module (MIRIM), which provides broad-band imaging in the 5 - 27 microns wavelength range for the James Webb Space Telescope. The imager has a 0"11 pixel scale and a total unobstructed view of 74"x113". The remainder of its nominal 113"x113" field is occupied by the coronagraphs and the low resolution spectrometer. We present the instrument optical and mechanical design. We show that the test data, as measured during the test campaigns undertaken at CEA-Saclay, at the Rutherford Appleton Laboratory, and at the NASA Goddard Space Flight Center, indicate that the instrument complies with its design requirements and goals. We also discuss the operational requirements (multiple dithers and exposures) needed for optimal scientific utilization of the MIRIM.
Sterile neutrinos in the electronvolt mass range are hinted at by a number of terrestrial neutrino experiments. However, such neutrinos are highly incompatible with data from the Cosmic Microwave Background and large scale structure. This paper discusses how charging sterile neutrinos under a new pseudoscalar interaction can reconcile eV sterile neutrinos with terrestrial neutrino data. We show that this model can reconcile eV sterile neutrinos in cosmology, providing a fit to all available data which is way better than the standard $\Lambda$CDM model with one additional fully thermalized sterile neutrino. In particular it also prefers a value of the Hubble parameter much closer to the locally measured value.
We compare the performance of two very different parallel gravitational $N$-body codes for astrophysical simulations on large GPU clusters, both pioneer in their own fields as well as in certain mutual scales - NBODY6++ and Bonsai. We carry out the benchmark of the two codes by analyzing their performance, accuracy and efficiency through the modeling of structure decomposition and timing measurements. We find that both codes are heavily optimized to leverage the computational potential of GPUs as their performance has approached half of the maximum single precision performance of the underlying GPU cards. With such performance we predict that a speed-up of $200-300$ can be achieved when up to 1k processors and GPUs are employed simultaneously. We discuss the quantitative information about comparisons of two codes, finding that in the same cases Bonsai adopts larger time steps as well as relative energy errors than NBODY6++, typically ranging from $10-50$ times larger, depending on the chosen parameters of the codes. While the two codes are built for different astrophysical applications, in specified conditions they may overlap in performance at certain physical scale, and thus allowing the user to choose from either one with finetuned parameters accordingly.
We investigate the variability behaviour of the broad Hb emission-line to
driving continuum variations in the best-studied AGN NGC 5548. For a particular
choice of BLR geometry, Hb surface emissivity based on photoionization models,
and using a scaled version of the 13 yr optical continuum light curve as a
proxy for the driving ionizing continuum, we explore several key factors that
determine the broad emission line luminosity L, characteristic size R(RW), and
variability amplitude (i.e., responsivity) eta, as well as the interplay
between them.
For fixed boundary models which extend as far as the hot-dust the predicted
delays for Hb are on average too long. However, the predicted variability
amplitude of Hb provides a remarkably good match to observations except during
low continuum states. We suggest that the continuum flux variations which drive
the redistribution in Hb surface emissivity F(r) do not on their own lead to
large enough changes in R(RW) or eta(eff). We thus investigate dust-bounded
BLRs for which the location of the effective outer boundary is modulated by the
continuum level and the dust-sublimation and dust-condensation timescales. We
find that in order to match the observed variability amplitude of broad Hb in
NGC 5548 a rather static outer boundary is preferred.
Intriguingly, we show that the most effective way of reducing the Hb delay,
while preserving its responsivity and equivalent width, is to invoke a smaller
value in the incident ionizing photon flux Phi(H) for a given ionizing
source--cloud radial distance r, than is normally inferred from the observed UV
continuum flux and typical models of the continuum SED.
Magnetic fields are fundamental to the dynamics of both accretion disks and the jets that they often drive. We review the basic physics of these phenomena, the past and current efforts to model them numerically with an emphasis on the jet-disk connection, and the observational constraints on the role of magnetic fields in the jets of active galaxies on all scales.
We study the $\rm{\gamma}$-ray luminosity and beaming effect for Fermi blazars. Our results are as follows. (i) There are significant correlations between $\rm{\gamma}$-ray luminosity and radio core luminosity, and between $\rm{\gamma}$-ray luminosity and $\rm{R_{v}}$, which suggests that the $\rm{\gamma}$-ray emissions have strong beaming effect. (ii) Using the $\rm{L_{ext}/M_{abs}}$ as an indicator of environment effects, we find that there have no significant correlation between $\rm{\gamma}$-ray luminosity and $\rm{L_{ext}/M_{abs}}$ for all sources when remove the effect of redshift. FSRQs considered alone also do not show a significant correlation, while BL Lacs still show a significant correlation when remove the effect of redshift. These results suggest that the $\rm{\gamma}$-ray emission may be affected by environment on the kiloparsec-scale for BL Lacs.
We examine the properties of the outflowing matter from an advective accretion disc around a spinning black hole. During accretion, rotating matter experiences centrifugal pressure supported shock transition that effectively produces a virtual barrier around the black hole in the form of post-shock corona (hereafter, PSC). Due to shock compression, PSC becomes hot and dense that eventually deflects a part of the inflowing matter as bipolar outflows because of the presence of extra thermal gradient force. In our approach, we study the outflow properties in terms of the inflow parameters, namely specific energy (${\mathcal E}$) and specific angular momentum ($\lambda$) considering the realistic outflow geometry around the rotating black holes. We find that spin of the black hole ($a_k$) plays an important role in deciding the outflow rate $R_{\dot m}$ (ratio of mass flux of outflow and inflow), in particular, $R_{\dot m}$ is directly correlated with $a_k$ for the same set of inflow parameters. It is found that a large range of the inflow parameters allows global accretion-ejection solutions and the effective area of the parameter space (${\mathcal E}$, $\lambda$) with and without outflow decreases with black hole spin ($a_k$). We compute the maximum outflow rate ($R^{max}_{\dot m}$) as function of black hole spin ($a_k$) and observe that $R^{max}_{\dot m}$ weakly depends on $a_k$ that lies in the range $\sim 10\%-18\%$ of the inflow rate for the adiabatic index $(\gamma)$ with $1.5 \ge \gamma \ge 4/3$. We present the observational implication of our approach while studying the steady/persistent Jet activities based on the accretion states of black holes. We discuss that our formalism seems to have the potential to explain the observed Jet kinetic power for several Galactic Black Hole sources (GBHs) and Active Galactic Nuclei (AGNs).
High resolution (R$\sim$22,500) spectra for 400 red clump giants, in four fields within $\rm -4.8^{\circ} \lesssim b \lesssim -3.4^{\circ}$ and $\rm -10^{\circ} \lesssim l \lesssim +10^{\circ}$, were obtained within the GIRAFFE Inner Bulge Survey (GIBS) project. To this sample we added another $\sim$ 400 stars in Baade's Window, observed with the identical instrumental configuration. We constructed the metallicity distributions for the entire sample, as well as for each field individually, in order to investigate the presence of gradients or field-to-field variations in the shape of the distributions. The metallicity distributions in the five fields are consistent with being drawn from a single parent population, indicating the absence of a gradient along the major axis of the Galactic bar. The global metallicity distribution is well fitted by two Gaussians. The metal poor component is rather broad, with a mean at $\rm <[Fe/H]>=-0.31$ dex and $\sigma=0.31$ dex. The metal-rich one is narrower, with mean $\rm <[Fe/H]>=+0.26$ and $\sigma=0.2$ dex. The [Mg/Fe] ratio follows a tight trend with [Fe/H], with enhancement with respect to solar in the metal-poor regime, similar to the one observed for giant stars in the local thick disc. [Ca/Fe] abundances follow a similar trend, but with a considerably larger scatter than [Mg/Fe]. A decrease in [Mg/Fe] is observed at $\rm [Fe/H]=-0.44$ dex. This \textit{knee} is in agreement with our previous bulge study of K-giants along the minor axis, but is 0.1 dex lower in metallicity than the one reported for the Bulge micro lensed dwarf and sub-giant stars. We found no variation in $\alpha$-element abundance distributions between different fields.
We develop a flexible set of action-based distribution functions (DFs) for stellar halos. The DFs have five free parameters, controlling the inner and outer density slope, break radius, flattening and anisotropy respectively. The DFs generate flattened stellar halos with a rapidly varying logarithmic slope in density, as well as a spherically aligned velocity ellipsoid with a long axis that points towards the Galactic centre - all attributes possessed by the stellar halo of the Milky Way. We use our action-based distribution function to model the blue horizontal branch stars extracted from the Sloan Digital Sky Survey as stellar halo tracers in a spherical Galactic potential. As the selection function is hard to model, we fix the density law from earlier studies and solve for the anisotropy and gravitational potential parameters. Our best fit model has a velocity anisotropy that becomes more radially anisotropic on moving outwards. It changes from $\beta \approx 0.4$ at Galactocentric radius of 15 kpc to $\approx 0.7$ at 60 kpc. This is a gentler increase than is typically found in simulations of stellar haloes built from the mutiple accretion of smaller systems. We find the potential corresponds to an almost flat rotation curve with amplitude of $\approx 200$ kms$^{-1}$ at these distances. This implies an enclosed mass of $\approx 4.5 \times 10^{11} M_\odot$ within a spherical shell of radius 50 kpc.
We study the soft excess variability of the Narrow Line Seyfert 1 galaxy IRAS 13224-3809. We considered all five archival XMM-Newton observations, and we applied the "flux-flux plots" method. We found that the flux-flux plots are highly affected by the choice of the light curves' time bin size, most probably due to the fast and large amplitude variations, and the intrinsic non-linear flux-flux relations in this source. We suggest that the smallest bin-size should be used in such cases. We constructed flux-flux plots in 11 energy bands below 1.7 keV, and we considered the 1.7-3 keV band, as being representative of the primary emission. The flux-flux plots are reasonably well fitted by a "power-law plus a constant" model. We detected significant positive constants in 3 out of 5 observations. The best-fit slopes are flatter than unity at energies below $\sim 0.9$ keV, where the soft excess is strongest. This suggests the presence of intrinsic spectral variability. A power-law like primary component, which is variable in flux and spectral slope (as $\Gamma\propto N_{\rm PL}^{0.1}$), and a soft-excess component which varies together with the primary continuum (as $F_{\rm excess}\propto F_{\rm primary}^{0.46}$) can broadly explain the flux-flux plots. In fact, it can result in the presence of positive "constants", even when a stable spectral component doe snot exist. Nevertheless, the possibility of a stable soft-band constant component cannot be ruled out, but its contribution to the observed 0.2-1 keV band flux should be less than $\sim 15$ %. The model constants in the flux-flux plots were consistent with zero in one observation, and negative, at energies below 1 keV, in another one. It is hard to explain these results with any spectral variability scenario, and they may signify the presence of a variable, warm absorber in the source.
Radial velocities measured with the 6-meter telescope are given for 5 faint dwarf galaxies. All of these galaxies are confirmed as very nearby objects. Two of them, KK135 (dIr) and UGC 1703 (dSph/dTr), are local isolated dwarfs, and the three others, UGCA 127sat (dIr), NGC 2683dw1 (dIr), and NGC891dwA (dTr), belong to companions of nearby massive spirals.
We investigated the polarization characteristics of a zebra pattern (ZP) in a type-IV solar radio burst observed with AMATERAS on 2011 June 21 for the purpose of evaluating the generation processes of ZP. Analyzing highly resolved spectral and polarization data revealed the frequency dependence of the degree of circular polarization and the delay between two polarized components for the first time. The degree of circular polarization was 50-70 percent right-handed and it varied little as a function of frequency. Cross-correlation analysis determined that the left-handed circularly polarized component was delayed by 50-70 ms relative to the right-handed component over the entire frequency range of the ZP and this delay increased with the frequency. We examined the obtained polarization characteristics by using pre-existing ZP models and concluded that the ZP was generated by the double plasma resonance process. Our results suggest that the ZP emission was originally generated in a completely polarized state in the O-mode and was partly converted into the X-mode near the source. Subsequently, the difference between the group velocities of the O-mode and X-mode caused the temporal delay.
Wright et al. 2014 have embarked on a search for advanced Karadashev Type III civilisations via the compilation of a sample of sources with extreme mid-IR emission and colours. In this scenario, the mid-IR emission is then primarily associated with waste heat energy by-products. I apply the Mid-IR radio correlation to this $\hat{G}$ sample (Griffith et al. 2015). I demonstrate that the mid-IR and radio luminosities are correlated for the sample with $q_{22}=1.35\pm0.42 $. By comparison, the First Look Survey (FLS) has $q_{22}=0.87\pm0.27$. The fact that the G-HAT sample largely follows the Mid-IR radio correlation, strongly suggests the vast majority of these sources are associated with galaxies in which natural astrophysical processes are dominant. This simple application of the mid-IR radio correlation can substantially reduce the number of false positives in the $\hat{G}$ catalogue, since galaxies occupied by advanced Kardashev Type III civilisations would be expected to exhibit very high values of $q$. Indeed I identify 9 outliers in the sample with $q_{22} > 2$ of which at least 3 have properties that are relatively well explained via standard astrophysical interpretations e.g. dust emission associated with nascent star formation and/or nuclear activity from a heavily obscured AGN. I also note that the comparison of resolved Mid-IR and radio images of galaxies on sub-galactic (kpc) scales can also be useful in identifying and recognising artificial mid-IR emission from less advanced intermediate Type II/III civilisations. Nevertheless, from the bulk properties of the $\hat{G}$ sample, I conclude that Kardashev Type-III civilisations are either very rare or do not exist in the local Universe.
We present a Python package LDTk that automates the calculation of custom stellar limb darkening (LD) profiles and model-specific limb darkening coefficients (LDC) using the library of PHOENIX-generated specific intensity spectra by Husser et al. (2013). The aim of the package is to facilitate analyses requiring custom generated limb darkening profiles, such as the studies of exoplanet transits--especially transmission spectroscopy, where the transit modelling is carried out for custom narrow passbands--eclipsing binaries (EBs), interferometry, and microlensing events. First, LDTk can be used to compute custom limb darkening profiles with uncertainties propagated from the uncertainties in the stellar parameter estimates. Second, LDTk can be used to estimate the limb-darkening-model specific coefficients with uncertainties for the most common limb-darkening models. Third, LDTk can be directly integrated into the log posterior computation of any pre-existing modelling code with minimal modifications. The last approach can be used to constrain the LD model parameter space directly by the LD profile, allowing for the marginalization over the LD parameter space without the need to approximate the constraint from the LD profile using a prior.
Ground and space-based sky surveys enable powerful cosmological probes based on measurements of galaxy properties and the distribution of galaxies in the Universe. These probes include weak lensing, baryon acoustic oscillations, abundance of galaxy clusters, and redshift space distortions; they are essential to improving our knowledge of the nature of dark energy. On the theory and modeling front, large-scale simulations of cosmic structure formation play an important role in interpreting the observations and in the challenging task of extracting cosmological physics at the needed precision. These simulations must cover a parameter range beyond the standard six cosmological parameters and need to be run at high mass and force resolution. One key simulation-based task is the generation of accurate theoretical predictions for observables, via the method of emulation. Using a new sampling technique, we explore an 8-dimensional parameter space including massive neutrinos and a variable dark energy equation of state. We construct trial emulators using two surrogate models (the linear power spectrum and an approximate halo mass function). The new sampling method allows us to build precision emulators from just 26 cosmological models and to increase the emulator accuracy by adding new sets of simulations in a prescribed way. This allows emulator fidelity to be systematically improved as new observational data becomes available and higher accuracy is required. Finally, using one LCDM cosmology as an example, we study the demands imposed on a simulation campaign to achieve the required statistics and accuracy when building emulators for dark energy investigations.
We present the results of near-infrared H- and K-band European Southern Observatory SINFONI integral field spectroscopy of the Seyfert galaxy NGC 1566. We investigate the central kpc of this nearby galaxy, concentrating on excitation conditions, morphology, and stellar content. NGC 1566 was selected from our NUGA (-south) sample and is a ringed, spiral galaxy with a stellar bar. We present emission and absorption line measurements in the central kpc of NGC 1566. Broad and narrow Br{\gamma} lines were detected. The detection of a broad Br{\gamma} component is a clear sign of a super-massive black hole in the center. Blackbody emission temperatures of ~1000 K are indicative of a hot dust component, the torus, in the nuclear region. The molecular hydrogen lines, hydrogen recombination lines, and [FeII] indicate that the excitation at the center is coming from an AGN. The central region is predominantly inhabited by molecular gas, dust, and an old K-M type giant stellar population. The molecular gas and stellar velocity maps both show a rotation pattern. The molecular gas velocity field shows a perturbation toward the center that is typical for bars or spiral density waves. The molecular gas species of warm H_2(1-0)S(1) and cold ^{12}CO(3-2) gas trace a nuclear gas disk of about 3" in radius with a nuclear spiral reaching toward the nucleus. From the equivalent width of H_2(1-0)S(1) a molecular ring with r<~3" can be inferred. This spiral seems to be an instrument that allows gas to fall toward the nucleus down to <50 pc scales. The excitation of molecular hydrogen in the nuclear gas disk is not clear but diagnostic diagrams show a distinction between the nuclear region and a <9 Myr old star forming region at the southwestern spiral arm. Possibly shocked gas is detected ~2" from the center, which is visible in dispersion maps of H$_2$(1-0)S(1) and ^{12}CO(3-2) and in the 0.87 mm continuum.
We adapted a multi-species escape model, developed for close-in extrasolar planets, to calculate the escape rates of CH4 and N2 from Pluto. In the absence of escape, CH4 should overtake N2 as the dominant species below the exobase. The CH4 profile depends strongly on the escape rate, however, and the typical escape rates predicted for Pluto lead to a nearly constant mixing ratio of less than 1 % below the exobase. In this case the CH4 escape rate is only 5-10 % of the N2 escape rate. Observations of the CH4 profile by the New Horizons/ALICE spectrograph can constrain the CH4 escape rate and provide a unique test for escape models.
Cosmic-ray scattering on magnetic turbulence leads to spatial diffusive propagation; if the scattering medium is moving, this will inevitably also cause changes in the momentum of the particles, so-called diffusive reacceleration. This can be described as diffusion in momentum space. Diffusive reacceleration has often been invoked to explain the peak observed in secondary-to-primary ratios at a few GeV, in particular Boron-to-Carbon. This avoids the necessity to postulate an ad-hoc break in the spatial diffusive coefficient, and has become almost a standard in modelling cosmic-ray spectra. However, at the levels invoked, the process implies a significant input of energy from the interstellar medium into cosmic rays, so that in such models interstellar space competes with the usual accelerators like supernova remnants. The questions arise: is reacceleration really occurring at the high level required to explain secondary-to-primary ratios? and are the energy requirements physically plausible? We address this issue using both analytical and numerical models of cosmic-ray propagation.
We analyze the shapes of cosmic string loops found in large-scale simulations of an expanding-universe string network. The simulation does not include gravitational back reaction, but we model that process by smoothing the loop using Lorentzian convolution. We find that loops at formation consist of generally straight segments separated by kinks. We do not see cusps or any cusp-like structure at the scale of the entire loop, although we do see very small regions of string that move with large Lorentz boosts. However, smoothing of the string almost always introduces two cusps on each loop. The smoothing process does not lead to any significant fragmentation of loops that were in non-self-intersecting trajectories before smoothing.
We derive and present a fast and accurate solution of the initial value problem for Keplerian motion in universal variables that does not use the Stumpff series. We find that it performs better than methods based on the Stumpff series.
We present a two-dimensional mapping of the gas flux distributions, as well as of the gas and stellar kinematics in the inner 220 pc of the Seyfert galaxy NGC 2110, using K-band integral field spectroscopy obtained with the Gemini NIFS at a spatial resolution of ~24pc and spectral resolution of ~40 km/s. The H2 emission extends over the whole field-of-view and is attributed to heating by X-rays from the AGN and/or by shocks, while the Brgamma emission is restricted to a bi-polar region extending along the South-East-North-West direction. The masses of the warm molecular gas and of the ionized gas are ~1.4x10^3 Msun and ~1.8x10^6 Msun, respectively. The stellar kinematics present velocity dispersions reaching 250km/s and a rotation pattern reaching an amplitude of 200 km/s. The gas velocity fields present a similar rotation pattern but also additional components that we attribute to inflows and outflows most clearly observed in the molecular gas emission. The inflows are observed beyond the inner 70 pc and are associated to a spiral arm seen in blueshift to the North-East and another in redshift to the South-West. We have estimated a mass inflow rate in warm molecular gas of ~4.6x10^-4 Msun/year. Within the inner 70 pc, another kinematic component is observed in the H2 emission that can be interpreted as due to a bipolar nuclear outflow oriented along the East-West direction, with a mass-outflow rate of ~4.3x10^-4 Msun/year in warm H2.
We use high resolution direct numerical simulations to show that helical turbulence can generate large-scale fields even in the presence of strong small-scale fields.During the kinematic stage, the unified large/small-scale dynamo grows fields with a shape-invariant eigenfunction, with most power peaked at small scales or large $k$. Nevertheless, the large-scale field can be clearly detected as an excess power at small $k$ in the negatively polarized component of the energy spectrum for a forcing with positively polarized waves. The strength of such kinematic large-scale field $\overline{B}$ relative to the total rms field $B_{rms}$ decreases with increasing magnetic Reynolds number, $Re_{M}$. However, as the Lorentz force becomes important, the field orders itself by saturating on successively larger scales. The magnetic power spectrum in the saturated state shows peaks at both the forcing wavenumber $k=k_f$, and at the box scale, $k=1$. The magnetic integral scale for the positively polarized waves, increases significantly from the kinematic stage to saturation. This implies that the small-scale field becomes as coherent as possible for a given forcing scale. Such an increase in the coherence scale of small-scale fields away from resistive scales averts the $Re_{M}$-dependent quenching of $\overline{B}/B_{rms}$. The $\overline{B}$, whose energy is measured in terms of the energy at $k=1$--$2$, grows from a value of the order of $4\%$ to about $40\%$ of $B_{rms}$ at saturation, aided in the final stages by helicity dissipation. Our results confirm that in helical turbulence, there is a single unified dynamo, with all scales initially growing together at one rate and, as the Lorentz force becomes important, successively larger scales saturate, with the largest scales continuing to grow (aided by small-scale magnetic helicity loss) as the small-scale field saturates.
C$_3$P$^-$ is analogous to the known interstellar anion C$_3$N$^-$ with phosphorus replacing the nitrogen in a simple step down the periodic table. In this work, it is shown that C$_3$P$^-$ is likely to possess a dipole-bound excited state. It has been hypothesized and observationally supported that dipole-bound excited states are an avenue through which anions could be formed in the interstellar medium. Additionally, C$_3$P$^-$ has a valence excited state that may lead to further stabilization of this molecule, and C$_3$P$^-$ has a larger dipole moment than neutral C$_3$P ($\sim 6$ D vs. $\sim 4$ D). As such, C$_3$P$^-$ is probably a more detectable astromolecule than even its corresponding neutral radical. Highly-accurate quantum chemical quartic force fields are also applied to C$_3$P$^-$ and its singly $^{13}$C substituted isotopologues in order to provide structures, vibrational frequencies, and spectroscopic constants that may aid in its detection.
The recently observed IceCube PeV events could be due to heavy dark matter (DM) decay. In this paper, we propose a simple DM model with extra $U(1)_X$ gauge symmetry and bridge it with standard model particles through heavy right-handed neutrino. The Dirac fermion DM $\chi$ with mass ~5 PeV can dominantly decay into a dark Higgs ($\phi$), the SM Higgs ($h$) and a neutrino ($\nu$). If the lifetime of $\chi$ is ~O($10^{28}$) sec, the resulting neutrino flux can fit data consistently. The neutrino flux from $\chi \rightarrow \phi h \nu$ in our model is softer than the one predicted from $\chi \rightarrow \nu h$, for example. We also discuss a possible mechanism to produce DM with the right relic abundance.
We consider the possibility of a gravitationally induced particle production through the mechanism of a nonminimal curvature-matter coupling. An interesting feature of this gravitational theory is that the divergence of the energy-momentum tensor is nonzero. As a first step in our study we reformulate the model in terms of an equivalent scalar-tensor theory, with two arbitrary potentials. By using the formalism of open thermodynamic systems, we interpret the energy balance equations in this gravitational theory from a thermodynamic point of view, as describing irreversible matter creation processes. The particle number creation rates, the creation pressure, and the entropy production rates are explicitly obtained as functions of the scalar field and its potentials, as well as of the matter Lagrangian. The temperature evolution laws of the newly created particles are also obtained. The cosmological implications of the model are briefly investigated, and it is shown that the late-time cosmic acceleration may be due to particle creation processes. Furthermore, it is also shown that due to the curvature--matter coupling, during the cosmological evolution a large amount of comoving entropy is also produced.
We present a simple model for electron transport in a possible layer of exotic nuclear clusters (in the so called nuclear pasta layer) between the crust and liquid core of a strongly magnetized neutron star. The electron transport there can be strongly anisotropic and gyrotropic. The anisotropy is produced by different electron effective collision frequencies along and across local symmetry axis in domains of exotic ordered nuclear clusters and by complicated effects of the magnetic field. We also calculate averaged kinetic coefficients in case local domains are freely oriented. Possible applications of the obtained results and open problems are outlined.
The post-Newtonian formulation of a general class of f(R) theories is set up to 3rd order approximation. It turns out that the information of a specific form of f(R) gravity is encoded in the Yukawa potential, which is contained in the perturbative expansion of the metric components. It is shown that the Yukawa potantial does appear in the 3rd order expression of the effective refraction index of light, although it is cancelled in the 2nd order expression. Therefore the f(R) theories are distinguishable from general relativity by gravitational lensing effect at the 3rd order post-Newtonian approximation. Our result opens the possibility to bring new insights into the issue of dark matter from f(R) gravity.
In this paper, we analyze a Bianchi Type I model with a scalar field in a
chaotic inflation potential, $V(\phi) = \frac{1}{2}\phi^2$ in the context of
stochastic eternal inflation. We use the typical slow-roll approximation in
combination with expansion-normalized variables in an orthonormal frame
approach to obtain a dynamical system which describes the dynamics of the shear
anisotropy and the inflaton field. We first show that the dynamics of the
inflaton field can be decoupled from the dynamics of the shear anisotropy. We
then use a fixed-points analysis in combination with global techniques from
topological dynamical systems theory to prove that the cosmological model under
consideration isotropizes irrespective of an inflationary epoch, which has also
described by other authors who have investigated a Bianchi Type I model under
similar configurations. We then show that for inflation to occur, the amount of
anisotropy must be very small.
We also give a description of the stochastic dynamics of the inflaton field
by using techniques from stochastic calculus. We show that the Klein-Gordon
equation becomes a stochastic differential equation with a highly nonlinear
drift term. In this case, the deceleration parameter itself becomes a random
variable, and we give details regarding when such a model can undergo
inflation. We finally derive the form of the long-term, stationary probability
distribution of the inflaton field, and show that it has the form of a
double-well potential. We then calculate the probability of inflation occurring
based on this approach. We conclude the paper by performing some numerical
simulations of the stochastic differential equation describing the dynamics of
the inflaton field. We conjecture that even in the case of stochastic eternal
inflation, one requires precise initial conditions for inflation to occur.
Links to: arXiv, form interface, find, astro-ph, recent, 1508, contact, help (Access key information)
In giant molecular clouds (GMCs), shocks driven by converging turbulent flows create high-density, strongly-magnetized regions that are locally sheetlike. In previous work, we showed that within these layers, dense filaments and embedded self-gravitating cores form by gathering material along the magnetic field lines. Here, we extend the parameter space of our three-dimensional, turbulent MHD core formation simulations. We confirm the anisotropic core formation model we previously proposed, and quantify the dependence of median core properties on the pre-shock inflow velocity and upstream magnetic field strength. Our results suggest that bound core properties are set by the total dynamic pressure (dominated by large-scale turbulence) and thermal sound speed c_s in GMCs, independent of magnetic field strength. For models with Mach number between 5 and 20, the median core masses and radii are comparable to the critical Bonnor-Ebert mass and radius defined using the dynamic pressure for P_ext. Our results correspond to M_core = 1.2 c_s^4/sqrt(G^3 rho_0 v_0^2) and R_core = 0.34 c_s^2/sqrt(G rho_0 v_0^2) for rho_0 and v_0 the large-scale mean density and velocity. For our parameter range, the median M_core ~ 0.1-1 M_sun, but a very high pressure cloud could have lower characteristic core mass. We find cores and filaments form simultaneously, and filament column densities are a factor ~2 greater than the surrounding cloud when cores first collapse. We also show that cores identified in our simulations have physical properties comparable to those observed in the Perseus cloud. Superthermal cores in our models are generally also magnetically supercritical, suggesting that the same may be true in observed clouds.
Studies of molecular clouds and young stars near the sun have provided invaluable insights into the process of star formation. Indeed, much of our physical understanding of this topic has been derived from such studies. Perhaps the two most fundamental problems confronting star formation research today are: 1) determining the origin of stellar mass and 2) deciphering the nature of the physical processes that control the star formation rate in molecular gas. As I will briefly outline here, observations and studies of local star forming regions are making particularly significant contributions toward the solution of both these important problems.
Strongly-magnetized, rapidly-rotating neutron stars are contenders for the central engines of both long-duration gamma-ray bursts (LGRBs) and hydrogen-poor super-luminous supernovae (SLSNe-I). Models for typical (~minute long) LGRBs invoke magnetars with high dipole magnetic fields (Bd > 1e15 G) and short spin-down times, while models for SLSNe-I invoke neutron stars with weaker fields and longer spin-down times of weeks. Here we identify a transition region in the space of Bd and birth period for which a magnetar can power both a long GRB and a luminous SN. In particular, we show that a 2 ms period magnetar with a spin-down time of ~1e4 s can explain the observations of both the ultra-long GRB 111209 and its associated luminous SN2011kl. For magnetars with longer spin down times, we predict even longer duration (~1e6 s) GRBs and brighter supernovae, a correlation that extends to Swift J2058+05 (commonly interpreted as a tidal disruption event). We further show that previous estimates of the maximum rotational energy of a proto-magnetar were too conservative and energies up to Emax ~1-2e53 erg are possible. The magnetar model can therefore comfortably accommodate the extreme energy requirements recently posed by the most luminous supernova ASASSN-15lh. The high ionization flux from a pulsar wind nebula powering ASASSN-15lh may lead to an "ionization break-out" X-ray burst over the coming months, which would be accompanied by an abrupt change in the optical spectrum. We conclude by briefly contrasting millisecond magnetar and black hole models for SLSNe and ultra-long GRBs.
We have recently introduced a novel statistical measure of dark matter clustering in phase space, the particle phase space average density ($P^2SAD$). In a two-paper series, we studied the structure of $P^2SAD$ in the Milky-Way-size Aquarius haloes, constructed a physically motivated model to describe it, and illustrated its potential as a powerful tool to predict signals sensitive to the nanostructure of dark matter haloes. In this letter, we report a remarkable universality of the clustering of dark matter in phase space as measured by $P^2SAD$ within the subhaloes of host haloes across different environments covering a range from dwarf-size to cluster-size haloes ($10^{10}-10^{15}$ M$_\odot$). Simulations show that the universality of $P^2SAD$ holds for more than 7 orders of magnitude, over a 2D phase space, covering over 3 orders of magnitude in distance/velocity, with a simple functional form that can be described by our model. Invoking the universality of $P^2SAD$, we can accurately predict the non-linear power spectrum of dark matter at small scales all the way down to the decoupling mass limit of cold dark matter particles. As an application, we compute the subhalo boost to the annihilation of dark matter in a wide range of host halo masses.
From stellar spectra, a variety of physical properties of stars can be derived. In particular, the chemical composition of stellar atmospheres can be inferred from absorption line analyses. These provide key information on large scales, such as the formation of our Galaxy, down to the small-scale nucleosynthesis processes that take place in stars and supernovae. By extending the observed wavelength range toward bluer wavelengths, we optimize such studies to also include critical absorption lines in metal-poor stars, and allow for studies of heavy elements (Z>38) whose formation processes remain poorly constrained. In this context, spectrographs optimized for observing blue wavelength ranges are essential, since many absorption lines at redder wavelengths are too weak to be detected in metal-poor stars. This means that some elements cannot be studied in the visual-redder regions, and important scientific tracers and science cases are lost. The present era of large public surveys will target millions of stars. Here we describe the requirements driving the design of the forthcoming survey instrument 4MOST, a multi-object spectrograph commissioned for the ESO VISTA 4m-telescope. We focus here on high-density, wide-area survey of stars and the science that can be achieved with high-resolution stellar spectroscopy. Scientific and technical requirements that governed the design are described along with a thorough line blending analysis. For the high-resolution spectrograph, we find that a sampling of >2.5 (pixels per resolving element), spectral resolution of 18000 or higher, and a wavelength range covering 393-436 nm, is the most well-balanced solution for the instrument. A spectrograph with these characteristics will enable accurate abundance analysis (+/-0.1 dex) in the blue and allow us to confront the outlined scientific questions. (abridged)
We use a semi-analytical model for the substructure of dark matter haloes to assess the too-big-to-fail (TBTF) problem. The model accurately reproduces the average subhalo mass and velocity functions, as well as their halo-to-halo variance, in N-body simulations. We construct thousands of realizations of Milky Way (MW) size host haloes, allowing us to investigate the TBTF problem with unprecedented statistical power. We examine the dependence on host halo mass and cosmology, and explicitly demonstrate that a reliable assessment of TBTF requires large samples of hundreds of host haloes. We argue that previous statistics used to address TBTF suffer from the look-elsewhere effect and/or disregard certain aspects of the data on the MW satellite population. We devise a new statistic that is not hampered by these shortcomings, and, using only data on the 9 known MW satellite galaxies with $V_{\rm max}>15{\rm kms}^{-1}$, demonstrate that $1.4^{+3.3}_{-1.1}\%$ of MW-size host haloes have a subhalo population in statistical agreement with that of the MW. However, when using data on the MW satellite galaxies down to $V_{\rm max}=8{\rm kms}^{-1}$, this MW consistent fraction plummets to $<5\times10^{-4}$ (at 68% CL). Hence, if it turns out that the inventory of MW satellite galaxies is complete down to 8km/s, then the maximum circular velocities of MW satellites are utterly inconsistent with $\Lambda$CDM predictions, unless baryonic effects can drastically increase the spread in $V_{\rm max}$ values of satellite galaxies compared to that of their subhaloes.
We study the structure and kinematics of the OH-streamer and the +80 km/s cloud and their interactions with the circumnuclear disk (CND) and with other molecular clouds in the vicinity of the Galactic centre (GC), and we map OH absorption at about 6" resolution at R $\le$ 10 pc from the GC, with about 9 km/s velocity resolution. The VLA was used to map OH line absorption at the 1665 and 1667 MHz lambda doublet main lines towards the Sagittarius A complex. Strong OH absorption was found in the OH-streamer, the southern streamer (SS), the +20, +50, and +80 km/s molecular clouds, the molecular belt, the CND, the expanding molecular ring (EMR), and the high negative velocity gas (HNVG). The OH-streamer was found to comprise three parts, head, middle, and tail, and to interact with the SS/+20, +80 km/s clouds and the CND. Optical depths and column densities have been calculated for the OH-streamer and the +80 km/s cloud. The OH-streamer, the SS, the +20 and +80 km/s clouds, and the CND are intimately related in position and velocity space. The OH-streamer was found to be a clumpy object stretching in projection from the inner radius of the CND at about 1.8 pc from Sgr A* towards and partly engulfing Sgr A*. As a side result of our data, a possible link between the near side of the EMR and the CND's southwest lobe was found. Additionally, we found OH absorption against all four of the previously known Compact HII Regions A-D, located east of Sgr A East, indicating their close association with the +50 km/s cloud.
The perfect 30-min cadence of the full-frame images from the Transiting Exoplanet Survey Satellite (TESS) will impose a hard Nyquist limit of 24 d$^{-1}$ ($\approx 278$ $\mu$Hz). This will be problematic for asteroseismology of stars with oscillation frequencies at or around that Nyquist limit, which will have insurmountable Nyquist ambiguities. TESS does offer some observing slots at shorter cadences, but these will be limited in number and competitive, while the full frame images will be the main data product for many types of variable stars. We show that the Nyquist ambiguities can be alleviated if, when TESS resumes observations after a downlink, integrations are not resumed at perfect cadence with those before the downlink. The time spent idling before integrations are resumed need only be around five minutes for satisfactory results, and observing time can be recouped from the downlink event if the telescope does not wait for a return to perfect cadence before resuming integrations. The importance of imperfect cadence after downlink is discussed in light of phase coverage of transit events.
Relativistic jets naturally occur in astrophysical systems that involve accretion onto compact objects, such as core collapse of massive stars in gamma-ray bursts (GRBs) and accretion onto supermassive black holes in active galactic nuclei (AGN). It is generally accepted that these jets are powered electromagnetically, by the magnetised rotation of a central compact object. However, how they produce the observed emission and survive the propagation for many orders of magnitude in distance without being disrupted by current-driven non-axisymmetric instabilities is the subject of active debate. We carry out time-dependent 3D relativistic magnetohydrodynamic simulations of relativistic, Poynting flux dominated jets. The jets are launched self-consistently by the rotation of a strongly magnetised central compact object. This determines the natural degree of azimuthal magnetic field winding, a crucial factor that controls jet stability. We find that the jets are susceptible to two types of instability: (i) a global, external kink mode that grows on long time scales and causes the jets to bodily bend sideways. Whereas this mode does not cause jet disruption over the simulated distances, it substantially reduces jet propagation speed. We show, via an analytic model, that the growth of the external kink mode depends on the slope of the ambient medium density profile. In flat density distributions characteristic of galactic cores, an AGN jet may stall, whereas in stellar envelopes the external kink weakens as the jet propagates outward; (ii) a local, internal kink mode that grows over short time scales and causes small-angle magnetic reconnection and conversion of about half of jet electromagnetic energy flux into heat. Based on the robustness and energetics of the internal kink mode, we suggest that this instability is the main dissipation mechanism responsible for powering GRB prompt emission.
Classical Cepheids are key probes of both stellar astrophysics and cosmology as standard candles and pulsating variable stars. It is important to understand Cepheids in unprecedented detail in preparation for upcoming GAIA, JWST and extremely-large telescope observations. Cepheid eclipsing binary stars are ideal tools for achieving this goal, however there are currently only three known systems. One of those systems, OGLE-LMC-CEP1812, raises new questions about the evolution of classical Cepheids because of an apparent age discrepancy between the Cepheid and its red giant companion. We show that the Cepheid component is actually the product of a stellar merger of two main sequence stars that has since evolved across the Hertzsprung gap of the HR diagram. This post-merger product appears younger than the companion, hence the apparent age discrepancy is resolved. We discuss this idea and consequences for understanding Cepheid evolution.
We report on initial results from a campaign to obtain optical imaging of a sample of Ultra Compact High Velocity Clouds (UCHVCs) discovered by the ALFALFA neutral hydrogen (HI) survey. UCHVCs are sources with velocities and sizes consistent with their being low-mass dwarf galaxies in the Local Volume, but without optical counterparts in existing catalogs. We are using the WIYN 3.5-m telescope and pODI camera to image these objects and search for an associated stellar population. In this paper, we present our observational strategy and method for searching for resolved stellar counterparts to the UCHVCs. We combine careful photometric measurements, a color-magnitude filter, and spatial smoothing techniques to search for stellar overdensities in the g- and i-band images. We also run statistical tests to quantify the likelihood that whatever overdensities we find are real and not chance superpositions of sources. We demonstrate the method by applying it to two data sets: WIYN imaging of Leo P, a UCHVC discovered by ALFALFA and subsequently shown to be a low-mass star-forming dwarf galaxy in the Local Volume, and WIYN imaging of AGC198606, an HI cloud identified by ALFALFA that is near in position and velocity to the Local Group dwarf Leo T. Applying the search method to the Leo P data yields an unambiguous detection (>99% confidence) of the galaxy's stellar population. Applying our method to the AGC198606 imaging yields a possible detection (92% confidence) of an optical counterpart located ~2.5 arc minutes away from the centroid of AGC198606's HI distribution and within the HI disk. We use the optical data to estimate a distance to the stellar counterpart between 373 and 393 kpc, with an absolute magnitude M_i = -4.67+/-0.09. Combining the WIYN data with our previous estimate of the HI mass of AGC198606 from WSRT imaging yields an HI-to-stellar mass ratio of ~45-110.
The performance of a wide-field adaptive optics system depends on input design parameters. Here we investigate the performance of a multi-conjugate adaptive optics system design for the European Extremely Large Telescope, using an end-to-end Monte-Carlo adaptive optics simulation tool, DASP. We consider parameters such as the number of laser guide stars, sodium layer depth, wavefront sensor pixel scale, number of deformable mirrors, mirror conjugation and actuator pitch. We provide potential areas where costs savings can be made, and investigate trade-offs between performance and cost. We conclude that a 6 laser guide star system using 3 DMs seems to be a sweet spot for performance and cost compromise.
Gas-phase abundances in HII regions of two spiral galaxies, NGC7793 and NGC4945, have been studied to determine their radial metallicity gradients. We used the strong-line method to derive oxygen abundances from spectra acquired with GMOS-S, the multi-object spectrograph on the 8m- Gemini South telescope. We found that NGC7793 has a well-defined gas-phase radial oxygen gradient of -0.321 $\pm$ 0.112 dex R$_{\rm 25}^{-1}$ (or -0.054 $\pm$ 0.019 dex kpc$^{-1}$) in the galactocentric range 0.17$<$R$_{\rm G}$/R$_{\rm 25}$ $<$ 0.82, not dissimilar from gradients calculated with direct abundance methods in galaxies of similar mass and morphology. We also determined a shallow radial oxygen gradient in NGC4945, -0.253 $\pm$ 0.149 dex R$_{\rm 25}^{-1}$ (or -0.019 $\pm$ 0.011 dex kpc$^{-1}$) for 0.04$<$R$_{\rm G}$/R$_{\rm 25}$ $<$ 0.51, where the larger relative uncertainty derives mostly from the larger inclination of this galaxy. NGC7793 and NGC4945 have been selected for this study because they are similar, in mass and morphology, to M33 and the Milky Way, respectively. Since at zeroth order we expect the radial metallicity gradients to depend on mass and galaxy type, we compared our galaxies in the framework of radial metallicity models best suited for M33 and the Galaxy. We found a good agreement between M33 and NGC7793, pointing toward similar evolution for the two galaxies. We notice instead differences between NGC4945 and the radial metallicity gradient model that best fits the Milky Way. We found that these differences are likely related to the presence of an AGN combined with a bar in the central regions of NGC4945, and to its interacting environment.
We report on a detailed analysis of the so-called ~1 Hz quasi-periodic oscillation (QPO) in the eclipsing and dipping neutron-star low-mass X-ray binary EXO 0748-676. This type of QPO has previously been shown to have a geometric origin. Our study focuses on the evolution of the QPO as the source moves through the color-color diagram, in which it traces out an atoll-source-like track. The QPO frequency increases from ~0.4 Hz in the hard state to ~25 Hz as the source approaches the soft state. Combining power spectra based on QPO frequency reveals additional features that strongly resemble those seen in non-dipping/eclipsing atoll sources. We show that the low-frequency QPOs in atoll sources and the ~1 Hz QPO in EXO 0748-676 follow similar relations with respect to the noise components in their power spectra. We conclude that the frequencies of both types of QPOs are likely set by (the same) precession of a misaligned inner accretion disk. For high-inclination systems, like EXO 0748-676, this results in modulations of the neutron-star emission due to obscuration or scattering, while for lower-inclination systems the modulations likely arise from relativistic Doppler boosting and light-bending effects.
We study the orbital evolution of hot Jupiters due to the excitation and damping of tidally driven $g$-modes within solar-type host stars. Linearly resonant $g$-modes (the dynamical tide) are driven to such large amplitudes in the stellar core that they excite a sea of other $g$-modes through weakly nonlinear interactions. By solving the dynamics of large networks of nonlinearly coupled modes, we show that the nonlinear dissipation rate of the dynamical tide is several orders of magnitude larger than the linear dissipation rate. As a result, we find that the orbits of planets with mass $M_p > 0.5M_J$ and period $P < 2\textrm{ days}$ decay on timescales that are small compared to the main-sequence lifetime of their solar-type hosts. This corresponds to stellar tidal quality factors $Q_\ast^\prime \simeq 10^5-10^6$ for this range of $M_p$ and $P$. Our results imply that there are $\simeq10$ currently known exoplanetary systems, including WASP-19b and HAT-P-36-b, with orbital decay timescales shorter than a Gyr. Rapid, tide induced orbital decay may explain the observed paucity of planets with $M_p > M_J$ and $P< 2\textrm{ days}$ around solar-type hosts and could generate detectable transit-timing variations in the near future.
Any non-spherical distribution of density inside planets and stars gives rise to a non-spherical external gravity and change of shape. If part or all of the observed zonal flows at the cloud deck of giant planets represent deep interior dynamics, then the density perturbations associated with the deep zonal flows could generate gravitational signals detectable by the planned Juno mission and the Cassini Proximal Orbits. It is currently debated whether the thermal wind equation (TWE) can be used to calculate the gravity field associated with deep zonal flows. Here we present a critical comparison between the Euler equation and the thermal wind equation. Our analysis shows that the applicability of the TWE in calculating the gravity moments depends crucially on retaining the non-sphericity of the background density and gravity. Only when the background non-sphericity of the planet is taken into account, the TWE makes accurate enough prediction (with a few tens of percent errors) for the high-degree gravity moments associated with deep zonal flows. Since the TWE is derived from the curl of the Euler equation and is a local relation, it necessarily says nothing about any density perturbations that contribute irrotational terms to the Euler equation and that has a non-local origin. However, the predicted corrections from these density contributions to the low harmonic degree gravity moments are not discernible from insignificant changes in interior models while the corrections at high harmonic degree are very small, tens of percent or less.
We present polarization observations of two Galactic plane fields centered on Galactic coordinates (l,b)=(0 deg,0 deg) and (329 deg, 0 deg) at Q- (43 GHz) and W-band (95 GHz), covering between 301 and 539 square degrees depending on frequency and field. These measurements were made with the QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hours of observations. The resulting maps represent the deepest large-area Galactic polarization observations published to date at the relevant frequencies with instrumental rms noise varying between 1.8 and 2.8 uK deg, 2.3-6 times deeper than corresponding WMAP and Planck maps. The angular resolution is 27.3' and 12.8' FWHM at Q- and W-band, respectively. We find excellent agreement between the QUIET and WMAP maps over the entire fields, and no compelling evidence for significant residual instrumental systematic errors in either experiment, whereas the Planck 44 GHz map deviates from these in a manner consistent with reported systematic uncertainties for this channel. We combine QUIET and WMAP data to compute inverse-variance-weighted average maps, effectively retaining small angular scales from QUIET and large angular scales from WMAP. From these combined maps, we derive constraints on several important astrophysical quantities, including a robust detection of polarized synchrotron spectral index steepening of ~0.2 off the plane, as well as the Faraday rotation measure toward the Galactic center (RM=-4000 +/- 200 rad m^-2), all of which are consistent with previously published results. Both the raw QUIET and the co-added QUIET+WMAP maps are made publicly available together with all necessary ancillary information.
With the hypothesis that cosmic string loops act as seeds for globular clusters in mind, we study the role that velocities of these strings will play in determining the mass distribution of globular clusters. Loops with high enough velocities will not form compact and roughly spherical objects and can hence not be the seeds for globular clusters. We compute the expected number density and mass function of globular clusters as a function of both the string tension and the peak loop velocity, and compare the results with the observational data on the mass distribution of globular clusters in our Milky Way. We determine the critical peak string loop velocity above which the agreement between the string loop model for the origin of globular clusters (neglecting loop velocities) and observational data is lost.
The Kepler mission has discovered more than 4000 exoplanet candidates. Many are in systems with tightly packed inner planets. Inside-Out Planet Formation (IOPF) has been proposed to explain these systems. It involves sequential in situ planet formation at the local pressure maximum of a retreating dead zone inner boundary (DZIB). Pebbles accumulate at this pressure trap, which builds up a ring, and then a planet. The planet is expected to grow until it opens a gap, which helps to both truncate pebble accretion and induce DZIB retreat that sets the location of formation of the next planet. This simple scenario may be modified if the planet migrates significantly from its formation location. Thus planet-disk interactions play a crucial role in the IOPF scenario. We present numerical simulations that first assess migration of planets of various masses that are forming at the DZIB of an active accretion disk, where the effective viscosity rapidly increases in the radially inward direction. We find that the disk's torques on the planet tend to trap the planet at a location very close to the initial pressure maximum where it formed. We then study gap opening by these planets to assess at what mass a significant gap is created. Finally we present a simple model for DZIB retreat due to penetration of X-rays from the star to the disk midplane. Overall, these simulations help to quantify both the mass scale of first,"Vulcan," planet formation and the orbital separation to the location of second planet formation.
Using all-sky maps obtained with COBE/DIRBE, we reanalyzed the diffuse sky brightness at 1.25 and 2.2 um, which consists of zodiacal light, diffuse Galactic light (DGL), integrated starlight (ISL), and isotropic emission including the extragalactic background light. Our new analysis including an improved estimate of the DGL and the ISL with the 2MASS data showed that deviations of the isotropic emission from isotropy were less than 10% in the entire sky at high Galactic latitude (|b|>35). The result of our analysis revealed a significantly large isotropic component at 1.25 and 2.2 um with intensities of 60.15 +/- 16.14 and 27.68 +/- 6.21 nWm-2sr-1, respectively. This intensity is larger than the integrated galaxy light, upper limits from gamma-ray observation, and potential contribution from exotic sources (i.e., Population III stars, intrahalo light, direct collapse black holes, and dark stars). We therefore conclude that the excess light may originate from the local universe; the Milky Way and/or the solar system.
We have observed high-dispersion echelle spectra of main-sequence stars in five nearby young associations -- Argus, Carina-Near, Hercules-Lyra, Orion and Subgroup B4 -- and derived abundances for elements ranging from Na to Eu. These are the first chemical abundance measurements for two of the five associations, while the remaining three associations are analysed more extensively in our study. Our results support the presence of chemical homogeneity among association members with a typical star-to-star abundance scatter of about 0.06 dex or less over many elements. The five associations show log$\epsilon$(Li) consistent with their age and share a solar chemical composition for all elements with the exception of Ba. We find that all the heavy elements (Y, Zr, La, Ce, Nd, Sm and Eu) exhibit solar ratios, i.e., [X/Fe] $\simeq$ 0, while Ba is overabundant by about 0.2-0.3 dex. The origin of the overabundance of Ba is a puzzle. Within the formulation of the s-process, it is difficult to create a higher Ba abundance without a similar increase in the s-process contributions to other heavy elements (La-Sm). Given that Ba is represented by strong lines of Ba II and La-Sm are represented by rather weak ionized lines, the suggestion, as previously made by other studies, is that the Ba abundance may be systematically overestimated by standard methods of abundance analysis perhaps because the upper reaches of the stellar atmospheres are poorly represented by standard model atmospheres. A novel attempt to analyse the Ba I line at 5535 \AA\ gives a solar Ba abundance for stars with effective temperatures hotter than about 5800 K but increasingly subsolar Ba abundances for cooler stars with apparent Ba deficiencies of 0.5 dex at 5100 K. This trend with temperature may signal a serious non-LTE effect on the Ba I line.
We present spectroscopic tests of MIR to FIR extinction laws in IRDC G028.36+00.07, a potential site of massive star and star cluster formation. Lim & Tan (2014) developed methods of FIR extinction mapping of this source using ${\it Spitzer}$-MIPS ${\rm 24\mu m}$ and ${\it Herschel}$-PACS ${\rm 70\mu m}$ images, and by comparing to MIR ${\it Spitzer}$-IRAC $3$--${\rm 8\mu m}$ extinction maps, found tentative evidence for grain growth in the highest mass surface density regions. Here we present results of spectroscopic infrared extinction (SIREX) mapping using ${\it Spitzer}$-IRS (14 to ${\rm 38\mu m}$) data of the same IRDC. These methods allow us to first measure the SED of the diffuse Galactic ISM that is in the foreground of the IRDC. We then carry out our primary investigation of measuring the MIR to FIR opacity law and searching for potential variations as a function of mass surface density within the IRDC. We find relatively flat, featureless MIR-FIR opacity laws that lack the $\sim{\rm 12\mu m}$ and $\sim{\rm 35\mu m}$ features associated with the thick water ice mantle models of Ossenkopf & Henning (1994). Their thin ice mantle models and the coagulating aggregate dust models of Ormel et al. (2011) are a generally better match to the observed opacity laws. We also find evidence for generally flatter MIR to FIR extinction laws as mass surface density increases, strengthening the evidence for grain and ice mantle growth in higher density regions.
Bimolecular diffusion coefficients are important parameters used by atmospheric models to calculate altitude profiles of minor constituents in an atmosphere. Unfortunately, laboratory measurements of these coefficients were never conducted at temperature conditions relevant to the atmosphere of Titan. Here we conduct a detailed uncertainty analysis of the bimolecular diffusion coefficient parameters as applied to Titan's upper atmosphere to provide a better understanding of the impact of uncertainty for this parameter on models. Because temperature and pressure conditions are much lower than the laboratory conditions in which bimolecular diffusion parameters were measured, we apply a Bayesian framework, a problem-agnostic framework, to determine parameter estimates and associated uncertainties. We solve the Bayesian calibration problem using the open-source QUESO library which also performs a propagation of uncertainties in the calibrated parameters to temperature and pressure conditions observed in Titan's upper atmosphere. Our results show that, after propagating uncertainty through the Massman model, the uncertainty in molecular diffusion is highly correlated to temperature and we observe no noticeable correlation with pressure. We propagate the calibrated molecular diffusion estimate and associated uncertainty to obtain an estimate with uncertainty due to bimolecular diffusion for the methane molar fraction as a function of altitude. Results show that the uncertainty in methane abundance due to molecular diffusion is in general small compared to eddy diffusion and the chemical kinetics description. However, methane abundance is most sensitive to uncertainty in molecular diffusion above 1200 km where the errors are nontrivial and could have important implications for scientific research based on diffusion models in this altitude range.
We present HI 21cm emission observations of the z ~ 0.00632 sub-damped Lyman-alpha absorber (sub-DLA) towards PG1216+069 made using the Arecibo Telescope and the Very Large Array (VLA). The Arecibo 21cm spectrum corresponds to an HI mass of ~ 3.2x10^7 solar masses, two orders of magnitude smaller than that of a typical spiral galaxy. This is surprising since in the local Universe the cross-section for absorption at high HI column densities is expected to be dominated by spirals. The 21cm emission detected in the VLA spectral cube has a low signal-to-noise ratio, and represents only half the total flux seen at Arecibo. Emission from three other sources is detected in the VLA observations, with only one of these sources having an optical counterpart. This group of HI sources appears to be part of complex "W", believed to lie in the background of the Virgo cluster. While several HI cloud complexes have been found in and around the Virgo cluster, it is unclear whether the ram pressure and galaxy harassment processes that are believed to be responsible for the creation of such clouds in a cluster environment are relevant at the location of this cloud complex. The extremely low metallicity of the gas, ~ 1/40 solar, also makes it unlikely that the sub-DLA consists of material that has been stripped from a galaxy. Thus, while our results have significantly improved our understanding of the host of this sub-DLA, the origin of the gas cloud remains a mystery
Tidal Disruption of stars by supermassive central black holes from dense rotating star clusters is modelled by high-accuracy direct N-body simulation. As in a previous paper on spherical star clusters we study the time evolution of the stellar tidal disruption rate and the origin of tidally disrupted stars, now according to several classes of orbits which only occur in axisymmetric systems (short axis tube and saucer). Compared with that in spherical systems, we found a higher TD rate in axisymmetric systems. The enhancement can be explained by an enlarged loss-cone in phase space which is raised from the fact that total angular momentum $\bf J$ is not conserved. As in the case of spherical systems, the distribution of the last apocenter distance of tidally accreted stars peaks at the classical critical radius. However, the angular distribution of the origin of the accreted stars reveals interesting features. Inside the influence radius of the supermassive black hole the angular distribution of disrupted stars has a conspicuous bimodal structure with a local minimum near the equatorial plane. Outside the influence radius this dependence is weak. We show that the bimodal structure of orbital parameters can be explained by the presence of two families of regular orbits, namely short axis tube and saucer orbits. Also the consequences of our results for the loss cone in axisymmetric galactic nuclei are presented.
Using six-dimesional phase-space information from the Fourth Data release of the Radial Velocity Experiment (RAVE) over the range of Galactic longitude 240$^{\circ}< l <$ 360$^{\circ}$ and $V_{LSR} < -239$ kms$^{-1}$, we have computed orbits for 329 RAVE stars that were originally selected as chemically and kinematically related to $\omega$ Centauri. The orbits were integrated in a Milky-Way-like axisymmetric Galactic potential, ignoring the effects of the dynamical evolution of $\omega$ Centauri due to the tidal effects of the Galaxy disk on the cluster along time. We also ignored secular changes in the Milky Way potential over time. In a Monte Carlo scheme, and under the assumption that the stars may have been ejected with velocities greater than the escape velocity ($V_{rel}>V_{esc,0}$) from the cluster, we identified 15 stars as having close encounters with $\omega$ Centauri: (\textit{i}) 8 stars with relative velocities $V_{rel}< 200 $ kms$^{-1}$ may have been ejected $\sim$ 200 Myr ago from $\omega$ Centauri; (\textit{ii}) other group of 7 stars were identified with high relative velocity $V_{rel}> 200 $ kms$^{-1}$ during close encounters, and seems unlikely that they have been ejected from $\omega$ Centauri. We also confirm the link between J131340.4-484714 as potential member of $\omega$ Centauri, and probably ejected $\sim$ 2.0 Myr ago, with a relative velocity $V_{rel}\sim80$ kms$^{-1}$.
The relative cosmic variance ($\sigma_v$) is a fundamental source of uncertainty in pencil-beam surveys and, as a particular case of count-in-cell statistics, can be used to estimate the bias between galaxies and their underlying dark-matter distribution. Our goal is to test the significance of the clustering information encoded in the $\sigma_v$ measured in the ALHAMBRA survey. We measure the cosmic variance of several galaxy populations selected with $B-$band luminosity at $0.35 \leq z < 1.05$ as the intrinsic dispersion in the number density distribution derived from the 48 ALHAMBRA subfields. We compare the observational $\sigma_v$ with the cosmic variance of the dark matter expected from the theory, $\sigma_{v,{\rm dm}}$. This provides an estimation of the galaxy bias $b$. The galaxy bias from the cosmic variance is in excellent agreement with the bias estimated by two-point correlation function analysis in ALHAMBRA. This holds for different redshift bins, for red and blue subsamples, and for several $B-$band luminosity selections. We find that $b$ increases with the $B-$band luminosity and the redshift, as expected from previous work. Moreover, red galaxies have a larger bias than blue galaxies, with a relative bias of $b_{\rm rel} = 1.4 \pm 0.2$. Our results demonstrate that the cosmic variance measured in ALHAMBRA is due to the clustering of galaxies and can be used to characterise the $\sigma_v$ affecting pencil-beam surveys. In addition, it can also be used to estimate the galaxy bias $b$ from a method independent of correlation functions.
This paper investigates the potential to detect tau neutrinos in the energy range of 1-1000 PeV searching for very inclined showers with imaging Cherenkov telescopes. A neutrino induced tau lepton escaping from the Earth may decay and initiate an air shower which can be detected by a fluorescence or Cherenkov telescope. We present here a study of the detection potential of Earth-skimming neutrinos taking into account neutrino interactions in the Earth crust, local matter distributions at various detector sites, the development of tau-induced showers in air and the detection of Cherenkov photons with IACTs. We analysed simulated shower images on the camera focal plane and implemented generic reconstruction chains based on Hillas parameters. We find that present IACTs can distinguish air showers induced by tau neutrinos from the background of hadronic showers in the PeV-EeV energy range. We present the neutrino trigger efficiency obtained for a few configurations being considered for the next-generation Cherenkov telescopes, i.e. the Cherenkov Telescope Array. Finally, for a few representative neutrino spectra expected from astrophysical sources, we compare the expected event rates at running IACTs to what expected for the dedicated IceCube neutrino telescope.
The ultra-high energy cosmic rays observed at the Earth are most likely accelerated in extra-galactic sources. For the typical luminosities invoked for such sources, the electric current associated to the flux of cosmic rays that leave them is large. The associated plasma instabilities create magnetic fluctuations that can efficiently scatter particles. We argue that this phenomenon forces cosmic rays to be self-confined in the source proximity for energies $E<E_{\rm cut}$, where $E_{\rm cut}\approx 10^{7} L_{44}^{2/3}$ GeV for low background magnetic fields ($B_{0}\ll nG$). For larger values of $B_{0}$, cosmic rays are confined close to their sources for energies $E<E_{\rm cut}\approx 2\times 10^{8} \lambda_{10} L_{44}^{1/4} B_{-10}^{1/2}$ GeV, where $B_{-10}$ is the field in units of $0.1$ nG, $\lambda_{10}$ is its coherence lengths in units of 10 Mpc and $L_{44}$ is the source luminosity in units of $10^{44}$ erg/s.
The mechanism through which meter-sized boulders grow to km-sized planetesimals in protoplanetary discs is a subject of active research, since it is critical for planet formation. To avoid spiralling into the protostar due to aerodynamic drag, objects must rapidly grow from cm-sized pebbles, which are tightly coupled to the gas, to large boulders of 1-100m in diameter. It is already well known that over-densities in the gaseous component of the disc provide potential sites for the collection of solids, and that significant density structures in the gaseous component of the disc (e.g., spiral density waves) can trap solids efficiently enough for the solid component of the disc to undergo further gravitational collapse due to their own self-gravity. In this work, we employ the PENCIL CODE to conduct local shearing sheet simulations of massive self-gravitating protoplanetary discs, to study the effect of anticyclonic transient vortices, or eddies, on the evolution of solids in these discs. We find that these types of structures are extremely efficient at concentrating small and intermediate-sized dust particles with friction times comparable to, or less than, the local orbital period of the disc. This can lead to significant over-densities in the solid component of the disc, with density enhancements comparable to, and even higher, than those within spiral density waves; increasing the rate of gravitational collapse of solids into bound structures.
We study a two component dark matter candidate inspired by the Minimal Walking Technicolor model. Dark matter consists of a dominant SIMP-like dark atom component made of bound states between primordial helium nuclei and a doubly charged technilepton, and a small WIMP-like component made of another dark atom bound state between a doubly charged technibaryon and a technilepton. This scenario is consistent with direct search experimental findings because the dominant SIMP component interacts too strongly to reach the depths of current detectors with sufficient energy to recoil and the WIMP-like component is too small to cause significant amount of events. In this context a metastable technibaryon that decays to $e^+e^+$, $\mu^+ \mu^+$ and $\tau^+ \tau^+$ can in principle explain the observed positron excess by AMS-02 and PAMELA, while being consistent with the photon flux observed by FERMI/LAT. We scan the parameters of the model and we find the best possible fit to the latest experimental data. We find that there is a small range of parameter space that this scenario can be realised under certain conditions regarding the cosmic ray propagation and the final state radiation. This range of parameters fall inside the region where the current run of LHC can probe, and therefore it will soon be possible to either verify or exclude conclusively this model of dark matter.
Diffuse radio emission in the form of radio halos and relics has been found in a number of merging galaxy clusters. These structures indicate that shock and turbulence associated with the merger accelerate electrons to relativistic energies. We report the discovery of a radio relic + radio halo system in PSZ1 G108.18-11.53 (z=0.335). This cluster hosts the second most powerful double radio relic system ever discovered. We observed PSZ1 G108.18-11.53 with the Giant Meterwave Radio Telescope (GMRT) and the Westerbork Synthesis Radio Telescope (WSRT). We obtained radio maps at 147, 323, 607 and 1380 MHz. We also observed the cluster with the Keck telescope, obtaining the spectroscopic redshift for 42 cluster members. From the injection index we obtained the Mach number of the shocks generating the two radio relics. For the southern shock we found M = 2.33^{+0.19}_{-0.26}, while the northern shock Mach number goes from M = 2.20^{+0.07}_{-0.14} in the north part down to M = 2.00^{+0.03}_{-0.08} in the southern region. If the relation between the injection index and the Mach number predicted by diffusive shock acceleration (DSA) theory holds, this is the first observational evidence for a gradient in the Mach number along a galaxy cluster merger shock.
KMOS (K-Band Multi Object Spectrograph) is a novel integral field spectrograph installed in the VLT's ANTU unit. The instrument offers an ability to observe 24 2.8"$\times$2.8" sub-fields positionable within a 7.2' patrol field, each sub-field producing a spectrum with a 14$\times$14-pixel spatial resolution. The main science drivers for KMOS are the study of galaxies, star formation, and molecular clouds, but its ability to simultaneously measure spectra of multiple stars makes KMOS an interesting instrument for exoplanet atmosphere characterization via transmission spectroscopy. We set to test whether transmission spectroscopy is practical with KMOS, and what are the conditions required to achieve the photometric precision needed, based on observations of a partial transit of WASP-19b, and full transits of GJ 1214b and HD 209458b. Our analysis uses the simultaneously observed comparison stars to reduce the effects from instrumental and atmospheric sources, and Gaussian processes to model the residual systematics. We show that KMOS can, in theory, deliver the photometric precision required for transmission spectroscopy. However, this is shown to require a) pre-imaging to ensure accurate centering and b) a very stable night with optimal observing conditions (seeing $\sim$0.8"). Combining these two factors with the need to observe several transits, each with a sufficient out-of-transit baseline (and with the fact that similar or better precision can be reached with telescopes and instruments with smaller pressure,) we conclude that transmission spectroscopy is not the optimal science case to take advantage of the abilities offered by KMOS and VLT.
Supernova remnants (SNRs) in the Galaxy are an important source of energy injection into the interstellar medium, and also of cosmic rays. Currently there are 294 known SNRs in the Galaxy, and their distribution with Galactocentric radius is of interest for various studies. Here I discuss some of the statistics of Galactic SNRs, including the observational selection effects that apply, and difficulties in obtaining distances for individual remnants from the `Sigma-D' relation. Comparison of the observed Galactic longitude distribution of a sample of bright Galactic SNRs -- which are not strongly affected by selection effects -- with those expected from models is used to constrain the Galactic distribution of SNRs. The best-fitting power-law/exponential model is more concentrated towards the Galactic centre than the widely used distribution obtained by Case & Bhattacharya (1998).
We report the detection of 48 millisecond pulsars (MSPs) out of 75 observed thus far using the LOFAR in the frequency range 110-188 MHz. We have also detected three MSPs out of nine observed in the frequency range 38-77 MHz. This is the largest sample of MSPs ever observed at these low frequencies, and half of the detected MSPs were observed for the first time at frequencies below 200 MHz. We present the average pulse profiles of the detected MSPs, their effective pulse widths and flux densities, and compare these with higher observing frequencies. The LOFAR pulse profiles will be publicly available via the EPN Database of Pulsar Profiles. We also present average values of dispersion measures (DM) and discuss DM and profile variations. About 35% of the MSPs show strong narrow profiles, another 25% exhibit scattered profiles, and the rest are only weakly detected. A qualitative comparison of the LOFAR MSP profiles with those at higher radio frequencies shows constant separation between profile components. Similarly, the profile widths are consistent with those observed at higher frequencies, unless scattering dominates at the lowest frequencies. This is very different from what is observed for normal pulsars and suggests a compact emission region in the MSP magnetosphere. The amplitude ratio of the profile components, on the other hand, can dramatically change towards low frequencies, often with the trailing component becoming dominant. As demonstrated by Dyks et al. (2010) this can be caused by aberration and retardation. This data set enables high-precision studies of pulse profile evolution with frequency, dispersion, Faraday rotation, and scattering in the interstellar medium. Characterizing and correcting these systematic effects may improve pulsar-timing precision at higher observing frequencies, where pulsar timing array projects aim to directly detect gravitational waves.
Based on our newly developed methods and the XMM-Newton large program of SN1006, we extract and analyze the spectra from 3596 tessellated regions of this SNR each with 0.3-8 keV counts $>10^4$. For the first time, we map out multiple physical parameters, such as the temperature ($kT$), electron density ($n_e$), ionization parameter ($n_et$), ionization age ($t_{ion}$), metal abundances, as well as the radio-to-X-ray slope ($\alpha$) and cutoff frequency ($\nu_{cutoff}$) of the synchrotron emission. We construct probability distribution functions of $kT$ and $n_et$, and model them with several Gaussians, in order to characterize the average thermal and ionization states of such an extended source. We construct equivalent width (EW) maps based on continuum interpolation with the spectral model of each regions. We then compare the EW maps of OVII, OVIII, OVII K$\delta-\zeta$, Ne, Mg, SiXIII, SiXIV, and S lines constructed with this method to those constructed with linear interpolation. We further extract spectra from larger regions to confirm the features revealed by parameter and EW maps, which are often not directly detectable on X-ray intensity images. For example, O abundance is consistent with solar across the SNR, except for a low-abundance hole in the center. This "O Hole" has enhanced OVII K$\delta-\zeta$ and Fe emissions, indicating recently reverse shocked ejecta, but also has the highest $n_et$, indicating forward shocked ISM. Therefore, a multi-temperature model is needed to decompose these components. The asymmetric metal distributions suggest there is either an asymmetric explosion of the SN or an asymmetric distribution of the ISM.
For line-driven winds from hot, luminous OB stars, we examine the subtle but
important role of diffuse, scattered radiation in determining both the topology
of steady-state solutions and intrinsic variability in the transonic wind base.
We use a smooth source function formalism to obtain nonlocal, integral
expressions for the direct and diffuse components of the line-force that
account for deviations from the usual localized, Sobolev forms. As the
scattering source function is reduced, we find the solution topology in the
transonic region transitions from X-type, with a unique wind solution, to a
nodal type, characterized by a degenerate family of solutions.
Specifically, in the idealized case of an optically thin source function and
a uniformly bright stellar disk, the unique X-type solution proves to be a
stable attractor to which time-dependent numerical radiation-hydrodynamical
simulations relax. But in models where the scattering strength is only modestly
reduced, the topology instead turns nodal, with the associated solution
degeneracy now manifest by intrinsic structure and variability that persist
down to the photospheric wind base. This highlights the potentially crucial
role of diffuse radiation for the dynamics and variability of line-driven
winds, and seriously challenges the use of Sobolev theory in the transonic wind
region. Since such Sobolev-based models are commonly used in broad applications
like stellar evolution and feedback, this prompts development of new wind
models, not only for further quantifying the intrinsic variability found here,
but also for computing new theoretical predictions of global properties like
velocity laws and mass-loss rates.
An active dipole antenna has been designed to measure transient electric field induced by ultra high energy cosmic rays for the CODALEMA experiment. The main requirements for this detector, composed of a low noise preamplifier placed close to a dipole antenna, are a wide bandwidth ranging from 100 kHz to 100 MHz and a good sensitivity on the whole spectrum.
We present a statistical framework for estimating global navigation satellite system (GNSS) non-ionospheric differential time delay bias. The biases are estimated by examining differences of measured line integrated electron densities (TEC) that are scaled to equivalent vertical integrated densities. The spatio-temporal variability, instrumentation dependent errors, and errors due to inaccurate ionospheric altitude profile assumptions are modeled as structure functions. These structure functions determine how the TEC differences are weighted in the linear least-squares minimization procedure, which is used to produce the bias estimates. A method for automatic detection and removal of outlier measurements that do not fit into a model of receiver bias is also described. The same statistical framework can be used for a single receiver station, but it also scales to a large global network of receivers. In addition to the Global Positioning System (GPS), the method is also applicable to other dual frequency GNSS systems, such as GLONASS (Globalnaya Navigazionnaya Sputnikovaya Sistema). The use of the framework is demonstrated in practice through several examples. A specific implementation of the methods presented here are used to compute GPS receiver biases for measurements in the MIT Haystack Madrigal distributed database system. Results of the new algorithm are compared with the current MIT Haystack Observatory MAPGPS bias determination algorithm. The new method is found to produce estimates of receiver bias that have reduced day-to-day variability and more consistent coincident vertical TEC values.
We determine accurate positions of the main satellites of Uranus: Miranda,
Ariel, Umbriel, Titania, and Oberon. Positions of Uranus, as derived from those
of these satellites, are also determined. The observational period spans from
1992 to 2011. All runs were made at the Pico dos Dias Observatory, Brazil.
We used the software called Platform for Reduction of Astronomical Images
Automatically (PRAIA) to minimise (digital coronography) the influence of the
scattered light of Uranus on the astrometric measurements and to determine
accurate positions of the main satellites. The positions of Uranus were then
indirectly determined by computing the mean differences between the observed
and ephemeris positions of these satellites. A series of numerical filters was
applied to filter out spurious data. These filters are mostly based on the
comparison between the positions of Oberon with those of the other satellites
and on the offsets as given by the differences between the observed and
ephemeris positions of all satellites.
We have, for the overall offsets of the five satellites, -29 (+/-63) mas in
right ascension and -27 (+/-46) mas in declination. For the overall difference
between the offsets of Oberon and those of the other satellites, we have +3
(+/-30) mas in right ascension and -2 (+/-28) mas in declination. Ephemeris
positions for the satellites were determined from DE432+ura111. Comparisons
using other modern ephemerides for the solar system -INPOP13c- and for the
motion of the satellites -NOE-7-2013- were also made. They confirm that the
largest contribution to the offsets we find comes from the motion of the
barycenter of the Uranus system around the barycenter of the solar system, as
given by the planetary ephemerides. Catalogues with the observed positions are
provided.
We present Rc-band surface photometry for 170 of the 203 galaxies in GHASP,
Gassendi H-Alpha survey of SPirals, a sample of late-type galaxies for which
high-resolution Fabry-Perot H{\alpha} maps have previously been obtained. Our
data set is constructed by new Rc-band observations taken at the Observatoire
de Haute-Provence (OHP), supplemented with Sloan Digital Sky Survey (SDSS)
archival data, obtained with the purpose of deriving homogeneous photometric
profiles and parameters. Our results include Rc-band surface brightness
profiles for 170 galaxies and $ugriz$ profiles for 108 of these objects. We
catalogue several parameters of general interest for further reference, such as
total magnitude, effective radius and isophotal parameters -- magnitude,
position angle, ellipticity and inclination. We also perform a structural
decomposition of the surface brightness profiles using a multi-component method
in order to separate disks from bulges and bars, and to observe the main
scaling relations involving luminosities, sizes and maximum velocities.
We determine the Rc-band Tully Fisher relation using maximum velocities
derived solely from H$\alpha$ rotation curves for a sample of 80 galaxies,
resulting in a slope of $-8.1 \pm 0.5$, zero point of $-3.0 \pm 1.0$ and an
estimated intrinsic scatter of $0.28 \pm 0.07$. We note that, different from
the TF-relation in the near-infrared derived for the same sample, no change in
the slope of the relation is seen at the low-mass end (for galaxies with
$V_{max} < 125$ km/s). We suggest that this different behaviour of the Tully
Fisher relation (with the optical relation being described by a single
power-law while the near-infrared by two) may be caused by differences in the
stellar mass to light ratio for galaxies with $V_{max} < 125$ km/s.
We explore the possibility to develop a new axion helioscope type, sensitive to the higher axion mass region favored by axion models. We propose to use a low background large volume TPC immersed in an intense magnetic field. Contrary to traditional tracking helioscopes, this detection technique takes advantage of the capability to directly detect the photons converted on the buffer gas which defines the axion mass sensitivity region, and does not require pointing the magnet to the Sun. The operation flexibility of a TPC to be used with different gas mixtures (He, Ne, Xe, etc) and pressures (from 10 mbar to 10 bar) will allow to enhance sensitivity for axion masses from few meV to several eV. We present different helioscope data taking scenarios, considering detection efficiency and axion absorption probability, and show the sensitivities reachable with this technique to be few $\times$ 10$^{-11}\,$GeV$^{-1}$ for a 5$\,$T$\,$m$^3$ scale TPC. We show that a few years program taking data with such setup would allow to probe the KSVZ axion model for axion masses above 100 meV.
The IceCube Collaboration initially reported the detection of 37 extraterrestrial neutrinos in the TeV - PeV energy range. The reconstructed neutrino events were obtained during three consecutive years of data taking, from 2010 to 2013. Although these events have been discussed to have an extragalactic origin, they have not correlated to any known source. Recently, the IceCube Collaboration reported a neutrino-induced muon event with energy of $2.6\pm0.3$ PeV which corresponds to the highest event ever detected. Neither the reconstructed direction of this event (J2000.0), detected on June 11 2014 at R.A.=110$^\circ$.34, Dec.=11$^\circ$.48 matches with any familiar source. Long gamma-ray bursts (lGRBs) are usually associated with the core collapse of massive stars leading a relativistic-collimated jets inside the star with high-energy neutrino production. These neutrinos have been linked to the 37 events previously detected by IceCube detector. In this work, we explore the conditions so that the highest neutrino recently detected could be generated by proton-photon and proton-hadron interactions in internal shocks inside GRB progenitor star and then be detected in IceCube detector. We found that not only low-luminosity but also high-luminosity GRBs could produce this neutrino in progenitors such as Wolf-Rayet and blue super giant stars.
We revisit the linear MRI in a cylindrical model of an accretion disk and uncover a number of attractive results overlooked in previous treatments. In particular, we elucidate the connection between local axisymmetric modes and global modes, and show that a local channel flow corresponds to the evanescent part of a global mode. In addition, we find that the global problem reproduces the local dispersion relation without approximation, a result that helps explain the success the local analysis enjoys in predicting global growth rates. MRI channel flows are nonlinear solutions to the governing equations in the local shearing box. However, only a small subset of MRI modes share the same property in global disk models, providing further evidence that the prominence of channels in local boxes is artificial. Finally, we verify our results via direct numerical simulations with the Godunov code RAMSES.
Extremely broad emission wings at H$\beta$ and H$\alpha$ have been found in VFTS data for five very luminous BA supergiants in or near 30 Doradus in the Large Magellanic Cloud. The profiles of both lines are extremely asymmetrical, which we have found to be caused by very broad diffuse interstellar bands (DIBs) in the longward wing of H$\beta$ and the shortward wing of H$\alpha$. These DIBs are well known to interstellar but not to many stellar specialists, so that the asymmetries may be mistaken for intrinsic features. The broad emission wings are generally ascribed to electron scattering, although we note difficulties for that interpretation in some objects. Such profiles are known in some Galactic hyper/supergiants and are also seen in both active and quiescent Luminous Blue Variables. No prior or current LBV activity is known in these 30 Dor stars, although a generic relationship to LBVs is not excluded; subject to further observational and theoretical investigation, it is possible that these very luminous supergiants are approaching the LBV stage for the first time. Their locations in the HRD and presumed evolutionary tracks are consistent with that possibility. The available evidence for spectroscopic variations of these objects is reviewed, while recent photometric monitoring does not reveal variability. A search for circumstellar nebulae has been conducted, with an indeterminate result for one of them.
Spectra of the cellular photospheric flows are determined from full-disk Doppler velocity observations acquired by the Helioseismic and Magnetic Imager (HMI) instrument on the Solar Dynamics Observatory (SDO) spacecraft. Three different analysis methods are used to separately determine spectral coefficients representing the poloidal flows, the toroidal flows, and the radial flows. The amplitudes of these spectral coefficients are constrained by simulated data analyzed with the same procedures as the HMI data. We find that the total velocity spectrum rises smoothly to a peak at a wavenumber of about 120 (wavelength of about 35 Mm), which is typical of supergranules. The spectrum levels off out to wavenumbers of about 400, and then rises again to a peak at a wavenumber of about 3500 (wavelength of about 1200 km), which is typical of granules. The velocity spectrum is dominated by the poloidal flow component (horizontal flows with divergence but no curl) at wavenumbers above 30. The toroidal flow component (horizontal flows with curl but no divergence) dominates at wavenumbers less than 30. The radial flow velocity is only about 3\% of the total flow velocity at the lowest wavenumbers, but increases in strength to become about 50\% at wavenumbers near 4000. The spectrum compares well with the spectrum of giant cell flows at the lowest wavenumbers and with the spectrum of granulation from a 3D radiative-hydrodynamic simulation at the highest wavenumbers.
We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95\% upper limit on the strain amplitude (at a frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times 10^{-15}$. For a broken power law model, we place priors on the strain amplitude derived from simulations of Sesana (2013) and McWilliams et al. (2014). We find that the data favor a broken power law to a pure power law with odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models, respectively. The McWilliams model is essentially ruled out by the data, and the Sesana model is in tension with the data under the assumption of a pure power law. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. We then place the most stringent limits so far on the energy density of relic GWs, $\Omega_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where $m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB, $\Omega_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a conservative limit of $G\mu<3.3\times 10^{-8}$ - a factor of 4 better than the joint Planck and high-$l$ CMB data from other experiments.
Dark Matter (DM) may be a thermal relic that annihilates into heavier states in the early Universe. This Forbidden DM framework accommodates a wide range of DM masses from keV to weak scales. An exponential hierarchy between the DM mass and the weak scale follows from the exponential suppression of the thermally averaged cross section. Stringent constraints from the cosmic microwave background are evaded because annihilations turn off at late times. We provide an example where DM annihilates into dark photons, which is testable through large DM self-interactions and direct detection.
We investigate the quantum mechanical transitions, induced by the combined effect of Gravitational wave (GW) and noncommutative (NC) structure of space, among the states of a 2-dimensional harmonic oscillator. The phonon modes excited by the passing GW within the resonant bar-detector are formally identical to forced harmonic oscillator and they represent a length variation of roughly the same order of magnitude as the characteristic length-scale of spatial noncommutativity estimated from the phenomenological upper bound of the NC parameter. This motivates our present work. We employ a number of different GW wave-forms that are typically expected from possible astronomical sources. We find that the transition probablities are quite sensitive to the nature of polarization of the GW. We further elaborate on the particular type of sources of GW radiation which can induce transitions that can be used as effective probe of the spatial noncommutative structure.
Links to: arXiv, form interface, find, astro-ph, recent, 1508, contact, help (Access key information)
Sugars of extraterrestrial origin have been observed in the interstellar medium (ISM), in at least one comet spectrum, and in several carbonaceous chondritic meteorites that have been recovered from the surface of the Earth. The origins of these sugars within the meteorites have been debated. To explore the possibility that sugars could be generated during shock events, this paper reports on the results of the first laboratory impact experiments wherein glycolaldehyde, found in the ISM, as well as glycolaldehyde mixed with montmorillonite clay, have been subjected to reverberated shocks from ~5 to >25 GPa. New biologically relevant molecules, including threose, erythrose and ethylene glycol, were identified in the resulting samples. These results show that sugar molecules can not only survive but also become more complex during impact delivery to planetary bodies.
The primary challenge of rocket propulsion is the burden of needing to accelerate the spacecraft's own fuel, resulting in only a logarithmic gain in maximum speed as propellant is added to the spacecraft. Light sails offer an attractive alternative in which fuel is not carried by the spacecraft, with acceleration being provided by an external source of light. By artificially illuminating the spacecraft with beamed radiation, speeds are only limited by the area of the sail, heat resistance of its material, and power use of the accelerating apparatus. In this paper, we show that leakage from a light sail propulsion apparatus in operation around a solar system analogue would be detectable. To demonstrate this, we model the launch and arrival of a microwave beam-driven light sail constructed for transit between planets in orbit around a single star, and find an optimal beam frequency on the order of tens of GHz. Leakage from these beams yields transients with flux densities of 0.1 Jy and durations of seconds at 100 pc. Because most travel within a planetary system would be conducted between the habitable worlds within that system, multiply-transiting exoplanetary systems offer the greatest chance of detection, especially when the planets are in projected conjunction as viewed from Earth. If interplanetary travel via beam-driven light sails is commonly employed in our galaxy, this activity could be revealed by radio follow-up of nearby transiting exoplanetary systems. The expected signal properties define a new strategy in the search for extraterrestrial intelligence (SETI).
The Two-Micron All-Sky Survey (2MASS) has mapped out the low-redshift Universe down to $K_S\sim14$ mag. As its near-infrared photometry primarily probes the featureless Rayleigh-Jeans tail of galaxy spectral energy distributions, colour-based redshift estimation is rather uninformative. Until now, redshift estimates for this dataset have relied on optical follow-up suffering from selection biases. Here we use the newly-developed technique of clustering-based redshift estimation to infer the redshift distribution of the 2MASS sources regardless of their optical properties. We characterise redshift distributions of objects from the Extended Source Catalog as a function of near-infrared colours and brightness and report some observed trends. We also apply the clustering redshift technique to dropout populations, sources with non-detections in one or more near-infrared bands, and present their redshift distributions. Combining all extended sources, we show that the redshift distribution of this sample extends up to $z\sim0.3$. We perform a similar analysis with the Point Source Catalog and show that it can be separated into stellar and extragalactic contributions with galaxies reaching $z\sim0.7$. We estimate that the Point Source Catalog contains 1.6 million extragalactic objects: as many as in the Extended Source Catalog but probing a cosmic volume ten times larger.
We spectroscopically survey the galaxy cluster XMM-LSS J02182-05102 (hereafter IRC 0218) using LRIS (optical) and MOSFIRE (near-infrared) on Keck I as part of the ZFIRE survey. IRC 0218 has a narrow redshift range of $1.612<z_{\rm spec}<1.635$ defined by 33 members of which 20 are at R$_{\rm proj}<1$ Mpc. The cluster redshift and velocity dispersion are $z_{\rm cl}=1.6233\pm0.0003$ and $\sigma_{\rm cl}=254\pm50$ km s$^{-1}$. We reach NIR line sensitivities of $\sim0.3\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$ that, combined with multi-wavelength photometry, provide extinction-corrected H$\alpha$ star formation rates (SFR), gas phase metallicities from [NII]/H$\alpha$, and stellar masses. We measure an integrated H$\alpha$ SFR of $\sim325{\rm M}_{\odot}$ yr$^{-1}$ (26 members; R$_{\rm proj}<2$ Mpc) and show that the elevated star formation in the cluster core (R$_{\rm proj}<0.25$ Mpc) is driven by the concentration of star-forming members, but the average SFR per H$\alpha$-detected galaxy is half that of members at R$_{\rm proj}\sim1$ Mpc. However, we do not detect any environmental imprint when comparing attenuation and gas phase metallicities: the cluster galaxies show similar trends with M$_{\star}$ as to the field, e.g. more massive galaxies have larger stellar attenuation. IRC 0218's gas phase metallicity-M$_{\star}$ relation (MZR) is offset to lower metallicities relative to $z\sim0$ and has a slope of $0.13\pm0.10$. Comparing the MZR in IRC 0218 to the COSMOS cluster at $z=2.1$ shows no evolution ($\Delta t\sim1$ Gyr): the MZR for both galaxy clusters are remarkably consistent with each other and virtually identical to several field surveys at $z\sim2$.
The Eclipse Mapping Method is an indirect imaging technique that transforms the shape of the eclipse light curve into a map of the surface brightness distribution of the occulted regions. Three decades of application of this technique to the investigation of the structure, the spectrum and the time evolution of accretion discs around white dwarfs in cataclysmic variables have enriched our understanding of these accretion devices with a wealth of details such as (but not limited to) moving heating/cooling waves during outbursts in dwarf novae, tidally-induced spiral shocks of emitting gas with sub-Keplerian velocities, elliptical precessing discs associated to superhumps, and measurements of the radial run of the disc viscosity through the mapping of the disc flickering sources. This chapter reviews the principles of the method, discusses its performance, limitations, useful error propagation procedures, as well as highlights a selection of applications aimed at showing the possible scientific problems that have been and may be addresses with it.
We describe the design and performance of the Medium Resolution Spectrometer (MRS) for the JWST-MIRI instrument. The MRS incorporates four coaxial spectral channels in a compact opto-mechanical layout that generates spectral images over fields of view up to 7.7 X 7.7 arcseconds in extent and at spectral resolving powers ranging from 1,300 to 3,700. Each channel includes an all-reflective integral field unit (IFU): an 'image slicer' that reformats the input field for presentation to a grating spectrometer. Two 1024 X 1024 focal plane arrays record the output spectral images with an instantaneous spectral coverage of approximately one third of the full wavelength range of each channel. The full 5 to 28.5 micron spectrum is then obtained by making three exposures using gratings and pass-band-determining filters that are selected using just two three-position mechanisms. The expected on-orbit optical performance is presented, based on testing of the MIRI Flight Model and including spectral and spatial coverage and resolution. The point spread function of the reconstructed images is shown to be diffraction limited and the optical transmission is shown to be consistent with the design expectations.
Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.
Long gamma-ray bursts (GRBs), among the most energetic events in the Universe, are explosions of massive and short-lived stars, so they pinpoint locations of recent star formation. However, several GRB host galaxies have recently been found to be deficient in molecular gas (H2), believed to be the fuel of star formation. Moreover, optical spectroscopy of GRB afterglows implies that the molecular phase constitutes only a small fraction of the gas along the GRB line-of-sight. Here we report the first ever 21 cm line observations of GRB host galaxies, using the Australia Telescope Compact Array, implying high levels of atomic hydrogen (HI), which suggests that the connection between atomic gas and star formation is stronger than previously thought, with star formation being potentially directly fuelled by atomic gas (or with very efficient HI-to-H2 conversion and rapid exhaustion of molecular gas), as has been theoretically shown to be possible. This can happen in low metallicity gas near the onset of star formation, because cooling of gas (necessary for star formation) is faster than the HI-to-H2 conversion. Indeed, large atomic gas reservoirs, together with low molecular gas masses, stellar and dust masses are consistent with GRB hosts being preferentially galaxies which have very recently started a star formation episode after accreting metal-poor gas from the intergalactic medium. This provides a natural route for forming GRBs in low-metallicity environments. The gas inflow scenario is also consistent with the existence of the companion HI object with no optical counterpart ~19 kpc from the GRB 060505 host, and with the fact that the HI centroids of the GRB 980425 and 060505 hosts do not coincide with optical centres of these galaxies, but are located close to the GRB positions.
First high-altitude observations of gravity wave (GW)-induced CO$_2$ density perturbations in the Martian thermosphere retrieved from NASA's NGIMS instrument on board the MAVEN satellite are presented and interpreted using the extended GW parameterization of Yi\u{g}it et al. [2008] and the Mars Climate Database as an input. Observed relative density perturbations between 180-220 km of 20-40 % demonstrate appreciable local time, latitude, and altitude variations. Modeling for the spatiotemporal conditions of the MAVEN observations suggests that GWs can directly propagate from the lower atmosphere to the thermosphere, produce appreciable dynamical effects, and likely contribute to the observed fluctuations. Modeled effects are somewhat smaller than the observed but their highly variable nature is in qualitative agreement with observations. Possible reasons for discrepancies between modeling and measurements are discussed.
We present two supernovae (SNe) discovered with the Hubble Space Telescope (HST) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), an HST multi-cycle treasury program. We classify both objects as Type Ia SNe and find redshifts of z = 1.80+-0.02 and 2.26 +0.02 -0.10, the latter of which is the highest redshift Type Ia SN yet seen. Using light curve fitting we determine luminosity distances and find that both objects are consistent with a standard Lambda-CDM cosmological model. These SNe were observed using the HST Wide Field Camera 3 infrared detector (WFC3-IR), with imaging in both wide- and medium-band filters. We demonstrate that the classification and redshift estimates are significantly improved by the inclusion of single-epoch medium-band observations. This medium-band imaging approximates a very low resolution spectrum (lambda/delta lambda ~ 100) which can isolate broad spectral absorption features that differentiate Type Ia SNe from their most common core collapse cousins. This medium-band method is also insensitive to dust extinction and (unlike grism spectroscopy) it is not affected by contamination from the SN host galaxy or other nearby sources. As such, it can provide a more efficient - though less precise - alternative to IR spectroscopy for high-z SNe.
The Cygnus OB2 Association is one of the nearest and largest collections of massive stars in the Galaxy. Situated at the heart of the "Cygnus X" complex of star-forming regions and molecular clouds, its distance has proven elusive owing to the ambiguous nature of kinematic distances along this $\ell\simeq80$ degree sightline and the heavy, patchy extinction. In an effort to refine the three-dimensional geometry of key Cygnus~X constituents, we have measured distances to four eclipsing double-lined OB-type spectroscopic binaries that are probable members of Cyg~OB2. We find distances of $1.33\pm0.17$, $1.32\pm0.07$, $1.44\pm0.18$, and $1.32\pm0.13$ kpc toward MT91~372, MT91~696, CPR2002~A36, and Schulte~3 respectively. We adopt a weighted average distance of 1.33$\pm$0.06~kpc. This agrees well with spectrophotometric estimates for the Association as a whole and with parallax measurements of protostellar masers in the surrounding interstellar clouds, thereby linking the ongoing star formation in these clouds with Cyg~OB2. We also identify Schulte 3C (O9.5V), a 4" visual companion to the 4.75 day binary Schulte~3(A+B), as a previously unrecognized Association member.
Magnetars are a special kind of neutron stars. There may also be accreting magnetars. From the studies of isolated magnetars, it is known that a neutron star with a strong dipole field only is not a magnetar. Super-slow X-ray pulsars may just be accreting high magnetic field neutron stars. The ultra-luminous X-ray pulsar NuSTAR J095551+6940.8 may be an accreting magnetar. It may be an accreting low magnetic field magnetar with multipole field of 10^14 G and dipole field of 10^12 G. This point of view is consistent with the study of isolated magnetars. An ultra-luminous X-ray pulsar phase in the binary evolution may result in massive millisecond pulsars.
The revolution of helio- and asteroseismology provides access to the detailed properties of stellar interiors by studying the star's oscillation modes. Among them, gravity (g) modes are formed by constructive interferences between progressive internal gravity waves (IGWs), propagating in stellar radiative zones. Our new 3D nonlinear simulations of the interior of a solar-like star allows us to study the excitation, propagation, and dissipation of these waves. The aim of this article is to clarify our understanding of the behavior of IGWs in a 3D radiative zone and to provide a clear overview of their properties. We use a method of frequency filtering that reveals the path of {individual} gravity waves of different frequencies in the radiative zone. We are able to identify the region of propagation of different waves in 2D and 3D, to compare them to the linear raytracing theory and to distinguish between propagative and standing waves (g modes). We also show that the energy carried by waves is distributed in different planes in the sphere, depending on their azimuthal wave number. We are able to isolate individual IGWs from a complex spectrum and to study their propagation in space and time. In particular, we highlight in this paper the necessity of studying the propagation of waves in 3D spherical geometry, since the distribution of their energy is not equipartitioned in the sphere.
r-Process nucleosynthesis in material ejected during neutron star mergers may lead to radioactively powered transients called kilonovae. The timescale and peak luminosity of these transients depend on the composition of the material after nuclear burning ceases, which determines the local heating rate from nuclear decays and the opacity. Kasen et al. (2013, ApJ, 774, 25) and Tanaka & Hotokezaka (2013, ApJ, 775, 113) pointed out that lanthanides can drastically increase the opacity in these outflows. We use the new general-purpose nuclear reaction network SkyNet to carry out a parameter study of r-process nucleosynthesis for a range of initial electron fractions $Y_e$, initial specific entropies $s$, and expansion timescales $\tau$. We find that the ejecta is lanthanide-free for $Y_e \gtrsim 0.22 - 0.30$, depending on $s$ and $\tau$. The heating rate is insensitive to $s$ and $\tau$, but certain, larger values of $Y_e$ lead to reduced heating rates, due to individual nuclides dominating the heating. With a simplified gray radiative transport scheme in spherical symmetry, we estimate the luminosity, time, and effective temperature at the peak of the light curves. We find that the luminosity peaks much earlier, at about a day in the lanthanide-free cases compared to a week in the lanthanide-rich cases. The heating rate does not change much as the ejecta becomes lanthanide-free with increasing $Y_e$, but the light curve peak becomes about an order of magnitude brighter because it peaks much earlier when the heating rate is larger. We also provide parametric fits for the heating rates between 0.1 and 100 days, and we provide a simple fit in $Y_e$, $s$, and $\tau$ to estimate whether the ejecta is lanthanide-rich or not.
Using a three-dimensional magnetohydrodynamic model, we simulate the magnetic reconnection in a single current sheet. We assume a finite guide field, a random perturbation on the velocity field and uniform resistivity. Our model enhances the reconnection rate relative to the classical Sweet-Parker model in the same configuration. The efficiency of magnetic energy conversion is increased by interactions between the multiple tearing layers coexisting in the global current sheet. This interaction, which forms a positive-feedback system, arises from coupling of the inflow and outflow regions in different layers across the current sheet. The coupling accelerates the elementary reconnection events, thereby enhancing the global reconnection rate. The reconnection establishes flux tubes along each tearing layer. Slow-mode shocks gradually form along the outer boundaries of these tubes, further accelerating the magnetic energy conversion.Such positive-feedback system is absent in two-dimensional simulation, three-dimensional reconnection without a guide field and a reconnection under a single perturbation mode. We refer to our model as the "shock-evoking positive-feedback" model.
Can a white dwarf, accreting hydrogen-rich matter from a non-degenerate companion star, ever exceed the Chandrasekhar mass and explode as a type Ia supernova? We explore the range of accretion rates that allow a white dwarf (WD) to secularly grow in mass, and derive limits on the accretion rate and on the initial mass that will allow it to reach $1.4M_\odot$ --- the Chandrasekhar mass. We follow the evolution through a long series of hydrogen flashes, during which a thick helium shell accumulates. This determines the effective helium mass accretion rate for long-term, self-consistent evolutionary runs with helium flashes. We find that net mass accumulation always occurs despite helium flashes. Although the amount of mass lost during the first few helium shell flashes is a significant fraction of that accumulated prior to the flash, that fraction decreases with repeated helium shell flashes. Eventually no mass is ejected at all during subsequent flashes. This unexpected result occurs because of continual heating of the WD interior by the helium shell flashes near its surface. The effect of heating is to lower the electron degeneracy throughout the WD, and especially in the outer layers. This key result yields helium burning that is quasi-steady state, instead of explosive. We thus find a remarkably large parameter space within which long-term, self-consistent simulations show that a WD can grow in mass and reach the Chandrasekhar limit, despite its helium flashes.
A tidal radius is a distance from a satellite orbiting in a host potential beyond which its material is stripped by the tidal force. We derive a revised expression for the tidal radius of a rotating satellite which properly takes into account the possibility of prograde and retrograde orbits of stars. Besides the eccentricity of the satellite orbit, the tidal radius depends also on the ratio of the satellite internal angular velocity to the orbital angular velocity. We compare our formula to the results of two $N$-body simulations of dwarf galaxies orbiting a Milky Way-like host on a prograde and retrograde orbit. The tidal radius for the retrograde case is larger than for the prograde. We introduce a kinematic radius separating stars still orbiting the dwarf galaxy from those already stripped and following the potential of the host galaxy. We find that the tidal radius matches very well the kinematic radius.
We present the third data release from the Australia Telescope Large Area Survey (ATLAS). These data combine the observations at 1.4 GHz before and after upgrades to the Australia Telescope Compact Array reaching a sensitivity of 14 microJy/beam in 3.6 deg^2 over the Chandra Deep Field South (CDFS) and of 17 microJy/beam in 2.7 deg^2 over the European Large Area ISO Survey South 1 (ELAIS-S1). We used a variety of array configurations to maximise the uv coverage resulting in a resolution of 16 by 7 arcsec in CDFS and of 12 by 8 arcsec in ELAIS-S1. After correcting for peak bias and bandwidth smearing, we find a total of 3034 radio source components above 5 sigma in CDFS, of which 514 (17 per cent) are considered to be extended. The number of components detected above 5 sigma in ELAIS-S1 is 2084, of which 392 (19 per cent) are classified as extended. The catalogues include reliable spectral indices (delta alpha < 0.2) between 1.40 and 1.71 GHz for ~350 of the brightest components.
We studied the outburst properties of the hyper-luminous X-ray source ESO 243-49 HLX-1, using the full set of Swift monitoring observations. We quantified the increase in the waiting time, recurrence time and e-folding rise timescale along the outburst sequence, and the corresponding decrease in outburst duration, total radiated energy, and e-folding decay timescale, which confirms previous findings. HLX-1 spends less and less time in outburst and more and more time in quiescence, but its peak luminosity remains approximately constant. We compared the HLX-1 outburst properties with those of bright Galactic low-mass X-ray binary transients (LMXBTs). Our spectral analysis strengthens the similarity between state transitions in HLX-1 and those in Galactic LMXBTs. We also found that HLX-1 follows the nearly linear correlations between the hard-to-soft state transition luminosity and the peak luminosity, and between the rate of change of X-ray luminosity during the rise phase and the peak luminosity, which indicates that the occurrence of the hard-to-soft state transition of HLX-1 is similar to those of Galactic LMXBTs during outbursts. We found that HLX-1 does not follow the correlations between total radiated energy and peak luminosity, and between total radiated energy and e-folding rise/decay timescales we had previously identified in Galactic LMXBTs. HLX-1 would follow those correlations if the distance were several hundreds of kpc. But invoking a much closer distance for HLX-1 is not a viable solution to this problem, as it introduces other more serious inconsistencies with the observations.
The estimation of cosmological constraints from observations of the large scale structure of the Universe, such as the power spectrum or the correlation function, requires the knowledge of the inverse of the associated covariance matrix, namely the precision matrix, $\mathbf{\Psi}$. In most analyses, $\mathbf{\Psi}$ is estimated from a limited set of mock catalogues. Depending on how many mocks are used, this estimation has an associated error which must be propagated into the final cosmological constraints. For future surveys such as Euclid and DESI, the control of this additional uncertainty requires a prohibitively large number of mock catalogues. In this work we test a novel technique for the estimation of the precision matrix, the covariance tapering method, in the context of baryon acoustic oscillation measurements. Even though this technique was originally devised as a way to speed up maximum likelihood estimations, our results show that it also reduces the impact of noisy precision matrix estimates on the derived confidence intervals, without introducing biases on the target parameters. The application of this technique can help future surveys to reach their true constraining power using a significantly smaller number of mock catalogues.
Hereafter we describe the activities of the $Grand \, Sud-Ouest$ Data Centre operated for INSU/CNRS by the OMP-IRAP and the Universit\'e Paul Sabatier (Toulouse), in a collaboration with the OASU-LAB (Bordeaux) and OREME-LUPM (Montpellier).
Current Solar Energetic Particle (SEP) propagation models describe the effects of IMF turbulence on SEPs as diffusion, using a Fokker-Planck (FP) equation. However, these models cannot explain the observed fast access of SEPs to regions widely separated in longitude within the heliosphere, across the average magnetic field, without use of unrealistically strong cross-field diffusion. Here, we show that accounting for the SEP propagation along field-lines that meander due to plasma turbulence provides an explanation for the wide SEP events. Using a model that includes the non-diffusive effects caused by field-line meandering, we reproduce the observed longitudinal extent of SEP peak fluxes characterised by a Gaussian profile with sigma=30-50 degrees, while current diffusion theory can only explain extents of 10 degrees with realistic diffusion coefficients. Our model also reproduces the timing of SEP arrival at distant longitudes, which cannot be explained using the diffusion model.
The Large Synoptic Survey Telescope (LSST) will be the largest time-domain photometric survey ever. In order to maximize the LSST science yield for a broad array of transient stellar phenomena, it is necessary to optimize the survey cadence, coverage, and depth via quantitative metrics that are specifically designed to characterize the time-domain behavior of various types of stellar transients. In this paper we present three such metrics built on the LSST Metric Analysis Framework (MAF) model (Jones et al. 2014). Two of the metrics quantify the ability of LSST to detect non-periodic and/or non-recurring transient events, and the ability of LSST to reliably measure periodic signals of various timescales. The third metric provides a way to quantify the range of stellar parameters in the stellar populations that LSST will probe. We provide example uses of these metrics and discuss some implications based on these metrics for optimization of the LSST survey for observations of stellar variables and transients.
Many spiral galaxy haloes show stellar streams with various morphologies when observed with deep images. The origin of these tidal features is discussed, either coming from a satellite infall or caused by residuals of an ancient, gas-rich major merger. By modelling the formation of the peculiar features observed in the NGC 4013 halo, we investigate their origin. By using GADGET -2 with implemented gas cooling, star formation, and feedback, we have modelled the overall NGC 4013 galaxy and its associated halo features. A gas-rich major merger occurring 2.7-4.6 Gyr ago succeeds in reproducing the NGC 4013 galaxy properties, including all the faint stellar features, strong gas warp, boxy-shaped halo and vertical 3.6 mum luminosity distribution. High gas fractions in the progenitors are sufficient to reproduce the observed thin and thick discs, with a small bulge fraction, as observed. A major merger is able to reproduce the overall NGC 4013 system, including the warp strength, the red colour and the high stellar mass density of the loop, while a minor merger model cannot. Because the gas-rich model suffices to create a pseudo-bulge with a small fraction of the light, NGC 4013 is perhaps the archetype of a late-type galaxy formed by a relatively recent merger. Then late type, pseudo-bulge spirals are not mandatorily made through secular evolution, and the NGC 4013 properties also illustrate that strong warps in isolated galaxies may well occur at a late phase of a gas-rich major merger.
The mass-loss process in Mira stars probably occurs in an asymmetric way where dust can form in inhomogeneous circumstellar molecular clumps. Following asymmetries along the pulsation cycle can give us clues about these mass-loss processes. We imaged the Mira star X Hya and its environnement at different epochs to follow the evolution of the morphology in the continuum and in the molecular bands. We observed X Hya with AMBER in J-H-K at low resolution at two epochs. We modelled squared visibilities with geometrical and physical models. We also present imaging reconstruction results obtained with MiRA and based on the physical a priori images. We report on the angular scale change of X Hya between the two epochs. 1D CODEX profiles allowed us to understand and model the spectral variation of squared visibilities and constrain the stellar parameters. Reconstructed model-dependent images enabled us to reproduce closure phase signals and the azimuthal dependence of squared visibilities. They show evidence for material inhomogeneities located in the immediate environment of the star.
Our aim is to precisely measure the physical parameters of the eclipsing binary IO Aqr and derive a distance to this system by applying a surface brightness - colour relation. Our motivation is to combine these parameters with future precise distance determinations from the GAIA space mission to derive precise surface brightness - colour relations for stars. We extensively used photometry from the Super-WASP and ASAS projects and precise radial velocities obtained from HARPS and CORALIE high-resolution spectra. We analysed light curves with the code JKTEBOP and radial velocity curves with the Wilson-Devinney program. We found that IO Aqr is a hierarchical triple system consisting of a double-lined short-period (P=2.37 d) spectroscopic binary and a low-luminosity and low-mass companion star orbiting the binary with a period of ~25000 d (~70 yr) on a very eccentric orbit. We derive high-precision (better than 1%) physical parameters of the inner binary, which is composed of two slightly evolved main-sequence stars (F5 V-IV + F6 V-IV) with masses of M1=1.569+/-0.004 and M2=1.655+/-0.004 M_sun and radii R1=2.19+/-0.02 and R2=2.49+/-0.02 R_sun. The companion is most probably a late K-type dwarf with mass ~0.6 M_sun. The distance to the system resulting from applying a (V-K) surface brightness - colour relation is 255+/-6(stat.)+/-6(sys.) pc, which agrees well with the Hipparcos value of 270+/-73 pc, but is more precise by a factor of eight.
Accreting white dwarfs (WDs) with non-degenerate companions are expected to emit in soft X-rays and the UV, if accreted H-rich material burns stably. They are an important component of the unresolved emission of elliptical galaxies, and their combined ionizing luminosity may significantly influence the optical line emission from warm ISM. In an earlier paper we modeled populations of accreting WDs, first generating WD with main-sequence, Hertzsprung gap and red giant companions with the population synthesis code \textsc{BSE}, and then following their evolution with a grid of evolutionary tracks computed with \textsc{MESA}. Now we use these results to estimate the soft X-ray (0.3-0.7keV), H- and He II-ionizing luminosities of nuclear burning WDs and the number of super-soft X-ray sources for galaxies with different star formation histories. For the starburst case, these quantities peak at $\sim 1$ Gyr and decline by $\sim 1-3$ orders of magnitude by the age of 10 Gyr. For stellar ages of $\sim$~10 Gyr, predictions of our model are consistent with soft X-ray luminosities observed by Chandra in nearby elliptical galaxies and He II 4686$\AA/\rm{H}{\beta}$ line ratio measured in stacked SDSS spectra of retired galaxies, the latter characterising the strength and hardness of the UV radiation field. However, the soft X-ray luminosity and He~II~4686$\AA/\rm{H}{\beta}$ ratio are significantly overpredicted for stellar ages of $\lesssim 4-8$ Gyr. We discuss various possibilities to resolve this discrepancy and tentatively conclude that it may be resolved by a modification of the typically used criteria of dynamically unstable mass loss for giant stars.
We present the full catalog of Young Stellar Objects (YSOs) identified in the 18 molecular clouds surveyed by the Spitzer Space Telescope "cores to disks" (c2d) and "Gould Belt" (GB) Legacy surveys. Using standard techniques developed by the c2d project, we identify 3239 candidate YSOs in the 18 clouds, 2966 of which survive visual inspection and form our final catalog of YSOs in the Gould Belt. We compile extinction corrected SEDs for all 2966 YSOs and calculate and tabulate the infrared spectral index, bolometric luminosity, and bolometric temperature for each object. We find that 326 (11%), 210 (7%), 1248 (42%), and 1182 (40%) are classified as Class 0+I, Flat-spectrum, Class II, and Class III, respectively, and show that the Class III sample suffers from an overall contamination rate by background AGB stars between 25% and 90%. Adopting standard assumptions, we derive durations of 0.40-0.78 Myr for Class 0+I YSOs and 0.26-0.50 Myr for Flat-spectrum YSOs, where the ranges encompass uncertainties in the adopted assumptions. Including information from (sub)millimeter wavelengths, one-third of the Class 0+I sample is classified as Class 0, leading to durations of 0.13-0.26 Myr (Class 0) and 0.27-0.52 Myr (Class I). We revisit infrared color-color diagrams used in the literature to classify YSOs and propose minor revisions to classification boundaries in these diagrams. Finally, we show that the bolometric temperature is a poor discriminator between Class II and Class III YSOs.
We report the results of a search for large velocity width, low-intensity line wings - a commonly used signature of molecular outflows - in four low redshift (ultra)luminous infrared galaxies (U/LIRGs) that appear to be dominated by star formation. The targets were drawn from a sample of fourteen such galaxies presented in Chung et al. (2011), who showed the stacked CO spectrum of the sample to exhibit 1000 km/s-wide line wings. We obtained sensitive, wide bandwidth imaging of our targets using the IRAM Plateau de Bure Interferometer. We detect each target at very high significance but do not find the claimed line wings in these four targets. Instead, we constrain the flux in the line wings to be only a few percent. Casting our results as mass outflow rates following Cicone et al. (2014) we show them to be consistent with a picture in which very high mass loading factors preferentially occur in systems with high AGN contributions to their bolometric luminosity. We identify one of our targets, IRAS05083 (VII Zw 31), as a candidate molecular outflow.
A complete, flux density limited sample of 96 faint ($> 0.5$ mJy) radio sources is selected from the 10C survey at 15.7 GHz in the Lockman Hole. We have matched this sample to a range of multi-wavelength catalogues, including SERVS, SWIRE, UKIDSS and optical data; multi-wavelength counterparts are found for 80 of the 96 sources and spectroscopic redshifts are available for 24 sources. Photometric reshifts are estimated for the sources with multi-wavelength data available; the median redshift of the sample is 0.91 with an interquartile range of 0.84. Radio-to-optical ratios show that at least 94 per cent of the sample are radio loud, indicating that the 10C sample is dominated by radio galaxies. This is in contrast to samples selected at lower frequencies, where radio-quiet AGN and starforming galaxies are present in significant numbers at these flux density levels. All six radio-quiet sources have rising radio spectra, suggesting that they are dominated by AGN emission. These results confirm the conclusions of Paper I that the faint, flat-spectrum sources which are found to dominate the 10C sample below $\sim 1$ mJy are the cores of radio galaxies. The properties of the 10C sample are compared to the SKADS Simulated Skies; a population of low-redshift starforming galaxies predicted by the simulation is not found in the observed sample.
We investigate the migration of low-mass planets ($5 M_{\oplus}$ and $20 M_{\oplus}$) in accretion discs threaded with a magnetic field using 2D MHD code in polar coordinates. We observed that, in the case of a strong azimuthal magnetic field where the plasma parameter is $\beta\sim 1-2$, density waves at the magnetic resonances exert a positive torque on the planet and may slow down or reverse its migration. However, when the magnetic field is weaker (i.e., the plasma parameter $\beta$ is relatively large), then non-axisymmetric density waves excited by the planet lead to growth of the radial component of the field and, subsequently, to development of the magneto-rotational instability, such that the disc becomes turbulent. Migration in a turbulent disc is stochastic, and the migration direction may change as such. To understand migration in a turbulent disc, both the interaction between a planet and individual turbulent cells, as well as the interaction between a planet and ordered density waves, have been investigated.
The Ulysses spacecraft provided the first opportunity to identify and study Interstellar Dust (ISD) in-situ in the Solar System between 1992 and 2007. Here we present the first comprehensive analysis of the ISD component in the entire Ulysses dust data set. We analysed several parameters of the ISD flow in a time-resolved fashion: flux, flow direction, mass index, and flow width. The general picture is in agreement with a time-dependent focussing/defocussing of the charged dust particles due to long-term variations of the solar magnetic field throughout a solar magnetic cycle of 22 years. In addition, we confirm a shift in dust direction of $50^{\circ} \pm 7^{\circ}$ in 2005, along with a steep, size-dependent increase in flux by a factor of 4 within 8 months. To date, this is difficult to interpret and has to be examined in more detail by new dynamical simulations. This work is part of a series of three papers. This paper concentrates on the time-dependent flux and direction of the ISD. In a companion paper (Kr\"uger et al., 2015) we analyse the overall mass distribution of the ISD measured by Ulysses, and a third paper discusses the results of modelling the flow of the ISD as seen by Ulysses (Sterken et al., 2015).
We present a new observable, position-dependent power spectrum, to measure the large-scale structure bispectrum in the squeezed configuration, where one wavenumber is much smaller than the other two. The squeezed-limit bispectrum measures how the small-scale power spectrum is modulated by a long-wavelength overdensity, which is due to gravitational evolution and possibly inflationary physics. We divide a survey into small subvolumes, compute the local power spectrum and the mean overdensity in each subvolume, and measure the correlation between them. The correlation measures the integral of the bispectrum, which is dominated by squeezed configurations if the scale of the local power spectrum is much smaller than the subvolume size. We use the separate universe approach to model how the small-scale power spectrum is affected by a long-wavelength overdensity gravitationally. This models the nonlinearity of the bispectrum better than the perturbation theory approach. Not only the new observable is easy to interpret, but it sidesteps the complexity of the full bispectrum estimation as both power spectrum and mean overdensity are easier to estimate than the full bispectrum. We report on the first measurement of the position-dependent correlation function from the SDSS-III BOSS DR10 CMASS sample. We detect the bispectrum of the CMASS sample, and constrain their nonlinear bias combining with anisotropic clustering and weak lensing. We finally study the response of the small-scale power spectrum to 1-3 long-wavelength overdensities. We compare the separate universe approach to separate universe simulations to unprecedented accuracy. We test the standard perturbation theory (SPT) hypothesis that the nonlinear n-point function is fully predicted by the linear power spectrum at the same time. We find discrepancies on small scales, which suggest that SPT fails even if it is calculated to all orders.
We report on deep XMM-Newton observations of the vertex filament in the southern giant lobe of the Fanaroff-Riley class I radio galaxy Centaurus A. We find no X-ray excess from the filament region and place a 3 sigma upper limit on the 1 keV flux density of the filament of 9.6 nJy. This directly constrains the electron density and magnetic field strength in the filament. For the first time in an individual filament, we show that the excess in synchrotron emissivity cannot be produced purely by excess electrons: the filament magnetic field strength must be higher than in the giant lobes as a whole, and close to or above the equipartition value for the filament. The filaments are not significantly overpressured with respect to the surrounding lobe with a pressure provided by relativistic electrons.
We investigate whether the far-UV continuum of nearby radio galaxies reveals
evidence for the presence of star forming or non-stellar components. If a UV
excess due to an extra radiation component exists we compare this with other
properties such as radio power, optical spectral type and the strength of the
emission lines. We also discuss the possible correlation between the
ultra-violet flux, IR properties and central black hole mass. We use two sampes
of low luminosity radio galaxies with comparable redshifts ($z < 0.2$).
Spectral Energy Distributions are constructed using a number of on-line
databases: GALEX, SDSS, 2MASS, and WISE. The parameter $XUV$ is introduced,
which measures the excess slope of the UV continuum between 4500 and 2000 \AA,
with respect to the UV radiation produced by the underlying old galaxy
component.
We find that the UV excess is usually small or absent in low luminosity
sources, but sets in abruptly at the transition radio power above which we find
mostly FRII sources. $XUV$ behaves very similarly to the strength of the
optical emission lines (in particular $H\alpha$). Below $P_{1.4 GHz} < 10^{24}$
WHz$^{-1}$ $XUV$ is close to zero. $XUV$ correlates strongly with the $H\alpha$
line strength, but only in sources with strong $H\alpha$ emission. There is a
strong correlation between $XUV$ and the slope of the mid-IR, as measured by
the WISE bands in the interval 3.4 to 22 $\mu$m, in the sense that sources with
a strong UV excess also have stronger IR emission. There is an inverse
correlation between $XUV$ and central black hole mass: strong UV excess objects
have, on average, $M_{BH}$ about 2-3 times less massive than those without UV
excess. Low luminosity radio galaxies tend to be more massive and contain more
massive black holes.
An analysis of absorption profiles of the $\lambda$6614 diffuse interstellar band recorded along the lines-of-sight towards HD 179406 (20 Aql) and HD 147889 is described. The difference in band shape is attributed to the degree of internal excitation of the carrier, which is principally due to vibrational hot bands although an electronic component may also be present. The results are discussed with respect to other models of diffuse band spectral line shape.
As my own work on the Sun's magnetic field started exactly 50 years ago at Crimea in the USSR, I have been a participant in the field during nearly half the time span since Hale's discovery in 1908 of magnetic fields in sunspots. The present historical account is accompanied by photos from my personal slide collection, which show a number of the leading personalities who advanced the field in different areas: measurement techniques, from photographic to photoelectric and imaging methods in spectro-polarimetry; theoretical foundations of MHD and the origin of cosmic magnetic fields (birth of dynamo theory); the quest for increased angular resolution from national projects to international consortia (for instruments both on ground and in space); introduction of the Hanle effect in astrophysics and the Second Solar Spectrum as its playground; small-scale nature of the field, the fundamental resolution limit, and transcending it by resolution-independent diagnostics.
It has been proposed that a galaxy's nova rate might be enhanced by the production of nova progenitor binaries in the dense cores of its globular clusters (GCs). To explore this idea, relative nova rates in three Virgo elliptical galaxies, M87, M49 and M84, which have significantly different GC specific frequencies ($S_{N}$) of 14, 3.6, and 1.6, respectively, were measured over the course of 4 epochs spanning a period of 14 months. To simplify the analysis, observations of the nearly equidistant galaxies were made on the same nights, with the same integration times, and through the same filter (H$\alpha$), so that the relative numbers of novae discovered would reflect the relative nova rates. At the conclusion of our survey we found a total of 27 novae associated with M87, 37 with M49, and 19 with M84. After correcting for survey completeness, we found annual nova rates of $154^{+23}_{-19}$, $189^{+26}_{-22}$, and $95^{+15}_{-14}$, for M87, M49, and M84, respectively, corresponding to $K$-band luminosity-specific nova rates of $3.8\pm1.0$, $3.4\pm0.6$, and $3.0\pm0.6$ novae per year per $10^{10}~L_{K,\odot}$. The overall results of our study suggest that a galaxy's nova rate simply scales with its luminosity, and is insensitive to its GC specific frequency. Two novae, one in M87 and one in M84, were found to be spatially coincident with known GCs. After correcting for the mass fraction in GCs, we estimate that novae are likely enhanced relative to the field by at least an order of magnitude in the GC systems of luminous Virgo ellipticals.
List of proceedings from the HAWC Collaboration presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.
We study black hole solutions at first order in the Hartle-Thorne slow-rotation approximation in Horndeski gravity theories. We derive the equations of motion including also cases where the scalar depends linearly on time. In the Hartle-Thorne formalism, all first-order rotational corrections are described by a single frame-dragging function. We show that the frame dragging function is exactly the same as in general relativity for all known BH solutions in shift-symmetric Horndeski theories, with the exception of theories with a linear coupling to the Gauss-Bonnet invariant. Our results extend previous no-hair theorems for a broad class of Horndeski gravity theories.
Dark Matter particles with inelastic interactions are ubiquitous in extensions of the Standard Model, yet remain challenging to fully probe with existing strategies. We propose a series of powerful searches at hadron and lepton colliders that are sensitive to inelastic dark matter dynamics. In representative models, we find that the LHC and BaBar could offer strong sensitivity to the thermal-relic dark matter parameter space for dark matter masses between ~100 MeV-100 GeV and fractional mass-splittings above the percent level; future searches at Belle II with a dedicated monophoton trigger could also offer sensitivity to thermal-relic scenarios with masses below a few GeV. Thermal scenarios with either larger masses or splittings are largely ruled out; lower masses remain viable yet may be accessible with other search strategies.
An interesting cosmological history was proposed by Prigogine {\it et al.} who considered the Universe as a thermodynamically open system. This scenario is characterized by a process of matter creation, which corresponds to an irreversible energy flow from the gravitational field to the pressureless matter fluid. Here, we show that the gravitationally induced particle production may arise from a nonminimal curvature-matter coupling. By considering the equivalent scalar-tensor theory, the cosmological implications of the model are discussed. As all known natural systems tend to a state of thermodynamic equilibrium, and assuming the universe is not different in this respect, we also discuss the conditions to attain the equilibrium state.
We develop a novel sampling theorem for functions defined on the three-dimensional rotation group SO(3) by associating the rotation group with the three-torus through a periodic extension. Our sampling theorem requires $4L^3$ samples to capture all of the information content of a signal band-limited at $L$, reducing the number of required samples by a factor of two compared to other equiangular sampling theorems. We present fast algorithms to compute the associated Fourier transform on the rotation group, the so-called Wigner transform, which scale as $O(L^4)$, compared to the naive scaling of $O(L^6)$. For the common case of a low directional band-limit $N$, complexity is reduced to $O(N L^3)$. Our fast algorithms will be of direct use in speeding up the computation of directional wavelet transforms on the sphere. We make our SO3 code implementing these algorithms publicly available.
We systematically examine the properties of null geodesics around an electrically charged, asymptotically flat black hole in Eddington-inspired Born-Infeld gravity, varying the electric charge of black hole and the coupling constant in the theory. We find that the radius of the unstable circular orbit for massless particle decreases with the coupling constant, if the value of the electrical charge is fixed. Additionally, we consider the strong gravitational lensing around such a black hole. We show that the deflection angle, the position angle of the relativistic images, and the magnification due to the light bending in strong gravitational field are quite sensitive to the parameters determining the black hole solution. Thus, through the accurate observations associated with the strong gravitational lensing, it might be possible to reveal the gravitational theory in a strong field regime.
We examine the spatial extrema (local maxima, minima and saddle points) of the covariant scalars (density, Hubble expansion, spatial curvature and eigenvalues of the shear and electric Weyl tensors) of the quasi-spherical Szekeres dust models. Sufficient conditions are obtained for the existence of distributions of multiple extrema in spatial comoving locations that can be prescribed through initial conditions. These distributions evolve without shell crossing singularities at least for ever expanding models (with or without cosmological constant) in the full evolution range where the models are valid. By considering the local maxima and minima of the density, our results allow for setting up elaborated networks of "pancake" shaped evolving cold dark matter over-densities and density voids whose spatial distribution and amplitudes can be controlled from initial data compatible with standard early Universe initial conditions. We believe that these results have an enormous range of potential application to provide fully relativistic non-perturbative modelling of cosmic structure at all scales.
Pitch-angle scattering rates for cosmic-ray particles in magnetohydrodynamic (MHD) simulations with imbalanced turbulence are calculated for fully evolving electromagnetic turbulence. We compare with theoretical predictions derived from the quasilinear theory of cosmic-ray diffusion for an idealized slab spectrum and demonstrate how cross helicity affects the shape of the pitch-angle diffusion coefficient. Additional simulations in evolving magnetic fields or static field configurations provide evidence that the scattering anisotropy in imbalanced turbulence is not primarily due to coherence with propagating Alfven waves, but an effect of the spatial structure of electric fields in cross-helical MHD turbulence.
We derive the relation between the amplitudes of short-baseline appearance and disappearance oscillations in 3+$N_{s}$ neutrino mixing schemes which is the origin of the appearance-disappearance tension that is found from the analysis of the existing data in any 3+$N_{s}$ neutrino mixing scheme. We illustrate the power of the relation to reveal the appearance-disappearance tension in the cases of 3+1 and 3+2 mixing using the results of a global fit of short-baseline neutrino oscillation data.
Loop quantum cosmology (LQC) provides promising resolutions to the trans-Planckian issue and initial singularity arising in the inflationary models of general relativity. In general, due to different quantization approaches, LQC involves two types of quantum corrections, the holonomy and inverse-volume, to both of the cosmological background evolution and perturbations. In this paper, using {\em the third-order uniform asymptotic approximations}, we derive explicitly the observational quantities of the slow-roll inflation in the framework of LQC with these quantum corrections. We calculate the power spectra, spectral indices, and running of the spectral indices for both scalar and tensor perturbations, whereby the tensor-to-scalar ratio is obtained. We expand all the observables at the time when the inflationary mode crosses the Hubble horizon. As the upper error bounds for the uniform asymptotic approximation at the third-order are $\lesssim 0.15\%$, these results represent the most accurate results obtained so far in the literature. It is also shown that with the inverse-volume corrections, both scalar and tensor spectra exhibit a deviation from the usual shape at large scales. Then, using the Planck, BAO and SN data we obtain new constraints on quantum gravitational effects from LQC corrections, and find that such effects could be within the detection of the forthcoming experiments.
The celebrated Weinberg theorem in cosmological perturbation theory states that there always exist two adiabatic scalar modes in which the comoving curvature perturbation is conserved on super-horizon scales. In particular, when the perturbations are generated from a single source, such as in single field models of inflation, both of the two allowed independent solutions are adiabatic and conserved on super-horizon scales. There are few known examples in literature which violate this theorem. We revisit the theorem and specify the loopholes in some technical assumptions which violate the theorem in models of non-attractor inflation, fluid inflation, solid inflation and in the model of pseudo conformal universe.
Links to: arXiv, form interface, find, astro-ph, recent, 1508, contact, help (Access key information)