We report the discovery of Lyman-alpha emission (Ly$\alpha$) in the bright galaxy EGSY-2008532660 (hereafter EGSY8p7) using the MOSFIRE spectrograph at the Keck Observatory. First reported by Roberts-Borsani et al. (2015), it was selected for spectroscopic observations because of its photometric redshift ($z_{phot}=8.57^{+0.22}_{-0.43}$), apparent brightness (H$_{160}=25.26\pm0.09$) and red Spitzer/IRAC [3.6]-[4.5] color indicative of contamination by strong oxygen emission in the [4.5] band. With a total integration of $\sim4.3$ hours, our data reveal an emission line at $\simeq11776$ {\AA} which we argue is likely Ly$\alpha$ at a redshift $z_{spec}=8.68$, in good agreement with the photometric estimate. The line was detected independently on two nights using different slit orientations and its detection significance is $\sim7.5\sigma$. An overlapping sky line contributes significantly to the uncertainty on the total line flux but not the overall significance. By direct addition and a Gaussian fit, we estimate a 95\% confidence range of 1.0 - 2.5 $\times 10^{-17}$ ergs cm$^{-2}$ sec$^{-1}$, corresponding to a rest-frame equivalent width of 17 - 42 \AA\ . EGSY8p7 is the most distant galaxy confirmed spectroscopically to date, and the third luminous source in the EGS field beyond $z_{phot}\gtrsim7.5$ with detectable Ly$\alpha$ emission viewed at a time when the intergalactic medium is expected to be fairly neutral. Although the reionization process was probably patchy, we discuss whether luminous sources with prominent IRAC color excesses may harbor harder ionizing spectra than the dominant fainter population thereby creating earlier ionized bubbles. Further spectroscopic follow-up of such bright sources promises important insight into the early formation of galaxies.
Infrared avalanche photodiode arrays represent a panacea for many branches of astronomy by enabling extremely low-noise, high-speed and even photon-counting measurements at near-infrared wavelengths. We recently demonstrated the use of an early engineering-grade infrared avalanche photodiode array that achieves a correlated double sampling read noise of 0.73 e- in the lab, and a total noise of 2.52 e- on sky, and supports simultaneous high-speed imaging and tip-tilt wavefront sensing with the Robo-AO visible-light laser adaptive optics system at the Palomar Observatory 1.5-m telescope. We report here on the improved image quality achieved simultaneously at visible and infrared wavelengths by using the array as part of an image stabilization control-loop with adaptive-optics sharpened guide stars. We also discuss a newly enabled survey of nearby late M-dwarf multiplicity as well as future uses of this technology in other adaptive optics and high-contrast imaging applications.
We characterize the astrometric distortion at the edges of thick, fully-depleted CCDs in the lab using a bench-top simulation of LSST observing. By illuminating an array of forty thousand pinholes (30mu m diameter) at the object plane of a f/1.2 optical reimager, thousands of PSFs can be imaged over a 4Kx4K pixel CCD. Each high purity silicon pixel, 10mu m square by 100mu m deep, can then be individually characterized through a series of sub-pixel dithers in the X/Y plane. The unique character [response, position, shape] of each pixel as a function of flux, wavelength, back side bias, etc. can be investigated. We measure the magnitude and onset of astrometric error at the edges of the detector as a test of the experimental setup, using a LSST prototype CCD. We show that this astrometric error at the edge is sourced from non-uniformities in the electric field lines that define pixel boundaries. This edge distortion must be corrected in order to optimize the science output of weak gravitational lensing and large scale structure measurements for the LSST.
SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter halos in order to better understand the abundance and evolution of the primary source of ionizing radiation. A series of high-resolution N-body simulations is run to quantify the abundance of dark matter halos as a function of mass $M$, accretion rate $\dot{M}$, and redshift $z$. A new fit for the halo mass function $dn/dM$ is $\approx 20\%$ more accurate at the high-mass end where bright galaxies are expected to reside. A novel approach is used to fit the halo accretion rate function $dn/d\dot{M}$ in terms of the halo mass function. Abundance matching against the observed galaxy luminosity function is used to estimate the luminosity-mass relation and the luminosity-accretion-rate relation. The inferred star formation efficiency is not monotonic with $M$ nor $\dot{M}$, but reaches a maximum value at a characteristic mass $\sim 2 \times 10^{11}\ M_\odot$ and a characteristic accretion rate $\sim 6 \times 10^2\ M_\odot/{\rm yr}$ at $z \approx 6$. We find a universal EoR luminosity-accretion-rate relation, which is used to construct a fiducial model for the galaxy luminosity function. The Schechter parameters evolve such that $\phi_\star$ decreases, $M_\star$ is more positive (fainter), and $\alpha$ is more negative (steeper) at higher redshifts. We forecast for the upcoming James Webb Space Telescope and show that with apparent magnitude limit $m_{\rm AB} \approx 31\ (32)$, it can observe $\gtrsim 11\ (24)$ unlensed galaxies per square degree per unit redshift at least down to $M_\star$ at $z \lesssim 13\ (14)$.
We present a novel approach to estimating the intensity mapping signal of any CO rotational line emitted during the Epoch of Reionization (EoR). Our approach is based on large velocity gradient (LVG) modeling, a radiative transfer modeling technique that generates the full CO spectral line energy distribution (SLED) for a specified gas kinetic temperature, volume density, velocity gradient, molecular abundance, and column density. These parameters, which drive the physics of CO transitions and ultimately dictate the shape and amplitude of the CO SLED, can be linked to the global properties of the host galaxy, mainly the star formation rate (SFR) and the SFR surface density. By further employing an empirically derived SFR-M relation for high redshift galaxies, we can express the LVG parameters, and thus the specific intensity of any CO rotational transition, as functions of the host halo mass M and redshift z. Integrating over the range of halo masses expected to host CO-luminous galaxies, i.e. M >= 10^8 M{_\odot}, we predict a mean CO(1-0) brightness temperature ranging from ~1 {\mu}K at z = 6 to ~ 0.2 {\mu}K at z = 10 in the case where the duty cycles of star formation and CO luminous activity are assumed to be 0.1 (f_{UV} = f_{duty} = 0.1). In this model, the CO emission signal remains strong for higher rotational levels, with < T_{CO} > ~ 0.3 and 0.1 {\mu}K for the CO J = 10->9 transition at z = 6 and 10 respectively. If instead we adopt duty cycles of unity, the estimated CO(1-0) brightness temperature declines to < T_{CO}>~ 0.6 {\mu}K at z = 6 and ~0.03 {\mu}K at z =10 respectively; the correspondingly reduced signal strengths of the higher J lines make detection of these transitions at high significance less likely in the f_{UV} = f_{duty} = 1 model.
We carried out the largest ($>3.5\times10^5$ Mpc$^3$, 26 deg$^2$) H$\alpha$ narrow band survey to date at $z\sim0.2$ in the SA22, W2 and XMMLSS extragalactic fields. Our survey covers a large enough volume to overcome cosmic variance and to sample bright and rare H$\alpha$ emitters up to an observed luminosity of $\sim10^{42.4}$ erg s$^{-1}$, equivalent to $\sim11 M_\odot$ yr$^{-1}$. Using our sample of $220$ sources brighter than $>10^{41.4}$ erg s$^{-1}$ ($>1 M_\odot$ yr$^{-1}$), we derive H$\alpha$ luminosity functions, which are well described by a Schechter function with $\phi^* = 10^{-2.85\pm0.03}$ Mpc$^{-3}$ and $L^*_{H\alpha} = 10^{41.71\pm0.02}$ erg s$^{-1}$ (with a fixed faint end slope $\alpha=-1.35$). We find that surveys probing smaller volumes ($\sim3\times10^4$ Mpc$^3$) are heavily affected by cosmic variance, which can lead to errors of over $100$ per cent in the characteristic density and luminosity of the H$\alpha$ luminosity function. We derive a star formation rate density of $\rho_\mathrm{SFRD} = 0.0094\pm0.0008$ $M_\odot$ yr$^{-1}$, in agreement with the redshift-dependent H$\alpha$ parametrisation from Sobral et al. (2013). The two-point correlation function is described by a single power law $\omega(\theta) = (0.159\pm0.012) \theta^{(-0.75\pm0.05)}$, corresponding to a clustering length of $r_0 = 3.3\pm0.8$ Mpc/h. We find that the most luminous H$\alpha$ emitters at $z\sim0.2$ are more strongly clustered than the relatively fainter ones. The $L^*_{H\alpha}$ H$\alpha$ emitters at $z\sim0.2$ in our sample reside in $\sim10^{12.5-13.5}$ $M_\odot$ dark matter haloes. This implies that the most star forming galaxies always reside in relatively massive haloes or group-like environments and that the typical host halo mass of star-forming galaxies is independent of redshift if scaled by $L_\mathrm{H\alpha}/L^*_{H\alpha}(z)$, as proposed by Sobral et al. (2010).
We perform a detailed modelling of the post-outburst surface emission of the low magnetic field magnetar SGR 0418+5729. The dipolar magnetic field of this source, B=6x10^12 G estimated from its spin-down rate, is in the observed range of magnetic fields for normal pulsars. The source is further characterized by a high pulse fraction and a single-peak profile. Using synthetic temperature distribution profiles, and fully accounting for the general-relativistic effects of light deflection and gravitational redshift, we generate synthetic X-ray spectra and pulse profiles that we fit to the observations. We find that asymmetric and symmetric surface temperature distributions can reproduce equally well the observed pulse profiles and spectra of SGR 0418. Nonetheless, the modelling allows us to place constraints on the system geometry (i.e. the angles $\psi$ and $\xi$ that the rotation axis makes with the line of sight and the dipolar axis, respectively), as well as on the spot size and temperature contrast on the neutron star surface. After performing an analysis iterating between the pulse profile and spectra, as done in similar previous works, we further employed, for the first time in this context, a Markov-Chain Monte-Carlo approach to extract constraints on the model parameters from the pulse profiles and spectra, simultaneously. We find that, to reproduce the observed spectrum and flux modulation: (a) the angles must be restricted to $65\deg < \psi+\xi < 125\deg$ or $235\deg < \psi+\xi <295\deg$; (b) the temperature contrast between the poles and the equator must be at least a factor of $\sim6$, and (c) the size of the hottest region ranges between 0.2-0.7 km (including uncertainties on the source distance). Last, we interpret our findings within the context of internal and external heating models.
We present a study of the X-ray flaring activity of Sgr A* during all the 150 XMM-Newton and Chandra observations pointed at the Milky Way center over the last 15 years. This includes the latest XMM-Newton and Chandra campaigns devoted to monitoring the closest approach of the very red Br-Gamma emitting object called G2. The entire dataset analysed extends from September 1999 through November 2014. We employed a Bayesian block analysis to investigate any possible variations in the characteristics (frequency, energetics, peak intensity, duration) of the flaring events that Sgr A* has exhibited since their discovery in 2001. We observe that the total bright-or-very bright flare luminosity of Sgr A* increased between 2013-2014 by a factor of 2-3 (~3.5 sigma significance). We also observe an increase (~99.9% significance) from 0.27+-0.04 to 2.5+-1.0 day^-1 of the bright-or-very bright flaring rate of Sgr A*, starting in late summer 2014, which happens to be about six months after G2's peri-center passage. This might indicate that clustering is a general property of bright flares and that it is associated with a stationary noise process producing flares not uniformly distributed in time (similar to what is observed in other quiescent black holes). If so, the variation in flaring properties would be revealed only now because of the increased monitoring frequency. Alternatively, this may be the first sign of an excess accretion activity induced by the close passage of G2. More observations are necessary to distinguish between these two hypotheses.
This work introduces a new fitting formalism for isophotes which enables more accurate modelling of galaxies with non-elliptical shapes, such as disk galaxies viewed edge-on or galaxies with X-shaped/peanut bulges. Within this scheme, the angular parameter which defines quasi-elliptical isophotes is transformed from the commonly used, but inappropriate, polar co-ordinate to the `eccentric anomaly'. This provides a superior description of deviations from ellipticity, better capturing the true isophotal shape. Furthermore, this makes it possible to accurately recover both the surface brightness profile, using the correct azimuthally-averaged isophote, and the two-dimensional model of any galaxy: the hitherto ubiquitous, but artificial, cross-like features in residual images are completely removed. The formalism has been implemented into the IRAF tasks $Ellipse$ and $Bmodel$ to create the new tasks `$Isofit$', and `$Cmodel$'. The new tools are demonstrated here with application to five galaxies, chosen to be representative case-studies for several areas where this technique makes it possible to gain new scientific insight. Specifically: properly quantifying boxy/disky isophotes via the fourth harmonic order in edge-on galaxies, quantifying X-shaped/peanut bulges, higher-order Fourier moments for modelling bars in disks, and complex isophote shapes. Higher order (n > 4) harmonics now become meaningful and may correlate with structural properties, as boxyness/diskyness is known to do. This work also illustrates how the accurate construction, and subtraction, of a model from a galaxy image facilitates the identification and recovery of over-lapping sources such as globular clusters and the optical counterparts of X-ray sources.
This paper develops a novel numerical method for obtaining structures of rapidly rotating stars based on a self-consistent field scheme. The solution is obtained iteratively. Both rapidly rotating barotropic and baroclinic equilibrium states are calculated self-consistently using this method. Two types of rotating baroclinic stars are investigated by changing the isentropic surfaces inside the star. Solution sequences of these are calculated systematically and critical rotation models beyond which no rotating equilibrium state exists are also obtained. All of these rotating baroclinic stars satisfy necessarily the Bjerknes-Rosseland rules. Self-consistent solutions of baro-clinic stars with shellular-type rotation are successfully obtained where the isentropic surfaces are oblate and the surface temperature is hotter at the poles than at the equator if it is assumed that the star is an ideal gas star. These are the first self-consistent and systematic solutions of rapidly rotating baroclinic stars with shellular-type rotations. Since they satisfy the stability criterion due to their rapid rotation, these rotating baroclinic stars would be dynamically stable. This novel numerical method and the solutions of the rapidly rotating baroclinic stars will be useful for investigating stellar evolution with rapid rotations.
We present extensive observations of the Type Ib/c SN2013ge from -13 to +457 days, including spectra and Swift UV-optical photometry beginning 2-4 days post explosion. This makes SN2013ge one of the best observed normal Type Ib/c SN at early times, when the light curve is particularly sensitive to the progenitor configuration and mixing of radioactive elements. These early observations reveal two distinct light curve components in the UV bands. The first component rises over 4-5 days and is visible for a week post-explosion. Spectra of the first component have a blue continuum and show a plethora of high velocity (~14,000 km/s) but narrow (~3500 km/s) features, indicating that the line forming region is restricted. The explosion parameters estimated for the bulk explosion are standard for Type Ib/c SN, while detailed analysis of optical and NIR spectra identify weak He features at early times, and nebular spectra show evidence for mixing and asymmetry in the bulk ejecta. In addition, SN2013ge exploded in a low metallicity environment and we have obtained some of the deepest radio and X-ray limits for a Type Ib/c SN to date that constrain the progenitor mass-loss rate. We are left with two distinct progenitor scenarios for SN2013ge depending on our interpretation of the early emission. If the first component is cooling envelope emission, then the progenitor of SN2013ge possessed a low-mass extended envelope. Alternatively, if the first component is due to outwardly mixed Ni-56 then our observations are consistent with the asymmetric ejection of a small amount of mass ahead of the bulk explosion. Current models for the collision of a SN shock with a binary companion cannot reproduce both the timescale and luminosity of the early emission in SN2013ge. Finally, we find that the spectra of the first component of SN2013ge are similar to those of the rapidly-declining SN2002bj.
We report the detection of the Galactic nuclear disc in line-of-sight kinematics of stars, measured with infrared spectroscopy from APOGEE. The nuclear disc is found to have a rotation velocity V ~ 120km/s comparable to the gas disc. The current data suggest that this disc is kinematically quite cold and has a small vertical extent of order 50pc. The stellar kinematics suggest a truncation radius of the stellar disc at a galactocentric radius R ~ 150pc, and provide tentative evidence for an overdensity at the position of the ring found in the molecular gas disc.
Correlations between the accretion luminosity and emission line luminosities (L_acc and L_line) of pre-main sequence (PMS) stars have been published for many different spectral lines, which are used to estimate accretion rates. Despite the origin of those correlations is unknown, this could be attributed to direct or indirect physical relations between the emission line formation and the accretion mechanism. This work shows that all (near-UV/optical/near-IR) L_acc-L_line correlations are the result of the fact that the accretion luminosity and the stellar luminosity (L_star) are correlated, and are not necessarily related with the physical origin of the line. Synthetic and observational data are used to illustrate how the L_acc-L_line correlations depend on the L_acc-L_star relationship. We conclude that because PMS stars show the L_acc-L_star correlation immediately implies that L_acc also correlates with the luminosity of all emission lines, for which the L_acc-L_line correlations alone do not prove any physical connection with accretion but can only be used with practical purposes to roughly estimate accretion rates. When looking for correlations with possible physical meaning, we suggest that L_acc/L_star and L_line/L_star should be used instead of L_acc and L_line. Finally, the finding that L_acc has a steeper dependence on L_star for T-Tauri stars than for intermediate-mass Herbig Ae/Be stars is also discussed. That is explained from the magnetospheric accretion scenario and the different photospheric properties in the near-UV.
This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlation functions of CMB temperature and polarization. They use the hybrid approach employed previously: pixel-based at low multipoles, $\ell$, and a Gaussian approximation to the distribution of cross-power spectra at higher $\ell$. The main improvements are the use of more and better processed data and of Planck polarization data, and more detailed foreground and instrumental models. More than doubling the data allows further checks and enhanced immunity to systematics. Progress in foreground modelling enables a larger sky fraction, contributing to enhanced precision. Improvements in processing and instrumental models further reduce uncertainties. Extensive tests establish robustness and accuracy, from temperature, from polarization, and from their combination, and show that the {\Lambda}CDM model continues to offer a very good fit. We further validate the likelihood against specific extensions to this baseline, such as the effective number of neutrino species. For this first detailed analysis of Planck polarization, we concentrate at high $\ell$ on E modes. At low $\ell$ we use temperature at all Planck frequencies along with a subset of polarization. These data take advantage of Planck's wide frequency range to improve the separation of CMB and foregrounds. Within the baseline cosmology this requires a reionization optical depth $\tau=0.078\pm0.019$, significantly lower than without high-frequency data for explicit dust monitoring. At high $\ell$ we detect residual errors in E, typically at the {\mu}K$^2$ level; we thus recommend temperature alone as the high-$\ell$ baseline. Nevertheless, Planck high-$\ell$ polarization spectra are already good enough to allow a separate high-accuracy determination of the {\Lambda}CDM parameters, consistent with those established from temperature alone.
The inverse Compton catastrophe is defined as a dramatic rise in the luminosity of inverse Compton scattered photons. It is described by a non-linear loop of radiative processes that sets in for high values of the electron compactness and is responsible for the efficient transfer of energy from electrons to photons, predominantly through inverse Compton scatterings. We search for the conditions that drive a magnetized non-thermal source to the inverse Compton catastrophe regime and study its multi-wavelength (MW) photon spectrum. We develop a generic analytical framework and use numerical calculations as a backup to the analytical predictions. We find that the escaping radiation from a source in the Compton catastrophe regime bears some unique features. The MW photon spectrum is a broken power law with a break at $\sim m_e c^2$ due to the onset of the Klein-Nishina suppression. The spectral index below the break energy depends on the electron and magnetic compactnesses logarithmically, while it is independent of the electron power-law index ($s$). The maximum radiating power emerges typically in the $\gamma$-ray regime, at energies $\sim m_e c^2$ ($\sim \gamma_{\max} m_e c^2$ ) for $s>2$ ($s\lesssim 2$), where $\gamma_{\max}$ is the maximum Lorentz factor of the injected electron distribution. We apply the principles of the inverse Compton catastrophe to blazars and $\gamma$-ray bursts using the analytical framework we developed, and show how these can be used to impose robust constraints on the source parameters.
The GAPS experiment is foreseen to carry out a dark matter search by
measuring low-energy cosmic-ray antideuterons and antiprotons with a novel
detection approach. It will provide a new avenue to access a wide range of
different dark matter models and masses from about 10GeV to 1TeV. The
theoretically predicted antideuteron flux resulting from secondary interactions
of primary cosmic rays is very low. Well-motivated theories beyond the Standard
Model contain viable dark matter candidates, which could lead to a significant
enhancement of the antideuteron flux due to annihilation or decay of dark
matter particles. This flux contribution is believed to be especially large at
low energies, which leads to a high discovery potential for GAPS. The GAPS
low-energy antiproton search will provide some of the most stringent
constraints on ~30GeV dark matter, will provide the best limits on primordial
black hole evaporation on galactic length scales, and explore new discovery
space in cosmic-ray physics.
GAPS is designed to achieve its goals via long duration balloon flights at
high altitude in Antarctica. The detector itself will consist of 10 planes of
Si(Li) solid state detectors and a surrounding time-of-flight system.
Antideuterons and antiprotons will be slowed down in the Si(Li) material,
replace a shell electron and form an excited exotic atom. The atom will be
deexcited by characteristic X-ray transitions and will end its life by the
formation of an annihilation pion/proton star. This unique event structure will
deliver a nearly background free detection possibility.
This review paper summarizes the state-of-the-art in laboratory based interstellar ice chemistry. The focus is on atom addition reactions, illustrating how water, carbon dioxide and methanol can form in the solid state at astronomically relevant temperatures, and also the formation of more complex species such as hydroxylamine, an important prebiotic molecule, and glycolaldehyde, the smallest sugar, is discussed. These reactions are particularly relevant during the dark ages of star and planet formation, i.e., when the role of UV light is restricted. A quantitative characterization of such processes is only possible through dedicated laboratory studies, i.e., under full control of a large set of parameters such as temperature, atom-flux, and ice morphology. The resulting numbers, physical and chemical constants, e.g., barrier heights, reaction rates and branching ratios, provide information on the molecular processes at work and are needed as input for astrochemical models, in order to bridge the timescales typical for a laboratory setting to those needed to understand the evolutionary stages of the interstellar medium. Details of the experiments as well as the astrochemical impact of the results are discussed.
Gas motions in galaxy clusters play important roles in determining the properties of the intracluster medium (ICM) and constraining cosmological parameters using X-ray and Sunyaev-Zel'dovich effect observations of galaxy clusters. The upcoming ASTRO-H mission, equipped with high-resolution X-ray spectrometer, will make the first direct measurements of gas motions in galaxy clusters through measurements of Doppler shifting and broadening of emission lines. However, the physical interpretation of the data will be challenging due to the complex thermal and velocity structures of the ICM. In this work, we investigate how well we can measure bulk and turbulent gas motions in the ICM with ASTRO-H, by analyzing mock ASTRO-H simulations of galaxy clusters extracted from cosmological hydrodynamic simulations. We assess how photon counts, spectral fitting methods, multiphase ICM structure, deprojections, and region selection affect the measurements of gas motions. We show that while ASTRO-H is capable of recovering the underlying spherically averaged velocity profiles to within 20% with reasonable amount of photon counts (>~200) in the 6.7 keV Fe XXV line complex, there are considerable azimuthal variations in the ICM velocities, even in dynamically relaxed systems, which must be taken into account when interpreting data and developing observing strategies. Finally, we show that ASTRO-H should enable direct measurements of the hydrostatic mass bias with an accuracy of <~5%, by accounting for both rotational and random velocities from Doppler shifts and broadening of emission lines. Our results are broadly applicable for future X-ray missions, such as Athena+ and SMART-X.
Using the relation between distance modulus (m-M) and redshift (z), deduced from Friedman-Robertson-Walker (FRW) metric and assuming different values of deceleration parameter (q0). We constrained the Hubble parameter (h). The estimates of the Hubble parameters we obtained using the median values of the data obtained from NASA Extragalactic Database (NED), are: h=0.7+/-0.3 for q0=0, h=0.6+/-0.3, for q0=1 and h=0.8+/-0.3, for q0=-1. The corresponding age ({\tau}) and size (R) of the observable universe were also estimated as: {\tau}=15+/-1 Gyrs, R=(5+/-2)x10^3 Mpc, {\tau}=18+/-1 Gyrs, R=(6+/-2)x10^3 Mpc and {\tau}=13+/-1 Gyrs, R=(4+/-2)x10^3 Mpc for q0=0, q0=1 and q0=-1 respectively.
The Next Generation Virgo Cluster Survey has recently determined the luminosity function of galaxies in the core of the Virgo cluster down to unprecedented magnitude and surface brightness limits. Comparing simulations of cluster formation to the derived central stellar mass function, we attempt to estimate the stellar-to-halo-mass ratio (SHMR) for dwarf galaxies, as it would have been before they fell into the cluster. This approach ignores several details and complications, e.g., the contribution of ongoing star formation to the present-day stellar mass of cluster members, and the effects of adiabatic contraction and/or violent feedback on the subhalo and cluster potentials. The final results are startlingly simple, however; we find that the trends in the SHMR determined previously for bright galaxies appear to extend down in a scale-invariant way to the faintest objects detected in the survey. These results extend measurements of the formation efficiency of field galaxies by two decades in halo mass, or five decades in stellar mass, down to some of the least massive dwarf galaxies known, with stellar masses of $\sim 10^5 M_\odot$.
We use the Jansson-Farrar JF12 magnetic field configuration in the context of point source searches by correlating the Telescope Array ultra-high energy cosmic ray data and the IceCube-40 neutrino candidates. As expected, we have found no correlations, thus, we devote this paper to the study of the effect of different magnetic field hypotheses on the minimum neutrino source flux strength required for a $5\sigma$ discovery and the derived $90\%$ CL upper limits. In this study we present a comparison between the JF12 field, that includes a combination of regular and random field components, and the standard turbulent magnetic field used in previous correlation analyses. For a wider perspective, we also incorporate in our comparison the cases of no magnetic field and the JF12 regular component alone and consider different power law indices $\alpha=2,\alpha=2.3$ for the neutrino point source flux. Collaterally, a novel parameterisation of the JF12 random component is introduced. We have observed that the discovery potential for a point source search is sensitive to changes in the magnetic field assumptions, being the difference between the models mentioned before between 15\% and 26\%.
We analyze the Suzaku XIS data of the central region of supernova remnant G332.5-5.6. The X-ray data are well described by a single non-equilibrium ionization thermal model, {\tt vnei}, with an absorbing hydrogen column density of 1.4$^{+0.4}_{-0.1}$ $\times$ 10$^{21}$ cm$^{-2}$. The plasma is characterized by an electron temperature of 0.49$^{+0.08}_{-0.06}$ keV with subsolar abundances for O (0.58$^{+0.06}_{-0.05}$ solar value) and Fe (0.72$^{+0.06}_{-0.05}$ solar value) and slightly overabundance for Mg (1.23$^{+0.14}_{-0.14}$ solar value). It seems that the central X-ray emission originates from projection effect or evaporation of residual clouds inside G332.5-5.6. We estimate a distance of 3.0 $\pm$ 0.8 kpc for G332.5-5.6 based on the extinction-distance relation. G332.5-5.6 has an age of 7 - 9 kyr.
Recently, two independent groups found very different results when measuring the central velocity dispersion of the galactic globular cluster NGC 6388 with different methods. While L\"utzgendorf et al. (2011) found a rising profile and a high central velocity dispersion (23.3 km/s), measurements obtained by Lanzoni et al. (2013) showed a value 40% lower. The value of the central velocity dispersion has a serious impact on the mass and possible presence of an intermediate-mass black hole at the center of NGC 6388. We use a photometric catalog of NGC 6388 to create a simulated SINFONI and ARGUS dataset. The construction of the IFU data cube is done with different observing conditions reproducing the conditions reported for the original observations as closely as possible. In addition, we produce an N-body realization of a 10^6 M_SUN stellar cluster with the same photometric properties as NGC 6388 to account for unresolved stars. We find that the individual radial velocities, i.e. the measurements from the simulated SINFONI data, are systematically biased towards lower velocity dispersions. The reason is that due to the wings in the point spread function the velocities get biased towards the mean cluster velocity. This study shows that even with AO supported observations, individual radial velocities in crowded fields are likely to be biased. The ARGUS observations do not show this kind of bias but were found to have larger uncertainties than previously obtained. We find a bias towards higher velocity dispersions in the ARGUS pointing when fixing the extreme velocities of the three brightest stars but find those variations are within the determined uncertainties. We rerun Jeans models and fit the kinematic profile with the new uncertainties. This yields a BH mass of M_BH = (2.8 +- 0.4) x 10^4 M_SUN and M/L ratio M/L = (1.6 +- 0.1) M_SUN/L_SUN, consistent with our previous results.
We find that variations in the UV emissions of cool M-dwarf stars have a potentially large impact upon atmospheric biosignatures in simulations of Earth-like exoplanets i.e. planets with Earths development, and biomass and a molecular nitrogen-oxygen dominated atmosphere. Starting with an assumed black-body stellar emission for an M7 class dwarf star, the stellar UV irradiation was increased stepwise and the resulting climate-photochemical response of the planetary atmosphere was calculated. Results suggest a Goldilocks effect with respect to the spectral detection of ozone. At weak UV levels, the ozone column was weak (due to weaker production from the Chapman mechanism) hence its spectral detection was challenging. At strong UV levels, ozone formation is stronger but its associated stratospheric heating leads to a weakening in temperature gradients between the stratosphere and troposphere, which results in weakened spectral bands. Also, increased UV levels can lead to enhanced abundances of hydrogen oxides which oppose the ozone formation effect. At intermediate UV (i.e. with x10 the stellar UV radiative flux of black body Planck curves corresponding to spectral class M7) the conditions are just right for spectral detection. Results suggest that the planetary O3 profile is sensitive to the UV output of the star from about(200-350) nm. We also investigated the effect of increasing the top-of-atmosphere incoming Lyman-alpha radiation but this had only a minimal effect on the biosignatures since it was efficiently absorbed in the uppermost planetary atmospheric layer, mainly by abundant methane. Earlier studies have suggested that the planetary methane is an important stratospheric heater which critically affects the vertical temperature gradient, hence the strength of spectral emission bands.
I compute the average polarisation asymmetry from the Klein-Nishina differential cross section on free electrons at rest. As expected from the expression for the asymmetry, the average asymmetry is found to decrease like the inverse of the incident photon energy asymptotically at high energy. I then compute a simple estimator of the polarisation fraction that makes optimal use of all the kinematic information present in an event final state, by the use of "moments" method, and I compare its statistical power to that of a simple fit of the azimuthal distribution. In contrast to polarimetry with pair creation, for which I obtained an improvement by a factor of larger than two in a previous work, here for Compton scattering the improvement is only of 10-20 %.
It has long been proposed that low frequency QPOs in stellar mass black holes or their equivalents in super massive black holes are results of resonances between infall and cooling time scales. We explicitly compute these two time scales in a generic situation to show that resonances are easily achieved. During an outburst of a transient black hole candidate (BHC), the accretion rate of the Keplerian disk as well as the geometry of the Comptonizing cloud change very rapidly. During some period, resonance condition between the cooling time scale (predominantly by Comptonization) and the infall time scale of the Comptonizing cloud is roughly satisfied. This leads to low frequency quasi-periodic oscillations (LFQPOs) of the Compton cloud and the consequent oscillation of hard X-rays. In this paper, we explicitly follow the BHC H 1743-322 during its 2010 outburst. We compute Compton cooling time and infall time on several days and show that QPOs take place when these two roughly agree within ~50%, i.e., the resonance condition is generally satisfied. We also confirm that for the sharper LFQPOs (i.e., higher Q-factors) the ratio of two time scales is very close to 1.
Discovered more than 40 years ago, impulsive solar energetic particle (SEP) events are still poorly understood. The enormous abundance enhancement of the rare 3He isotope is the most striking feature of these events, though large enhancements in heavy and ultra-heavy nuclei are also observed. Recurrent 3He-rich SEPs in impulsive events have only been observed for limited time periods, up to a few days which is typically the time that a single stationary spacecraft is magnetically connected to the source active regions on the Sun. With the launch of the two STEREO spacecraft we now have the possibility of longer connection time to solar active regions. We examined the evolution of source regions showing repeated 3He-rich SEP emissions for relatively long time periods. We found that recurrent 3He-rich SEPs in these long-lived sources occur after the emergence of magnetic flux.
The last solar minimum was, by recent standards, unusually deep and long. We are now close to the maximum of the subsequent solar cycle, which is relatively weak. In this article we make comparisons between different global (unresolved) measures of the Sun's magnetic activity, to investigate how they are responding to this weak-activity epoch. We focus on helioseismic data, which are sensitive to conditions, including the characteristics of the magnetic field, in the solar interior. Also considered are measures of the magnetic field in the photosphere (sunspot number and sunspot area), the chromosphere and corona (10.7cm radio flux and 530.3nm green coronal index), and two measures of the Sun's magnetic activity closer to Earth (the interplanetary magnetic field and the galactic cosmic-ray intensity). Scaled versions of the activity proxies diverge from the helioseismic data around 2000, indicating a change in relationship between the proxies. The degree of divergence varies from proxy to proxy with sunspot area and 10.7cm flux showing only small deviations, while sunspot number, coronal index, and the two interplanetary proxies show much larger departures. In Cycle 24 the deviations in the solar proxies and the helioseismic data decrease, raising the possibility that the deviations observed in Cycle 23 are just symptomatic of a 22-year Hale cycle. However, the deviations in the helioseismic data and the interplanetary proxies increase in Cycle 24. Interestingly the divergence in the solar proxies and the helioseismic data are not reflected in the shorter-term variations (often referred to as quasi-biennial oscillations) observed on top of the dominant 11-year solar cycle. However, despite being highly correlated in Cycle 22, the short-term variations in the interplanetary proxies show very little correlation with the helioseismic data during Cycles 23 and 24.
Aims. We aim to develop a chemical model that contains a consistent
description of spin-state chemistry in reactions involving chemical species
with multiple deuterons. We apply the model to the specific case of deuterated
ammonia, to derive values for the various spin-state ratios.
Methods. We apply symmetry rules in the complete scrambling assumption to
calculate branching ratio tables for reactions between chemical species that
include multiple protons and/or deuterons. Reaction sets for both gas-phase and
grain-surface chemistry are generated using an automated routine that forms all
possible spin-state variants of any given reaction with up to six H/D atoms.
Single-point and modified Bonnor-Ebert models are used to study the density and
temperature dependence of ammonia and its isotopologs, and the associated
spin-state ratios.
Results. We find that the spin-state ratios of the ammonia isotopologs are,
at late times, very different from their statistical values. The ratios are
rather insensitive to variations in the density, but present strong temperature
dependence. We derive high peak values ($\sim$ 0.1) for the deuterium fraction
in ammonia, in agreement with previous (gas-phase) models. The deuterium
fractionation is strongest at high density, corresponding to a high degree of
depletion, and also presents temperature dependence. We find that in the
temperature range 5 to 20 K, the deuterium fractionation peaks at $\sim$ 15 K
while most of the ortho/para (and meta/para for $\rm ND_3$) ratios present a
minimum at 10 K (ortho/para $\rm NH_2D$ has instead a maximum at this
temperature).
Conclusions. Owing to the density and temperature dependence found in the
abundances and spin-state ratios of ammonia and its isotopologs, it is evident
that observations of ammonia and its deuterated forms can provide important
constraints on the physical structure of molecular clouds.
Whether the new line at ~3.5 keV, recently detected in different samples of galaxy clusters, Andromeda galaxy and central part of our Galaxy, is due to Potassium emission lines, is now unclear. By using the latest astrophysical atomic emission line database AtomDB v. 3.0.2, we show that the most prospective method to directly check its Potassium origin will be the study of K XIX emission line complex at ~3.7 keV with future X-ray imaging spectrometers such as Soft X-ray spectometer on-board Astro-H mission or microcalorimeter on-board Micro-X sounding rocket experiment. To further reduce the remaining (factor ~3-5) uncertainty of the 3.7/3.5 keV ratio one should perform more precise modeling including removal of significant spatial inhomogeneities, detailed treatment of background components, and further extension of the modeled energy range.
The currently operating space missions, as well as those that will be launched in the near future, (will) deliver high-quality data for millions of stellar objects. Since the majority of stellar astrophysical applications still (at least partly) rely on spectroscopic data, an efficient tool for the analysis of medium- to high-resolution spectroscopy is needed. We aim at developing an efficient software package for the analysis of medium- to high-resolution spectroscopy of single stars and those in binary systems. The major requirements are that the code has a high performance, represents the state-of-the-art analysis tool, and provides accurate determinations of atmospheric parameters and chemical compositions for different types of stars. We use the method of atmosphere models and spectrum synthesis, which is one of the most commonly used approaches for the analysis of stellar spectra. Our Grid Search in Stellar Parameters (GSSP) code makes use of the OpenMPI implementation, which makes it possible to run in parallel mode. The method is first tested on the simulated data and is then applied to the spectra of real stellar objects. The majority of test runs on the simulated data were successful in the sense that we could recover the initially assumed sets of atmospheric parameters. We experimentally find the limits in signal-to-noise ratios of the input spectra, below which the final set of parameters gets significantly affected by the noise. Application of the GSSP package to the spectra of three Kepler stars, KIC11285625, KIC6352430, and KIC4931738, was also largely successful. We found an overall agreement of the final sets of the fundamental parameters with the original studies. For KIC6352430, we found that dependence of the light dilution factor on wavelength cannot be ignored, as it has significant impact on the determination of the atmospheric parameters of this binary system.
Extreme deconvolution (XD) of broad-band photometric data can both separate stars from quasars and generate probability density functions for quasar redshifts, while incorporating flux uncertainties and missing data. Mid-infrared photometric colors are now widely used to identify hot dust intrinsic to quasars, and the release of all-sky WISE data has led to a dramatic increase in the number of IR-selected quasars. Using forced-photometry on public WISE data at the locations of SDSS point sources, we incorporate this all-sky data into the training of the XDQSOz models originally developed to select quasars from optical photometry. The combination of WISE and SDSS information is far more powerful than SDSS alone, particularly at $z>2$. The use of SDSS$+$WISE photometry is comparable to the use of SDSS$+$ultraviolet$+$near-IR data. We release a new public catalogue of 5,537,436 (total; 3,874,639 weighted by probability) potential quasars with probability $P_{\textrm{QSO}} > 0.2$. The catalogue includes redshift probabilities for all objects. We also release an updated version of the publicly available set of codes to calculate quasar and redshift probabilities for various combinations of data. Finally, we demonstrate that this method of selecting quasars using WISE data is both more complete and efficient than simple WISE color-cuts, especially at high redshift. Our fits verify that above $z \sim 3$ WISE colors become bluer than the standard cuts applied to select quasars. Currently, the analysis is limited to quasars with optical counterparts, and thus cannot be used to find highly obscured quasars that WISE color-cuts identify in significant numbers.
The bulk chemical compositions of planets are uncertain, even for major elements such as Mg and Si. This is due to the fact that the samples available for study all originate from relatively shallow depths. Comparison of the stable isotope compositions of planets and meteorites can help overcome this limitation. Specifically, the non-chondritic Si isotope composition of the Earth's mantle was interpreted to reflect the presence of Si in the core, which can also explain its low density relative to pure Fe-Ni alloy. However, we have found that angrite meteorites display a heavy Si isotope composition similar to the lunar and terrestrial mantles. Because core formation in the angrite parent-body (APB) occurred under oxidizing conditions at relatively low pressure and temperature, significant incorporation of Si in the core is ruled out as an explanation for this heavy Si isotope signature. Instead, we show that equilibrium isotopic fractionation between gaseous SiO and solid forsterite at 1370 K in the solar nebula could have produced the observed Si isotope variations. Nebular fractionation of forsterite should be accompanied by correlated variations between the Si isotopic composition and Mg/Si ratio following a slope of 1, which is observed in meteorites. Consideration of this nebular process leads to a revised Si concentration in the Earth's core of 3.6 (+6.0/-3.6) wt% and provides estimates of Mg/Si ratios of bulk planetary bodies.
Magnetars are the strongest magnets in the present universe and the combination of extreme magnetic field, gravity and density makes them unique laboratories to probe current physical theories (from quantum electrodynamics to general relativity) in the strong field limit. Magnetars are observed as peculiar, burst--active X-ray pulsars, the Anomalous X-ray Pulsars (AXPs) and the Soft Gamma Repeaters (SGRs); the latter emitted also three "giant flares," extremely powerful events during which luminosities can reach up to 10^47 erg/s for about one second. The last five years have witnessed an explosion in magnetar research which has led, among other things, to the discovery of transient, or "outbursting," and "low-field" magnetars. Substantial progress has been made also on the theoretical side. Quite detailed models for explaining the magnetars' persistent X-ray emission, the properties of the bursts, the flux evolution in transient sources have been developed and confronted with observations. New insight on neutron star asteroseismology has been gained through improved models of magnetar oscillations. The long-debated issue of magnetic field decay in neutron stars has been addressed, and its importance recognized in relation to the evolution of magnetars and to the links among magnetars and other families of isolated neutron stars. The aim of this paper is to present a comprehensive overview in which the observational results are discussed in the light of the most up-to-date theoretical models and their implications. This addresses not only the particular case of magnetar sources, but the more fundamental issue of how physics in strong magnetic fields can be constrained by the observations of these unique sources.
Spectral analysis is a powerful tool to investigate stellar properties and it
has been widely used for decades now. However, the methods considered to
perform this kind of analysis are mostly based on iteration among a few
diagnostic lines to determine the stellar parameters. While these methods are
often simple and fast, they can lead to errors and large uncertainties due to
the required assumptions.
Here we present a method based on Bayesian statistics to find simultaneously
the best combination of effective temperature, surface gravity, projected
rotational velocity, and microturbulence velocity, using all the available
spectral lines. Different tests are discussed to demonstrate the strength of
our method, which we apply to 54 mid-resolution spectra of field and cluster B
stars obtained at the Observatoire du Mont-M\'egantic. We compare our results
with those found in the literature. Differences are seen which are well
explained by the different methods used. We conclude that the B-star
microturbulence velocities are often underestimated. We also confirm the trend
that B stars in clusters are on average faster rotators than field B stars.
We report on the analysis of the data collected by Swift, INTEGRAL and RXTE of the Black Hole Candidate (BHC) 4U 1630-47 during 3 consecutive outbursts occurred in 2006, 2008 and 2010, respectively. We show that, although a similar spectral and temporal behaviour in the energy range between 2-10 keV, these 3 outbursts present pronounced differences above 20 keV. In fact, the 2010 outburst extends at high energies without any detectable cut-off until 150-200 keV, while the other two previous outbursts, occurred on 2006 and 2008, are not detected at all above 20 keV. Moreover, the 2008 outburst does not show any detectable hard state in its final phases and even during the 2010 outburst, the final hard state shows some peculiarities rarely observed in other BHC. We also investigate on the peculiar huge variation of 4U 1630-47 hydrogen column density (N$_{H}$) reported in the literature using the Swift/XRT data. In fact this instrument is one of the most suitable for this purpose thanks to its lower energy coverage.
Gaia's astrometric solution aims to determine at least five parameters for each star, together with appropriate estimates of their uncertainties and correlations. This requires at least five distinct observations per star. In the early data reductions the number of observations may be insufficient for a five-parameter solution, and even after the full mission many stars will remain under-observed, including faint stars at the detection limit and transient objects. In such cases it is reasonable to determine only the two position parameters. Their formal uncertainties would however grossly underestimate the actual errors, due to the neglected parallax and proper motion. We aim to develop a recipe to calculate sensible formal uncertainties that can be used in all cases of under-observed stars. Prior information about the typical ranges of stellar parallaxes and proper motions is incorporated in the astrometric solution by means of Bayes' rule. Numerical simulations based on the Gaia Universe Model Snapshot (GUMS) are used to investigate how the prior influences the actual errors and formal uncertainties when different amounts of Gaia observations are available. We develop a criterion for the optimum choice of priors, apply it to a wide range of cases, and derive a global approximation of the optimum prior as a function of magnitude and galactic coordinates. The feasibility of the Bayesian approach is demonstrated through global astrometric solutions of simulated Gaia observations. With an appropriate prior it is possible to derive sensible positions with realistic error estimates for any number of available observations. Even though this recipe works also for well-observed stars it should not be used where a good five-parameter astrometric solution can be obtained without a prior. Parallaxes and proper motions from a solution using priors are always biased and should not be used.
We have performed a WISE (Wide-Field Infrared Survey Explorer) based study to identify and characterize young stellar objects (YSOs) in 12x12 degree Perseus OB2 association. Spectral energy distribution (SED) slope in range of 3.4-12 micron and a 5sigma selection criteria were used to select our initial sample. Further manual inspection reduced our final catalog to 156 known and 119 YSO candidate. The spatial distribution of newly found YSOs all over the field shows an older generation of star formation which most of its massive members have evolved into main sequence stars. In contrast, the majority of younger members lie within the Perseus molecular cloud and currently active star forming clusters such as NGC1333 and IC348. We also identified additional 66 point sources which passed YSO selection criteria but are likely AGB stars. However their spatial distribution suggests that they may contain a fraction of the YSOs. Comparing our results with the commonly used color-color selections, we found that while color selection method fails in picking up bright but evolved weak disks, our SED fitting method can identify such sources, including transitional disks. In addition we have less contamination with background sources such as galaxies, but in a price of loosing fainter (Jmag > 12) YSOs. Finally we employed a Bayesian Monte Carlo SED fitting method to determine the characteristics of each YSO candidate. Distribution of SED slopes and model driven age and mass confirms separated YSO populations with suggested three age groups of younger than 1 Myr old, 1-5 Myr old, and older than 5 Myrs which agrees with the age of Per OB2 association and currently star forming sites within the cloud.
We present a new set of analytic models for the expansion of HII regions powered by UV photoionisation from massive stars and compare them to a new suite of radiative magnetohydrodynamic simulations of self-gravitating molecular clouds. To perform these simulations we use RAMSES-RT, an Eulerian adaptive mesh magnetohydrodynamics code with radiative transfer of UV photons. We find two models that fit the simulation results well, and give a physically-motivated criterion for determining which of the models should be used. In one model, the ionisation front is only resisted by the ram pressure from the external medium, which we model as a spherically symmetric power law model. In the other, the front stalls at an equilibrium radius at which the ram pressure from accretion or turbulence in the cloud is balanced by the thermal pressure of the photoionised gas. If this stalling radius is larger than the radius of the cloud, the ionisation front can escape the cloud and expand freely in most directions. Otherwise, the front stalls at the predicted radius. We also measure the response of Jeans unstable gas to the HII regions to predict the impact of UV radiation on star formation in the cloud. We find that the mass in unstable gas can be explained by a model in which the clouds are evaporated by UV photons, suggesting that the net feedback on star formation should be negative.
We introduce a new mechanism for generating magnetic fields in the recombination era. This Harrison-like mechanism utilizes vorticity in baryons that is sourced through the Bose-Einstein condensate of axions via gravitational interactions. The magnetic fields generated are on the galactic scales $\sim 10\,{\rm kpc}$ and have a magnitude of the order of $B\sim10^{-23}\,{\rm G}$ today. The field has a greater magnitude than those generated from other mechanisms relying on second order perturbation theory, and is sufficient to provide a seed for battery mechanisms.
The temperature of a circumbinary disc edge should undulate due to variations in illumination as a function of binary orbital phase. We explore circumbinary disc temperature variations as a source of broad-band infrared light curve variability. Approximating the wall of a circumbinary disc edge as a wide optically thick cylinder with surface temperature dependent on its illumination, we find that a binary comprised of 1 M$_\odot$ and 0.5 M$_\odot$ pre-main sequence stars in a $\sim$15.5 day period, would exhibit the largest amplitude variations of $\sim$9% at 3.77 and 4.68 {\mu}m as seen by a distant observer. The amplitude of variations and shape of the light curve is sensitive to the luminosity and mass ratios of the stars in the binary, the radius of the circumbinary disc clearing, the binary separation, and the orbital inclination. The light curve variations are smooth and very red with a non-sinusoidal shape for most of the parameter space explored. Possible morphologies include a single peak with a flat region, two peaks of different heights or a single dip.
Recently, a new mechanism to generate a naturally small electroweak scale has been proposed. It exploits the coupling of the Higgs to an axion-like field and a long era in the early universe where the axion unchains a dynamical screening of the Higgs mass. We present a new realization of this idea with the new feature that it leaves no signs of new physics up to a rather large scale, 10^9 GeV, except for two very light and weakly coupled axion-like states. One of the scalars can be a viable Dark Matter candidate. Such a cosmological Higgs-axion interplay could be tested with a number of experimental strategies.
We demonstrate that current and planned underground neutrino experiments could offer a powerful probe of few-MeV dark matter when combined with a nearby high-intensity low-to-medium energy electron accelerator. This experimental setup, an underground beam-dump experiment, is capable of decisively testing the thermal freeze-out mechanism for several natural dark matter scenarios in this mass range. We present the sensitivity reach in terms of the mass-coupling parameter space of existing and planned detectors, such as Super-K, SNO+, and JUNO, in conjunction with a hypothetical 100 MeV energy accelerator. This setup can also greatly extend the sensitivity of direct searches for new light weakly-coupled force-carriers independently of their connection to dark matter.
The precise measurement of cosmic-ray antiparticles serves as important means for identifying the nature of dark matter. Recent years showed that identifying the nature of dark matter with cosmic-ray positrons and higher energy antiprotons is difficult, and has lead to a significantly increased interest in cosmic-ray antideuteron searches. Antideuterons may also be generated in dark matter annihilations or decays, offering a potential breakthrough in unexplored phase space for dark matter. Low-energy antideuterons are an important approach because the flux from dark matter interactions exceeds the background flux by more than two orders of magnitude in the low-energy range for a wide variety of models. This review is based on the "dbar14 - dedicated cosmic-ray antideuteron workshop", which brought together theorists and experimentalists in the field to discuss the current status, perspectives, and challenges for cosmic-ray antideuteron searches and discusses the motivation for antideuteron searches, the theoretical and experimental uncertainties of antideuteron production and propagation in our Galaxy, as well as give an experimental cosmic-ray antideuteron search status update. This report is a condensed summary of the article "Review of the theoretical and experimental status of dark matter identification with cosmic-ray antideuteron" (arXiv:1505.07785).
We present a definition of unsigned magnification in gravitational lensing valid on arbitrary convex normal neighborhoods of time oriented Lorentzian manifolds. This definition is a function defined at any two points along a null geodesic that lie in a convex normal neighborhood, and foregoes the usual notions of lens and source planes in gravitational lensing. Rather, it makes essential use of the van Vleck determinant, which we present via the exponential map, and Etherington's definition of luminosity distance for arbitrary spacetimes. We then specialize our definition to spacetimes, like Schwarzschild's, in which the lens is compact and isolated, and show that our magnification function is monotonically increasing along any geodesic contained within a convex normal neighborhood.
We review some recent proposals for relativistic models of dark matter in the context of bimetric gravity. The aim is to solve the problems of cold dark matter (CDM) at galactic scales, and to reproduce the phenomenology of the modified Newtonian dynamics (MOND), while still being in agreement with the standard cosmological model $\Lambda$-CDM at large scales. In this context a promising alternative is dipolar dark matter (DDM) in which two different species of dark matter particles are separately coupled to the two metrics of bigravity and are linked together by an internal vector field. The phenomenology of MOND then results from a mechanism of gravitational polarization. Probably the best formulation of the model is within the framework of recently developed massive bigravity theories. Then the gravitational sector of the model is safe by construction, but a ghostly degree of freedom in the decoupling limit is still present in the dark matter sector. Future work should analyse the cosmological solutions of the model and check the post-Newtonian parameters in the solar system.
In this study, we propose a new approach to constrain the coupling of axion-like particles (ALPs) to photons. One intriguing property of these ALPs is their mixing with photons within magnetic fields. This mixing allows photons propagating in magnetic fields to convert into ALPs and \textit{vice versa}. Plasma effects can lead to resonant conversion, further enhancing the conversion probability. For suitable ALP masses, this resonant conversion can occur for cosmic microwave background photons transversing galaxy clusters which would distort the CMB spectrum in the direction of galaxy clusters. We compare the predicted distortion with recent measurements of the thermal Sunyaev-Zeldovich Compton parameter to obtain upper limits on the coupling between photons and ALPs. The constraints apply to the mass range of approximately $2\cdot 10^{-14}$ eV $ \lesssim m_\phi \lesssim 3\cdot 10^{-12}$ eV in which resonant photon-ALP conversions can occur. Using simple galaxy cluster models, we obtain new limits for this mass range, which are up to two orders of magnitude stronger than existing ones. These limits also reduce the available parameter space for explaining the soft X-ray excess of the Coma Cluster by photon-ALP conversion.
The reliability of the first detection is one of the most interesting challenges for the gravitational wave community. To increase the detection confidence, the LIGO and Virgo collaboration have already started coincident observations between gravitational waves detectors and other astronomical instruments, like electromagnetic or neutrino detectors. This can be done in two directions: searching for gravitational waves triggered by the electromagnetic informations, or pointing the electromagnetic telescopes to the sky position given in real time by the gravitational wave analysis. The success of the latter case depends strongly on the analysis speed of gravitational wave pipelines to analyze data and extract any gravitational wave candidate with as much information as possible. In this paper we discuss the case of the coherent Waveburst pipeline, the main pipeline used in the past scientific LIGO-Virgo analyses for the search of gravitational wave transients, reporting the capability of making an all-sky and all-time analysis and the analysis speed performance.
Links to: arXiv, form interface, find, astro-ph, recent, 1507, contact, help (Access key information)
We report the discovery and early evolution of ASASSN-15lh, the most luminous supernova ever found. At redshift z=0.2326, ASASSN-15lh reached an absolute magnitude of M_{u,AB} ~ -23.5 and bolometric luminosity L_bol ~ 2.2x10^45 ergs/s, which is >~ 2 times more luminous than any previously known supernova. Its spectra match the hydrogen-poor sub-class of super-luminous supernovae (SLSNe-I), whose energy sources and progenitors are poorly understood. In contrast to known SLSNe-I, most of which reside in star-forming, dwarf galaxies, its host appears to be a luminous galaxy (M_V ~ -22; M_K ~ -25.1) with little star formation. In the two months since its first detection, ASASSN-15lh has radiated ~7.5x10^51 ergs, challenging the popular magnetar model for the engine of SLSNe-I.
We introduce the MUSCEL Program (MUltiwavelength observations of the Structure, Chemistry and Evolution of LSB galaxies), a project aimed at determining the star-formation histories of low surface brightness galaxies. MUSCEL utilizes ground-based optical spectra and space-based UV and IR photometry to fully constrain the star-formation histories of our targets with the aim of shedding light on the processes that led low surface brightness galaxies down a different evolutionary path from that followed by high surface brightness galaxies, such as our Milky Way. Here we present the spatially-resolved optical spectra of UGC 628, observed with the VIRUS-P IFU at the 2.7-m Harlen J. Smith Telescope at the McDonald Observatory, and utilize emission-line diagnostics to determine the rate and distribution of star formation as well as the gas-phase metallicity and metallicity gradient. We find highly clustered star formation throughout UGC 628, excluding the core regions, and a log(O/H) metallicity around -4.2, with more metal rich regions near the edges of the galactic disk. Based on the emission-line diagnostics alone, the current mode of star formation, slow and concentrated in the outer disk, appears to have dominated for quite some time, although there are clear signs of a much older stellar population formed in a more standard inside-out fashion.
In this talk we present a novel framework that unifies the stunning success of MOND on galactic scales with the triumph of the LambdaCDM model on cosmological scales. This is achieved through the rich and well-studied physics of superfluidity. The dark matter and MOND components have a common origin, representing different phases of a single underlying substance. In galaxies, dark matter thermalizes and condenses to form a superfluid phase. The superfluid phonons couple to baryonic matter particles and mediate a MOND-like force. Our framework naturally distinguishes between galaxies (where MOND is successful) and galaxy clusters (where MOND is not): dark matter has a higher temperature in clusters, and hence is in a mixture of superfluid and normal phase. The rich and well-studied physics of superfluidity leads to a number of striking observational signatures, which we briefly discuss. Remarkably the critical temperature and equation of state of the dark matter superfluid are similar to those of known cold atom systems. Identifying a precise cold atom analogue would give important insights on the microphysical interactions underlying DM superfluidity. Tantalizingly, it might open the possibility of simulating the properties and dynamics of galaxies in laboratory experiments.
We present the first spatially resolved polarized scattered light H-band detection of the DoAr 28 transitional disk. Our two epochs of imagery detect the scattered light disk from our effective inner working angle of 0.10" (13 AU) out to 0.50" (65 AU). This inner working angle is interior to the location of the system's gap inferred by previous studies using SED modeling (15 AU). We detected a candidate point source companion 1.08" northwest of the system; however, our second epoch of imagery strongly suggests that this object is a background star. We constructed a grid of Monte Carlo Radiative Transfer models of the system, and our best fit models utilize a modestly inclined (50 deg), 0.01 Msun disk that has a partially depleted inner gap from the dust sublimation radius out to ~8 AU. Subtracting this best fit, axi-symmetric model from our polarized intensity data reveals evidence for two small asymmetries in the disk, which could be attributable to variety of mechanisms.
Mass-to-light versus colour relations (MLCRs), derived from stellar population synthesis models, are widely used to estimate galaxy stellar masses (M$_*$) yet a detailed investigation of their inherent biases and limitations is still lacking. We quantify several potential sources of uncertainty, using optical and near-infrared (NIR) photometry for a representative sample of nearby galaxies from the Virgo cluster. Our method for combining multi-band photometry with MLCRs yields robust stellar masses, while errors in M$_*$ decrease as more bands are simultaneously considered. The prior assumptions in one's stellar population modelling dominate the error budget, creating a colour-dependent bias of up to 0.6 dex if NIR fluxes are used (0.3 dex otherwise). This matches the systematic errors associated with the method of spectral energy distribution (SED) fitting, indicating that MLCRs do not suffer from much additional bias. Moreover, MLCRs and SED fitting yield similar degrees of random error ($\sim$0.1-0.14 dex) when applied to mock galaxies and, on average, equivalent masses for real galaxies with M$_* \sim$ 10$^{8-11}$ M$_{\odot}$. The use of integrated photometry introduces additional uncertainty in M$_*$ measurements, at the level of 0.05-0.07 dex. We argue that using MLCRs, instead of time-consuming SED fits, is justified in cases with complex model parameter spaces (involving, for instance, multi-parameter star formation histories) and/or for large datasets. Spatially-resolved methods for measuring M$_*$ should be applied for small sample sizes and/or when accuracies less than 0.1 dex are required. An Appendix provides our MLCR transformations for ten colour permutations of the $grizH$ filter set.
We present results on the SFR-$M_*$ relation (i.e., the "main sequence") among star-forming galaxies at $1.37\leq z \leq2.61$ using the MOSFIRE Deep Evolution Field (MOSDEF) survey. Based on a sample of 261 star-forming galaxies with observations of H$\alpha$ and H$\beta$ emission lines, we have estimated robust dust-corrected instantaneous star-formation rates (SFRs) over a large dynamic range in stellar mass ($\sim 10^{9.0}-10^{11.5}M_\odot$). We find a tight correlation between SFR(H$\alpha$) and $M_*$ with an intrinsic scatter of 0.36 dex, 0.05 dex larger than that of UV-based SFRs. This increased scatter is consistent with predictions from numerical simulations of 0.03 - 0.1 dex, and is attributed to H$\alpha$ more accurately tracing SFR variations. The slope of the $\log(\text{SFR})-\log(M_*)$ relation, using SFR(H$\alpha$), at $1.4< z<2.6$ and over the stellar mass range of $10^{9.5}$ to $10^{11.5}M_\odot$ is $0.65\pm 0.09$. We find that different assumptions for the dust correction, such as using the stellar $E(B-V)$ with a Calzetti et al. (2000) attenuation curve, as well as the sample biases against red and dusty star-forming galaxies at large masses, could yield steeper slopes. Moreover, not correcting the Balmer emission line fluxes for the underlying Balmer absorption results in overestimating the dust extinction of H$\alpha$ and SFR(H$\alpha$) at the high-mass end by 2.1 (2.5) at $10^{10.6} M_\odot$ ($10^{11.1} M_\odot$) and artificially increases the slope of the main-sequence. The shallower main-sequence slope found here compared to that of galaxy evolution simulations may be indicative of different feedback processes governing the low- and/or high-mass end of the main sequence.
We examine the radio properties of the Brightest Cluster Galaxies (BCGs) in a large sample of X-ray selected galaxy clusters comprising the Brightest Cluster Sample (BCS), the extended BCS (eBCS) and ROSAT-ESO Flux Limited X-ray (REFLEX) cluster catalogues. We have multi-frequency radio observations of the BCG using a variety of data from the Australia Telescope Compact Array (ATCA), Jansky Very Large Array (VLA) and Very Long Baseline Array (VLBA) telescopes. The radio spectral energy distributions (SEDs) of these objects are decomposed into a component attributed to on-going accretion by the active galactic nuclei (AGN) that we refer to as the 'core', and a more diffuse, ageing component we refer to as the 'non-core'. These BCGs are matched to previous studies to determine whether they exhibit emission lines (principally H-alpha), indicative of the presence of a strong cooling cluster core. We consider how the radio properties of the BCGs vary with cluster environmental factors. Line emitting BCGs are shown to generally host more powerful radio sources, exhibiting the presence of a strong, distinguishable core component in about 60% of cases. This core component more strongly correlates with the BCG's [OIII]5007A line emission. For BCGs in line-emitting clusters, the X-ray cavity power correlates with both the extended and core radio emission, suggestive of steady fuelling of the AGN over bubble-rise time-scales in these clusters.
By exploiting two sets of high-resolution images obtained with HST ACS/WFC over a baseline of ~10 years we have measured relative proper motions of ~70,000 stars in the stellar system Terzan 5. The results confirm the membership of the three sub-populations with different iron abudances discovered in the system. The orbit of the system has been derived from a first estimate of its absolute proper motion, obtained by using bulge stars as reference. The results of the integration of this orbit within an axisymmetric Galactic model exclude any external accretion origin for this cluster. Terzan 5 is known to have chemistry similar to the Galactic bulge; our findings support a kinematic link between the cluster and the bulge, further strengthening the possibility that Terzan 5 is the fossil remnant of one of the pristine clumps that originated the bulge.
We introduce a new project to understand helium reionization using fully coupled $N$-body, hydrodynamics, and radiative transfer simulations. This project aims to capture correctly the thermal history of the intergalactic medium (IGM) as a result of reionization and make predictions about the Lyman-$\alpha$ forest and baryon temperature-density relation. The dominant sources of radiation for this transition are quasars, so modeling the source population accurately is very important for making reliable predictions. In this first paper, we present a new method for populating dark matter halos with quasars. Our set of quasar models include two different light curves, a lightbulb (simple on/off) and symmetric exponential model, and luminosity-dependent quasar lifetimes. Our method self-consistently reproduces an input quasar luminosity function (QLF) given a halo catalog from an $N$-body simulation, and propagates quasars through the merger history of halo hosts. After calibrating quasar clustering using measurements from BOSS, we find that the characteristic mass of quasar hosts is $M_h \sim 2.5 \times 10^{12} M_\odot$ $h^{-1}$ for the lightbulb model, and $M_h \sim 2.3 \times 10^{12} M_\odot$ $h^{-1}$ for the exponential model. In the exponential model, the peak quasar luminosity for a given halo mass is larger than that in the lightbulb model, typically by a factor of 1.5-2. The effective lifetime for quasars in the lightbulb model is 59 Myr, and in the exponential case, the effective time constant is about 15 Myr. We include semi-analytic calculations of helium reionization, and discuss how to include these quasars as sources of ionizing radiation for full hydrodynamics with radiative transfer simulations in order to study helium reionization.
We consider the high radio frequency (15 GHz - 353 GHz) properties and variability of 35 Brightest Cluster Galaxies (BCGs). These are the most core-dominated sources drawn from a parent sample of more than 700 X-ray selected clusters, thus allowing us to relate our results to the general population. We find that >6.0% of our parent sample (>15.1% if only cool-core clusters are considered) contain a radio-source at 150 GHz of at least 3mJy (~1x10^23 W/Hz at our median redshift of z~0.13). Furthermore, >3.4% of the BCGs in our parent sample contain a peaked component (Gigahertz Peaked Spectrum, GPS) in their spectra that peaks above 2 GHz, increasing to >8.5% if only cool-core clusters are considered. We see little evidence for strong variability at 15 GHz on short (week-month) time-scales although we see variations greater than 20% at 150 GHz over 6-month times-frames for 4 of the 23 sources with multi-epoch observations. Much more prevalent is long-term (year-decade time-scale) variability, with average annual amplitude variations greater than 1% at 15 GHz being commonplace. There is a weak trend towards higher variability as the peak of the GPS-like component occurs at higher frequency. We demonstrate the complexity that is seen in the radio spectra of BCGs and discuss the potentially significant implications of these high-peaking components for Sunyaev-Zel'dovich cluster searches.
We derive the transfer equations for polarized radiation in the atmospheres
of compact stars, which take into account a frequency redistribution of
radiation within and near a cyclotron line core. The equations are valid in the
magnetic fields up to $10^{13}$ G and can be used for numerical modeling of a
cyclotron line formation in the warm magnetospheric plasmas of compact stars.
We present two forms of such equations. The first form, for the intensities
of ordinary and extraordinary modes, is applicable for the compact stars with a
moderate magnetic field strength up to $10^{10}-10^{11}$ G. The second form,
for the Stokes parameters, is more complex, but applicable even if a linear
mode coupling takes place somewhere in the scattering-dominated atmosphere.
Analysing dispersion properties of a magnetized plasma, we show that the linear
mode coupling is possible for a wide range of parameters and originates from a
partial cancellation of the plasma and vacuum contributions to the refraction
indices.
An exceptional $\gamma$-ray outburst from 3C 279 is detected by {\it Fermi}-Large Area Telescope (LAT) in 2015 June. In the energy range of 0.1$-$300 GeV, the highest flux measured is (39.1$\pm$2.5) $\times$ 10$^{-6}$ \phflux, which is the highest $\gamma$-ray flux ever detected from 3C 279, exceeding the previous historically brightest flare observed by {\it EGRET} in 1996. The high activity period consists of three major flares with the last one being the brightest. All but one flares show a faster rise and slower decay pattern and at the peak of the activity, the $\gamma$-ray spectrum is found to show a clear signature of break/curvature. The obtained spectral parameters hint for the peak of the inverse Compton emission to lie in the LAT energy range (around $\sim$1 GeV) which is in contrast to that seen during the 2013 December and 2014 April $\gamma$-ray flares of 3C 279. From the $\gamma\gamma$ pair opacity arguments, the minimum Doppler factor is estimated to be 14 and the location of the $\gamma$-ray emitting region is found to be either at the outer edge of the broad line region or farther out from it.
We run hydrodynamical simulations of a 2D isothermal non self-gravitating inviscid gas flowing in a rigidly rotating externally imposed potential formed by only two components: a monopole and a quadrupole. We explore systematically the effects of varying the quadrupole while keeping fixed the monopole and discuss the consequences for the interpretation of longitude-velocity diagrams in the Milky Way. We find that the gas flow can constrain the quadrupole of the potential and the characteristics of the bar that generates it. The exponential scale length of the bar must be at least $1.5\rm\, kpc$. The strength of the bar is also constrained. Our global interpretation favours a pattern speed of $\Omega=40\,\rm km s^{-1} {kpc}^{-1}$. We find that for most observational features, there exist a value of the parameters that matches each individual feature well, but is difficult to reproduce all the important features at once. Due to the intractably high number of parameters involved in the general problem, quantitative fitting methods that can run automatic searches in parameter space are necessary.
Refined constraints on chameleon theories are calculated for atom-interferometry experiments, using a numerical approach consisting in solving for a four-region model the static and spherically symmetric Klein-Gordon equation for the chameleon field. By modeling not only the test mass and the vacuum chamber but also its walls and the exterior environment, the method allows to probe new effects on the scalar field profile and the induced acceleration of atoms. In the case of a weakly perturbing test mass, the effect of the wall is to enhance the field profile and to lower the acceleration inside the chamber by up to one order of magnitude. In the thin-shell regime, significant deviations from the analytical estimations are found, even when measurements are realized in the immediate vicinity of the test mass. Close to the vacuum chamber wall, the acceleration becomes negative and potentially measurable. This prediction could be used to discriminate between fifth-force effects and systematic experimental uncertainties, by doing the experiment at several key positions inside the vacuum chamber. The influence of the wall thickness and density is also studied. For the chameleon potential $V(\phi) = \Lambda^{4+\alpha} / \phi^\alpha$ and a coupling function $A(\phi) = \exp(\phi /M)$, one finds $M \gtrsim 7 \times 10^{16}$ GeV, independently of the power-law index. For $V(\phi) = \Lambda^4 (1+ \Lambda/ \phi)$ one finds $M \gtrsim 4 \times 10^{16}$ GeV. Future experiments able to measure an acceleration $a \sim 10^{-11} \mathrm{m/s^2}$ would probe the chameleon parameter space up to the Planck scale. Our method can easily be extended to constrain other models with a screening mechanism, such as symmetron, dilaton and f(R) theories.
The unknown nature of dark energy motivates continued cosmological tests of large-scale gravitational physics. We present a new consistency check based on the relative amplitude of non-relativistic galaxy peculiar motions, measured via redshift-space distortion, and the relativistic deflection of light by those same galaxies traced by galaxy-galaxy lensing. We take advantage of the latest generation of deep, overlapping imaging and spectroscopic datasets, combining the Red Cluster Sequence Lensing Survey (RCSLenS), the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the WiggleZ Dark Energy Survey and the Baryon Oscillation Spectroscopic Survey (BOSS). We quantify the results using the "gravitational slip" statistic E_G, which we estimate as 0.48 +/- 0.10 at z=0.32 and 0.30 +/- 0.07 at z=0.57, the latter constituting the highest redshift at which this quantity has been determined. These measurements are consistent with the predictions of General Relativity, for a perturbed Friedmann-Robertson-Walker metric in a Universe dominated by a cosmological constant, which are E_G = 0.41 and 0.36 at these respective redshifts. The combination of redshift-space distortion and gravitational lensing data from current and future galaxy surveys will offer increasingly stringent tests of fundamental cosmology.
We present high-precision photometric light curves of five O-type stars observed with the refurbished {\it Kepler\/} satellite during its Campaign 0. For one of the stars, we also assembled high-resolution ground-based spectroscopy with the {\sc hermes} spectrograph attached to the 1.2-m Mercator telescope. The stars EPIC202060097 (O9.5V) and EPIC202060098 (O7V) exhibit monoperiodic variability due to rotational modulation with an amplitude of 5.6 mmag and 9.3 mmag and a rotation period of 2.63 d and 5.03 d, respectively. EPIC202060091 (O9V) and EPIC202060093 (O9V:pe) reveal variability at low frequency but the cause is unclear. EPIC202060092 (O9V:p) is discovered to be a spectroscopic binary with at least one multiperiodic $\beta\,$Cep-type pulsator whose detected mode frequencies occur in the range $[0.11,6.99]$ d$^{-1}$ and have amplitudes between 0.8 and 2.0 mmag. Its pulsation spectrum is shown to be fully compatible with the ones predicted by core-hydrogen burning O-star models. Despite the short duration of some 33\,d and the limited data quality with a precision near 100 $\mu$mag of these first K2 data, the diversity of possible causes for O-star variability already revealed from campaigns of similar duration by the MOST and CoRoT satellites is confirmed with {\it Kepler}. We provide an overview of O-star space photometry and give arguments why future K2 monitoring during Campaigns 11 and 13 at short cadence, accompanied by time-resolved high-precision high-resolution spectroscopy opens up the possibility of in-depth O-star seismology.
Context. Edge-on galaxies can offer important insights in galaxy evolution as
they are the only systems where the distribution of the different components
can be studied both radially and vertically. The HEROES project was designed to
investigate the interplay between the gas, dust, stars and dark matter (DM) in
a sample of 7 massive edge-on spiral galaxies.
Aims. In this second HEROES paper we present an analysis of the atomic gas
content of 6 out of 7 galaxies in our sample. The remaining galaxy was recently
analysed according to the same strategy. The primary aim of this work is to
constrain the surface density distribution, the rotation curve and the geometry
of the gas disks in a homogeneous way. In addition we identify peculiar
features and signs of recent interactions.
Methods. We construct detailed tilted-ring models of the atomic gas disks
based on new GMRT 21-cm observations of NGC 973 and UGC 4277 and re-reduced
archival HI data of NGC 5907, NGC 5529, IC 2531 and NGC 4217. Potential
degeneracies between different models are resolved by requiring a good
agreement with the data in various representations of the data cubes.
Results. From our modelling we find that all but one galaxy are warped along
the major axis. In addition, we identify warps along the line of sight in three
galaxies. A flaring gas layer is required to reproduce the data only for one
galaxy, but (moderate) flares cannot be ruled for the other galaxies either. A
coplanar ring-like structure is detected outside the main disk of NGC 4217,
which we suggest could be the remnant of a recent minor merger event. We also
find evidence for a radial inflow of 15 +- 5 km/s in the disk of NGC 5529,
which might be related to the ongoing interaction with two nearby companions.
(Abridged)
We analyze a set of volume limited samples from SDSS DR12 to quantify the degree of inhomogeneity at different length scales using Shannon entropy. We find that the galaxy distributions exhibit a higher degree of inhomogeneity as compared to a Poisson point process at all length scales. Our analysis indicates that signatures of inhomogeneities in the galaxy distributions persist at least upto a length scale of $120 \, h^{-1}\, {\rm Mpc}$. The galaxy distributions appear to be homogeneous on a scale of $140 \, h^{-1}\, {\rm Mpc}$ and beyond. Analyzing a set of mock galaxy samples from a semi analytic galaxy catalogue from the Millennium simulation we find a scale of transition to homogeneity at $\sim 100 \, h^{-1}\, {\rm Mpc}$.
It is well known that, for bright gamma-ray pulsars with high statistics above a few GeV, the phase averaged spectral energy distribution (SED) is harder than a simple exponential cutoff above the break. We perform phase-resolved spectral analyses of bright gamma-ray pulsars and demonstrate that, even over narrow phase ranges, the SEDs of gamma-ray pulsars above the break energy are harder than a simple exponential cutoff. We argue within a radiation-reaction limited curvature framework that this is indicative of non-stationary emission or emission from multiple zones. Further, we address a common problem faced when fitting hard spectral tails with a power-law times a sub-exponential function. Namely, that the sub-exponent parameter does not describe any parameters of physical models of pulsar emission. We introduce a simple analytical fit function to solve this problem.
DG CVn is a binary system in which one of the components is an M type dwarf ultra fast rotator, only three of which are known in the solar neighborhood. Observations of DG CVn by the Swift satellite and several ground-based observatories during its super-flare event on 2014 allowed us to perform a complete hard X-ray - optical follow-up of a super-flare from the red-dwarf star. The observations support the fact that the super-flare can be explained by the presence of (a) large active region(s) on the surface of the star. Such activity is similar to the most extreme solar flaring events. This points towards a plausible extrapolation between the behaviour from the most active red-dwarf stars and the processes occurring in the Sun.
Basic ideas of muon tracker technique for the solution of primary cosmic ray (PCR) composition problem in the energy range $10^{17}-10^{18}$ eV are presented. The approach uses MC simulation data made with CORSIKA6.990 for "Pamir-XXI" site conditions. Similar technology can certainly be developed for other observation levels and interaction models. One can probably extend it to much higher primary energies.
Nonhelical shear dynamos are studied with a particular focus on the possibility of coherent dynamo action. The primary results -- serving as a follow up to the results of Squire & Bhattacharjee [arXiv:1506.04109 (2015)] -- pertain to the "magnetic shear-current effect" as a viable mechanism to drive large-scale magnetic field generation. This effect raises the interesting possibility that the saturated state of the small-scale dynamo could drive large-scale dynamo action, and is likely to be important in the unstratified regions of accretion disk turbulence. In this paper, the effect is studied at low Reynolds numbers, removing the complications of small-scale dynamo excitation and aiding analysis by enabling the use of quasi-linear statistical simulation methods. In addition to the magnetically driven dynamo, new results on the kinematic nonhelical shear dynamo are presented. These illustrate the relationship between coherent and incoherent driving in such dynamos, demonstrating the importance of rotation in determining the relative dominance of each mechanism.
We present a stacking analysis of the complete sample of Early Type Galaxies (ETGs) in the \textit{Chandra} COSMOS (C-COSMOS) survey, to explore the nature of the X-ray luminosityin the redshift and stellar luminosity ranges \(0<z<1.5\) and \({10}^{9}<L_K/L_{\astrosun}<{10}^{13}\). Using established scaling relations, we subtract the contribution of X-ray binary populations, to estimate the combined emission of hot ISM and AGN. To discriminate between the relative importance of these two components, we (1) compare our results with the relation observed in the local universe \(L_{X,gas}\propto L_K^{4.5}\) for hot gaseous halos emission in ETGs, and (2) evaluate the spectral signature of each stacked bin. We find two regimes where the non-stellar X-ray emission is hard, consistent with AGN emission. First, there is evidence of hard, absorbed X-ray emission in stacked bins including relatively high z (\(\sim 1.2\)) ETGs with average high X-ray luminosity (\(L_{X-LMXB}\gtrsim 6\times{10}^{42}\mbox{ erg}/\mbox{s}\)). These luminosities are consistent with the presence of of highly absorbed "hidden" AGNs in these ETGs, which are not visible in their optical-IR spectra and spectral energy distributions. Second, confirming the early indication from our C-COSMOS study of X-ray detected ETGs, we find significantly enhanced X-ray luminosity in lower stellar mass ETGs (\(L_K\lesssim{10}^{11}L_{\astrosun}\)), relative to the local \(L_{X,gas}\propto L_K^{4.5}\) relation. The stacked spectra of these ETGs also suggest X-ray emission {harder than expected from gaseous hot halos}. This emission is consistent with inefficient accretion \({10}^{-5}-{10}^{-4}\dot{M}_{Edd}\) onto \(M_{BH}\sim {10}^{6}-{10}^{8}\,M_{\astrosun}\).
We report the discovery of a bright blue quasar: SDSS J022218.03-062511.1. This object was discovered spectroscopically while searching for hot white dwarfs that may be used as calibration sources for large sky surveys such as the Dark Energy Survey or the Large Synoptic Survey Telescope project. We present the calibrated spectrum, spectral line shifts and report a redshift of z = 0.521 +/- 0.0015 and a rest-frame g-band luminosity of 8.71 X 10^11 L(Sun).
In this paper we present multi-band period-luminosity (P-L) relations for fundamental mode Cepheids in the SMC. The optical VI-band mean magnitudes for these SMC Cepheids were taken from the third phase of the Optical Gravitational Lensing Experiment (OGLE-III) catalog. We also matched the OGLE-III SMC Cepheids to 2MASS and SAGE-SMC catalog to derive mean magnitudes in the JHK-bands and the four {\it Spitzer} IRAC bands, respectively. All photometry was corrected for extinction by adopting the Zaritsky's extinction map. Cepheids with periods smaller than $\sim2.5$ days were removed from the sample. In addition to the extinction corrected P-L relations in nine filters from optical to infrared, we also derived the extinction-free Wesenheit function for these Cepheids. We tested the nonlinearity of these SMC P-L relations (except the $8.0\mu\mathrm{m}$-band P-L relation) at 10 days: none of the P-L relations show statistically significant evidence of nonlinearity. When compared to the P-L relations in the LMC, the t-test results revealed that there is a difference between the SMC/LMC P-L slopes only in the V- and J-band. Further, we found excellent agreement between the SMC/LMC Wesenheit P-L slope. The difference in LMC and SMC Period-Wesenheit relation LMC and SMC zero points was found to be $\Delta \mu=0.483\pm0.015$ mag. This amounts to a difference in distance modulus between the LMC and SMC.
Sun-like stellar oscillations are excited by turbulent convection and have been discovered in some 500 main sequence and sub-giant stars and in more than 12,000 red giant stars. When such stars are near gravitational wave sources, low-order quadrupole acoustic modes are also excited above the experimental threshold of detectability, and they can be observed, in principle, in the acoustic spectra of these stars. Such stars form a set of natural detectors to search for gravitational waves over a large spectral frequency range, from $10^{-7}$ Hz to $10^{-2}$ Hz. In particular, these stars can probe the $10^{-6}$ Hz -- $10^{-4}$ Hz spectral window which cannot be probed by current conventional gravitational wave detectors, such as SKA and eLISA. The PLATO stellar seismic mission will achieve photospheric velocity amplitude accuracy of $~ {\rm cm/s}$. For a gravitational wave search, we will need to achieve accuracies of the order of $10^{-2}{\rm cm/s}$, i.e., at least one generation beyond PLATO. However, we have found that multi-body stellar systems have the ideal setup for this type of gravitational wave search. This is the case for triple stellar systems formed by a compact binary and an oscillating star. Continuous monitoring of the oscillation spectra of these stars to a distance of up to a kpc could lead to the discovery of gravitational waves originating in our galaxy or even elsewhere in the universe. Moreover, unlike experimental detectors, this observational network of stars will allow us to study the progression of gravitational waves throughout space.
A pair of curved shocks in a collisionless plasma is examined with a two-dimensional particle-in-cell (PIC) simulation. The shocks are created by the collision of two electron-ion clouds at a speed that exceeds everywhere the threshold speed for shock formation. A variation of the collision speed along the initially planar collision boundary, which is comparable to the ion acoustic speed, yields a curvature of the shock that increases with time. The spatially varying Mach number of the shocks results in a variation of the downstream density in the direction along the shock boundary. This variation is eventually equilibrated by the thermal diffusion of ions. The pair of shocks is stable for tens of inverse ion plasma frequencies. The angle between the mean flow velocity vector of the inflowing upstream plasma and the shock's electrostatic field increases steadily during this time. The disalignment of both vectors gives rise to a rotational electron flow, which yields the growth of magnetic field patches that are coherent over tens of electron skin depths.
Using photometry and high resolution spectroscopy we investigate for the first time the physical connection between the open clusters NGC 5617 and Trumpler 22. Based on new CCD photometry we report their spatial proximity and common age of ~70 Myr. Based on high resolution spectra collected using the HERMES and UCLES spectrographs on the Anglo-Australian telescope, we present radial velocities and abundances for Fe, Na, Mg, Al, Si, Ca and Ni. The measured radial velocities are -38.63 +/-2.25 km/s for NGC 5617 and -38.46 +/-2.08 km/s for Trumpler 22. The mean metallicity of NGC 5617 was found to be [Fe/H] =-0.18 +/-0.02 and for Trumpler 22 was found to be [Fe/H] = -0.17 +/-0.04. The two clusters share similar abundances across the other elements, indicative of a common chemical enrichment history of these clusters. Together with common motions and ages we confirm that NGC 5617 and Trumpler 22 are a primordial binary cluster pair in the Milky Way.
We analyze the X-ray spectra of the atoll 4U~1705-44 when the source undergoes the island-banana state transition. We use the Rossi X-ray Timing Explorer (RXTE) and BeppoSAX observations for this analysis. We demonstrate that the broad-band energy spectral distributions for all evolutinary states can be fitted by a model, consisting two Comptonized components. One arises from the seed photons coming from a neutron star (NS) atmosphere at a temperature kT_{s1}<1.5 keV (herein Comptb) and a second resulting from the seed photons of T_{s2}~1.1-1.3 keV coming from the disk (herein Comptb2). We found that we needed to add a low-temperature blackbody and an iron-line ({Gaussian}) component to the model in order to obtain high-quality fits. The data analysis using this model indicates that the power-law photon index Gamma_{1} of our model is always about 2, independently of the spectral state. Another parameter, Gamma_{2} demonstrates a two-phase behavior depending on the spectral state. Gamma_{2} is quasi-constant at Gamma_{2}~ 2 when the electron temperature kT^{(2)}_e<80 keV and Gamma_{2} is less than 2, in the range of 1.3<Gamma_{2}<2, when kT^{(2)}_e>80 keV. This phase is similar to that was previously found in the Z-source Sco X-1. We interpret the decreasing index phase using a model in which a super-Eddington radiation pressure from the neutron star causes an expansion of the Compton cloud similar to that found previously in Sco~X-1 during the Flaring branch.
We present a dynamical analysis of the galaxy cluster AC114 based on a catalogue of 524 velocities. Of these, 169 (32%) are newly obtained at ESO (Chile) with the VLT and the VIMOS spectrograph. Data on individual galaxies are presented and the accuracy of the measured velocities is discussed. Dynamical properties of the cluster are derived. We obtain an improved mean redshift value z= 0.31665 +/- 0.0008 and velocity dispersion \sigma= 1893+73-82 \kms. A large velocity dispersion within the core radius and the shape of the infall pattern suggests that this part of the cluster is in a radial phase of relaxation with a very elongated radial filament spanning 12000 \kms. A radial foreground structure is detected within the central 0.5/h Mpc radius, recognizable as a redshift group at the same central redshift value. We analyze the color distribution for this archetype Butcher-Oemler galaxy cluster and identify the separate red and blue galaxy sequences. The latter subset contains 44% of confirmed members of the cluster, reaching magnitudes as faint as R_{f}= 21.1 (1.0 magnitude fainter than previous studies). We derive a mass M_{200}= (4.3 \pm 0.7) x 10^15 Msun/h. In a subsequent paper we will utilize the spectral data presented here to explore the mass-metallicity relation for this intermediate redshift cluster.
The \textit{Spitzer}/Infrared Spectrograph spectra of three spectroscopically anomalous galaxies (IRAS~F10398+1455, IRAS~F21013-0739 and SDSS~J0808+3948) are modeled in terms of a mixture of warm and cold silicate dust, and warm and cold carbon dust. Their unique infrared (IR) emission spectra are characterized by a steep $\simali$5--8$\mum$ emission continuum, strong emission bands from polycyclic aromatic hydrocarbon (PAH) molecules, and prominent silicate emission. The steep $\simali$5--8$\mum$ emission continuum and strong PAH emission features suggest the dominance of starbursts, while the silicate emission is indicative of significant heating from active galactic nuclei (AGNs). With warm and cold silicate dust of various compositions ("astronomical silicate," amorphous olivine, or amorphous pyroxene) combined with warm and cold carbon dust (amorphous carbon, or graphite), we are able to closely reproduce the observed IR emission of these %spectroscopically anomalous galaxies. We find that the dust temperature is the primary cause in regulating the steep $\sim$5--8$\mum$ continuum and silicate emission, insensitive to the exact silicate or carbon dust mineralogy and grain size $a$ as long as $a\simlt1\mum$. More specifically, the temperature of the $\simali$5--8$\mum$ continuum emitter (which is essentially carbon dust) of these galaxies is $\sim$250--400$\K$, much lower than that of typical quasars which is $\sim$640$\K$. Moreover, it appears that larger dust grains are preferred in quasars. The lower dust temperature and smaller grain sizes inferred for these three galaxies compared with that of quasars could be due to the fact that they may harbor a young/weak AGN which is not maturely developed yet.
We have studied the relationship between the solar-wind speed $[V]$ and the coronal magnetic-field properties (a flux expansion factor [$f$] and photospheric magnetic-field strength [$B_{\mathrm{S}}$]) at all latitudes using data of interplanetary scintillation and solar magnetic field obtained for 24 years from 1986 to 2009. Using a cross-correlation analyses, we verified that $V$ is inversely proportional to $f$ and found that $V$ tends to increase with $B_{\mathrm{S}}$ if $f$ is the same. As a consequence, we find that $V$ has extremely good linear correlation with $B_{\mathrm{S}}/f$. However, this linear relation of $V$ and $B_{\mathrm{S}}/f$ cannot be used for predicting the solar-wind velocity without information on the solar-wind mass flux. We discuss why the inverse relation between $V$ and $f$ has been successfully used for solar-wind velocity prediction, even though it does not explicitly include the mass flux and magnetic-field strength, which are important physical parameters for solar-wind acceleration.
We report the first detection of linearly polarized emission at an observing wavelength of 350~$\mu$m from the radio-loud Active Galactic Nucleus (AGN) 3C~279. We conducted polarization observations for 3C~279 using the SHARP polarimeter in the Caltech Submillimeter Observatory (CSO) on 2014 March 13 and 14. For the first time, we detected the linear polarization with the degree of polarization of 13.3\%$\pm$3.4\% {\bf($3.9\sigma$)} and the Electric Vector Position Angle (EVPA) of 34.7$^\circ\pm5.6^\circ$. We also observed 3C~279 simultaneously at 22, 43, and 86~GHz in dual polarization with the Korean VLBI Network (KVN) on 2014 March 6 (single dish) and imaged in milliarcsecond (mas) scales at 22, 43, 86, and 129~GHz on March 22 (VLBI). We found that the degree of linear polarization increases from 10\% to 13\% at 22~GHz to 350~$\mu$m and the EVPAs at all observing frequencies are parallel within $<10^\circ$ to the direction of the jet at mas scale, implying that the integrated magnetic fields are perpendicular to the jet in the innermost regions. We also found that the Faraday rotation measures RM are in a range of $-6.5\times10^2 \sim -2.7\times10^3$~rad~m$^{-2}$ between 22-86~GHz, and are scaled as a function of wavelength: $|{\rm RM}|\propto\lambda^{-2.2}$. These results indicate that the mm and sub-mm polarization emission are generated in the compact jet within 1~mas scale and affected by a Faraday screen in or in the close proximity of the jet.
We report on the development and construction of the high-purity germanium spectrometer setup GIOVE (Germanium Inner Outer Veto), recently built and now operated at the shallow underground laboratory of the Max-Planck-Institut f\"ur Kernphysik, Heidelberg. Particular attention was paid to the design of a novel passive and active shield, aiming at efficient rejection of environmental and muon induced radiation backgrounds. The achieved sensitivity level of <100 {\mu}Bq/kg for primordial radionuclides from U and Th in typical {\gamma} ray sample screening measurements is unique among instruments located at comparably shallow depths and can compete with instruments at far deeper underground sites.
We report results of investigation of amplitude calibration for very long baseline interferometry (VLBI) observations with Korean VLBI Network (KVN). Amplitude correction factors are estimated based on comparison of KVN observations at 22~GHz correlated by Daejeon hardware correlator and DiFX software correlator in Korea Astronomy and Space Science Institute (KASI) with Very Long Baseline Array (VLBA) observations at 22~GHz by DiFX software correlator in National Radio Astronomy Observatory (NRAO). We used the observations for compact radio sources, 3C~454.3 and NRAO~512, which are almost unresolved for baselines in a range of 350-477~km. Visibility data of the sources obtained with similar baselines at KVN and VLBA are selected, fringe-fitted, calibrated, and compared for their amplitudes. We found that visibility amplitudes of KVN observations should be corrected by factors of 1.10 and 1.35 when correlated by DiFX and Daejeon correlators, respectively. These correction factors are attributed to the combination of two steps of 2-bit quantization in KVN observing systems and characteristics of Daejeon correlator.
The Large and Small Magellanic Clouds (LMC and SMC, respectively) are observed to have characteristic dust extinction curves that are quite different from those of the Galaxy (e.g., strength of the 2175 A bump). Although the dust composition and size distribution of the Magellanic Clouds (MCs) that can self-consistently explain their observed extinction curves have been already proposed, it remain unclear whether and how the required dust properties can be achieved in the formation histories of the MCs. We therefore investigate the time evolution of the dust properties of the MCs and thereby derive their extinction curves using one-zone chemical evolution models with formation and evolution of small and large silicate and carbonaceous dust grains and dusty winds associated with starburst events. We find that the observed SMC extinction curve without a conspicuous 2175 A bump can be reproduced well by our SMC model, if the small carbon grains can be selectively lost through the dust wind during the latest starburst about 0.2 Gyr ago. We also find that the LMC extinction curve with a weak 2175 A bump can be reproduced by our LMC model with less efficient removal of dust through dust wind. We discuss possible physical reasons for different dust wind efficiencies between silicate and graphite and among galaxies.
van Dokkum and Conroy revisited the unexpectedly strong Na I lines at 8200 A found in some giant elliptical galaxies and interpreted it as evidence for unusually bottom-heavy initial mass function. Jeong et al. later found a large population of galaxies showing equally-extraordinary Na D doublet absorption lines at 5900 A (Na D excess objects: NEOs) and showed that their origins can be different for different types of galaxies. While a Na D excess seems to be related with the interstellar medium (ISM) in late-type galaxies, smooth-looking early-type NEOs show little or no dust extinction and hence no compelling sign of ISM contributions. To further test this finding, we measured the doppler components in the Na D lines. We hypothesized that ISM would have a better (albeit not definite) chance of showing a blueshift doppler departure from the bulk of the stellar population due to outflow caused by either star formation or AGN activities. Many of the late-type NEOs clearly show blueshift in their Na D lines, which is consistent with the former interpretation that the Na D excess found in them is related with star formation-caused gas outflow. On the contrary, smooth-looking early-type NEOs do not show any notable doppler component, which is also consistent with the interpretation of Jeong et al. that the Na D excess in early-type NEOs is likely not related with ISM activities but is purely stellar in origin.
Baldwin (1982) wrote that "the distribution of sources in the radio luminosity, P, overall physical size, D, diagram" could be considered as "the radio astronomer's H-R diagram". However, unlike the case of stars, not only the intrinsic properties of the jets, but also those of the host galaxy and the intergalactic medium are relevant to explain the evolutionary tracks of radio radio sources. In this contribution I review the current status of our understanding of the evolution of radio sources from a theoretical and numerical perspective, using the P-D diagram as a framework. An excess of compact (linear size < 10 kpc) sources could be explained by low-power jets being decelerated within the host galaxy, as shown by recent numerical simulations. These decelerated jets could also explain the population of the radio sources that have been recently classified as FR0. I will discuss the possible tracks that radio sources may follow within this diagram, and some of the physical processes that can explain the different tracks.
We present multilevel radiative transfer modeling of the scattering polarization observed in the solar O I infrared triplet around 777 nm. We demonstrate that the scattering polarization pattern observed on the solar disk forms in the chromosphere, far above the photospheric region where the bulk of the emergent intensity profiles originates. We study the sensitivity of the polarization pattern to the thermal structure of the solar atmosphere and to the presence of weak magnetic fields (0.01 - 100 G) through the Hanle effect, showing that the scattering polarization signals of the oxygen infrared triplet encode information on the magnetism of the solar chromosphere.
Stellar conditions leading to a possible semi-convective mixing are discussed in three relevant cases: (1) low mass MS stars in which the CNO cycle takes progressively the lead over the PP chain due to the increase in temperature as core hydrogen burning proceeds, (2) massive MS stars which experience a large contri- bution of the radiation pressure to the total pressure and (3) core helium burning stars for which the production of carbon in the core increases the opacity. A short discussion of semi-convection in terms of instability of non radial modes follows.
We present a nucleosynthesis sensitivity study for the $\gamma$-process in a Supernova type II model within the NuGrid research platform. The simulations aimed at identifying the relevant local production and destruction rates for the p-nuclei of molybdenum and at determining the sensitivity of the final abundances to these rates. We show that local destruction rates strongly determine the abundance of $^{92}$Mo and $^{94}$Mo, and quantify the impact.
We discuss properties of oscillatory convective modes in low-mass red giants, and compare them with observed properties of the long secondary periods (LSPs) of semi-regular red giant variables. Oscillatory convective modes are very nonadiabatic g$^{-}$ modes and they are present in luminous stars, such as red giants with $\log L/{\rm L}_\odot \ga 3$. Finite amplitudes for these modes are confined to the outermost nonadiabatic layers, where the radiative energy flux is more important than the convective energy flux. The periods of oscillatory convection modes increase with luminosity, and the growth times are comparable to the oscillation periods. The LSPs of red giants in the Large Magellanic Cloud (LMC) are observed to lie on a distinct period-luminosity sequence called sequence D. This sequence D period-luminosity relation is roughly consistent with the predictions for dipole oscillatory convective modes in AGB models if we adopt a mixing length of 1.2 pressure scale height ($\alpha = 1.2$). However, the effective temperature of the red-giant sequence of the LMC is consistent to models with $\alpha=1.9$, which predict periods too short by a factor of two.
The use of acoustic holography in the high-frequency $p$-mode spectrum can resolve the source distributions of enhanced acoustic emissions within halo structures surrounding active regions. In doing so, statistical methods can then be applied to ascertain relationships with the magnetic field. This is the focus of this study. The mechanism responsible for the detected enhancement of acoustic sources around solar active regions has not yet been explained. Furthermore the relationship between the magnetic field and enhanced acoustic emission has not yet been comprehensively examined. We have used vector magnetograms from the \Helioseismic and Magnetic Imager (HMI) on-board the Solar Dynamics Observatory (SDO) to image the magnetic-field properties in the halo. We have studied the acoustic morphology of an active region, with a complex halo and "glories," and we have linked some acoustic properties to the magnetic-field configuration. In particular, we find that acoustic sources are significantly enhanced in regions of intermediate field strength with inclinations no different from the distributions found in the quiet Sun. Additionally we have identified a transition region between the active region and the halo, in which the acoustic source power is hindered by inclined fields of intermediate field strength. Finally, we have compared the results of acoustic emission maps, calculated from holography, and the commonly used local acoustic maps, finding that the two types of maps have similar properties with respect to the magnetic field but lack spatial correlation when examining the highest-powered regions.
The dynamical evolution of the refractive index of the tensor modes of the geometry produces a specific class of power spectra characterized by a blue (i.e. slightly increasing) slope which is directly determined by the competition of the slow-roll parameter and of the rate of variation of the refractive index. Throughout the conventional stages of the inflationary and post-inflationary evolution, the microwave background anisotropies measurements, the pulsar timing limits and the big-bang nucleosythesis constraints set stringent bounds on the refractive index and on its rate of variation. Within the physically allowed region of the parameter space the cosmic background of relic gravitons leads to a potentially large signal for the ground based detectors (in their advanced version) and for the proposed space-borne interferometers. Conversely, the lack of direct detection of the signal will set a qualitatively new bound on the dynamical variation of the refractive index.
This paper contains the result of the elaboration of informations about Solar Irradiation, Geological, Meteorological and Climatic from the point of view of the quantitative data and social interaction recorded in Taviano (LE) over 24 years. These data are compared to check local variations, long term trends, and correlation with mean annual temperature. The ultimate goal of this work is to understand long term climatic changes in this geographic area. The classes of event considerated are hydrogeological phenomena, sun irradiation, seismic, volcanic, meteorological and climatological event. Only event occurred between 1990 and 2014 are considerated. The analysis is performed using a statistical approach. A particular care is used to minimize any effect due to prejudices in case of lack of data. Finally, we calculate the annual average from the monthly ones. Data on this paper don't come from a complete census of phenomena; they are considered enough representative of the accepted vulnerability level at the beginning of this study.
Bar-induced gas inflows towards the galaxy centres are recognized as a key agent for the secular evolution of galaxies. One immediate consequence is the accumulation of gas in the centre of galaxies where it can form stars and alter the chemical and physical properties. We use a sample of nearby face--on disc galaxies with available SDSS spectra to study whether the properties of the ionised gas in the central parts (radii <~0.6-2.1 kpc) of barred galaxies are altered by the presence of a bar, and whether the bar effect is related to bar and/or parent galaxy properties. The distributions of all parameters analysed are different for barred and unbarred galaxies, except for the R23 metallicity tracer and the oxygen abundance (from photoionisation models). The median values point towards (marginally) larger dust content, star formation rate per unit area, electron density and ionisation parameter in the centres of barred galaxies than in the unbarred counterpart. The most remarkable barred/unbarred difference appears in the [NII]6583/Ha line ratio, which is on average ~25% larger in barred galaxies, due to a larger N/O in the centres of these galaxies. We observe an enhancement of the central gas differences in later-type galaxies or galaxies with less massive bulges. However the bar seems to have a lower impact on the central gas properties for galaxies with more massive bulges (M_bulge > 10^10 M_sun) or galaxies with total stellar mass above ~ 10^10.8 M_sun. In conclusion, we find observational evidence that the presence of a galactic bar affects the central ionised gas properties of disc galaxies, where the most striking effect is an enhancement in the N/O abundance ratio, which can be qualitatively interpreted as due to a different origin or evolutionary processes for less and more massive bulges, with the gaseous phase of the former having currently a closer relation with bars.
The Radiation Assessment Detector (RAD), on board Mars Science Laboratory's (MSL) rover Curiosity, measures the {energy spectra} of both energetic charged and neutral particles along with the radiation dose rate at the surface of Mars. With these first-ever measurements on the Martian surface, RAD observed several effects influencing the galactic cosmic ray (GCR) induced surface radiation dose concurrently: [a] short-term diurnal variations of the Martian atmospheric pressure caused by daily thermal tides, [b] long-term seasonal pressure changes in the Martian atmosphere, and [c] the modulation of the primary GCR flux by the heliospheric magnetic field, which correlates with long-term solar activity and the rotation of the Sun. The RAD surface dose measurements, along with the surface pressure data and the solar modulation factor, are analysed and fitted to empirical models which quantitatively demonstrate} how the long-term influences ([b] and [c]) are related to the measured dose rates. {Correspondingly we can estimate dose rate and dose equivalents under different solar modulations and different atmospheric conditions, thus allowing empirical predictions of the Martian surface radiation environment.
After seven years of science operation, the Fermi mission has brought great advances in the study of Gamma-ray Bursts (GRBs). Over 1600 GRBs have been detected by the Gamma-ray Burst Monitor, and more than 100 of these are also detected by the Large Area Telescope above 30 MeV. We will give an overview of these observations, presenting the common properties in the GRB temporal and spectral behavior at high energies. We will also highlight the unique characteristics of some individual bursts. The main physical implications of these results will be discussed, along with open questions regarding GRB modeling in their prompt and temporally-extended emission phases.
Revised spectroscopic parameters for the HF molecule and a new CN line list in the 2.3 mu region have been recently available, allowing a revision of the F content in AGB stars. AGB carbon stars are the only observationally confirmed sources of fluorine. Nowadays there is not a consensus on the relevance of AGB stars in its Galactic chemical evolution. The aim of this article is to better constrain the contribution of these stars with a more accurate estimate of their fluorine abundances. Using new spectroscopic tools and LTE spectral synthesis, we redetermine fluorine abundances from several HF lines in the K-band in a sample of Galactic and extragalactic AGB carbon stars of spectral types N, J and SC spanning a wide range of metallicities. On average, the new derived fluorine abundances are systematically lower by 0.33 dex with respect to previous determinations. This may derive from a combination of the lower excitation energies of the HF lines and the larger macroturbulence parameters used here as well as from the new adopted CN line list. Yet, theoretical nucleosynthesis models in AGB stars agree with the new fluorine determinations at solar metallicities. At low metallicities, an agreement between theory and observations can be found by handling in a different way the radiative/convective interface at the base of the convective envelope. New fluorine spectroscopic measurements agree with theoretical models at low and at solar metallicity. Despite this, complementary sources are needed to explain its observed abundance in the solar neighbourhood.
Cosmic rays up to at least PeV energies are usually described in the framework of an elementary scenario that involves acceleration by objects that are located in the disk of the Milky Way, such as supernova remnants or massive star-forming regions, and then diffusive propagation throughout the Galaxy. Details of the propagation process are so far inferred mainly from the composition of cosmic rays measured near the Earth and then extrapolated to the whole Galaxy. The details of the propagation in the Galactic halo and the escape into the intergalactic medium remain uncertain. The densities of cosmic rays in specific locations can be traced via the gamma rays they produce in inelastic collisions with clouds of interstellar gas. Therefore, we analyze 73 months of Fermi-LAT data from 300 MeV to 10 GeV in the direction of several high- and intermediate-velocity clouds that are located in the halo of the Milky Way. These clouds are supposed to be free of internal sources of cosmic rays and hence any gamma-ray emission from them samples the large-scale distribution of Galactic cosmic rays. We evaluate for the first time the gamma-ray emissivity per hydrogen atom up to ~7 kpc above the Galactic disk. The emissivity is found to decrease with distance from the disk, which provides direct evidence that cosmic rays at the relevant energies originate therein. Furthermore, the emissivity of one of the targets, the upper intermediate-velocity Arch, hints at a 50% decline of the cosmic-ray intensity within 2 kpc from the disk.
The dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter-dominated objects known. Due to their proximity, high dark matter content, and lack of astrophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most promising targets for the indirect detection of dark matter via gamma rays. Here we report on gamma-ray observations of Milky Way dwarf spheroidal satellite galaxies based on 6 years of Fermi Large Area Telescope data processed with the new Pass 8 reconstruction and event-level analysis. None of the dwarf galaxies are significantly detected in gamma rays, and we present upper limits on the dark matter annihilation cross section from a combined analysis of the 15 most promising dwarf galaxies. The constraints derived are among the strongest to date using gamma rays and lie below the canonical thermal relic cross section for WIMPs of mass $\lesssim 100~GeV$ annihilating via the $b \bar b$ and $\tau^{+}\tau^{-}$ channels.
Recent developments of the nuclear emulsion technology led to the production of films with nanometric silver halide grains suitable to track low energy nuclear recoils with submicrometric length. This improvement opens the way to a directional Dark Matter detection, thus providing an innovative and complementary approach to the on-going WIMP searches. An important background source for these searches is represented by neutron-induced nuclear recoils that can mimic the WIMP signal. In this paper we provide an estimation of the contribution to this background from the intrinsic radioactive contamination of nuclear emulsions. We also report the induced background as a function of the read-out threshold, by using a GEANT4 simulation of the nuclear emulsion, showing that it amounts to about 0.02 neutrons per year per kilogram, fully compatible with the design of a 10 kg$\times$year exposure.
In X-ray binary star systems consisting of a compact object that accretes material from an orbiting secondary star, there is no straightforward means to decide if the compact object is a black hole or a neutron star. To assist this process we develop a Bayesian statistical model which makes use of the fact that X-ray binary systems appear to cluster based on their compact object type when viewed from a 3-dimensional coordinate system derived from X-ray spectral data, where the first coordinate is the ratio of counts in mid to low energy band (color 1), the second coordinate is the ratio of counts in high to low energy band (color 2), and the third coordinate is the sum of counts in all three bands. Precisely, we use this model to estimate the probabilities that an X-ray binary system contains a black hole, non-pulsing neutron star or pulsing neutron star. In particular we utilize a latent variable model in which the latent variables follow a Gaussian process prior distribution, and hence we are able to induce the spatial correlation we believe exists between systems of the same type. The utility of this approach is evidenced by the accurate prediction of system types using Rossi X-ray Timing Explorer All Sky Monitor data, but it is not flawless. In particular, non-pulsing neutron systems containing "bursters" which are close to the boundary demarcating systems containing black holes tend to be classified as black hole systems. As a byproduct of our analyses, we provide the astronomer with public R code that can be used to predict the compact object type of X-ray binaries given training data.
Strongly-interacting matter in the form of nuggets of nuclear-density material are not currently excluded as dark matter candidates in the ten gram to hundred kiloton mass range. A recent variation on quark nugget dark matter models postulates that a first-order imbalance between matter and antimatter at the quark-gluon phase transition in the early universe could lead to most of the dark matter bound into heavy (baryon number $B \sim 10^{25}$) anti-quark nuggets in the current epoch, explaining both the dark matter preponderance and the matter-antimatter asymmetry. Interactions of these massive objects with normal matter in the Earth and Sun will lead to annihilation and an associated neutrino flux in the $\sim 20-50$ MeV range. We calculate these fluxes for anti-quark nuggets of sufficient number density to account for the dark matter and find that current neutrino flux limits from Super-Kamiokande exclude these objects as major dark matter candidates at a high confidence level over a very wide range of masses. Anti-quark nuggets in the previously allowed mass range cannot account for more than $\sim 1/5$ of the dark matter flux.
Multiple tracers of the cosmic density field, with different bias, number and luminosity evolution, can be used to measure the large-scale properties of the Universe. We show how an optimal combination of tracers can be used to detect general-relativistic effects in the observed density of sources. We forecast for the detectability of these effects, as well as measurements of primordial non-Gaussianity and large-scale lensing magnification with current and upcoming large-scale structure experiments. In particular we quantify the significance of these detections in the short term with experiments such as the Dark Energy Survey (DES), and in the long term with the Large Synoptic Survey Telescope (LSST) and the Square Kilometre Array (SKA). We review the main observational challenges that must be overcome to carry out these measurements.
Host star metallicities have been used to infer observational constraints on planet formation throughout the history of the exoplanet field. The giant planet metallicity correlation has now been widely accepted, but questions remain as to whether the metallicity correlation extends to the small terrestrial-sized planets. Here, we report metallicities for a sample of 518 stars in the Kepler field that have no detected transiting planets and compare their metallicity distribution to a sample of stars that hosts small planets (Rp < 1.7 R_Earth). Importantly, both samples have been analyzed in a homogeneous manner using the same set of tools (Stellar Parameters Classification tool; SPC). We find the average metallicity of the sample of stars without detected transiting planets to be [m/H]_SNTP,dwarf = -0.02 +- 0.02 dex and the sample of stars hosting small planets to be [m/H]_STP = -0.02 +- 0.02 dex. The average metallicities of the two samples are indistinguishable within the uncertainties, and the two-sample Kolmogorov-Smirnov test yields a p-value of 0.68 (0.41 sigma), indicating a failure to reject the null hypothesis that the two samples are drawn from the same parent population. We conclude that the homogeneous analysis of the data presented here support the hypothesis that stars hosting small planets have a metallicity similar to stars with no known transiting planets in the same area of the sky.
While supernova remnants (SNRs) are widely thought to be powerful cosmic-ray accelerators, indirect evidence comes from a small number of well-studied cases. Here we systematically determine the gamma-ray emission detected by the Fermi Large Area Telescope (LAT) from all known Galactic SNRs, disentangling them from the sea of cosmic-ray generated photons in the Galactic plane. Using LAT data we have characterized the 1-100 GeV emission in 279 regions containing SNRs, accounting for systematic uncertainties caused by source confusion and instrumental response. We have also developed a method to explore some systematic effects on SNR properties caused by the modeling of the interstellar emission (IEM). The IEM contributes substantially to gamma-ray emission in the regions where SNRs are located. To explore the systematics we consider different model construction methods, different model input parameters, and independently fit the model components to the gamma-ray data. We will describe this analysis method in detail. In the First Fermi LAT SNR Catalog there are 30 sources classified as SNRs, using spatial overlap with the radio position. For all the remaining regions we evaluated upper limits on SNRs' emission. In this work we will present a study of the aggregate characteristics of SNRs, such as comparisons between GeV and radio sizes as well as fluxes and spectral indexes and with TeV.
Context. We have conducted a 3D MHD simulation of the solar corona above an active region in full scale and high resolution, which shows coronal loops, and plasma flows within them, similar to observations. Aims. We want to find the connection between the photospheric energy input by field-line braiding with the coronal energy conversion by Ohmic dissipation of induced currents. Methods. To this end we compare the coronal energy input and dissipation within our simulation domain above different fields of view, e.g. for a small loops system in the active region (AR) core. We also choose an ensemble of field lines to compare, e.g., the magnetic energy input to the heating per particle along these field lines. Results. We find an enhanced Ohmic dissipation of currents in the corona above areas that also have enhanced upwards-directed Poynting flux. These regions coincide with the regions where hot coronal loops within the AR core are observed. The coronal density plays a role in estimating the coronal temperature due to the generated heat input. A minimum flux density of about 200 Gauss is needed in the photosphere to heat a field line to coronal temperatures of about 1 MK. Conclusions. This suggests that the field-line braiding mechanism provides the coronal energy input and that the Ohmic dissipation of induced currents dominates the coronal heating mechanism.
Light gauge bosons can lead to resonant interactions between high-energy astrophysical neutrinos and the cosmic neutrino background. We study this possibility in detail, considering the ability of IceCube to probe such scenarios. We find the most dramatic effects in models with a very light $Z'$ ($m_{Z'} \lesssim 10$ MeV), which can induce a significant absorption feature at $E_{\nu} \sim\,$5-10$\,{\rm TeV} \times (m_{Z'}/{\rm MeV})^2$. In the case of the inverted hierarchy and a small sum of neutrino masses, such a light $Z'$ can result in a broad and deep spectral feature at $\sim\,$0.1-10$\,{\rm PeV} \times (m_{Z'}/{\rm MeV})^2$. Current IceCube data already excludes this case for a $Z'$ lighter than a few MeV and couplings greater than $g\sim10^{-4}$. We emphasize that the ratio of neutrino flavors observed by IceCube can be used to further increase their sensitivity to $Z'$ models and to other exotic physics scenarios.
The observed value of the Higgs mass indicates the possibility that there is no supersymmetry below the Planck scale and that the Higgs can play the role of the inflaton. We examine the general structure of the saddle point inflation in string-inspired theory without supersymmetry. We point out that the string scale is fixed to be around the GUT scale $\sim10^{16}$GeV in order to realize successful inflation. We find that the inflaton can be naturally identified with the Higgs field.
50% of the heavy element abundances are produced via slow neutron capture reactions in different stellar scenarios. The underlying nucleosynthesis models need the input of neutron capture cross sections. One of the fundamental signatures for active nucleosynthesis in our galaxy is the observation of long-lived radioactive isotopes, such as $^{60}$Fe with a half-life of $2.60\times10^6$ yr. To reproduce this $\gamma$-activity in the universe, the nucleosynthesis of $^{60}$Fe has to be understood reliably. A $^{60}$Fe sample produced at the Paul-Scherrer-Institut was activated with thermal and epithermal neutrons at the research reactor at the Johannes Gutenberg-Universit\"at Mainz. The thermal neutron capture cross section has been measured for the first time to $\sigma_{\text{th}}=0.226 \ (^{+0.044}_{-0.049})$ b. An upper limit of $\sigma_{\text{RI}} < 0.50$ b could be determined for the resonance integral. An extrapolation towards the astrophysicaly interesting energy regime between $kT$=10 keV and 100 keV illustrates that the s-wave part of the direct capture component can be neglected.
Neutron capture cross sections of unstable isotopes are important for neutron-induced nucleosynthesis as well as for technological applications. A combination of a radioactive beam facility, an ion storage ring and a high flux reactor would allow a direct measurement of neutron induced reactions over a wide energy range on isotopes with half lives down to minutes. The idea is to measure neutron-induced reactions on radioactive ions in inverse kinematics. This means, the radioactive ions will pass through a neutron target. In order to efficiently use the rare nuclides as well as to enhance the luminosity, the exotic nuclides can be stored in an ion storage ring. The neutron target can be the core of a research reactor, where one of the central fuel elements is replaced by the evacuated beam pipe of the storage ring. Using particle detectors and Schottky spectroscopy, most of the important neutron-induced reactions, such as (n,$\gamma$), (n,p), (n,$\alpha$), (n,2n), or (n,f), could be investigated.
With the standard deviation for the logarithm of the re-scaled range $\langle |F(t+\tau)-F(t)|\rangle$ of simulated fractal Brownian motions $F(t)$ given in a previous paper \cite{q14}, the method of least squares is adopted to determine the slope, $S$, and intercept, $I$, of the log$(\langle |F(t+\tau)-F(t)|\rangle)$ vs $\rm{log}(\tau)$ plot to investigate the limitation of this procedure. It is found that the reduced $\chi^2$ of the fitting decreases with the increase of the Hurst index, $H$ (the expectation value of $S$), which may be attributed to the correlation among the re-scaled ranges. Similarly, it is found that the errors of the fitting parameters $S$ and $I$ are usually smaller than their corresponding standard deviations. These results show the limitation of using the simple least square method to determine the dimension of a fractal time series. Nevertheless, they may be used to reinterpret the fitting results of the least square method to determine the dimension of fractal Brownian motions more self-consistently. The currency exchange rate between Euro and Dollar is used as an example to demonstrate this procedure and a fractal dimension of 1.511 is obtained for spans greater than 30 transactions.
Effects of the completely unknown symmetry (isovector) potential of the \D on the total and differential \rpi in heavy-ion collisions at beam energies from 100 to 1000 MeV/A are explored within an isospin-dependent transport model IBUU. The effects are found to be negligible at beam energies above the pion production threshold due to the very short lifetimes of less than 2 fm/c for $\Delta$ resonances with masses around $m_{\Delta}=1232$ MeV, leaving the $\pi^-/\pi^+$ ratios of especially the energetic pions still a reliable probe of the high-density behavior of nuclear symmetry energy $E_{sym}(\rho)$. However, as the beam energy becomes deeply sub-threshold for pion production, effects of the $\Delta$ symmetry potential becomes appreciable especially on the \rpi of low-energy pions from the decays of low-mass $\Delta$ resonances which have lived long enough to be affected by their mean-field potentials, providing a useful tool to study the symmetry potential and spectroscopy of $\Delta$ resonances in neutron-rich nuclear matter. Interestingly though, even at the deeply sub-threshold beam energies, the differential \rpi of energetic pions remains sensitive to the \esym at supra saturation densities with little influence from the uncertain symmetry potential of the $\Delta$ resonance.
Links to: arXiv, form interface, find, astro-ph, recent, 1507, contact, help (Access key information)
At present, there is a great worldwide interest in the development of technologies that allow information about the X-ray emission from pulsating cosmic sources to be used to obtain navigation solutions for deep-space spacecraft. In this paper, we illustrate the technique for determining the spatial position of a spacecraft based on the already existing data from the RXTE X-ray space observatory. We show that the spacecraft position toward the Crab pulsar can be determined using an X-ray detector with an effective area of about 0.6 sq.m in the energy range 3-15 keV with an accuracy up to 730 m in a signal integration time of 1000 s. Extending the energy range to 1 keV (the efficiency of the RXTE/PCA spectrometer decreases dramatically at energies below 3 keV) will allow a spacecraft position accuracy of 400-450 m to be achieved at the same effective area and up to 300-350 m when using detectors with an effective area of ~1 sq.m in the energy range 1-10 keV.
How complex organic - and potentially prebiotic - molecules are formed in regions of low- and high-mass star-formation remains a central question in astrochemistry. In particular, with just a few sources studied in detail, it is unclear what role environment plays in complex molecule formation. In this light, a comparison of relative abundances of related species between sources might be useful to explain observed differences. We seek to measure the relative abundance between three important complex organic molecules, ethylene glycol ((CH$_2$OH)$_2$), glycolaldehyde (CH$_2$OHCHO) and methyl formate (HCOOCH$_3$), toward high-mass protostars and thereby provide additional constraints on their formation pathways. We use IRAM 30-m single dish observations of the three species toward two high-mass star-forming regions - W51/e2 and G34.3+0.2 - and report a tentative detection of (CH2OH)2 toward both sources. Assuming that (CH$_2$OH)$_2$, CH$_2$OHCHO and HCOOCH$_3$ spatially coexist, relative abundance ratios, HCOOCH$_3$/(CH$_2$OH)$_2$, of 31 and 35 are derived for G34.3+0.2 and W51/e2, respectively. CH$_2$OHCHO is not detected, but the data provide lower limits to the HCOOCH$_3$/CH$_2$OHCHO abundance ratios of $\ge$193 for G34.3+0.2 and $\ge$550 for W51/e2. A comparison of these results to measurements from various sources in the literature indicates that the source luminosities may be correlated with the HCOOCH$_3$/(CH$_2$OH)$_2$ and HCOOCH$_3$/CH$_2$OHCHO ratios. This apparent correlation may be a consequence of the relative timescales each source spend at different temperatures-ranges in their evolution. Furthermore, we obtain lower limits to the ratio of (CH$_2$OH)$_2$/CH2OHCHO for G34.3+0.2 ($\ge$6) and W51/e2 ($\ge$16). This result confirms that a high (CH$_2$OH)$_2$/CH$_2$OHCHO abundance ratio is not a specific property of comets, as previously speculated.
We present observations of Swift J1112.2-8238, and identify it as a candidate relativistic tidal disruption flare (rTDF). The outburst was first detected by Swift/BAT in June 2011 as an unknown, long-lived (order of days) $\gamma$-ray transient source. We show that its position is consistent with the nucleus of a faint galaxy for which we establish a likely redshift of $z=0.89$ based on a single emission line that we interpret as the blended [OII]$\lambda3727$ doublet. At this redshift, the peak X/$\gamma$-ray luminosity exceeded $10^{47}$ ergs s$^{-1}$, while a spatially coincident optical transient source had $i^{\prime} \sim 22$ (M$_g \sim -21.4$ at $z=0.89$) during early observations, $\sim 20$ days after the Swift trigger. These properties place Swift J1112.2-8238 in a very similar region of parameter space to the two previously identified members of this class, Swift J1644+57 and Swift J2058+0516. As with those events the high-energy emission shows evidence for variability over the first few days, while late time observations, almost 3 years post-outburst, demonstrate that it has now switched off. Swift J1112.2-8238 brings the total number of such events observed by Swift to three, interestingly all detected by Swift over a $\sim$3 month period ($<3\%$ of its total lifetime as of March 2015). While this suggests the possibility that further examples may be uncovered by detailed searches of the BAT archives, the lack of any prime candidates in the years since 2011 means these events are undoubtedly rare.
Precise measurements of the baryon acoustic oscillation (BAO) scale as a standard ruler in the clustering pattern of large-scale structure is a central goal of current and future galaxy surveys. The BAO peak may be sharpened using the technique of density-field reconstruction, in which the bulk displacements of galaxies are estimated using a Zel'dovitch approximation. We use numerical simulations to demonstrate how the accuracy of this approximation depends strongly on local environment, and how this information may be used to construct an improved BAO measurement through environmental re-weighting and using higher-order perturbation theory. We outline further applications of the displacement field for testing cosmological models.
The dispersion of the star-formation main sequence (SFMS) reflects the diversity of star formation histories and variation in star formation rates (SFRs) in star-forming galaxies (SFGs) with similar stellar masses ($M^\ast$). We examine the dispersion of local SFMS using a complete sample of Sloan Digital Sky Survey galaxies at 0.01$<z<$0.03 with $\log(M^\ast/M_\odot)>$8.8. The SFRs are estimated from H$\alpha$ in combination with 22$\mu m$ observation from WISE. The catalog of bulge+disk decomposition from Simard et al. (2011) is available for the sample galaxies. We measure the dispersion of specific SFR (SSFR) as a function of $M^*$. We confirm that the dispersion increases with $M^*$ from 0.37$\pm0.01$dex at $\log(M^\ast/M_\odot)<$9.6 to 0.51$\pm0.02$dex at $\log(M^\ast/M_\odot)>$10.2. Despite star formation is mostly associated with disks, the dispersion of disk SSFR still increases with $M^*$. We conclude that the presence of bulges/bars is likely responsible for the large dispersion of SSFR in massive SFGs while low-mass SFGs are mostly disk-dominated and thus with small dispersion. Our results suggest that star formation on galactic scales is dramatically affected by central dense structures through both enhancing and/or quenching processes; while lower-mass SFGs tend to have less bursty star formation histories. However, the dispersion of SSFR becomes significantly smaller and remains constant when only disk-dominated SFGs are counted. This finding implies that the impact of stochastic stellar feedback on star formation is likely to follow the same pattern in all disk galaxies, showing no correlation with halo potential.
We present the first optical (590--890 nm) imaging polarimetry observations of the pre-transitional protoplanetary disk around the young solar analog LkCa 15, addressing a number of open questions raised by previous studies. We detect the previously unseen far side of the disk gap, confirm the highly eccentric scattered-light gap shape that was postulated from near-infrared imaging, at odds with the symmetric gap inferred from millimeter interferometry. Furthermore, we resolve the inner disk for the first time and trace it out to 30 AU. This new source of scattered light may contribute to the near-infrared interferometric signal attributed to the protoplanet candidate LkCa 15 b, which lies embedded in the outer regions of the inner disk. Finally, we present a new model for the system architecture of LkCa 15 that ties these new findings together. These observations were taken during science verification of SPHERE ZIMPOL and demonstrate this facility's performance for faint guide stars under adverse observing conditions.
We present 50 nights of polarimetric observations of HD 189733 in $B$ band using the POLISH2 aperture-integrated polarimeter at the Lick Observatory Shane 3-m telescope. This instrument, commissioned in 2011, is designed to search for Rayleigh scattering from short-period exoplanets due to the polarized nature of scattered light. Since these planets are spatially unresolvable from their host stars, the relative contribution of the planet-to-total system polarization is expected to vary with an amplitude of order 10 parts per million (ppm) over the course of the orbit. Non-zero and also variable at the 10 ppm level, the inherent polarization of the Lick 3-m telescope limits the accuracy of our measurements and currently inhibits conclusive detection of scattered light from this exoplanet. However, the amplitude of observed variability conservatively sets a $3 \sigma$ upper limit to the planet-induced polarization of the system of 58 ppm in $B$ band, which is consistent with a previous upper limit from the POLISH instrument at the Palomar Observatory 5-m telescope (Wiktorowicz 2009). A physically-motivated Rayleigh scattering model, which includes the depolarizing effects of multiple scattering, is used to conservatively set a $3 \sigma$ upper limit to the geometric albedo of HD 189733b of $A_g < 0.37$. This value is consistent with the value $A_g = 0.226 \pm 0.091$ derived from occultation observations with HST STIS (Evans et al. 2013), but it is inconsistent with the large $A_g = 0.61 \pm 0.12$ albedo reported by (Berdyugina et al. 2011).
We use a non-detection in $\nu = 1.4\,$GHz Green Bank Telescope observations of the ultra-faint dwarf spheroidal galaxy Segue I, which could be immersed in a non-negligible halo magnetic field of the Milky Way, to place bounds on particle dark matter properties. We model the galaxy using an Einasto dark matter profile, and compute the expected synchrotron flux from dark matter annihilation as a function of the magnetic field strength $B$, diffusion coefficient $D_0$, and particle mass $m_\chi$ for different annihilation channels. The data strongly disfavor annihilations to $e^+e^-$ for $m_\chi \lesssim 50\,$GeV, but are not sensitive to the $b \bar b$ channel. Adopting a fiducial $B \sim 2\,\mu$G inferred from Segue I's proximity to the Milky Way, our models of annihilation to $\tau^+\tau^-$ with $m_\chi = 30\,$GeV require an intermediate value of $D_0$ for consistency with the data. The most compelling limits are obtained for WIMP annihilation to $\mu^+\mu^-$: we exclude $m_\chi \lesssim 30\,$GeV$\,\rightarrow\mu^+\mu^-$ at 95% confidence, unless $D_0$ exceeds the Milky Way value or $B$ is significantly smaller than we have assumed.
We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter density slope, $\alpha$, on the ratio between stellar-to-halo mass. We show that this relation holds approximately at all redshifts (with an intrinsic scatter of ~0.18 in $\alpha$). This implies that in practically all haloes the shape of their inner density profile changes quite substantially over cosmic time, as they grow in stellar and total mass. Thus, depending on their final stellar-to-halo mass ratio, haloes can either form and keep a substantial density core (size~1 kpc), or form and then destroy the core and re-contract the halo, going back to a cuspy profile, which is even steeper than CDM predictions for massive galaxies (~1e12 Msun). We show that results from the NIHAO suite are in good agreement with recent observational measurements of $\alpha$ in dwarf galaxies. Overall our results suggest that the notion of a universal density profile for dark matter haloes is no longer valid in the presence of galaxy formation.
The number of magnetic stars detected among massive stars is small; nevertheless, the role played by the magnetic field in stellar evolution cannot be disregarded. Links between line profile variability, enhancements/depletions of surface chemical abundances, and magnetic fields have been identified for low-mass B-stars, but for the O-type domain this is almost unexplored. Based on FORS2 and HARPS spectropolarimetric data, we present the first detection of a magnetic field in HD54879, a single slowly rotating O9.7 V star. Using two independent and different techniques we obtained the firm detection of a surface average longitudinal magnetic field with a maximum amplitude of about 600 G, in modulus. A quantitative spectroscopic analysis of the star with the stellar atmosphere code FASTWIND results in an effective temperature and a surface gravity of 33000$\pm1000$ K and 4.0$\pm0.1$ dex. The abundances of carbon, nitrogen, oxygen, silicon, and magnesium are found to be slightly lower than solar, but compatible within the errors. We investigate line-profile variability in HD54879 by complementing our spectra with spectroscopic data from other recent OB-star surveys. The photospheric lines remain constant in shape between 2009 and 2014, although H$\alpha$ shows a variable emission. The H$\alpha$ emission is too strong for a standard O9.7 V and is probably linked to the magnetic field and the presence of circumstellar material. Its normal chemical composition and the absence of photospheric line profile variations make HD54879 the most strongly magnetic, non-variable single O-star detected to date.
Dynamical scattering of binaries and triple systems of stars, planets, and compact objects may produce highly inclined triple systems subject to Kozai-Lidov (KL) oscillations, potentially leading to collisions, mergers, Type Ia supernovae, and other phenomena. We present the results of more than 400 million gravitational scattering experiments of binary-binary, triple-single, and triple-binary scattering. We compute the cross sections for all possible outcomes and explore their dependencies on incoming velocity, mass, semi-major axis, and eccentricity, including analytic fits and discussion of the velocity dependence. For the production of new triple systems by scattering we find that compact triples are preferred, with ratios of outer to inner semi-major axes of ~few--100, flat or quasi-thermal eccentricity distributions, and flat distributions in cosine of the mutual inclination. Dynamically formed triples are thus subject to strong KL oscillations, the "eccentric Kozai mechanism," and non-secular effects. For single and binary flyby encounters with triple systems, we compute the cumulative cross section for changes to the mutual inclination, eccentricity, and semi-major axis ratio. We apply these results to scattering events in the field, open clusters, and globular clusters, and explore the implications for Type Ia supernovae via collisions and mergers, stellar collisions, and the lifetime and dynamical isolation of triple systems undergoing KL oscillations. An Appendix provides an analysis of the velocity dependence of the collision cross section in binary-single scattering.
This paper uses the equations of motion that govern the dynamics of galaxies in the Local Volume to measure the (total) masses of the Large Magellanic Cloud (LMC), Milky Way (MW) and Andromeda (M31) galaxies. We implement a Bayesian technique devised by Pe\~narrubia et al. (2014) to simultaneously fit point-mass orbits to published distances and velocities of galaxies within 3~Mpc as well as to the relative position and velocity vectors of the Milky Way-Andromeda pair derived from HST observations. Our analysis returns a LMC mass $M_{LMC}=0.25_{-0.08}^{+0.09}\times 10^{12}M_\odot$ at a 68% confidence level. The masses of the Milky Way, $M_{MW}=1.04_{-0.23}^{+0.26}\times 10^{12}M_\odot$, and Andromeda, $M_{M31}=1.33_{-0.33}^{+0.39}\times 10^{12}M_\odot$, are consistent with previous estimates that neglect the impact of the LMC on the observed Hubble flow. Such a large LMC mass is indicative of an extended dark matter halo and supports the scenario where this galaxy is just past its first pericentric approach. Consequently, our results suggest that the LMC may induce significant perturbations on the Galactic potential.
Using 3D global hydro simulations coupled with radiative transfer calculations, we study the appearance of density waves induced by giant planets in direct imaging observations at near infrared wavelengths. We find that a 6 MJ planet in a typical disk around a 1 M_sun star can produce prominent and detectable spiral arms both interior and exterior to its orbit. The inner arms have (1) two well separated arms in roughly m=2 symmetry, (2) exhibit ~10-15 degrees pitch angles, (3) ~180-270 degrees extension in the azimuthal direction, and (4) ~150% surface brightness enhancement, all broadly consistent with observed spiral arms in the SAO 206462 and MWC 758 systems. The outer arms cannot explain observations as they are too tightly wound given typical disk scale height. We confirm previous results that the outer density waves excited by a 1 MJ planet exhibit low contrast in the IR and are practically not detectable. We also find that 3D effects of the waves are important. Compared to isothermal models, density waves in adiabatic disks exhibit weaker contrast in surface density but stronger contrast in scattered light images, due to a more pronounced vertical structure in the former caused by shock heating. To drive observed pairs of arms with an external companion on a circular orbit, a massive planet, possibly a brown dwarf, is needed at around [r~0.7", PA~10 degrees] (position angle PA from north to east) in SAO 206462 and [r~0.6 , PA~10 degrees] in MWC 758. Their existence may be confirmed by direct imaging planet searches.
Recent direct imaging observations have revealed spiral structure in protoplanetary disks. Previous studies have suggested that planet-induced spiral arms cannot explain these spiral patterns, as 1) the pitch angle of the spiral arm is larger in observations than that predicted by the linear density wave theory, 2) the contrast of the spiral arm is higher in observations than in synthetic observations based on two dimensional planet-disk simulations. We have carried out three dimensional (3-D) hydrodynamical simulations to study spiral wakes/shocks excited by young planets. We find that, in contrast with linear theory, the pitch angle of spiral arms does depend on the planet mass, which can be explained by the non-linear density wave theory. The more massive is the planet, the larger pitch angle the spiral arm has. A secondary spiral arm, especially for the inner arms, is also excited by the planet. The more massive is the planet, the larger is the separation in the azimuthal direction between the primary and secondary arms. We also find that although the arms in the outer disk do not exhibit much vertical motion, the inner arms have significant vertical motion, which boosts the density perturbation at the disk atmosphere by more than a factor of 10 compared with that at the disk midplane. Combining hydrodynamical models with Monte-Carlo radiative transfer calculations, we find that the inner spiral arms are considerably more prominent in synthetic near-IR images using full 3-D hydrodynamical models than images based on 2-D models, indicating the need to model observations with full 3-D hydrodynamics. Overall, spiral arms (especially inner arms) excited by planetary-mass objects are prominent features that are observable by current near-IR imaging facilities, and the shape of the spiral arms informs us not only about the position but also about the mass of the companion.
A gauge singlet scalar with non-minimal coupling to gravity can drive inflation and later freeze out to become cold dark matter. We explore this idea by revisiting inflation in the singlet direction (S-inflation) and Higgs Portal Dark Matter in light of the Higgs discovery, limits from LUX and observations by Planck. We show that large regions of parameter space remain viable, so that successful inflation is possible and the dark matter relic abundance can be reproduced. Moreover, the scalar singlet can stabilise the electroweak vacuum and at the same time overcome the problem of unitarity-violation during inflation encountered by Higgs Inflation, provided the singlet is a real scalar. The 2-$\sigma$ Planck upper bound on $n_s$ imposes that the singlet mass is below 2 TeV, so that almost the entire allowed parameter range can be probed by XENON1T.
To measure the stellar and orbital properties of the metal-poor RS CVn binary o Draconis (o Dra), we directly detect the companion using interferometric observations obtained with the Michigan InfraRed Combiner at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array. The H-band flux ratio between the primary and secondary stars is the highest confirmed flux ratio (370 +/- 40) observed with long-baseline optical interferometry. These detections are combined with radial velocity data of both the primary and secondary stars, including new data obtained with the Tillinghast Reflector Echelle Spectrograph on the Tillinghast Reflector at the Fred Lawrence Whipple Observatory and the 2-m Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We determine an orbit from which we find model-independent masses and ages of the components (M_A = 1.35 +\- 0.05 M_Sun, M_B = 0.99 +\- 0.02 M_Sun, system age = 3.0 -\+ 0.5 Gyr). An average of a 23-year light curve of o Dra from the Tennessee State University Automated Photometric Telescope folded over the orbital period newly reveals eclipses and the quasi-sinusoidal signature of ellipsoidal variations. The modeled light curve for our system's stellar and orbital parameters confirm these ellipsoidal variations due to the primary star partially filling its Roche lobe potential, suggesting most of the photometric variations are not due to stellar activity (starspots). Measuring gravity darkening from the average light curve gives a best-fit of beta = 0.07 +\- 0.03, a value consistent with conventional theory for convective envelope stars. The primary star also exhibits an anomalously short rotation period, which, when taken with other system parameters, suggests the star likely engulfed a low-mass companion that had recently spun-up the star.
We present a chemical abundance distribution study in 14 $\alpha$, odd-Z, even-Z, light, and Fe-peak elements of approximately 3200 intermediate metallicity giant stars from the APOGEE survey. The main aim of our analysis is to explore the Galactic disk-halo transition region between -1.20 $<$ [Fe/H] $<$ -0.55 as a means to study chemical difference (and similarities) between these components. In this paper, we show that there is an $\alpha$-poor and $\alpha$-rich sequence within both the metal-poor and intermediate metallicity regions. Using the Galactic rest-frame radial velocity and spatial positions, we further separate our sample into the canonical Galactic components. We then studied the abundances ratios, of Mg, Ti, Si, Ca, O, S, Al, C+N, Na, Ni, Mn, V, and K for each of the components and found the following: (1) the $\alpha$-poor halo subgroup is chemically distinct in the $\alpha$-elements (particularly O, Mg, and S), Al, C+N, and Ni from the $\alpha$-rich halo, consistent with the literature confirming the existence of an $\alpha$-poor accreted halo population; (2) the canonical thick disk and halo are not chemically distinct in all elements indicating a smooth transition between the thick disk and halo; (3) a subsample of the $\alpha$-poor stars at metallicities as low as [Fe/H] $\sim$ -0.85 dex are chemically and dynamically consistent with the thin disk indicating that the thin disk may extend to lower metallicities than previously thought, and (4) that the location of the most metal-poor thin disk stars are consistent with a negative radial metallicity gradient. Finally, we used our analysis to suggest a new set of chemical abundance planes ([$\alpha$/Fe], [C+N/Fe], [Al/Fe], and [Mg/Mn]) that may be able to chemically label the Galactic components in a clean and efficient way independent of kinematics.
In this work, I explore an empirically motivated model for investigating the relationship between galaxy stellar masses, star formation rates and their halo masses and mass accretion histories. The core statistical quantity in this model is the stellar mass assembly distribution, $P(dM_{*}/dt|\mathbf{X},a)$, which specifies the probability density distribution of stellar mass assembly rates given a set of halo properties $\mathbf{X}$ and epoch $a$. Predictions from this model are obtained by integrating the stellar mass assembly distribution (SMAD) over halo merger trees, easily obtained from modern, high-resolution $N$-body simulations. Further properties of the galaxies hosted by the halos can be obtained by post-processing the stellar mass assembly histories with stellar population synthesis models. In my particular example implementation of this model, I use the \citet{behroozi13a} constraint on the median stellar mass assembly rates of halos as a function of their mass and redshift to construct an example parameterization of $P(dM_{*}/dt|\mathbf{X},a)$. This SMAD is then integrated over individual halo mass accretion histories from $N$-body merger trees starting at z = 4, using simple rules to account for merging halos. I find that this a simple model can reproduce qualitatively the bimodal features of the low-redshift galaxy population, including the qualitative split in the two-point clustering as a function of specific star formation rate. These results indicate that models which directly couple halo and galaxy growth through simple efficiency functions can naturally predict the star formation rate bimodality in higher-order statistics of the galaxy field, such as its two-point correlations or galactic conformity signals.
We present Sedonu, a new open source, steady-state, special relativistic Monte Carlo (MC) neutrino transport code, available at bitbucket.org/srichers/sedonu. The code calculates the energy- and angle-dependent neutrino distribution function on fluid backgrounds of any number of spatial dimensions, calculates the rates of change of fluid internal energy and electron fraction, and solves for the equilibrium fluid temperature and electron fraction. We apply this method to snapshots from two dimensional simulations of accretion disks left behind by binary neutron star mergers, varying the input physics and comparing to the results obtained with a leakage scheme for the case of a central black hole and a central hypermassive neutron star. Neutrinos are guided away from the densest regions of the disk and escape preferentially around 45 degrees from the equatorial plane. Neutrino heating is strengthened by MC transport a few scale heights above the disk midplane near the innermost stable circular orbit, potentially leading to a stronger neutrino-driven wind. Neutrino cooling in the dense midplane of the disk is stronger when using MC transport, leading to a globally higher cooling rate by a factor of a few and a larger leptonization rate by an order of magnitude. We calculate neutrino pair annihilation rates and estimate that an energy of 2.8e46 erg is deposited within 45 degrees of the symmetry axis over 300 ms when a central BH is present. Similarly, 1.9e48 erg is deposited over 3 s when an HMNS sits at the center, but neither estimate is likely to be sufficient to drive a GRB jet.
Interpretation of the oscillation spectrum of the slowly pulsating B-type star HD21071 is presented. We show that non-rotating models cannot account for the two highest amplitude frequencies and taking into account the effects of rotation is necessary. Rotating seismic models are constructed using various chemical compositions, opacity data, core overshooting parameters and rotational velocities. There are prospects for seismic modelling of SPB stars, even if no asymptotic pattern is observed in their oscillation spectra, provided an unambiguous mode identification is doable and the effects of rotation are properly included.
A growing body of evidence suggests that multiple dynamo mechanisms can drive magnetic variability on different timescales, not only in the Sun but also in other stars. Many solar activity proxies exhibit a quasi-biennial ($\sim$2 year) variation, which is superimposed upon the dominant 11 year cycle. A well-characterized stellar sample suggests at least two different relationships between rotation period and cycle period, with some stars exhibiting long and short cycles simultaneously. Within this sample, the solar cycle periods are typical of a more rapidly rotating star, implying that the Sun might be in a transitional state or that it has an unusual evolutionary history. In this work, we present new and archival observations of dual magnetic cycles in the young solar analog HD 30495, an $\sim$1 Gyr-old G1.5V star with a rotation period near 11 days. This star falls squarely on the relationships established by the broader stellar sample, with short-period variations at $\sim$1.7 years and a long cycle of $\sim$12 years. We measure three individual long-period cycles and find durations ranging from 9.6-15.5 years. We find the short-term variability to be intermittent, but present throughout the majority of the time series, though its occurrence and amplitude are uncorrelated with the longer cycle. These essentially solar-like variations occur in a Sun-like star with more rapid rotation, though surface differential rotation measurements leave open the possibility of a solar equivalence.
The variation of the solar diameter in time and in position angle has implications in astrophysics and in general relativity, as the long series of studies attest. The Transits of Venus in 2004 and 2012 have been carefully studied because of the rarity of the phenomenon and its historical importance due the AU measure and to the discovery of Venus atmosphere. The characterization of Venus atmosphere and the measure of the solar diameter to the milliarcsecond level of precision have been studied also from satellite images. The results of the solar diameter measurements made with the observations in Athens (2004) and at the Huairou Solar Observing Station in China (2012) are presented. The topic of the oblateness of the Sun at sunset and its intrinsic value is drafted to introduce the general public to the relativistic relevance of measuring the solar figure, in the occasion of the International Year of Light 2015.
Galactic cosmic ray (CRs) sources, classically proposed to be Supernova Remnants (SNRs), must meet the energetic particle content required by direct measurements of high energy CRs. Indirect gamma-ray measurements of SNRs with the Fermi Large Area Telescope (LAT) have now shown directly that at least three SNRs accelerate protons. With the first Fermi LAT SNR Catalog, we have systematically characterized the GeV gamma-rays emitted by 279 SNRs known primarily from radio surveys. We present these sources in a multiwavelength context, including studies of correlations between GeV and radio size, flux, and index, TeV index, and age and environment tracers, in order to better understand effects of evolution and environment on the GeV emission. We show that previously sufficient models of SNRs' GeV emission no longer adequately describe the data. To address the question of CR origins, we also examine the SNRs' maximal CR contribution assuming the GeV emission arises solely from proton interactions. Improved breadth and quality of multiwavelength data, including distances and local densities, and more, higher resolution gamma-ray data with correspondingly improved Galactic diffuse models will strengthen this constraint.
Measuring the solar diameter at all position angles gives the complete figure of the Sun. Their asphericities have implications in classical physics and general relativity, and the behavior of the optical systems used in the direct measurements is to be known accurately. A solar filter is a plane-parallel glass with given absorption, and here we study the departures from the parallelism of the faces of a crystal slab 5 mm thick, because of static deformations. These deformations are rescaled to the filter's dimensions. Related to the Solar Disk Sextant experiment and to the Reflecting Heliometer of Rio de Janeiro a simplified model of the influences of the inclination between the external and the internal surfaces of a glass solar filter, is discussed.
We use a large sample of isolated dark matter halo pairs drawn from cosmological N-body simulations to identify candidate systems whose kinematics match that of the Local Group of Galaxies (LG). We find, in agreement with the "timing argument" and earlier work, that the separation and approach velocity of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the pair of ~ 5*10^12 M_sun. A mass this large, however, is difficult to reconcile with the small relative tangential velocity of the pair, as well as with the small deceleration from the Hubble flow observed for the most distant LG members. Halo pairs that match these three criteria have average masses a factor of ~2 times smaller than suggested by the timing argument, but with large dispersion, spanning more than a decade in mass. Guided by these results, we have selected 12 halo pairs with total mass in the range 1.6-3.6 *10^12 M_sun for the APOSTLE project (A Project Of Simulations of The Local Environment), a suite of resimulations at various numerical resolution levels (reaching up to ~10^4 M_sun per gas particle) that use the hydrodynamical code and subgrid physics developed for the EAGLE project. These simulations reproduce, by construction, the main kinematics of the MW-M31 pair, and produce satellite populations whose overall number, luminosities, and kinematics are in good agreement with observations of the MW and M31 companions. These diagnostics are sensitive to the total mass assumed for the MW-M31 pair; indeed, the LG satellite population would be quite difficult to reproduce for pair masses as high as indicated by the timing argument. The APOSTLE candidate systems thus provide an excellent testbed to confront directly many of the predictions of the Lambda-CDM cosmology with observations of our local Universe.
The light curves of some luminous supernovae are suspected to be powered by the spindown energy of a rapidly rotating magnetar. Here we describe a possible signature of the central engine: a burst of shock breakout emission occurring several days after the supernova explosion. The energy input from the magnetar inflates a high-pressure bubble that drives a shock through the pre-exploded supernova ejecta. If the magnetar is powerful enough, that shock will near the ejecta surface and become radiative. At the time of shock breakout, the ejecta will have expanded to a large radius (~10^{14} cm) so that the radiation released is at optical/ultraviolet wavelengths (T ~ 20,000 K) and lasts for several days. The luminosity and timescale of this magnetar driven shock breakout are similar to the first peak observed recently in the double-peaked light curve of SN-LSQ14BDQ. However, for a large region of model parameter space, the breakout emission is predicted to be dimmer than the diffusive luminosity from direct magnetar heating. A distinct double peaked light curve may only be conspicuous in events with a slowly rising diffusive light curve, due perhaps to a large ejecta mass and/or inefficient thermalization of the magnetar spindown energy at early times. Otherwise, the breakout may be detectable as a kink in the early luminosity or color evolution, or as a small but abrupt rise in the photospheric velocity. A similar breakout signature may accompany other central engines in supernovae, such as a black hole accreting fallback material.
Despite almost all being acquired as photons, astronomical data from different instruments and at different stages in its life may exist in different formats to serve different purposes. Beyond the data itself, descriptive information is associated with it as metadata, either included in the data format or in a larger multi-format data structure. Those formats may be used for the acquisition, processing, exchange, and archiving of data. It has been useful to use similar formats, or even a single standard to ease interaction with data in its various stages using familiar tools. Knowledge of the evolution and advantages of present standards is useful before we discuss the future of how astronomical data is formatted. The evolution of the use of world coordinates in FITS is presented as an example.
We study rotating thermal convection in spherical shells as prototype for flow in the cores of terrestrial planets, gas planets or in stars. We base our analysis on a set of about 450 direct numerical simulations of the (magneto)hydrodynamic equations under the Boussinesq approximation. The Ekman number ranges from $10^{-3}$ to $10^{-6}$. Four sets of simulations are considered: non-magnetic simulations and dynamo simulations with either free-slip or no-slip flow boundary conditions. The non-magnetic setup with free-slip boundaries generates the strongest zonal flows. Both non-magnetic simulations with no-slip flow boundary conditions and self-consistent dynamos with free-slip boundaries have drastically reduced zonal-flows. Suppression of shear leads to a substantial gain in heat-transfer efficiency, increasing by a factor of 3 in some cases. Such efficiency enhancement occurs as long as the convection is significantly influenced by rotation. At higher convective driving the heat-transfer efficiency trends towards that of the classical non-rotating Rayleigh-B\'enard system. Analysis of the latitudinal distribution of heat flow reveals that the shear is most effective at suppressing heat-transfer at low-latitudes. We also explore the influence of the magnetic field on the {\em non-zonal} flow components to test the idea of `magnetostrophic' convection. Comparing non-magnetic with dynamo simulations for no-slip boundaries, which reduce the influence of zonal flow, we find that at an Ekman number of $10^{-5}$ the flow in polar regions is significantly affected by the presence of the magnetic field. Near the equator, changes remain subtle. In a few simulations at an Ekman number of $10^{-6}$ we see dramatic changes in flow structure at {\em all} latitudes due to the presence of the magnetic field, hinting that we are approaching a magnetostrophic regime at such low Ekman numbers.
We report the discovery of one extremely metal-poor (EMP; [Fe/H]<-3) and one ultra metal-poor (UMP; [Fe/H]<-4) star selected from the SDSS/SEGUE survey. These stars were identified as EMP candidates based on their medium-resolution (R~2,000) spectra, and were followed-up with high-resolution (R~35,000) spectroscopy with the Magellan-Clay Telescope. Their derived chemical abundances exhibit good agreement with those of stars with similar metallicities. We also provide new insights on the formation of the UMP stars, based on comparison with a new set of theoretical models of supernovae nucleosynthesis. The models were matched with 20 UMP stars found in the literature, together with one of the program stars (SDSS J1204+1201), with [Fe/H]=-4.34. From fitting their abundances, we find that the supernovae progenitors, for stars where carbon and nitrogen are measured, had masses ranging from 20.5 M_sun to 28 M_sun and explosion energies from 0.3 to 0.9x10^51 erg. These results are highly sensitive to the carbon and nitrogen abundance determinations, which is one of the main drivers for future high-resolution follow-up of UMP candidates. In addition, we are able to reproduce the different CNO abundance patterns found in UMP stars with a single progenitor type, by varying its mass and explosion energy.
We use a new, improved version of the HI Parkes All-Sky Survey to search for HI emission from nine new, ultra-faint Milky Way satellite galaxy candidates recently discovered in data from the Dark Energy Survey. None of the candidates is detected in HI, implying upper limits for their HI masses of typically several hundred to a few thousand solar masses. The resulting upper limits on M_HI / L_V and M_HI / M_star suggest that at least some of the new galaxy candidates are HI deficient. This finding is consistent with the general HI deficiency of satellite galaxies located within the Milky Way's virial radius and supports the hypothesis that gas is being removed from satellites by tidal and ram-pressure forces during perigalactic passages. In addition, some of the objects may be embedded in, and interacting with, the extended neutral and ionised gas filaments of the Magellanic Stream.
In this review we describe recent observational and theoretical developments in our understanding of pulsar winds and pulsar-wind nebulae (PWNe). We put special emphasis on the results from observations of well-characterized PWNe of various types (e.g., torus-jet and bowshock-tail), the most recent MHD modeling efforts, and the status of the flaring Crab PWN puzzle.
The integrated Sachs-Wolfe (ISW) effect provides us the information of the time evolution of gravitational potential. The cross-correlation between the cosmic microwave background (CMB) and the large scale structure (LSS) is known as a promising way to extract the ISW effect. Compared to CMB, the matter fluctuation can grow non-linearly and this is represented in the gravitational potential. Compared to the linear ISW effect, this non-linear ISW effect known as the Rees-Sciama (RS) effect shows the unique behavior by changing the anti-correlated cross correlation between the CMB and the mass tracer into the positively correlated cross correlation. We show that the dependence of this flipping scale on dark energy models and it might be used as a new method to investigate dark energy models.
Measurements of pulsar flux densities are of great importance for understanding the pulsar emission mechanism and for predictions of pulsar survey yields and the pulsar population at large. Typically these flux densities are determined from phase-averaged "pulse profiles", but this method has limited applicability at low frequencies because the observed pulses can easily be spread out by interstellar effects like scattering or dispersion, leading to a non-pulsed continuum component that is necessarily ignored in this type of analysis. In particular for the class of the millisecond pulsars (MSPs) at frequencies below 200MHz, such interstellar effects can seriously compromise de- tectability and measured flux densities. In this paper we investigate MSP spectra based on a complementary approach, namely through investigation of archival con- tinuum imaging data. Even though these images lose sensitivity to pulsars since the on-pulse emission is averaged with off-pulse noise, they are insensitive to effects from scattering and provide a reliable way to determine the flux density and spectral indices of MSPs based on both pulsed and unpulsed components. Using the 74MHz VLSSr as well as the 325MHz WENSS and 1.4GHz NVSS catalogues, we investigate the imaging flux densities of MSPs and evaluate the likelihood of spectral turn-overs in this population. We determine three new MSP spectral indices and identify six new MSPs with likely spectral turn-overs.
New insights into the formation of interstellar formamide, a species of great relevance in prebiotic chemistry, are provided by electronic structure and kinetic calculations for the reaction NH2 + H2CO -> NH2CHO + H. Contrarily to what previously suggested, this reaction is essentially barrierless and can, therefore, occur under the low temperature conditions of interstellar objects thus providing a facile formation route of formamide. The rate coefficient parameters for the reaction channel leading to NH2CHO + H have been calculated to be A = 2.6x10^{-12} cm^3 s^{-1}, beta = -2.1 and gamma = 26.9 K in the range of temperatures 10-300 K. Including these new kinetic data in a refined astrochemical model, we show that the proposed mechanism can well reproduce the abundances of formamide observed in two very different interstellar objects: the cold envelope of the Sun-like protostar IRAS16293-2422 and the molecular shock L1551-B2. Therefore, the major conclusion of this Letter is that there is no need to invoke grain-surface chemistry to explain the presence of formamide provided that its precursors, NH2 and H2CO, are available in the gas-phase.
Filamentary structures are common in molecular clouds. Explaining how they fragment to dense cores is a missing step in understanding their role in star formation. We perform a case study of whether low-mass filaments are close-to hydrostatic prior to their fragmentation, and whether their fragmentation agrees with gravitational fragmentation models. For this, we study the 6.5 pc long Musca molecular cloud that is an ideal candidate for a filament at an early stage of fragmentation. We employ dust extinction mapping in conjunction with near-infrared data from the NEWFIRM instrument, and 870 um dust continuum emission data from the LABOCA instrument, to estimate column densities. We use the data to identify fragments from the cloud and to determine the radial density distribution of its filamentary part. We compare the cloud's morphology with 13CO and C18O line emission observed with the APEX/SHeFI instrument. The Musca cloud is pronouncedly fragmented at its ends, but harbours a remarkably well-defined, 1.6 pc long filament in its Center region. The line mass of the filament is 21-31 Ms pc^-1 and FWHM 0.07 pc. Its radial profile can be fitted with a Plummer profile that has the power-index of 2.6 \pm 11%, flatter than that of an infinite hydrostatic filament. The profile can also be fitted with a hydrostatic cylinder truncated by external pressure. These models imply a central density of 5-10 x 10^4 cm^-3. The fragments in the cloud have a mean separation of 0.4 pc, in agreement with gravitational fragmentation. These properties, together with the subsonic and velocity-coherent nature of the cloud, suggest a scenario in which an initially hydrostatic cloud is currently gravitationally fragmenting. The fragmentation has started a few tenths of a Myr ago from the cloud ends, leaving its center yet relatively non-fragmented, possibly because of gravitational focusing in a finite geometry.
A nonlinear kinetic theory, combining cosmic ray (CR) acceleration in supernova remnants (SNRs) with their gas dynamics, is used to re-examine the nonthermal properties of the remnant of SN 1987A for an extended evolutionary period of 5-50 yr. This spherically symmetric model is approximately applied to the different features of the SNR which consist of (i) a blue supergiant wind and bubble, and (ii) of the swept-up red supergiant (RSG) wind structures in the form of an HII region, an equatorial ring (ER), and an hourglass region. The RSG wind involves a mass loss rate that decreases significantly with elevation above and below the equatorial plane. The model adapts recent three-dimensional hydrodynamical simulations by Potter et al. (2014) which use a significantly smaller ionized mass of the ER than assumed in the earlier studies by the present authors. The SNR shock has recently swept up the ER which is the densest region in the immediate circumstellar environment. Therefore the expected gamma-ray energy flux density at TeV-energies at the current epoch has already reached its maximal value $\sim 10^{-13}$ erg cm$^{-2}$s$^{-1}$. This flux should decrease by a factor of about two over the next ten years.
We use data at 131, 171, and 304 A from the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO) to search for hot flux ropes in 141 M-class and X-class solar flares that occurred at solar longitudes equal to or larger than 50 degrees. Half of the flares were associated with coronal mass ejections (CMEs). The goal of our survey is to assess the frequency of hot flux ropes in large flares irrespective of their formation time relative to the onset of eruptions. The flux ropes were identified in 131 A images using morphological criteria and their high temperatures were confirmed by their absence in the cooler 171 and 304 A passbands. We found hot flux ropes in 45 of our events (32% of the flares); 11 of them were associated with confined flares while the remaining 34 were associated with eruptive flares. Therefore almost half (49%) of the eruptive events involved a hot flux rope configuration. The use of supplementary Hinode X-Ray Telescope (XRT) data indicates that these percentages should be considered as lower limits of the actual rates of occurrence of hot flux ropes in large flares.
We have observed the prompt emission of GRB100418A, from its beginning by the MAXI/SSC (0.7-7 keV) on board the International Space Station followed by the Swift/XRT (0.3-10 keV) observation. The light curve can be fitted by a combination of a power law component and an exponential component (decay constant is $31.6\pm 1.6$). The X-ray spectrum is well expressed by the Band function with $E_{\rm p}\leq$8.3 keV. This is the brightest GRB showing a very low value of $E_{\rm p}$. It is also consistent with the Yonetoku-relation ($E_{\rm p}$-$L_{\rm p}$) while it is not clear with the Amati-relation ($E_{\rm p}$-$E_{\rm iso}$).
We study selected properties of Solar Energetic Particle (SEP) events as inferred from their associated radio emissions. We used a catalogue of 115 SEP events that consists of entries of proton intensity enhancements at one AU, with complete coverage over solar cycle 23, based on high-energy (~68 MeV) protons from SOHO/ERNE and we calculated the proton release time at the Sun using velocity dispersion analysis (VDA). After an initial rejection of cases with unrealistic VDA path lengths, we assembled composite radio spectra for the remaining events using data from ground-based and space-borne radio-spectrographs. For every event we registered the associated radio emissions and we divided the events in groups according to their associated radio emissions. The proton release was found to be most often accompanied by both type III and II radio bursts, but a good association percentage was also registered in cases accompanied by type IIIs only. The worst association was found for the cases with type II only association. These radio association percentages support the idea that both flare- and shock-resident particle release processes are observed in high-energy proton events. In cases of type III-associated events we extended our study to the timings between the type III radio emission, the proton release, and the electron release as inferred from VDA based on Wind/3DP 20-646 keV data. Typically, the protons are released after the start of the associated type III bursts and simultaneously or before the release of energetic electrons. For the cases with type II radio association we found that the distribution of the proton release heights had a maximum at ~2.5 Rs. Most (69%) of the flares associated to our SEP events were located at the western hemisphere, with a peak within the well-connected region of 50-60 deg western longitude.
We discuss the properties of H$\alpha$ emission stars across the sample of 22035 spectra from the Gaia-ESO Survey internal data release, observed with the GIRAFFE instrument and largely belonging to stars in young open clusters. Automated fits using two independent Gaussian profiles and a third component that accounts for the nebular emission allow us to discern distinct morphological types of H$\alpha$ line profiles with the introduction of a simplified classification scheme. All in all we find 3765 stars with intrinsic emission and sort their spectra into eight distinct morphological categories: single--component emission, emission blend, sharp emission peaks, double emission, P-Cygni, inverted P-Cygni, self--absorption, and emission in absorption. We have more than one observation for 1430 stars in our sample, thus allowing a quantitative discussion of the degree of variability of H$\alpha$ emission profiles, which is expected for young, active objects. We present a catalogue of stars with properties of their H$\alpha$ emission line profiles, morphological classification, analysis of variability with time and the supplementary information from the SIMBAD, VizieR, and ADS databases. The records in SIMBAD indicate the presence of H$\alpha$ emission for roughly 25% of all stars in our catalogue, while at least 305 of them have already been more thoroughly investigated according to the references in ADS. The most frequently identified morphological categories in our sample of spectra are emission blend (23%), emission in absorption (22%), and self--absorption (16%). Objects with repeated observations demonstrate that our classification into discrete categories is generally stable through time, but categories P-Cygni and self--absorption seem less stable, which is the consequence of discrete classification rules, as well as of the fundamental change in profile shape.
Rare event search experiments using liquid xenon as target and detection medium require ultra-low background levels to fully exploit their physics potential. Cosmic ray induced activation of the detector components and, even more importantly, of the xenon itself during production, transportation and storage at the Earth's surface, might result in the production of radioactive isotopes with long half-lives, with a possible impact on the expected background. We present the first dedicated study on the cosmogenic activation of xenon after 345 days of exposure to cosmic rays at the Jungfraujoch research station at 3470m above sea level, complemented by a study of copper which has been activated simultaneously. We have directly observed the production of 7Be, 101Rh, 125Sb, 126I and 127Xe in xenon, out of which only 125Sb could potentially lead to background for a multi-ton scale dark matter search. The production rates for five out of eight studied radioactive isotopes in copper are in agreement with the only existing dedicated activation measurement, while we observe lower rates for the remaining ones. The specific saturation activities for both samples are also compared to predictions obtained with commonly used software packages, where we observe some underpredictions, especially for xenon activation.
The Star Formation Rate (SFR) is one of the main parameters used to analyze the evolution of galaxies through time. The need for recovering the light reprocessed by dust commonly requires the use of low spatial resolution far-infrared data. Recombination-line luminosities provide an alternative, although uncertain dust-extinction corrections based on narrow-band imaging or long-slit spectroscopy have traditionally posed a limit to their applicability. Integral Field Spectroscopy (IFS) is clearly the way to overcome such limitation. We obtain integrated H{\alpha}, ultraviolet (UV) and infrared (IR)-based SFR measurements for 272 galaxies from the CALIFA survey at 0.005 < z < 0.03 using single-band and hybrid tracers. We provide updated calibrations, both global and split by properties (including stellar mass and morphological type), referred to H{\alpha}. The extinction-corrected H{\alpha} luminosity agrees with the updated hybrid SFR estimators based on either UV or H{\alpha} plus IR luminosity over the full range of SFRs (0.03-20 M$_{\odot}$ yr$^{-1}$). The coefficient that weights the amount of energy produced by newly-born stars that is reprocessed by dust on the hybrid tracers, a$_{IR}$, shows a large dispersion. However, it does not became increasingly small at high attenuations, as expected if significant highly-obscured H$\alpha$ emission would be missed. Lenticulars, early-type spirals and type-2 AGN host galaxies show smaller coefficients due to the contribution of optical photons and AGN to dust heating. In the Local Universe the H{\alpha} luminosity derived from IFS observations can be used to measure SFR, at least in statistically-significant, optically-selected galaxy samples. The analysis of the SFR calibrations by galaxies properties could be potentially used by other works to study the impact of different selection criteria in the SFR values derived.
Gamma-ray bursts (GRBs), which have isotropic energy up to $10^{54}$ erg, would be the ideal tool to study the properties of early universe: including dark energy, star formation rate, and the metal enrichment history of the Universe. We will briefly review the progress on the field of GRB cosmology. Meanwhile, X-ray flares, which may have important clues to the central engine, are common phenomena in the GRB afterglows. We present statistical results of X-ray flares, i.e., energy, duration time and waiting time distributions, and compare the results with solar flares. The similarity between the two kinds of flares are found, which may indicates that the physical mechanism of GRB X-ray flares is magnetic reconnection.
Photon imaging for MeV gammas has serious difficulties due to huge backgrounds and unclearness in images, which are originated from incompleteness in determining the physical parameters of Compton scattering in detection, e.g., lack of the directional information of the recoil electrons. The recent major mission/instrument in the MeV band, CGRO/COMPTEL, which was Compton Camera (CC), detected mere $\sim30$ persistent sources. It is in stark contrast with $\sim$2000 sources in the GeV band. Here we report the performance of an Electron-Tracking Compton Camera (ETCC), and prove that it has a good potential to break through this stagnation in MeV gamma-ray astronomy. The ETCC provides all the parameters of Compton-scattering by measuring 3-D recoil electron tracks; then the SPD (Scatter Plane Deviation) lost in CCs is recovered. The energy loss rate (dE/dx), which CCs cannot measure, is also obtained, and is found to be indeed helpful to reduce the background under conditions similar to space. Accordingly the significance in gamma detection is improved severalfold. On the other hand, SPD is essential to determine the Point Spread Function quantitatively. The SPD resolution is improved close to the theoretical limit for multiple scattering of recoil electrons. With such a well-determined PSF, we demonstrate for the first time that it is possible to provide reliable sensitivity in Compton imaging without utilising an optimization algorithm. As such, this study highlights the fundamental weak-points of CCs. In contrast we demonstrate the possibility of ETCC reaching the sensitivity below $1\times10^{-12}$ erg cm$^{-2}$ s$^{-1}$ at 1 MeV.
Radio giant pulses provide a unique opportunity to study the pulsar radio emission mechanism in exquisite detail. Previous studies have revealed a wide range of properties and phenomena, including extraordinarily high brightness temperatures, sub-nanosecond emission features, and banded dynamic spectra. New measurements of giant pulse characteristics can help guide and test theoretical emission models. To this end, an extensive observation campaign has begun which will provide more than 500 hours on the Crab with a 34-meter antenna located in California, USA. The observations are being done as part of an educational outreach program called the Goldstone-Apple Valley Radio Telescope (GAVRT). This antenna has a novel wide bandwidth receiver which provides up to 8 GHz of instantaneous bandwidth in the range of 2.5 to 14 GHz. These observations will provide detailed information about the variability, amplitude distribution, and detailed frequency structure of radio giant pulses. In addition, a database of pulses from these observations and others of the Crab pulsar is being created which will simplify multiwavelength correlation analysis.
The adiabatic component of perturbations is damped during the kinetic decoupling due to the collision with relativistic component on sub-horizon scales. However the isocurvature part is free from the damping and could be large enough to make a substantial contribution to the formation of small scale structure. We explicitly study the weakly interacting massive particles as dark matter with an early matter dominated period before radiation domination and show that the isocurvature perturbation is generated during the phase transition and leaves imprint in the observable signatures for the small scale structure.
We introduce "anamorphic" cosmology, an approach for explaining the smoothness and flatness of the universe on large scales and the generation of a nearly scale-invariant spectrum of adiabatic density perturbations. The defining feature is a smoothing phase that acts like a contracting universe based on some Weyl frame-invariant criteria and an expanding universe based on other frame-invariant criteria. An advantage of the contracting aspects is that it is possible to avoid the multiverse and measure problems that arise in inflationary models. Unlike ekpyrotic models, anamorphic models can be constructed using only a single field and can generate a nearly scale-invariant spectrum of tensor perturbations. Anamorphic models also differ from pre-big bang and matter bounce models that do not explain the smoothness. We present some examples of cosmological models that incorporate an anamorphic smoothing phase.
The success of the first measurement of the light bending by the solar gravitational field is due to the particular stellar field during the Eddington's 1919 total eclipse of the Sun, near the Hyades, giving the opportunity to measure the gravitational bending of the light to the astronomers in two expeditions in Brazil, Sobral, and on the Principe Island in the Atlantic Ocean. The geometrical properties of this field and another field in Leo are discussed in view of repeating this experiment of General Relativity with SOHO satellite data in the context of the International Year of Light 2015.
Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are predicted to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. Here we summarize the results of the first systematic search for gamma-ray emission in Fermi-LAT data in the energy range from 100 MeV to 300 GeV from a large sample of SNe exploding in dense CSM. We search for a gamma-ray excess at the position of 147 SNe Type IIn in a one year time window after the optical peak time. In addition we combine the closest and optically brightest sources of our sample in a joint likelihood analysis in three different time windows (3, 6 and 12 months). No excess gamma-ray emission is found and limits on the gamma-ray luminosity and the ratio of gamma-ray to optical luminosity are presented.
Galaxy shapes are subject to distortions due to the tidal field of the Universe. The cross-correlation of galaxy lensing with the lensing of the Cosmic Microwave Background (CMB) cannot easily be separated from the cross-correlation of galaxy intrinsic shapes with CMB lensing. Previous work suggested that the intrinsic alignment contamination can be $15\%$ of this cross-spectrum for the CFHT Stripe 82 (CS82) and Atacama Cosmology Telescope surveys. Here we re-examine these estimates using up-to-date observational constraints of intrinsic alignments at a redshift more similar to that of CS82 galaxies. We find a $\approx$ $10\%$ contamination of the cross-spectrum from red galaxies, with $\approx$ $3\%$ uncertainty due to uncertainties in the redshift distribution of source galaxies and the modelling of the spectral energy distribution. Blue galaxies are consistent with being unaligned, but could contaminate the cross-spectrum by an additional $9.5\%$ within current $95\%$ confidence levels. While our fiducial estimate of alignment contamination is similar to previous work, our work suggests that the relevance of alignments for CMB lensing-galaxy lensing cross-correlation remains largely unconstrained. Little information is currently available about alignments at $z>1.2$. We consider the upper limiting case where all $z>1.2$ galaxies are aligned with the same strength as low redshift luminous red galaxies, finding as much as $\approx$ $60\%$ contamination.
A gamma-ray burst (GRB) is a strong and fast gamma-ray emission from the explosion of stellar systems (massive stars or coalescing binary compact stellar remnants), happening at any possible redshift, and detected by space missions. Although GRBs are the most energetic events after the Big Bang, systematic search (started after the first localization in 1997) led to only 374 spectroscopic redshift measurements. For less than half, the host galaxy is detected and studied in some detail. Despite the small number of known hosts, their impact on our understanding of galaxy formation and evolution is immense. These galaxies offer the opportunity to explore regions which are observationally hostile, due to the presence of gas and dust, or the large distances reached. The typical long-duration GRB host galaxy at low redshift is small, star-forming and metal poor, whereas, at intermediate redshift, many hosts are massive, dusty and chemically evolved. Going even farther in the past of the Universe, at z > 5, long-GRB hosts have never been identified, even with the deepest NIR space observations, meaning that these galaxies are very small (stellar mass < 10^7 M_sun). We considered the possibility that some high-z GRBs occurred in primordial globular clusters, systems that evolved drastically since the beginning, but would have back then the characteristics necessary to host a GRB. At that time, the fraction of stellar mass contained in proto globular clusters might have been orders of magnitude higher than today. Plus, these objects contained in the past many massive fast rotating binary systems, which are also regarded as a favorable situation for GRBs. The common factor for all long GRBs at any redshift is the stellar progenitor: it is a very massive rare/short-lived star, present in young regions, whose redshift evolution is closely related to the star-formation history of the Universe.
Cosmic structures leave an imprint on the microwave background radiation through the integrated Sachs-Wolfe effect. We construct a template map of the linear signal using the SDSS-III Baryon Acoustic Oscillation Survey at redshift 0.43 < z < 0.65. We verify the imprint of this map on the Planck CMB temperature map at the 97% confidence level and show consistency with the density-temperature cross-correlation measurement. Using this ISW reconstruction as a template we investigate the presence of ISW sources and further examine the properties of the Granett-Neyrinck-Szapudi supervoid and supercluster catalogue. We characterise the three-dimensional density profiles of these structures for the first time and demonstrate that they are significant structures. Model fits demonstrate that the supervoids are elongated along the line-of-sight and we suggest that this special orientation may be picked out by the void-finding algorithm in photometric redshift space. We measure the mean temperature profiles in Planck maps from public void and cluster catalogues. In an attempt to maximise the stacked ISW signal we construct a new catalogue of super-structures based upon local peaks and troughs of the gravitational potential. However, we do not find a significant correlation between these structures and the CMB temperature.
{\it Kepler} satellite photometry and phase-resolved spectroscopy of the ultracompact AM CVn type binary SDSS J190817.07+394036.4 are presented. The average spectra reveal a variety of weak metal lines of different species, including silicon, sulphur and magnesium as well as many lines of nitrogen, beside the strong absorption lines of neutral helium. The phase-folded spectra and the Doppler tomograms reveal an S-wave in emission in the core of the He I 4471 \AA\,absorption line at a period of $P_{\rm orb}=1085.7\pm2.8$\,sec identifying this as the orbital period of the system. The Si II, Mg II and the core of some He I lines show an S-wave in absorption with a phase offset of $170\pm15^\circ$ compared to the S-wave in emission. The N II, Si III and some helium lines do not show any phase variability at all. The spectroscopic orbital period is in excellent agreement with a period at $P_{\rm orb}=1085.108(9)$\,sec detected in the three year {\it Kepler} lightcurve. A Fourier analysis of the Q6 to Q17 short cadence data obtained by {\it Kepler} revealed a large number of frequencies above the noise level where the majority shows a large variability in frequency and amplitude. In an O-C analysis we measured a $\vert\dot{P}\vert\sim1.0\,$x$\,10^{-8}\,$s\,s$^{-1}$ for some of the strongest variations and set a limit for the orbital period to be $\vert\dot{P}\vert<10^{-10}$s\,s$^{-1}$. The shape of the phase folded lightcurve on the orbital period indicates the motion of the bright spot. Models of the system were constructed to see whether the phases of the radial velocity curves and the lightcurve variation can be combined to a coherent picture. However, from the measured phases neither the absorption nor the emission can be explained to originate in the bright spot.
We study the magnetic field generation in a neutron star within the model based on the magnetic field instability in the nuclear matter owing to the electron-nucleon parity violating interaction. We suggest that the growing magnetic field takes the energy from thermal background fermions in the neutron star matter. The system of kinetic equations for the spectra of the magnetic helicity density and magnetic energy density as well as the chiral imbalance are solved numerically accounting for this energy source. We obtain that, for the initial conditions corresponding to a typical neutron star, the large scale magnetic field $\sim 10^{15}\thinspace\text{G}$ is generated during $(10^4-10^5)\thinspace\text{yr}$. We suggest that the proposed model describes strong magnetic fields observed in magnetars.
We study the properties of satellites in the environment of massive star-forming galaxies at z~1.8 in the COSMOS field, using a sample of 215 galaxies on the main sequence of star formation with an average mass of 10^11 Msun. At z>1.5, these galaxies typically trace halos of mass >10^13 Msun. We use optical-near-infrared photometry to estimate stellar masses and star formation rates (SFR) of centrals and satellites down to ~6*10^9 Msun. We stack data around 215 central galaxies to statistically detect their satellite halos, finding an average of ~3 galaxies in excess of the background density. We fit the radial profiles of satellites with simple beta-models, and compare their integrated properties to model predictions. We find that the total stellar mass of satellites amounts to 68% of the central galaxy, while SED modeling and far-infrared photometry consistently show their total SFR to be 25-35% of the central's rate. We also see significant variation in the specific SFR of satellites within the halo with, in particular, a sharp decrease at <100 kpc. After considering different potential explanations, we conclude that this is likely an environmental signature of the hot inner halo. This effect can be explained in the first order by a simple free-fall scenario, suggesting that these low-mass environments can shut down star formation in satellites on relatively short timescales of ~0.3 Gyr.
Far away from any sunspot, a bright flare erupted on November 1st, 2014, with onset at 4:44 UT and a duration of around three hours, causing a C2.7-class flare. The blast was associated with the sudden disappearance of a large dark solar filament. The rest of the filament flew out into space, forming the core of a massive CME. Despite the location of the explosion over the sun's southeastern region (near the eastern edge of the sun) not be geoeffective, a radiation storm, that is, solar energetic particles (SEP) started to reach the Earth around 14:00 UT, reaching the condition of an S1 (minor) radiation storm level on Nov. 2th. In coincidence with onset of the S1 radiation storm (SEP above 5 MeV), the Tupi telescopes located at $22^090'$S; $43^020'$W, within the South Atlantic Anomaly (SAA) detected a muon enhancement caused by relativistic protons from this solar blast. In addition an increase in the particle intensity was found also at South Pole neutron monitor. This means that there was a transverse propagation to the interplanetary magnetic field of energetic solar particles. However, we show that perpendicular diffusion alone cannot explain these observations, it is necessary a combination with further processes as a very high speed, at least of a fraction the CME shocks, close to the ecliptic plane.
Both astronomical observations of the interaction of Type II supernova remnants (SNR) with dense interstellar clouds as well as cosmochemical studies of the abundances of daughter products of short-lived radioisotopes (SLRIs) formed by supernova nucleosynthesis support the hypothesis that the Solar Systems SLRIs may have been derived from a supernova. This paper continues a series devoted to examining whether such a shock wave could have triggered the dynamical collapse of a dense, presolar cloud core and simultaneously injected sufficient abundances of SLRIs to explain the cosmochemical evidence. Here we examine the effects of shock waves striking clouds whose spin axes are oriented perpendicular, rather than parallel, to the direction of propagation of the shock front. The models start with 2.2 solar mass cloud cores and shock speeds of 20 or 40 km/sec. Central protostars and protoplanetary disks form in all models, though with disk spin axes aligned somewhat randomly. The disks derive most of their angular momentum not from the initial cloud rotation, but from the Rayleigh-Taylor fingers that also inject shock wave SLRIs. Injection efficiencies, fi, the fraction of the incident shock wave material injected into the collapsing cloud core, are 0.04 - 0.1 in these models, similar to when the rotation axis is parallel to the shock propagation direction. Evidently altering the rotation axis orientation has only a minor effect on the outcome, strengthening the case for this scenario as an explanation for the Solar Systems SLRIs.
Double-Fourier interferometry is the most viable path to sub-arcsecond spatial resolution for future astronomical instruments that will observe the universe at far-infrared wavelengths. The double transform spatio-spectral interferometry couples pupil plane beam combination with detector arrays to enable imaging spectroscopy of wide fields, that will be key to accomplishing top-level science goals. The wide field of view and the necessity for these instruments to fly above the opaque atmosphere create unique characteristics and requirements compared to instruments on ground-based telescopes. In this paper, we discuss some characteristics of single-baseline spatio-spectral interferometers. We investigate the impact of intensity and optical path difference noise on the interferogram and the spectral signal-to-noise ratio. We apply our findings to the special case of the Balloon Experimental Twin Telescope for Infrared Interferometry (BETTII), a balloon payload that will be a first application of this technique at far-infrared wavelengths on a flying platform.
Detailed characterization of an extrasolar planet's atmosphere provides the best hope for distinguishing the makeup of its outer layers, and the only hope for understanding the interplay between initial composition, chemistry, dynamics & circulation, and disequilibrium processes. In recent years, some areas have seen rapid progress while developments in others have come more slowly and/or have been hotly contested. This article gives an observer's perspective on the current understanding of extrasolar planet atmospheres prior to the considerable advances expected from the next generation of observing facilities. Atmospheric processes of both transiting and directly-imaged planets are discussed, including molecular and atomic abundances, cloud properties, thermal structure, and planetary energy budgets. In the future we can expect a continuing and accelerating stream of new discoveries, which will fuel the ongoing exoplanet revolution for many years to come.
We calculate the upper bound on the reheating temperature given the non-observation of gravitational waves if the number of efolds during inflation are the minimum number required to address the horizon problem as formulated in terms of entropy. This bound is valid for canonical single field slow roll inflation with a generic potential. Our bound numerically is $T_{\text{reh}}\lesssim1.7\times10^{13}$ GeV, which is a factor of 428 less the usual bound one obtains from the non-observation of gravitational waves alone. If inflation lasted much longer than the minimum number of required efolds, our bound relaxes to coincide with the usual bound. We discuss the relevance for studies of primordial black holes.
We investigate if the hemispherical asymmetry in the CMB is produced from parity-violating excited initial condition. We show that in the limit where the deviations from the Bunch-Davies vacuum is large and the scale of new physics is maximally separated from the inflationary Hubble parameter, the primordial power spectrum is modulated only by dipole and quadrupole terms. Requiring the dipole contribution in the power spectrum accounts for the observed power asymmetry, $A=0.07\pm0.022$, we show that the amount of quadrupole terms is roughly equal to $A^2$, which is still consistent with the bounds from the CMB. The mean local bispectrum which gets enhanced for the excited initial states is within the $1\sigma$ bound of Planck 2015 results, $f_{\rm NL}\simeq 4.17$, but reachable by future CMB experiments. The amplitude of the local non-gaussianity modulates around this mean value, approximately depending on the angle that the short wavelength mode makes with the preferred direction. The amount of variation maximizes for the configurations that are coplanar with the preferred direction. For counterclockwise oriented configurations, maximum and minimum values for the non-gaussianity will occur for the ones in which the short wavelength mode is, respectively, antiparallel and parallel with the preferred direction. The difference of non-gaussianity between these two configurations is as large as $\simeq 1.2$ which can be used to distinguish this scenario from other scenarios that try to explain the observed hemispherical asymmetry. Such modulation in non-gaussianity is minimized for the configurations that are in the equator plane orthogonal to the preferred direction.
We investigate the inflationary implications of extensions of Poincare symmetry. The simplest constructions with local scale invariance lead to universal predictions: the spectral index is $n_s = 1-2/N$, in excellent agreement with Planck data, while the tensor-to-scalar ratio is determined by a free parameter to $r = 12 \alpha / N^2$. For the special value $\alpha=1$ one finds symmetry enhancement to the full conformal group. We show that these findings hold both for two-derivative scalar-tensor theories as well as higher-derivative gravity. Therefore scale invariance underlies a promising set of inflationary models.
We define a generalization of scalar fields with non-canonical kinetic term which we call exotic k-essence or briefly, exotik. These fields are generated by the global description of cosmological models with two interactive fluids in the dark sector and under certain conditions, they correspond to usual k-essences. The formalism is applied to the cases of constant potential and of inverse square potential and also we develop the purely exotik version for the modified holographic Ricci type of dark energy (MHR), where the equations of state are not constant. With the kinetic function $F=1+mx$ and the inverse square potential we recover, through the interaction term, the identification between k-essences and quintessences of exponential potential, already known for Friedmann-Robertson-Walker and Bianchi type I geometries. Worked examples are shown that include the self-interacting MHR and also models with crossing of the phantom divide line (PDL).
The majority of the matter in the universe is still unidentified and under investigation by both direct and indirect means. Many experiments searching for the recoil of dark-matter particles off target nuclei in underground laboratories have established increasingly strong constraints on the mass and scattering cross sections of weakly interacting particles, and some have even seen hints at a possible signal. Other experiments search for a possible mixing of photons with light scalar or pseudo-scalar particles that could also constitute dark matter. Furthermore, annihilation or decay of dark matter can contribute to charged cosmic rays, photons at all energies, and neutrinos. Many existing and future ground-based and satellite experiments are sensitive to such signals. Finally, data from the Large Hadron Collider at CERN are scrutinized for missing energy as a signature of new weakly interacting particles that may be related to dark matter. In this review article we summarize the status of the field with an emphasis on the complementarity between direct detection in dedicated laboratory experiments, indirect detection in the cosmic radiation, and searches at particle accelerators.
Neutrons produced by the carbon fusion reaction 12C(12C,n)23Mg play an important role in stellar nucleosynthesis. However, past studies have shown large discrepancies between experimental data and theory, leading to an uncertain cross section extrapolation at astrophysical energies. We present the first direct measurement that extends deep into the astrophysical energy range along with a new and improved extrapolation technique based on experimental data from the mirror reaction 12C(12C,p)23Na. The new reaction rate has been determined with a well-defined uncertainty that exceeds the precision required by astrophysics models. Using our constrained rate, we find that 12C(12C,n)23Mg is crucial to the production of Na and Al in Pop-III Pair Instability Supernovae. It also plays a non-negligible role in the production of weak s-process elements as well as in the production of the important galactic gamma-ray emitter 60Fe.
Links to: arXiv, form interface, find, astro-ph, recent, 1507, contact, help (Access key information)
The major mass fraction of the envelope of hot luminous stars is radiatively stable. However, the partial ionisation of hydrogen, helium and iron gives rise to extended sub-surface convection zones in all of them. In this work, we investigate the effect of the pressure induced by the turbulent motion in these zones based on the mixing length theory, and search for observable consequences. We find that the turbulent pressure fraction can amount up to ~5% in OB supergiants, and to ~30% in cooler supergiants. The resulting structural changes are, however, not significantly affecting the evolutionary tracks compared to previous calculations. Instead, a comparison of macroturbulent velocities derived from high quality spectra of OB stars with the turbulent pressure fraction obtained in corresponding stellar models reveals a strong correlation of these two quantities. We discuss a possible physical connection, and conclude that turbulent pressure fluctuations may drive high-order oscillations, which - as conjectured earlier - manifest themselves as macroturbulence in the photospheres of hot luminous stars.
We report on an informal survey about the use of software in the worldwide astronomical community. The survey was carried out between December 2014 and February 2015, collecting responses from 1142 astronomers, spanning all career levels. We find that all participants use software in their research. The vast majority of participants, 90%, write at least some of their own software. Even though writing software is so wide-spread among the survey participants, only 8% of them report that they have received substantial training in software development. Another 49% of the participants have received "little" training. The remaining 43% have received no training. We also find that astronomers' software stack is fairly narrow. The 10 most popular tools among astronomers are (from most to least popular): Python, shell scripting, IDL, C/C++, Fortran, IRAF, spreadsheets, HTML/CSS, SQL and Supermongo. Across all participants the most common programing language is Python ($67\pm 2\%$), followed by IDL ($44\pm 2\%$), C/C++ ($37\pm 2\%$) and Fortran ($28\pm 2\%$). IRAF is used frequently by $24\pm 1\%$ of participants. We show that all trends are largely independent of career stage, area of research and geographic location.
Evidence for an extraterrestrial flux of high-energy neutrinos has now been found in multiple searches with the IceCube detector. The first solid evidence was provided by a search for neutrino events with deposited energies $\gtrsim30$~TeV and interaction vertices inside the instrumented volume. Recent analyses suggest that the extraterrestrial flux extends to lower energies and is also visible with throughgoing, $\nu_\mu$-induced tracks from the Northern hemisphere. Here, we combine the results from six different IceCube searches for astrophysical neutrinos in a maximum-likelihood analysis. The combined event sample features high-statistics samples of shower-like and track-like events. The data are fit in up to three observables: energy, zenith angle and event topology. Assuming the astrophysical neutrino flux to be isotropic and to consist of equal flavors at Earth, the all-flavor spectrum with neutrino energies between 25~TeV and 2.8~PeV is well described by an unbroken power law with best-fit spectral index $-2.50\pm0.09$ and a flux at 100~TeV of $\left(6.7_{-1.2}^{+1.1}\right)\cdot10^{-18}\,\mathrm{GeV}^{-1}\mathrm{s}^{-1}\mathrm{sr}^{-1}\mathrm{cm}^{-2}$. Under the same assumptions, an unbroken power law with index $-2$ is disfavored with a significance of 3.8~$\sigma$ ($p=0.0066\%$) with respect to the best fit. This significance is reduced to 2.1~$\sigma$ ($p=1.7\%$) if instead we compare the best fit to a spectrum with index $-2$ that has an exponential cut-off at high energies. Allowing the electron neutrino flux to deviate from the other two flavors, we find a $\nu_e$~fraction of $0.18\pm0.11$ at Earth. The sole production of electron neutrinos, which would be characteristic of neutron-decay dominated sources, is rejected with a significance of 3.6~$\sigma$ ($p=0.014\%$).
We present a weak-lensing analysis of the merging {\em Frontier Fields} (FF) cluster Abell~2744 using new Subaru/Suprime-Cam imaging. The wide-field lensing mass distribution reveals this cluster is comprised of four distinct substructures. Simultaneously modeling the two-dimensional reduced shear field using a combination of a Navarro--Frenk--White (NFW) model for the main core and truncated NFW models for the subhalos, we determine their masses and locations. The total mass of the system is constrained as $M_\mathrm{200c} = (2.06\pm0.42)\times10^{15}\,M_\odot$. The most massive clump is the southern component with $M_\mathrm{200c} = (7.7\pm3.4)\times10^{14}\,M_\odot$, followed by the western substructure ($M_\mathrm{200c} = (4.5\pm2.0)\times10^{14}\,M_\odot$) and two smaller substructures to the northeast ($M_\mathrm{200c} = (2.8\pm1.6)\times10^{14}\,M_\odot$) and northwest ($M_\mathrm{200c} = (1.9\pm1.2)\times10^{14}\,M_\odot$). The presence of the four substructures supports the picture of multiple mergers. Using a composite of hydrodynamical binary simulations we explain this complicated system without the need for a "slingshot" effect to produce the northwest X-ray interloper, as previously proposed. The locations of the substructures appear to be offset from both the gas ($87^{+34}_{-28}$ arcsec, 90\% CL) and the galaxies ($72^{+34}_{-53}$ arcsec, 90\% CL) in the case of the northwestern and western subhalos. To confirm or refute these findings, high resolution space-based observations extending beyond the current FF limited coverage to the west and northwestern area are essential.
In quiescent low-mass X-ray binaries (qLMXBs) containing neutron stars, the origin of the thermal X-ray component may be either release of heat from the core of the neutron star, or continuing low-level accretion. In general, heat from the core should be stable on timescales $<10^4$ years, while continuing accretion may produce variations on a range of timescales. While some quiescent neutron stars (e.g. Cen X-4, Aql X-1) have shown variations in their thermal components on a range of timescales, several others, particularly those in globular clusters with no detectable nonthermal hard X-rays (fit with a powerlaw), have shown no measurable variations. Here, we constrain the spectral variations of 12 low mass X-ray binaries in 3 globular clusters over $\sim10$ years. We find no evidence of variations in 10 cases, with limits on temperature variations below 11% for the 7 qLMXBs without powerlaw components, and limits on variations below 20% for 3 other qLMXBs that do show non-thermal emission. However, in 2 qLMXBs showing powerlaw components in their spectra (NGC 6440 CX 1 & Terzan 5 CX 12) we find marginal evidence for a 10% decline in temperature, suggesting the presence of continuing low-level accretion. This work adds to the evidence that the thermal X-ray component in quiescent neutron stars without powerlaw components can be explained by heat deposited in the core during outbursts. Finally, we also investigate the correlation between hydrogen column density (N$_H$) and optical extinction (A$_V$) using our sample and current models of interstellar X-ray absorption, finding $N_H ({\rm cm}^{-2}) = (2.81\pm0.13)\times10^{21} A_V$.
We compare the dark matter halos' structural parameters derived for four Milky Way dwarf spheroidal galaxies to those of subhalos found in cosmological $N$-body simulations. We confirm that estimates of the mass at a single fixed radius are fully consistent with the observations. However, when a second structural parameter such as the logarithmic slope of the dark halo density profile measured close to the half-light radius is included in the comparison, we find little to no overlap between the satellites and the subhalos. Typically the right mass subhalos have steeper profiles at these radii than measurements of the dSph suggest. Using energy arguments we explore if it is possible to solve this discrepancy by invoking baryonic effects. Assuming that feedback from supernovae can lead to a reshaping of the halos, we compute the required efficiency and find entirely plausible values for a significant fraction of the subhalos and even as low as 0.1%. This implies that care must be taken not to exaggerate the effect of supernovae feedback as this could make the halos too shallow. These results could be used to calibrate and possibly constrain feedback recipes in hydrodynamical simulations.
Magnetohydrodynamic accretion disk simulations suggest that much of the energy liberated by the magnetorotational instability (MRI) can be channeled into large-scale toroidal magnetic fields through dynamo action. Under certain conditions, this field can dominate over gas and radiation pressure in providing vertical support against gravity, even close to the midplane. Using a simple model for the creation of this field, its buoyant rise, and its coupling to the gas, we show how disks could be driven into this magnetically dominated state and deduce the resulting vertical pressure and density profiles. Applying an established criterion for MRI to operate in the presence of a toroidal field, we show that magnetically supported disks can have two distinct MRI-active regions, separated by a "dead zone" where local MRI is suppressed, but where magnetic energy continues to flow upward from the dynamo region below. We suggest that the relative strengths of the MRI zones, and the local poloidal flux, determine the spectral states of X-ray binaries. Specifically, "intermediate" and "hard" accretion states occur when MRI is triggered in the hot, upper zone of the corona, while disks in "soft" states do not develop the upper MRI zone. We discuss the conditions under which various transitions should take place and speculate on the relationship of dynamo activity to the various types of quasi-periodic oscillations that sometimes appear in the hard spectral components. The model also explains why luminous accretion disks in the "soft" state show no signs of the thermal/viscous instability predicted by standard alpha models.
Through our HARPS radial velocity survey for planets around solar twin stars, we have identified a promising Jupiter twin candidate around the star HIP11915. We characterize this Keplerian signal and investigate its potential origins in stellar activity. Our analysis indicates that HIP11915 hosts a Jupiter-mass planet with a 3600-day orbital period and low eccentricity. Although we cannot definitively rule out an activity cycle interpretation, we find that a planet interpretation is more likely based on a joint analysis of RV and activity index data. The challenges of long-period radial velocity signals addressed in this paper are critical for the ongoing discovery of Jupiter-like exoplanets. If planetary in nature, the signal investigated here represents a very close analog to the solar system in terms of both Sun-like host star and Jupiter-like planet.
We present Ha maps at 1kpc spatial resolution for star-forming galaxies at z~1, made possible by the WFC3 grism on HST. Employing this capability over all five 3D-HST/CANDELS fields provides a sample of 2676 galaxies. By creating deep stacked Halpha (Ha) images, we reach surface brightness limits of 1x10^-18\erg\s\cm^2\arcsec^2, allowing us to map the distribution of ionized gas out to >10kpc for typical L* galaxies at this epoch. We find that the spatial extent of the Ha distribution increases with stellar mass as r(Ha)[kpc]=1.5(Mstars/10^10Msun)^0.23. Furthermore, the Ha emission is more extended than the stellar continuum emission, consistent with inside-out assembly of galactic disks. This effect, however, is mass dependent with r(Ha)/r(stars)=1.1(M/10^10Msun)^0.054, such that at low masses r(Ha)~r(stars). We map the Ha distribution as a function of SFR(IR+UV) and find evidence for `coherent star formation' across the SFR-M plane: above the main sequence, Ha is enhanced at all radii; below the main sequence, Ha is depressed at all radii. This suggests that at all masses the physical processes driving the enhancement or suppression of star formation act throughout the disks of galaxies. It also confirms that the scatter in the star forming main sequence is real and caused by variations in the star formation rate at fixed mass. At high masses (10^10.5<M/Msun<10^11), above the main sequence, Ha is particularly enhanced in the center, plausibly building bulges and/or supermassive black holes. Below the main sequence, the star forming disks are more compact and a strong central dip in the EW(Ha), and the inferred specific star formation rate, appears. Importantly though, across the entirety of the SFR-M plane, the absolute star formation rate as traced by Ha is always centrally peaked, even in galaxies below the main sequence.
Wolf-Rayet (WR) stars have a severe impact on their environments owing to their strong ionizing radiation fields and powerful stellar winds. Since these winds are considered to be driven by radiation pressure, it is theoretically expected that the degree of the wind mass-loss depends on the initial metallicity of WR stars. Following our comprehensive studies of WR stars in the Milky Way, M31, and the LMC, we derive stellar parameters and mass-loss rates for all seven putatively single WN stars known in the SMC. Based on these data, we discuss the impact of a low-metallicity environment on the mass loss and evolution of WR stars. The quantitative analysis of the WN stars is performed with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The physical properties of our program stars are obtained from fitting synthetic spectra to multi-band observations. In all SMC WN stars, a considerable surface hydrogen abundance is detectable. The majority of these objects have stellar temperatures exceeding 75 kK, while their luminosities range from 10^5.5 to 10^6.1 Lsun. The WN stars in the SMC exhibit on average lower mass-loss rates and weaker winds than their counterparts in the Milky Way, M31, and the LMC. By comparing the mass-loss rates derived for WN stars in different Local Group galaxies, we conclude that a clear dependence of the wind mass-loss on the initial metallicity is evident, supporting the current paradigm that WR winds are driven by radiation. A metallicity effect on the evolution of massive stars is obvious from the HRD positions of the SMC WN stars at high temperatures and high luminosities. Standard evolution tracks are not able to reproduce these parameters and the observed surface hydrogen abundances. Homogeneous evolution might provide a better explanation for their evolutionary past.
Results from the IceCube Neutrino Observatory have recently provided compelling evidence for the existence of a high energy astrophysical neutrino flux utilizing a dominantly Southern Hemisphere dataset consisting primarily of nu_e and nu_tau charged current and neutral current (cascade) neutrino interactions. In the analysis presented here, a data sample of approximately 35,000 muon neutrinos from the Northern sky was extracted from data taken during 659.5 days of livetime recorded between May 2010 and May 2012. While this sample is composed primarily of neutrinos produced by cosmic ray interactions in the Earth's atmosphere, the highest energy events are inconsistent with a hypothesis of solely terrestrial origin at 3.7 sigma significance. These neutrinos can, however, be explained by an astrophysical flux per neutrino flavor at a level of Phi(E_nu) = 9.9^{+3.9}_{-3.4} times 10^{-19} GeV^{-1} cm^{-2} sr^{-1} s^{-1} ({E_nu / 100 TeV})^{-2}, consistent with IceCube's Southern Hemisphere dominated result. Additionally, a fit for an astrophysical flux with an arbitrary spectral index was performed. We find a spectral index of 2.2^{+0.2}_{-0.2}, which is also in good agreement with the Southern Hemisphere result.
Chondrule formation remains one of the most elusive early Solar System events. Here, we take the novel approach of employing numerical simulations to investigate chondrule origin beyond purely cosmochemical methods. We model the transport of generically-produced chondrules and dust in a 1D viscous protoplanetary disk model, in order to constrain the chondrule formation events. For a single formation event we are able to match analytical predictions of the memory chondrule and dust populations retain of each other (complementarity), finding that a large mass accretion rate ($\gtrsim 10^{-7}$~M$_\odot$~yr$^{-1}$) allows for delays on the order of the disk's viscous timescale between chondrule formation and chondrite accretion. Further, we find older disks to be severely diminished of chondrules, with accretion rates $\lesssim 10^{-9}$~M$_\odot$~yr$^{-1}$ for nominal parameters. We then characterize the distribution of chondrule origins in both space and time, as functions of disk parameters and chondrule formation rates, in runs with continuous chondrule formation and both static and evolving disks. Our data suggest that these can account for the observed diversity between distinct chondrite classes, if some diversity in accretion time is allowed for.
Starburst regions with multiple powerful winds of young massive stars and supernova remnants are favorable sites for high-energy cosmic ray acceleration. A supernova shock colliding with a fast wind from a compact cluster of young stars allows the acceleration of protons to energies well above the standard limits of diffusive shock acceleration in an isolated SN. The proton spectrum in such a wind-supernova PeV accelerator is hard with a large flux in the high-energy-end of the spectrum producing copious gamma-rays and neutrinos in inelastic nuclear collisions. We argue that SN shocks in the Westerlund 1 cluster in the Milky Way may accelerate protons to about 40 PeV. Once accelerated, these CRs will diffuse into surrounding dense clouds and produce neutrinos with fluxes sufficient to explain a fraction of the events detected by IceCube Observatory from the inner Galaxy.
$\gamma$-ray radiation from pulsars is usually thought to be mostly produced by the synchro-curvature losses of accelerated particles. Here we present a systematic study of all currently reported, good-quality Fermi-LAT pulsar spectral data. We do so by applying a model which follows the particle dynamics and consistently computes the emission of synchro-curvature radiation. By fitting observational data on a case by case basis, we are able to obtain constraints about the parallel electric field, the typical lengthscale over which particles emit the bulk of the detected radiation, and the number of involved particles. The model copes well with data of several dozens of millisecond and young pulsars. By correlating the inferred model parameters with the observed timing properties, some trends are discovered. First, a non-negligible part of the radiation comes from the loss of perpendicular momentum soon after pair creation. Second, the electric field strongly correlates with both the inverse of the emission lengthscale and the magnetic field at light cylinder, thus ruling out models with high-energy photon production close to the surface. These correlations unify young and millisecond pulsars under the same physical scenario, and predict that magnetars are intrinsically $\gamma$-ray quiet via syncrhro-curvature processes, since magnetospheric particles are not accelerated enough to emit a detectable $\gamma$-ray flux.
In this paper we model the observed absorption measure distribution (AMD) in Mrk 509 with a single-zone absorber in pressure equilibrium. AMD is usually constructed from observations of narrow absorption lines in radio-quiet active galaxies with warm absorbers. We study the properties of the warm absorber in Mrk 509 using recently published broad-band spectral energy distribution observed with different instruments. This spectrum is an input in radiative transfer computations with full photoionisation treatment using titan code. We show that the simplest way to fully reproduce the shape of AMD is to assume that the warm absorber is a single zone under constant total pressure. With this assumption we found theoretical AMD which matches the observed one determined on the basis of 600 ks RGS XMM-Newton spectrum of Mrk 509. Our model puts strong constraints that the density of the warm absorber should be high enough to produce strong opacity jumps which are responsible for observed AMD dips.
We present results from multi-wavelength simultaneous X-ray and radio observations of the black hole X-ray binary V404 Cyg in quiescence. Our coverage with NuSTAR provides the very first opportunity to study the X-ray spectrum of V404 Cyg at energies above 10 keV. The unabsorbed broad-band (0.3-30 keV) quiescent luminosity of the source is 8.9$\times$10$^{32}$ erg s$^{-1}$ for a distance of 2.4 kpc. The source shows clear variability on short time scales in radio, soft X-ray and hard X-ray bands in the form of multiple flares. The broad-band X-ray spectra obtained from XMM-Newton and NuSTAR can be characterized with a power-law model having photon index {\Gamma}=2.13$\pm$0.07 (90% confidence errors); however, residuals at high energies indicate spectral curvature significant at a 3{\sigma} confidence level with e-folding energy of the cutoff to be 19$^{+19}_{-7}$ keV. Such curvature can be explained using synchrotron emission from the base of a jet outflow. Radio observations using the JVLA reveal that the spectral index evolves on very fast time-scales, switching between optically thick and thin synchrotron emission, possibly due to instabilities in the compact jet or stochastic instabilities in accretion rate. We explore different scenarios to explain this fast variability.
We report on the energetics of the 2012 March 12 polar coronal mass ejection (CME) originating from a southern latitude of ~60o. The polar CME is similar to low-latitude CMEs in almost all respects: three-part morphology, post eruption arcade (PEA), CME and filament kinematics, CME mass and kinetic energy, and the relative thermal energy content of the PEA. From polarized brightness images, we estimate the CME mass, which is close to the average mass of low-latitude CMEs. The CME kinetic energy (3.3x1030 erg) is also typical of the general population of CMEs. From photospheric magnetograms, we estimate the free energy (1.8x1031 erg) in the polar crown source region, which we find is sufficient to power the CME and the PEA. About 19% of the free energy went into the CME kinetic energy. We compute the thermal energy content of the PEA (2.3x1029 erg) and find it to be a small fraction (6.8%) of the CME kinetic energy. This fraction is remarkably similar to that in active region CMEs associated with major flares. We also show that the 2012 March 12 is one among scores of polar CMEs observed during the maximum phase of cycle 24. The cycle 24 polar crown prominence eruptions have the same rate of association with CMEs as those from low-latitudes. This investigation supports the view that all CMEs are magnetically propelled from closed field regions, irrespective of their location on the Sun (polar crown filament regions, quiescent filament regions or active regions).
Isolated galaxies in low-density regions are significant in the sense that they are least affected by the hierarchical pattern of galaxy growth, and interactions with perturbers, at least for the last few Gyr. To form a comprehensive picture of the star formation history of isolated galaxies, we constructed a catalog of isolated galaxies and their comparison sample in relatively denser environments. The galaxies are drawn from the SDSS DR7 in the redshift range of $0.025<z<0.044$. We performed a visual inspection and classified their morphology following the Hubble classification scheme. For the spectroscopic study, we make use of the OSSY catalog. We confirm most of the earlier understanding on isolated galaxies. The most remarkable additional results are as follows. Isolated galaxies are dominantly late type with the morphology distribution (E: S0: S: Irr) = (9.9: 11.3: 77.6: 1.2)\%. The frequency of elliptical galaxies among isolated galaxies is only a third of that of the comparison sample. Most of the photometric and spectroscopic properties are surprisingly similar between isolated and comparison samples. However, early-type isolated galaxies are less massive by 50\% and younger (by H$\beta$) by 20\% than their counterparts in the comparison sample. This can be explained as a result of different merger and star formation histories for differing environments in the hierarchical merger paradigm. We provide an on-line catalog for the list and properties of our sample galaxies.
Many jets are detected at X-ray wavelengths in the Sun's polar regions, and the ejected plasma along the jets has been suggested to contribute mass to the fast solar wind. From in-situ measurements in the magnetosphere, it has been found that the fast solar wind has photospheric abundances while the slow solar wind has coronal abundances. Therefore, we investigated the abundances of polar jets to determine whether they are the same as that of the fast solar wind. For this study, we selected 22 jets in the polar region observed by Hinode/EIS (EUV Imaging Spectrometer) and XRT (X-Ray Telescope) simultaneously on 2007 November 1-3. We calculated the First Ionization Potential (FIP) bias factor from the ratio of the intensity between high (S) and low (Si, Fe) FIP elements using the EIS spectra. The values of the FIP bias factors for the polar jets are around 0.7-1.9, and 75$\%$ of the values are in the range of 0.7-1.5, which indicates that they have photospheric abundances similar to the fast solar wind. The results are consistent with the reconnection jet model where photospheric plasma emerges and is rapidly ejected into the fast wind.
We analyze the Cosmic Microwave Background Radiation (CMBR) temperature and polarization data in order to extract the signal of correlation between l and l+1 multipoles in the multipole ranges, 2-64, 30-64 and 30-100. Such a correlation is predicted by the dipole modulation model proposed on the basis of the observed hemispherical anisotropy in temperature field. An anisotropic or inhomogeneous model of primordial power spectrum which leads to such correlations in temperature field also predicts similar correlations in CMBR polarization. Our results for the case of temperature using the latest PLANCK data agree with those obtained by earlier analysis. We also find a very strong signal of correlation in the polarization data. Surprisingly, however, the preferred direction in the case of polarization points in the direction close to the CMBR dipole which is very different from the corresponding direction in the case of temperature.
The paradigm for GRB prompt emission is changing. Since early in the CGRO era, the empirical Band function has been considered a good description of the keV-MeV spectra although its shape is very often inconsistent with the predictions of the pure synchrotron emission scenarios. We have recently established a new observational model analyzing data of the NASA Fermi Gamma-ray Space Telescope. In this model, GRB prompt emission is a combination of three main emission components: (i) a thermal-like component that we interpreted so far as the jet photosphere emission, (ii) a non-thermal component that we interpreted so far as synchrotron radiation, and (iii) an additional non-thermal (cutoff) power-law most likely of inverse Compton origin. In this article we reanalyze some of the bright GRBs observed with CGRO/BATSE with the new model, namely GRBs 941017, 970111 and 990123. We conclude that BATSE data are fully consistent with the recent results obtained with Fermi: some bright BATSE GRBs exhibit three separate components during the prompt phase with similar spectral parameters as these reported from Fermi data. In addition, the analysis of the BATSE GRBs with the new model results in a relation between the time-resolved energy flux of the non-thermal component and its corresponding $\nu$F$_\nu$ peak energy (i.e., F$_{i}^{nTh}$-E$_{peak,i}^{nTh}$) that has a similar index as the one initially derived from Fermi data. For GRBs with known redshift (z) this results in a possible universal relation between the luminosity of the non-thermal component and its corresponding $\nu$F$_\nu$ peak energy in the rest frame (i.e., L$_{i}^{nTh}$-E$_{peak,i}^{NT,rest}$). We estimated z for GRBs 941017 and 970111 using GRB 990123--with z=1.61--as a reference. The estimated z for GRB 941017 is typical for long GRBs and the estimated z for GRB 970111 is right in the range of the expected values for this burst.
We investigate the location of the radio jet bases ("radio cores") of blazars in radio images, and their stationarity by means of dense very long baseline interferometry (VLBI) observations. In order to measure the position of a radio core, we conducted 12 epoch astrometric observation of the blazar Markarian 421 with the VLBI Exploration of Radio Astrometry at 22 GHz immediately after a large X-ray flare, which occurred in the middle of 2011 September. For the first time, we find that the radio core is not stationary but rather changes its location toward 0.5 mas downstream. This angular scale corresponds to the de-projected length of a scale of $10^5$ Schwarzschild radii (Rs) at the distance of Markarian~421. This radio-core wandering may be a new type of manifestation associated with the phenomena of large X-ray flares.
Using star-forming galaxies sample in the nearby Universe (0.02<z<0.10) selected from the SDSS (DR7) and GALEX all-sky survey (GR5), we present a new empirical calibration for predicting dust extinction of galaxies from H-alpha-to-FUV flux ratio. We find that the H-alpha dust extinction (A(Ha)) derived with H-alpha/H-beta ratio (Balmer decrement) increases with increasing H-alpha/UV ratio as expected, but there remains a considerable scatter around the relation, which is largely dependent on stellar mass and/or H-alpha equivalent width (EW(Ha)). At fixed H-alpha/UV ratio, galaxies with higher stellar mass (or galaxies with lower EW(Ha)) tend to be more highly obscured by dust. We quantify this trend and establish an empirical calibration for predicting A(Ha) with a combination of H-alpha/UV ratio, stellar mass and EW(Ha), with which we can successfully reduce the systematic uncertainties accompanying the simple H-alpha/UV approach by ~15-30%. The new recipes proposed in this study will provide a convenient tool for predicting dust extinction level of galaxies particularly when Balmer decrement is not available. By comparing A(Ha) (derived with Balmer decrement) and A(UV) (derived with IR/UV luminosity ratio) for a subsample of galaxies for which AKARI FIR photometry is available, we demonstrate that more massive galaxies tend to have higher extra extinction towards the nebular regions compared to the stellar continuum light. Considering recent studies reporting smaller extra extinction towards nebular regions for high-redshift galaxies, we argue that the dust geometry within high-redshift galaxies resemble more like low-mass galaxies in the nearby Universe.
Galactic archaeology is the study of the history of star formation and chemical evolution in the Milky Way, based on present-day stellar populations. Studies of young stars are a key anchor point for Galactic archaeology, since quantities like the initial mass function and the star formation rate can be studied directly in young clusters and star forming regions. Conversely, massive spectroscopic Galactic archaeology surveys can be used as a data source for young star studies.
SDSS J082053.53+000843.4 belongs to the HW Vir family of short period binary systems and was first identified by Geier in 2011. Whilst three subsequent papers have focused on the morphology of this system, little has been published on the system period and its constancy. Here we provide the first published times of minima together with a revised ephemeris and a binary period, 0.0962404(1)d, that is four orders of magnitude more precise than previous declared values.
The underlying mechanisms driving the quenching of dwarf-mass satellite galaxies remain poorly constrained, but recent studies suggest they are particularly inefficient for those satellites with stellar mass 10$^{\rm 9}$ M$_{\odot}$. We investigate the characteristic evolution of these systems with chemodynamical simulations and idealised models of their tidal/hydrodynamic interactions within the 10$^{\rm 13-13.5}$ M$_{\odot}$ group-mass hosts in which they are preferentially quenched. Our fiducial simulations highlight the role played by secular star formation and stellar bars, and demonstrate a transition from a gas-rich to passive, HI-deficient state (i.e. $\Delta$SFR$\le$-1, def$_{\rm HI}$$\ge$0.5) within 6 Gyr of first infall. Furthermore, in the 8-10 Gyr in which these systems have typically been resident within group hosts, the bulge-to-total ratio of an initially bulgeless disc can increase to 0.3$<$B/T$<$0.4, its specific angular momentum $\lambda_{\rm R}$ reduce to $\sim$0.5, and strong bisymmetries formed. Ultimately, this scenario yields satellites resembling dwarf S0s, a result that holds for a variety of infall inclinations/harassments albeit with broad scatter. The key assumptions here lie in the rapid removal of the satellite's gaseous halo upon virial infall, and the satellite's local intra-group medium density being defined by the host's spherically-averaged profile. We demonstrate how quenching can be greatly enhanced if the satellite lies in an overdensity, consistent with recent cosmological-scale simulations but contrasting with observationally-inferred quenching mechanisms/timescales; an appraisal of these results with respect to the apparent preferential formation of dS0s/S0s in groups is also given.
We analyse the magnetic activity characteristics of the planet hosting Sun-like star, HD 1237, using HARPS spectro-polarimetric time-series data. We find evidence of rotational modulation of the magnetic longitudinal field measurements consistent with our ZDI analysis, with a period of 7 days. We investigate the effect of customising the LSD mask to the line depths of the observed spectrum and find that it has a minimal effect on shape of the extracted Stokes V profile but does result in a small increase in the S/N ($\sim$ 7%). We find that using a Milne-Eddington solution to describe the local line profile provides a better fit to the LSD profiles in this slowly rotating star, which also impacts the recovered ZDI field distribution. We also introduce a fit-stopping criterion based on the information content (entropy) of the ZDI maps solution set. The recovered magnetic field maps show a strong (+90 G) ring-like azimuthal field distribution and a complex radial field dominating at mid latitudes ($\sim$45 degrees). Similar magnetic field maps are recovered from data acquired five months apart. Future work will investigate how this surface magnetic field distribution impacts the coronal magnetic field and extended environment around this planet-hosting star.
We present self-consistent triaxial stellar systems that have analytic
distribution functions (DFs) expressed in terms of the actions. These provide
triaxial density profiles with cores or cusps at the centre. They are the first
self-consistent triaxial models with analytic DFs suitable for modelling giant
ellipticals and dark haloes. Specifically, we study triaxial models that
reproduce the Hernquist profile from Williams & Evans (2015), as well as
flattened isochrones of the form proposed by Binney (2014). We explore the
kinematics and orbital structure of these models in some detail. The models
typically become more radially anisotropic on moving outwards, have velocity
ellipsoids aligned in Cartesian coordinates in the centre and aligned in
spherical polar coordinates in the outer parts.
In projection, the ellipticity of the isophotes and the position angle of the
major axis of our models generally changes with radius. So, a natural
application is to elliptical galaxies that exhibit isophote twisting. As
triaxial St\"ackel models do not show isophote twists, our DFs are the first to
generate mass density distributions that do exhibit this phenomenon, typically
with a gradient of $\approx 10^\circ$/effective radius, which is comparable to
the data.
Triaxiality is a natural consequence of models that are susceptible to the
radial orbit instability. We show how a family of spherical models with
anisotropy profiles that transition from isotropic at the centre to radially
anisotropic becomes unstable when the outer anisotropy is made sufficiently
radial. Models with a larger outer anisotropy can be constructed but are found
to be triaxial. We argue that the onset of the radial orbit instability can be
identified with the transition point when adiabatic relaxation yields strongly
triaxial rather than weakly spherical endpoints.
We used a new generation of asymptotic giant branch (AGB) stellar models that include dust formation in the stellar winds to find the links between evolutionary models and the observed properties of a homogeneous sample of Large Magellanic Cloud (LMC) planetary nebulae (PNe). Comparison between the evolutionary yields of elements such as CNO and the corresponding observed chemical abundances is a powerful tool to shed light on evolutionary processes such as hot bottom burning (HBB) and third dredge-up (TDU). We found that the occurrence of HBB is needed to interpret the nitrogen-enriched (log(N/H)+12>8) PNe. In particular, N-rich PNe with the lowest carbon content are nicely reproduced by AGB models of mass M >=6 Mo, whose surface chemistry reflects the pure effects of HBB. PNe with log(N/H)+12<7.5 correspond to ejecta of stars that have not experienced HBB, with initial mass below about 3 Mo. Some of these stars show very large carbon abundances, owing to the many TDU episodes experienced. We found from our LMC PN sample that there is a threshold to the amount of carbon accumulated at AGB surfaces, log(C/H)+12<9. Confirmation of this constraint would indicate that, after the C-star stage is reached,AGBs experience only a few thermal pulses, which suggests a rapid loss of the external mantle, probably owing to the effects of radiation pressure on carbonaceous dust particles present in the circumstellar envelope. The implications of these findings for AGB evolution theories and the need to extend the PN sample currently available are discussed.
Microwave Kinetic Inductance Detectors (MKIDs) are the most attractive radiation detectors for far-infrared and sub-mm astronomy: They combine ultimate sensitivity with the possibility to create very large detector arrays, in excess of 10 000 pixels. This is possible by reading-out the arrays using RF frequency division multiplexing, which allows multiplexing ratios in excess of 1000 pixels per readout line. We describe a novel readout system for large arrays of MKIDs, operating in a 2 GHz band in the 4-8 GHz range. The readout, which is a combination of a digital front- and back-end and an analog up- and down-converter system, can read out up to 4000 detectors simultaneously with 1 kHz datarate. The system achieves a readout noise power spectral density of -98 dBc/Hz while reading 1000 carriers simultaneously, which scales linear with the number of carriers. We demonstrate that 4000 state-of-the-art Aluminium-NbTiN MKIDs can be read out without deteriorating their intrinsic performance.
We use images acquired with the Hubble Space Telescope Wide Field Camera 3 and new models to probe the Horizontal Branch (HB) population of the We use images acquired with the Hubble Space Telescope Wide Field Camera 3 and new models to probe the horizontal branch (HB) population of the Galactic globular cluster (GC) NGC 2419. A detailed analysis of the composite HB highlights three populations:(1) the blue luminous HB, hosting standard helium stars (Y=0.25) with a very small spread of mass, (2) a small population of stars with intermediate helium content (0.26<Y<=0.29), and (3) the well-populated extreme HB. We can fit the last group with models having high helium abundance (Y \sim 0.36), half of which (the hottest part, 'blue hook' stars) are identified as possible 'late flash mixed stars'. The initial helium abundance of this extreme population is in nice agreement with the predicted helium abundance in the ejecta of massive asymptotic giant branch (AGB) stars of the same metallicity as NGC 2419. This result further supports the hypothesis that second-generation stars in GCs formed from the ashes of intermediate-mass AGB stars. We find that the distribution in magnitude of the blue hook stars is larger than that predicted by theoretical models. We discuss the possible uncertainties in the magnitude scales and different attempts to model this group of stars. Finally, we suggest that consistency can be better achieved if we assume core masses larger than predicted by our models. This may be possible if the progenitors were fast rotators on the main sequence. If further study confirms this interpretation, a fast initial rotation would be a strong signature of the peculiarity of extreme second-generation stars in GCs.
We present a detailed study based on infrared photometry of all Galactic RV Tauri stars from the General Catalogue of Variable Stars (GCVS). RV Tauri stars are the brightest among the population II Cepheids. They are thought to evolve away from the asymptotic giant branch (AGB) towards the white dwarf domain. IRAS detected several RV Tauri stars because of their large IR excesses and it was found that they occupy a specific region in the [12] - [25], [25] - [60] IRAS two-colour diagram. We used the all sky survey of WISE to extend these studies and compare the infrared properties of all RV Tauri stars in the GCVS with a selected sample of post-AGB objects with the goal to place the RV Tauri pulsators in the context of post-AGB evolution. Moreover, we correlated the IR properties of both the RV Tauri stars and the comparison sample with other observables like binarity and the presence of a photospheric chemical anomaly called depletion. We find that Galactic RV Tauri stars display a range of infrared properties and we differentiate between disc sources, objects with no IR excess and objects for which the spectral energy distribution (SED) is uncertain. We obtain a clear correlation between disc sources and binarity. RV Tauri stars with a variable mean magnitude are exclusively found among the disc sources. We also find evidence for disc evolution among the binaries. Furthermore our studies show that the presence of a disc seems to be a necessary but not sufficient condition for the depletion process to become efficient.
We investigate the influence of different analytical parameterizations and fit functions for the local star formation rate in AMR simulations of an isolated disk galaxy with the Nyx code. Such parameterizations express the star formation efficiency as function of the local turbulent Mach number and viral parameter. By employing the method of adaptively refined large eddy simulations, we are able to evaluate these physical parameters from the numerically unresolved turbulent energy associated with the grid scale. We consider both single and multi free-fall variants of star formation laws proposed by Padoan & Nordlund, Hennebelle & Chabrier, and Krumholz & McKee. We find that the global star formation rate and the relation between the local star formation rate and the gas column density is reproduced in agreement with observational constraints by all multi free-fall models of star formation. Some models with obsolete calibration or a single free-fall time scale, however, result in an overly clumpy disk that does not even remotely resemble the structure of observed spirals.
We show how the range of application of the principal component analysis-based inversion method of Paletou et al. (2015) can be extended to late-type stars data. Besides being an extension of its original application domain, for FGK stars, we also used synthetic spectra for our learning database. We discuss our results on effective temperatures against previous evaluations made available from Vizier and Simbad services at CDS.
C$_2$H is a representative hydrocarbon that is abundant and ubiquitous in the interstellar medium (ISM). To study its chemical properties, we present Submillimeter Array (SMA) observations of the C$_2$H $N=3-2$ and HC$_3$N $J=30-29$ transitions and the 1.1 mm continuum emission toward four OB cluster-forming regions, AFGL 490, ON 1, W33 Main, and G10.6-0.4, which cover a bolometric luminosity range of $\sim$10$^3$--10$^6$ $L_{\odot}$. We found that on large scales, the C$_2$H emission traces the dense molecular envelope. However, for all observed sources, the peaks of C$_2$H emission are offset by several times times 10$^4$ AU from the peaks of 1.1 mm continuum emission, where the most luminous stars are located. By comparing the distribution and profiles of C$_2$H hyperfine lines and the 1.1 mm continuum emission, we find that the C$_2$H column density (and abundance) around the 1.1 mm continuum peaks is lower than those in the ambient gas envelope. Chemical models suggest that C$_2$H might be transformed to other species owing to increased temperature and density; thus, its reduced abundance could be the signpost of the heated molecular gas in the $\sim$10$^4$ AU vicinity around the embedded high-mass stars. Our results support such theoretical prediction for centrally embedded $\sim10^3$--$10^6L_{\odot}$ OB star-forming cores, while future higher-resolution observations are required to examine the C$_2$H transformation around the localized sites of high-mass star formation.
We apply our recently suggested statistical approach to thermal evolution of isolated neutron stars and accreting quasistationary neutron stars in X-ray transients for constraining the position and relative broadening alpha of the direct Urca threshold of powerful neutrino emission in neutron star cores. We show that most likely explanation of observations corresponds to alpha = 0.08 - 0.10 and to the neutron star mass, at which the direct Urca process is open, M_D = (1.6 - 1.8) M_sun.
Non-local thermodynamic equilibrium (NLTE) effects in diagnostically important solar Fe I lines are important due to the strong sensitivity of Fe I to ionizing UV radiation, which may lead to a considerable under-population of the Fe I levels in the solar atmosphere and, therefore, to a sizeable weakening of Fe I lines. Such NLTE effects may be intensified or weakened by horizontal radiative transfer (RT) in a three-dimensionally (3-D) structured atmosphere. We analyze the influence of horizontal RT on commonly used Fe I lines in a snapshot of a 3-D radiation magneto-hydrodynamic (MHD) simulation of a plage region. NLTE- and horizontal RT effects occur with considerable strength (up to 50% in line depth or equivalent width) in the analyzed snapshot. As they may have either sign and both signs occur with approximately the same frequency and strength, the net effects are small when considering spatially averaged quantities. The situation in the plage atmosphere turns out to be rather complex. Horizontal transfer leads to line-weakening relative to 1-D NLTE transfer near the boundaries of kG magnetic elements. Around the centers of these elements, however, we find an often significant line-strengthening. This behavior is in contrast to that expected from previous 3-D RT computations in idealized flux-tube models, which display only a line weakening. The origin of this unexpected behavior lies in the fact that magnetic elements are surrounded by dense and relatively cool down-flowing gas, which forms the walls of the magnetic elements. The continuum in these dense walls is often formed in colder gas than in the central part of the magnetic elements. Consequently, the central parts of the magnetic element experience a sub-average UV-irradiation leading to the observed 3-D NLTE line strengthening.
We gauge the impact of spacecraft-induced effects on the inferred variability properties of the light curve of the Seyfert 1 AGN Zw 229-15 observed by \Kepler. We compare the light curve of Zw 229-15 obtained from the Kepler MAST database with a re-processed light curve constructed from raw pixel data (Williams & Carini, 2015). We use the first-order structure function, $SF(\delta t)$, to fit both light curves to the damped power-law PSD of Kasliwal, Vogeley & Richards, 2015. On short timescales, we find a steeper log-PSD slope ($\gamma = 2.90$ to within $10$ percent) for the re-processed light curve as compared to the light curve found on MAST ($\gamma = 2.65$ to within $10$ percent)---both inconsistent with a damped random walk which requires $\gamma = 2$. The log-PSD slope inferred for the re-processed light curve is consistent with previous results (Carini & Ryle, 2012, Williams & Carini, 2015) that study the same re-processed light curve. The turnover timescale is almost identical for both light curves ($27.1$ and $27.5$~d for the reprocessed and MAST database light curves). Based on the obvious visual difference between the two versions of the light curve and on the PSD model fits, we conclude that there remain significant levels of spacecraft-induced effects in the standard pipeline reduction of the Kepler data. Re-processing the light curves will change the model inferenced from the data but is unlikely to change the overall scientific conclusion reached by Kasliwal et al. 2015---not all AGN light curves are consistent with the DRW.
On 2014 Dec. 9.61, the All-Sky Automated Survey for SuperNovae (ASAS-SN or "Assassin") discovered ASASSN-14lp just $\sim2$ days after first light using a global array of 14-cm diameter telescopes. ASASSN-14lp went on to become a bright supernova ($V = 11.94$ mag), second only to SN 2014J for the year. We present prediscovery photometry (with a detection less than a day after first light) and ultraviolet through near-infrared photometric and spectroscopic data covering the rise and fall of ASASSN-14lp for more than 100 days. We find that ASASSN-14lp had a broad light curve ($\Delta m_{15}(B) = 0.796 \pm 0.001_{\textrm{stat}}$), a $B$-band maximum at $2457015.823 \pm 0.030_{\textrm{stat}}$, a rise time of $16.94^{+ 0.11 }_{- 0.11 }$ days, and moderate host--galaxy extinction ($E(B-V)_{\textrm{host}} = 0.329 \pm 0.001_{\textrm{stat}}$). Using ASASSN-14lp we derive a distance modulus for NGC 4666 of $\mu = 30.834 \pm 0.003_{\textrm{stat}} \pm 0.16_{\textrm{syst}}$ corresponding to a distance of $14.68 \pm 0.02_{\textrm{stat}} \pm 1.15_{\textrm{syst}}$ Mpc. However, a tip of the red giant branch distance to the host galaxy should be measured to allow ASASSN-14lp to be added to the calibrating sample of Type Ia supernovae. Finally, using our early-time photometric and spectroscopic data along with our derived light curve properties, we rule out red giant secondaries with limits on the radius of a non-degenerate companion as small as $0.34 \rm{R}_\odot$ for favorable viewing angles and estimates of the explosion time.
The space mission PLATO will usher in a new era of exoplanetary science by expanding our current inventory of transiting systems and constraining host star ages, which are currently highly uncertain. This capability might allow PLATO to detect changes in planetary system architecture with time, particularly because planetary scattering due to Lagrange instability may be triggered long after the system was formed. Here, we utilize previously published instability timescale prescriptions to determine PLATO's capability to detect a trend of decreasing planet frequency with age for systems with equal-mass planets. For two-planet systems, our results demonstrate that PLATO may detect a trend for planet masses which are at least as massive as super-Earths. For systems with three or more planets, we link their initial compactness to potentially detectable frequency trends in order to aid future investigations when these populations will be better characterized.
4U 1538-52, an absorbed high mass X-ray binary with an orbital period of 3.73 days, shows moderate orbital intensity modulations with a low level of counts during the eclipse. Several models have been proposed to explain the accretion at different orbital phases by a spherically symmetric stellar wind from the companion. The aim of this work is to study both the light curve and orbital phase spectroscopy of this source in the long term. Particularly, the folded light curve and the changes of the spectral parameters with orbital phase to analyse the stellar wind of QV Nor, the mass donor of this binary system. We used all the observations made from the Gas Slit Camera on board MAXI of 4U 1538-52 covering many orbits continuously. We obtained the good interval times for every orbital phase range which were the input to extract our data. We estimated the orbital period of the system and then folded the light curves and we fitted the X-ray spectra with the same model for every orbital phase spectrum. We also extracted the averaged spectrum of all the MAXI data available. The MAXI spectra in the 2-20 keV energy range were fitted with an absorbed Comptonization of cool photons on hot electrons. We found a strong orbital dependence of the absorption column density but neither the fluorescence iron emission line nor low energy excess were needed to fit the MAXI spectra. The variation of the spectral parameters over the binary orbit were used to examine the mode of accretion onto the neutron star in 4U 1538-52. We deduce a best value of $\dot{M}/v_\infty=0.65\times 10^{-9}$ $M_{\odot} \, yr^{-1}/(km \, s^{-1})$ for QV Nor.
Radio-loud active galactic nuclei are among the most powerful objects in the
universe. In these objects, most of the emission comes from relativistic jets
getting their power from the accretion of matter onto supermassive black holes.
However, despite the number of studies, a jet's acceleration to relativistic
speeds is still poorly understood.
It is widely known that jets contain relativistic particles that emit
radiation through several physical processes, one of them being the inverse
Compton scattering of photons coming from external sources. In the case of a
plasma composed of electrons and positrons continuously heated by the
turbulence, inverse Compton scattering can lead to relativistic bulk motions
through the Compton rocket effect. We investigate this process and compute the
resulting bulk Lorentz factor in the complex photon field of an AGN composed of
several external photon sources.
We consider various sources here: the accretion disk, the dusty torus, and
the broad line region. We take their geometry and anisotropy carefully into
account in order to numerically compute the bulk Lorentz factor of the jet at
every altitude.
The study, made for a broad range of parameters, shows interesting and
unexpected behaviors of the bulk Lorentz factor, exhibiting acceleration and
deceleration zones in the jet. We investigate the patterns of the bulk Lorentz
factor along the jet depending on the source sizes and on the observation angle
and we finally show that these patterns can induce variability in the AGN
emission with timescales going from hours to months.
We investigate a small-scale ($\approx$ 1.5 Mm along the slit), supersonic downflow of about 90 km s$^{-1}$ in the transition region above the light-bridged sunspot umbra in AR 11836. The observations were obtained with the Interface Region Spectrograph (IRIS) on 2013 September 2, from 16:40 to 17:59 UT. The downflow shows up as red-shifted "satellite" lines of the Si IV and O IV transition region lines and is remarkably steady over the observing period of nearly 80 min. The downflow is not visible in the chromospheric lines, which only show an intensity enhancement at the location of the downflow. The density inferred from the line ratio of the red-shifted satellites of the O IV lines ($N_\mathrm{e} = 10^{10.6\pm0.25} \mathrm{cm}^{-3}$) is only a factor 2 smaller than the one inferred from the main components ($N_\mathrm{e} = 10^{10.95\pm0.20} \mathrm{cm}^{-3}$). Consequently, this implies a substantial mass flux ($\approx 5 \times 10^{-7}$ g cm$^{-2}$ s$^{-1}$), which would evacuate the overlying corona on time scales of the order of 10 s. We interpret these findings as evidence of a stationary termination shock of a supersonic siphon flow in a cool loop rooted in the central umbra of the spot.
We present a derivation of the radiometer equation based on the original references and fundamental statistical concepts. We then perform numerical simulations of white noise to illustrate the radiometer equation in action. Finally, we generate 1/f and 1/f^2 noise, demonstrate that it is non-stationary, and use it to simulate the effect of gain fluctuations on radiometer performance.
We have made an estimation of the synchrotron peak frequency ($\nu_{peak}^{s}$) for six very low synchrotron peaked (VLSP) blazars. These objects were selected as VLSP candidates (with the $\nu_{peak}^{s} \leq 10^{13}$ Hz) from the archival data. We have build spectral energy distribution using quasi-simultaneous observations at the Zeiss-1000 and RATAN-600 telescopes of the Special astrophysical observatory of RAS and made an estimation of the $\nu_{peak}^{s}$. We confirmed classification as a VLSP for three sources (PKS 0446+11, [HB89] 1308+326 and 3C 345), for three other blazars we have calculated $\nu_{peak}^{s}>10^{13}$ Hz.
We address the highly debated issue of constraining the gamma-ray emission region in blazars from cross-correlation analysis using discrete correlation function between radio and gamma-ray light curves. The significance of the correlations is evaluated using two different approaches: simulating light curves and mixed source correlations. The cross-correlation analysis yielded 26 sources with significant correlations. In most of the sources, the gamma-ray peaks lead the radio with time lags in the range +20 and +690 days, whereas in sources 1633+382 and 3C 345 we find the radio emission to lead the gamma rays by -15 and -40 days, respectively. Apart from the individual source study, we stacked the correlations of all sources and also those based on sub-samples. The time lag from the stacked correlation is +80 days for the whole sample and the distance travelled by the emission region corresponds to 7 pc. We also compared the start times of activity in radio and gamma rays of the correlated flares using Bayesian block representation. This shows that most of the flares at both wavebands start at almost the same time, implying a co-spatial origin of the activity. The correlated sources show more flares and are brighter in both bands than the uncorrelated ones.
The IAU Commission 4 Working Group on Standardizing Access to Ephemerides recommends the use of the Spacecraft and Planet Kernel (SPK) format as a standard format for the position ephemerides of planets and other natural solar system bodies, and the use of the Planetary Constants Kernel (PCK) format for the orientation of these bodies. It further recommends that other supporting data be stored in a text PCK. These formats were developed for use by the SPICE Toolkit by the Navigation and Ancillary Information Facility of NASA's Jet Propulsion Laboratory (JPL). The CALCEPH library developed by the Institut de mecanique celeste de calcul des ephemerides (IMCCE) is also able to make use of these files. High accuracy ephemerides available in files conforming to the SPK and PCK formats include: the Development Ephemerides (DE) from JPL, Integrateur Numerique Planetaire de l'Observatoire de Paris (INPOP) from IMCCE, and the Ephemerides Planets and the Moon (EPM), developed by the Institute for Applied Astronomy (IAA). The bulk of this report is a description of the portion of PCK and SPK formats required for these ephemerides. New SPK and PCK data types, both called Type 20: Chebyshev (Velocity Only), have been added. Other changes to the specification are (i) a new object identification number for coordinate time ephemerides and (ii) a set of three new data types that use the TCB rather than the TDB time scale for the ephemerides, but are otherwise identical to their TDB versions.
We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected ellipticity components of galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by galaxy type. We then compare our observations to measurements taken from the Millennium Simulation, assuming different models of galaxy-halo misalignment. Using the Millennium Simulation we verify that the statistical estimator used removes contamination from cosmic shear. We also detect an additional signal in the simulation, which we interpret as the impact of intrinsic shape-shear alignments between the lenses and their large-scale structure environment. These alignments are likely to have caused some of the previous observational constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find $f_\mathrm{h}=-0.04 \pm 0.25$ for early-type galaxies, which is consistent with current models for the galaxy-halo misalignment predicting $f_\mathrm{h}\simeq 0.20$. For late-type galaxies we measure $f_\mathrm{h}=0.69_{-0.36}^{+0.37}$ from CFHTLenS. This can be compared to the simulated results which yield $f_\mathrm{h}\simeq 0.02$ for misaligned late-type models.
Fermi-LAT >30 MeV observations have increased the number of detected solar flares by almost a factor of 10 with respect to previous space observations. These sample both the impulsive and long duration phases of GOES M and X class flares. Of particular interest is the recent detections of three solar flares whose position behind the limb was confirmed by the STEREO-B spacecraft. While gamma-ray emission up to tens of MeV resulting from proton interactions has been detected before from occulted solar flares, the significance of these particular events lies in the fact that these are the first detections of >100 MeV gamma-ray emission from footpoint-occulted flares. We will present the Fermi-LAT, RHESSI and STEREO observations of these flares and discuss the various emission scenarios for these sources and implications for the particle acceleration mechanisms.
We report on CO (J = 2 - 1) mapping with the IRAM 30-m HERA receiver array of CGCG 97-079, an irregular galaxy in the merging galaxy cluster Abell 1367 (z = 0.022). We find that $\sim$ 80% of the detected CO (J = 2 - 1) is projected within a 16 arcsec$^{2}$ (6.5 kpc$^{2}$) region to the north and west of the optical/NIR centre, with the intensity maximum offset $\sim 10$ arcsec (4 kpc) NW of the optical/NIR centre and $\sim$ 7 arcsec (3 kpc) south-east of the HI intensity maximum. Evolutionary synthesis models indicate CGCG 97-079 experienced a burst of star formation $\sim$ 10$^8$ yr ago, most likely triggered by a tidal interaction with CGCG 97-073. For CGCG 97-079 we deduce an infall velocity to the cluster of $\sim$ 1000 km s$^{-1}$ and moderate ram pressure (P$_\mathrm{ram} \sim 10^{-11}$ dyn cm$^{-2}$). The observed offset in CGCG 97-079 of the highest density HI and CO (J = 2 - 1) from the stellar components has not previously been observed in galaxies currently undergoing ram pressure stripping, although previous detailed studies of gas morphology and kinematics during ram pressure stripping were restricted to significantly more massive galaxies with deeper gravitational potential wells. We conclude the observed cold gas density maxima offsets are most likely the result of ram pressure and/or the high-speed tidal interaction with CGCG 97-073. However ram pressure stripping is likely to be playing a major role in the perturbation of lower density gas.
We present a time-variability study of young stellar objects in the cluster IRAS 20050+2720, performed at 3.6 and 4.5 micron with the Spitzer Space Telescope; this study is part of the Young Stellar Object VARiability project (YSOVAR). We have collected light curves for 181 cluster members over 40 days. We find a high variability fraction among embedded cluster members of ca. 70%, whereas young stars without a detectable disk display variability less often (in ca. 50% of the cases) and with lower amplitudes. We detect periodic variability for 33 sources with periods primarily in the range of 2-6 days. Practically all embedded periodic sources display additional variability on top of their periodicity. Furthermore, we analyze the slopes of the tracks that our sources span in the color-magnitude diagram (CMD). We find that sources with long variability time scales tend to display CMD slopes that are at least partially influenced by accretion processes, while sources with short variability time scales tend to display extinction-dominated slopes. We find a tentative trend of X-ray detected cluster members to vary on longer time scales than the X-ray undetected members.
Accretion flows around black holes generally result in mass-outflows that exhibit irregular behavior quite often. Using 2D time-dependent hydrodynamical calculations, we show that the mass-outflow is unstable in the cases of thick accretion flows such as the low angular momentum accretion flow and the advection-dominated accretion flow. For the low angular momentum flow, the inward accreting matter on the equatorial plane interacts with the outflowing gas along the rotational axis and the centrifugally supported oblique shock is formed at the interface of both the flows, when the viscosity parameter $\alpha$ is as small as $\alpha \le 10^{-3}$. The hot and rarefied blobs, which result in the eruptive mass-outflow, are generated in the inner shocked region and grow up toward the outer boundary. The advection-dominated accretion flow attains finally in the form of a torus disc with the inner edge of the disc at $3R_{\rm g} \le r \le 6R_{\rm g}$ and the center at $6R_ {\rm g} \le r \le 10R_{\rm g}$, and a series of hot blobs is intermittently formed near the inner edge of the torus and grows up along the outer surface of the torus. As a result, the luminosity and the mass-outflow rate are modulated irregularly where the luminosity is enhanced by 10-40% and the mass-outflow rate is increaed by a factor of few up to ten. We interpret the unstable nature of the outflow to be due to the Kelvin-Helmholtz instability, examining the Richardson number for the Kelvin-Helmholtz criterion in the inner region of the flow. We propose that the flare phenomena of Sgr A* may be induced by the unstable mass-outflow as is found in this work.
One of the puzzles associated with tidal disruption event candidates (TDEs) is that there is a dichotomy between the color temperatures of ${\rm few}\times 10^4$ K for TDEs discovered with optical and UV telescopes, and the color temperatures of ${\rm few}\times 10^5 - 10^6$ K for TDEs discovered with X-ray satellites. Here we propose that high-temperature TDEs are produced when the tidal debris of a disrupted star self-intersects relatively close to the SMBH, in contrast to the more distant self-intersection that leads to lower color temperatures. In particular, we note from simple ballistic considerations that greater apsidal precession in an orbit is the key to closer self-intersection. Thus larger values of $\beta$, the ratio of the tidal radius to the pericenter distance of the initial orbit, are more likely to lead to high temperatures. For a given star and $\beta$, apsidal precession also increases for larger black hole masses, but larger black hole masses imply a lower temperature at a fixed Eddington ratio. Thus the expected dependence of the temperature on the mass of the black hole is non-monotonic. We find that in order to produce a soft X-ray temperature TDE, a deeply plunging stellar orbit with $\beta> 3$ is needed and a black hole mass of $\lesssim 5\times 10^6 M_\odot$ is favored. Although observations of TDEs are comparatively scarce and are likely dominated by selection effects, it is encouraging that both predictions are consistent with current data.
Non-thermal pressure in galaxy clusters leads to underestimation of the mass of galaxy clusters based on hydrostatic equilibrium with thermal gas pressure. This occurs even for dynamically relaxed clusters that are used for calibrating the mass-observable scaling relations. We show that the analytical model for non-thermal pressure developed in Shi & Komatsu 2014 can correct for this so-called `hydrostatic mass bias', if most of the non-thermal pressure comes from bulk and turbulent motions of gas in the intracluster medium. Our correction works for the sample average irrespective of the mass estimation method, or the dynamical state of the clusters. This makes it possible to correct for the bias in the hydrostatic mass estimates from X-ray surface brightness and the Sunyaev-Zel'dovich observations that will be available for clusters in a wide range of redshifts and dynamical states.
Recent progress in the understanding of how externally driven magnetic reconnection evolves is organized in terms of parameter space diagrams. These diagrams are constructed using four pivotal dimensionless parameters: the Lundquist number $S$, the magnetic Prandtl number $P_m$, the amplitude of the boundary perturbation $\hat \Psi_0$, and the perturbation wave number $\hat k$. This new representation highlights the parameters regions of a given system in which the magnetic reconnection process is expected to be distinguished by a specific evolution. Contrary to previously proposed phase diagrams, the diagrams introduced here take into account the dynamical evolution of the reconnection process and are able to predict slow or fast reconnection regimes for the same values of $S$ and $P_m$, depending on the parameters that characterize the external drive, never considered so far. These features are important to understand the onset and evolution of magnetic reconnection in diverse physical systems
Science is an inherently quantitative endeavor, and general education science courses are taken by a majority of college students. As such, they are a powerful venue for advancing students' skills and attitudes toward mathematics. This article reports on the development and validation of the Quantitative Reasoning for College Science (QuaRCS) Assessment, a numeracy assessment instrument designed for college-level general education science students. It has been administered to more than four thousand students over eight semesters of refinement. We show that the QuaRCS is able to distinguish varying levels of quantitative literacy and present performance statistics for both individual items and the instrument as a whole. Responses from a survey of forty-eight Astronomy and Mathematics educators show that these two groups share views regarding which quantitative skills are most important in the contexts of science literacy and educated citizenship, and the skills assessed with the QuaRCS are drawn from these rankings. The fully-developed QuaRCS assessment was administered to nearly two thousand students in nineteen general education science courses and one STEM major course in early 2015, and results reveal that the instrument is valid for both populations.
The ANTARES experiment consists of an array of photomultipliers distributed along 12 lines and located deep underwater in the Mediterranean Sea. It searches for astrophysical neutrinos collecting the Cherenkov light induced by the charged particles, mainly muons, produced in neutrino interactions around the detector. Since at energies of $\sim$10 TeV the muon and the incident neutrino are almost collinear, it is possible to use the ANTARES detector as a neutrino telescope and identify a source of neutrinos in the sky starting from a precise reconstruction of the muon trajectory. To get this result, the arrival times of the Cherenkov photons must be accurately measured. A to perform time calibrations with the precision required to have optimal performances of the instrument is described. The reconstructed tracks of the atmospheric muons in the ANTARES detector are used to determine the relative time offsets between photomultipliers. Currently, this method is used to obtain the time calibration constants for photomultipliers on different lines at a precision level of 0.5 ns. It has also been validated for calibrating photomultipliers on the same line, using a system of LEDs and laser light devices.
We have measured electron impact ionization (EII) for Fe 7+ from the ionization threshold up to 1200 eV. The measurements were performed using the TSR heavy ion storage ring. The ions were stored long enough prior to measurement to remove most metastables, resulting in a beam of 94% ground state ions. Comparing with the previously recommended atomic data, we find that the Arnaud & Raymond (1992) cross section is up to about 40\% larger than our measurement, with the largest discrepancies below about 400~eV. The cross section of Dere (2007) agrees to within 10%, which is about the magnitude of the experimental uncertainties. The remaining discrepancies between measurement and the most recent theory are likely due to shortcomings in the theoretical treatment of the excitation-autoionization contribution.
We analytically calculate the influence of a plasma on the shadow of a black hole (or of another compact object). We restrict to spherically symmetric and static situations, where the shadow is circular. The plasma is assumed to be non-magnetized and pressure-less. We derive the general formulas for a spherically symmetric plasma density on an unspecified spherically symmetric and static spacetime. The formalism applies not only to black holes but also, e.g., to wormholes. As examples for the underlying spacetime model, we consider the Schwarzschild spacetime and the Ellis wormhole. In particular, we treat the case that the plasma is in radial free fall from infinity onto a Schwarzschild black hole. The perspectives of actually observing the influence of a plasma on the shadows of supermassive black holes are discussed.
The so called $f(X)$ hybrid metric-Palatini gravity presents a unique viable generalisation of the $f(R)$ theories within the metric-affine formalism. Here the cosmology of the $f(X)$ theories is studied using the dynamical system approach. The method consists of formulating the propagation equation in terms of suitable (expansion-normalised) variables as an autonomous system. The fixed points of the system then represent exact cosmological solutions described by power-law or de Sitter expansion. The formalism is applied to two classes of $f(X)$ models, revealing both standard cosmological fixed points and new accelerating solutions that can be attractors in the phase space. In addition, the fixed point with vanishing expansion rate are considered with special care in order to characterise the stability of Einstein static spaces and bouncing solutions.
The long-standing model-independent annual modulation effect measured by the DAMA Collaboration is examined in the framework of asymmetric mirror dark matter interacting with target nuclei in the detector via the kinetic mixing between mirror and ordinary photons. The allowed physical ranges for the kinetic mixing parameter are obtained taking into account various existing uncertainties in nuclear and particle physics quantities as well as in the density and velocity distributions of dark matter.
Links to: arXiv, form interface, find, astro-ph, recent, 1507, contact, help (Access key information)
We describe the fall of Annama meteorite occurred in the remote Kola Peninsula (Russia) close to Finnish border on April 19, 2014 (local time). The fireball was instrumentally observed by the Finnish Fireball Network. From these observations the strewnfield was computed and two first meteorites were found only a few hundred meters from the predicted landing site on May 29th and May 30th 2014, so that the meteorite (an H4-5 chondrite) experienced only minimal terrestrial alteration. The accuracy of the observations allowed a precise geocentric radiant to be obtained, and the heliocentric orbit for the progenitor meteoroid to be calculated. Backward integrations of the orbits of selected near-Earth asteroids and the Annama meteoroid showed that they rapidly diverged so that the Annama meteorites are unlikely related to them. The only exception seems to be the recently discovered 2014UR116 that shows a plausible dynamic relationship. Instead, analysis of the heliocentric orbit of the meteoroid suggests that the delivery of Annama onto an Earth-crossing Apollo type orbit occurred via the 4:1 mean motion resonance with Jupiter or the nu6 secular resonance, dynamic mechanisms that are responsible for delivering to Earth most meteorites studied so far.
The field of astrobiology has made huge strides in understanding the habitable zones around stars (Stellar Habitable Zones) where life can begin, sustain its existence and evolve into complex forms. A few studies have extended this idea by modelling galactic-scale habitable zones (Galactic Habitable Zones) for our Milky Way and specific elliptical galaxies. However, estimating the habitability for galaxies spanning a wide range of physical properties has so far remained an outstanding issue. Here, we present a "cosmobiological" framework that allows us to sift through the entire galaxy population in the local Universe and answer the question "Which type of galaxy is most likely to host complex life in the cosmos"? Interestingly, the three key astrophysical criteria governing habitability (total mass in stars, total metal mass and ongoing star formation rate) are found to be intricately linked through the "fundamental metallicity relation" as shown by SDSS (Sloan Digital Sky Survey) observations of more than a hundred thousand galaxies in the local Universe. Using this relation we show that metal-rich, shapeless giant elliptical galaxies at least twice as massive as the Milky Way (with a tenth of its star formation rate) can potentially host ten thousand times as many habitable (earth-like) planets, making them the most probable "cradles of life" in the Universe.
We study a two-parameter extension of the cosmological standard model $\Lambda$CDM in which cold dark matter interacts with a new form of dark radiation. The two parameters correspond to the energy density in the dark radiation $\Delta N_\mathrm{eff}$ and the interaction strength between the dark matter and dark radiation fluids. The interactions give rise to a very weak "dark matter drag" which damps the growth of matter density perturbations throughout radiation domination, allowing to reconcile the tension between predictions of large scale structure from the CMB and direct measurements of $\sigma_8$. We perform a precision fit to Planck CMB data, BAO, large scale structure, and direct measurements of the expansion rate of the universe today. Our model lowers the $\chi$-squared relative to $\Lambda$CDM by about 11, corresponding to a preference for non-zero dark matter drag by more than $3 \sigma$. Particle physics models which naturally produce a dark matter drag of the required form include the recently proposed non-Abelian dark matter model in which the dark radiation corresponds to massless dark gluons.
We present a method to identify Ultra Faint Dwarf Galaxy (UFDG) candidates in the halo of the Milky Way using the future Gaia catalogue and we explore its detection limits and completeness. The method is based on the Wavelet Transform and searches for over-densities in the combined space of sky coordinates and proper motions, using kinematics in the search for the first time. We test the method with a Gaia mock catalogue that has the Gaia Universe Model Snapshot (GUMS) as a background, and use a library of around 30 000 UFDGs simulated as Plummer spheres with a single stellar population. For the UFDGs we use a wide range of structural and orbital parameters that go beyond the range spanned by real systems, where some UFDGs may remain undetected. We characterize the detection limits as function of the number of observable stars by Gaia in the UFDGs with respect to that of the background and their apparent sizes in the sky and proper motion planes. We find that the addition of proper motions in the search improves considerably the detections compared to a photometric survey at the same magnitude limit. Our experiments suggest that Gaia will be able to detect UFDGs that are similar to some of the known UFDGs even if the limit of Gaia is around 2 magnitudes brighter than that of SDSS, with the advantage of having a full-sky catalogue. We also see that Gaia could even find some UFDGs that have lower surface brightness than the SDSS limit.
Current and future large redshift surveys, as the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), will use Emission-Line Galaxies (ELG) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range $0.6 < z < 1.7$. With current data, we explore the halo-galaxy connection by measuring three clustering properties of $g$-selected ELGs as matter tracers in the redshift range $0.6 < z < 1$: (i) the redshift-space two-point correlation function using spectroscopic redshifts from the BOSS ELG sample and VIPERS; (ii) the angular two-point correlation function on the footprint of the CFHT-LS; (iii) the galaxy-galaxy lensing signal around the ELGs using the CFHTLenS. We interpret these observations by mapping them onto the latest high-resolution MultiDark Planck N-body simulation, using a novel (Sub)Halo-Abundance Matching technique that accounts for the ELG incompleteness. ELGs at $z\sim0.8$ live in halos of $(1\pm 0.5)\times10^{12}\,h^{-1}$M$_{\odot}$ and 22.5$\pm2.5$% of them are satellites belonging to a larger halo. The halo occupation distribution of ELGs indicates that we are sampling the galaxies in which stars form in the most efficient way, according to their stellar-to-halo mass ratio.
NGC 4395 is a bulgeless spiral galaxy, harboring one of the nearest known type 1 Seyfert nuclei. Although there is no consensus on the mass of its central engine, several estimates suggest it to be one of the lightest massive black holes (MBHs) known. We present the first direct dynamical measurement of the mass of this MBH from a combination of two-dimensional gas kinematic data, obtained with the adaptive optics assisted near infrared integral field spectrograph Gemini/NIFS, and high-resolution multiband photometric data from Hubble Space Telescope's Wide Field Camera 3 (HST/WFC3). We use the photometric data to model the shape and stellar mass-to-light ratio (M/L) of the nuclear star cluster. From the Gemini/NIFS observations, we derive the kinematics of warm molecular hydrogen gas as traced by emission through the H$_2$ 1--0 S(1) transition. These kinematics show a clear rotational signal, with a position angle orthogonal to NGC 4395's radio jet. Our best fitting tilted ring models of the kinematics of the molecular hydrogen gas contain a black hole with mass $M=4_{-3}^{+8}\times 10^5$ M$_\odot$ (3$\sigma$ uncertainties) embedded in a nuclear star cluster of mass $M=2 \times 10^6$ M$_\odot$. Our black hole mass measurement is in excellent agreement with the reverberation mapping mass estimate of Peterson et al. (2005), but shows some tension with other mass measurement methods based on accretion signals.
The origin of high velocity cool gas seen in galactic winds remains unknown. Following Wang (1995), we argue that rapid radiative cooling in initially hot (10^7-10^8 K) thermally-driven outflows can produce fast neutral atomic and photoionized cool gas. Outflows with hot gas mass-loading factor relative to star formation rate of beta > 0.5 cool on scales ranging from the size of the host to tens of kpc. We provide scalings for the cooling radius r_cool, density, column density, emission measure, radiative efficiency, and cool gas velocity. At r_cool, the gas produces X-ray and then UV/optical line emission at velocities of hundreds to thousands of km/s with a total power bounded from above by the energy injection rate 0.01 L_star if the flow is powered by steady-state star formation with luminosity L_star. The wind is thermally and convectively unstable at and beyond r_cool. Thermal instability can amplify density fluctuations by a factor of ~100, potentially leading to a multi-phase medium. Cooled winds can decelerate in the extended gravitational potential of galaxies and may explain the prevalence of cool gas in galactic halos. We forward a picture of winds whereby cool clouds are initially accelerated from the host by ram pressure of the hot flow, but are rapidly shredded and incorporated into the hot flow by the Kelvin-Helmholtz instability. This increases the hot wind mass loading, seeding radiative and thermal instability and cool gas rebirth. We show that if the cooled wind re-shocks as it sweeps up the circumgalactic medium that its cooling time is short, thus depositing cool gas far out into the halo. Finally, we show that conduction can dominate energy transport in low-beta hot galactic winds, leading to much flatter temperature profiles compared to the nominal expectation from adiabaticity, potentially consistent with X-ray observations of some local starbursts. (Abridged)
Cosmic voids in the large-scale structure of the Universe affect the peculiar motions of objects in their vicinity. Although these motions are difficult to observe directly, the clustering pattern of their surrounding tracers in redshift space is influenced in a unique way. This allows to investigate the interplay between densities and velocities around voids, which is solely dictated by the laws of gravity. With the help of N-body simulations and derived mock-galaxy catalogs we calculate the average density fluctuations inside and outside voids identified with a watershed algorithm in redshift space and compare the results with the expectation from general relativity and the LCDM model of cosmology. We find that simple linear-theory predictions work remarkably well in describing the dynamics of voids even on relatively small scales. Adopting a Bayesian inference framework, we determine the full posterior probability distribution of our model parameters and forecast the achievable accuracy on measurements of the growth rate of structure and the geometric distortion through the Alcock-Paczy\'nski effect. Their relative uncertainties in galaxy surveys with number densities comparable to the SDSS MAIN (CMASS) sample that probe a volume of $1h^{-3}{\rm Gpc}^3$ yield $\sigma_{f/b}/(f/b)\sim40\%$ ($60\%$) and $\sigma_{D_AH}/D_AH\sim5\%$ ($8\%$), respectively. The presented method is highly model independent; its viability lies in the underlying assumption of statistical isotropy of the Universe.
Measurements of the galaxy power spectrum contain a wealth of information about the Universe. Its optimal extraction is vital if we are to truly understand the micro-physical nature of dark matter and dark energy. In Smith & Marian (2015) we generalized the power spectrum methodology of Feldman et al. (1994) to take into account the key tenets of galaxy formation: galaxies form and reside exclusively in dark matter haloes; a given dark matter halo may host galaxies of various luminosities; galaxies inherit the large-scale bias associated with their host halo. In this paradigm we derived the optimal weighting and reconstruction scheme for maximizing the signal-to-noise on a given band power estimate. For a future all-sky flux-limited galaxy redshift survey of depth b_J ~22, we now demonstrate that the optimal weighting scheme does indeed provide improved S/N at the level of ~20% when compared to Feldman et al. (1994) and ~60% relative to Percival et al. (2003), for scales of order k~0.5 Mpc/h. Using a Fisher matrix approach, we show that the cosmological information yield is also increased relative to these alternate methods -- especially the primordial power spectrum amplitude and dark energy equation of state.
We investigate the general interaction between an eccentric planet and a coplanar debris disc of the same mass, using analytical theory and n-body simulations. Such an interaction could result from a planet-planet scattering or merging event. We show that when the planet mass is comparable to that of the disc, the former is often circularised with little change to its semimajor axis. The secular effect of such a planet can cause debris to apsidally anti-align with the planet's orbit (the opposite of what may be naively expected), leading to the counter-intuitive result that a low-mass planet may clear a larger region of debris than a higher-mass body would. The interaction generally results in a double-ringed debris disc, which is comparable to those observed in HD 107146 and HD 92945. As an example we apply our results to HD 107146, and show that the disc's morphology and surface brightness profile can be well-reproduced if the disc is interacting with an eccentric planet of comparable mass (~10-100 Earth masses). This hypothetical planet had a pre-interaction semimajor axis of 30 or 40 au (similar to its present-day value) and an eccentricity of 0.4 or 0.5 (which would since have reduced to ~0.1). Thus the planet (if it exists) presently resides near the inner edge of the disc, rather than between the two debris peaks as may otherwise be expected. Finally we show that disc self-gravity can be important in this mass regime and, whilst it would not affect these results significantly, it should be considered when probing the interaction between a debris disc and a planet.
We present an IR-monitoring survey with the $Spitzer$ Space Telescope of the
star forming region GGD 12-15. Over 1000 objects were monitored including about
350 objects within the central 5 arcminutes which is found to be especially
dense in cluster members. The monitoring took place over 38 days and is part of
the Young Stellar Object VARiability (YSOVAR) project. The region was also the
subject of a contemporaneous 67ks $Chandra$ observation. The field includes 119
previously identified pre-main sequence star candidates. X-rays are detected
from 164 objects, 90 of which are identified with cluster members. Overall, we
find that about half the objects in the central 5 arcminutes are young stellar
objects based on a combination of their spectral energy distribution, IR
variability and X-ray emission.
Most of the stars with IR excess relative to a photosphere show large
amplitude (>0.1 mag) mid-IR variability. There are 39 periodic sources, all but
one of these is found to be a cluster member. Almost half of the periodic
sources do not show IR excesses. Overall, more than 85% of the Class I, flat
spectrum, and Class II sources are found to vary. The amplitude of the
variability is larger in more embedded young stellar objects. Most of the
Class~I/II objects exhibit redder colors in a fainter state, compatible with
time-variable extinction. A few become bluer when fainter, which can be
explained with significant changes in the structure of the inner disk. A search
for changes in the IR due to X-ray events is carried out, but the low number of
flares prevented an analysis of the direct impact of X-ray flares on the IR
lightcurves. However, we find that X-ray detected Class II sources have longer
timescales for change in the mid-IR than a similar set of non-X-ray detected
Class IIs.
We have coupled a fast, parametrized star cluster evolution code to a Markov Chain Monte Carlo code to determine the distribution of probable initial conditions of observed star clusters, which may serve as a starting point for future $N$-body calculations. In this paper we validate our method by applying it to a set of star clusters which have been studied in detail numerically with $N$-body simulations and Monte Carlo methods: the Galactic globular clusters M4, 47 Tucanae, NGC 6397, M22, $\omega$ Centauri, Palomar 14 and Palomar 4, the Galactic open cluster M67, and the M31 globular cluster G1. For each cluster we derive a distribution of initial conditions that, after evolution up to the cluster's current age, evolves to the currently observed conditions. We find that there is a connection between the morphology of the distribution of initial conditions and the dynamical age of a cluster and that a degeneracy in the initial half-mass radius towards small radii is present for clusters which have undergone a core collapse during their evolution. We find that the results of our method are in agreement with $N$-body and Monte Carlo studies for the majority of clusters. We conclude that our method is able to find reliable posteriors for the determined initial mass and half-mass radius for observed star clusters, and thus forms an suitable starting point for modeling an observed cluster\rq{}s evolution.
We present a unified computational framework which can be used to describe impulsive flares on the Sun and on dMe stars. The models assume that the flare impulsive phase is caused by a beam of charged particles that is accelerated in the corona and propagates downward depositing energy and momentum along the way. This rapidly heats the lower stellar atmosphere causing it to explosively expand and dramatically brighten. Our models consist of flux tubes that extend from the sub-photosphere into the corona. We simulate how flare-accelerated charged particles propagate down one-dimensional flux tubes and heat the stellar atmosphere using the Fokker-Planck kinetic theory. Detailed radiative transfer is included so that model predictions can be directly compared with observations. The flux of flare-accelerated particles drives return currents which additionally heat the stellar atmosphere. These effects are also included in our models. We examine the impact of the flare-accelerated particle beams on model solar and dMe stellar atmospheres and perform parameter studies varying the injected particle energy spectra. We find the atmospheric response is strongly dependent on the accelerated particle cutoff energy and spectral index.
Cosmological constraints from galaxy clusters rely on accurate measurements of the mass and internal structure of clusters. An important source of systematic uncertainty in cluster mass and structure measurements is the secure selection of background galaxies that are gravitationally lensed by clusters. This issue has been shown to be particular severe for faint blue galaxies. We therefore explore the selection of faint blue background galaxies, by reference to photometric redshift catalogs derived from the COSMOS survey and our own observations of massive galaxy clusters at z~0.2. We show that methods relying on photometric redshifts of galaxies in/behind clusters based on observations through five filters, and on deep 30-band COSMOS photometric redshifts are both inadequate to identify safely faint blue background galaxies. This is due to the small number of filters used by the former, and absence of massive galaxy clusters at redshifts of interest in the latter. We therefore develop a pragmatic method to combine both sets of photometric redshifts to select a population of blue galaxies based purely on photometric analysis. This sample yields stacked weak-lensing results consistent with our previously published results based on red galaxies. We also show that the stacked clustercentric number density profile of these faint blue galaxies is consistent with expectations from consideration of the lens magnification signal of the clusters. Indeed, the observed number density of blue background galaxies changes by ~10-30 per cent across the radial range over which other surveys assume it to be flat.
We present a comprehensive analysis of strong-lensing, weak-lensing shear and magnification data for a sample of 16 X-ray-regular and 4 high-magnification galaxy clusters selected from the CLASH survey. Our analysis combines constraints from 16-band HST observations and wide-field multi-color imaging taken primarily with Subaru/Suprime-Cam. We reconstruct surface mass density profiles of individual clusters from a joint analysis of the full lensing constraints, and determine masses and concentrations for all clusters. We find internal consistency of the ensemble mass calibration to be $\le 5\% \pm 6\%$ by comparison with the CLASH weak-lensing-only measurements of Umetsu et al. For the X-ray-selected subsample, we examine the concentration-mass relation and its intrinsic scatter using a Bayesian regression approach. Our model yields a mean concentration of $c|_{z=0.34} = 3.95 \pm 0.35$ at $M_{200c} \simeq 14\times 10^{14}M_\odot$ and an intrinsic scatter of $\sigma(\ln c_{200c}) = 0.13 \pm 0.06$, in excellent agreement with LCDM predictions when the CLASH selection function based on X-ray morphological regularity and the projection effects are taken into account. We also derive an ensemble-averaged surface mass density profile for the X-ray-selected subsample by stacking their individual profiles. The stacked mass profile is well described by a family of density profiles predicted for cuspy dark-matter-dominated halos, namely, the NFW, Einasto, and DARKexp models, whereas the single power-law, cored isothermal and Burkert density profiles are disfavored by the data. We show that cuspy halo models that include the two-halo term provide improved agreement with the data. For the NFW halo model, we measure a mean concentration of $c_{200c} = 3.79^{+0.30}_{-0.28}$ at $M_{200c} = 14.1^{+1.0}_{-1.0}\times 10^{14}M_\odot$, demonstrating robust consistency between complementary analysis methods.
Gas can be violently stripped from their galaxy disks in rich clusters, and be dispersed over 100kpc-scale tails or plumes. Young stars have been observed in these tails, suggesting they are formed in situ. This will contribute to the intracluster light, in addition to tidal stripping of old stars. We want to quantify the efficiency of intracluster star formation. We present CO(1--0) and CO(2--1) observations, made with the IRAM-30m telescope, towards the ram-pressure stripped tail northeast of NGC4388 in Virgo. HII regions found all along the tails, together with dust patches have been targeted. We detect molecular gas in 4 positions along the tail, with masses between 7x10$^5$ to 2x10$^6$ M$_\odot$. Given the large distance from the NGC 4388 galaxy, the molecular clouds must have formed in situ, from the HI gas plume. We compute the relation between surface densities of star formation and molecular gas in these regions, and find that the star formation has very low efficiency. The corresponding depletion time of the molecular gas can be up to 500 Gyr and more. Since this value exceeds a by far Hubble time, this gas will not be converted into stars, and will stay in a gaseous phase to join the intracluster medium.
Intensity oscillations of coronal bright points (BPs) have been studied for past several years. It has been known for a while that these BPs are closed magnetic loop like structures. However, initiation of such intensity oscillations is still an enigma. There have been many suggestions to explain these oscillations, but modeling of such BPs have not been explored so far. Using a multithreaded nanoflare heated loop model we study the behavior of such BPs in this work. We compute typical loop lengths of BPs using potential field line extrapolation of available data (Chandrashekhar et al. 2013), and set this as the length of our simulated loops. We produce intensity like observables through forward modeling and analyze the intensity time series using wavelet analysis, as was done by previous observers. The result reveals similar intensity oscillation periods reported in past observations. It is suggested these oscillations are actually shock wave propagations along the loop. We also show that if one considers different background subtractions, one can extract adiabatic standing modes from the intensity time series data as well, both from the observed and simulated data.
The SKA will be transformational for many areas of science, but in particular for the study of neutron stars and their usage as tools for fundamental physics in the form of radio pulsars. Since the last science case for the SKA, numerous and unexpected advances have been made broadening the science goals even further. With the design of SKA Phase 1 being finalised, it is time to confront the new knowledge in this field, with the prospects promised by this exciting new telescope. While technically challenging, we can build our expectations on recent discoveries and technical developments that have reinforced our previous science goals.
Rotational shear layers at the boundary between radiative and convective zones, tachoclines, play a key role in the dynamo process of magnetic field generation in the Sun and solar-like stars. We present two sets of global simulations of rotating turbulent convection and dynamo. The first set considers a stellar convective envelope only; the second one, aiming at the formation of a tachocline, considers also the upper part of the radiative zone. Our results indicate that dynamo properties like the growth rate, the saturation energy and mode depend on the Rossby (Ro) number. The models with slow rotation in the first set of simulations reproduce remarkably well the solar differential rotation in the convection zone. Depending on the value of Ro either oscillatory (with ~2 yr period) or steady dynamo solutions are obtained. The models in the second set naturally develop a tachocline which, in turn, leads to the generation of strong mean magnetic field (~1 Tesla). Since the field is also deposited into the stable deeper layer, its evolutionary time-scale is much longer than in the models without a tachocline. Surprisingly, the magnetic field in the upper turbulent convection zone evolves in the same time scale as the deep field. These models result in either an oscillatory dynamo with ~30 yr period or in a steady dynamo depending on the value of Ro. Finally, we study the evolution of the magnetic field in terms of the mean-field dynamo source terms with the dynamo coefficients computed using the FOSA approximation. For the oscillatory models without a tachocline the evolution of the fields seems to be consistent with a pattern of dynamo waves propagating according to the Parker-Yoshimura sign rule. In the models with tachocline the dynamics is more complex and could involve other transport mechanisms due to the longer time-scale.
NGC 1851 is an intriguing Galactic globular cluster, with multiple stellar evolutionary sequences, light and heavy element abundance variations and indications of a surrounding stellar halo. We present the first results of a spectroscopic study of red giant stars within and outside of the tidal radius of this cluster. Our results identify nine probable new cluster members (inside the tidal radius) with heliocentric radial velocities consistent with that of NGC 1851. We also identify, based on their radial velocities, four probable extratidal cluster halo stars at distances up to ~3.1 times the tidal radius, which are supportive of previous findings that NGC 1851 is surrounded by an extended stellar halo. Proper motions were available for 12 of these 13 stars and all are consistent with that of NGC 1851. Apart from the cluster members and cluster halo stars, our observed radial velocity distribution agrees with the expected distribution from a Besancon disk/N-body stellar halo Milky Way model generated by the Galaxia code, suggesting that no other structures at different radial velocities are present in our field. The metallicities of these stars are estimated using equivalent width measurements of the near infrared calcium triplet absorption lines and are found, within the limitations of this method, to be consistent with that of NGC 1851. In addition we recover 110 red giant cluster members from previous studies based on their radial velocities and identify three stars with unusually high radial velocities.
We report a discovery of 6 massive galaxies with both extremely large Lya equivalent width and evolved stellar population at z ~ 3. These MAssive Extremely STrong Lya emitting Objects (MAESTLOs) have been discovered in our large-volume systematic survey for strong Lya emitters (LAEs) with twelve optical intermediate-band data taken with Subaru/Suprime-Cam in the COSMOS field. Based on the SED fitting analysis for these LAEs, it is found that these MAESTLOs have (1) large rest-frame equivalent width of EW_0(Lya) ~ 100--300 A, (2) M_star ~ 10^10.5--10^11.1 M_sun, and (3) relatively low specific star formation rates of SFR/M_star ~ 0.03--1 Gyr^-1. Three of the 6 MAESTLOs have extended Ly$\alpha$ emission with a radius of several kpc although they show very compact morphology in the HST/ACS images, which correspond to the rest-frame UV continuum. Since the MAESTLOs do not show any evidence for AGNs, the observed extended Lya emission is likely to be caused by star formation process including the superwind activity. We suggest that this new class of LAEs, MAESTLOs, provides a missing link from star-forming to passively evolving galaxies at the peak era of the cosmic star-formation history.
The modification of star formation (SF) in galaxy interactions is a complex process, with SF observed to be both enhanced in major mergers and suppressed in minor pair interactions. Such changes likely to arise on short timescales and be directly related to the galaxy-galaxy interaction time. Here we investigate the link between dynamical phase and direct measures of SF on different timescales for pair galaxies, targeting numerous star-formation rate (SFR) indicators and comparing to pair separation, individual galaxy mass and pair mass ratio. We split our sample into the higher (primary) and lower (secondary) mass galaxies in each pair and find that SF is indeed enhanced in all primary galaxies but suppressed in secondaries of minor mergers. We find that changes in SF of primaries is consistent in both major and minor mergers, suggesting that SF in the more massive galaxy is agnostic to pair mass ratio. We also find that SF is enhanced/suppressed more strongly for short-time duration SFR indicators (e.g. H-alpha), highlighting recent changes to SF in these galaxies, which are likely to be induced by the interaction. We propose a scenario where the lower mass galaxy has its SF suppressed by gas heating or stripping, while the higher mass galaxy has its SF enhanced, potentially by tidal gas turbulence and shocks. This is consistent with the seemingly contradictory observations for both SF suppression and enhancement in close pairs.
In a search with the Parkes radio telescope of 56 unidentified Fermi-LAT gamma-ray sources, we have detected 11 millisecond pulsars (MSPs), 10 of them discoveries, of which five were reported in Kerr et al. (2012). We did not detect radio pulsations from another six pulsars now known in these sources. We describe the completed survey, which included multiple observations of many targets done to minimize the impact of interstellar scintillation, acceleration effects in binary systems, and eclipses. We consider that 23 of the 39 remaining sources may still be viable pulsar candidates. We present timing solutions and polarimetry for five of the MSPs, and gamma-ray pulsations for PSR J1903-7051 (pulsations for five others were reported in the second Fermi-LAT catalog of gamma-ray pulsars). Two of the new MSPs are isolated and five are in >1 d circular orbits with 0.2-0.3 Msun presumed white dwarf companions. PSR J0955-6150, in a 24 d orbit with a ~0.25 Msun companion but eccentricity of 0.11, belongs to a recently identified class of eccentric MSPs. PSR J1036-8317 is in an 8 hr binary with a >0.14 Msun companion that is probably a white dwarf. PSR J1946-5403 is in a 3 hr orbit with a >0.02 Msun companion with no evidence of radio eclipses.
Context. MWC 656 is the recently discovered first binary system case composed
of a Be-type star and an accreting black hole. Its low X-ray luminosity
indicates that the system is in a quiescent X-ray state.
Aims. The aim of our investigation is to establish if the MWC 656 system has
detectable radio emission and if the radio characteristics are consistent with
those of quiescent black hole systems.
Methods.We used three archived VLA data sets, one hour each, at 3 GHz and
seven new VLA observations, two hours each, at 10 GHz to produce very high
sensitivity images, down to $\sim$1$\,\mu$Jy.
Results.We detected the source twice in the new observations: in the first
VLA run, at periastron passage, with a flux density of 14.2$\,\pm\,$2.9 $\mu$Jy
and by combining all together the other six VLA runs, with a flux density of
$3.7 \pm 1.4$ $\mu$Jy. The resulting combined map of the archived observations
has the sensitivity of $1 \sigma = 6.6\, \mu Jy$ but no radio emission is there
detected.
Conclusions. The radio and X-ray luminosities agree with the behaviour of
accreting binary black holes in the hard and quiescent state. In particular,\
\mwc{} in the $L_X$, $L_R$ plane occupies the same region as A0620$-$00 and XTE
J1118+480, the faintest known black holes up to now.
We present weak-lensing mass measurements of 50 X-ray luminous galaxy clusters at $0.15\le z\le0.3$, based on high quality observations with Suprime-Cam mounted on the 8.2-m Subaru telescope. We pay close attention to possible systematic biases, aiming to control them at the $\lt4$ per cent level. The dominant source of systematic bias in weak-lensing measurements of the mass of individual galaxy clusters is contamination of background galaxy catalogues by faint cluster and foreground galaxies. We extend our conservative method for selecting background galaxies with $(V-i')$ colours redder than the red sequence of cluster members to use a colour-cut that depends on cluster-centric radius. This allows us to define background galaxy samples that suffer $\le1$ per cent contamination, and comprise $13$ galaxies per square arcminute. Thanks to the purity of our background galaxy catalogue, the largest systematic in our measurement is a shape measurement bias of $3$ per cent, that we measure using custom-made simulations that probe weak shears upto $g=0.3$. Our individual cluster mass and concentration measurements are in excellent agreement with predictions of the mass-concentration relation. Equally, our stacked shear profile is in excellent agreement with the Navarro Frenk and White profile. Our new LoCuSS mass measurements are consistent with the CCCP and CLASH surveys, and in tension with the Weighing the Giants (WtG) at $\sim2\sigma$ significance.
Whilst most galaxy properties scale with galaxy mass, similar scaling relations for angular momentum are harder to demonstrate. A lognormal (LN) density distribution for disc mass provides a good overall fit to the observational data for disc rotation curves for a wide variety of galaxy types and luminosities. In this paper, the total angular momentum J and energy $\vert{}$E$\vert{}$ were computed for 38 disc galaxies from the published rotation curves and plotted against the derived disc masses, with best fit slopes of 1.683$\pm{}$0.018 and 1.643$\pm{}$0.038 respectively, using a theoretical model with a LN density profile. The derived mean disc spin parameter was $\lambda{}$=0.423$\pm{}$0.014. Using the rotation curve parameters V$_{max}$ and R$_{max}$ as surrogates for the virial velocity and radius, the virial mass estimator $M_{disc}\propto{}R_{max}V_{max}^2$ was also generated, with a log-log slope of 1.024$\pm{}$0.014 for the 38 galaxies, and a proportionality constant ${\lambda{}}^*=1.47\pm{}0.20\times{}{10}^5\ M_{sun\ }{kpc}^{-1}{km}^{-2}\ s^2$. This relationship was close to the theoretical slope of 1, and had less scatter than the corresponding Tully Fisher relation, $M\propto{}{\left(V_{rot}\right)}^{\alpha{}}$, suggesting that the virial mass estimator may provide an alternative method to determine disc masses.
We present a multi-wavelength study of a newly discovered compact group (CG), SDSS J0959+1259, based data from XMM-Newton, SDSS and the Calar Alto optical imager BUSCA. With a maximum velocity offset of 500 km s$^{-1}$, a mean redshift of 0.035, and a mean spatial extension of 480 kpc, this CG is exceptional in having the highest concentration of nuclear activity in the local Universe, established with a sensitivity limit L$_{X}>4\times $10$^{40}$ erg s$^{-1}$ in 2--10 keV band and R-band magnitude $M_R < -19$. The group is composed of two type-2 Seyferts, one type-1 Seyfert, two LINERs and three star forming galaxies. Given the high X-ray luminosity of LINERs which reaches $\sim 10^{41}$ erg s$^{-1}$, it is likely that they are also accretion driven, bringing the number of active nuclei in this group to to 5 out of 8 (AGN fraction of 60\%). The distorted shape of one member of the CG suggests that strong interactions are taking place among its galaxies through tidal forces. Therefore, this system represents a case study for physical mechanisms that trigger nuclear activity and star formation in CGs.
We present spectral and timing analysis of NuSTAR observations of the accreting X-ray pulsar 2RXP J130159.6-635806. The source was serendipitously observed during a campaign focused on the gamma-ray binary PSR B1259-63 and was later targeted for a dedicated observation. The spectrum has a typical shape for accreting X-ray pulsars, consisting of a simple power law with an exponential cutoff starting at ~7 keV with a folding energy of E_fold=~18 keV. There is also an indication of the presence of a 6.4 keV iron line in the spectrum at the ~3 sigma significance level. NuSTAR measurements of the pulsation period reveal that the pulsar has undergone a strong and steady spin-up for the last 20 years. The pulsed fraction is estimated to be ~80%, and is constant with energy up to 40 keV. The power density spectrum shows a break towards higher frequencies relative to the current spin period. This, together with steady persistent luminosity, points to a long-term mass accretion rate high enough to bring the pulsar out of spin equilibrium.
Recent studies have suggested the possibility of significantly obscuring super-soft X-ray sources in relatively modest amounts of local matter lost from the binaries themselves. If correct, then this would have explained the paucity of observed super-soft X-ray sources and would have significance for the search for single-degenerate type Ia supernova progenitors. We point out that earlier studies of circumbinary obscuration ignored photo-ionisations of the gas by the emission from the super-soft X-ray source. We revisit the problem using a full, self-consistent calculation of the ionisation state of the circumbinary material photo-ionised by the radiation of the central source. Our results show that the circumstellar mass-loss rates required for obcuration of super-soft X-ray sources is about an order of magnitude larger than those reported in earlier studies, for comparable model parameters. While this does not entrirely rule out the possibility of circumstellar material obscuring super-soft X-ray sources, it makes it unlikely that this effect alone can account for the majority of the missing super-soft X-ray sources. We discuss the observational appearance of hypothetical obscured nuclear burning white dwarfs and show that they have signatures making them distinct from photo-ionised nebulae around super-soft X-ray sources imbedded in the low density ISM.
HI intensity mapping is an emerging tool to probe dark energy. Observations of the redshifted HI signal will be contaminated by instrumental noise, atmospheric and Galactic foregrounds. The latter is expected to be four orders of magnitude brighter than the HI emission we wish to detect. We present a simulation of single-dish observations including an instrumental noise model with 1/f and white noise, and sky emission with a diffuse Galactic foreground and HI emission. We consider two foreground cleaning methods: spectral parametric fitting and principal component analysis. For a smooth frequency spectrum of the foreground and instrumental effects, we find that the parametric fitting method provides residuals that are still contaminated by foreground and 1/f noise, but the principal component analysis can remove this contamination down to the thermal noise level. This method is robust for a range of different models of foreground and noise, and so constitutes a promising way to recover the HI signal from the data. However, it induces a leakage of the cosmological signal into the subtracted foreground of around 5%. The efficiency of the component separation methods depends heavily on the smoothness of the frequency spectrum of the foreground and the 1/f noise. We find that as, long as the spectral variations over the band are slow compared to the channel width, the foreground cleaning method still works.
The internal angular momentum distribution of a star is key to determine its evolution. Fortunately, the stellar internal rotation can be probed through studies of rotationally-split non-radial oscillation modes. In particular, detection of non-radial gravity modes (g modes) in massive young stars has become feasible recently thanks to the Kepler space mission. Our aim is to derive the internal rotation profile of the Kepler B8V star KIC 10526294 through asteroseismology. We interpret the observed rotational splittings of its dipole g modes using four different approaches based on the best seismic models of the star and their rotational kernels. We show that these kernels can resolve differential rotation the radiative envelope if a smooth rotational profile is assumed and the observational errors are small. Based on Kepler data, we find that the rotation rate near the core-envelope boundary is well constrained to $163\pm89$ nHz. The seismic data are consistent with rigid rotation but a profile with counter-rotation within the envelope has a statistical advantage over constant rotation. Our study should be repeated for other massive stars with a variety of stellar parameters in order to deduce the physical conditions that determine the internal rotation profile of young massive stars, with the aim to improve the input physics of their models.
The spin-down of a neutron star, e.g. due to magneto-dipole losses, results in compression of the stellar matter and induces nuclear reactions at phase transitions between different nuclear species in the crust. We show that this mechanism is effective in heating recycled pulsars, in which the previous accretion process has already been compressing the crust, so it is not in nuclear equilibrium. We calculate the corresponding emissivity and confront it with available observations, showing that it might account for the likely thermal ultraviolet emission of PSR J0437-4715.
It is shown that the unusual thermodynamic properties of matter within the region of two-phase coexistence in hybrid stars result in a change of the standard condition for beginning of convection. In particular, the thermal flux transported by convection may be directed towards the stellar center. We discuss favorable circumstances leading to such an effect of "inverse convection" and its possible influence on the thermal evolution of hybrid stars.
Molecular clouds consist typically of 3/4 H2, 1/4 He and traces of heavier
elements. In an earlier work we showed that at very low temperatures and high
densities, H2 can be in a phase transition leading to the formation of ice
clumps as large as comets, or even planets. However, He has very different
chemical properties and no phase transition is expected before H2 in dense ISM
conditions. The gravitational stability of fluid mixtures has been studied
before, but not including a phase transition.
We study the gravitational stability of binary fluid mixtures with special
emphasis if one component is in a phase transition. The results are aimed at
applications in molecular cloud conditions.
We study the gravitational stability of van der Waals fluid mixtures using
linearised analysis and examine virial equilibrium conditions using the
Lennard-Jones inter-molecular potential. Then, combining the Lennard-Jones and
gravitational potentials, the non-linear dynamics of fluid mixtures are studied
using the molecular dynamics code LAMMPS.
Besides the classical ideal-gas Jeans instability criterion, a fluid mixture
is always gravitationally unstable if it is in a phase transition. In unstable
situations the species can separate: in some conditions He precipitates faster
than H2, while in other conditions the converse occurs. Also, for an initial
gas phase collapse the geometry is essential: contrary to spherical or
filamentary collapses, sheet-like collapses starting below 15 K allow to easily
reach H2 condensation conditions because then it is the fastest, and both the
increase of heating and opacity are limited.
Depending on density, temperature and mass, either rocky H2 planetoids, or
gaseous He planetoids form. H2 planetoids are favoured by high density, low
temperature and low mass, while He planetoids need more mass and can form at
temperature well above the critical one.
Next-generation cosmological surveys will probe ever larger volumes of the Universe, including the largest scales, near and beyond the horizon. On these scales, the galaxy power spectrum carries signatures of local primordial non-Gaussianity (PNG) and horizon-scale General Relativistic (GR) effects. But cosmic variance severely limits detection of horizon-scale effects. In order to beat down cosmic variance, we can combine surveys via the multi-tracer technique. This method benefits from large bias differences between two tracers of the underlying dark matter distribution, which suggests a multi-wavelength combination of large volume surveys that are planned on a similar time-scale. We show that the combination of two contemporaneous surveys, a large neutral hydrogen intensity mapping survey in SKA Phase\,1 and a Euclid-like photometric survey, will provide unprecedented constraints on PNG as well as detection of the GR effects. We forecast that the error on local PNG will break through the cosmic-variance limit on cosmic microwave background surveys, and achieve $\sigma(f_{\rm NL})\simeq1.37-0.48$, depending on assumed priors and the final bias and source counts. Moreover it should make the first measurements of GR effects with $\sim7\%$ accuracy which are more robust to the assumed fiducial model.
Recent work by Bulbul et al. and Boyarsky et al. has suggested that a line feature at approx. 3.5 keV in the X-ray spectra of galaxy clusters and individual galaxies seen with XMM-Newton is due to the decay of sterile neutrinos, a dark matter candidate. This identification has been criticized by Jeltema and Profumo on the grounds that model spectra suggest that atomic transitions in helium-like potassium (K XVIII) and chlorine (Cl XVI) are more likely to be the emitters. Here it is pointed out that the K XVIII lines have been observed in numerous solar flare spectra at high spectral resolution with the RESIK crystal spectrometer and also appear in Chandra HETG spectra of the coronally active star sigma Gem. In addition, the solar flare spectra at least indicate a mean coronal potassium abundance which is a factor of between 9 and 11 higher than the solar photospheric abundance. This fact, together with the low statistical quality of the XMM-Newton spectra, completely accounts for the approx. 3.5 keV feature and there is therefore no need to invoke a sterile neutrino interpretation of the observed line feature at 3.5 keV.
In this article, we argue that models based on machine learning (ML) can be very effective in estimating the non-linear matter power spectrum ($P(k)$). We employ the prediction ability of the supervised ML algorithms to build an estimator for the $P(k)$. The estimator is trained on a set of cosmological models, and redshifts for which the $P(k)$ is known, and it learns to predict $P(k)$ for any other set. We review three ML algorithms -- Random Forest, Gradient Boosting Machines, and K-Nearest Neighbours -- and investigate their prime parameters to optimize the prediction accuracy of the estimator. We also compute an optimal size of the training set, which is realistic enough, and still yields high accuracy. We find that, employing the optimal values of the internal parameters, a set of $50-100$ cosmological models is enough to train the estimator that can predict the $P(k)$ for a wide range of cosmological models, and redshifts. Using this configuration, we build a blackbox -- Supervised Estimator for Matter Power Spectrum (SEMPS) -- that computes the $P(k)$ to 5-10$\%$ accuracy up to $k\sim 10 h^{-1}{\rm Mpc}$ with respect to the reference model (cosmic emulator). We also compare the estimations of SEMPS to that of the Halofit, and find that for the $k$-range where the cosmic variance is low, SEMPS estimates are better than that of the Halofit. The predictions of the SEMPS are instantaneous in the sense that it can evaluate up to 500 $P(k)$ in less than one second, which makes it ideal for many applications like visualisations, weak lensing, emulations, likelihood analysis etc.. As a supplement to this article, we provide a publicly available software package.
The solar system has changed dramatically since its birth, and so did our understanding of it. A considerable research effort has been invested in the past decade in an attempt to reconstruct the solar system history, including the earliest stages some 4.5 billion years ago. The results indicate how several processes, such as planetary migration and dynamical instabilities, acted to relax the orbital spacing of the outer planets, and provided the needed perturbation to explain the present planetary orbits that are not precisely circular and coplanar. Here we highlight this work and illustrate the key results in a computer simulation that unifies several recently developed theories. The emerging view represents another step away from the initial perception of the solar system as part of unchanging heavens.
The recent measurements of the masses of the pulsar J00737-3039B and of the companion J1756-2251 and pulsars PSR J1614-2230, PSR J0348-0432 demonstrate the existence of compact stars with masses in a broad range from 1.2 to 2 $M_\odot$. To fulfill the constraint $M_{\rm max}>2M_{\odot}$ and to demonstrate the possibility of cooling scenarios for purely hadronic and further for hybrid stars we exploit the stiff DD2 hadronic equation of state producing a maximum neutron star mass $M\simeq 2.43 M_{\odot}$. We show that the "nuclear medium cooling" scenario for neutron stars comfortably explains the whole set of cooling curves just by a variation of the star masses without the necessity for the occurrence of the direct Urca reaction. To describe the cooling data with the very stiff DD2 equation of state we select a proton gap profile from those exploited in the literature and allow for a variation of the effective pion gap controlling the efficiency of the medium modified Urca process. Fast cooling of young neutron stars like it is seen in the data for Cas A is explained with the DD2 equation of state when the following conditions are provided: the presence of an efficient medium modified Urca process, and a large proton gap at densities $n\le 2n_0$ vanishing for $n\ge (2.5 - 3) n_0$, where $n_0$ is the saturation nuclear density.
Many observed CEMP stars are found in binary systems and show enhanced abundances of $s$-elements. The origin of the chemical abundances of these CEMP-$s$ stars is believed to be accretion in the past of enriched material from a primary star in the AGB phase. We investigate the mechanism of mass transfer and the process of nucleosynthesis in low-metallicity AGB stars by modelling the binary systems in which the observed CEMP-$s$ stars were formed. For this purpose we compare a sample of $67$ CEMP-$s$ stars with a grid of binary stars generated by our binary evolution and nucleosynthesis model. We classify our sample CEMP-$s$ stars in three groups based on the observed abundance of europium. In CEMP$-s/r$ stars the europium-to-iron ratio is more than ten times higher than in the Sun, whereas it is lower than this threshold in CEMP$-s/nr$ stars. No measurement of europium is currently available for CEMP-$s/ur$ stars. On average our models reproduce well the abundances observed in CEMP-$s/nr$ stars, whereas in CEMP-$s/r$ stars and CEMP-$s/ur$ stars the abundances of the light-$s$ elements are systematically overpredicted by our models and in CEMP-$s/r$ stars the abundances of the heavy-$s$ elements are underestimated. In all stars our modelled abundances of sodium overestimate the observations. This discrepancy is reduced only in models that underestimate the abundances of most of the $s$-elements. Furthermore, the abundance of lead is underpredicted in most of our model stars. These results point to the limitations of our AGB nucleosynthesis model, particularly in the predictions of the element-to-element ratios. Finally, in our models CEMP-$s$ stars are typically formed in wide systems with periods above 10000 days, while most of the observed CEMP-$s$ stars are found in relatively close orbits with periods below 5000 days.
We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.
We investigate roles of Alfvenic waves in the weakly-ionized atmosphere of hot Jupiters by carrying out non-ideal magnetohydrodynamic (MHD) simulations with Ohmic diffusion in one-dimensional magnetic flux tubes. Turbulence at the surface excites Alfven waves and they propagate upward to drive hot (~ 10^4 K) outflows. The magnetic diffusion plays an important role in the dissipation of the Alfvenic waves in the weakly ionized atmosphere of hot Jupiters. The mass-loss rate of the spontaneously driven planetary wind is considerably reduced, in comparison with that obtained from ideal MHD simulations because the Alfvenic waves are severely damped at low altitudes in the atmosphere, whereas the wave heating is still important in the heating of the upper atmosphere. Dependence on the surface temperature, planetary radius, and velocity dispersion at the surface is also investigated. We find an inversion phenomenon of the transmitted wave energy flux; the energy flux carried by Alfven waves in the upper atmosphere has a nonmonotonic correlation with the input energy flux from the surface in a certain range of the surface temperature because the resistivity is determined by the global physical properties of the atmosphere in a complicated manner. We also point out that the heating and mass loss are expected only in limited zones if the open magnetic field is confined in the limited regions.
The inner region of the Milky Way is one of the most interesting and complex regions of the gamma-ray sky. The intense interstellar emission and resolved point sources, as well as potential contributions by other sources such as unresolved source populations and dark matter, complicate the interpretation of the data. In this paper the Fermi LAT team analysis of a 15x15 degree region about the Galactic centre is described. The methodology for point-source detection and treatment of the interstellar emission is given. In general, the bulk of the gamma-ray emission from this region is attributable to a combination of these two contributions. However, low-intensity residual emission remains and its characterisation is discussed.
We study the stability properties of Rayleigh unstable flows both in the purely hydrodynamic and magnetohydrodynamic (MHD) regimes for two different values of the shear $q=2.1, 4.2$ ($q = - d\ln\Omega / d\ln r$) and compare it with the Keplerian case $q=1.5$. The Rayleigh stability criterion states that hydrodynamic shear flows are stable for $q<2$ but with a weak magnetic field, Rayleigh stable flows are unstable to the magnetorotational instability (MRI). The shearing box approximation is not particularly suited for the $q>2$ regime as the volume averaged velocities ($k=0$ mode) are unstable in this regime but the advantage of using a pseudospectral code is that the $k=0$ mode is conserved. We find that the $q>2$ regime is unstable to turbulence both in the hydrodynamic and in the MHD limit (with an initially weak magnetic field). In the $q>2$ regime, the velocity fluctuations dominate the magnetic fluctuations whereas in the $q<2$ regime the magnetic fluctuations dominate. This highlights two different paths to MHD turbulence implied by the two regimes, suggesting that in the $q>2$ regime the instability produces primarily velocity fluctuations that cause magnetic fluctuations, with the causality reversed for the $q<2$ MRI unstable regime. We also find that the magnetic field correlation is increasingly localized as the shear is increased in the Rayleigh unstable regime, a trend not present for the MRI unstable regime.
A detailed analysis is presented of 33 RR Lyrae stars in Pisces observed with the Kepler space telescope over the 8.9-day long K2 Two-Wheel Concept Engineering Test. The sample includes not only fundamental-mode and first overtone (RRab and RRc) stars but the first two double-mode (RRd) stars that Kepler detected and the only modulated first-overtone star ever observed from space so far. The precision of the extracted K2 light curves made it possible to detect low-amplitude additional modes in all subtypes. All RRd and non-modulated RRc stars show the additional mode at P_X/P_1~0.61 that was detected in previous space-based photometric measurements. A periodicity longer than the fundamental mode was tentatively identified in one RRab star that might belong to a gravity mode. We determined the photometric [Fe/H] values for all fundamental-mode stars and provide the preliminary results of our efforts to fit the double-mode stars with non-linear hydrodynamic pulsation models. The results from this short test run indicate that the K2 mission will be, and has started to be, an ideal tool to expand our knowledge about RR Lyrae stars. As a by-product of the target search and analysis, we identified 165 bona-fide double-mode RR Lyrae stars from the Catalina Sky Survey observations throughout the sky, 130 of which are new discoveries.
The study is devoted to radial pulsations in Population I Cepheids with masses from 5.4 to 6 Msun. Solution of the equations of radiation hydrodynamics and turbulent convection for nonlinear stellar pulsations was obtained with initial conditions as the core helium burning models of computed evolutionary sequences. For each value of the initial mass we considered stellar evolution from the zero age main sequence to central helium exhaustion with three values of the convective overshoot parameter: aov=0.1, 0.15 and 0.2. Models for the Cepheid SZ Tau with pulsation period 3.149 day can be constructed only for aov=0.1 and aov=0.15. The star is at the evolutionary stage of the second crossing of the instability strip and pulsates in the first overtone just near the boundary between the fundamental mode and the first overtone. The mass, radius and age of the star are in the ranges 5.46 <= M/Msun <= 5.75, 41.5 <= R/Rsun <= 42.3 and 6.9e7 yr <= t <= 8.0e7 yr, respectively. Predicted period change rates are of -0.4 s/yr.
The Starobinsky model of inflation, consistent with Planck 2015, has a peculiar form of the action, which contains the leading Einstein term R, the R^2 term with a huge coefficient, and negligible higher order terms. We propose an explanation of this form based on compactification of extra dimensions. Once tuning of order 10^(-4) is accepted to suppress the linear term R, we no more have to suppress higher order terms, which give nontrivial corrections to the Starobinsky model. We show our predictions of the spectral index, its running, and the tensor-to-scalar ratio. Finally, we discuss quantum gravity may appear at the scale greater than O(10^16) GeV.
A $^{238}$U projectile beam was used to create cadmium isotopes via abrasion-fission at 410 MeV/u in a beryllium target at the entrance of the in-flight separator FRS at GSI. The fission fragments were separated with the FRS and injected into the isochronous storage ring ESR for mass measurements. The Isochronous Mass Spectrometry (IMS) was performed under two different experimental conditions, with and without B$\rho$-tagging at the dispersive central focal plane of the FRS. In the experiment with B$\rho$-tagging the magnetic rigidity of the injected fragments was determined by an accuracy of $2\times 10^{-4}$. A new method of data analysis, using a correlation matrix for the combined data set from both experiments, has provided mass values for 25 different isotopes for the first time. The high selectivity and sensitivity of the experiment and analysis has given access even to rare isotopes detected with a few atoms per week. In this letter we present for the $^{129,130,131}$Cd isotopes mass values directly measured for the first time. The Cd results clearly show a very pronounced shell effect at $N=82$ which is in agreement with the conclusion from $\gamma$-ray spectroscopy of $^{130}$Cd and confirms the assumptions of modern shell-model calculations.
The flatness of the inflaton potential and lightness of the Higgs could have the common origin of the breaking of a global symmetry. This scenario provides a unified framework of Goldstone Inflation and Composite Higgs, where the inflaton and the Higgs both have a pseudo--Goldstone boson nature. The inflaton reheats the Universe via decays to the Higgs and subsequent secondary production of other SM particles via the top and massive vector bosons. We find that inflationary predictions and perturbative reheating conditions are consistent with CMB data for sub--Planckian values of the fields, as well as opening up the possibility of inflation at the TeV scale. We explore this exciting possibility, leading to an interplay between collider data and cosmological constraints.
The axions/axion like particles (ALPs) may constitute a major part of dark matter. Recently people find that dark matter axions can thermalize and form a Bose-Einstein condensate with a long correlation length. For the ALPs the thermalization scenario is similar. We find that for the linear regime of perturbation the ALPs are different from ordinary point like dark matter particles with additional terms in the first order velocity equation. The differences are especially compelling for string theory originated lighter ALPs. Also, axions/ALPs with a long correlation length can be thermalized due to gravitational interaction therefore alter the entropy of large scale. We propose that it can be a mechanism to explain the anomalies of Cosmic Microwave Background (CMB) multipoles if the mass of ALPs is order of $10^{-29}{\rm eV}$.
Based on mean field calculations with Skyrme interactions, we extract a constraint on the isovector effective mass in nuclear matter at saturation density $\rho_0$, i.e., $m_{v}^{\ast}(\rho_0)=(0.77\pm0.03) m$ by combining the experimental data of the centroid energy of the isovector giant dipole resonance (IVGDR) and the electric dipole polarizability in $^{208}$Pb. Meanwhile, the isoscalar effective mass at $\rho_0$ is determined to be $m_{s}^{\ast}(\rho_0)=(0.91\pm0.05) m$ by analyzing the measured excitation energy of the isoscalar giant quadrupole resonance (ISGQR) in $^{208}$Pb. From the constrained $m_{s}^{\ast}(\rho_0)$ and $m_{v}^{\ast}(\rho_0)$, we obtain the isospin splitting of nucleon effective mass in asymmetric nuclear matter of isospin asymmetry $\delta$ at $\rho_0$ as $(m_n^{\ast}(\rho_0,\delta)-m_p^{\ast}(\rho_0,\delta))/m = \Delta m^*_1(\rho_0) \delta + O(\delta^3)$ with the linear isospin splitting coefficient $\Delta m^*_1(\rho_0) = 0.33\pm0.16$. Furthermore, the constraints on $m_{v}^{\ast}(\rho)$, $m_{s}^{\ast}(\rho)$ and $\Delta m^*_1(\rho)$ at other densities are obtained from the similar analyses and we find that the $\Delta m^*_1(\rho)$ increases with the density.
Macroscopic irreversible processes emerge from fundamental physical laws of reversible character. The source of the local irreversibility seems to be not in the laws themselves but in the initial and boundary conditions of the equations that represent the laws. In this work we propose that the asymmetric screening of currents by black hole event horizons determines, locally, a preferred direction for the flux of electromagnetic energy. We study the growth of black hole event horizons due to the cosmological expansion and accretion of cosmic microwave background radiation, for different cosmological models. We propose generalized McVittie comoving metrics and integrate the rate of accretion of cosmic microwave background radiation onto a supermassive black hole over cosmic time. We find that for flat, open, and closed Friedmann cosmological models, the ratio of the total area of the black hole event horizons with respect to the area of a radial comoving space-like hypersurface is always larger than one. Since accretion of cosmic radiation set an absolute lower limit to the total matter accreted by black holes, this implies that the causal past and future are not symmetric for any spacetime event. The asymmetry causes a net Poynting flux in the global future direction; the latter is in turn related to the ever increasing thermodynamic entropy. Thus, we expose a connection between four different "time arrows": cosmological, electromagnetic, gravitational, and thermodynamic.
Using 137,562 quasars in the redshift range $2.1\leq z\leq3.5$ from the Data Release 11 (DR11) of the Baryon Oscillation Spectroscopic Survey (BOSS) of Sloan Digital Sky Survey (SDSS)-III, the BOSS-SDSS collaboration estimated the expansion rate $H(z=2.34)=222\pm7$ km/s/Mpc of Universe, and reported that this value is in tension with the predictions of flat $\Lambda$CDM model at around 2.5$\sigma$ level. In this letter, we briefly describe some attempts made in the literature to relieve the tension, and show that the tension can naturally be alleviated in non-flat $\Lambda$CDM model with positive curvature. However, this idea confronts with the inflation paradigm which predicts almost a spatially flat Universe. Nevertheless, the theoretical consistency of the non-flat $\Lambda$CDM model with the new result from BOSS deserves attention of the community.
We show that a new class of helical phase inflation models can be simply realized in minimal supergravity, wherein the inflaton is the phase component of a complex field and its potential admits a deformed helicoid structure. We find a new unique complex-valued index $\chi$ that characterizes almost the entire region of the $n_s-r$ plane favored by new Planck observations. Continuously varying the index $\chi$, predictions interpolate from quadratic/natural inflation parameterized by a phase/axion decay constant to Starobinsky-like inflation parameterized by the $\alpha$-parameter. We demonstrate that the simple supergravity construction realizing Starobinsky-like inflation can be obtained from a more microscopic model by integrating out heavy fields, and that the flat phase direction for slow-roll inflation is protected by a mildly broken global $U(1)$ symmetry. %, which is mildly broken at the inflation energy scale. We study the geometrical origin of the index $\chi$, and find that it corresponds to a linear constraint relating \kah moduli. We argue that such a linear constraint is a natural result of moduli stabilization in Type \MyRoman{2} orientifold compactifications on Calabi-Yau threefolds with geometric and non-geometric fluxes. Possible choices for the index $\chi$ are discrete points on the complex plane that relate to the distribution of supersymmetric Minkowski vacua on moduli space. More precise observations of the inflationary epoch in the future may provide a better estimation of the index $\chi$. Since $\chi$ is determined by the fluxes and vacuum expectation values of complex structure moduli, such observations would characterize the geometry of the internal space as well.
Links to: arXiv, form interface, find, astro-ph, recent, 1507, contact, help (Access key information)