We present kinematic and photometric evidence for an accretion event in the halo of the cD galaxy M87 in the last Gyr. Using velocities for ~300 planetary nebulas (PNs) in the M87 halo, we identify a chevron-like substructure in the PN phase-space. We implement a probabilistic Gaussian mixture model to identify the PNs that belong to the chevron. From analysis of deep V-band images of M87, we find that the region with the highest density of PNs associated to the chevron, is a crown-shaped substructure in the optical light. We assign a total of N_(PN,sub)=54 to the substructure, which extends over ~50 kpc along the major axis where we also observe radial variations of the ellipticity profile and a colour gradient. The substructure has highest surface brightness in a 20kpc x 60kpc region around 70 kpc in radius. In this region, it causes an increase in surface brightness by >60%. The accretion event is consistent with a progenitor galaxy with a V-band luminosity of L=2.8\pm1.0 x 10^9 L_(sun,V), a colour of (B-V)=0.76\pm0.05, and a stellar mass of M=6.4\pm2.3 x 10^9 M_sun. The accretion of this progenitor galaxy has caused an important modification of the outer halo of M87 in the last Gyr. By itself it is strong evidence that the galaxy's cD halo is growing through the accretion of smaller galaxies as predicted by hierarchical galaxy evolution models.
The density variance - Mach number relation of the turbulent interstellar medium is relevant for theoretical models of the star formation rate, efficiency, and the initial mass function of stars. Here we use high-resolution hydrodynamical simulations with grid resolutions of up to 1024^3 cells to model compressible turbulence in a regime similar to the observed interstellar medium. We use Fyris Alpha, a shock-capturing code employing a high-order Godunov scheme to track large density variations induced by shocks. We investigate the robustness of the standard relation between the logarithmic density variance (sigma_s^2) and the sonic Mach number (M) of isothermal interstellar turbulence, in the non-isothermal regime. Specifically, we test ideal gases with diatomic molecular (gamma = 7/5) and monatomic (gamma = 5/3) adiabatic indices. A periodic cube of gas is stirred with purely solenoidal forcing at low wavenumbers, leading to a fully-developed turbulent medium. We find that as the gas heats in adiabatic compressions, it evolves along the relationship in the density variance - Mach number plane, but deviates significantly from the standard expression for isothermal gases. Our main result is a new density variance - Mach number relation that takes the adiabatic index into account: sigma_s^2 = ln {1+b^2*M^[(5*gamma+1)/3]} and provides good fits for b*M <= 1. A theoretical model based on the Rankine-Hugoniot shock jump conditions is derived, sigma_s^2 = ln {1+(gamma+1)*b^2*M^2/[(gamma-1)*b^2*M^2+2]}, and provides good fits also for b*M > 1. We conclude that this new relation for adiabatic turbulence may introduce important corrections to the standard relation, if the gas is not isothermal.
Motivated by the recent discovery of several dwarf galaxies near the Large Magellanic Cloud (LMC), we study the accretion of massive satellites onto Milky Way (MW)/M31-like halos using the ELVIS suite of N-body simulations. We identify 25 surviving subhalos near the expected mass of the LMC, and investigate the lower-mass satellites that were associated with these subhalos before they fell into the MW/M31 halos. Typically, 7% of the overall z=0 satellite population of MW/M31 halos were in a surviving LMC-group prior to falling into the MW/M31 halo. This fraction, however, can vary between 1% and 25%, being higher for groups with higher-mass and/or more recent infall times. Groups of satellites disperse rapidly in phase space after infall, and their distances and velocities relative to the group center become statistically similar to the overall satellite population after 4-8 Gyr. We quantify the likelihood that satellites were associated with an LMC-mass group as a function of both distance and velocity relative to the LMC at z=0. The close proximity in distance of the nine Dark Energy Survey candidate dwarf galaxies to the LMC suggest that ~2-4 are likely associated with the LMC. Furthermore, if several of these dwarfs nearby to the LMC are genuine members, then the LMC-group probably fell into the MW very recently, <2 Gyr ago. If the connection with the LMC is established with the help of the follow-up velocity measurements, these "satellites of satellites" represent prime candidates to study the affects of group pre-processing on lower mass dwarfs.
We calculate the colours and luminosities of redshift z = 0.1 galaxies from the EAGLE simulation suite using the GALAXEV population synthesis models. We take into account obscuration by dust in birth clouds and diffuse ISM using a two-component screen model, following the prescription of Charlot and Fall. We compare models in which the dust optical depth is constant to models where it depends on gas metallicity, gas fraction and orientation. The colours of EAGLE galaxies for the more sophisticated models are in broad agreement with those of observed galaxies. In particular, EAGLE produces a red sequence of passive galaxies and a blue cloud of star forming galaxies, with approximately the correct fraction of galaxies in each population and with g-r colours within 0.1 magnitudes of those observed. Luminosity functions from UV to NIR wavelengths differ from observations at a level comparable to systematic shifts resulting from a choice between Petrosian and Kron photometric apertures. Despite the generally good agreement there are clear discrepancies with observations. The blue cloud of EAGLE galaxies extends to somewhat higher luminosities than in the data, consistent with the modest underestimate of the passive fraction in massive EAGLE galaxies. There is also a moderate excess of bright blue galaxies compared to observations. The overall level of agreement with the observed colour distribution suggests that EAGLE galaxies at z = 0.1 have ages, metallicities and levels of obscuration that are comparable to those of observed galaxies.
We show that the extended main sequence turnoffs seen in intermediate age Large Magellanic Cloud (LMC) clusters, often attributed to age spreads of several hundred Myr, may be easily accounted for by variable stellar rotation in a coeval population. We compute synthetic photometry for grids of rotating stellar evolution models and interpolate them to produce isochrones at a variety of rotation rates and orientations. An extended main sequence turnoff naturally appears in color-magnitude diagrams at ages just under 1 Gyr, peaks in extent between ~1 and 1.5 Gyr, and gradually disappears at around 2 Gyr in age. We then fit our interpolated isochrones by eye to four LMC clusters with very extended main sequence turnoffs: NGC 1783, 1806, 1846, and 1987. In each case, stellar populations with a single age and metallicity can comfortably account for the observed extent of the turnoff region.
[abridged] The ALESS survey has followed-up a sample of 122 sub-millimeter sources in the Extended Chandra Deep Field South at 870um with ALMA, allowing to pinpoint the positions of sub-millimeter galaxies (SMGs) to 0.3'' and to find their precise counterparts at different wavelengths. This enabled the first compilation of the multi-wavelength spectral energy distributions (SEDs) of a statistically reliable survey of SMGs. In this paper, we present a new calibration of the MAGPHYS modelling code that is optimized to fit these UV-to-radio SEDs of z>1 star-forming galaxies using an energy balance technique to connect the emission from stellar populations, dust attenuation and dust emission in a physically consistent way. We derive statistically and physically robust estimates of the photometric redshifts and physical parameters for the ALESS SMGs. We find that they have a median stellar mass $M_\ast=(8.9\pm0.1)\times10^{10} M_\odot$, SFR$=280\pm70 M_\odot$/yr, overall V-band dust attenuation $A_V=1.9\pm0.2$ mag, dust mass $M_\rm{dust}=(5.6\pm1.0)\times10^8 M_\odot$, and average dust temperature Tdust~40 K. The average intrinsic SED of the ALESS SMGs resembles that of local ULIRGs in the IR range, but the stellar emission of our average SMG is brighter and bluer, indicating lower dust attenuation, possibly because they are more extended. We explore how the average SEDs vary with different parameters, and we provide a new set of SMG templates. To put the ALESS SMGs into context, we compare their stellar masses and SFRs with those of less actively star-forming galaxies at the same redshifts. At z~2, about half of the SMGs lie above the star-forming main sequence, while half are at the high-mass end of the sequence. At higher redshifts (z~3.5), the SMGs tend to have higher SFR and Mstar, but the fraction of SMGs that lie significantly above the main sequence decreases to less than a third.
We investigate beta-interactions of free nucleons and their impact on the electron fraction (Y_e) and r-process nucleosynthesis in ejecta characteristic of binary neutron star mergers (BNSMs). For that we employ trajectories from a relativistic BNSM model to represent the density-temperature evolutions in our parametric study. In the high-density environment, positron captures decrease the neutron richness at the high temperatures predicted by the hydrodynamic simulation. Circumventing the complexities of modelling three-dimensional neutrino transport, (anti)neutrino captures are parameterized in terms of prescribed neutrino luminosities and mean energies, guided by published results and assumed as constant in time. Depending sensitively on the adopted neutrino-antineutrino luminosity ratio, neutrino processes increase Y_e to values between 0.25 and 0.40, still allowing for a successful r-process compatible with the observed solar abundance distribution and a significant fraction of the ejecta consisting of r-process nuclei. If the electron neutrino luminosities and mean energies are relatively large compared to the antineutrino properties, the mean Y_e might reach values >0.40 so that neutrino captures seriously compromise the success of the r-process. In this case, the r-abundances remain compatible with the solar distribution, but the total amount of ejected r-material is reduced to a few percent, because the production of iron-peak elements is favored. Proper neutrino physics, in particular also neutrino absorption, have to be included in BNSM simulations before final conclusions can be drawn concerning r-processing in this environment and concerning observational consequences like kilonovae, whose peak brightness and color temperature are sensitive to the composition-dependent opacity of the ejecta.
We present the discovery of a transiting exoplanet candidate in the K2 Field-1 with an orbital period of 9.1457 hours: EPIC 201637175b. The highly variable transit depths, ranging from $\sim$0\% to 1.3\%, are suggestive of a planet that is disintegrating via the emission of dusty effluents. We characterize the host star as an M-dwarf with $T_{\rm eff} \simeq 3800$. We have obtained ground-based transit measurements with several 1-m class telescopes and with the GTC. These observations (1) improve the transit ephemeris; (2) confirm the variable nature of the transit depths; (3) indicate variations in the transit shapes; and (4) demonstrate clearly that at least on one occasion the transit depths were significantly wavelength dependent. The latter three effects tend to indicate extinction of starlight by dust rather than by any combination of solid bodies. The K2 observations yield a folded light curve with lower time resolution but with substantially better statistical precision compared with the ground-based observations. We detect a significant ``bump' just after the transit egress, and a less significant bump just prior to transit ingress. We interpret these bumps in the context of a planet that is not only likely streaming a dust tail behind it, but also has a more prominent leading dust trail that precedes it. This effect is modeled in terms of dust grains that can escape to beyond the planet's Hill sphere and effectively undergo `Roche lobe overflow', even though the planet's surface is likely underfilling its Roche lobe by a factor of 2.
KIC 2835289 is a triple stellar system that consists of an inner ellipsoidal variable with an orbital period of ~0.86 days and an outer star that eclipses the inner pair every ~750 days. Two eclipse events were observed by the Kepler mission, but they do not fully constrain our photodynamical models. The next eclipse event will occur on May 14 UT, and we solicit community involvement to follow it up from the ground. All details are available in the attached call. Please contact the authors to join the follow-up campaign.
$\delta$-sunspots, with highly complex magnetic structures, are very productive in energetic eruptive events, such as X-class flares and homologous eruptions. We here study the formation of such complex magnetic structures by numerical simulations of magnetic flux emergence from the convection zone into the corona in an active-region-scale domain. In our simulation, two pairs of bipolar sunspots form on the surface, originating from two buoyant segments of a single subsurface twisted flux rope, following the approach of Toriumi et al. (2014). Expansion and rotation of the emerging fields in the two bipoles drive the two opposite polarities into each other with apparent rotating motion, producing a compact $\delta$-sunspot with a sharp polarity inversion line. The formation of the $\delta$-sunspot in such a realistic-scale domain produces emerging patterns similar to those formed in observations, e.g. the inverted polarity against Hale's law, the curvilinear motion of the spot, strong transverse field with highly sheared magnetic and velocity fields at the PIL. Strong current builds up at the PIL, giving rise to reconnection, which produces a complex coronal magnetic connectivity with non-potential fields in the Delta-spot overlaid by more relaxed fields connecting the two polarities at the two ends.
Ensemble modeling of CMEs provides a probabilistic forecast of CME arrival time which includes an estimation of arrival time uncertainty from the spread and distribution of predictions and forecast confidence in the likelihood of CME arrival. The real-time ensemble modeling of CME propagation uses the WSA-ENLIL+Cone model installed at the CCMC and executed in real-time. The current implementation evaluates the sensitivity of WSA-ENLIL+Cone model simulations of CME propagation to initial CME parameters. We discuss the results of real-time ensemble simulations for a total of 35 CME events between January 2013 - July 2014. For the 17 events where the CME was predicted to arrive at Earth, the mean absolute arrival time prediction error was 12.3 hours, which is comparable to the errors reported in other studies. For predictions of CME arrival at Earth the correct rejection rate is 62% and the false-alarm rate is 38%. The arrival time was within the range of the ensemble arrival predictions for 8 out of 17 events. The Brier Score for CME arrival predictions is 0.15 (where 1 is a perfect forecast), indicating that on average, the predicted likelihood of CME arrival is fairly accurate. The reliability of ensemble CME arrival predictions is heavily dependent on the initial distribution of CME input parameters, particularly the median and spread. Preliminary analysis of the probabilistic forecasts suggests undervariability, indicating that these ensembles do not sample a wide enough spread in CME input parameters. Prediction errors can also arise from ambient model parameters, the accuracy of the solar wind background derived from coronal maps, or other model limitations. Finally, predictions of the Kp geomagnetic index differ from observed values by less than one for 11 out of 17 of the ensembles and Kp prediction errors computed from the mean predicted Kp show a mean absolute error of 1.3.
Coherent alignments of galaxy shapes, often called "intrinsic alignments" (IA), are the most significant source of astrophysical uncertainty in weak lensing measurements. We develop the tidal alignment model of IA and demonstrate its success in describing observational data. We also describe a technique to separate IA from galaxy-galaxy lensing measurements. Applying this technique to luminous red galaxy lenses in the Sloan Digital Sky Survey, we constrain potential IA contamination from associated sources to be below a few percent.
In this work, we analyse coordinated observations spanning chromospheric, TR and coronal temperatures at very high resolution which reveal essential characteristics of thermally unstable plasmas. Coronal rain is found to be a highly multi-thermal phenomenon with a high degree of co-spatiality in the multi-wavelength emission. EUV darkening and quasi-periodic intensity variations are found to be strongly correlated to coronal rain showers. Progressive cooling of coronal rain is observed, leading to a height dependence of the emission. A fast-slow two-step catastrophic cooling progression is found, which may reflect the transition to optically thick plasma states. The intermittent and clumpy appearance of coronal rain at coronal heights becomes more continuous and persistent at chromospheric heights just before impact, mainly due to a funnel effect from the observed expansion of the magnetic field. Strong density inhomogeneities on spatial scales of 0.2"-0.5" are found, in which TR to chromospheric temperature transition occurs at the lowest detectable scales. The shape of the distribution of coronal rain widths is found to be independent of temperature with peaks close to the resolution limit of each telescope, ranging from 0.2" to 0.8". However we find a sharp increase of clump numbers at the coolest wavelengths and especially at higher resolution, suggesting that the bulk of the rain distribution remains undetected. Rain clumps appear organised in strands in both chromospheric and TR temperatures, suggesting an important role of thermal instability in the shaping of fundamental loop substructure. We further find structure reminiscent of the MHD thermal mode. Rain core densities are estimated to vary between 2x10^{10} cm^{-3} and 2.5x10^{11} cm^{-3} leading to significant downward mass fluxes per loop of 1-5x10^{9} g s^{-1}, suggesting a major role in the chromosphere-corona mass cycle.
Cosmological simulations still lack numerical resolution or physical processes to simulate dwarf galaxies in sufficient details. Accurate numerical simulations of individual dwarf galaxies are thus still in demand. We aim at (i) studying in detail the coupling between stars and gas in a galaxy, exploiting the so-called stellar hydrodynamical approach, and (ii) studying the chemo-dynamical evolution of individual galaxies starting from self-consistently calculated initial gas distributions. We present a novel chemo-dynamical code in which the dynamics of gas is computed using the usual hydrodynamics equations, while the dynamics of stars is described by the stellar hydrodynamics approach, which solves for the first three moments of the collisionless Boltzmann equation. The feedback from stellar winds and dying stars is followed in detail. In particular, a novel and detailed approach has been developed to trace the aging of various stellar populations, which enables an accurate calculation of the stellar feedback depending on the stellar age. We build initial equilibrium models of dwarf galaxies that take gas self-gravity into account and present different levels of rotational support. Models with high rotational support develop prominent bipolar outflows; a newly-born stellar population in these models is preferentially concentrated to the galactic midplane. Models with little rotational support blow away a large fraction of the gas and the resulting stellar distribution is extended and diffuse. The stellar dynamics turns out to be a crucial aspect of galaxy evolution. If we artificially suppress stellar dynamics, supernova explosions occur in a medium heated and diluted by the previous activity of stellar winds, thus artificially enhancing the stellar feedback (abridged).
In the era of rapidly increasing amounts of time series data, classification
of variable objects has become the main objective of time-domain astronomy.
Classification of irregularly sampled time series is particularly difficult
because the data can not be represented naturally as a plain vector which
directly can be fed into a classifier. In the literature, various statistical
features derived from time series serve as a representation. Typically, the
usefulness of the derived features is judged in an empirical fashion according
to their predictive power.
In this work, an alternative to the feature-based approach is investigated.
In this new representation the time series is described by a density model.
Similarity between each pair of time series is quantified by the distance
between their respective models. The density model captures all the information
available, also including measurement errors. Hence, we view this model as a
generalisation to the static features which directly can be derived, e.g., as
moments from the density.
In the numerical experiments, we use data from the OGLE and ASAS surveys and
demonstrate that the proposed representation performs up to par with the best
currently used feature-based approaches. While the density representation
preserves all static information present in the observational data, the
features are only a less complete description. The density representation is an
upper boundary in terms of information made available to the classifier.
Consequently, the predictive power of the proposed classification depends on
the choice of similarity measure and classifier, only. We therefore expect that
the proposed method yields performance close to an optimal classifier. Due to
its principled nature, we advocate that this new approach of representing time
series has potential in tasks beyond classification, e.g., unsupervised
learning.
A number of possible mechanisms have been suggested to generate density inhomogeneities in the early Universe which could survive until the onset of primordial nucleosynthesis and generate neutron-rich regions. In this work we are not concerned with how the inhomogeneities were generated but we want to focus on the effect of such inhomogeneities on primordial nucleosynthesis. One of the proposed signatures of inhomogeneity, the synthesis of very heavy elements by neutron capture in a primordial r-process, was analyzed for varying baryon to photon ratios $\eta$ and fluctuation length scales $L$.
Baryon density inhomogeneities in the early universe can give rise to a floor of heavy elements (up to $A\approx 270$) produced in a primordial r-process with fission cycling. A parameter study with variation of the global baryon to photon ratio $\eta$, and under inclusion of neutron diffusion effects was performed. New results concerning the dependence of the results on nuclear physics parameters are presented.
Water maser emission at 22 GHz is a useful probe to study the transition between the nearly spherical mass-loss in the AGB to a collimated one in the post-AGB phase. In their turn, collimated jets in the post-AGB phase could determine the shape of planetary nebulae (PNe) once photoionization starts. We intend to find new cases of post-AGB stars and PNe with water maser emission, including water fountains or water-maser-emitting PNe. We observed water maser emission in a sample of 133 objects, with a significant fraction being post-AGB and young PN candidate sources with strong obscuration. We detected this emission in 15 of them, of which seven are reported here for the first time. We identified three water fountain candidates: IRAS 17291-2147, with a total velocity spread of ~96 km/s in its water maser components and two sources (IRAS 17021-3109 and IRAS 17348-2906) that show water maser emission outside the velocity range covered by OH masers. We have also identified IRAS 17393-2727 as a possible new water-maser-emitting PN. The detection rate is higher in obscured objects (14%) than in those with optical counterparts (7%), consistent with previous results. Water maser emission seems to be common in objects that are bipolar in the near-IR (43% detection rate). The water maser spectra of water fountain candidates like IRAS 17291-2147 show significantly less maser components than others (e.g., IRAS 18113-2503). We speculate that most post-AGBs may show water maser emission with wide enough velocity spread (> 100 km/s) when observed with enough sensitivity and/or for long enough periods of time. Therefore, it may be necessary to single out a special group of "water fountains", probably defined by their high maser luminosities. We also suggest that the presence of both water and OH masers in a PN is a better tracer of its youth, rather than the presence of just one of these species.
I present an exact and explicit solution to the scalar (Stokes flux intensity) radio interferometer imaging equation on a spherical surface which is valid also for non-coplanar interferometer configurations. This imaging equation is comparable to $w$-term imaging algorithms, but by using a spherical rather than a Cartesian formulation this term has no special significance. The solution presented also allows direct identification of the scalar (spin 0 weighted) spherical harmonics on the sky. The method should be of interest for future multi-spacecraft interferometers, wide-field imaging with non-coplanar arrays, and CMB spherical harmonic measurements using interferometers.
Knowledge of the geometry of pulsational modes is a prerequisite for seismic modelling of stars. In the case of slowly pulsating B-type (SPB) pulsators, the simple zero-rotation approach so far used for mode identification is usually not valid because pulsational frequencies are often of the order of the rotational frequency. Moreover, this approach allows us to determine only the spherical harmonic degree, $\ell$, while the azimuthal order, $m$, is beyond its reach. On the other hand, because of the density of oscillation spectra of SPB stars, knowledge of $m$ is indispensable if one wants to assign the radial order, $n$, to the observed frequency peaks. Including the effects of rotation via the traditional approximation, we perform identification of the mode angular numbers ($\ell,~m$) for 31 SPB stars with available multicolour time series photometry. Simultaneously, constraints on the rotational velocity, $V_{\rm rot}$, and the inclination angle, $i$, are determined assuming uniform rotation and a constant value of $V_{\rm rot}\sin i$. Dependence of the results on the adopted model is tested using HD\,21071 as an example. Despite some model uncertainties and limitations of the method, our studies show the correct approach to identifying the low-frequency oscillation modes.
We develop an efficient method based on the linear regression algorithm to probe the cosmological CPT violation using the CMB polarisation data. We validate this method using simulated CMB data and apply it to recent CMB observations. We find that a combined data sample of BICEP1 and BOOMERanG 2003 favours a nonzero isotropic rotation angle at $2.3\sigma$ confidence level, ie, $\Delta\alpha=-3.3 \pm1.4$ deg (68% CL) with systematics included.
The Rosetta mission and its exquisite measurements have revived the debate on whether comets are pristine planetesimals or collisionally evolved objects. We investigate the collisional evolution experienced by the precursors of current comet nuclei during the early stages of the Solar System, in the context of the so-called "Nice Model". We consider two environments for the collisional evolution: (1) the trans-planetary planetesimal disk, from the time of gas removal until the disk was dispersed by the migration of the ice giants, and (2) the dispersing disk during the time that the scattered disk was formed. Simulations have been performed, using different methods in the two cases, to find the number of destructive collisions typically experienced by a comet nucleus of 2km radius. In the widely accepted scenario, where the dispersal of the planetesimal disk occurred at the time of the Late Heavy Bombardment about 4Gy ago, comet-sized planetesimals have a very small chance to survive against destructive collisions in the disk. On the extreme assumption that the disk was dispersed directly upon gas removal, there is a chance for a significant fraction of the planetesimals to remain intact. However, these survivors would still bear the marks of many non-destructive impacts. Thus, the Nice Model of Solar System evolution predicts that typical km-sized comet nuclei are predominantly fragments resulting from collisions experienced by larger parent bodies. An important goal for further research is to investigate, whether the observed properties of comet nuclei are compatible with such a collisional origin.
We present predictions for the clustering of galaxies selected by their total infra-red luminosity ($L_{\rm IR}$), and their emission at far infra-red (FIR) and sub-millimetre (sub-mm) wavelengths. We combine a new version of the GALFORM semi-analytic model of galaxy formation, implemented in a Millennium-style $N$-body simulation utilising the WMAP7 cosmology, with a self-consistent model for calculating the absorption and re-emission of stellar radiation by dust. In the model, galaxies selected at 850 $\mu$m predominantly reside in dark matter halos of mass ~$10^{11.5}-10^{12}$ $h^{-1}$ M$_{\odot}$, independent of redshift (for $0.2<z<4$) or flux (for $0.25<S_{850\mu\rm m}<4$ mJy). Around the peak of their redshift distribution ($z$~2.5) the brightest galaxies ($S_{850 \mu\rm m}$>4 mJy) exhibit a correlation length of $r_{0}=5.5_{-0.5}^{+0.3}$ $h^{-1}$ Mpc, consistent with observations. We show further that these galaxies evolve into $z=0$ descendants with stellar mass ~$10^{11}$ $h^{-1}$ M$_{\odot}$ occupying halos which span a broad range in mass ~$10^{12}-10^{14}$ $h^{-1}$ M$_{\odot}$. The FIR emissivity at shorter wavelengths (250, 350 and 500 $\mu$m) in our model is dominated by galaxies in the same halo mass range, again independent of redshift (for $0.5<z<5$). We compare our predictions for the angular power spectrum of Cosmic Infra-red Background (CIB) anisotropies at these wavelengths with recent observations, and find that the model agrees with the observed power to within a factor of ~2 over all scales and wavelengths, an improvement over earlier versions of the model. Simulating sub-mm imaging at 850 $\mu$m, we show that source confusion due to the coarse angular resolution of single-dish telescopes at this wavelength can significantly bias angular clustering measurements, severely complicating the interpretation of such observations.
We present a study of the molecular CO gas and mid/far infrared radiation arising from the environment surrounding the Wolf-Rayet (W-R) star 130. We use the multi-wavelength data to analyze the properties of the dense gas and dust, and its possible spatial correlation with that of Young Stellar Objects (YSOs). We use CO J=1-0 data from the FCRAO survey as tracer of the molecular gas, and mid/far infrared data from the recent WISE and Herschel space surveys to study the dust continuum radiation and to identify a population of associated candidate YSOs. The spatial distribution of the molecular gas shows a ring-like structure very similar to that observed in the HI gas, and over the same velocity interval. The relative spatial distribution of the HI and CO components is consistent with a photo-dissociation region. We have identified and characterized four main and distinct molecular clouds that create this structure. Cold dust is coincident with the dense gas shown in the CO measurements. We have found several cYSOs that lie along the regions with the highest gas column density, and suggest that they are spatially correlated with the shell. These are indicative of regions of star formation induced by the strong wind and ionization of the WR star.
We describe and execute a novel approach to observationally estimate the lifetimes of giant molecular clouds (GMCs). We focus on the cloud population between the two main spiral arms in M51 (the inter-arm region) where cloud destruction via shear and star formation feedback dominates over formation processes. By monitoring the change in GMC number densities and properties from one side of the inter-arm to the other, we estimate the lifetime as a fraction of the inter-arm travel time. We find that GMC lifetimes in M51's inter-arm are finite and short, 20 to 30 Myr. Such short lifetimes suggest that cloud evolution is influenced by environment, in which processes can disrupt GMCs after a few free-fall times. Over most of the region under investigation shear appears to regulate the lifetime. As the shear timescale increases with galactocentric radius, we expect cloud destruction to switch primarily to star formation feedback at larger radii. We identify a transition from shear- to feedback-dominated disruption through a change in the behavior of the GMC number density. The signature suggests that shear is more efficient at completely dispersing clouds, whereas feedback transforms the population, e.g. by fragmenting high mass clouds into lower mass pieces. Compared to the characteristic timescale for molecular hydrogen in M51, our short lifetimes suggest that gas can remain molecular while clouds disperse and reassemble. We propose that galaxy dynamics regulates the cycling of molecular material from diffuse to bound (and ultimately star-forming) objects, contributing to long observed molecular depletion times in normal disk galaxies. We also speculate that, in more extreme environments such as elliptical galaxies and concentrated galaxy centers, star formation can be suppressed when the shear timescale becomes so short that some clouds can not survive to collapse and form stars.
We describe the search for Lyman-break galaxies (LBGs) near the sub-millimeter bright starburst galaxy HFLS3 at $z$$=$6.34 and a study on the environment of this massive galaxy during the end of reionization. We performed two independent selections of LBGs on images obtained with the Gran Telescopio Canarias (GTC) and the Hubble Space Telescope (HST) by combining non-detections in bands blueward of the Lyman-break and color selection. A total of 10 objects fulfilling the LBG selection criteria at $z$$>$5.5 were selected over the 4.54 and 55.5 arcmin$^2$ covered by our HST and GTC images, respectively. The photometric redshift, UV luminosity, and the star-formation rate of these sources were estimated with models of their spectral energy distribution. These $z$$\sim$6 candidates have physical properties and number densities in agreement with previous results. The UV luminosity function of this field at $z$$\sim$6 shows no strong evidence for an overdensity of relatively bright objects (m$_{F105W}$$<$25.9) associated with HFLS3. A Voronoi tessellation analysis also did not allow a detection of an overdensity around HFLS3. However we identified three faint objects at less than three arcseconds from HFLS3 with color consistent with those expected for $z$$\sim$6 galaxies. Deeper data are needed to confirm their redshifts and to study their association with HFLS3 and the galaxy merger that may be responsible for the massive starburst.
The neutron-drip transition in the dense matter constituting the interior of neutron stars generally refers to the appearance of unbound neutrons as the matter density reaches some threshold density $\rho_\textrm{drip}$. This transition has been mainly studied under the cold catalyzed matter hypothesis. However, this assumption is unrealistic for accreting neutron stars. After examining the physical processes that are thought to be allowed in both accreting and nonaccreting neutron stars, suitable conditions for the onset of neutron drip are derived and general analytical expressions for the neutron drip density and pressure are obtained. Moreover, we show that the neutron-drip transition occurs at lower density and pressure than those predicted within the mean-nucleus approximation. This transition is studied numerically for various initial composition of the ashes from X-ray bursts and superbursts using microscopic nuclear mass models.
We attempt to determine the molecular composition of disks around young low-mass stars in the $\rho$ Oph region and to compare our results with a similar study performed in the Taurus-Auriga region. We used the IRAM 30 m telescope to perform a sensitive search for CN N=2-1 in 29 T Tauri stars located in the $\rho$ Oph and upper Scorpius regions. $^{13}$CO J=2-1 is observed simultaneously to provide an indication of the level of confusion with the surrounding molecular cloud. The bandpass also contains two transitions of ortho-H$_2$CO, one of SO, and the C$^{17}$O J=2-1 line, which provides complementary information on the nature of the emission. Contamination by molecular cloud in $^{13}$CO and even C$^{17}$O is ubiquitous. The CN detection rate appears to be lower than for the Taurus region, with only four sources being detected (three are attributable to disks). H$_2$CO emission is found more frequently, but appears in general to be due to the surrounding cloud. The weaker emission than in Taurus may suggest that the average disk size in the $\rho$ Oph region is smaller than in the Taurus cloud. Chemical modeling shows that the somewhat higher expected disk temperatures in $\rho$ Oph play a direct role in decreasing the CN abundance. Warmer dust temperatures contribute to convert CN into less volatile forms. In such a young region, CN is no longer a simple, sensitive tracer of disks, and observations with other tracers and at high enough resolution with ALMA are required to probe the gas disk population.
Aims. We aim to describe the pre-main sequence and main-sequence evolution of
X-ray and extreme-ultaviolet radiation of a solar mass star based on its
rotational evolution starting with a realistic range of initial rotation rates.
Methods. We derive evolutionary tracks of X-ray radiation based on a
rotational evolution model for solar mass stars and the rotation-activity
relation. We compare these tracks to X-ray luminosity distributions of stars in
clusters with different ages.
Results. We find agreement between the evolutionary tracks derived from
rotation and the X-ray luminosity distributions from observations. Depending on
the initial rotation rate, a star might remain at the X-ray saturation level
for very different time periods, approximately from 10 Myr to 300 Myr for slow
and fast rotators, respectively.
Conclusions. Rotational evolution with a spread of initial conditions leads
to a particularly wide distribution of possible X-ray luminosities in the age
range of 20 to 500 Myrs, before rotational convergence and therefore X-ray
luminosity convergence sets in. This age range is crucial for the evolution of
young planetary atmospheres and may thus lead to very different planetary
evolution histories.
A site selection of potential observatory locations in Turkey have been carried out by using Multi-Criteria Decision Analysis (MCDA) coupled with Geographical Information Systems (GIS) and satellite imagery which in turn reduced cost and time and increased the accuracy of the final outcome. The layers of cloud cover, digital elevation model, artificial lights, precipitable water vapor, aerosol optical thickness and wind speed were studied in the GIS system. In conclusion of MCDA, the most suitable regions were found to be located in a strip crossing from southwest to northeast including also a diverted region in southeast of Turkey. These regions are thus our prime candidate locations for future on-site testing. In addition to this major outcome, this study has also been applied to locations of major observatories sites. Since no goal is set for \textit{the best}, the results of this study is limited with a list of positions. Therefore, the list has to be further confirmed with on-site tests. A national funding has been awarded to produce a prototype of an on-site test unit (to measure both astronomical and meteorological parameters) which might be used in this list of locations.
Binary mergers (NSMs) of double neutron star (and black hole-neutron star) systems are suggested to be major sites of r-process elements in the Galaxy by recent hydrodynamical and nucleosynthesis studies. It has been pointed out, however, that the estimated long lifetimes of neutron star binaries are in conflict with the presence of r-process-enhanced halo stars at metallicities as low as [Fe/H] ~ -3. To resolve this problem, we examine the role of NSMs in the early Galactic chemical evolution on the assumption that the Galactic halo was formed from merging sub-halos. We present simple models for the chemical evolution of sub-halos with total final stellar masses between 10^4 M_solar and 2 x 10^8 M_solar. Typical lifetimes of compact binaries are assumed to be 100 Myr (for 95% of their population) and 1 Myr (for 5%), according to recent binary population synthesis studies. The resulting metallcities of sub-halos and their ensemble are consistent with the observed mass-metallicity relation of dwarf galaxies in the Local Group, and the metallicity distribution of the Galactic halo, respectively. We find that the r-process abundance ratios [r/Fe] start increasing at [Fe/H] <= -3 if the star formation efficiencies are smaller for less massive sub-halos. In addition, the sub-solar [r/Fe] values (observed as [Ba/Fe] ~ -1.5 for [Fe/H] < -3) are explained by the contribution from the short-lived (~1 Myr) binaries. Our results indicate that NSMs may have a substantial contribution to the r-process element abundances throughout the Galactic history.
We present the results of approximately three years of observations of Planck Sunyaev-Zeldovich (SZ) sources with telescopes at the Canary Islands observatories, as part of the general optical follow-up programme undertaken by the Planck collaboration. In total, 78 SZ sources are discussed. Deep imaging observations were obtained for most of those sources; spectroscopic observations in either in long-slit or multi-object modes were obtained for many. We found optical counterparts for 73 of the 78 candidates. This sample includes 53 spectroscopic redshifts determinations, 20 of them obtained with a multi-object spectroscopic mode. The sample contains new redshifts for 27 Planck clusters that were not included in the first Planck SZ source catalogue (PSZ1).
We study the non-perturbative dynamics of the Standard Model (SM) after inflation, in the regime where the SM is decoupled from (or weakly coupled to) the inflationary sector. We use classical lattice simulations in an expanding box in (3+1) dimensions, modeling the SM gauge interactions with both global and Abelian-Higgs analogue scenarios. We consider different post-inflationary expansion rates. During inflation, the Higgs forms a condensate, which starts oscillating soon after inflation ends. Via non-perturbative effects, the oscillations lead to a fast decay of the Higgs into the SM species, transferring most of the energy into $Z$ and $W^{\pm}$ bosons. All species are initially excited far away from equilibrium, but their interactions lead them into a stationary stage, with exact equipartition among the different energy components. From there on the system eventually reaches equilibrium. We have characterized in detail, in the different expansion histories considered, the evolution of the Higgs and of its dominant decay products, until equipartition is established.
The aim of this letter is to discuss the virtual identity of two recent tidal theories: the creep tide theory of Ferraz-Mello (Cel. Mech. Dyn. Astron. 116, 109, 2013) and the Maxwell model developed by Correia et al. (Astron. Astrophys. 571, A50, 2014). It includes the discussion of the basic equations of the theories, which, in both cases, include an elastic and an anelastic component, and shows that the basic equations of the two theories are equivalent and differ by only a numerical factor in the anelastic tide. It also includes a discussion of the lags: the lag of the full tide (geodetic), dominated by the elastic component, and the phase of the anelastic tide. In rotating rocky bodies not trapped in a spin-orbit resonance (e.g., the Earth) the geodetic lag is close to zero and the phase of the anelastic tide is close to 90 degrees. The results obtained from combining tidal solutions from satellite tracking data and from Topex/Poseidon satellite altimeter data, by Ray et al., are extended to determine the phase of the Earth's anelastic tide as $ \sigma_0=89.80\pm 0.05$ degrees.
We present here the first spectroscopic and photometric analysis of the
double-lined eclipsing binary containing the classical, first-overtone Cepheid
OGLE-LMC-CEP-2532 (MACHO 81.8997.87). The system has an orbital period of 800
days and the Cepheid is pulsating with a period of 2.035 days.
Using spectroscopic data from three high-class telescopes and photometry from
three surveys spanning 7500 days we are able to derive the dynamical masses for
both stars with an accuracy better than 3%. This makes the Cepheid in this
system one of a few classical Cepheids with an accurate dynamical mass
determination (M_1=3.90 +/- 0.10 M_sun). The companion is probably slightly
less massive (3.82 +/- 0.10 M_sun), but may have the same mass within errors
(M_2/M_1= 0.981 +/- 0.015). The system has an age of about 185 million years
and the Cepheid is in a more advanced evolutionary stage.
For the first time precise parameters are derived for both stars in this
system. Due to the lack of the secondary eclipse for many years not much was
known about the Cepheid's companion. In our analysis we used extra information
from the pulsations and the orbital solution from the radial velocity curve.
The best model predicts a grazing secondary eclipse shallower than 1 mmag,
hence undetectable in the data, about 370 days after the primary eclipse.
The dynamical mass obtained here is the most accurate known for a
first-overtone Cepheid and will contribute to the solution of the Cepheid mass
discrepancy problem.
One of the most intriguing results from the gamma-ray instruments in orbit
has been the detection of powerful flares from the Crab Nebula. These flares
challenge our understanding of pulsar wind nebulae and models for particle
acceleration. We report on the portion of a multiwavelength campaign using
Keck, HST, and Chandra concentrating on a small emitting region, the Crab's
inner knot, located a fraction of an arcsecond from the pulsar.
We find that the knot's radial size, tangential size, peak flux, and the
ratio of the flux to that of the pulsar are correlated with the projected
distance of the knot from the pulsar. A new approach, using singular value
decomposition for analyzing time series of images, was introduced yielding
results consistent with the more traditional methods while some uncertainties
were substantially reduced.
We exploit the characterization of the knot to discuss constraints on
standard shock-model parameters that may be inferred from our observations
assuming the inner knot lies near to the shocked surface. These include
inferences as to wind magnetization, shock shape parameters such as incident
angle and poloidal radius of curvature, as well as the IR/optical emitting
particle enthalpy fraction. We find that while the standard shock model gives
good agreement with observation in many respects, there remain two puzzles: (a)
The observed angular size of the knot relative to the pulsar--knot separation
is much smaller than expected; (b) The variable, yet high degree of
polarization reported is difficult to reconcile with a highly relativistic
downstream flow.
We study the effects of energy transport in the Sun by asymmetric dark matter with momentum and velocity-dependent interactions, with an eye to solving the decade-old Solar Abundance Problem. We study effective theories where the dark matter-nucleon scattering cross-section goes as $v_{\rm rel}^{2n}$ and $q^{2n}$ with $n = -1, 0, 1 $ or $2$, where $v_{\rm rel}$ is the dark matter-nucleon relative velocity and $q$ is the momentum exchanged in the collision. Such cross-sections can arise generically as leading terms from the most basic nonstandard DM-quark operators. We employ a high-precision solar simulation code to study the impact on solar neutrino rates, the sound speed profile, convective zone depth, surface helium abundance and small frequency separations. We find that the majority of models that improve agreement with the observed sound speed profile and depth of the convection zone also reduce neutrino fluxes beyond the level that can be reasonably accommodated by measurement and theory errors. However, a few specific points in parameter space yield a significant overall improvement. A 3-5 GeV DM particle with $\sigma_{SI} \propto q^2$ is particularly appealing, yielding more than a $6\sigma$ improvement with respect to standard solar models, while being allowed by direct detection and collider limits. We provide full analytical capture expressions for $q$- and $v_{\rm rel}$-dependent scattering, as well as complete likelihood tables for all models.
Scalar fields appear in many theories beyond the Standard Model of particle physics. In the early universe, they are exposed to extreme conditions, including high temperature and rapid cosmic expansion. Understanding their behavior in this environment is crucial to understand the implications for cosmology. We calculate the finite temperature effective action for the field expectation value in two particularly important cases, for damped oscillations near the ground state and for scalar fields with a flat potential. We find that the behavior in both cases can in good approximation be described by a complex valued effective potential that yields Markovian equations of motion. Near the potential minimum, we recover the solution to the well-known Langevin equation. For large field values we find a very different behavior, and our result for the damping coefficient significantly differs from the expressions given in the literature. We illustrate our results in a simple scalar model, for which we give analytic approximations for the effective potential and damping coefficient. We also provide various expressions for loop integrals at finite temperature that are useful for future calculations in other models.
We review recent developments on cosmology in extended teleparallel gravity, called "$F(T)$ gravity" with $T$ the torsion scalar in teleparallelism. We explore various cosmological aspects of $F(T)$ gravity including the evolution of the equation of state for the universe, finite-time future singularities, thermodynamics, and four-dimensional effective $F(T)$ gravity theories coming from the higher-dimensional Kaluza-Klein (KK) and Randall-Sundrum (RS) theories.
We present the results from the first ensemble prediction model for major solar flares (M and X classes). Using the probabilistic forecasts from three models hosted at the Community Coordinated Modeling Center (NASA-GSFC) and the NOAA forecasts, we developed an ensemble forecast by linearly combining the flaring probabilities from all four methods. Performance-based combination weights were calculated using a Monte Carlo-type algorithm by applying a decision threshold $P_{th}$ to the combined probabilities and maximizing the Heidke Skill Score (HSS). Using the probabilities and events time series from 13 recent solar active regions (2012 - 2014), we found that a linear combination of probabilities can improve both probabilistic and categorical forecasts. Combination weights vary with the applied threshold and none of the tested individual forecasting models seem to provide more accurate predictions than the others for all values of $P_{th}$. According to the maximum values of HSS, a performance-based weights calculated by averaging over the sample, performed similarly to a equally weighted model. The values $P_{th}$ for which the ensemble forecast performs the best are 25 \% for M-class flares and 15 \% for X-class flares. When the human-adjusted probabilities from NOAA are excluded from the ensemble, the ensemble performance in terms of the Heidke score, is reduced.
A new perspective on quintessence cosmology with an exponential potential $V(\phi)=V_{0}\exp(-\lambda\phi)$ is introduced. Different from the traditional phase space analysis, we analyze the cosmological parameters' space. We make a comparison of these two analyses and demonstrate many advantages of this new perspective. The whole evolution history of the universe is available in the cosmological parameters' space analysis. Associating the observational data, the current state of universe could be pinpointed in phase diagram, which could offer more information. We find that the stable attractor solutions are the final states of cosmic evolution instead of current state, while the intermediate evolutionary state is quite uncertain. We also numerically present valid regions of the initial conditions of the autonomous system. Using the numerical method, we obtain some results in an agreement with the phase space analysis. Moreover, the numerical method is appropriate to any dimension.
The first measurement on the antiproton to proton ratio made by the AMS-02 collaboration agrees with the expection from conventional cosmic-ray secondaries in the kinetic energy range $\sim 10-100$ GeV, which can be turned into stringent upper limits on the dark matter (DM) annihilation cross sections above $\sim 300$ GeV. Using the GALPROP code, we derive the upper limits in various propagation models and DM profiles. We show that in the "conventional" propagation model, for the $q\bar q$, $b\bar b$, and $WW$ final states, the constraints can be more stringent than that derived from the recent Ferm-LAT gamma-ray data on the dwarf spheroidal satellite galaxies. Making use of the typical minimal, median and maximal models obtained from a previous global fit, we show that the variation of the upper limits is around a factor of five.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
Highest resolution imaging in astronomy is achieved by interferometry, connecting telescopes over increasingly longer distances, and at successively shorter wavelengths. Here, we present the first diffraction-limited images in visual light, produced by an array of independent optical telescopes, connected electronically only, with no optical links between them. With an array of small telescopes, second-order optical coherence of the sources is measured through intensity interferometry over 180 baselines between pairs of telescopes, and two-dimensional images reconstructed. The technique aims at diffraction-limited optical aperture synthesis over kilometre-long baselines to reach resolutions showing details on stellar surfaces and perhaps even the silhouettes of transiting exoplanets. Intensity interferometry circumvents problems of atmospheric turbulence that constrain ordinary interferometry. Since the electronic signal can be copied, many baselines can be built up between dispersed telescopes, and over long distances. Using arrays of air Cherenkov telescopes, this should enable the optical equivalent of interferometric arrays currently operating at radio wavelengths.
We investigate the properties of halo gas using three cosmological `zoom-in' simulations of realistic Milky Way-galaxy analogs with varying sub-grid physics. In all three cases, the mass of hot ($T > 10^6$ K) halo gas is $\sim 1\%$ of the host's virial mass. Hot halos extend to 140 kpc from the galactic center and are surrounded by a bubble of warm-hot ($T = 10^5 - 10^6$K) gas that extends to the virial radius. Simulated halos agree well outside 20-30 kpc with the $\beta$-model of Miller $\&$ Bregman (2014) based on OVII absorption and OVIII emission measurements. Warm-hot and hot gas contribute up to $80\%$ of the total gas reservoir, and contain nearly the same amount of baryons as the stellar component. The mass of warm-hot and hot components falls into the range estimated for $L^*$ galaxies. With key observational constraints on the density of the Milky Way corona being satisfied, we show that concealing of the ubiquitous warm-hot baryons, along with the ejection of just $20-30 \%$ of the diffuse gas out of the potential wells by supernova-driven outflows, can solve the "missing baryon problem". The recovered baryon fraction within 3 virial radii is $90\%$ of the universal value. With a characteristic density of $\sim 10^{-4}$ cm$^{-3}$ at $50-80$ kpc, diffuse coronae meet the requirement for fast and complete ram-pressure stripping of the gas reservoirs in dwarf galaxy satellites, which signals the importance of satellite accretion in the assembly of halos and explains naturally how dSphs lost their gas soon after infall.
Einstein's theory of General Relativity (GR) is tested accurately within the local universe i.e., the Solar System, but this leaves open the possibility that it is not a good description at the largest scales in the Universe. The standard model of cosmology assumes GR as the theory to describe gravity on all scales. In 1998, astronomers made the surprising discovery that the expansion of the Universe is accelerating, not slowing down. This late-time acceleration of the Universe has become the most challenging problem in theoretical physics. Within the framework of GR, the acceleration would originate from an unknown dark energy. Alternatively, it could be that there is no dark energy and GR itself is in error on cosmological scales. The standard model of cosmology is based on a huge extrapolation of our limited knowledge of gravity. This discovery of the late time acceleration of the Universe may require us to revise the theory of gravity and the standard model of cosmology based on GR. We will review recent progress in constructing modified gravity models as an alternative to dark energy and developing cosmological tests of gravity.
Several hundred young stars lie in the innermost parsec of our Galaxy. The super-massive black hole (SMBH) might capture planets orbiting these stars, and bring them onto nearly radial orbits. The same fate might occur to planetary embryos (PEs), i.e. protoplanets born from gravitational instabilities in protoplanetary disks. In this paper, we investigate the emission properties of rogue planets and PEs in the Galactic center. In particular, we study the effects of photoevaporation, caused by the ultraviolet background. Rogue planets can hardly be detected by current or forthcoming facilities, unless they are tidally disrupted and accrete onto the SMBH. In contrast, photoevaporation of PEs (especially if the PE is being tidally stripped) might lead to a recombination rate as high as ~10^45 s^-1, corresponding to a Brackett-gamma luminosity ~10^31 erg s^-1, very similar to the observed luminosity of the dusty object G2. We critically discuss the possibility that G2 is a rogue PE, and the major uncertainties of this model.
We confirm the planetary nature of KOI-372b (aka Kepler object of interest K00372.01), a giant transiting exoplanet orbiting a solar-analog G2V star. The mass of KOI-372b and the eccentricity of its orbit were accurately derived thanks to a series of precise radial velocity measurements obtained with the CAFE spectrograph mounted on the CAHA 2.2-m telescope. A simultaneous fit of the radial-velocity data and Kepler photometry revealed that KOI-372b is a dense Jupiter-like planet with a mass of Mp=3.25 Mjup and a radius of Rp=0.882 Rjup. KOI-372b is moving on a quite eccentric orbit, e=0.172, making a complete revolution around its parent star in 125.6 days. The semi-major axis of the orbit is 0.4937 au, implying that the planet is close to its habitable zone (roughly 0.5 au from it). By analysing the mid-transit times of the 12 transit events of KOI-372b recorded by the Kepler spacecraft, we found a clear transit time variation, which is attributable to the presence of a planet c in a wider orbit. We estimated that KOI-372c has a mass between 0.13 and 0.31 Mjup, also revolving on an eccentric orbit (e=0.17-0.24) in roughly 460 days, at a mean distance of 1.2 au from the host star, within the boundaries of its habitable zone. The analysis of the CAFE spectra revealed a relatively high photospheric lithium content, A(Li)=2.48 dex, suggesting that the parent star is relatively young. From a gyrochronological analysis, we estimate that the age of this planetary system is 1.0 Gyr.
Many core-collapse supernova progenitors are presumed to be in binary systems. If a star explodes in a binary system, the early supernova light curve can be brightened by the collision of the supernova ejecta with the companion star. The early brightening can be observed when the observer is in the direction of the hole created by the collision. Based on a population synthesis model, we estimate the fractions of core-collapse supernovae in which the light-curve brightening by the collision can be observed. We find that 0.19% of core-collapse supernova light curves can be observed with the collisional brightening. Type Ibc supernova light curves are more likely to be brightened by the collision (0.53%) because of the high fraction of the progenitors being in binary systems and their proximity to the companion stars. Type II and IIb supernova light curves are less affected (~1e-3% and ~1e-2%, respectively). Although the early, slow light-curve declines of some Type IIb and Ibc supernovae are argued to be caused by the collision with the companion star (e.g. SN 2008D), the small expected fraction, as well as the unrealistically small separation required, disfavour the argument. The future transient survey by the Large Synoptic Survey Telescope is expected to detect ~10 Type Ibc supernovae with the early collisional brightening per year, and they will be able to provide information on supernova progenitors in binary systems.
Peculiar velocity measurements are the only tool available in the low-redshift Universe for mapping the large-scale distribution of matter and can thus be used to constrain cosmology. Using redshifts from the 2M++ redshift compilation, we reconstruct the density of galaxies within 200 Mpc/h, allowing for the first time good sampling of important superclusters such as the Shapley Concentration. We compare the predicted peculiar velocities from 2M++ to Tully-Fisher and SNe peculiar velocities. We find a value of $\beta^* \equiv \Omega_{\rm{m}}^{0.55}/b^* = 0.431 \pm 0.021$, suggesting $\Omega_{\rm{m}}^{0.55}\sigma_{\rm{8,lin}} = 0.401 \pm 0.024$, in good agreement with other probes. The predicted peculiar velocity of the Local Group arising from the 2M++ volume alone is $540 \pm 40$ km/s, towards $l = 268 \pm 4, b= 38 \pm 6$, only $10^\circ$ out of alignment with the Cosmic Microwave Background dipole. To account for velocity contributions arising from sources outside the 2M++ volume, we fit simultaneously for $\beta^*$ and an external bulk flow in our analysis. We find that an external bulk flow is preferred at the 5.1$\sigma$ level, and the best fit has a velocity of $159\pm23$ km/s towards $l=304 \pm 11, b = 6 \pm 13$. Finally, the predicted bulk flow of a 50 Mpc/h Gaussian-weighted volume centred on the Local Group is $230 \pm 30$ km/s, in the direction $l=293\pm 8, b = 14 \pm 10$, in agreement with predictions from $\Lambda$CDM.
Geological activity is thought to be important for the origin of life and for maintaining planetary habitability. We show that transient sulfate aerosols could be a signature of exoplanet volcanism, and therefore a geologically active world. A detection of transient aerosols, if linked to volcanism, could thus aid in habitability evaluations of the exoplanet. On Earth, subduction-induced explosive eruptions inject SO2 directly into the stratosphere, leading to the formation of sulfate aerosols with lifetimes of months to years. We demonstrate that the rapid increase and gradual decrease in sulfate aerosol loading associated with these eruptions may be detectable in transit transmission spectra with future large-aperture telescopes, such as the James Webb Space Telescope (JWST) and European Extremely-Large Telescope (E-ELT) for a planetary system at a distance of 10 pc, assuming an Earth-like atmosphere, bulk composition, and size. Specifically, we find that a S/N of 12.1 and 7.1 could be achieved with E-ELT (assuming photon-limited noise) for an Earth-analog orbiting a Sun-like star and M5V star, respectively, even without multiple transits binned together. We propose that the detection of this transient signal would strongly suggest an exoplanet volcanic eruption, if potential false positives such as dust storms or bolide impacts can be ruled out. Furthermore, because scenarios exist in which O2 can form abiotically in the absence of volcanic activity, a detection of transient aerosols that can be linked to volcanism, along with a detection of O2, would be a more robust biosignature than O2 alone.
We present a new likelihood-ratio ranking statistic for use in searches for gravitational waves from the inspiral and merger of compact object binaries. Expanding on previous work, the ranking statistic incorporates a model for the correlations in the signal-to-noise ratios with which signals will be seen in a network of ground-based antennas while retaining an algebraic procedure for mapping ranking statistic values to false-alarm probability. Additionally, the ranking statistic enables the implementation of a rigorous signal rate estimation technique. We implement the ranking statistic and demonstrate its use including signal rate estimation in an analysis of a simulated signal population in simulated noise.
We study the details of the DAMA/LIBRA results and compare those with the recent published DM Ice results of ICE Cube. In various recent papers, it was shown that the 40K peak on DAMA/LIBRA data leaves no room for a Dark Matter signal in the bulk of the data. Using Information Theory for the different types of detection environments, we show that annual variation calculations and the DM Ice data reinforce the claims that the DAMA/LIBRA detector is not observing Dark Matter WIMPs.
Proving the magnetic configuration of solar spicules has hitherto been difficult due to the lack of spatial resolution and image stability during off-limb ground-based observations. We report spectropolarimetric observations of spicules taken in the He I 1083 nm spectral region with the Tenerife Infrared Polarimeter II at the German Vacuum Tower Telescope of the Observatorio del Teide (Tenerife; Canary Islands; Spain). The data provide the variation with geometrical height of the Stokes I, Q, U, and V profiles whose encoded information allows the determination of the magnetic field vector by means of the HAZEL inversion code. The inferred results show that the average magnetic field strength at the base of solar spicules is about 80 gauss and then it decreases rapidly with height to about 30 gauss at a height of 3000 km above the visible solar surface. Moreover, the magnetic field vector is close to vertical at the base of the chromosphere and has mid inclinations (about 50 degree) above 2 Mm height.
The use of Reduced Proper Motion in identifying isolated white dwarfs has long been used as a proxy for the absolute magnitude in a population with known kinematics. This, however, introduces a proper motion detection limit on top of the existing photometric limit. How the survey volume is hampered by this extra parameter is discussed in Hambly et al. 2012. In this work, we discuss some robust outlier rejection methods in order to minimise the proper motion limit and hence maximise the survey volume. The generalised volume, corrected for the distance of the Sun from the Galactic Plane, is integrated explicitly.
The traditional Schmidt density estimator has been proven to be unbiased and effective in a magnitude limited sample. Previously, efforts have been made to generalise it for populations with non-uniform density and proper motion limited cases. This work shows that the then good assumptions for a proper motion limited sample are no longer sufficient to cope with modern data. Populations with larger differences in the kinematics as compared to the Local Standard of Rest are most severely affected. We show that this systematic bias can be removed by treating the discovery fraction inseparable from the generalised maximum volume integrand. The treatment can be applied to any proper motion limited sample with good knowledge of the kinematics. This work demonstrates the method through application to a mock catalogue of a white dwarf-only solar neighbourhood for various scenarios and compared against the traditional treatment using a survey with Pan-STARRS-like characteristics.
We present SymPix, a special-purpose spherical grid optimized for efficient sampling of rotationally invariant linear operators. This grid is conceptually similar to the Gauss-Legendre (GL) grid, aligning sample points with iso-latitude rings located on Legendre polynomial zeros. Unlike the GL grid, however, the number of grid points per ring varies as a function of latitude, avoiding expensive over-sampling near the poles and ensuring nearly equal sky area per grid point. The ratio between the number of grid points in two neighbouring rings is required to be a low-order rational number (3, 2, 1, 4/3, 5/4 or 6/5) to maintain a high degree of symmetries. Our main motivation for this grid is to solve linear systems using multi-grid methods, and to construct efficient preconditioners through pixel-space sampling of the linear operator in question. The GL grid is not suitable for these purposes due to its massive over-sampling near the poles, leading to nearly degenerate linear systems, while HEALPix, another commonly used spherical grid, exhibits few symmetries, and is therefore computationally inefficient for these purposes. As a benchmark and representative example, we compute a preconditioner for a linear system with both HEALPix and SymPix that involves the operator $D + B^T N^{-1} B$, where $B$ and $D$ may be described as both local and rotationally invariant operators, and $N$ is diagonal in pixel domain. For a bandwidth limit of $\ell_\text{max}=3000$, we find that SymPix, due to its higher number of internal symmetries, yields average speed-ups of 360 and 23 for $B^T N^{-1} B$ and $D$, respectively, relative to HEALPix.
High-precision frequencies of acoustic modes in red giant stars are now available thanks to the long observing length and high-quality of the light curves provided by the NASA Kepler mission, thus allowing to probe the interior of evolved cool low-mass stars with unprecedented level of detail. We characterize the acoustic signature of the helium second ionization zone in a sample of 18 low-mass low-luminosity red giants by exploiting new mode frequency measurements derived from more than four years of Kepler observations. We analyze the second frequency differences of radial acoustic modes in all the stars of the sample by using the Bayesian code Diamonds. We find clear acoustic glitches due to the signature of helium second ionization in all the stars of the sample. We measure the acoustic depth and the characteristic width of the acoustic glitches with a precision level on average around $\sim$2% and $\sim$8%, respectively. We find good agreement with theoretical predictions and existing measurements from the literature. Lastly, we derive the amplitude of the glitch signal at $\nu_\mathrm{max}$ for the second differences and for the frequencies with an average precision of $\sim$6%, obtaining values in the range 0.14-0.24 $\mu$Hz, and 0.08-0.33 $\mu$Hz, respectively, which can be used to investigate the helium abundance in the stars.
High-resolution (0.5 arcsec) CO(2-1) observations performed with the Atacama Large Millimetre/submillimetre Array have been used to trace the kinematics of the molecular gas in the Seyfert 2 galaxy{IC~5063}. Although one of the most radio-loud Seyfert galaxy, IC~5063 is a relatively weak radio source (P_1.4GHz = 3 x 10^23 W Hz^-1). The data reveal that the kinematics of the gas is very complex. A fast outflow of molecular gas extends along the entire radio jet (~ 1 kpc), with the highest outflow velocities about 0.5 kpc from the nucleus, at the location of the brighter hot-spot in the W lobe. All the observed characteristics can be described by a scenario of a radio plasma jet expanding into a clumpy medium, interacting directly with the clouds and inflating a cocoon that drives a lateral outflow into the interstellar medium. This suggests that most of the observed cold molecular outflow is due to fast cooling of the gas after the passage of a shock and that it is the end product of the cooling process.
Most type-Ic core-collapse supernovae (CCSNe) produce $^{56}$Ni and neutron stars (NSs) or black holes (BHs). The dipole radiation of nascent NSs has usually been neglected in explaining supernovae (SNe) with peak absolute magnitude $M_{\rm peak}$ in any band are $\gtrsim -19.5$~mag, while the $^{56}$Ni can be neglected in fitting most type-Ic superluminous supernovae (SLSNe Ic) whose $M_{\rm peak}$ in any band are $\lesssim -21$~mag, since the luminosity from a magnetar (highly magnetized NS) can outshine that from a moderate amount of $^{56}$Ni. For luminous SNe Ic with $-21 \lesssim M_{\rm peak}\lesssim -19.5$~mag, however, both contributions from $^{56}$Ni and NSs cannot be neglected without serious modeling, since they are not SLSNe and the $^{56}$Ni mass could be up to $\sim 0.5 M_{\odot}$. In this paper we propose a unified model that contain contributions from both $^{56}$Ni and a nascent NS. We select three luminous SNe Ic-BL, SN~2010ay, SN~2006nx, and SN~14475, and show that, if these SNe are powered by $^{56}$Ni, the ratio of $M_{\rm Ni}$ to $M_{\rm ej}$ are unrealistic. Alternatively, we invoke the magnetar model and the hybrid ($^{56}$Ni + NS) model and find that they can fit the observations, indicating that our models are valid and necessary for luminous SNe Ic. Owing to the lack of late-time photometric data, we cannot break the parameter degeneracy and thus distinguish among the model parameters, but we can expect that future multi-epoch observations of luminous SNe can provide stringent constraints on $^{56}$Ni yields and the parameters of putative magnetars.
We investigate the rotational evolution of young stars through Monte Carlo simulations. We simulate 280,000 stars, each of which is assigned a mass, a rotational period, and a mass accretion rate. The mass accretion rate depends on mass and time, following power-laws indices 1.4 and -1.5, respectively. A mass-dependent accretion threshold is defined below which a star is considered as diskless, which results in a distribution of disk lifetimes that matches observations. Stars are evolved at constant angular spin rate while accreting and at constant angular momentum when they become diskless. We recover the bimodal period distribution seen in several young clusters. The short period peak consists mostly of diskless stars and the long period one is mainly populated by accreting stars. Both distributions present a long tail towards long periods and a population of slowly rotating diskless stars is observed at all ages. We reproduce the observed correlations between disk fraction and spin rate, as well as between IR excess and rotational period. The period-mass relation we derive from the simulations exhibits the same global trend as observed in young clusters only if we release the disk locking assumption for the lowest mass stars. We find that the time evolution of median specific angular momentum follows a power law index of -0.65 for accreting stars and of -0.53 for diskless stars, a shallower slope that results from a wide distribution of disk lifetimes. Using observationally-documented distributions of disk lifetimes, mass accretion rates, and initial rotation periods, and evolving an initial population from 1 to 12 Myr, we reproduce the main characteristics of pre-main sequence angular momentum evolution, which supports the disk locking hypothesis. (abridged)
Relations between star formation rates along the spiral arms and the velocities of gas inflow into the arms in grand-design galaxy NGC 628 were studied. We found that the radial distribution of average star formation rate in individual star formation regions in regular spiral arms correlates with the velocity of gas inflow into the spiral arms. Both distributions have maxima at a galactocentric distance of 4.5-5 kpc. There are no correlations between the radial distributions of average star formation rate in star formation regions in spiral arms and outside spiral arms in the main disc. We also did not find a correlation between the radial distribution of average star formation rate in star formation regions in spiral arms and HI column density.
Solar flares and CMEs have a broad range of magnitudes. This review discusses the possibility of "extreme events," defined as those with magnitudes greater than have been seen in the existing historical record. For most quantitative measures, this direct information does not extend more than a century and a half into the recent past. The magnitude distributions (occurrence frequencies) of solar events (flares/CMEs) typically decrease with the parameter measured or inferred (peak flux, mass, energy etc. Flare radiation fluxes tend to follow a power law slightly flatter than $S^{-2}$, where S represents a peak flux; solar particle events (SPEs) follow a still flatter power law up to a limiting magnitude, and then appear to roll over to a steeper distribution, which may take an exponential form or follow a broken power law. This inference comes from the terrestrial $^{14}$C record and from the depth dependence of various radioisotope proxies in the lunar regolith and in meteorites. Recently major new observational results have impacted our use of the relatively limited historical record in new ways: the detection of actual events in the $^{14}$C tree-ring records, and the systematic observations of flares and "superflares" by the Kepler spacecraft. I discuss how these new findings may affect our understanding of the distribution function expected for extreme solar events.
Excess emission over expected diffuse astrophysical backgrounds in the direction of the Galactic center region has been claimed at various wavelengths, from radio to gamma rays. Among particle models advocated to explain such observations, several invoke interactions between dark matter particles and ordinary matter, such as cosmic rays, interstellar gas or free electrons. Depending on the specific interstellar matter particles' species and energy, such models predict distinct morphological features. In this study we make detailed predictions for the morphology of models where the relevant electromagnetic emission is proportional to the product of the dark matter density profile and the density of interstellar matter or cosmic rays. We compare the predicted latitudinal and longitudinal distributions with observations, and provide the associated set of relevant spatial templates.
Spiral arms that emerge from the ends of a galactic bar are important in interpreting observations of our and external galaxies. It is therefore important to understand the physical mechanism that causes them. We find that these spiral arms can be understood as kinematic density waves generated by librations around underlying ballistic closed orbits. This is even true in the case of a strong bar, provided the librations are around the appropriate closed orbits and not around the circular orbits that form the basis of the epicycle approximation. An important consequence is that it is a potential's orbital structure that determines whether a bar should be classified as weak or strong, and not crude estimates of the potential's deviation from axisymmetry.
We present a new method to constrain the grain size in protoplanetary disks with polarization observations at millimeter wavelengths. If dust grains are grown to the size comparable to the wavelengths, the dust grains are expected to have a large scattering opacity and thus the continuum emission is expected to be polarized due to self-scattering. We perform 3D radiative transfer calculations to estimate the polarization degree for the protoplanetary disks having radial Gaussian-like dust surface density distributions, which have been recently discovered. The maximum grain size is set to be $100 {\rm~\mu m}$ and the observing wavelength to be 870 ${\rm \mu m}$. We find that the polarization degree is as high as 2.5% with a subarcsec spatial resolution, which is likely to be detected with near-future ALMA observations. The emission is polarized due to scattering of anisotropic continuum emission. The map of the polarization degree shows a double peaked distribution and the polarization vectors are in the radial direction in the inner ring and in the azimuthal direction in the outer ring. We also find the wavelength dependence of the polarization degree: the polarized emission is strongest if dust grains have a maximum size of $a_{\rm max}\sim\lambda/2\pi$, where $\lambda$ is the observing wavelength. Hence, multi-wave and spatially resolved polarization observations toward protoplanetary disks enable us to put a constraint on the grain size. The constraint on the grain size from polarization observations is independent of or may be even stronger than that from the opacity index.
Planet formation is one explanation for the partial clearing of dust observed in the disks of some T Tauri stars. Indeed studies using state-of-the-art high angular resolution techniques have very recently begun to observe planetary companions in these so-called transitional disks. The goal of this work is to use spectra of the transitional disk object LkCa 15 obtained with X-Shooter on the Very Large Telescope to investigate the possibility of using spectro-astrometry to detect planetary companions to T Tauri stars. It is argued that an accreting planet should contribute to the total emission of accretion tracers such as H$\alpha$ and therefore planetary companions could be detected with spectro-astrometry in the same way as it has been used to detect stellar companions to young stars. A probable planetary-mass companion was recently detected in the disk of LkCa 15. Therefore, it is an ideal target for this pilot study. We studied several key accretion lines in the wavelength range 300 nm to 2.2 $\mu$m with spectro-astrometry. While no spectro-astrometric signal is measured for any emission lines the accuracy achieved in the technique is used to place an upper limit on the contribution of the planet to the flux of the H$\alpha$, Pa$\gamma$, and Pa$\beta$ lines. The derived upper limits on the flux allows an upper limit of the mass accretion rate, log($\dot{M}_{acc}$) = -8.9 to -9.3 for the mass of the companion between 6 M$_{Jup}$ and 15 M$_{Jup}$, respectively, to be estimated (with some assumptions).
The equation describing the secular diffusion of a self-gravitating collisionless system induced by an exterior perturbation is derived while assuming that the timescale corresponding to secular evolution is much larger than that corresponding to the natural frequencies of the system. Its two dimensional formulation for a tepid galactic disc is also derived using the epicyclic approximation. Its WKB limit is found while assuming that only tightly wound transient spirals are sustained by the disc. It yields a simple quadrature for the diffusion coefficients which provides a straightforward understanding of the loci of maximal diffusion within the disc.
The main orbital signatures of the secular evolution of an isolated self-gravitating stellar Mestel disc are recovered using a dressed Fokker-Planck formalism in angle-action variables. The shot-noise-driven formation of narrow ridges of resonant orbits is recovered in the WKB limit of tightly wound transient spirals, for a tepid Toomre-stable tapered disc. The relative effect of the bulge, the halo, the disc temperature and the spectral properties of the shot noise are investigated in turn. For such galactic discs all elements seem to impact the locus and direction of the ridge. For instance, when the halo mass is decreased, we observe a transition between a regime of heating in the inner regions of the disc through the inner Lindblad resonance to a regime of radial migration of quasi-circular orbits via the corotation resonance in the outer part of the disc. The dressed secular formalism captures both the nature of collisionless systems (via their natural frequencies and susceptibility), and their nurture via the structure of the external perturbing power spectrum. Hence it provides the ideal framework in which to study their long term evolution.
We consider a sample of ten GRBs with long lasting ($\gtrsim10^2\rm sec$) emission detected by Fermi/LAT and for which X-ray data around $1\,$day are also available. We assume that both the X-rays and the GeV emission are produced by electrons accelerated at the external shock, and show that the X-ray and the LAT fluxes lead to very different estimates of the initial kinetic energy of the blast wave. The energy estimated from LAT is on average $\sim50$ times larger than the one estimated from X-rays. We model the data (accounting also for optical detections around $1\,$day, if available) to unveil the reason for this discrepancy and find that good modelling within the external shock scenario is always possible and leads to two possibilities: either the X-ray emitting electrons (unlike the GeV emitting electrons) are in the slow cooling regime or ii) the X-ray synchrotron flux is strongly suppressed by Compton cooling, whereas, due to the Klein-Nishina suppression, this effect is much smaller at GeV energies. In both cases the X-ray flux is no longer a robust proxy for the blast wave kinetic energy. On average, both cases require weak magnetic fields ($10^{-6}\lesssim \epsilon_B \lesssim 10^{-3}$) and relatively large isotropic equivalent kinetic blast wave energies, in the range $10^{53}\rm erg$$<E_{0,kin}<10^{55}\rm erg$. These energies are larger than those estimated from the X-ray flux alone, and imply smaller inferred values of the prompt efficiency mechanism, reducing the efficiency requirements on the still uncertain mechanism responsible for prompt emission.
Fluctuations in a stellar system's gravitational field cause the orbits of stars to evolve. The resulting evolution of the system can be computed with the orbit-averaged Fokker-Planck equation once the diffusion tensor is known. We present the formalism that enables one to compute the diffusion tensor from a given source of noise in the gravitational field when the system's dynamical response to that noise is included. In the case of a cool stellar disc we are able to reduce the computation of the diffusion tensor to a one-dimensional integral. We implement this formula for a tapered Mestel disc that is exposed to shot noise and find that we are able to explain analytically the principal features of a numerical simulation of such a disc. In particular the formation of narrow ridges of enhanced density in action space is recovered. As the disc's value of Toomre's $Q$ is reduced and the disc becomes more responsive, there is a transition from a regime of heating in the inner regions of the disc through the inner Lindblad resonance to one of radial migration of near-circular orbits via the corotation resonance in the intermediate regions of the disc. The formalism developed here provides the ideal framework in which to study the long-term evolution of all kinds of stellar discs.
We examine the heating of the intra-cluster medium (ICM) of cooling flow clusters of galaxies by jet-inflated bubbles and conclude that mixing of hot bubble gas with the ICM is the dominate heating process. We use the PLUTO hydrodynamical code in full 3D to properly account for the inflation of the bubbles and to the multiple vortices induced by the jets and bubbles. The vortices mix some hot shocked jet gas with the ICM. For the parameters used the mixing process accounts for approximately 80% of the energy transferred from the jets to the ICM. Only about 20% of the transferred energy is channelled to the kinetic energy of the ICM. Part of this develops as ICM turbulence. We conclude that turbulent heating plays a smaller role than mixing. Heating by shocks is less efficient even.
A series of optical and one near-infrared nebular spectra covering the first year of the Type Ia supernova SN 2011fe are presented and modelled. The density profile that proved best for the early optical/ultraviolet spectra, "rho-11fe", was extended to lower velocities to include the regions that emit at nebular epochs. Model rho-11fe is intermediate between the fast deflagration model W7 and a low-energy delayed-detonation. Good fits to the nebular spectra are obtained if the innermost ejecta are dominated by neutron-rich, stable Fe-group species, which contribute to cooling but not to heating. The correct thermal balance can thus be reached for the strongest [FeII] and [FeIII] lines to be reproduced with the observed ratio. The 56Ni mass thus obtained is 0.47 +/- 0.05 Mo. The bulk of 56Ni has an outermost velocity of ~8500 km/s. The mass of stable iron is 0.23 +/- 0.03 Mo. Stable Ni has low abundance, ~10^{-2} Mo. This is sufficient to reproduce an observed emission line near 7400 A. A sub-Chandrasekhar explosion model with mass 1.02 Mo and no central stable Fe does not reproduce the observed line ratios. A mock model where neutron-rich Fe-group species are located above 56Ni following recent suggestions is also shown to yield spectra that are less compatible with the observations. The densities and abundances in the inner layers obtained from the nebular analysis, combined with those of the outer layers previously obtained, are used to compute a synthetic bolometric light curve, which compares favourably with the light curve of SN 2011fe.
The Blandford-Znajek (BZ) mechanism describes a process extracting rotation energy from a spinning black hole (BH) via magnetic field lines penetrating the event horizon of central BH. We report, for the first time, a general analytic approach to study force-free jets launched by the BZ mechanism, and its three immediate applications: (1) we present a high-order split monopole perturbation solution to the BZ mechanism, which accurately pins down the energy extraction rate $\dot E$ and well describes the structure of BH magnetosphere for all range of BH spins ($0\leq a\leq 1$); (2) the approach yields an exact constraint for the monopole field configuration in the Kerr spacetime, $I = \Omega (1-A_\phi^2)$, where $A_\phi$ is the $\phi-$component of electromagnetic field potential, $\Omega$ is the angular velocity of magnetic field lines and $I$ is the poloidal electric current. The constraint is of particular importance to benchmark the accuracy of numerical simulations; (3) we prove the uniqueness of solutions to BZ mechanism, i.e., whenever there is a Schwarzschild metric solution there exists a unique Kerr metric solution approaching it at far distance.
We present a new method for deriving the stellar birth function (SBF) of resolved stellar populations. The SBF (stars born per unit mass, time, and metallicity) is the combination of the initial mass function (IMF), the star-formation history (SFH), and the metallicity distribution function (MDF). The framework of our analysis is that of Poisson Point Processes (PPPs), a class of statistical models suitable when dealing with points (stars) in a multidimensional space (the measurement space of multiple photometric bands). The theory of PPPs easily accommodates the modeling of measurement errors as well as that of incompleteness. Compared to most of the tools used to study resolved stellar populations, our method avoids binning stars in the color-magnitude diagram and uses the entirety of the information (i.e., the whole likelihood function) for each data point; the proper combination of the individual likelihoods allows the computation of the posterior probability for the global population parameters. This includes unknowns such as the IMF slope and combination of SFH and MDF, which are rarely solved for simultaneously in the literature, however entangled and correlated they might be. Our method also allows proper inclusion of nuisance parameters, such as distance and extinction distributions. The aim of this paper, is to assess the validity of this new approach under a range of assumptions, using only simulated data. Forthcoming work will show applications to real data. Although it has a broad scope of possible applications, we have developed this method to study multi-band HST observations of the Milky Way Bulge. Therefore we will focus on simulations with characteristics similar to those of the Galactic Bulge.
We developed a model for young giant exoplanets (Exoplanet
Radiative-convective Equilibrium Model or Exo-REM). Input parameters are
planet's surface gravity (g), effective temperature (Teff ) and elemental
composition. Under the additional assumption of thermochemical equilibrium, the
model predicts the equilibrium temperature profile and mixing ratio profiles of
the most important gases. Opacity sources include the H$_2$-He
collision-induced absorption and molecular lines from H$_2$O, CO, CH$_4$
(updated with the Exomol linelist), NH$_3$, VO, TiO, Na and K. Absorption by
iron and silicate cloud particles is added above the expected condensation
levels with a fixed scale height and a given optical depth at some reference
wavelength. Scattering was not included at this stage.
We applied Exo-REM to photometric and spectral observations of the planet
beta Pictoris b obtained in a series of near IR filters. We derived Teff = 1550
$\pm$ 150 K, log(g) = 3.5 $\pm$ 1, and a radius R = 1.76 $\pm$ 0.24 R Jup
(2-$\sigma$ error bars). These values are comparable to those found in the
literature, although with more conservative error bars, but consistent with the
model accuracy. We finally investigated the precision to which the above
parameters can be constrained from SPHERE measurements using different sets of
near IR filters as well as near low resolution spectroscopy.
A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.
Planet formation studies are often focused on solar-type stars, implicitly considering our Sun as reference point. This approach overlooks, however, that Herbig Ae/Be stars are in some sense much better targets to study planet formation processes empirically, with their disks generally being larger, brighter and simply easier to observe across a large wavelength range. In addition, massive gas giant planets have been found on wide orbits around early type stars, triggering the question if these objects did indeed form there and, if so, by what process. In the following I briefly review what we currently know about the occurrence rate of planets around intermediate mass stars, before discussing recent results from Herbig Ae/Be stars in the context of planet formation. The main emphasis is put on spatially resolved polarized light images of potentially planet forming disks and how these images - in combination with other data - can be used to empirically constrain (parts of) the planet formation process. Of particular interest are two objects, HD100546 and HD169142, where, in addition to intriguing morphological structures in the disks, direct observational evidence for (very) young planets has been reported. I conclude with an outlook, what further progress we can expect in the very near future with the next generation of high-contrast imagers at 8-m class telescopes and their synergies with ALMA.
A report is made on the luminosity and pulse-period evolution of the Be binary X-ray pulsar, GX 304$-$1, during a series of outbursts from 2009 to 2013 observed by MAXI/GSC, RXTE/PCA, and Fermi/GBM. In total, twelve outbursts repeated by $\sim$ 132.2 days were observed, which is consistent with the X-ray periodicity of this object observed in the 1970s. These 12 outbursts, together with those in the 1970s, were found to all recur with a well defined period of 132.189$\pm$0.02 d, which can be identified with the orbital period. The pulse period at $\sim 275$ s, obtained from the RXTE/PCA and Fermi/GBM data, apparently exhibited a periodic modulation synchronized with the outburst period, suggesting the pulsar orbital motion, which is superposed on a secular spin-up trend throughout the entire active phase. The observed pulse-period changes were successfully represented by a model composed of the binary orbital modulation and pulsar spin up caused by mass accretion through an accretion disk. The orbital elements obtained from the best-fit model, including the projected orbital semi-major axis $a_{\rm x}\sin i \simeq 500-600$ light-s and an eccentricity $e \simeq 0.5$, are typical of Be binary X-ray pulsars.
The Planck satellite detectors are calibrated in the 2015 release using the "orbital dipole", which is the time-dependent dipole generated by the Doppler effect due to the motion of the satellite around the Sun. Such an effect has also relativistic time-dependent corrections of relative magnitude 10^(-3), due to coupling with the "solar dipole" (the motion of the Sun compared to the CMB rest frame), which are included in the data calibration by the Planck collaboration. We point out that such corrections are subject to a frequency-dependent multiplicative factor. This factor differs from unity especially at the highest frequencies, relevant for the HFI instrument. Since currently Planck calibration errors are dominated by systematics, to the point that polarization data is currently unreliable at large scales, such a correction can in principle be highly relevant for future data releases.
We report the results of two multi-chord stellar occultations by the dwarf planet (1) Ceres that were observed from Brazil on 2010 August 17, and from the USA on 2013 October 25. Four positive detections were obtained for the 2010 occultation, and nine for the 2013 occultation. Elliptical models were adjusted to the observed chords to obtain Ceres' size and shape. Two limb fitting solutions were studied for each event. The first one is a nominal solution with an indeterminate polar aspect angle. The second one was constrained by the pole coordinates as given by Drummond et al. Assuming a Maclaurin spheroid, we determine an equatorial diameter of 972 $\pm$ 6 km and an apparent oblateness of 0.08 $\pm$ 0.03 as our best solution. These results are compared to all available size and shape determinations for Ceres made so far, and shall be confirmed by the NASA's Dawn space mission.
Aims. Every 5.5 years eta Car's light curve and spectrum change remarkably across all observed wavelength bands. We compare the recent spectroscopic event in mid-2014 to the events in 2003 and 2009 and investigate long-term trends. Methods. Eta Car was observed with HST STIS, VLT UVES, and CTIO 1.5m CHIRON for a period of more than two years in 2012-2015. Archival observations with these instruments cover three orbital cycles. Results. Important spectroscopic diagnostics show significant changes in 2014 compared to previous events. While the timing of the first HeII 4686 flash was remarkably similar to previous events, the HeII equivalent widths were slightly larger and the line flux increased compared to 2003. The second HeII peak occurred at about the same phase as in 2009, but was stronger. The HeI line flux grew in 2009-2014 compared to 1998-2003. On the other hand, Halpha and FeII lines show the smallest emission strengths ever observed. Conclusions. The basic character of the spectroscopic events has changed in the past 2-3 cycles; ionizing UV radiation dramatically weakened during each pre-2014 event but not in 2014. The strengthening of HeI emission and the weakening of the lower-excitation wind features in our direct line of sight implies a substantial change in the physical parameters of the emitting regions. The polar spectrum at FOS4 shows less changes in the broad wind emission lines, which may be explained by the latitude-dependent wind structure of eta Car. The quick and strong recovery of the HeII emission in 2014 supports a scenario, in which the wind-wind shock may not have completely collapsed as was proposed for previous events. All this may be the consequence of just one elementary change, namely a strong decrease in the primary's mass-loss rate.
We present the first fully calibrated H$_2$, 1-0 S(1) image of the entire 30 Doradus nebula. The observations were conducted using the NOAO Extremely Wide-Field Infrared Imager on the CTIO 4-meter Blanco Telescope. Together with a NEWFIRM Br$\gamma$ image of 30 Doradus, our data reveal the morphologies of the warm molecular gas and ionized gas in 30 Doradus. The brightest H$_2$-emitting area, which extends from the northeast to the southwest of R136, is a photodissociation region viewed face-on, while many clumps and pillar features located at the outer shells of 30 Doradus are photodissociation regions viewed edge-on. Based on the morphologies of H$_2$, Br$\gamma$, $^{12}$CO, and 8$\mu$m emission, the H$_2$ to Br$\gamma$ line ratio and Cloudy models, we find that the H$_2$ emission is formed inside the photodissociation regions of 30 Doradus, 2 - 3 pc to the ionization front of the HII region, in a relatively low-density environment $<$ 10$^4$ cm$^{-3}$. Comparisons with Br$\gamma$, 8$\mu$m, and CO emission indicate that H$_2$ emission is due to fluorescence, and provide no evidence for shock excited emission of this line.
Millisecond pulsars occur abundantly in globular clusters. They are expected to be responsible for several spectral components in the radio through gamma-ray waveband (e.g., involving synchrotron and inverse Compton emission), as have been seen by Radio Telescope Effelsberg, Chandra X-ray Observatory, Fermi Large Area Telescope, and the High Energy Stereoscopic System (H.E.S.S.) in the case of Terzan 5 (with fewer spectral components seen for other globular clusters). H.E.S.S. has recently performed a stacking analysis involving 15 non-detected globular clusters and obtained quite constraining average flux upper limits above 230 GeV. We present a model that assumes millisecond pulsars as sources of relativistic particles and predicts multi-wavelength emission from globular clusters. We apply this model to the population of clusters mentioned above to predict the average spectrum and compare this to the H.E.S.S. upper limits. Such comparison allows us to test whether the model is viable, leading to possible constraints on various average cluster parameters within this framework.
O-C diagram is a useful technique to analyse the period changes of a pulsator by using the maximum (or minimum) value points which have been obtained from the historical data. But if an object is a double-mode or multi-mode pulsator, the extreme value points are the results of all the modes other than just the fundamental mode. We discussed these situations and give out some criteria to judge whether the O-C diagram can be used in these situations.
We report the identification of elongated (triaxial or prolate) galaxies in cosmological simulations at $z\simeq2$. These are preferentially low-mass galaxies ($M_s \le 10^{9.5} \ M_\odot$), residing in dark-matter (DM) haloes with strongly elongated inner parts, a common feature of high-redshift DM haloes in the $\Lambda$CDM cosmology. Feedback slows formation of stars at the centres of these halos, so that a dominant and prolate DM distribution gives rise to galaxies elongated along the DM major axis. As galaxies grow in stellar mass, stars dominate the total mass within the galaxy half-mass radius, making stars and DM rounder and more oblate. A large population of elongated galaxies produces a very asymmetric distribution of projected axis ratios, as observed in high-z galaxy surveys. This indicates that the majority of the galaxies at high redshifts are not discs or spheroids but rather galaxies with elongated morphologies.
Recent high-resolution observations have shown stellar winds to harbour complexities which strongly deviate from spherical symmetry, generally assumed as standard wind model. One such morphology is the archimedean spiral, generally believed to be formed by binary interactions, which has been directly observed in multiple sources. We seek to investigate the manifestation in the observables of spiral structures embedded in the spherical outflows of cool stars. We aim to provide an intuitive bedrock with which upcoming ALMA data can be compared and interpreted. By means of an extended parameter study, we model rotational CO emission from the stellar outflow of asymptotic giant branch stars. To this end, we develop a simplified analytical parametrised description of a 3D spiral structure. This model is embedded into a spherical wind, and fed into the 3D radiative transfer code LIME, which produces 3D intensity maps throughout velocity space. Subsequently, we investigate the spectral signature of rotational transitions of CO of the models, as well as the spatial aspect of this emission by means of wide-slit PV diagrams. Additionally, the potential for misinterpretation of the 3D data in a 1D context is quantified. Finally, we simulate ALMA observations to explore the impact of interefrometric noise and artifacts on the emission signatures. The spectral signatures of the CO rotational transition v=0 J=3-2 are very efficient at concealing the dual nature of the outflow. Only a select few parameter combinations allow for the spectral lines to disclose the presence of the spiral structure. The inability to disentangle the spiral from the spherical signal can result in an incorrect interpretation in a 1D context. Consequently, erroneous mass loss rates would be calculated...
(Abridged) Neutrino interactions beyond the standard model may affect the cosmological evolution and can be constrained through observations. We consider the possibility that neutrinos possess secret scalar or pseudoscalar interactions mediated by the Nambu-Goldstone boson of a still unknown spontaneously broken global $U(1)$ symmetry, as in, e.g. , Majoron models. In such scenarios, neutrinos still decouple at $T\simeq 1$ MeV, but become tightly coupled again ('recouple') at later stages of the cosmological evolution. We use available observations of CMB anisotropies, including Planck 2013 and the joint BICEP2/Planck 2015 data, to derive constraints on the quantity $\gamma_{\nu \nu}^4$, parameterizing the neutrino collision rate due to (pseudo)scalar interactions. We consider both a minimal extension of the standard $\Lambda$CDM model, and scenarios with extra relativistic species or non-vanishing tensors. We find a typical constraint $\gamma_{\nu \nu}^4 < 0.9\times 10^{-27}$ (95\% C.L.), implying an upper limit on the redshift $z_{rec}$ of neutrino recoupling $< 8500$. In the framework of Majoron models, the upper limit on $\gamma_{\nu \nu}$ roughly translates on a constraint $g < 8.2\times 10^{-7}$ on the Majoron-neutrino coupling constant $g$. In general, the data show a weak ($\sim 1\sigma$) but intriguing preference for non-zero values of $\gamma_{\nu \nu}^4$, with best fits in the range $\gamma_{\nu \nu}^4 = (0.15 - 0.35)\times 10^{-27}$, depending on the particular dataset. This is more evident when either observations from ACT and SPT are included, or the possibility of non-vanishing tensor modes is considered. In particular, for the minimal model $\Lambda$CDM +$\gamma_{\nu \nu}$ and including the Planck 2013, ACT and SPT data, we report $\gamma_{\nu \nu}^4=( 0.45^{+0.15}_{-0.38} )\times10^{-27}$ ($200 < z_{rec} < 5700$) at 68\% confidence level.
In dense parts of interstellar clouds (> 10^5 cm^-3), dust & gas are expected to be in thermal equilibrium, being coupled via collisions. However, previous studies have shown that the temperatures of the dust & gas may remain decoupled even at higher densities. We study in detail the temperatures of dust & gas in the photon-dominated region S 140, especially around the deeply embedded infrared sources IRS 1-3 and at the ionization front. We derive the dust temperature and column density by combining Herschel PACS continuum observations with SOFIA observations at 37 $\mu$m and SCUBA at 450 $\mu$m. We model these observations using greybody fits and the DUSTY radiative transfer code. For the gas part we use RADEX to model the CO 1-0, CO 2-1, 13CO 1-0 and C18O 1-0 emission lines mapped with the IRAM-30m over a 4' field. Around IRS 1-3, we use HIFI observations of single-points and cuts in CO 9-8, 13CO 10-9 and C18O 9-8 to constrain the amount of warm gas, using the best fitting dust model derived with DUSTY as input to the non-local radiative transfer model RATRAN. We find that the gas temperature around the infrared sources varies between 35 and 55K and that the gas is systematically warmer than the dust by ~5-15K despite the high gas density. In addition we observe an increase of the gas temperature from 30-35K in the surrounding up to 40-45K towards the ionization front, most likely due to the UV radiation from the external star. Furthermore, detailed models of the temperature structure close to IRS 1 show that the gas is warmer and/or denser than what we model. Finally, modelling of the dust emission from the sub-mm peak SMM 1 constrains its luminosity to a few ~10^2 Lo. We conclude that the gas heating in the S 140 region is very efficient even at high densities, most likely due to the deep UV penetration from the embedded sources in a clumpy medium and/or oblique shocks.
The secular evolution of an infinitely thin tepid isolated galactic disc made
of a finite number of particles is described using the inhomogeneous
Balescu-Lenard equation. Assuming that only tightly wound transient spirals are
present in the disc, a WKB approximation provides a simple and tractable
quadrature for the corresponding drift and diffusion coefficients. It provides
insight into the physical processes at work during the secular diffusion of a
self-gravitating discrete disc and makes quantitative predictions on the
initial variations of the distribution function in action space.
When applied to the secular evolution of an isolated stationary
self-gravitating Mestel disc, this formalism predicts initially the importance
of the corotation resonance in the inner regions of the disc leading to a
regime involving radial migration and heating. It predicts in particular the
formation of a "ridge like" feature in action space, in agreement with
simulations, but over-estimates the timescale involved in its appearance. Swing
amplification is likely to resolve this discrepancy.
In astrophysics, the inhomogeneous Balescu-Lenard equation and its WKB limit
may also describe the secular diffusion of giant molecular clouds in galactic
discs, the secular migration and segregation of planetesimals in
proto-planetary discs, or even the long-term evolution of population of stars
within the Galactic center.
The discovery of millisecond pulsars switching between states powered either by the rotation of their magnetic field or by the accretion of matter, has recently proved the tight link shared by millisecond radio pulsars and neutron stars in low-mass X-ray binaries. Transitional millisecond pulsars also show an enigmatic intermediate state in which the neutron star is surrounded by an accretion disk, it emits coherent X-ray pulsations, but is sub-luminous in X-rays with respect to accreting neutron stars, and is brighter in gamma-rays than millisecond pulsars in the rotation-powered state. Here, we model the X-ray and gamma-ray emission observed from PSR J1023+0038 in such a state based on the assumption that most of the disk in-flow is propelled away by the rapidly rotating neutron star magnetosphere, and that electrons can be accelerated to energies of a few GeV at the turbulent disk-magnetosphere boundary. We show that the synchrotron and self-synchrotron Compton emission coming from such a region, together with the hard disk emission typical of low states of accreting compact objects, is able to explain the radiation observed in the X-ray and gamma-ray band. The average emission observed from PSR J1023+0038 is modelled by a disk in-flow with a rate of $(1-3)\times10^{-11} M_{\odot}/yr$, truncated at a radius ranging between 30 and 45 km, compatible with the hypothesis of a propelling magnetosphere. We compare the results we obtained with models that rather assume that a rotation-powered pulsar is turned on, showing how the spin down power released in similar scenarios is hardly able to account for the magnitude of the observed emission.
In the last decade ground-based Imaging Atmospheric Cherenkov Telescopes have discovered roughly 30 pulsar wind nebulae at energies above 100 GeV. We present first results from a leptonic emission code that models the spectral energy density of a pulsar wind nebula by solving the Fokker-Planck transport equation and calculating inverse Compton and synchrotron emissivities. Although models such as these have been developed before, most of them model the geometry of a pulsar wind nebula as that of a single sphere. We have created a time-dependent, multi-zone model to investigate changes in the particle spectrum as the particles diffuse through the pulsar wind nebula, as well as predict the radiation spectrum at different positions in the nebula. We calibrate our new model against a more basic previous model and fit the observed spectrum of G0.9+0.1, incorporating data from the High Energy Stereoscopic System as well as radio and X-ray experiments.
The accreting millisecond X-ray pulsar SAX J1808.4--3658 shows peculiar low luminosity states known as "reflares" after the end of the main outburst. During this phase the X-ray luminosity of the source varies by up to three orders of magnitude in less than 1-2 days. The lowest X-ray luminosity observed reaches a value of ~1e32 erg/s, only a factor of a few brighter than its typical quiescent level. We investigate the 2008 and 2005 reflaring state of SAX J1808.4-3658 to determine whether there is any evidence for a change in the accretion flow with respect to the main outburst. We perform a multiwavelength photometric and spectral study of the 2005 and 2008 reflares with data collected during an observational campaign covering the near-infrared, optical, ultra-violet and X-ray band. We find that the NIR/optical/UV emission, expected to some from the outer accretion disk shows variations in luminosity which are 1--2 orders of magnitude shallower than in X-rays. The X-ray spectral state observed during the reflares does not change substantially with X-ray luminosity indicating a rather stable configuration of the accretion flow. We investigate the most likely configuration of the innermost regions of the accretion flow and we infer an accretion disk truncated at or near the co-rotation radius. We interpret these findings as due to either a strong outflow (due to a propeller effect) or a trapped disk (with limited/no outflow) in the inner regions of the accretion flow.
We investigate features of the lateral distribution function (LDF) of the radio signal emitted by cosmic ray air-showers with primary energies $> 0.1$~EeV and its connection to air-shower parameters such as energy and shower maximum using CoREAS simulations made for the configuration of the Tunka-Rex antenna array. Taking into account all significant contributions to the total radio emission, such as by the geomagnetic effect, the charge excess, and the atmospheric refraction we parameterize the radio LDF. This parameterization is two-dimensional and has several free parameters. The large number of free parameters is not suitable for experiments of sparse arrays operating at low SNR (signal-to-noise ratios). Thus, exploiting symmetries, we decrease the number of free parameters and reduce the LDF to a simple one-dimensional function. The remaining parameters can be fit with a small number of points, i.e. as few as the signal from three antennas above detection threshold. Finally, we present a method for the reconstruction of air-shower parameters, in particular, energy and $X_{\mathrm{max}}$ (shower maximum), which can be reached with a theoretical accuracy of better than 15\% and 30~g/cm$^2$, respectively.
We present optical integral field spectroscopy of the circum-nuclear gas of the Seyfert 2 galaxy NGC 1386. The data cover the central 7$^{\prime\prime} \times 9^{\prime\prime}$ (530 $\times$ 680 pc) at a spatial resolution of 0.9" (68 pc), and the spectral range 5700-7000 \AA\ at a resolution of 66 km s$^{-1}$. The line emission is dominated by a bright central component, with two lobes extending $\approx$ 3$^{\prime\prime}$ north and south of the nucleus. We identify three main kinematic components. The first has low velocity dispersion ($\bar \sigma \approx $ 90 km s$^{-1}$), extends over the whole field-of-view, and has a velocity field consistent with gas rotating in the galaxy disk. We interpret the lobes as resulting from photoionization of disk gas in regions where the AGN radiation cones intercept the disk. The second has higher velocity dispersion ($\bar \sigma \approx$ 200 km s$^{-1}$) and is observed in the inner 150 pc around the continuum peak. This component is double peaked, with redshifted and blueshifted components separated by $\approx$ 500 km s$^{-1}$. Together with previous HST imaging, these features suggest the presence of a bipolar outflow for which we estimate a mass outflow rate of $\mathrm{\dot M} \gtrsim $ 0.1 M$_{\odot}$ yr$^{-1}$. The third component is revealed by velocity residuals associated with enhanced velocity dispersion and suggests that outflow and/or rotation is occurring approximately in the equatorial plane of the torus. A second system of velocity residuals may indicate the presence of streaming motions along dusty spirals in the disk.
Continuous and precise space-based photometry has made it possible to measure the orbital frequency modulation of pulsating stars in binary systems with extremely high precision over long time spans. Frequency modulation caused by binary orbital motion manifests itself as a multiplet with equal spacing of the orbital frequency in the Fourier transform. The amplitudes and phases of the peaks in these multiplets reflect the orbital properties, hence the orbital parameters can be extracted by analysing such precise photometric data alone. We derive analytically the theoretical relations between the multiplet properties and the orbital parameters, and present a method for determining these parameters, including the eccentricity and the argument of periapsis, from a quintuplet or a higher order multiplet. This is achievable with the photometry alone, without spectroscopic radial velocity measurements. We apply this method to Kepler mission data of KIC8264492, KIC9651065, and KIC10990452, each of which is shown to have an eccentricity exceeding 0.5. Radial velocity curves are also derived from the Kepler photometric data. We demonstrate that the results are in good agreement with those obtained by another technique based on the analysis of the pulsation phases.
We present the latest developments to the phase modulation method for finding binaries among pulsating stars. We demonstrate how the orbital elements of a pulsating binary star can be obtained analytically, that is, without converting time delays to radial velocities by numerical differentiation. Using the time delays directly offers greater precision, and allows the parameters of much smaller orbits to be derived. The method is applied to KIC9651065, KIC10990452, and KIC8264492, and a set of the orbital parameters is obtained for each system. Radial velocity curves for these stars are deduced from the orbital elements thus obtained.
Using Hubble Space Telescope Advanced Camera for Surveys (HST/ACS) and Wide Field Camera 3 (WFC3) observations from the Panchromatic Hubble Andromeda Treasury (PHAT), we present new period-luminosity relations for Cepheid variables in M31. Cepheid from several ground-based studies are identified in the PHAT pho- tometry to derive new Period-Luminosity and Wesenheit Period-Luminosity relations in the NIR and visual filters. We derive a distance modulus to M31 of 24.51+/-0.08 in the IR bands and 24.32+/-0.09 in the visual bands, including the first PL relations in the F475W and F814W filters for M31. Our derived visual and IR distance moduli dis- agree at slightly more than a 1-{\sigma} level. Differences in the Period-Luminosity relations between ground-based and HST observations are investigated for a subset of Cepheids. We find a significant discrepancy between ground-based and HST Period-Luminosity relations with the same Cepheids, suggesting adverse effects from photometric contam- ination in ground-based Cepheid observations. Additionally, a statistically significant radial trend in the PL relation is found which does not appear to be explained by metallicity.
We observed the Fermi-discovered gamma-ray binary 1FGL J1018.6-5856 at 20 epochs over 50 days using the CHIRON spectrograph, obtaining spectra at R~25,000 covering 4090-8908A. The average spectrum confirms an O6 V((f)) spectral type and extinction E(B-V) = 1.35+/-0.04. Variable absorption line equivalent widths suggest substantial contamination by wind line features. The limited S/N ratio hindered accurate continuum definition and prevented measurement of a high quality radial velocity curve. Nevertheless, the best data indicate a radial velocity amplitude <40 km/s for the He II lines and substantially lower for H I. We argue that this indicates a most likely compact object mass <2.2Msun. While black hole solutions are not excluded, a neutron star source of the gamma-ray emission seems preferred.
Near infrared images from the COBE satellite presented the first clear evidence that our Milky Way galaxy contains a boxy shaped bulge. Recent years have witnessed a gradual paradigm shift in the formation and evolution of the Galactic bulge. Bulges were commonly believed to form in the dynamical violence of galaxy mergers. However, it has become increasingly clear that the main body of the Milky Way bulge is not a classical bulge made by previous major mergers, instead it appears to be a bar seen somewhat end-on. The Milky Way bar can form naturally from a precursor disk and thicken vertically by the internal firehose/buckling instability, giving rise to the boxy appearance. This picture is supported by many lines of evidence, including the asymmetric parallelogram shape, the strong cylindrical rotation (i.e., nearly constant rotation regardless of the height above the disk plane), the existence of an intriguing X-shaped structure in the bulge, and perhaps the metallicity gradients. We review the major theoretical models and techniques to understand the Milky Way bulge. Despite the progresses in recent theoretical attempts, a complete bulge formation model that explains the full kinematics and metallicity distribution is still not fully understood. Upcoming large surveys are expected to shed new light on the formation history of the Galactic bulge.
The apsidal precession frequency in a fixed gravitational potential increases with the radial range of the orbit (eccentricity). Although the frequency increase is modest it can have important implications for wave dynamics in galaxy discs, which have not been previously explored in detail. One of the most interesting consequences is that for a given pattern frequency, each Lindblad resonance does not exist in isolation, but rather is the parent of a continuous sequence of resonant radii, a Lindblad Zone, with each radius in this zone characterized by a specific eccentricity. In the epicyclic approximation the precession or epicyclic frequency does not depend on epicycle size, and this phenomenon is not captured. A better approximation for eccentric orbits is provided by p-ellipse curves (Struck 2006), which do exhibit this effect. Here the p-ellipse approximation and precession-eccentricity relation are used as tools for finding the resonant radii generated from various Lindblad parent resonances. Simple, idealized examples, in flat rotation curve and near solid-body discs, are used to show that ensembles of eccentric resonant orbits excited in Lindblad Zones can provide a backbone for generating a variety of (kinematic) bars and spiral waves. In cases balancing radius-dependent circular frequencies and eccentricity-dependent precession, a range of resonant orbits can maintain their form in the pattern frame, and do not wind up. Eccentric resonance orbits require a strong perturbation to excite them, and may be produced mostly in galaxy interactions or by strong internal disturbances.
Earth and orbital based radar observations of asteroids provide a unique opportunity to characterize surface roughness and the dielectric properties of their surfaces, as well as potentially explore some of their shallow subsurface physical properties. If the dielectric and topographic properties of asteroid's surfaces are defined, we can constrain their textural characteristics as well as potential subsurface volatile enrichment using the observed radar backscatter. To achieve this objective, we establish the first dielectric model of asteroid Vesta for the case of a dry, volatile-poor regolith -- employing an analogy to the dielectric properties of lunar regolith, and adjusted for the surface densities and temperatures deduced from Dawn's Visible and InfraRed mapping spectrometer (VIR). Our model suggests that the dielectric constant at the surface of Vesta is relatively constant, ranging from 2.0 to 2.1 from the night- to day-side of Vesta, while the loss tangent shows slight variation as a function of diurnal temperature, ranging from 0.011 to 0.014. We estimate the surface porosity to be ~55% in the upper meter of the regolith, as derived from VIR observations. This is ~20% higher than previous estimation of porosity derived from previous Earth-based X- and S-band radar observation. We suggest that the radar backscattering properties of asteroid Vesta will be mainly driven by the changes in surface roughness rather than potential dielectric variations in the upper regolith in the X- and S-band.
We present a novel determination of the astrophysical uncertainties associated to the secondary antiproton flux originating from cosmic-ray spallation on the interstellar gas. We select a set of propagation models compatible with the recent B/C data from PAMELA, and find those providing minimal and maximal antiproton fluxes in different energy ranges. We use this result to determine the most conservative bounds on relevant Dark Matter (DM) annihilation channels: We find that the recent claim of a DM interpretation of a gamma-ray excess in the Galactic Center region cannot be ruled out by current antiproton data. Finally, we discuss the impact of the recently released preliminary data from AMS-02. In particular, we provide a reference model compatible with proton, helium and B/C spectra from this experiment. Remarkably, the main propagation parameters of this model are in perfect agreement with the best fit presented in our earlier statistical analyses. We also show that the antiproton-to-proton ratio does not exhibit any significant anomaly at high energy with respect to our predictions.
We show that the event excess observed by the IceCube collaboration at TeV--PeV energies, usually interpreted as evidence for astrophysical neutrinos, can be explained alternatively by the scattering of highly boosted dark matter particles. Specifically, we consider a scenario where a $\sim 4$ PeV scalar dark matter particle $\phi$ can decay to a much lighter dark fermion $\chi$, which in turn scatters off nuclei in the IceCube detector. Besides these events, which are exclusively shower-like, the model also predicts a secondary population of events at $\mathcal{O}(100\ \text{TeV})$ originating from the 3-body decay $\phi \to \chi \bar\chi a$, where $a$ is a pseudoscalar which mediates dark matter--Standard Model interactions and whose decay products include neutrinos. This secondary population also includes track-like events, and both populations together provide an excellent fit to the IceCube data. We then argue that a relic abundance of light Dark Matter particles $\chi$, which may constitute a subdominant component of the Dark Matter in the Universe, can have exactly the right properties to explain the observed excess in GeV gamma rays from the galactic center region. Our boosted Dark Matter scenario also predicts fluxes of $\mathcal{O}(10)$ TeV positrons and $\mathcal{O}(100\ \text{TeV})$ photons from 3-body cascade decays of the heavy Dark Matter particle $\phi$, and we show how these can be used to constrain parts of the viable parameter space of the model. Direct detection limits are weak due to the pseudoscalar couplings of $\chi$. Accelerator constraints on the pseudoscalar mediator $a$ lead to the conclusion that the preferred mass of $a$ is $\gtrsim 10$ GeV and that large coupling to $b$ quarks but suppressed or vanishing coupling to leptons are preferred.
In a cyclic entropy model in which the extroverse is jettisoned at turnaround with a Come Back Empty assumption, we address matching of the contraction scale factor $f(t_T)a(t)$ to the expansion scale factor $a(t)$, where $f(t_T)$ is the ratio at turnaround of the introverse to extroverse radii. Such matching is necessary for infinite cyclicity and fixes the CBE period at $\sim 2.6$ Ty. We then compare such a CBE model with alternative cyclic cosmologies including Penrose Conformal Cyclic Cosmology of period $10^{100}$ y, and speculate that the CBE model may be related to the CCC model by a highly nontrivial isomorphism.
We present a method for determining surface flows from solar images based upon optical flow techniques. We apply the method to sets of images obtained by a variety of solar imagers to assess its performance. The {\tt opflow3d} procedure is shown to extract accurate velocity estimates when provided perfect test data and quickly generates results consistent with completely distinct methods when applied on global scales. We also validate it in detail by comparing it to an established method when applied to high-resolution datasets and find that it provides comparable results without the need to tune, filter or otherwise preprocess the images before its application.
In this paper, based on the works of Capozziello et al., we have studied the Noether symmetry approach in the cosmological model with scalar and gauge fields proposed recently by Soda et al. The correct Noether symmetries and related Lie algebra are given according to the minisuperspace quantum cosmological model. The Wheeler-De Witt (WDW) equation is presented after quantization and the classical trajectories are then obtained in the semi-classical limit. The oscillating features of the wave function in the cosmic evolution recover the so-called Hartle criterion, and the selection rule in minisuperspace quantum cosmology is strengthened. Then we have realized now the proposition that Noether symmetries select classical universes.
We have constructed a spherically symmetric structure model in a cosmological background filled with perfect fluid with non-vanishing pressure as an exact solution of Einstein equations using the Lema\^{i}tre solution. To study its local and quasi-local characteristics including the novel features of its central black hole, we have suggested an algorithm to integrate the equations numerically. The result shows intriguing effects of the pressure inside the structure. The evolution of the central black hole within the FRW universe, its decoupling from the expanding parts of the model, the structure of its space-like apparent horizon, the limiting case of the dynamical horizon tending to a slowly evolving horizon, and the decreasing mass in-fall to the black hole is also studied. We have also calculated the redshift of a light emitted from nearby the cosmological structure to an observer in the FRW background and have shown that it contains both the local gravitational and the cosmological redshift with some observational consequences. It has been shown that this type of cosmological black holes have the flexibility to match with the NFW dark matter density profile.
We study the dynamics of a plasma of charged relativistic fermions at very high temperature $T\gg m$, where $m$ is the fermion mass, coupled to the electromagnetic field. In particular, we derive a magneto-hydrodynamical description of the evolution of such a plasma. We show that, as compared to conventional MHD for a plasma of non-relativistic particles, the hydrodynamical description of the relativistic plasma involves new degrees of freedom described by a pseudo-scalar field originating in a local asymmetry in the densities of left-handed and right-handed fermions. This field can be interpreted as an effective axion field. Taking into account the chiral anomaly we present dynamical equations for the evolution of this field, as well as of other fields appearing in the MHD description of the plasma. Due to its non-linear coupling to helical magnetic fields, the axion field significantly affects the dynamics of a magnetized plasma and can give rise to a novel type of inverse cascade.
In this paper we propose and extensively study mimetic $f({\cal G})$ modified
gravity models, with various scenarios of cosmological evolution, with or
without extra matter fluids. The easiest formulation is based on the use of
Lagrange multiplier constraint.
In certain versions of this theory, it is possible to realize accelerated
expansion of the Universe or even unified evolution which includes inflation
with dark energy, and at the same time in the same theoretical framework, dark
matter is described by the theory. This is achieved by the re-parametrization
of the metric tensor, which introduces a new degree of freedom in the
cosmological equations and leads to appearance of the mimetic "dark matter"
component. In the context of mimetic $f({\cal G})$ theory, we also provide some
quite general reconstruction schemes, which enable us to find which $f({\cal
G})$ gravity generates a specific cosmological evolution. In addition, we also
provide the general reconstruction technique for Lagrange multiplier $f({\cal
G})$ gravity. All our results are accompanied by illustrative examples, with
special emphasis on bouncing cosmologies.
We study primordial tensor power-spectra generated during inflation in bimetric gravity. More precisely, we examine a homogeneous expanding spacetime in a minimal bimetric model with an inflaton and calculate tensor perturbations on the homogeneous background under slow-roll approximation. In terms of the mass eigenstates, only the power-spectrum of the massless state remains constant and both the power-spectrum of the massive state and the cross power-spectrum rapidly decay during inflation. The amplitude of the physical power-spectrum is suppressed due to the flavor mixing. All power-spectra in the flavor eigenstates coincide with each other up to the first order of the slow-roll parameter.
We present results from a numerical forward model to evaluate one-dimensional reduced power spectral densities (PSD) from arbitrary energy distributions in $\mathbf{k}$-space. In this model, we can separately calculate the diagonal elements of the spectral tensor for incompressible axisymmetric turbulence with vanishing helicity. Given a critically balanced turbulent cascade with $k_\|\sim k_\perp^\alpha$ and $\alpha<1$, we explore the implications on the reduced PSD as a function of frequency. The spectra are obtained under the assumption of Taylor's hypothesis. We further investigate the functional dependence of the spectral index $\kappa$ on the field-to-flow angle $\theta$ between plasma flow and background magnetic field from MHD to electron kinetic scales. We show that critically balanced turbulence asymptotically develops toward $\theta$-independent spectra with a slope corresponding to the perpendicular cascade. This occurs at a transition frequency $f_{2D}(L,\alpha,\theta)$, which is analytically estimated and depends on outer scale $L$, critical balance exponent $\alpha$ and field-to-flow angle $\theta$. We discuss anisotropic damping terms acting on the $\mathbf{k}$-space distribution of energy and their effects on the PSD. Further, we show that the spectral anisotropies $\kappa(\theta)$ as found by Horbury et al. (2008) and Chen et al. (2010) in the solar wind are in accordance with a damped critically balanced cascade of kinetic Alfv\'en waves. We also model power spectra obtained by von Papen & Saur (2014) in Saturn's plasma sheet and find that the change of spectral indices inside $9\,R_\mathrm{s}$ can be explained by damping on electron scales.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
Recent $N$-body simulations have shown that Einasto radial profiles provide the most accurate description of dark matter halos. Predictions based on the traditional NFW functional form may fail to describe the structural properties of cosmic objects at the percent level required by precision cosmology. We computed the systematic errors expected for weak lensing analyses of clusters of galaxies if one wrongly models the lens properties. Even though the NFW fits of observed tangential shear profiles can be excellent, viral masses and concentrations of very massive halos ($>\sim10^{15}M_\odot/h$) can be over- and underestimated by $\sim 10$ per cent, respectively. Misfitting effects also steepen the observed mass-concentration relation, in a way similar to that seen in multiwavelength observations of galaxy groups and clusters. Einasto lenses can be distinguished from NFW halos either with deep observations of very massive structures ($>\sim10^{15}M_\odot/h$) or by stacking the shear profiles of thousands of group-sized lenses ($\sim 10^{14}M_\odot/h$).
Observations of warped discs can give insight into the nature of angular momentum transport in accretion discs. Only a few objects are known to show strong periodicity on long timescales, but when such periodicity is present it is often attributed to precession of the accretion disc. The X-ray binary Hercules X-1/HZ Herculis (Her X-1) is one of the best examples of such periodicity and has been linked to disc precession since it was first observed. By using the current best-fitting models to Her X-1, which invoke precession driven by radiation warping, I place a constraint on the effective viscosities that act in a warped disc. These effective viscosities almost certainly arise due to turbulence induced by the magneto-rotational instability. The constraints derived here are in agreement with analytical and numerical investigations into the nature of magneto-hydrodynamic disc turbulence, but at odds with some recent global simulations.
While observations of large-scale structure and the cosmic microwave background (CMB) provide strong constraints on the amplitude of the primordial power spectrum (PPS) on scales larger than 10 Mpc, the amplitude of the power spectrum on sub-galactic length scales is much more poorly constrained. We study early structure formation in a cosmological model with a blue-tilted PPS. We assume that the standard scale-invariant PPS is modified at small length scales as $P(k) \sim k^{m_{\rm s}}$ with $m_{\rm s} > 1$. We run a series of cosmological hydrodynamic simulations to examine the dependence of the formation epoch and the characteristic mass of primordial stars on the tilt of the PPS. In models with $m_{\rm s} > 1$, star-forming gas clouds are formed at $z > 100$, when formation of hydrogen molecules is inefficient because the intense CMB radiation destroys chemical intermediates. Without efficient coolant, the gas clouds gravitationally contract while keeping a high temperature. The protostars formed in such "hot" clouds grow very rapidly by accretion to become extremely massive stars that may leave massive black holes with a few hundred solar-masses at $z > 100$. The shape of the PPS critically affects the properties and the formation epoch of the first generation of stars. Future experiments of the CMB polarization and the spectrum distortion may provide important information on the nature of the first stars and their formation epoch, and hence on the shape of the small-scale power spectrum.
We compare the structure of star-forming molecular clouds in different regions of Orion A to determine how the column density probability distribution function (N-PDF) varies with environmental conditions such as the fraction of young protostars. A correlation between the N-PDF slope and Class 0 protostar fraction has been previously observed in a low-mass star-formation region (Perseus) by Sadavoy; here we test if a similar correlation is observed in a high-mass star-forming region. We use Herschel data to derive a column density map of Orion A. We use the Herschel Orion Protostar Survey catalog for accurate identification and classification of the Orion A young stellar object (YSO) content, including the short-lived Class 0 protostars (with a $\sim$ 0.14 Myr lifetime). We divide Orion A into eight independent 13.5 pc$^2$ regions; in each region we fit the N-PDF distribution with a power-law, and we measure the fraction of Class 0 protostars. We use a maximum likelihood method to measure the N-PDF power-law index without binning. We find that the Class 0 fraction is higher in regions with flatter column density distributions. We test the effects of incompleteness, YSO misclassification, resolution, and pixel-scale. We show that these effects cannot account for the observed trend. Our observations demonstrate an association between the slope of the power-law N-PDF and the Class 0 fractions within Orion A. Various interpretations are discussed including timescales based on the Class 0 protostar fraction assuming a constant star-formation rate. The observed relation suggests that the N-PDF can be related to an "evolutionary state" of the gas. If universal, such a relation permits an evaluation of the evolutionary state from the N-PDF power-law index at much greater distances than those accesible with protostar counts. (abridged)
We investigate the relationship between star formation (SF) and level of relaxation in a sample of 379 galaxy clusters at z < 0.2. We use data from the Sloan Digital Sky Survey to measure cluster membership and level of relaxation, and to select star-forming galaxies based on mid-infrared emission detected with the Wide-Field Infrared Survey Explorer. For galaxies with absolute magnitudes M_r < -19.5, we find an inverse correlation between SF fraction and cluster relaxation: as a cluster becomes less relaxed, its SF fraction increases. Furthermore, in general, the subtracted SF fraction in all unrelaxed clusters (0.117 +/- 0.003) is higher than that in all relaxed clusters (0.097 +/- 0.005). We verify the validity of our SF calculation methods and membership criteria through analysis of previous work. Our results agree with previous findings that a weak correlation exists between cluster SF and dynamical state, possibly because unrelaxed clusters are less evolved relative to relaxed clusters.
Recently it has been shown that a large fraction of the dwarf satellite galaxies orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and kinematically coherent planar structure. The presence of such a structure seems to challenge the current Cold Dark Matter paradigm of structure formation, which predicts a more uniform distribution of satellites around central objects. We show that it is possible to obtain a thin, extended, rotating plane of satellites resembling the one in Andromeda in cosmological collisionless simulations based on the Cold Dark Matter model. Our new high resolution simulations show a correlation between the formation time of the dark matter halo and the thickness of the plane of satellites. Our simulations have a high incidence of satellite planes as thin, extended, and as rich as the one in Andromeda and with a very coherent kinematic structure when we select high concentration/early forming halos. By tracking the formation of the satellites in the plane we show that they have been mainly accreted onto the main object along thin dark matter filaments at high redshift. Our results show that the presence of a thin, extended, rotating plane of satellites is not a challenge for the Cold Dark Matter paradigm, but actually supports one of the predictions of this paradigm related to the presence of filaments of dark matter around galaxies at high redshift.
We study the evolution of degenerate electron cores primarily composed of the carbon burning products oxygen, neon, and magnesium (hereafter ONeMg cores) that are undergoing compression. Electron capture reactions on A=20 and A=24 isotopes reduce the electron fraction and heat the core. We develop and use a new capability of the Modules for Experiments in Stellar Astrophysics (MESA) stellar evolution code that provides a highly accurate implementation of these key reactions. These new accurate rates and the ability of MESA to perform extremely small spatial zoning demonstrates a thermal runaway in the core triggered by the temperature and density sensitivity of the Ne-20 electron capture reactions. Both analytics and numerics show that this thermal runaway does not trigger core convection, but rather leads to a centrally concentrated (r < km) thermal runaway that will subsequently launch an oxygen deflagration wave from the center of the star. We use MESA to perform a parameter study that quantifies the influence of the magnesium mass fraction, the central temperature, the compression rate, and uncertainties in the electron capture reaction rates on the ONeMg core evolution. This allows us to establish a lower limit on the central density at which the oxygen deflagration wave initiates of $\rho_c > 8.5 \times 10^9\, \textrm{g cm}^{-3}$. Based on previous work and order-of-magnitude calculations, we expect objects which ignite oxygen at or above these densities to collapse and form a neutron star. Calculations such as these are an important step in producing more realistic progenitor models for studies of the signature of accretion-induced collapse.
We present a study of 15 new brown dwarfs belonging to the $\sim7$ Myr old 25 Orionis group and Orion OB1a sub-association with spectral types between M6 and M9 and estimated masses between $\sim0.07$M$_\odot$ and $\sim0.01$ M$_\odot$. By comparing them through a Bayesian method with low mass stars ($0.8\lesssim$ M/M$_\odot\lesssim0.1$) from previous works in the 25 Orionis group, we found statistically significant differences in the number fraction of classical T Tauri stars, weak T Tauri stars, class II, evolved discs and purely photospheric emitters at both sides of the sub-stellar mass limit. Particularly we found a fraction of $3.9^{+2.4}_{-1.6}~\%$ low mass stars classified as CTTS and class II or evolved discs, against a fraction of $33.3^{+10.8}_{-9.8}~\%$ in the sub-stellar mass domain. Our results support the suggested scenario in which the dissipation of discs is less efficient for decreasing mass of the central object.
Context. Z\,CMa is a binary composed of an embedded Herbig Be and an FU Ori
class star separated by $\sim100$ au. Observational evidence indicate a complex
environment in which each star has a circumstellar disk and drives a jet, and
the whole system is embedded in a large dusty envelope.
Aims. We aim to probe the circumbinary environment of Z\,CMa in the inner 400
au in scattered light.
Methods. We use high contrast imaging polarimetry with VLT/NaCo at $H$ and
$K_s$ bands.
Results. The central binary is resolved in both bands. The polarized images
show three bright and complex structures: a common dust envelope, a sharp
extended feature previously reported in direct light, and an intriguing bright
clump located $0\farcs3$ south of the binary, which appears spatially connected
to the sharp extended feature.
Conclusions.We detect orbital motion when compared to previous observations,
and report a new outburst driven by the Herbig star. Our observations reveal
the complex inner environment of Z\,CMa with unprecedented detail and contrast.
Transition disc systems are young stars that appear to be on the verge of dispersing their protoplanetary discs. We explore the nature of these systems by comparing the stellar accretion rates and disc masses of transition discs and normal T Tauri stars in Taurus and Ophiuchus. After controlling for the known dependencies of stellar accretion rate and disc mass and on age, stellar accretion rate on stellar mass, and disc mass on the presence of stellar or sub-stellar companions, we find that the normal T Tauri stars show a trend of stellar accretion rate increasing with disc mass. The transition discs tend to have higher average disc masses than normal T Tauri stars as well as lower accretion rates than normal T Tauri stars of the same disc mass. These results are most consistent with the interpretation that the transition discs have formed objects massive enough to alter the accretion flow, i.e., single or multiple giant planets. Several Ophiuchus T Tauri stars that are not known transition disc systems also have very low accretion rates for their disc masses. We speculate on the possible nature of these sources.
With the loss of a second reaction wheel, resulting in the inability to point continuously and stably at the same field of view, the NASA Kepler satellite recently entered a new mode of observation known as the K2 mission. The data from this redesigned mission present a specific challenge; the targets systematically drift in position on a ~6 hour time scale, inducing a significant instrumental signal in the photometric time series --- this greatly impacts the ability to detect planetary signals and perform asteroseismic analysis. Here we detail our version of a reduction pipeline for K2 target pixel data, which automatically: defines masks for all targets in a given frame; extracts the target's flux- and position time series; corrects the time series based on the apparent movement on the CCD (either in 1D or 2D) combined with the correction of instrumental and/or planetary signals via the KASOC filter (Handberg & Lund 2014), thus rendering the time series ready for asteroseismic analysis; computes power spectra for all targets, and identifies potential contaminations between targets. From a test of our pipeline on a sample of targets from the K2 campaign 0, the recovery of data for multiple targets increases the amount of potential light curves by a factor ${\geq}10$. Our pipeline could be applied to the upcoming TESS (Ricker et al. 2014) and PLATO 2.0 (Rauer et al. 2013) missions.
We present a new and computationally efficient method for characterizing very low mass companions using low resolution ($R\sim$30) near-infrared ($YJH$) spectra from high contrast imaging campaigns with integral field spectrograph (IFS) units. We conduct a detailed quantitative comparison of the efficacy of this method through tests on simulated data comparable in spectral coverage and resolution to the currently operating direct imaging systems around the world. In particular, we simulate Project 1640 data as an example of the use, accuracy, and precision of this technique. We present results from comparing simulated spectra of M, L, and T dwarfs with a large and finely-sampled grid of synthetic spectra using Markov Chain Monte Carlo techniques. We determine the precision and accuracy of effective temperature and surface gravity inferred from fits to PHOENIX dusty and cond, which we find reproduce the low-resolution spectra of all objects within the adopted flux uncertainties. Uncertainties in effective temperature decrease from $\pm$100-500 K for M dwarfs to as small as $\pm$30 K for some L and T spectral types. Surface gravity is constrained to within 0.2-0.4 dex for mid-L through T dwarfs, but uncertainties are as large as 1.0 dex or more for M dwarfs. Results for effective temperature from low-resolution $YJH$ spectra generally match predictions from published spectral type-temperature relationships except for L-T transition objects and young objects. Single-band spectra (i.e., narrower wavelength coverage) result in larger uncertainties and often discrepant results, suggesting that high contrast IFS observing campaigns can compensate for low spectral resolution by expanding the wavelength coverage for reliable characterization of detected companions. [Abstract truncated]
We present a model describing the evolution of Fanaroff-Riley type I and II radio AGN, and the transition between these classes. We quantify galaxy environments using a semi-analytic galaxy formation model, and apply our model to a volume-limited low redshift ($0.03 \leqslant z \leqslant 0.1$) sample of observed AGN to determine the distribution of jet powers and active lifetimes at the present epoch. Radio sources in massive galaxies are found to remain active for longer, spend less time in the quiescent phase, and inject more energy into their hosts than their less massive counterparts. The jet power is independent of the host stellar mass within uncertainties, consistent with maintenance-mode AGN feedback paradigm. The environments of these AGN are in or close to long-term heating-cooling balance. We also examine the properties of high- and low-excitation radio galaxy sub-populations. The HERGs are younger than LERGs by an order of magnitude, whilst their jet powers are greater by a factor of four. The Eddington-scaled accretion rates and jet production efficiencies of these populations are consistent with LERGs being powered by radiatively inefficient advection dominated accretion flows (ADAFs), while HERGs are fed by a radiatively efficient accretion mechanism.
We present 13CO and C18O (1-0), (2-1), and (3-2) maps towards the core-forming Perseus B1-E clump using observations from the James Clerk Maxwell Telescope (JCMT), Submillimeter Telescope (SMT) of the Arizona Radio Observatory, and IRAM 30 m telescope. We find that the 13CO and C18O line emission both have very complex velocity structures, indicative of multiple velocity components within the ambient gas. The (1-0) transitions reveal a radial velocity gradient across B1-E of 1 km/s/pc that increases from north-west to south-east, whereas the majority of the Perseus cloud has a radial velocity gradient increasing from south-west to north-east. In contrast, we see no evidence of a velocity gradient associated with the denser Herschel-identified substructures in B1-E. Additionally, the denser substructures have much lower systemic motions than the ambient clump material, which indicates that they are likely decoupled from the large-scale gas. Nevertheless, these substructures themselves have broad line widths (0.4 km/s) similar to that of the C18O gas in the clump, which suggests they inherited their kinematic properties from the larger-scale, moderately dense gas. Finally, we find evidence of C18O depletion only toward one substructure, B1-E2, which is also the only object with narrow (transonic) line widths. We suggest that as prestellar cores form, their chemical and kinematic properties are linked in evolution, such that these objects must first dissipate their turbulence before they deplete in CO.
We use hydrodynamic simulations to study the interaction of realistic active galactic nucleus (AGN) feedback mechanisms (accretion-disk winds & Compton heating) with a multi-phase interstellar medium (ISM). Our ISM model includes radiative cooling and explicit stellar feedback from multiple processes. We simulate radii ~0.1-100 pc around an isolated (non-merging) black hole. These are the scales where the accretion rate onto the black hole is determined and where AGN-powered winds and radiation couple to the ISM. Our primary results include: (1) The black hole accretion rate on these scales is determined by exchange of angular momentum between gas and stars in gravitational instabilities. This produces accretion rates of ~0.03-1 Msun/yr, sufficient to power a luminous AGN. (2) The gas disk in the galactic nucleus undergoes an initial burst of star formation followed by several Myrs where stellar feedback suppresses the star formation rate per dynamical time. (3) AGN winds injected at small radii with momentum fluxes ~L/c couple efficiently to the ISM and have a dramatic effect on the ISM properties in the central ~100 pc. AGN winds suppress the nuclear star formation rate by a factor of ~10-30 and the black hole accretion rate by a factor of ~3-30. They increase the total outflow rate from the galactic nucleus by a factor of ~10. The latter is broadly consistent with observational evidence for galaxy-scale atomic and molecular outflows driven by AGN rather than star formation. (4) In simulations that include AGN feedback, the predicted column density distribution towards the black hole is reasonably consistent with observations, whereas absent AGN feedback, the black hole is isotropically obscured and there are not enough optically-thin sight lines to explain observed Type I AGN. A 'torus-like' geometry arises self-consistently because AGN feedback evacuates the gas in the polar regions.
Gravitational lensing is a potentially powerful tool for elucidating the origin of gamma-ray emission from distant sources. Cosmic lenses magnify the emission from distance sources and produce time delays between mirage images. Gravitationally-induced time delays depend on the position of the emitting regions in the source plane. The Fermi/LAT satellite continuously monitors the entire sky and detects gamma-ray flares, including those from gravitationally-lensed blazars. Therefore, temporal resolution at gamma-ray energies can be used to measure these time delays, which, in turn, can be used to resolve the origin of the gamma-ray flares spatially. We provide a guide to the application and Monte Carlo simulation of three techniques for analyzing these unresolved light curves: the Autocorrelation Function, the Double Power Spectrum, and the Maximum Peak Method. We apply these methods to derive time delays from the gamma-ray light curve of the gravitationally-lensed blazar PKS 1830-211. The result of temporal analysis combined with the properties of the lens from radio observations yield an improvement in spatial resolution at gamma-ray energies by a factor of 10000. We analyze four active periods. For two of these periods, the emission is consistent with origination from the core and for the other two, the data suggest that the emission region is displaced from the core by more that ~1.5 kpc. For the core emission, the gamma-ray time delays, $23\pm0.5$ days and $19.7\pm1.2$ days, are consistent with the radio time delay $26^{+4}_{-5}$ days.
The [CII] 158$\mu$m fine-structure line is known to trace regions of active star formation and is the main coolant of the cold, neutral atomic medium. In this \textit{Letter}, we report a strong detection of the [CII] line in the host galaxy of the brightest quasar known at $z>6.5$, the Pan-STARRS1 selected quasar PSO J036.5078+03.0498 (hereafter P036+03), using the IRAM NOEMA millimeter interferometer. Its [CII] and total far-infrared luminosities are $(5.8 \pm 0.7) \times 10^9 \,L_\odot$ and $(7.6\pm1.5) \times 10^{12}\,L_\odot$, respectively. This results in a $L_{[CII]} /L_{TIR}$ ratio of $\sim 0.8\times 10^{-3}$, which is at the high end for those found for active galaxies, though it is lower than the average found in typical main sequence galaxies at $z\sim 0$. We also report a tentative additional line which we identify as a blended emission from the $3_{22} - 3_{13}$ and $5_{23} - 4_{32}$ H$_2$O transitions. If confirmed, this would be the most distant detection of water emission to date. P036+03 rivals the current prototypical luminous J1148+5251 quasar at $z=6.42$, in both rest-frame UV and [CII] luminosities. Given its brightness and because it is visible from both hemispheres (unlike J1148+5251), P036+03 has the potential of becoming an important laboratory for the study of star formation and of the interstellar medium only $\sim 800\,$Myr after the Big Bang.
Time-dependent three-dimensional magnetohydrodynamics (MHD) simulation modules are implemented at the Joint Science Operation Center (JSOC) of Solar Dynamics Observatory (SDO). The modules regularly produce three-dimensional data of the time-relaxed minimum-energy state of the solar corona using global solar-surface magnetic-field maps created from Helioseismic Magnetic Imager (HMI) full-disk magnetogram data. With the assumption of polytropic gas with specific heat ratio of 1.05, three types of simulation products are currently generated: i) simulation data with medium spatial resolution using the definitive calibrated synoptic map of the magnetic field with a cadence of one Carrington rotation, ii) data with low spatial resolution using the definitive version of the synchronic frame format of the magnetic field, with a cadence of one day, and iii) low-resolution data using near-real-time (NRT) synchronic format of the magnetic field on daily basis. The MHD data available in the JSOC database are three-dimensional, covering heliocentric distances from 1.025 to 4.975 solar radii, and contain all eight MHD variables: the plasma density, temperature and three components of motion velocity, and three components of the magnetic field. This article describes details of the MHD simulations as well as the production of the input magnetic-field maps, and details of the products available at the JSOC database interface. In order to assess the merits and limits of the model, we show the simulated data in early 2011 and compare with the actual coronal features observed by the Atmospheric Imaging Assembly (AIA) and the near-Earth in-situ data.
We propose a method for observing transiting exoplanets with near-infrared high-resolution spectrometers. We aim to create a robust data analysis method for recovering atmospheric transmission spectra from transiting exoplanets over a wide wavelength range in the near infrared. By using an inverse method approach, combined with stellar models and telluric transmission spectra, the method recovers the transiting exoplanet's atmospheric transmittance at high precision over a wide wavelength range. We describe our method and have tested it by simulating observations. This method is capable of recovering transmission spectra of high enough accuracy to identify absorption features from molecules such as O2, CH4, CO2, and H2O. This accuracy is achievable for Jupiter-size exoplanetsat S/N that can be reached for 8m class telescopes using high-resolution spectrometers (R>20 000) during a single transit, and for Earth-size planets and super-Earths transiting late K or M dwarf stars at S/N reachable during observations of less than 10 transits. We also analyse potential error sources to show the robustness of the method. Detection and characterization of atmospheres of both Jupiter-size planets and smaller rocky planets looks promising using this set-up.
We report here the highest resolution near-IR imaging to date of the HD 141569A disc taken as part of the NICI Science Campaign. We recover 4 main features in the NICI images of the HD 141569 disc discovered in previous HST imaging: 1) an inner ring / spiral feature. Once deprojected, this feature does not appear circular. 2) an outer ring which is considerably brighter on the western side compared to the eastern side, but looks fairly circular in the deprojected image. 3) an additional arc-like feature between the inner and outer ring only evident on the east side. In the deprojected image, this feature appears to complete the circle of the west side inner ring and 4) an evacuated cavity from 175 AU inwards. Compared to the previous HST imaging with relatively large coronagraphic inner working angles (IWA), the NICI coronagraph allows imaging down to an IWA of 0.3". Thus, the inner edge of the inner ring/spiral feature is well resolved and we do not find any additional disc structures within 175 AU. We note some additional asymmetries in this system. Specifically, while the outer ring structure looks circular in this deprojection, the inner bright ring looks rather elliptical. This suggests that a single deprojection angle is not appropriate for this system and that there may be an offset in inclination between the two ring / spiral features. We find an offset of 4+-2 AU between the inner ring and the star center, potentially pointing to unseen inner companions.
We have developed a model to predict the post-collision brightness increase of sub-catastrophic collisions between asteroids and to evaluate the likelihood of a survey detecting these events. It is based on the cratering scaling laws of Holsapple and Housen (2007) and models the ejecta expansion following an impact as occurring in discrete shells each with their own velocity. We estimate the magnitude change between a series of target/impactor pairs, assuming it is given by the increase in reflecting surface area within a photometric aperture due to the resulting ejecta. As expected the photometric signal increases with impactor size, but we find also that the photometric signature decreases rapidly as the target asteroid diameter increases, due to gravitational fallback. We have used the model results to make an estimate of the impactor diameter for the (596) Scheila collision of D=49-65m depending on the impactor taxonomy, which is broadly consistent with previous estimates. We varied both the strength regime (highly porous and sand/cohesive soil) and the taxonomic type (S-, C- and D-type) to examine the effect on the magnitude change, finding that it is significant at early stages but has only a small effect on the overall lifetime of the photometric signal. Combining the results of this model with the collision frequency estimates of Bottke et al. (2005), we find that low-cadence surveys of approximately one visit per lunation will be insensitive to impacts on asteroids with D<20km if relying on photometric detections.
The effect of metallicity on the observed light of Type Ia supernovae (SNe Ia) could lead to systematic errors as the absolute magnitudes of local and distant SNe Ia are compared to measure luminosity distances and determine cosmological parameters. The UV light may be especially sensitive to metallicity, though different modeling methods disagree as to the magnitude, wavelength dependence, and even the sign of the effect. The outer density structure, ^56 Ni, and to a lesser degree asphericity, also impact the UV. We compute synthetic photometry of various metallicity-dependent models and compare to UV/optical photometry from the Swift Ultra-Violet/Optical Telescope. We find that the scatter in the mid-UV to near-UV colors is larger than predicted by changes in metallicity alone and is not consistent with reddening. We demonstrate that a recently employed method to determine relative abundances using UV spectra can be done using UVOT photometry, but we warn that accurate results require an accurate model of the cause of the variations. The abundance of UV photometry now available should provide constraints on models that typically rely on UV spectroscopy for constraining metallicity, density, and other parameters. Nevertheless, UV spectroscopy for a variety of SN explosions is still needed to guide the creation of accurate models. A better understanding of the influences affecting the UV is important for using SNe Ia as cosmological probes, as the UV light may test whether SNe Ia are significantly affected by evolutionary effects.
This paper presents an extensive overview of known and proposed RR Lyrae type stars in binary systems. The aim is to revise and extend the list with new Galactic field systems. We utilized maxima timings for eleven RRab type stars with suspicious behaviour available in the GEOS database, and determined maxima timings on the basis of data from various sky surveys and own observations. This significantly extended amount of suitable maxima timings. We applied non-linear least-squares method to model proposed Light Time Effect (LiTE) in O-C diagrams, and determined orbital parameters for studied systems. In contrast to recent findings, our analysis showed preference for decades-long periods instead of periods in the order of years. Secondary components were found to be low mass objects predominantly. However, for two of sample stars, RZ Cet and AT Ser, the mass of the companion of more than one solar mass suggests that it is a neutron star or a black hole. We found that semi-major axes of proposed orbits are between 1 and 20 astronomical units. Because studied stars belong to the brightest, and therefore the closest RR Lyraes, maximal angular distances between components during orbit in particular systems should be at least between 1 and 13 mas. This significantly improves the chance to detect and resolve both stars using current telescopes and interferometers. However, our interpretation of the O-C diagrams as a consequence of the LiTE should be considered as preliminary since we do not have reliable spectroscopic measurements. On the other hand our models give a prediction of the period and radial velocity evolution in years to come which should be sufficient for plausible proof of the binarity.
We present the SLoWPoKES-II catalog of low-mass visual binaries identified from the Sloan Digital Sky Survey by matching photometric distances. The candidate pairs are vetted by comparing the stellar information. The candidate pairs are vetted by comparing the stellar density at their respective Galactic positions to Monte Carlo realizations of a simulated Milky Way. In this way, we are able to identify large numbers of bona fide wide binaries without need of proper motions. 105,537 visual binaries with angular separations of $\sim$1-20", were identified, each with a probability of chance alignment of $\lesssim$5%. This is the largest catalog of bona fide wide binaries to date, and it contains a diversity of systems---in mass, mass ratios, binary separations, metallicity, and evolutionary states---that should facilitate follow-up studies to characterize the properties of M dwarfs and white dwarfs. There is a subtle but definitive suggestion of multiple populations in the physical separation distribution, supporting earlier findings. We suggest that wide binaries are comprised of multiple populations, most likely representing different formation modes. There are 141 M7 or later wide binary candidates, representing a 7-fold increase in the number currently known. These binaries are too wide to have been formed via the ejection mechanism. Finally, we found a 6% of spectroscopically confirmed M dwarfs are not included in the SDSS STAR catalog; they are misclassified as extended sources due to the presence of nearby or partially resolved companion. The SLoWPoKES-II catalog is publicly available to the entire community on the world wide web via the Filtergraph data visualization portal.
We present the rest-frame optical spectral energy distribution and stellar masses of six Herschel- selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies were first identified with Herschel/SPIRE imaging data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The targets were observed with Spitzer/IRAC at 3.6 and 4.5um. Due to the spatial resolution of the IRAC observations at the level of 2 arcseconds, the lensing features of a background DSFG in the near-infrared are blended with the flux from the foreground lensing galaxy in the IRAC imaging data. We make use of higher resolution Hubble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles of the foreground lensing galaxy (or galaxies) as a way to model the foreground components, in order to successfully disentangle the foreground lens and background source flux densities in the IRAC images. The flux density measurements at 3.6 and 4.5um, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical spectral energy distribution of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs to establish the stellar mass, dust mass, star-formation rate, visual extinction, and other parameters for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 8 x 10^10 to 4 x 10^11 Msun and star-formation rates of around 100 Msun yr-1. This puts these lensed sub-millimeter systems well above the SFR-M* relation observed for normal star-forming galaxies at similar redshifts. The high values of SFR inferred for these systems are consistent with a major merger-driven scenario for star formation.
The Herschel Multi-tiered Extragalactic Survey (HerMES) has identified large numbers of dusty star-forming galaxies (DSFGs) over a wide range in redshift. A detailed understanding of these DSFGs is hampered by the poor spatial resolution of Herschel. We present 870um 0.45" imaging obtained in Cycle 0 with the Atacama Large Millimeter/submillimeter Array (ALMA) of a sample of 29 HerMES DSFGs. The ALMA imaging reveals that these DSFGs comprise a total of 62 sources (down to the 5-sigma limit in our ALMA sample; sigma~0.2 mJy). Optical imaging indicates that 36 of the ALMA sources experience a significant flux boost from gravitational lensing (mu>1.1), but only 6 are strongly lensed and show multiple images. We introduce and make use of uvmcmcfit, a general purpose and publicly available Markov chain Monte Carlo visibility plane analysis tool to analyze the source properties. Combined with our previous work on brighter Herschel sources, the lens models presented here tentatively favor intrinsic number counts for DSFGs with a break near 8 mJy at 880um and a steep fall off at higher flux densities. Nearly 70% of the Herschel sources break down into multiple ALMA counterparts, consistent with previous research indicating that the multiplicity rate is high in bright sources discovered in single-dish sub-mm or FIR surveys. The ALMA counterparts to our Herschel targets are located significantly closer to each other than ALMA counterparts to sources found in the LABOCA ECDFS Submillimeter Survey. Theoretical models underpredict the excess number of sources with small separations seen in our ALMA sample. The high multiplicity rate and low projected separations between sources seen in our sample argue in favor of interactions and mergers plausibly driving both the prodigious emission from the brightest DSFGs as well as the sharp downturn above S_880 = 8 mJy.
The Chinese Small Telescope ARray (CSTAR) has observed an area around the Celestial South Pole at Dome A since 2008. About $20,000$ light curves in the i band were obtained lasting from March to July, 2008. The photometric precision achieves about 4 mmag at i = 7.5 and 20 mmag at i = 12 within a 30 s exposure time. These light curves are analyzed using Lomb--Scargle, Phase Dispersion Minimization, and Box Least Squares methods to search for periodic signals. False positives may appear as a variable signature caused by contaminating stars and the observation mode of CSTAR. Therefore the period and position of each variable candidate are checked to eliminate false positives. Eclipsing binaries are removed by visual inspection, frequency spectrum analysis and locally linear embedding technique. We identify 53 eclipsing binaries in the field of view of CSTAR, containing 24 detached binaries, 8 semi-detached binaries, 18 contact binaries, and 3 ellipsoidal variables. To derive the parameters of these binaries, we use the Eclipsing Binaries via Artificial Intelligence (EBAI) method. The primary and the secondary eclipse timing variations (ETVs) for semi-detached and contact systems are analyzed. Correlated primary and secondary ETVs confirmed by false alarm tests may indicate an unseen perturbing companion. Through ETV analysis, we identify two triple systems (CSTAR J084612.64-883342.9 and CSTAR J220502.55-895206.7). The orbital parameters of the third body in CSTAR J220502.55-895206.7 are derived using a simple dynamical model.
Transient black hole candidates are interesting objects to study in X-rays as these sources show rapid evolutions in their spectral and temporal properties. In this paper, we study the spectral properties of the Galactic transient X-ray binary MAXI~J1659-152 during its very first outburst after discovery with the archival data of RXTE Proportional Counter Array instruments. We make a detailed study of the evolution of accretion flow dynamics during its 2010 outburst through spectral analysis using the Chakrabarti-Titarchuk two-component advective flow (TCAF) model as an additive table model in XSPEC. Accretion flow parameters (Keplerian disk and sub-Keplerian halo rates, shock location and shock strength) are extracted from our spectral fits with TCAF. We studied variations of these fit parameters during the entire outburst as it passed through three spectral classes: hard, hard-intermediate, and soft-intermediate. We compared our TCAF fitted results with standard combined disk black body (DBB) and power-law (PL) model fitted results and found that variations of disk rate with DBB flux and halo rate with PL flux are generally similar in nature. There appears to be an absence of the soft state unlike what is seen in other similar sources.
Recent observations by the Large Area Telescope (LAT) onboard the Fermi satellite have revealed bright gamma-ray emission from middle-aged supernova remnants (SNRs) inside our Galaxy. These remnants, which also possess bright non-thermal radio shells, are often found to be interacting directly with surrounding gas clouds. We explore the non-thermal emission mechanism at these dynamically evolved SNRs by constructing a hydrodynamical model. Two scenarios of particle acceleration, either a re-acceleration of Galactic cosmic rays (CRs) or an efficient nonlinear diffusive shock acceleration (NLDSA) of particles injected from downstream, are considered. Using parameters inferred from observations, our models are contrasted with the observed spectra of SNR W44. For the re-acceleration case, we predict a significant enhancement of radio and GeV emission as the SNR undergoes a transition into the radiative phase. If sufficiently strong magnetic turbulence is present in the molecular cloud, the re-acceleration scenario can explain the observed broadband spectral properties. The NLDSA scenario also succeeds in explaining the $\gamma$-ray spectrum but fails to reproduce the radio spectral index. Efficient NLDSA also results in a significant post-shock non-thermal pressure that limits the compression during cooling and prevents the formation of a prominent dense shell. Some other interesting differences between the two models in hydrodynamical behavior and resulting spectral features are illustrated.
Solar flares typically have an impulsive phase that followed by a gradual phase as best seen in soft X-ray emissions. A recent discovery based on the EUV Variability Experiment (EVE) observations onboard the Solar Dynamics Observatory (SDO) reveals that some flares exhibit a second large peak separated from the first main phase peak by tens of minutes to hours, which is coined as the flare's EUV late phase. In this paper, we address the origin of the EUV late phase by analyzing in detail two late phase flares, an M2.9 flare on 2010 October 16 and an M1.4 flare on 2011 February 18, using multi-passband imaging observations from the Atmospheric Imaing Assembly (AIA) onboard SDO. We find that: (1) the late phase emission originates from a different magnetic loop system, which is much larger and higher than the main phase loop system. (2) The two loop systems have different thermal evolution. While the late phase loop arcade reaches its peak brightness progressively at a later time spanning for more than one hour from high to low temperatures, the main phase loop arcade reaches its peak brightness at almost the same time (within several minutes) in all temperatures. (3) Nevertheless, the two loop systems seem to be connected magnetically, forming an asymmetric magnetic quadruple configuration. (4) Further, the footpoint brightenings in UV wavelengths show a systematic delay of about one minute from the main flare region to the remote footpoint of the late phase arcade system. We argue that the EUV late phase is the result of a long-lasting cooling process in the larger magnetic arcade system.
The second peak in the Fe XVI 33.5 nm line irradiance observed during solar flares by Extreme ultraviolet Variability Experiment (EVE) is known as Extreme UltraViolet (EUV) late phase. Our previous paper (Liu et al. 2013) found that the main emissions in the late phase are originated from large-scale loop arcades that are closely connected to but different from the post flare loops (PFLs), and we also proposed that a long cooling process without additional heating could explain the late phase. In this paper, we define the extremely large late phase because it not only has a bigger peak in the warm 33.5 irradiance profile, but also releases more EUV radiative energy than the main phase. Through detailedly inspecting the EUV images from three point-of-view, it is found that, besides the later phase loop arcades, the more contribution of the extremely large late phase is from a hot structure that fails to erupt. This hot structure is identified as a flux rope, which is quickly energized by the flare reconnection and later on continuously produces the thermal energy during the gradual phase. Together with the late-phase loop arcades, the fail to erupt flux rope with the additional heating create the extremely large EUV late phase.
Measurements of solar wind turbulence reveal the ubiquity of discontinuities. In this study, we investigate how the discontinuities, especially rotational discontinuities (RDs), are formed in magnetohydrodynamic (MHD) turbulence. In a simulation of the decaying compressive three-dimensional (3-D) MHD turbulence with an imposed uniform background magnetic field, we detect RDs with sharp field rotations and little variations of magnetic field intensity as well as mass density. At the same time, in the de Hoffman-Teller (HT) frame, the plasma velocity is nearly in agreement with the Alfv\'{e}n speed, and is field-aligned on both sides of the discontinuity. We take one of the identified RDs to analyze in details its 3-D structure and temporal evolution. By checking the magnetic field and plasma parameters, we find that the identified RD evolves from the steepening of the Alfv\'{e}n wave with moderate amplitude, and that steepening is caused by the nonuniformity of the Alfv\'{e}n speed in the ambient turbulence.
Grand rotation curves (GRC) within ~400 kpc of M31 and the Milky Way were constructed by combining disk rotation velocities and radial velocities of satellite galaxies and globular clusters. The GRC for the Milky Way was revised using the most recent Solar rotation velocity. The derived GRCs were deconvolved into a de Vaucouleurs bulge, exponential disk, and a dark halo with the Navarro-Frenk-White (NFW) density profile by the least chi-squares fitting. Comparison of the best-fit parameters revealed similarity of the disks and bulges of the two galaxies, whereas the dark halo mass of M31 was found to be twice the Galaxy's. We show that the NFW model may be a realistic approximation to the observed dark halos in these two giant spirals.
In this letter we present results from intra-night monitoring in three colors (BRI) of the blazar S4 0954+65 during its recent (2015 February) unprecedented high state. We find violent variations on very short time scales, reaching magnitude change levels of 0.1-0.2 mag/h. On some occasions, changes of ~0.1 mag are observed even within ~10 min. During the night of 14.02.2015 an exponential drop of ~0.7 magnitudes is detected for about 5 hours. Cross-correlation between the light curves does not reveal any detectable wavelength-dependent time delays, larger than ~5 min. Color changes "bluer-when-brighter" are observed on longer time scales. Possible variability mechanisms to explain the observations are discussed and a preference to the geometrical one is given.
Context. The abundances of the three main isotopes of oxygen are altered in the course of the CNO-cycle. When the first dredge-up mixes the burning products to the surface, the nucleosynthesis processes can be probed by measuring oxygen isotopic ratios. Aims. By measuring 16O/17O and 16O/18O in red giants of known mass we compare the isotope ratios with predictions from stellar and galactic evolution modelling. Methods. Oxygen isotopic ratios were derived from the K-band spectra of six red giants. The sample red giants are open cluster members with known masses of between 1.8 and 4.5 Msun . The abundance determination employs synthetic spectra calculated with the COMARCS code. The effect of uncertainties in the nuclear reaction rates, the mixing length, and of a change in the initial abundance of the oxygen isotopes was determined by a set of nucleosynthesis and mixing models using the FUNS code. Results. The observed 16O/17O ratios are in good agreement with the model results, even if the measured values do not present clear evidence of a variation with the stellar mass. The observed 16O/18O ratios are clearly lower than the predictions from our reference model. Variations in nuclear reaction rates and mixing length parameter both have only a very weak effect on the predicted values. The 12C/13C ratios of the K giants studied implies the absence of extra-mixing in these objects. Conclusions. A comparison with galactic chemical evolution models indicates that the 16O/18O abundance ratio underwent a faster decrease than predicted. To explain the observed ratios, the most likely scenario is a higher initial 18O abundance combined with a lower initial 16 O abundance. Comparing the measured 18 O/17 O ratio with the corresponding value for the ISM points towards an initial enhancement of 17O as well. Limitations imposed by the observations prevent this from being a conclusive result.
Context. Recent observations of solar prominences show the presence of
turbulent flows that may be caused by Kelvin-Helmholtz instabilites (KHI).
However, the observed flow velocities are below the classical threshold for the
onset of KHI in fully ionized plasmas.
Aims. We investigate the effect of partial ionization on the onset of KHI in
dense and cool cylindrical magnetic flux tubes surrounded by a hotter and
lighter environment.
Methods. The linearized governing equations of a partially ionized two-fluid
plasma are used to describe the behavior of small-amplitude perturbations
superimposed on a magnetic tube with longitudinal mass flow. A normal mode
analysis is performed to obtain the dispersion relation for linear
incompressible waves. We focus on the appearance of unstable solutions and
study the dependence of their growth rates on various physical parameters. An
analytical approximation of the KHI linear growth rate for slow flows and
strong ion-neutral coupling is obtained. An application to solar prominence
threads is given.
Results. The presence of a neutral component in a plasma may contribute to
the onset of the KHI even for sub-Alfv\'enic longitudinal shear flows.
Collisions between ions and neutrals reduce the growth rates of the unstable
perturbations but cannot completely suppress the instability.
Conclusions. Turbulent flows in solar prominences with sub-Alfv\'enic flow
velocities may be interpreted as consequences of KHI in partially ionized
plasmas.
We have found two molecular clouds having velocities of 2\,km\,s$^{-1}$ and 14\,km\,s$^{-1}$ toward the super star cluster RCW\,38 by observations of CO ($J=$1--0 and 3--2) transitions. The two clouds are likely physically associated with the cluster as verified by the high intensity ratio of the $J$=3--2 emission to the $J$=1--0 emission, the bridging feature connecting the two clouds in velocity and the morphological correspondence with the infrared dust emission. Since the total mass of the clouds and the cluster is too small to gravitationally bind the velocity, we suggest that a collision happened by chance between the two clouds. We present a scenario that the collision triggered formation of the $\sim$20 candidate O stars which are localized, within $\sim$0.5\,pc of the cluster center in the 2\,km\,s$-1$ cloud, just toward the northern tip of the 14\,km\,s$^{-1}$ cloud. The other member low-mass stars are likely pre-existent prior to the collision since they are distributed outside the 14\,km\,s$^{-1}$ cloud. The formation timescale of the O stars is estimated to be $\sim3\times10^4$\,yrs ($=\sim$0.5\,pc/16\,km\,s$^{-1}$), implying a mass accretion rate $\sim10^{-3}$\,$M_\odot$\,yr$^{-1}$ for a 20\,$M_\odot$ star. This is the third super star cluster alongside of Westerlund 2 and NGC\,3603 where cloud-cloud collision triggered the cluster formation among the few youngest super star clusters that are associated with dusty nebulae.
Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.
We study the relation between accretion, black hole mass and jet power in AGN, by using a large group of blazars detected by the Fermi Large Area Telescope and radio galaxies. Our main results are as follows. (i) The jet power of FSRQs and FRII-HEG depends on the black hole mass, which suggests that the FSRQs and FRII-HEG are in Radiation-Pressure Dominated regime. The jet power of BL Lacs and FRI-LEG depends on the accretion, which suggests that the BL Lacs and FRI-LEG are in the Gas-Pressure Dominated regime. (ii) We find that most of FSRQs and BL Lacs have $\rm{P_{jet}>L_{BZ}^{max}}$, which suggests that the Blandford-Znajek mechanism is insufficient to explain the jet power of these objects. (iii) The FSRQs are roughly separated from BL Lacs by the Ledlow-Owen's dividing line in the $\rm{\log P_{jet}-\log M}$ plane, which supports the unified scheme of AGN. (iv) The FSRQs and BL Lacs have a clear division at $\rm{L_{bol}/L_{Edd}\sim0.01}$, and the distribution of Eddington ratios of BL Lacs and FSRQs exhibits a bimodal nature, which imply that the accretion mode of FSRQs may be different from that of BL Lacs. (v) We find a significant correlation between broad line luminosity and jet power, which supports a direct tight connection between jet power and accretion.
The formation of the Earth's core is a consequence of planetary accretion and processes in the Earth's interior. The mechanical process of planetary differentiation is likely to occur in large, if not global, magma oceans created by the collisions of planetary embryos. Metal-silicate segregation in magma oceans occurs rapidly and efficiently unlike grain scale percolation according to laboratory experiments and calculations. Geochemical models of the core formation process as planetary accretion proceeds are becoming increasingly realistic. Single stage and continuous core formation models have evolved into multi-stage models that are couple to the output of dynamical models of the giant impact phase of planet formation. The models that are most successful in matching the chemical composition of the Earth's mantle, based on experimentally-derived element partition coefficients, show that the temperature and pressure of metal-silicate equilibration must increase as a function of time and mass accreted and so must the oxygen fugacity of the equilibrating material. The latter can occur if silicon partitions into the core and through the late delivery of oxidized material. Coupled dynamical accretion and multi-stage core formation models predict the evolving mantle and core compositions of all the terrestrial planets simultaneously and also place strong constraints on the bulk compositions and oxidation states of primitive bodies in the protoplanetary disk.
We employed published rotation periods of {\it Kepler} field stars to test whether stars hosting planets tend to rotate more slowly than stars without known planets. Spectroscopic vsini observations of nearby stars with planets have indicated that they tend to have smaller visni values. We employ data for {\it Kepler} Objects of Interest (KOIs) from the first 16 quarters of its original mission; stellar parameters are based on the analysis of the first 17 quarters. We confirm that KOI stars rotate more slowly with much greater confidence than we had previously found for nearby stars with planets. Furthermore, we find that stars with planets of all types rotate more slowly, not just stars with giant planets.
We present a comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.
The alignments between galaxies, their underlying matter structures, and the cosmic web constitute vital ingredients for a comprehensive understanding of gravity, the nature of matter, and structure formation in the Universe. We provide an overview on the state of the art in the study of these alignment processes and their observational signatures, aimed at a non-specialist audience. The development of the field over the past one hundred years is briefly reviewed. We also discuss the impact of galaxy alignments on measurements of weak gravitational lensing, and discuss avenues for making theoretical and observational progress over the coming decade.
A new family classification, based on a catalog of proper elements with $\sim
384,000$ numbered asteroids and on new methods is available. For the $45$
dynamical families with $>250$ members identified in this classification, we
present an attempt to obtain statistically significant ages: we succeeded in
computing ages for $37$ collisional families. We used a rigorous method,
including a least squares fit of the two sides of a V-shape plot in the proper
semimajor axis, inverse diameter plane to determine the corresponding slopes,
an advanced error model for the uncertainties of asteroid diameters, an
iterative outlier rejection scheme and quality control. The best available
Yarkovsky measurement was used to estimate a calibration of the Yarkovsky
effect for each family. The results are presented separately for the families
originated in fragmentation or cratering events, for the young, compact
families and for the truncated, one-sided families. For all the computed ages
the corresponding uncertainties are provided. We found 2 cases where two
separate dynamical families form together a single V-shape with compatible
slopes, thus indicating a single collisional event. We have also found 3
examples of dynamical families containing multiple collisional families, plus a
dubious case. We have found 2 cases of families containing a conspicuous
subfamily, such that it is possible to measure the slope of a distinct V-shape,
thus the age of the secondary collision. We also provide data on the central
gaps appearing in some families.
The ages computed in this paper are obtained with a single and uniform
methodology, thus the ages of different families can be compared, providing a
first example of collisional chronology of the asteroid main belt.
Galaxy shapes are not randomly oriented, rather they are statistically aligned in a way that can depend on formation environment, history and galaxy type. Studying the alignment of galaxies can therefore deliver important information about the astrophysics of galaxy formation and evolution as well as the growth of structure in the Universe. In this review paper we summarise key measurements of intrinsic alignments, divided by galaxy type, scale and environment. We also cover the statistics and formalism necessary to understand the observations in the literature. With the emergence of weak gravitational lensing as a precision probe of cosmology, galaxy alignments took on an added importance because they can mimic cosmic shear, the effect of gravitational lensing by large-scale structure on observed galaxy shapes. This makes intrinsic alignments an important systematic effect in weak lensing studies. We quantify the impact of intrinsic alignments on cosmic shear surveys and finish by reviewing practical mitigation techniques which attempt to remove contamination by intrinsic alignments.
Core-collapse supernovae produce an intense burst of electron antineutrinos in the few-tens-of-MeV range. Several Large Liquid Scintillator-based Detectors (LLSD) are currently operated worldwide, being very effective for low energy antineutrino detection through the Inverse Beta Decay (IBD) process. In this article, we develop a procedure for the prompt extraction of the supernova location by revisiting the details of IBD kinematics over the broad energy range of supernova neutrinos. Combining all current scintillator-based detector, we show that one can locate a canonical supernova at 10 kpc with an accuracy of 45 degrees (68% C.L.). After the addition of the next generation of scintillator-based detectors, the accuracy could reach 12 degrees (68% C.L.), therefore reaching the performances of the large water Cerenkov neutrino detectors. We also discuss a possible improvement of the SuperNova Early Warning System (SNEWS) inter-experiment network with the implementation of a directionality information in each experiment. Finally, we discuss the possibility to constrain the neutrino energy spectrum as well as the mass of the newly born neutron star with the LLSD data
We present and discuss highly accurate photometry obtained through medium Stromgren y,b bands and narrow [OIII], Halpha bands covering 500 days of the evolution of Nova Del 2013 since its maximum brightness. This is by far the most complete study of any nova observed in such photometric systems. The nova behaviour in these photometric bands is very different from that observed with the more conventional broad bands like UBVRI or ugriz, providing unique information about extent and ionization of the ejecta, the onset of critical phases like the transition between optically thick and thin conditions, and re-ionization by the central super-soft X-ray source. The actual transmission profiles of the y, b, [OIII] and Halpha photometric filters have been accurately measured at different epochs and different illumination angles, to evaluate in detail their performance under exact operating conditions. The extreme smoothness of both the Halpha and [OIII] lightcurves argues for absence of large and abrupt discontinuities in the ejecta of Nova Del 2013. Should they exist, glitches in the lightcurves would have appeared when the ionization and/or recombination fronts overtook them. During the period of recorded very large variability (up to 100x over a single day) in the super-soft X-ray emission (from day +69 to +86 past V maximum), no glitch in excess of 1% was observed in the optical photometry, either in the continuum (Stromgren y) or in the lines ([OIII] and Halpha filters), or in a combination of the two (Stromgren b, Johnson B and V). Considering that the recombination time scale in the ejecta was one week at that time, this excludes global changes of the white dwarf as the source of the X-ray variability and supports instead clumpy ejecta passing through the line of sight to us as its origin.
We investigate the effect of equilateral-type primordial trispectrum on the halo/galaxy bispectrum. We consider three types of equilateral primordial trispectra which are generated by quartic operators naturally appeared in the effective field theory of inflation and can be characterized by three non-linearity parameters, $g_{\rm NL} ^{\dot{\sigma}^4}$, $g_{\rm NL} ^{\dot{\sigma}^2 (\partial \sigma)^2}$, and $g_{\rm NL} ^{(\partial \sigma)^4}$. Recently, constraints on these parameters have been investigated from CMB observations by using WMAP9 data. In order to consider the halo/galaxy bispectrum with the equilateral-type primordial trispectra, we adopt the integrated Perturbation Theory (iPT) in which the effects of primordial non-Gaussianity are wholly encapsulated in the linear primordial polyspectrum for the evaluation of the biased polyspectrum. We show the shapes of the halo/galaxy bispectrum with the equilateral-type primordial trispectra, and find that the primordial trispectrum characterized by $g_{\rm NL} ^{\dot{\sigma}^4}$ provides the same scale-dependence as the gravity-induced halo/galaxy bispectrum. Hence, it would be difficult to obtain the constraint on $g_{\rm NL} ^{\dot{\sigma}^4}$ from the observations of the halo/galaxy bispectrum. On the other hand, the primordial trispectra characterized by $g_{\rm NL} ^{\dot{\sigma}^2 (\partial \sigma)^2}$ and $g_{\rm NL} ^{(\partial \sigma)^4}$ provide the common scale-dependence which is different from that of the gravity-induced halo/galaxy bispectrum on large scales. Hence future observations of halo/galaxy bispectrum would give constraints on the non-linearity parameters, $g_{\rm NL} ^{\dot{\sigma}^2 (\partial \sigma)^2}$ and $g_{\rm NL} ^{(\partial \sigma)^4}$ independently from CMB observations and it is expected that these constraints can be comparable to ones obtained by CMB.
We have observed high-dispersion echelle spectra of red giant members in the five open clusters NGC 1342, NGC 1662, NGC 1912, NGC 2354 and NGC 2447 and determined their radial velocities and chemical compositions. These are the first chemical abundance measurements for all but NGC 2447. We combined our clusters from this and previous papers with a sample drawn from the literature for which we remeasured the chemical abundances to establish a common abundance scale. With this homogeneous sample of open clusters, we study the relative elemental abundances of stars in open clusters in comparison with field stars as a function of age and metallicity. We find a range of mild enrichment of heavy (Ba-Eu) elements in young open cluster giants over field stars of the same metallicity. Our analysis succinct that the youngest stellar generations in cluster might be under-represented by the solar neighbourhood field stars.
Cappellari (2008) presented a flexible and efficient method to model the stellar kinematics of anisotropic axisymmetric and spherical stellar systems. The spherical formalism could be used to model the line-of-sight velocity second moments allowing for essentially arbitrary radial variation in the anisotropy and general luminous and total density profiles. Here we generalize the spherical formalism by providing the expressions for all three components of the projected second moments, including the two proper motion components. A reference implementation is now included in the public JAM package available at this http URL
In this paper, we analyze the linear stability of a stellar accretion disk, having a stratified morphology. The study is performed in the framework of ideal magneto-hydrodynamics and therefore it results in a characterization of the linear unstable magneto-rotational modes. The peculiarity of the present scenario consists of adopting the magnetic flux function as the basic dynamical variable. Such a representation of the dynamics allows to make account of the co-rotation theorem as a fundamental feature of the ideal plasma equilibrium, evaluating its impact on the perturbation evolution too. According to the Alfvenic nature of the Magneto-rotational instability, we consider an incompressible plasma profile and perturbations propagating along the background magnetic field. Furthermore, we develop a local perturbation analysis, around fiducial coordinates of the background configuration and dealing with very small scale of the linear dynamics in comparison to the background inhomogeneity size. The main issue of the present study is that the condition for the emergence of unstable modes is the same in the stratified plasma disk, as in the case of a thin configuration. Such a feature is the result of the cancelation of the vertical derivative of the disk angular frequency from the dispersion relation, which implies that only the radial profile of the differential rotation is responsible for the trigger of growing modes.
The shapes of galaxies are not randomly oriented on the sky. During the galaxy formation and evolution process, environment has a strong influence, as tidal gravitational fields in large-scale structure tend to align the shapes and angular momenta of nearby galaxies. Additionally, events such as galaxy mergers affect the relative alignments of galaxies throughout their history. These "intrinsic galaxy alignments" are known to exist, but are still poorly understood. This review will offer a pedagogical introduction to the current theories that describe intrinsic galaxy alignments, including the apparent difference in intrinsic alignment between early- and late-type galaxies and the latest efforts to model them analytically. It will then describe the ongoing efforts to simulate intrinsic alignments using both $N$-body and hydrodynamic simulations. Due to the relative youth of this field, there is still much to be done to understand intrinsic galaxy alignments and this review summarises the current state of the field, providing a solid basis for future work.
Observations show that the lower thermosphere of Mars ($\sim$100--140 km) is up to 40 K colder than the current general circulation models (GCMs) can reproduce. Possible candidates for physical processes missing in the models are larger abundances of atomic oxygen facilitating stronger CO$_2$ radiative cooling, and thermal effects of gravity waves. Using two state-of-the-art Martian GCMs, the Laboratoire de M\'et\'eorologie Dynamique and Max Planck Institute models that self-consistently cover the atmosphere from the surface to the thermosphere, these physical mechanisms are investigated. Simulations demonstrate that the CO$_2$ radiative cooling with a sufficiently large atomic oxygen abundance, and the gravity wave-induced cooling can alone result in up to 40 K colder temperature in the lower thermosphere. Accounting for both mechanisms produce stronger cooling at high latitudes. However, radiative cooling effects peak above the mesopause, while gravity wave cooling rates continuously increase with height. Although both mechanisms act simultaneously, these peculiarities could help to further quantify their relative contributions from future observations.
The evolution of a circumstellar disk from its gas-rich protoplanetary to gas-poor debris stage is not well understood. It is apparent that disk-clearing progresses from the inside-out on a short time-scale, and photoevaporation models are frequently invoked to explain this process. However, the photoevaporation rates predicted by recent models differ by up to two orders of magnitude, resulting in uncertain time-scales for the final stages of disk clearing. The best candidates for studying this stage are weak line T Tauri stars (WTTS) with significant IR excess. We here aim to provide observational constraints on theories of disk-clearing by measuring the dust masses and CO content of a sample of such WTTS. We use ALMA band-6 to obtain continuum and $^{12}$CO(2-1) line fluxes for a sample of 24 WTTS stars with a known IR-excess. For these systems, we infer the dust mass from the continuum observations, and derive disk luminosities and ages to allow comparison with previously detected systems. We detect continuum emission in only 4 of 24 systems, and no $^{12}$CO(2-1) emission in any. For those systems without a continuum detection, the dust mass and fractional disk luminosity upper-limits suggest they are in the debris disk regime, making them some of the youngest debris disks known. Of those with a continuum detection, three are possible photoevaporating disks but photodissociation has likely reduced the CO abundance to below our detection limit. The low fraction of continuum detections implies that once accretion onto the star stops, the clearing of the majority of dust progresses very rapidly. Most WTTS with IR excess are not in transition but resemble debris disks. The dust in these disks is either primordial and survived the disk clearing, or is of second generation origin. In the latter case, the presence of giant planets within these systems might be the cause.
The existence of the cosmic neutrino background (CnuB) is a fundamental prediction of the standard Big Bang cosmology. Although current cosmological probes provide indirect observational evidence, the direct detection of the CnuB in a laboratory experiment is a great challenge to the present experimental techniques. We discuss the future prospects for the direct detection of the CnuB, with the emphasis on the method of captures on beta-decaying nuclei and the PTOLEMY project. Other possibilities using the electron-capture (EC) decaying nuclei, the annihilation of extremely high-energy cosmic neutrinos (EHEC\nus) at the Z-resonance, and the atomic de-excitation method are also discussed in this review.
While dark matter self-interactions may solve several problems with structure formation, so far only the effects of two-body scatterings of dark matter particles have been considered. We show that, if a subdominant component of dark matter is charged under an unbroken $U(1)$ gauge group, collective dark plasma effects need to be taken into account to understand its dynamics. Plasma instabilities can lead to collisionless dark matter shocks in galaxy cluster mergers which might have been already observed in the Abell 3827 and 520 clusters. As a concrete model we propose a thermally produced dark pair plasma of vectorlike fermions. In this scenario the interacting dark matter component is expected to be separated from the stars and the non-interacting dark matter halos in cluster collisions. In addition, the missing satellite problem is softened, while constraints from all other astrophysical and cosmological observations are avoided.
We investigate the relation between the annihilation of dark matter (DM) particles into lepton pairs and into 2-body final states including one or two photons. We parametrize the DM interactions with leptons in terms of contact interactions, and calculate the loop-level annihilation into monochromatic gamma rays, specifically computing the ratio of the DM annihilation cross sections into two gamma rays versus lepton pairs. While the loop-level processes are generically suppressed in comparison with the tree-level annihilation into leptons, we find that some choices for the mediator spin and coupling structure lead to large branching fractions into gamma-ray lines. This result has implications for a dark matter contribution to the AMS-02 positron excess. We also explore the possibility of mediators which are charged under a dark symmetry and find that, for these loop-level processes, an effective field theory description is accurate for DM masses up to about half the mediator mass.
We consider the sector of Horndeski's gravity characterized by a coupling between the kinetic scalar field term and the Einstein tensor. Our goal is to find realistic neutron star configurations in this framework. We show that, in a certain limit, there exist solutions that are identical to the Schwarzschild metric outside the star but change considerably inside, where the scalar field is not trivial. We study numerically the equations and find the region of the parameter space where neutron stars exist. We determine their internal pressure and mass-radius relation, and we compare them with standard general relativity models.
We present a unifying treatment of dark energy and modified gravity that allows distinct conformal-disformal couplings of matter species to the gravitational sector. In this very general approach, we derive the conditions to avoid ghost and gradient instabilities. We compute the equations of motion for background quantities and linear perturbations. We illustrate our formalism with two simple scenarios, where either cold dark matter or a relativistic fluid is nonminimally coupled. This extends previous studies of coupled dark energy to a much broader spectrum of gravitational theories.
The AMS-02 collaboration has recently reported the antiproton to proton ratio with improved accuracy. In view of uncertainties of the production and the propagation of the cosmic rays, the observed ratio is still consistent with the secondary astrophysical antiproton to proton ratio. However, it is nonetheless enticing to examine whether the observed spectrum can be explained by a strongly motivated dark matter, the wino dark matter. As we will show, we find that the antiproton flux from the wino annihilation can explain the observed spectrum well for its mass range 2.5-3 TeV. The fit to data becomes particularly well compared to the case without the annihilation for the thermal wino dark matter case with a mass about 3 TeV. The ratio is predicted to be quickly decreased at the energy several hundreds of GeV, if this possibility is true, and it will be confirmed or ruled out in near future when the AMS-02 experiment accumulates enough data at this higher energy region.
Cosmological alpha-attractors give a natural explanation for the spectral index n_s of inflation as measured by Planck while predicting a range for the tensor-to-scalar ratio r, consistent with all observations, to be measured more precisely in future detection of gravity waves. Their embedding into supergravity exploits the hyperbolic geometry of the Poincare disk or half-plane. These geometries are isometric under Mobius transformations, which include the shift symmetry of the inflaton field. We introduce a new Kahler potential frame that explicitly preserves this symmetry, enabling the inflaton to be light. Moreover, we include higher-order curvature deformations, which can stabilize a direction orthogonal to the inflationary trajectory. We illustrate this new framework by stabilizing the single superfield alpha-attractors.
The measurement of frequency shifts for light beams exchanged between two test masses nearly in free fall is at the heart of gravitational wave detection. It is envisaged that the derivative of the frequency shift is in fact limited by differential forces acting on those test masses. We calculate the derivative of the frequency shift with a fully covariant, gauge-independent and coordinate-free method. This method is general and does not require a congruence of nearby beams' null geodesics as done in previous work. We show that the derivative of the parallel transport is the only means by which gravitational effects shows up in the frequency shift. This contribution is given as an integral of the Riemann tensor --the only physical observable of curvature-- along the beam's geodesic. The remaining contributions are: the difference of velocities, the difference of non-gravitational forces, and finally fictitious forces, either locally at the test masses or non-locally integrated along the beam's geodesic. As an application relevant to gravitational wave detection, we work out the frequency shift in the local Lorentz frame of nearby geodesics.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
Correlations between the star formation rates (SFRs) of nearby galaxies (so-called galactic conformity) have been observed for projected separations up to 4 Mpc, an effect not predicted by current semi-analytic models. We investigate correlations between the mass accretion rates (dMvir/dt) of nearby halos as a potential physical origin for this effect. We find that pairs of host halos "know about" each others' assembly histories even when their present-day separation is greater than thirty times the virial radius of either halo. These distances are far too large for direct interaction between the halos to explain the correlation in their dMvir/dt. Instead, halo pairs at these distances reside in the same large-scale tidal environment, which regulates dMvir/dt for both halos. Larger halos are less affected by external forces, which naturally gives rise to a mass dependence of the halo conformity signal. SDSS measurements of galactic conformity exhibit a qualitatively similar dependence on stellar mass, including how the signal varies with distance. Based on the expectation that halo accretion and galaxy SFR are correlated, we predict the scale-, mass- and redshift-dependence of large-scale galactic conformity, finding that the signal should drop to undetectable levels by z > 1. These predictions are testable with current surveys to z ~ 1; confirmation would establish a strong correlation between dark matter halo accretion rate and central galaxy SFR.
Newborn black holes in collapsing massive stars can be accompanied by a fallback disk. The accretion rate is typically super-Eddington and strong disk outflows are expected. Such outflows could be directly observed in some failed explosions of compact (blue supergiants or Wolf-Rayet stars) progenitors, and may be more common than long-duration gamma-ray bursts. Using an analytical model, we show that the fallback disk outflows produce blue UV-optical transients with a peak bolometric luminosity of ~10^(42-43) erg s^-1 (peak R-band absolute AB magnitudes of -16 to -18) and an emission duration of ~ a few to ~ 10 days. The spectra are likely dominated intermediate mass elements, but will lack much radioactive nuclei and iron-group elements. The above properties are broadly consistent with some of the rapid blue transients detected by Pan-STARRS and PTF. This scenario can be distinguished from alternative models using radio observations within a few years after the optical peak.
Connecting galaxies with their descendants (or progenitors) at different redshifts can yield strong constraints on galaxy evolution. Observational studies have historically selected samples of galaxies using a physical quantity, such as stellar mass, either above a constant limit or at a constant cumulative number density. Investigation into the efficacy of these selection methods has not been fully explored. Using a set of four semi-analytical models based on the output of the Millennium Simulation, we find that selecting galaxies at a constant number density (in the range $-4.3 < \log\ n\ [\mathrm{Mpc}^{-3}\ h^{3}] < -3.0$) is superior to a constant stellar mass selected sample, although it still has significant limitations. Recovery of the average stellar mass, stellar mass density and average star-formation rate is highly dependent on the choice of number density but can all be recovered to within $<50\%$ at the commonly employed choice of $\log\ n\ [\mathrm{Mpc}^{-3}\ h^{3}] = -4.0$, corresponding to $\log M_\odot / h \sim 11.2$ at $z=0$, but this increases at lower mass limits. We show that there is a large scatter between the location of a given galaxy in a rank ordering based on stellar mass between different redshifts. We find that the inferred velocity dispersion may be a better tracer of galaxy properties, although further investigation is warranted into simulating this property. Finally, we find that over large redshift ranges selection at a constant number density is more effective in tracing the progenitors of modern galaxies than vice-versa.
We hypothesize that at least some of the recently discovered class of calcium-rich gap transients are tidal detonation events of white dwarfs (WDs) by black holes (BHs) or possibly neutron stars. We show that the properties of the calcium-rich gap transients agree well with the predictions of the tidal detonation model. Under the predictions of this model, we use a follow-up X-ray observation of one of these transients, SN 2012hn, to place weak upper limits on the detonator mass of this system that include all intermediate-mass BHs (IMBHs). As these transients are preferentially in the stellar haloes of galaxies, we discuss the possibility that these transients are tidal detonations of WDs caused by random flyby encounters with IMBHs in dwarf galaxies or globular clusters. This possibility has been already suggested in the literature but without connection to the calcium-rich gap transients. In order for the random flyby cross-section to be high enough, these events would have to be occurring inside these dense stellar associations. However, there is a lack of evidence for IMBHs in these systems, and recent observations have ruled out all but the very faintest dwarf galaxies and globular clusters for a few of these transients. Another possibility is that these are tidal detonations caused by three-body interactions, where a WD is perturbed toward the detonator in isolated multiple star systems. We highlight a number of ways this could occur, even in lower-mass systems with stellar-mass BHs or neutron stars. Finally, we outline several new observational tests of this scenario, which are feasible with current instrumentation.
We study the stellar content in the tidal tails of three nearby merging galaxies, NGC 520, NGC 2623, and NGC 3256, using BVI imaging taken with the Advanced Camera for Surveys on board the Hubble Space Telescope. The tidal tails in all three systems contain compact and fairly massive young star clusters, embedded in a sea of diffuse, unresolved stellar light. We compare the measured colors and luminosities with predictions from population synthesis models to estimate cluster ages and find that clusters began forming in tidal tails during or shortly after the formation of the tails themselves. We find a lack of very young clusters ($\le 10$ Myr old), implying that eventually star formation shuts off in the tails as the gas is used up or dispersed. There are a few clusters in each tail with estimated ages that are older than the modeled tails themselves, suggesting that these may have been stripped out from the original galaxy disks. The luminosity function of the tail clusters can be described by a single power-law, $dN/dL \propto L^\alpha$, with $-2.6 < \alpha < -2.0$. We find a stellar age gradient across some of the tidal tails, which we interpret as a superposition of 1) newly formed stars and clusters along the dense center of the tail and 2) a sea of broadly distributed, older stellar material ejected from the progenitor galaxies.
Detailed characterization of exoplanets has begun to yield measurements of their atmospheric properties that constrain the planets' origins and evolution. For example, past observations of the dayside emission spectrum of the hot Jupiter WASP-12b indicated that its atmosphere has a high carbon-to-oxygen ratio (C/O $>$ 1), suggesting it had a different formation pathway than is commonly assumed for giant planets. Here we report a precise near-infrared transmission spectrum for WASP-12b based on six transit observations with the Hubble Space Telescope/Wide Field Camera 3. We bin the data in 13 spectrophotometric light curves from 0.84 - 1.67 $\mu$m and measure the transit depths to a median precision of 51 ppm. We retrieve the atmospheric properties using the transmission spectrum and find strong evidence for water absorption (7$\sigma$ confidence). This detection marks the first high-confidence, spectroscopic identification of a molecule in the atmosphere of WASP-12b. The retrieved 1$\sigma$ water volume mixing ratio is between $10^{-5}-10^{-2}$, which is consistent with C/O $>$ 1 to within 2$\sigma$. However, we also introduce a new retrieval parameterization that fits for C/O and metallicity under the assumption of chemical equilibrium. With this approach, we constrain C/O to $0.5^{+0.2}_{-0.3}$ at $1\,\sigma$ and rule out a carbon-rich atmosphere composition (C/O$>1$) at $>3\sigma$ confidence. Further observations and modeling of the planet's global thermal structure and dynamics would aid in resolving the tension between our inferred C/O and previous constraints. Our findings highlight the importance of obtaining high-precision data with multiple observing techniques in order to obtain robust constraints on the chemistry and physics of exoplanet atmospheres.
We report the discovery of 28 promising and a total of 58 new lens candidates from the CFHT Legacy Survey (CFHTLS) based on about 11 million classifications performed by citizen scientists as part of the first Space Warps lens search. The goal of the blind lens search was to identify lenses missed by lens finding robots (the RingFinder on galaxy scales and ArcFinder on group/cluster scales), which have been previously used to mine the CFHTLS for lenses. We compare some properties of lens samples detected by these algorithms to the SpaceWarps sample and found that they are broadly similar. The image separation distribution calculated from the SpaceWarps discovered sample shows that our previous constraints on the average density profile of the lens population are robust. Space Warps recovers about 60% of the known sample and the new candidates show a richer variety compared to the lenses found by the two robots. We find that analyzing only those classifications which are performed by the high power volunteers, Space Warps can achieve a detection rate of up to 75% for the known lens sample. Training and calibration of the performance of citizen scientists is crucial for the success of Space Warps. We also present the SIMCT pipeline, used for generating a sample of realistic simulated lensed images in the CFHTLS, and a sample of duds and false positives used in the training. Such a training sample has a legacy value for testing future lens finding algorithms. We make our training sample publicly available.
We present results from Subaru/FMOS near-infrared (NIR) spectroscopy of 118 star-forming galaxies at $z\sim1.5$ in the Subaru Deep Field. These galaxies are selected as [OII]$\lambda$3727 emitters at $z\approx$ 1.47 and 1.62 from narrow-band imaging. We detect H$\alpha$ emission line in 115 galaxies, [OIII]$\lambda$5007 emission line in 45 galaxies, and H$\beta$, [NII]$\lambda$6584, and [SII]$\lambda\lambda$6716,6731 in 13, 16, and 6 galaxies, respectively. Including the [OII] emission line, we use the six strong nebular emission lines in the individual and composite rest-frame optical spectra to investigate physical conditions of the interstellar medium in star-forming galaxies at $z\sim$1.5. We find a tight correlation between H$\alpha$ and [OII], which suggests that [OII] can be a good star formation rate (SFR) indicator for galaxies at $z\sim1.5$. The line ratios of H$\alpha$/[OII] are consistent with those of local galaxies. We also find that [OII] emitters have strong [OIII] emission lines. The [OIII]/[OII] ratios are larger than normal star-forming galaxies in the local Universe, suggesting a higher ionization parameter. Less massive galaxies have larger [OIII]/[OII] ratios. With evidence that the electron density is consistent with local galaxies, the high ionization of galaxies at high redshifts may be attributed to a harder radiation field by a young stellar population and/or an increase in the number of ionizing photons from each massive star.
Many high-state non-magnetic cataclysmic variables (CVs) exhibit blue-shifted absorption or P-Cygni profiles associated with ultraviolet (UV) resonance lines. These features imply the existence of powerful accretion disk winds in CVs. Here, we use our Monte Carlo ionization and radiative transfer code to investigate whether disk wind models that produce realistic UV line profiles are also likely to generate observationally significant recombination line and continuum emission in the optical waveband. We also test whether outflows may be responsible for the single-peaked emission line profiles often seen in high-state CVs and for the weakness of the Balmer absorption edge (relative to simple models of optically thick accretion disks). We find that a standard disk wind model that is successful in reproducing the UV spectra of CVs also leaves a noticeable imprint on the optical spectrum, particularly for systems viewed at high inclination. The strongest optical wind-formed recombination lines are H$\alpha$ and He II $\lambda4686$. We demonstrate that a higher-density outflow model produces all the expected H and He lines and produces a recombination continuum that can fill in the Balmer jump at high inclinations. This model displays reasonable verisimilitude with the optical spectrum of RW Trianguli. No single-peaked emission is seen, although we observe a narrowing of the double-peaked emission lines from the base of the wind. Finally, we show that even denser models can produce a single-peaked H$\alpha$ line. On the basis of our results, we suggest that winds can modify, and perhaps even dominate, the line and continuum emission from CVs.
The boundaries of cold dark matter halos are commonly defined to enclose a density contrast $\Delta$ relative to a reference (mean or critical) density. We argue that a more physical boundary of halos is the radius at which accreted matter reaches its first orbital apocenter after turnaround. This splashback radius, $R_{sp}$, manifests itself as a sharp density drop in the halo outskirts, at a location that depends upon the mass accretion rate. We present calibrations of $R_{sp}$ and the enclosed mass, $M_{sp}$, as a function of the accretion rate and alternatively peak height. We find that $R_{sp}$ varies between $\approx0.8-1R_{200m}$ for rapidly accreting halos and $\approx1.5R_{200m}$ for slowly accreting halos. The extent of a halo and its associated environmental effects can thus extend well beyond the conventionally defined "virial" radius. We show that $M_{sp}$ and $R_{sp}$ evolve relatively strongly compared to other commonly used definitions. In particular, $M_{sp}$ evolves significantly even for the smallest dwarf-sized halos at $z=0$. We also contrast $M_{sp}$ with the mass enclosed within four scale radii of the halo density profile, $M_{<4rs}$, which characterizes the inner halo. During the early stages of halo assembly, $M_{sp}$ and $M_{<4rs}$ evolve similarly, but in the late stages $M_{<4rs}$ stops increasing while $M_{sp}$ continues to grow significantly. This illustrates that halos at low $z$ can have "quiet" interiors while continuing to accrete mass in their outskirts. We discuss potential observational estimates of the splashback radius and show that it may already have been detected in galaxy clusters.
We present basic properties of $\sim$3,300 emission line galaxies detected by the FastSound survey, which are mostly H$\alpha$ emitters at $z \sim$ 1.2-1.5 in the total area of about 20 deg$^2$, with the H$\alpha$ flux sensitivity limit of $\sim 1.6 \times 10^{-16} \rm erg \ cm^{-2} s^{-1}$ at 4.5 sigma. This paper presents the catalogs of the FastSound emission lines and galaxies, which will be open to the public in the near future. We also present basic properties of typical FastSound H$\alpha$ emitters, which have H$\alpha$ luminosities of $10^{41.8}$-$10^{43.3}$ erg/s, SFRs of 20--500 $M_\odot$/yr, and stellar masses of $10^{10.0}$--$10^{11.3}$ $M_\odot$. The 3D distribution maps for the four fields of CFHTLS W1--4 are presented, clearly showing large scale clustering of galaxies at the scale of $\sim$ 100--600 comoving Mpc. Based on 1,105 galaxies with detections of multiple emission lines, we estimate that contamination of non-H$\alpha$ lines is about 4% in the single-line emission galaxies, which are mostly [OIII]$\lambda$5007. This contamination fraction is also confirmed by the stacked spectrum of all the FastSound spectra, in which H$\alpha$, [NII]$\lambda \lambda$6548,6583, [SII]$\lambda \lambda$6717, 6731, and [OI]$\lambda \lambda$6300,6364 are seen.
We assess the fraction of the Milky Way's circumgalactic medium (CGM) eluding detection due to its velocity being similar to gas in the disk. This is achieved using synthetic observations of the CGM in a simulated MW-mass galaxy that shows similar CGM kinematics to the MW and external L$\sim$L$_*$ galaxies. As viewed by a mock observer at a location similar to the Sun, only 50$\%$ (by mass) of the gas moves at high velocity ($|v_{\rm LSR}|\geq$100 km s$^{-1}$ or $|v_{\rm DEV}|\geq$50 km s$^{-1}$) in the simulated CGM and would be observable. The low velocity gas is thermodynamically similar to the high velocity gas, indicating the 50$\%$ observable fraction is applicable to spectral lines from the radio to the ultraviolet. We apply the observable mass fraction (50$\%$) to current estimates of the MW's CGM, and find a corrected total mass of 2.8$\times$10$^{8} M_{\odot}$ for gas below 10$^6$K within $\sim15$ kpc (this excludes the Magellanic System). This is less than the total mass of the CGM extending out to $\sim$150 kpc in other L$\sim$L$_*$ galaxies. However, we find similar OVI column densities when the discrepancy in path length between the MW and external galaxies is considered. The coherent spatial and kinematic distribution of low velocity gas in the simulated CGM suggests that current HI observations of the MW's CGM may miss large low velocity HI complexes. In addition, current mass estimates of the MW's CGM based on high-velocity line observations with distance constraints may miss a non-negligible fraction of gas in the outer halo which can be obscured if it moves at a velocity similar to the gas in the lower halo.
PKS 0521-36 is an Active Galactic Nucleus (AGN) with uncertain classification. We investigate the properties of this source from radio to gamma rays. The broad emission lines in the optical and UV bands and steep radio spectrum indicate a possible classification as an intermediate object between broad-line radio galaxies (BLRG) and steep spectrum radio quasars (SSRQ). On pc-scales PKS 0521-36 shows a knotty structure similar to misaligned AGN. The core dominance and the gamma-ray properties are similar to those estimated for other SSRQ and BLRG detected in gamma rays, suggesting an intermediate viewing angle with respect to the observer. In this context the flaring activity detected from this source by Fermi-LAT between 2010 June and 2012 February is very intriguing. We discuss the gamma-ray emission of this source in the framework of the structured jet scenario, comparing the spectral energy distribution (SED) of the flaring state in 2010 June with that of a low state. We present three alternative models corresponding to three different choices of the viewing angles theta_v = 6, 15, and 20 deg. We obtain a good fit for the the first two cases, but the SED obtained with theta_v = 15 deg if observed at a small angle does not resemble that of a typical blazar since the synchrotron emission should dominate by a large factor (about 100) the inverse Compton component. This suggests that a viewing angle between 6 and 15 deg is preferred, with the rapid variability observed during gamma-ray flares favouring a smaller angle. However, we cannot rule out that PKS 0521-36 is the misaligned counterpart of a synchrotron-dominated blazar.
We cross-correlate cosmic microwave background (CMB) lensing and galaxy weak lensing maps using the Planck 2013 and 2015 data and the 154 deg^2 Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). This measurement probes large-scale structure at intermediate redshifts ~0.9, between the high- and low-redshift peaks of the CMB and CFHTLenS lensing kernels, respectively. Using the noise properties of these data sets and standard Planck 2015 LCDM cosmological parameters, we forecast a signal-to-noise ratio ~4.6 for the cross-correlation. We find that the noise level of our actual measurement agrees well with this estimate, but the amplitude of the signal lies well below the theoretical prediction. The best-fit amplitudes of our measured cross-correlations are $A_{2013}=0.48\pm0.26$ and $A_{2015}=0.44\pm0.22$ using the 2013 and 2015 Planck CMB lensing maps, respectively, where $A=1$ corresponds to the fiducial Planck 2015 LCDM prediction. Due to the low measured amplitude, the detection significance is moderate (~2$\sigma$) and the data are in tension with the theoretical prediction (~2-$2.5\sigma$). The tension is reduced somewhat when compared to predictions using WMAP9 parameters, for which we find $A_{2013}=0.56\pm0.30$ and $A_{2015}=0.52\pm0.26$. We consider various systematic effects, finding that photometric redshift uncertainties, contamination by intrinsic alignments, and effects due to the masking of galaxy clusters in the Planck 2015 CMB lensing reconstruction are able to help resolve the tension at a significant level (~10% each). An overall multiplicative bias in the CFHTLenS shear data could also play a role, which can be tested with existing data. We close with forecasts for measurements of this cross-correlation using ongoing and future weak lensing surveys, which will definitively test the significance of the tension in our results with respect to LCDM.
Details of various unknown physical processes during the cosmic dawn and the epoch of reionization can be extracted from observations of the redshifted 21-cm signal. These observations, however, will be affected by the evolution of the signal along the line-of-sight which is known as the "light-cone effect". We model this effect by post-processing a dark matter $N-$body simulation with a 1-D radiative transfer code. We find that the effect is much stronger and dramatic in presence of inhomogeneous heating and Ly$\alpha$ coupling compared to the case where these processes are not accounted for. One finds increase (decrease) in the coeval spherically averaged power spectrum up to a factor of 3 (0.6) at large scales ($k \sim 0.05\, \rm Mpc^{-1}$), though these numbers are highly dependent on the source model. Consequently, the peak and trough-like features seen in the evolution of the large-scale power spectrum can be smoothed out to a large extent if the width of the frequency bands used in the experiment is large. We argue that it is important to account for the light-cone effect for any 21-cm signal prediction during cosmic dawn.
We present a study of the spatial distribution and kinematics of star-forming galaxies in 30 massive clusters at 0.15<z<0.30, combining wide-field Spitzer 24um and GALEX NUV imaging with highly-complete spectroscopy of cluster members. The fraction (f_SF) of star-forming cluster galaxies rises steadily with cluster-centric radius, increasing fivefold by 2r200, but remains well below field values even at 3r200. This suppression of star formation at large radii cannot be reproduced by models in which star formation is quenched in infalling field galaxies only once they pass within r200 of the cluster, but is consistent with some of them being first pre-processed within galaxy groups. Despite the increasing f_SF-radius trend, the surface density of star-forming galaxies actually declines steadily with radius, falling ~15x from the core to 2r200. This requires star-formation to survive within recently accreted spirals for 2--3Gyr to build up the apparent over-density of star-forming galaxies within clusters. The velocity dispersion profile of the star-forming galaxy population shows a sharp peak of 1.44-sigma_v at 0.3r500, and is 10--35% higher than that of the inactive cluster members at all cluster-centric radii, while their velocity distribution shows a flat, top-hat profile within r500. All of these results are consistent with star-forming cluster galaxies being an infalling population, but one that must also survive ~0.5--2Gyr beyond passing within r200. By comparing the observed distribution of star-forming galaxies in the stacked caustic diagram with predictions from the Millennium simulation, we obtain a best-fit model in which SFRs decline exponentially on quenching time-scales of 1.73\pm0.25 Gyr upon accretion into the cluster.
We present an analysis of the impact of the tree rings seen in the candidate sensors of the Large Synoptic Survey Telescope (LSST) on galaxy-shape measurements. The tree rings are a consequence of transverse electric fields caused by circularly symmetric impurity gradients in the silicon of the sensors. They effectively modify the pixel area and shift the photogenerated charge around, displacing the observed photon positions. The displacement distribution generates distortions that cause spurious shears correlated with the tree-rings patterns, potentially biasing cosmic shear measurements. In this paper we quantify the amplitude of the spurious shear caused by the tree rings on the LSST candidate sensors, and calculate its 2-point correlation function. We find that 2-point correlation function of the spurious shear on an area equivalent to the LSST field of view is order of about $10^{-13}$, providing a negligible contribution to the 2-point correlation of the cosmic shear signal. Additional work is underway, and the final results and analysis will be published elsewhere (Okura et al. (2015), in prep.)
A sample of 111 Fermi blazars each with a well-established radio core luminosity, broad-line luminosity, bolometric luminosity and black hole mass has been compiled from the literatures.We present a significant correlation between radio core and broad-line emission luminosities that supports a close link between accretion processes and relativistic jets. Analysis reveals a relationship of $\rm{LogL_{BLR}\sim(0.81\pm0.06)LogL_{R}^{C}}$ which is consistant with theoretical predicted coefficient and supports that blazar jets are powered by energy extraction from a rapidly spinning Kerr black hole through the magnetic field provided by the accretion disk. Through studying the correlation between the intrinsic bolometric luminosity and the black hole mass, we find a relationship of $\rm{{Log}\frac{L_{in}}{L_{\odot}}=(0.95\pm0.26){Log}\frac{M}{M_{\odot}}+(3.53\pm2.24)}$ which supports mass-luminosity relation for Fermi blazars derived in this work is a powerlaw relation similar to that for main-sequence stars. Finally, evolutionary sequence of blazars is discussed.
We compiled a sample of 24 flat-spectrum radio quasars (FSRQs), 21 BL Lacs, 13 FRI and 12 FRII radio galaxies to study the jet power and broad-line luminosity relation, and the jet power and black hole mass relation. Furthermore, we obtain the histograms of key parameters. Our main results are as the follows:(i) We found that the FSRQs are roughly separated from BL Lacs by the Ledlow-Owen FRI/FRII dividing line in the $\rm{\log{P_{jet}}-\log{M}}$ plane. This result supports the unification model that BL Lacs unify with FRI radio galaxies and FSRQs unify with FRII radio galaxies.(ii)We found that the coefficient of the best-fit linear regression equation of $\rm{\log{L_{BLR}}-\log{P_{jet}}}$ relation is very close to 1 for our sample. The correlation between broad line luminosity and jet power is significant which supports that jet power has a close link with accretion.
It is thought that planetary mass companions may form through gravitational disk instabilities or core accretion. Identifying such objects in the process of formation would provide the most direct test for the competing formation theories. One of the most promising candidates for a planetary mass object still in formation is the third object in the FWTau system. We here present ALMA cycle 1 observations confirming the recently published 1.3 mm detection of a dust disk around this third object and present for the first time a clear detection of a single peak 12CO(2-1) line, providing direct evidence for the simultaneous existence of a gas disk. We perform radiative transfer modeling of the third object in FW Tau and find that current observations are consistent with a planetary mass object embedded in a disk which is externally irradiated by the binary companion and seen at an inclination of i<15 deg. However, we also find that a near edge-on disk around a more massive substellar object can explain the observations if cloud contamination causes the single peak shape of the 12CO(2-1) line. Although this possibility appears less likely, further observations with ALMA, aiming for the detection of less contaminated gas lines, are required to conclusively unveil the nature of the third object in FWTau.
Weak gravitational lensing enables us to search clusters without the conventional assumption on the relation between visible and dark matter. We explore a variety of statistics of clusters selected with cosmic shear measurement by utilizing both analytic models and large numerical simulations. We first develop a halo model to predict the abundance and the clustering of weak lensing selected clusters. Observational effects such as galaxy shape noise are included in our model. We then generate realistic mock weak lensing catalogs to test the accuracy of our analytic model. To this end, we perform full-sky ray-tracing simulations that allow us to have multiple realizations of a large continuous area. We model the masked regions on the sky using the actual positions of bright stars, and generate 200 mock weak lensing catalogs with sky coverage of $\sim$1000 squared degrees. We utilize the large set of mock catalogs to evaluate the covariance matrices between the local and non-local statistics. We show that our theoretical model agrees well with the ensemble average of statistics and their covariances calculated directly from the mock catalogues. With a typical selection threshold, ignoring shape noise correction causes overestimation of the clustering of weak lensing selected clusters with a level of about $10\%$, and shape noise correction boosts the cluster abundance by a factor of a few. We calculate the cross-covariances using the halo model with accounting for the effective reduction of the survey area due to masks. The covariance of the cosmic shear auto power spectrum is affected by the mode-coupling effect that originates from sky masking. Our model and the results can be readily used for cosmological analysis with ongoing and future weak lensing surveys.
Magnetic reconnection is an important physical process in various explosive phenomena in the universe. In the previous studies, it was found that fast re- connection takes place when the thickness of a current sheet becomes on the order of a microscopic length such as the ion larmor radius or the ion inertial length. In this study, we investigated the pinching process of a current sheet by the Lorentz force in a low-{\beta} plasma using one-dimensional magnetohydrodynam- ics (MHD) simulations. It is known that there is an exact self-similar solution for this problem that neglects gas pressure. We compared the non-linear MHD dynamics with the analytic self-similar solution. From the MHD simulations, we found that with the gas pressure included the implosion process deviates from the analytic self-similar solution as t {\rightarrow} t 0, where t 0 is the explosion time when the thickness of a current sheet of the analytic solution becomes 0. We also found a pair of MHD fast-mode shocks are generated and propagate after the formation of the pinched current sheet as t {\rightarrow} t 0 . On the basis of the Rankine-Hugoniot relations, we derived the scaling law of the physical quantities with respect to the initial plasma beta in the pinched current sheet. Our study could help us to estimate the physical quantities in the pinched current sheet formed in a low-\b{eta} plasma.
We report spatial fluctuation analysis of the sky brightness in near-infrared from observations toward the north ecliptic pole (NEP) by the AKARI at 2.4 and 3.2 micron. As a follow up study of our previous work on the Monitor field of AKARI, we used NEP deep survey data, which covered a circular area of about 0.4 square degrees, in order to extend fluctuation analysis at angular scales up to 1000". We found residual fluctuation over the estimated shot noise at larger angles than the angular scale of the Monitor field. The excess fluctuation of the NEP deep field smoothly connects with that of the Monitor field at angular scales with a few hundreds arcseconds and extends without any significant variation to larger angular scales up to 1000". By comparing excess fluctuations at two wavelengths, we confirm a blue spectrum feature similar to the result of the Monitor field. We find that the result of this study is consistent with Spitzer Space Telescope observations at 3.6 micron. The origin of the excess fluctuation in the near-infrared background remains to be answered, but we could exclude zodiacal light, diffuse Galactic light, and unresolved faint galaxies at low-redshift based on the comparison with mid- and far-infrared brightness, ground based near-infrared images.
The FORS2 instrument is one of the most widely used and productive instruments on the Very Large Telescope. This article reports on a project to improve the quality of the reduced FORS2 spectra that can be produced with the software provided by ESO. The result of this effort is that spectra of significantly higher quality can now be produced with substantially lower effort by the science user of the data.
Solar flares are an explosive phenomenon, where super-sonic flows and shocks are expected in and above the post-flare loops. To understand the dynamics of post-flare loops, a two-dimensional magnetohydrodynamic (2D MHD) simulation of a solar flare has been carried out. We found new shock structures in and above the post-flare loops, which were not resolved in the previous work by Yokoyama and Shibata 2001. To study the dynamics of flows along the reconnected magnetic field, kinematics and energetics of the plasma are investigated along selected field lines. It is found that shocks are crucial to determine the thermal and flow structures in the post-flare loops. On the basis of the 2D MHD simulation, we have developed a new post-flare loop model which we call the pseudo-2D MHD model. The model is based on the 1D MHD equations, where all the variables depend on one space dimension and all the three components of the magnetic and velocity fields are considered. Our pseudo-2D model includes many features of the multi-dimensional MHD processes related to magnetic reconnection (particularly MHD shocks), which the previous 1D hydrodynamic models are not able to include. We compare the shock formation and energetics of a specific field line in the 2D calculation with those in our pseudo-2D MHD model, and we found that they give similar results. This model will allow us to study the evolution of the post-flare loops in a wide parameter space without expensive computational cost and without neglecting important physics associated with magnetic reconnection.
The 13 Myr old star HD106906 is orbited by a debris disk of at least 0.067 M_Moon with an inner and outer radius of 20 AU and 120 AU, respectively, and by a planet at a distance of 650 AU. We use this curious combination of a close low-mass disk and a wide planet to motivate our simulations of this system. We study the parameter space of the initial conditions to quantify the mass loss from the debris disk and its lifetime under the influence of the planet. We find that when the planet orbits closer to the star than about 50 AU and with low inclination relative to the disk (less than about 10 degrees), more disk material is perturbed outside than inside the region constrained by observations on timescales shorter than 1 Myr. Considering the age of the system, such a short lifetime of the disk is incompatible with the timescale for planet--planet scattering which is one of the scenarios suggested to explain the wide separation of the planet. For some configurations when the planet's orbit is inclined with respect to the disk, the latter will start to wobble. We argue that this wobbling is caused by a mechanism similar to the Kozai--Lidov oscillations. We also observe various resonant structures (such as rings and spiral arms) induced in the disk by the planet.
We report the detection during the ALMA Cycle 0 of SiS rotational lines in high-vibrational states as well as SiO and SiC$_2$ lines in their ground vibrational state, towards IRC+10216. The spatial distribution of these molecules shows compact emission for SiS and a more extended emission for SiO and SiC$_2$ , and also proves the existence of an increase in the SiC$_2$ emission at the outer shells of the circumstellar envelope. We analyze the excitation conditions of the vibrationally excited SiS using the population diagram technique and we used a large velocity gradient model to compare with the observations. We found moderate discrepancies between the observations and the models that could be explained if SiS lines detected are optically thick. Additionally, the line profiles of the detected rotational lines in the high energy vibrational states show a decreasing linewidth with increasing energy levels. This may evidence that these lines could be excited only in the inner shells, i.e. the densest and hottest, of the circumstellar envelope of IRC+10216.
Context. Extremely reddened AGB stars lose mass at high rates of >10^-5 Msun/yr. This is the very last stage of AGB evolution, in which stars in the mass range 2.0--4.0 Msun (for solar metallicity) should have been converted to C stars already. The extremely reddened AGB stars in the Galactic bulge are however predominantly O-rich, implying that they might be either low-mass stars or stars at the upper end of the AGB mass range. Aims. To determine the mass range of the most reddened AGB stars in the Galactic bulge. Methods. Using Virtual Observatory tools, we constructed spectral energy distributions of a sample of 37 evolved stars in the Galactic bulge with extremely red IRAS colours. We fitted DUSTY models to the observational data to infer the bolometric fluxes. Applying individual corrections for interstellar extinction and adopting a common distance, we determined luminosities and mass-loss rates, and inferred the progenitor mass range from comparisons with AGB evolutionary models. Results. The observed spectral energy distributions are consistent with a classification as reddened AGB stars, except for two stars, which are proto-planetary nebula candidates. For the AGB stars, we found luminosities in the range 3000--30,000 Lsun and mass-loss rates 10^-5--3x10^-4 Msun/yr. The corresponding mass range is 1.1--6.0 Msun assuming solar metallicity. Conclusions. Contrary to the predictions of the evolutionary models, the luminosity distribution is continuous, with many O-rich AGB stars in the mass range in which they should have been converted into C stars already. We suspect that bulge AGB stars have higher than solar metallicity and therefore may avoid the conversion to C-rich. The presence of low-mass stars in the sample shows that their termination of the AGB evolution also occurs during a final phase of very high mass-loss rate, leading to optically thick circumstellar shells.
The gamma-ray observations of molecular clouds associated with supernova remnants are considered one of the most promising ways to search for a solution of the problem of cosmic ray origin. Here we briefly review the status of the field, with particular emphasis on the theoretical and phenomenological aspects of the problem.
We have conducted radioastronomical observations of 9 dark clouds with the IRAM 30m telescope. We present the first identification in space of the ketenyl radical (HCCO) toward the starless core Lupus-1A and the molecular cloud L483, and the detection of the related molecules ketene (H2CCO) and acetaldehyde (CH3CHO) in these two sources and 3 additional dark clouds. We also report the detection of the formyl radical (HCO) in the 9 targeted sources and of propylene (CH2CHCH3) in 4 of the observed sources, which extends significantly the number of dark clouds where these molecules are known to be present. We derive a beam-averaged column density of HCCO of 5e11 cm-2 in both Lupus-1A and L483, which means that the ketenyl radical is just 10 times less abundant than ketene in these sources. The non-negligible abundance of HCCO found implies that there must be a powerful formation mechanism able to counterbalance the efficient destruction of this radical through reactions with neutral atoms. The column densities derived for HCO, (0.5-2.7)e12 cm-2, and CH2CHCH3, (1.9-4-2)e13 cm-2, are remarkably uniform across the sources where these species are detected, confirming their ubiquity in dark clouds. Gas phase chemical models of cold dark clouds can reproduce the observed abundances of HCO, but cannot explain the presence of HCCO in Lupus-1A and L483 and the high abundances derived for propylene. The chemistry of cold dark clouds needs to be revised in the light of these new observational results.
Synchrotron diffuse radiation (SDR) emission is one of the major Galactic components, in the 100 MHz up to 100 GHz frequency range. Its spectrum and sky map provide valuable measure of the galactic cosmic ray electrons (GCRE) in the relevant energy range, as well as of the strength and structure of the Galactic magnetic fields (GMF), both regular and random ones. This emission is an astrophysical sky foreground for the study of the Cosmic Microwave Background (CMB), and the extragalactic microwave measurements, and it needs to be modelled as better as possible. In this regard, in order to get an accurate description of the SDR in the Galaxy, we use - for the first time in this context - 3-dimensional GCRE models obtained by running the DRAGON code. This allows us to account for a realistic spiral arm pattern of the source distribution, demanded to get a self-consistent treatment of all relevant energy losses influencing the final synchrotron spectrum.
We consider a higher order term in the $\delta N$ expansion for the CMB power asymmetry generated by a superhorizon isocurvature field fluctuation. The term can generate the asymmetry without requiring a large value of $f_{NL}$. Instead it produces a non-zero value of $g_{NL}$. A combination of constraints lead to an allowed region in $f_{NL}-g_{NL}$ space. To produce the asymmetry with this term without a large value of $f_{NL}$ we find that the isocurvature field needs to contribute less than the inflaton towards the power spectrum of the curvature perturbation.
We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields.
The search for life via characterization of earth-like planets in the habitable zone is one of the key scientific objectives in Astronomy. We describe a new phase-occulting (PO) interferometric nulling coronagraphy (NC) approach. The PO-NC approach employs beamwalk and freeform optical surfaces internal to the interferometer cavity to introduce a radially dependent plate scale difference between each interferometer arm (optical path) that nulls the central star at high contrast while transmitting the off-axis field. The design is readily implemented on segmented-mirror telescope architectures, utilizing a single nulling interferometer to achieve high throughput, a small inner working angle (IWA), sixth-order or higher starlight suppression, and full off-axis discovery space, a combination of features that other coronagraph designs generally must trade. Unlike previous NC approaches, the PO-NC approach does not require pupil shearing; this increases throughput and renders it less sensitive to on-axis common-mode telescope errors, permitting relief of the observatory stability required to achieve contrast levels of $\leq10^{-10}$. Observatory operations are also simplified by removing the need for multiple telescope rolls and shears to construct a high contrast image. The design goals for a PO nuller are similar to other coronagraphs intended for direct detection of habitable zone (HZ) exoEarth signal: contrasts on the order of $10^{-10}$ at an IWA of $\leq3\lambda/D$ over $\geq10$% bandpass with a large ($>10$~m) segmented aperture space-telescope operating in visible and near infrared bands. This work presents an introduction to the PO nulling coronagraphy approach based on its Visible Nulling Coronagraph (VNC) heritage and relation to the radial shearing interferometer.
We report the discovery of a new group of double-periodic RR Lyrae stars from the analysis of the OGLE-IV Galactic bulge photometry. In 11 stars identified in the OGLE catalog as first overtone pulsators (RRc stars) we detect additional longer period variability of low amplitude, in the mmag regime. One additional star of the same type is identified in a published analysis of the Kepler space photometry. The period ratio between the shorter first overtone period and a new, longer period lies in a narrow range around 0.686. Thus, the additional period is longer than the expected period of the undetected radial fundamental mode. The obvious conclusion that addition periodicity corresponds to a gravity or a mixed mode faces difficulties, however.
A new sample of stars, representative of the solar neighbourhood luminosity function, is constructed from the Hipparcos catalogue and the Fifth Catalogue of Nearby Stars. We have cross-matched to sources in the 2MASS catalogue so that for all stars individually determined Near Infrared photometry (NIR) is available on a homogeneous system (typically K_s). The spatial completeness of the sample has been carefully determined by statistical methods, and the NIR luminosity function of the stars has been derived by direct star counts. We find a local volume luminosity of 0.121 +/- 0.004 L_K_sun/(pc**3), corresponding to a volumetric mass-to-light ratio of M/L_K = 0.31 +/- 0.02 M_sun/L_K_sun, where giants contribute 80 per cent to the light but less than 2 per cent to the stellar mass. We derive the surface brightness of the solar cylinder with the help of a vertical disc model. We find a surface brightness of 99 L_K_sun/(pc**2) with an uncertainty of approximately 10 %. This corresponds to a mass-to-light ratio for the solar cylinder of M/L_K = 0.34 M_sun/L_K_sun. The mass-to-light ratio for the solar cylinder is only 10% larger than the local value despite the fact that the local population has a much larger contribution of young stars. It turns out that the effective scale heights of the lower main sequence carrying most of the mass is similar to that of the giants, which are dominating the NIR light. The corresponding colour for the solar cylinder is V-K=2.89 mag compared to the local value of V-K = 2.46 mag. An extrapolation of the local surface brightness to the whole Milky Way yields a total luminosity of M_K = -24.2 mag. The Milky Way falls in the range of K band Tully-Fisher (TF) relations from the literature.
We considered 18 solar flares observed between June 2010 and July 2012, in which high energy >100 MeV {\gamma}-emission was registered by the Large Area Telescope (LAT) aboard FermiGRO. We examined for these {\gamma}-events soft X-ray observations by GOES, hard X-ray observations by the Anti-Coincidence Shield of the SPectrometer aboard INTEGRAL (ACS SPI) and the Gamma-Ray burst Monitor (GBM) aboard FermiGRO. Hard X-ray and {\pi}0-decay {\gamma}-ray emissions are used as tracers of electron and proton acceleration, respectively. Bursts of hard X-ray were observed by ACS SPI during impulsive phase of 13 events. Bursts of hard X-ray >100 keV were not found during time intervals, when prolonged hard {\gamma}-emission was registered by LAT/FermiGRO. Those events showing prolonged high-energy gamma-ray emission not accompanied by >100 keV hard X-ray emission are interpreted as an indication of either different acceleration processes for protons and electrons or as the presence of a proton population accelerated during the impulsive phase of the flare and subsequently trapped by some magnetic structure. In-situ energetic particle measurements by GOES and STEREO (High Energy Telescope, HET) shows that five of these {\gamma}-events were not accompanied by SEP events at 1 AU, even when multi-point measurements including STEREO are taken into account. Therefore accelerated protons are not always released into the heliosphere. A longer delay between the maximum temperature and the maximum emission measure characterises flares with prolonged high energy {\gamma}-emission and solar proton events.
A simple method of the vertical muon energy spectrum simulations have been suggested. These calculations have been carried out in terms of various models of hadronic interactions. The most energetic $ \pi^\pm $-mesons and K$^\pm $-mesons produced in hadron interactions contribute mainly in to this energy spectrum of muons due to the very steep energy spectrum of the primary particles. So, some constraints on the hadronic interaction models may be set from a comparison of calculated results with the cosmic data on the vertical muon energy spectrum. This comparison showed that the most energetic secondary particles production is too high in case of the QGSJET II-04 model and rather low in case of the QGSJET II-03 model. These conclusion have been supported by the LHC data.
The MIMAC experiment is a $\mu$-TPC matrix project for directional dark matter search. Directional detection is a strategy based on the measurement of the WIMP flux anisotropy due to the solar system motion with respect to the dark matter halo. The main purpose of MIMAC project is the measurement of the energy and the direction of nuclear recoils in 3D produced by elastic scattering of WIMPs. Since June 2012 a bi-chamber prototype is operating at the Modane underground laboratory. In this paper, we report the first ionization energy and 3D track observations of nuclear recoils produced by the radon progeny. This measurement shows the capability of the MIMAC detector and opens the possibility to explore the low energy recoil directionality signature.
We present Atacama Large Millimeter Array (ALMA) detections of atomic carbon line and dust continuum emission in two UV-luminous galaxies at redshift 6. The far-infrared (FIR) luminosities of these galaxies are substantially lower than similar starbursts at later cosmic epochs, indicating an evolution in the dust properties with redshift, in agreement with the evolution seen in ultraviolet (UV) attenuation by dust. The [CII] to FIR ratios are found to be higher than at low redshift showing that [CII] should be readily detectable by ALMA within the reionization epoch. One of the two galaxies shows a complex merger nature with the less massive component dominating the UV emission and the more massive component dominating the FIR line and continuum. Using the interstellar atomic carbon line to derive the systemic redshifts we investigate the velocity of Lyman alpha emission emerging from high-z galaxies. In contrast to previous work, we find no evidence for decreasing Lyman alpha velocity shifts at high-redshift. We observe an increase in velocity shifts from z~ 2 to z~6, consistent with the effects of increased IGM absorption.
Nuclear star clusters are among the densest stellar systems known and are common in both early- and late-type galaxies. They exhibit scaling relations with their host galaxy which may be related to those of supermassive black holes. These may therefore help us to unravel the complex physical processes occurring at the centres of galaxies. The properties of nuclear stellar systems suggest that their formation requires both dissipational and dissipationless processes. They have stellar populations of different ages, from stars as old as their host galaxy to young stars formed in the last 100 Myr. Therefore star formation must be happening either directly in the nuclear star cluster or in its vicinity. The secular processes that fuel the formation of pseudobulges very likely also contributes to nuclear star cluster growth.
We examine the combined effects of winds and photoionizing radiation from O--type stars on embedded stellar clusters formed in model turbulent molecular clouds covering a range of masses and radii. We find that feedback is able to increase the quantities of dense gas present, but decreases the rate and efficiency of the conversion of gas to stars relative to control simulations in which feedback is absent. Star formation in these calculations often proceeds at a rate substantially slower than the freefall rate in the dense gas. This decoupling is due to the weakening of, and expulsion of gas from, the deepest parts of the clouds' potential wells where most of the star formation occurs in the control simulations. This results in large fractions of the stellar populations in the feedback simulation becoming dissociated from dense gas. However, where star formation \emph{does} occur in both control and feedback simulations, it does so in dense gas, so the correlation between star formation activity and dense gas is preserved. The overall dynamical effects of feedback on the \emph{clusters} are minimal, with only small fraction of stars becoming unbound, despite large quantities of gas being expelled from some clouds. This owes to the settling of the stars into virialised and stellar--dominated configurations before the onset of feedback. By contrast, the effects of feedback on the observable properties of the clusters -- their U--, B-- and V--band magnitudes -- are strong and sudden. The timescales on which the clusters become visible and unobscured are short compared with the timescales which the clouds are actually destroyed.
We present correlations between 9 CO transition ($J=4-3$ to $12-11$) lines and beam-matched far-infrared (FIR) luminosities ($L_{\mathrm{FIR},\,b}$) among 167 local galaxies, using Herschel SPIRE/FTS spectroscopic data and PACS photometry data. We adopt entire-galaxy FIR luminosities ($L_{\mathrm{FIR},\,e}$) from the {\it{IRAS}} Revised Bright Galaxy Sample and correct to $L_{\mathrm{FIR},\,b}$ using PACS images to match the varying FTS beams. All 9 correlations between $L'_{\mathrm{CO}}$ and $L_{\mathrm{FIR},\,b}$ are essentially linear, even for the highest transition $J=12-11$. This supports the notion that dense molecular gas ($n_{\mathrm{H}_2}\gtrsim10^{4-6}\,cm^{-3}$) linearly correlates with the star formation rate (SFR). We divide the entire sample into three subsamples and find that smaller sample size can induce large difference in the correlation slopes. We also derive an average CO spectral line energy distribution (SLED) for the entire sample and discuss the implied average molecular gas properties for these local galaxies. We further extend our sample to high-z galaxies with literature CO($J=5-4$) data as an example, including submillimeter galaxies (SMGs) and "normal" star-forming BzKs. BzKs have similar FIR/CO(5-4) ratios as that of local galaxies, well following the local correlation, whereas SMGs roughly distribute around or slightly above local correlation with large uncertainties. Finally, by using Galactic CO($J=10-9$) data as well as very limited high-z CO($J=10-9$) data, we verify that the CO($J=10-9$) -- FIR correlation successfully extends to Galactic YSOs, suggesting that linear correlations are most likely valid over nearly 15 orders of magnitude.
The Active Particle-induced X-ray Spectrometer (APXS) is one of the payloads on board the Yutu rover of Chang'E-3 mission. In order to assess the instrumental performance of APXS, a ground verification test was done for two unknown samples (basaltic rock, mixed powder sample). In this paper, the details of the experiment configurations and data analysis method are presented. The results show that the elemental abundance of major elements can be well determined by the APXS with relative deviations < 15 wt. % (detection distance = 30 mm, acquisition time = 30 min). The derived detection limit of each major element is inversely proportional to acquisition time and directly proportional to detection distance, suggesting that the appropriate distance should be < 50mm.
We study the vertical gradient in azimuthal velocity of spiral galaxy NGC 4244 in a thin disk model. With surface density accounting for the rotation curve, we model the gradient properties in the approximation of quasi-circular orbits and find the predictions to be consistent with the gradient properties inferred from measurements. This consistency may suggest that the mass distribution in this galaxy is flattened.
In this article we investigate the outer and inner mass distributions of the irregular galaxies UGC 4284 and UGC 11861, taking advantage of published HI and H{\alpha} high resolution rotation curves and constraining the stellar disk of both galaxies throughout stellar population synthesis studies. In addition we take into account the gas content of both galaxies deriving the HI+He rotation curve. The deduced baryonic rotation curves (star+gas) are inadequate to account for the total mass of UGC 4284 and UGC 11861, for that reason we examine the possibility of dark matter to explain the incongruity between the observed HI and H{\alpha} rotation curves of UGC 4284 and UGC 11861 and the derived baryonic rotation curves. We consider NFW, Burkert, DiCintio, Einasto, and the Stadel dark matter halos, to analyse the dark matter content of UGC 4284 and UGC 11861. The principal results of this work are that cored dark matter models better reproduce the dark matter H{\alpha} and HI rotation curves of UGC 11861 and the dark matter HI rotation curve of UGC 4284, while, the H{\alpha} rotation curve of UGC 4284 is better reproduced by a cuspy DiCintio DM model. In general, cored exponential two-parameters models Einasto and Stadel, give better fits than Burkert. This trend, as well as to confirm past results, presents for the first time a comparison between two different exponential dark matter models, Einasto and Stadel, in an attempt to better constrain the range of possible exponential dark matter models applied to real galaxies.
Scale-dependent halo bias due to local primordial non-Gaussianity provides a strong test of single-field inflation. While it is universally understood that single-field inflation predicts negligible scale-dependent bias compared to current observational uncertainties, there is still disagreement on the exact level of scale-dependent bias at a level that could strongly impact inferences made from future surveys. In this paper, we clarify this confusion and derive in various ways that there is exactly zero scale-dependent bias in single-field inflation. Much of the current confusion follows from the fact that single-field inflation does predict a mode coupling of matter perturbations at the level of $f_{NL}^{loc} \approx -5/3$, which naively would lead to scale-dependent bias. However, we show explicitly that this mode coupling cancels out when perturbations are evaluated at a fixed physical scale rather than fixed coordinate scale. Furthermore, we show how the absence of scale-dependent bias can be derived easily in any gauge. This result can then be incorporated into a complete description of the observed galaxy clustering, including the previously studied general relativistic terms, which are important at the same level as scale-dependent bias of order $f_{NL}^{loc} \sim 1$. This description will allow us to draw unbiased conclusions about inflation from future galaxy clustering data.
We present an analysis of Spitzer/IRAC primary transit and secondary eclipse lightcurves measured for HD209458b, using Gaussian process models to marginalise over the intrapixel sensitivity variations in the 3.6 micron and 4.5 micron channels and the ramp effect in the 5.8 micron and 8.0 micron channels. The main advantage of this approach is that we can account for a broad range of degeneracies between the planet signal and systematics without actually having to specify a deterministic functional form for the latter. Our results do not confirm a previous claim of water absorption in transmission. Instead, our results are more consistent with a featureless transmission spectrum, possibly due to a cloud deck obscuring molecular absorption bands. For the emission data, our values are not consistent with the thermal inversion in the dayside atmosphere that was originally inferred from these data. Instead, we agree with another re-analysis of these same data, which concluded a non-inverted atmosphere provides a better fit. We find that a solar-abundance clear-atmosphere model without a thermal inversion underpredicts the measured emission in the 4.5 micron channel, which may suggest the atmosphere is depleted in carbon monoxide. An acceptable fit to the emission data can be achieved by assuming that the planet radiates as an isothermal blackbody with a temperature of $1484\pm 18$ K.
The theory of cosmological fluctuations assumes that the pre-inflationary state of the universe was the quantum vacuum of a scalar field(s) coupled to gravity. The observed cosmic microwave background fluctuations are then interpreted as quantum fluctuations. Here we consider alternate interpretations of the classic calculations of scalar and tensor power spectra by replacing the quantum vacuum with a classical statistical distribution, and suggest a way of distinguishing the quantum from the classical alternatives. The possibility that the latter is governed by a fundamental length scale as in string theory is also explored.
We investigate the validity of the equivalence principle near horizons in string theory, analyzing the breakdown of effective field theory caused by longitudinal string spreading effects. An experiment is set up where a detector is thrown into a black hole a long time after an early infalling string. Light cone gauge calculations, taken at face value, indicate a detectable level of root-mean-square longitudinal spreading of the initial string as measured by the late infaller. This results from the large relative boost between the string and detector in the near horizon region, which develops automatically despite their modest initial energies outside the black hole and the weak curvature in the geometry. We subject this scenario to basic consistency checks, using these to obtain a relatively conservative criterion for its detectability. In a companion paper, we exhibit longitudinal nonlocality in well-defined gauge-invariant S-matrix calculations, obtaining results consistent with the predicted spreading albeit not in a direct analogue of the black hole process. We discuss applications of this effect to the firewall paradox, and estimate the time and distance scales it predicts for new physics near black hole and cosmological horizons.
The Square Kilometer Array (SKA) is an international effort to build the world's largest radio telescope, with one square kilometer collecting area. Besides its ambitious scientific objectives, such as probing the cosmic dawn and cradle of life, SKA also demands several revolutionary technological breakthroughs, with ultra-high precision synchronisation of the frequency references for thousands of antennas being one of them. In this report, aimed at applications to SKA, we demonstrate a frequency reference synchronization and dissemination scheme with the phase noise compensation function placed at the client site. Hence, one central hub can be linked to a large number of client sites, forming a star-shaped topology. As a performance test, the 100 MHz reference signal from a Hydrogen maser clock is disseminated and recovered at two remote sites. Phase noise characteristics of the recovered reference frequency signal coincides with that of the hydrogen-maser source and satisfies SKA requirement.
It has been pointed out that the null energy condition can be violated stably in some non-canonical scalar-field theories. This allows us to consider the Galilean Genesis scenario in which the universe starts expanding from Minkowski spacetime and hence is free from the initial singularity. We use this scenario to study the early-time completion of inflation, pushing forward the recent idea of Pirtskhalava et al. We present a generic form of the Lagrangian governing the background and perturbation dynamics in the Genesis phase, the subsequent inflationary phase, and the graceful exit from inflation, as opposed to employing the effective field theory approach. Our Lagrangian belongs to a more general class of scalar-tensor theories than the Horndeski theory and Gleyzes-Langlois-Piazza-Vernizzi generalization, but still has the same number of the propagating degrees of freedom, and thus can avoid Ostrogradski instabilities. We investigate the generation and evolution of primordial perturbations in this scenario and show that one can indeed construct a stable model of inflation preceded by (generalized) Galilean Genesis.
We argue a scenario motivated by the context of string landscape, where our universe is produced by a new vacuum bubble embedded in an old bubble and these bubble universes have not only different cosmological constants, but also their own different gravitational constants. We study these effects on the primordial curvature perturbations. In order to construct a model of varying gravitational constants, we use the Jordan-Brans-Dicke (JBD) theory where different expectation values of scalar fields produce difference of constants. In this system, we investigate the nucleation of bubble universe and dynamics of the wall separating two spacetimes. In particular, the primordial curvature perturbation on superhorizon scales can be affected by the wall trajectory as the boundary effect. We show the effect of gravitational constant in the exterior bubble universe can provide a peak like a bump feature at a large scale in a modulation of power spectrum.
The present Editorial introduces the Special Issue dedicated by the journal Universe to the General Theory of Relativity, the beautiful theory of gravitation of Einstein, a century after its birth. It reviews some of its key features in a historical perspective, and, in welcoming distinguished researchers from all over the world to contribute it, some of the main topics at the forefront of the current research are outlined.
Although f(R) modifications of late time cosmology is successful in explaining present cosmic acceleration, it is very difficult to simultaneously satisfy the fifth-force constraint. Even when the fifth-force constraint is satisfied, the effective scalar degree of freedom may move to a point (close to its minima) in the field space where the Ricci scalar diverges. We elucidate this point further with a specific example of f(R) gravity that incorporates several viable f(R) gravity models in the literature. In particular, we show that the nonlinear evolution of the scalar field in pressureless contracting dust can easily lead to the curvature singularity, making this theory unviable.
We derive the second order correction to the scalar and tensor spectral tilts for the inflationary models with non-minimally derivative coupling. The non-minimally kinetic coupling to Einstein tensor brings the energy scale in the inflationary models down to be sub-Planckian. In the high friction limit, the Lyth bound is modified with an extra suppression factor, so that the field excursion of the inflaton is sub-Planckian. The inflationary models with non-minimally derivative coupling are more consistent with observations.
We present a new approach to find accurate solutions to the Poisson equation, as obtained from the steady-state limit of a diffusion equation with strong source terms. For this purpose, we start from Boltzmann's kinetic theory and investigate the influence of higher order terms on the resulting macroscopic equations. By performing an appropriate expansion of the equilibrium distribution, we provide a method to remove the unnecessary terms up to a desired order and show that it is possible to find, with high level of accuracy, the steady-state solution of the diffusion equation for sizeable Knudsen numbers. In order to test our kinetic approach, we discretise the Boltzmann equation and solve the Poisson equation, spending up to six order of magnitude less computational time for a given precision than standard lattice Boltzmann methods.
CDMS II data from the 5-tower runs at the Soudan Underground Laboratory were reprocessed with an improved charge-pulse fitting algorithm. Two new analysis techniques to reject surface-event backgrounds were applied to the 612 kg days germanium-detector WIMP-search exposure. An extended analysis was also completed by decreasing the 10 keV analysis threshold to $\sim$5 keV, to increase sensitivity near a WIMP mass of 8 GeV/$c^2$. After unblinding, there were zero candidate events above a deposited energy of 10 keV and 6 events in the lower-threshold analysis. This yielded minimum WIMP-nucleon spin-independent scattering cross-section limits of $1.8 \times 10^{-44}$ and $1.18 \times 10 ^{-41}$ cm$^2$ at 90\% confidence for 60 and 8.6 GeV/$c^2$ WIMPs, respectively. This improves the previous CDMS II result by a factor of 2.4 (2.7) for 60 (8.6) GeV/$c^2$ WIMPs.
The low-mass X-ray binary Scorpius X-1 (Sco X-1) is potentially the most luminous source of continuous gravitational-wave radiation for interferometers such as LIGO and Virgo. For low-mass X-ray binaries this radiation would be sustained by active accretion of matter from its binary companion. With the Advanced Detector Era fast approaching, work is underway to develop an array of robust tools for maximizing the science and detection potential of Sco X-1. We describe the plans and progress of a project designed to compare the numerous independent search algorithms currently available. We employ a mock-data challenge in which the search pipelines are tested for their relative proficiencies in parameter estimation, computational efficiency, robust- ness, and most importantly, search sensitivity. The mock-data challenge data contains an ensemble of 50 Scorpius X-1 (Sco X-1) type signals, simulated within a frequency band of 50-1500 Hz. Simulated detector noise was generated assuming the expected best strain sensitivity of Advanced LIGO and Advanced VIRGO ($4 \times 10^{-24}$ Hz$^{-1/2}$). A distribution of signal amplitudes was then chosen so as to allow a useful comparison of search methodologies. A factor of 2 in strain separates the quietest detected signal, at $6.8 \times 10^{-26}$ strain, from the torque-balance limit at a spin frequency of 300 Hz, although this limit could range from $1.2 \times 10^{-25}$ (25 Hz) to $2.2 \times 10^{-26}$ (750 Hz) depending on the unknown frequency of Sco X-1. With future improvements to the search algorithms and using advanced detector data, our expectations for probing below the theoretical torque-balance strain limit are optimistic.
The parallel globe is an old, very simple and ingenious device that, when systematically employed in astronomy classes, becomes a teaching tool with great potential. Properly oriented according to the local meridian, this instrument allows us to follow the shadows in any region of the Earth that is illuminated by the Sun, as well as offering a clear view of the terminator, the fast-moving grey line that divides the day from the night on our planet. With knowledge of the shadows, it is possible to estimate the latitude of a site and to infer local solar time anywhere in the planet's sunlit hemisphere. Furthermore, by using the parallel globe we may understand simply the existence of regions in which objects sometimes do not cast shadows, and also other regions which, on the contrary, sometimes become "long-shadow" countries. In this work, we first review the device and the basics of its assembly and operation. In the second part, we describe in detail some activities targeted to facilitate its use in the classroom, which our research group has been developing during teacher training workshops.
Recently the AMS-02 experiment reported an excess of cosmic ray antiprotons over the expected astrophysical background. We interpret the excess as a signal from annihilating or decaying dark matter and find that the observed spectrum is well fitted by adding contributions from the annihilation or decay of dark matter with mass of O(TeV) or larger. Interestingly, Wino dark matter with mass of around 3 TeV, whose thermal relic abundance is consistent with present dark matter abundance, can explain the antiproton excess. We also discuss the implications for the decaying gravitino dark matter with R-parity violation.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
We study the dynamical evolution of a stellar disk orbiting a massive black hole. We explore the role of two-body relaxation, mass segregation, stellar evolution and binary heating in affecting the disk evolution, and consider the impact of the nuclear cluster structure and the stellar-disk mass-function. We use analytic arguments and numerical calculations, and apply them to study the evolution of a stellar disk (similar to that observed in the Galactic center; GC), both on the short (few Myr) and longer (100 Myr) evolutionary timescales. We find the dominant processes affecting the disk evolution are two-body relaxation and mass segregation where as binary heating have only a little contribution. Massive stars play a dominant role in kinematically heating low mass stars, and driving them to high eccentricities/inclinations. Multi-mass models with realistic mass-functions for the disk stars show the disk structure to be mass stratified, with the most massive stars residing in thinner structures. Stellar evolution plays an important role in decreasing the number of massive stars with time, thereby leading to slower relaxation, where the remnant compact objects of these stars are excited to higher eccentricities/inclinations. At these later evolutionary stages dynamical heating by the nuclear cluster plays a progressively more important role. We conclude that the high eccentricities of the disk-stars in the Galactic Center suggest that the disk formed with initially high eccentricities, or that collective or secular processes dominate the disk evolution. Finally, we find that the disk structure is expected to keep a thin structure even after 100 Myrs. It therefore suggests earlier disks now containing only older, lower mass stars might still be observed in the Galactic center, unless destroyed/smeared by other non-two-body relaxation processes.
From pulsar scintillations we infer the presence of sheet-like structures in the ISM; it has been suggested that these are current sheets. Current sheets probably play an important role in heating the solar corona, and there is evidence for their presence in the solar wind. Such magnetic discontinuities have been found in numerical simulations with particular boundary conditions, as well as in simulations using an incompressible equation of state. Here, I investigate their formation under more general circumstances by means of topological considerations as well as numerical simulations of the relaxation of an arbitrary smoothly-varying magnetic field. The simulations are performed with a variety of parameters and boundary conditions: in low, high and of-order-unity plasma-$\beta$ regimes, with periodic and fixed boundaries, with and without a friction force, at various resolutions and with various diffusivities. Current sheets form, over a dynamical timescale, under {\it all} conditions explored. At higher resolution they are thinner, and there is a greater number of weaker current sheets. The magnetic field eventually relaxes into a smooth minimum energy state, the energy of which depends on the magnetic helicity, as well as on the nature of the boundaries.
The massive cluster MACSJ1149.5+2223 (z=0.544) displays five very large lensed images of a well resolved spiral galaxy at z_spect=1.491. It is within one of these images that the first example of a multiply-lensed supernova has been detected recently in deep Hubble Frontier Field imaging. The depth of this data also reveals many HII regions within the lensed spiral galaxy which we identify between the five counter-images. Here we expand the capability of our free-form method to incorporate these HII regions locally, with other reliable lensed galaxies added for a global solution. This improved accuracy allows us to estimate when the Refsdal supernova will appear within the other lensed images of the spiral galaxy to an accuracy of about 7%. We predict the reappearance of this supernova in one of the counter-images (RA=11:49:36.025, DEC=+22:23:48.11, J2000) on November first 2015 (plus minus 25 days), offering a unique opportunity to study the early phases of this supernova and to examine the consistency of the mass model and the cosmological model that have an impact on the time delay prediction.
Kozai-Lidov (KL) oscillations in hierarchical triple systems have found application to many astrophysical contexts, including planet formation, type Ia supernovae, and supermassive black hole dynamics. The period of these oscillations is known at the order-of-magnitude level, but dependences on the initial mutual inclination or inner eccentricity are not typically included. In this work I calculate the period of KL oscillations ($t_{\textrm{KL}}$) exactly in the test particle limit at quadrupole order (TPQ). I explore the parameter space of all hierarchical triples at TPQ and show that except for triples on the boundary between libration and rotation, the period of KL oscillations does not vary by more than a factor of a few. The exact period may be approximated to better than 2 per cent for triples with mutual inclinations between 60$^{\circ}$ and 120$^{\circ}$ and initial eccentricities less than $\sim$0.3. In addition, I derive an analytic expression for the period of octupole-order oscillations due to the `eccentric KL mechanism' (EKM). I show that the timescale for EKM oscillations is proportional to $\epsilon_{\textrm{oct}}^{-1/2}$, where $\epsilon_{\textrm{oct}}$ measures the strength of octupole perturbations relative to quadrupole perturbations.
We use public data for 105783 quasars from The Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) that include spectral monochromatic luminosities at 5100\AA, 3000\AA, and 1350\AA, and the corresponding observed broad-band ugriz, VRI (converted), JHK and WISE magnitudes, and derive broad-band-to-monochromatic luminosity ratios independent of a cosmological model. The ratios span the redshift range of z=0.1-4.9 and may serve as a proxy for measuring the bolometric luminosity, broad line region (BLR) radii and/or black hole masses, whenever flux-calibrated spectra are unavailable or the existing spectra have low signal-to-noise ratios. They are provided both in tabular and parametrized form.
We have succeeded in obtaining magnetised equilibrium states with differential rotation and differential toroidal magnetic fields. If an internal toroidal field of a proto-neutron star is wound up from the initial poloidal magnetic field by differential rotation, the distribution of the toroidal magnetic field is determined by the profile of this differential rotation. However, the distributions of the toroidal fields in all previous magnetised equilibrium studies do not represent the magnetic winding by the differential rotation of the star. In this paper, we investigate a formulation of a differential toroidal magnetic field that represents the magnetic field wound up by differential rotation. We have developed two functional forms of differential toroidal fields which correspond to a v-constant and a j-constant field in analogy to differential rotations. As the degree of the differential becomes very high, the toroidal magnetic field becomes highly localised and concentrated near the rotational axis. Such a differential toroidal magnetic field would suppress the low-T/|W| instability more efficiently even if the total magnetic field energy is much smaller than that of a non-differential toroidal magnetic field.
The detection of reflected light from an exoplanet is a difficult technical challenge at optical wavelengths. Even though this signal is expected to replicate the stellar signal, not only is it several orders of magnitude fainter, but it is also hidden among the stellar noise. We apply a variant of the cross-correlation technique to HARPS observations of 51 Peg to detect the reflected signal from planet 51 Peg b. Our method makes use of the cross-correlation function of a binary mask with high-resolution spectra to amplify the minute planetary signal that is present in the spectra by a factor proportional to the number of spectral lines when performing the cross correlation. The resulting cross-correlation functions are then normalized by a stellar template to remove the stellar signal. Carefully selected sections of the resulting normalized CCFs are stacked to increase the planetary signal further. The recovered signal allows probing several of the planetary properties, including its real mass and albedo. We detect evidence for the reflected signal from planet 51 Peg b at a significance of 3\sigma_noise. The detection of the signal permits us to infer a real mass of 0.46^+0.06_-0.01 M_Jup (assuming a stellar mass of 1.04\;M_Sun) for the planet and an orbital inclination of 80^+10_-19 degrees. The analysis of the data also allows us to infer a tentative value for the (radius-dependent) geometric albedo of the planet. The results suggest that 51Peg b may be an inflated hot Jupiter with a high albedo (e.g., an albedo of 0.5 yields a radius of 1.9 \pm 0.3 R_Jup for a signal amplitude of 6.0\pm0.4 x 10^-5). We confirm that the method we perfected can be used to retrieve an exoplanet's reflected signal, even with current observing facilities. The advent of next generation of observing facilities will yield new opportunities for this type of technique to probe deeper into exoplanets.
In the course of the TOPoS (Turn Off Primordial Stars) survey, aimed at discovering the lowest metallicity stars, we have found several carbon-enhanced metal-poor (CEMP) stars. We here present our analysis of six CEMP stars. Calcium and carbon are the only elements that can be measured in all six stars. The range is -5.0<=[Ca/H]< -2.1 and 7.12<=A(C)<=8.65. For star SDSS J1742+2531 we were able to detect three FeI lines from which we deduced [Fe/H]=-4.80, from four CaII lines we derived [Ca/H]=-4.56, and from synthesis of the G-band we derived A(C)=7.26. For SDSS J1035+0641 we were not able to detect any iron lines, yet we could place a robust (3sigma) upper limit of [Fe/H]< -5.0 and measure the Ca abundance, with [Ca/H]=-5.0, and carbon, A(C)=6.90. No lithium is detected in the spectrum of SDSS J1742+2531 or SDSS J1035+0641, which implies a robust upper limit of A(Li)<1.8 for both stars. Our measured carbon abundances confirm the bimodal distribution of carbon in CEMP stars, identifying a high-carbon band and a low-carbon band. We propose an interpretation of this bimodality according to which the stars on the high-carbon band are the result of mass transfer from an AGB companion, while the stars on the low-carbon band are genuine fossil records of a gas cloud that has also been enriched by a faint supernova (SN) providing carbon and the lighter elements. (Abridged)
We investigate the variability of exoplanetary radio emission using stellar magnetic maps and 3D field extrapolation techniques. We use a sample of hot Jupiter hosting stars, focusing on the HD 179949, HD 189733 and tau Boo systems. Our results indicate two time-scales over which radio emission variability may occur at magnetised hot Jupiters. The first is the synodic period of the star-planet system. The origin of variability on this time-scale is the relative motion between the planet and the interplanetary plasma that is co-rotating with the host star. The second time-scale is the length of the magnetic cycle. Variability on this time-scale is caused by evolution of the stellar field. At these systems, the magnitude of planetary radio emission is anticorrelated with the angular separation between the subplanetary point and the nearest magnetic pole. For the special case of tau Boo b, whose orbital period is tidally locked to the rotation period of its host star, variability only occurs on the time-scale of the magnetic cycle. The lack of radio variability on the synodic period at tau Boo b is not predicted by previous radio emission models, which do not account for the co-rotation of the interplanetary plasma at small distances from the star.
We present both the technical overview and main science drivers of the fourth
phase of the Optical Gravitational Lensing Experiment (hereafter OGLE-IV).
OGLE-IV is currently one of the largest sky variability surveys worldwide,
targeting the densest stellar regions of the sky. The survey covers over 3000
square degrees in the sky and monitors regularly over a billion sources.
The main targets include the inner Galactic Bulge and the Magellanic System.
Their photometry spans the range of $12<I<21$ mag and $13<I<21.7$ mag,
respectively. Supplementary shallower Galaxy Variability Survey covers the
extended Galactic bulge and 2/3 of the whole Galactic disk within the magnitude
range of $10<I<19$ mag. All OGLE-IV surveys provide photometry with
milli-magnitude accuracy at the bright end. The cadence of observations varies
from 19-60 minutes in the inner Galactic bulge to 1-3 days in the remaining
Galactic bulge fields, Magellanic System and the Galactic disk.
OGLE-IV provides the astronomical community with a number of real time
services. The Early Warning System (EWS) contains information on two thousand
gravitational microlensing events being discovered in real time annually, the
OGLE Transient Detection System (OTDS) delivers over 200 supernovae a year. We
also provide the real time photometry of unpredictable variables such as
optical counterparts to the X-ray sources and R CrB stars.
Hundreds of thousands new variable stars have already been discovered and
classified by the OGLE survey. The number of new detections will be at least
doubled during the current OGLE-IV phase. The survey was designed and optimized
primarily to conduct the second generation microlensing survey for exoplanets.
It has already contributed significantly to the increase of the discovery rate
of microlensing exoplanets and free-floating planets.
We obtain solutions to the coupled Schr\"odinger-Poisson equations. The solutions describe the evolution of cold dark matter density perturbations in an otherwise homogeneous expanding Friedmann universe. We discuss the relationships between descriptions of cold dark matter in terms of a pressureless fluid, in terms of a wavefunction, of a classical scalar field, and a quantum scalar field. We identify the regimes where the various descriptions coincide and where they differ.
We have used the Submillimeter Array to image, at ~1.5" resolution, C2H (3-2) emission from the circumstellar disk orbiting the nearby (D = 54 pc), ~8 Myr-old, ~0.8 Msun classical T Tauri star TW Hya. The SMA imaging reveals that the C2H emission exhibits a ring-like morphology. Based on a model in which the C2H column density follows a truncated radial power-law distribution, we find that the inner edge of the ring lies at ~45 AU, and that the ring extends to at least ~120 AU. Comparison with previous (single-dish) observations of C2H (4-3) emission indicates that the C2H molecules are subthermally excited and, hence, that the emission arises from the relatively warm, tenuous upper atmosphere of the disk. We propose that the C2H emission most likely traces particularly efficient photo-destruction of small grains and/or photodesorption and photodissociation of hydrocarbons derived from grain ice mantles in the surface layers of the outer disk. The presence of a C2H ring in the TW Hya disk hence likely serves as a marker of dust grain processing and radial and vertical grain size segregation within the disk.
Hydrocarbon organic material, as found in the interstellar medium, exists in complex mixtures of aromatic and aliphatic forms. It is considered to be originated from carbon enriched giant stars during their final stages of evolution, when very strong mass loss occurs in a few thousand years on their way to become planetary nebulae. We show here that the same organic compounds appear to be formed in previous stages of the evolution of giant stars. More specifically, during the first ascending giant branch K-type stars. According to our model this happens only when these stars are being abruptly enriched with lithium together with the formation of a circumstellar shell with a strong mass loss during just a few thousand years. This sudden mass loss is, on an average, a thousand times larger than that of normal Li-poor K giant stars. This shell would later be detached, specially when the star stops its Li enrichment and a rapid photospheric Li depletion occurs. In order to gain extra carbon-based material to form the organic hydrocarbons, and also to explain the presence of complex inorganic compounds in these stars, we propose an interaction of these strong winds with remaining asteroidal/cometary disks that already existed around these stars since they were dwarf A-type stars. The mechanism of interaction presented here is successful to explain the presence of inorganic compounds, however it is unable to produce new carbon free atoms to form the organic hydrocarbon compounds. Finally, we discuss some suggestions and speculations that can eventually help solving the long-standing puzzle of Li-rich giants.
Attention has focussed recently on models of inflation that involve a second or more fields with a mass near the inflationary Hubble parameter $H$, as may occur in supersymmetric theories if the supersymmetry-breaking scale is not far from $H$. Quasi-single-field (QSF) inflation is a relatively simple family of phenomenological models that serve as a proxy for theories with additional fields with masses $m\sim H$. Here we consider the tensor-scalar-scalar (tss) three-point function that arises in QSF inflation. Since QSF involves fields in addition to the inflaton, the consistency conditions between correlations that arise in single-clock inflation are not necessarily satisfied. As a result, the tensor-scalar-scalar correlation may be larger than in single-field inflation. This tss correlation gives rise to local departures from statistical isotropy, or in other words, a nontrivial four-point function. The presence of the tensor mode may moreover be inferred geometrically from the shape dependence of the four-point function. We estimate the size of a galaxy survey required to detect this tss correlation in QSF inflation as a function of $H$. Our study of primordial correlators which include gravitons in seeking imprints of additional fields with masses $m\sim H$ during inflation can be seen as complementary to the recent "cosmological collider physics" proposal.
Much of the dynamical structure of the Kuiper belt can be explained if Neptune migrated over several AU, and/or if Neptune was scattered to an eccentric orbit during planetary instability. An outstanding problem with the existing formation models is that the distribution of orbital inclinations they predicted is narrower than the one inferred from observations. Here we perform numerical simulations of Kuiper belt formation starting from an initial state with Neptune at 20 < a_{N,0} < 30 AU and a dynamically cold outer disk extending from beyond a_{N,0} to 30 AU. Neptune's orbit is migrated into the disk on an e-folding timescale 1 <= tau <= 100 Myr. A small fraction (~10^{-3}) of the disk planetesimals become implanted into the Kuiper belt in the simulations. By analyzing the orbital distribution of the implanted bodies in different cases we find that the inclination constraint implies that tau >= 10 Myr and a_{N,0} <= 25 AU. The models with tau < 10 Myr do not satisfy the inclination constraint, because there is not enough time for various dynamical processes to raise inclinations. The slow migration of Neptune is consistent with other Kuiper belt constraints, and with recently developed models of planetary instability/migration. Neptune's eccentricity and inclination are never large in these models (e_N<0.1, i_N<2 deg), as required to avoid excessive orbital excitation in the >40 AU region, where the Cold Classicals presumably formed.
In this study we describe an innovative method to determine potential sites for optical and infrared astronomical observations in the Andes region of northern South America. The method computes the Clear sky fraction (CSF) from Geostationary Observational Environmental Satellite (GOES) data for the years 2008-12 through a comparison with temperatures obtained from long-term records of weather stations and atmospheric temperature profiles from radiosonde. Criteria for sky clearance were established for two infrared GOES channels in order to determine potential sites in the Andes region of northern South-America. The method was validated using the reported observed hours at Observatorio Nacional de Llano del Hato in Venezuela. Separate CSF percentages were computed for dry and rainy seasons for both, photometric and spectroscopic night qualities. Twelve sites with five year averages of CSF for spectroscopic nights larger than 30% during the dry seasons were found to be suitable for astronomical observations. The best site with (220$\pm$42) spectroscopic clear nights per year is located in the Andes of Venezuela (70$^{\circ}$28'48"W, 9$^{\circ}$5'60"N) at an altitude of 3480 meters. Lower quality regions were found in Sierra Nevada de Santamarta and Serran\'ia del Perij\'a with (126$\pm$34) and (111$\pm$27) clear nights per year, respectively. Sites over the Andes are identified in Norte de Santander with (107$\pm$23) and in the north-east part of Boyac\'a with a mean of (94$\pm$13) clear nights per year. Two sites at low latitude located in Ecuador with more than 100 clear nights per year and with similar seasonal CSF percentages were also identified. Five year evolution suggest a possible correlation be tween the lowest percentages observed during the rainy seasons of 2010 and 2011 with positive values of the Southern Oscillation Index.
We perform smoothed particle hydrodynamics (SPH) simulations for merging binary carbon-oxygen (CO) white dwarfs (WDs) with masses of $1.1$ and $1.0$ $M_\odot$, until the merger remnant reaches a dynamically steady state. Using these results, we assess whether the binary could induce a thermonuclear explosion, and whether the explosion could be observed as a type Ia supernova (SN Ia). We investigate three explosion mechanisms: a helium-ignition following the dynamical merger (`helium-ignited violent merger model'), a carbon-ignition (`carbon-ignited violent merger model'), and an explosion following the formation of the Chandrasekhar mass WD (`Chandrasekhar mass model'). An explosion of the helium-ignited violent merger model is possible, while we predict that the resulting SN ejecta are highly asymmetric since its companion star is fully intact at the time of the explosion. The carbon-ignited violent merger model can also lead to an explosion. However, the envelope of the exploding WD spreads out to $\sim 0.1R_\odot$; it is much larger than that inferred for SN 2011fe ($< 0.1R_\odot $) while much smaller than that for SN 2014J ($\sim 1R_\odot$). For the particular combination of the WD masses studied in this work, the Chandrasekhar mass model is not successful to lead to an SN Ia explosion. Besides these assessments, we investigate the evolution of unbound materials ejected through the merging process (`merger ejecta'), assuming a case where the SN Ia explosion is not triggered by the helium- or carbon-ignition during the merger. The merger ejecta interact with the surrounding interstellar medium, and form a shell. The shell has a bolometric luminosity of more than $2 \times 10^{35}$ ergs$^{-1}$ lasting for $\sim 2 \times 10^4$ yr. If this is the case, Milky Way should harbor about $10$ such shells at any given time.
We present 850 micron observations of the 2-3 Myr cluster IC 348 in the Perseus molecular cloud using the SCUBA-2 camera on the James Clerk Maxwell Telescope. Our SCUBA-2 map has a diameter of 30 arcmin and contains ~370 cluster members, including ~200 objects with IR excesses. We detect a total of 13 discs. Assuming standard dust properties and a gas to dust mass ratio of 100, we derive disc masses ranging from 1.5 to 16 M_JUP . We also detect 8 Class 0/I protostars. We find that the most massive discs (M_Disc > 3 MJUP ; 850 micron fux > 10 mJy) in IC 348 tend to be transition objects according to the characteristic "dip" in their infrared Spectral Energy Distributions (SEDs). This trend is also seen in other regions. We speculate that this could be an initial conditions effect (e.g., more massive discs tend to form giant planets that result in transition disc SEDs) and/or a disc evolution effect (the formation of one or more massive planets results in both a transition disc SED and a reduction of the accretion rate, increasing the lifetime of the outer disc). A stacking analysis of the discs that remain undetected in our SCUBA-2 observations suggests that their median 850 micron flux should be ~1 mJy, corresponding to a disc mass ~0.3 M_JUP (gas plus dust) or ~1 M_Earth of dust. While the available data are not deep enough to allow a meaningful comparison of the disc luminosity functions between IC 348 and other young stellar clusters, our results imply that disc masses exceeding the Minimum Mass Solar Nebula are very rare (~1%) at the age of IC 348, specially around very low-mass stars.
Using the "Updated Nearby Galaxy Catalog", we consider different properties of companion galaxies around luminous hosts in the Local Volume. The data on stellar masses, linear diameters,surface brightnesses, HI-richness, specific star formation rate (sSFR), and morphological types are discussed for members of the nearest groups, including the Milky Way and M 31 groups, as a function of their separation from the hosts. Companion galaxies in groups tend to have lower stellar masses, smaller linear diameters and fainter mean surface brightnesses as the distance to their host decreases. The hydrogen-to-stellar mass ratio of the companions increases with their linear projected separation from the dominant luminous galaxy. This tendency is more expressed around the bulge-dominated hosts. While linear separation of the companions decreases, their mean sSFR becomes lower, accompanied with the increasing sSFR scatter. Typical linear projected separation of dSphs around the bulge-dominated hosts, 350 kpc, is substantially larger than that around the disk-dominated ones, 130 kpc. This difference probably indicates the presence of larger hot/warm gas haloes around the early-type host galaxies. The mean fraction of dSph (quenched) companions in 11 the nearest groups as a function of their projected separation R_p can be expressed as f(E) = 0.55 - 0.69 x R_p. The fraction of dSphs around the Milky Way and M 31 looks much higher than in other nearby groups because the quenching efficiency dramatically increases towards the ultra-low mass companions. We emphasize that the observed properties of the Local Group are not typical for other groups in the Local Volume due to the role of selection effects caused by our location inside the Local Group.
Filament eruptions often lead to coronal mass ejections (CMEs), which can affect critical technological systems in space and on the ground when they interact with the geo-magnetosphere in high speeds. Therefore, it is an important issue to investigate the acceleration mechanisms of CMEs in solar/space physics. Based on observations and simulations, the resistive magnetic reconnection and the ideal instability of magnetic flux rope have been proposed to accelerate CMEs. However, it remains elusive whether both of them play a comparable role during a particular eruption. It has been extremely difficult to separate their contributions as they often work in a close time sequence during one fast acceleration phase. Here we report an intriguing filament eruption event, which shows two apparently separated fast acceleration phases and provides us an excellent opportunity to address the issue. Through analyzing the correlations between velocity (acceleration) and soft (hard) X-ray profiles, we suggest that the instability and magnetic reconnection make a major contribution during the first and second fast acceleration phases, respectively. Further, we find that both processes have a comparable contribution to accelerate the filament in this event.
Starless molecular cores are natural laboratories for interstellar molecular chemistry research. The chemistry of ices in such objects was investigated with a three-phase (gas, surface, and mantle) model. We considered the center part of five starless cores, with their physical conditions derived from observations. The ice chemistry of oxygen, nitrogen, sulfur, and complex organic molecules (COMs) was analyzed. We found that an ice-depth dimension, measured, e.g., in monolayers, is essential for modeling of chemistry in interstellar ices. Particularly, the H2O:CO:CO2:N2:NH3 ice abundance ratio regulates the production and destruction of minor species. It is suggested that photodesorption during core collapse period is responsible for high abundance of interstellar H2O2 and O2H, and other species synthesized on the surface. The calculated abundances of COMs in ice were compared to observed gas-phase values. Smaller activation barriers for CO and H2CO hydrogenation may help explain the production of a number of COMs. The observed abundance of methyl formate HCOOCH3 could be reproduced with a 1kyr, 20K temperature spike. Possible desorption mechanisms, relevant for COMs, are gas turbulence (ice exposure to interstellar photons) or a weak shock within the cloud core (grain collisions). To reproduce the observed COM abundances with the present 0D model, 1-10% of ice mass needs to be sublimated. We estimate that the lifetime for starless cores likely does not exceed 1Myr. Taurus cores are likely to be younger than their counterparts in most other clouds.
We apply the statefinder hierarchy plus the fractional growth parameter to explore the extended Ricci dark energy (ERDE) model, in which there are two independent coefficients $\alpha$ and $\beta$. By adjusting them, we plot evolution trajectories of some typical parameters, including Hubble expansion rate $E$, deceleration parameter $q$, the third and fourth order hierarchy $S_3^{(1)}$ and $S_4^{(1)}$ and fractional growth parameter $\epsilon$, respectively, as well as several combinations of them. For the case of variable $\alpha$ and constant $\beta$, in the low-redshift region the evolution trajectories of $E$ are in high degeneracy and that of $q$ separate somewhat. However, the $\Lambda$CDM model is confounded with ERDE in both of these two cases. $S_3^{(1)}$ and $S_4^{(1)}$, especially the former, perform much better. They can differentiate well only varieties of cases within ERDE except $\Lambda$CDM in the low-redshift region. For high-redshift region, combinations $\{S_n^{(1)},\epsilon\}$ can break the degeneracy. Both of $\{S_3^{(1)},\epsilon\}$ and $\{S_4^{(1)},\epsilon\}$ have the ability to discriminate ERDE with $\alpha=1$ from $\Lambda$CDM, of which the degeneracy cannot be broken by all the before-mentioned parameters. For the case of variable $\beta$ and constant $\alpha$, $S_3^{(1)}(z)$ and $S_4^{(1)}(z)$ can only discriminate ERDE from $\Lambda$CDM. Nothing but pairs $\{S_3^{(1)},\epsilon\}$ and $\{S_4^{(1)},\epsilon\}$ can discriminate not only within ERDE but also ERDE from $\Lambda$CDM. Finally we find that $S_3^{(1)}$ is surprisingly a better choice to discriminate within ERDE itself, and ERDE from $\Lambda$CDM as well, rather than $S_4^{(1)}$.
The geometry of the cosmic web drives in part the spin acquisition of galaxies. This can be explained in a Lagrangian framework, by identifying the specific long-wavelength correlations within the primordial Gaussian random field which are relevant to spin acquisition. Tidal Torque Theory is revisited in the context of such anisotropic environments, biased by the presence of a filament within a wall. The point process of filament-type saddles represents it most efficiently. The constrained misalignment between the tidal and the inertia tensors in the vicinity of filament-type saddles simply explains the distribution of spin directions. This misalignment implies in particular an azimuthal orientation for the spins of more massive galaxies and a spin alignment with the filament for less massive galaxies. This prediction is found to be in qualitative agreement with measurements in Gaussian random fields and N-body simulations. It relates the transition mass to the geometry of the saddle, and accordingly predicts its measured scaling with the mass of non-linearity. Implications for galaxy formation and weak lensing are briefly discussed, as is the dual theory of spin alignments in walls.
The ambitious science goals of the Large Synoptic Survey Telescope (LSST) have motivated a search for new and unexpected sources of systematic error in the LSST camera. Flat-field images are a rich source of data on sensor anomalies, although such effects are typically dwarfed by shot noise in a single flat field. After combining many ($\sim 500$) such images into `ultraflats' to reduce the impact of shot noise, we perform photon transfer analysis on a pixel-by-pixel basis and observe no spatial structure in pixel linearity or gain at light levels of 100 ke$^-$ and below. At 125 ke$^-$, a columnar structure is observed in the gain map--we attribute this to a flux-dependent charge transfer inefficiency. We also probe small-scale variations in effective pixel size by analyzing pixel-neighbor correlations in ultraflat images, where we observe clear evidence of intrinsic variation in effective pixel size in an LSST prototype sensor near the $\sim .3\%$ level.
Context. Most stars are born in clusters, thus the protoplanetary discs surrounding the newly formed stars might be influenced by this environment. Isolated star-disc encounters have previously been studied, and it was shown that very close encounters are necessary to completely destroy discs. However, relatively distant encounters are still able to change the disc size considerably. Aims. We quantify the importance of disc-size reduction that is due to stellar encounters in an entire stellar population. Methods. We modelled young, massive clusters of different densities using the code Nbody6 to determine the statistics of stellar encounter parameters. In a second step, we used these parameters to investigate the effect of the environments on the disc size. For this purpose, we performed a numerical experiment with an artificial initial disc size of 105 AU. Results. We quantify to which degree the disc size is more sensitive to the cluster environment than to the disc mass or frequency. We show that in all investigated clusters a large portion of discs is significantly reduced in size. After 5 Myr, the fraction of discs smaller than 1000 AU in ONC-like clusters with an average number density of 60pc$^3$, the fraction of discs smaller than 1000 AU is 65%, while discs smaller than 100 AU make up 15%. These fractions increase to 84% and 39% for discs in denser clusters like IC 348 (500pc$^3$). Even in clusters with a density four times lower than in the ONC (15pc$^3$), about 43% of all discs are reduced to sizes below 1 000 AU and roughly 9% to sizes below 100 AU. Conclusions. For any disc in the ONC that initially was larger than 1 000 AU, the probability to be truncated to smaller disc sizes as a result of stellar encounters is quite high. Thus, among other effects, encounters are important in shaping discs and potentially forming planetary systems in stellar clusters.
We use Spitzer observations of the rich population of Asymptotic Giant Branch stars in the Large Magellanic Cloud (LMC) to test models describing the internal structure and nucleosynthesis of the most massive of these stars, i.e. those with initial mass above $\sim 4M_{\odot}$. To this aim, we compare Spitzer observations of LMC stars with the theoretical tracks of Asymptotic Giant Branch models, calculated with two of the most popular evolution codes, that are known to differ in particular for the treatment of convection. Although the physical evolution of the two models are significantly different, the properties of dust formed in their winds are surprisingly similar, as is their position in the colour-colour (CCD) and colour-magnitude (CMD) diagrams obtained with the Spitzer bands. This model independent result allows us to select a well defined region in the ($[3.6]-[4.5], [5.8]-[8.0]$) plane, populated by AGB stars experiencing Hot Bottom Burning, the progeny of stars with mass $M\sim 5.5M_{\odot}$. This result opens up an important test of the strength hot bottom burning using detailed near-IR (H and K bands) spectroscopic analysis of the oxygen-rich, high luminosity candidates found in the well defined region of the colour-colour plane. This test is possible because the two stellar evolution codes we use predict very different results for the surface chemistry, and the C/O ratio in particular, owing to their treatment of convection in the envelope and of convective boundaries during third dredge-up. The differences in surface chemistry are most apparent when the model stars reach the phase with the largest infrared emission.
CONTEXT: The Virgo direction has been observed at many wavelengths in the recent years, in particular in the ultraviolet with GALEX. The far ultraviolet (FUV) diffuse light detected by GALEX bears interesting information on the large scale distribution of Galactic dust, owing to the GALEX FUV band sensitivity and resolution. AIMS: We aim to characterise the ultraviolet large scale distribution of diffuse emission in the Virgo direction. A map of this emission may become useful for various studies by identifying regions where dust affects observations by either scattering light or absorbing radiation. METHODS: We construct mosaics of the FUV and near ultraviolet diffuse emission over a large sky region (RA 12 to 13 hours, DEC 0 to 20 degrees) surrounding the Virgo cluster, using all the GALEX available data in the area. We test for the first time the utilisation of the FUV diffuse light as a Galactic extinction E(B-V) tracer. RESULTS: The FUV diffuse light scattered on cirrus reveals details in their geometry. Despite a large dispersion, the FUV diffuse light correlates roughly with other Galactic dust tracers (coming from IRAS, Herschel, Planck), offering an opportunity to use the FUV emission to locate them in future studies with a better resolution (about 5 arcsec native resolution, 20 arcsec pixels maps presented in this paper) than several usual tracers. Estimating the Galactic dust extinction on the basis of this emission allows us to find a smaller dispersion in the NUV-i colour of background galaxies at a given E(B-V)than with other tracers. The diffuse light mosaics obtained in this work are made publicly available.
The high-frequency-peaked BL Lac (HBL) 1ES 0806+524 (z = 0.138) was discovered in VHE $\gamma$ rays in 2008. Until now, the broad-band spectrum of 1ES 0806+524 has been only poorly characterized, in particular at high energies. We analysed multiwavelength observations from $\gamma$ rays to radio performed from 2011 January to March, which were triggered by the high activity detected at optical frequencies. These observations constitute the most precise determination of the broad-band emission of 1ES 0806+524 to date. The stereoscopic MAGIC observations yielded a $\gamma$-ray signal above 250 GeV of $(3.7 \pm 0.7)$ per cent of the Crab Nebula flux with a statistical significance of 9.9 $\sigma$. The multiwavelength observations showed significant variability in essentially all energy bands, including a VHE $\gamma$-ray flare that lasted less than one night, which provided unprecedented evidence for short-term variability in 1ES 0806+524. The spectrum of this flare is well described by a power law with a photon index of $2.97 \pm 0.29$ between $\sim$150 GeV and 1 TeV and an integral flux of $(9.3 \pm 1.9)$ per cent of the Crab Nebula flux above 250 GeV. The spectrum during the non-flaring VHE activity is compatible with the only available VHE observation performed in 2008 with VERITAS when the source was in a low optical state. The broad-band spectral energy distribution can be described with a one-zone Synchrotron Self Compton model with parameters typical for HBLs, indicating that 1ES 0806+524 is not substantially different from the HBLs previously detected.
We report the presence of high significance diffuse radio emission from the Triangulum Australis cluster using observations made with the KAT-7 telescope and propose that this emission is a giant radio halo. We compare the radio power from this proposed halo with X-ray and SZ measurements and demonstrate that it is consistent with the established scaling relations for cluster haloes. By combining the X-ray and SZ data we calculate the ratio of non-thermal to thermal electron pressure within Triangulum Australis to be $X=0.658\pm0.054$. We use this ratio to constrain the maximum magnetic field strength within the halo region to be $B_{\rm max, halo} = 33.08\,\mu$G and compare this with the minimum field strength from equipartition of $B_{\rm min, halo} = 0.77(1+k)^{2/7}\,\mu$G to place limits on the range of allowed magnetic field strength within this cluster. We compare these values to those for more well-studied systems and discuss these results in the context of equipartition of non-thermal energy densities within clusters of galaxies.
We present precision radial velocities and stellar population parameters for 77 star clusters in the Local Group galaxy M33. Our GTC and WHT observations sample both young, massive clusters and known/candidate globular clusters, spanning ages ~ 10^6 - 10^10 yr, and metallicities, [M/H] ~-1.7 to solar. The cluster system exhibits an age-metallicity relation; the youngest clusters are the most metal-rich. When compared to HI data, clusters with [M/H] ~ -1.0 and younger than ~ 4 Gyr are clearly identified as a disc population. The clusters show evidence for strong time evolution in the disc radial metallicity gradient (d[M/H]dt / dR = 0.03 dex/kpc/Gyr). The oldest clusters have stronger, more negative gradients than the youngest clusters in M33. The clusters also show a clear age-velocity dispersion relation. The line of sight velocity dispersions of the clusters increases with age similar to Milky Way open clusters and stars. The general shape of the relation is reproduced by disc heating simulations, and the similarity between the relations in M33 and the Milky Way suggests that heating by substructure, and cooling of the ISM both play a role in shaping this relation. We identify 12 "classical" GCs, six of which are newly identified GC candidates. The GCs are more metal-rich than Milky Way halo clusters, and show weak rotation. The inner (R < 4.5 kpc) GCs exhibit a steep radial metallicity gradient (d[M/H]/dR = -0.29+-0.11 dex/kpc) and an exponential-like surface density profile. We argue that these inner GCs are thick disc rather than halo objects.
Approximate Bayesian Computation (ABC) enables parameter inference for complex physical systems in cases where the true likelihood function is unknown, unavailable, or computationally too expensive. It relies on the forward simulation of mock data and comparison between observed and synthetic catalogues. Here we present cosmoabc, a Python ABC sampler featuring a Population Monte Carlo (PMC) variation of the original ABC algorithm, which uses an adaptive importance sampling scheme. The code is very flexible and can be easily coupled to an external simulator, while allowing to incorporate arbitrary distance and prior functions. As an example of practical application, we coupled cosmoabc with the numcosmo library and demonstrate how it can be used to estimate posterior probability distributions over cosmological parameters based on measurements of galaxy clusters number counts without computing the likelihood function. cosmoabc is published under the GPLv3 license on PyPI and GitHub and documentation is available at this http URL
We describe Space Warps, a novel gravitational lens discovery service that yields samples of high purity and completeness through crowd-sourced visual inspection. Carefully produced colour composite images are displayed to volunteers via a classi- fication interface which records their estimates of the positions of candidate lensed features. Simulated lenses, and expert-classified images which lack lenses, are inserted into the image stream at random intervals; this training set is used to give the vol- unteers feedback on their performance, as well as to calibrate it in order to allow dynamical updates to the probability of any image they classify to contain a lens. Low probability systems are retired from the site periodically, concentrating the sample towards a set of candidates. Having divided 160 square degrees of Canada-France- Hawaii Telescope Legacy Survey (CFHTLS) imaging into some 430,000 overlapping 84 by 84 arcsecond tiles and displaying them on the site, we were joined by around 37,000 volunteers who contributed 11 million image classifications over the course of 8 months. The sample was reduced to 3368 Stage I candidates; these were then refined to yield a sample that we expect to be over 90% complete and 30% pure. We comment on the scalability of the Space Warps system to the wide field survey era, based on our finding that searches of 10$^5$ images can be performed by a crowd of 10$^5$ volunteers in 6 days.
For the past few years, we have observed the central half parsec of our Galaxy in the mid-infrared from 2.8 to 5.1 micron. Our aim is to improve our understanding of the direct environment of SgrA*, the supermassive blackhole at the centre of the Milky Way. This work is described in the present paper and by Moultaka et al. 2015 (submitted). Here, we focus on the study of the spatial distribution of the 12CO ice and gas-phase absorptions. We observed the central half parsec with ISAAC spectrograph located at the UT3/VLT ESO telescope in Chile. The slit was placed along 22 positions arranged parallel to each other to map the region. We built the first data cube in this wavelength range covering the central half parsec. The wavelength interval of the used M-band filter ranges from 4.6 to 5.1 micron. It hosts the P- and R- branches of the ro-vibrational transitions of the gaseous 12CO and 13CO, as well as the absorption band attributed to the 12CO ice at 4.675 micron. Using two calibrators, we could disentangle the local from the line-of-sight absorptions and provide a first-order estimate of the foreground extinction. We find residual ices and gase-phase CO that can be attributed to local absorptions due to material from the interstellar and/or the circumstellar medium of the central parsec. Our finding implies temperatures of the order of 10 to 60K which is in agreement with the presence of water ices in the region highlighted by Moultaka et al. (2004, 2005).
To understand the origin of Solar Energetic Particles (SEPs), we must study their injection time relative to other solar eruption manifestations. Traditionally the injection time is determined using the Velocity Dispersion Analysis (VDA) where a linear fit of the observed event onset times at 1 AU to the inverse velocities of SEPs is used to derive the injection time and path length of the first-arriving particles. VDA does not, however, take into account that the particles that produce a statistically observable onset at 1 AU have scattered in the interplanetary space. We use Monte Carlo test particle simulations of energetic protons to study the effect of particle scattering on the observable SEP event onset above pre-event background, and consequently on VDA results. We find that the VDA results are sensitive to the properties of the pre-event and event particle spectra as well as SEP injection and scattering parameters. In particular, a VDA-obtained path length that is close to the nominal Parker spiral length does not imply that the VDA injection time is correct. We study the delay to the observed onset caused by scattering of the particles and derive a simple estimate for the delay time by using the rate of intensity increase at the SEP onset as a parameter. We apply the correction to a magnetically well-connected SEP event of June 10 2000, and show it to improve both the path length and injection time estimates, while also increasing the error limits to better reflect the inherent uncertainties of VDA.
We present the results from a survey, designed to investigate the accretion process of massive young stellar objects (MYSOs) through near infrared narrow band imaging using the H$_2$ $\nu$=1-0 S(1) transition filter. A sample of 353 Massive Young Stellar Object (MYSO) candidates was selected from the Red MSX Source survey using photometric criteria at longer wavelengths (infrared and submillimeter) and chosen with positions throughout the Galactic Plane. Our survey was carried out at the SOAR Telescope in Chile and CFHT in Hawaii covering both hemispheres. The data reveal that extended H$_2$ emission is a good tracer of outflow activity, which is a signpost of accretion process on young massive stars. Almost half of the sample exhibit extended H$_2$ emission and 74 sources (21\%) have polar morphology, suggesting collimated outflows. The polar-like structures are more likely to appear on radio-quiet sources, indicating these structures occur during the pre-UCHII phase. We also found an important fraction of sources associated with fluorescent H$_2$ diffuse emission that could be due to a more evolved phase. The images also indicate only $\sim$23\% (80) of the sample is associated with extant (young) stellar clusters. These results support the scenario in which massive stars are formed by accretion disks, since the merging of low mass stars would not produce outflow structures.
The high velocity dispersion compact cloud CO-0.30-0.07 is a peculiar molecular clump discovered in the central moleculr zone of the Milky Way, which is characterized by its extremely broad velocity emissions ($\sim 145\ \rm{km s^{-1}}$) despite the absence of internal energy sources. We present new interferometric maps of the cloud in multiple molecular lines in frequency ranges of 265--269 GHz and 276--280 GHz obtained using the Sumbmillimeter Array, along with the single-dish images previously obtained with the ASTE 10-m telescope. The data show that the characteristic broad velocity emissions are predominantly confined in two parallel ridges running through the cloud center. The central ridges are tightly anti-correlated with each other in both space and velocity, thereby sharply dividing the entire cloud into two distinct velocity components (+15 km s$^{-1}$ and +55 km s$^{-1}$). This morphology is consistent with a model in which the two velocity components collide with a relative velocity of 40 $\mathrm{km s^{-1}}$ at the interface defined by the central ridges, although an alternative explanation with a highly inclined expanding-ring model is yet to be fully invalidated. We have also unexpectedly detected several compact clumps ($\lesssim 0.1\ $pc in radius) likely formed by shock compression. The clumps have several features in common with typical star-forming clouds: high densities ($10^{6.5-7.5}\ \mathrm{cm^{-3}}$), rich abundances of hot-core-type molecular species, and relatively narrow velocity widths apparently decoupled from the furious turbulence dominating the cloud. The cloud CO-0.30-0.07 is possibly at an early phase of star formation activity triggered by the shock impact.
RX J1140.1+0307 (hereafter RX1140) is a Narrow Line Seyfert 1 (NLS1) with one of the lowest black hole masses known in an AGN (M ${\le} 10^6$ M$_{\odot}$). We show results from two new {\it XMM-Newton} observations, showing soft 2-10~keV spectra, a strong excess at lower energies, and fast X-ray variability as is typical of this class. The soft excess can be equally well fit by either low temperature Comptonisation or highly smeared, ionised reflection models, but we use a covariance analysis of the fast X-ray variability as well as lag and coherence spectra to show that the low temperature Comptonisation model gives a better description of the break in variability properties between soft and hard X-rays. Both models also require an additional component at the softest energies, as expected from the accretion disc. However, this inner disc spectrum does not join smoothly onto the variable optical and far UV emission (which should be produced in the outer disc) unless the mass is underestimated by an order of magnitude. The variable optical and far UV emission instead suggests that $L/L_{Edd}\sim 10$ through the outer disc, in which case advection and/or wind losses are required to explain the observed broadband spectral energy distribution. However, the similarity of the X-ray properties of RX1140 to other simple NLS1 such as PG 1244+026, RE J1034+396 and RX~J0136 means it is likely that these are also super-Eddington sources. This means their spectral energy distribution cannot be used to determine black hole spin despite appearing to be disc dominated. It also means that the accretion geometry close to the black hole is unlikely to be a flat disc as assumed in the new X-ray reverberation mapping techniques.
Massive stars burn hydrogen through the CNO cycle during most of their evolution. When mixing is efficient, or when mass transfer in binary systems happens, chemically processed material is observed at the surface of O and B stars. ON stars show stronger lines of nitrogen than morphologically normal counterparts. Whether this corresponds to the presence of material processed through the CNO cycle or not is not known. Our goal is to answer this question. We perform a spectroscopic analysis of a sample of ON stars with atmosphere models. We determine the fundamental parameters as well as the He, C, N, and O surface abundances. We also measure the projected rotational velocities. We compare the properties of the ON stars to those of normal O stars. We show that ON stars are usually helium-rich. Their CNO surface abundances are fully consistent with predictions of nucleosynthesis. ON stars are more chemically evolved and rotate - on average - faster than normal O stars. Evolutionary models including rotation cannot account for the extreme enrichment observed among ON main sequence stars. Some ON stars are members of binary systems, but others are single stars as indicated by stable radial velocities. Hence, mass transfer is not a simple explanation for the observed chemical properties. We conclude that ON stars show extreme chemical enrichment at their surface, consistent with nucleosynthesis through the CNO cycle. Its origin is not clear at present.
The space probe 'New Horizon' was launched on 19th of January 2006 in order to study Pluto and its moons. Spacecraft will fly by Pluto as close as 12500 km in the middle of July 2015 and will get the most detailed images of Pluto and its moon until this moment. At the same time, observation obtained by the ground-based telescopes may also be helpful for the research of such distant system. Thereby, the Laboratory of observational astrometry of Pulkovo Observatory of RAS made a decision to reprocess observations obtained during last decade. More than 350 positional observations of Pluto - Charon system were carried out with the mirror astrograph ZA-320M at Pulkovo and Maksutov telescope MTM-500M near Kislovodsk. These observations were processed by means of software system APEX-II developed in Pulkovo observatory and numerical simulation was performed to calculate the differences between positions of photocenter and barycenter of Pluto - Charon system.
We present optical photometric and spectroscopic observations of supernova 2013ej. It is one of the brightest type II supernovae exploded in a nearby ($\sim 10$ Mpc) galaxy NGC 628. The light curve characteristics are similar to type II SNe, but with a relatively shorter ($ \sim85 $ day) and steeper ($ \sim1.7 $ mag (100 d)$^{-1} $ in V) plateau phase. The SN shows a large drop of 2.4 mag in V band brightness during plateau to nebular transition. The absolute ultraviolet (UV) light curves are identical to SN 2012aw, showing a similar UV plateau trend extending up to 85 days. The radioactive $^{56}$Ni mass estimated from the tail luminosity is $ 0.02 $M$_{\odot}$ which is significantly lower than typical type IIP SNe. The characteristics of spectral features and evolution of line velocities indicate that SN 2013ej is a type II event. However, light curve characteristics and some spectroscopic features provide strong support in classifying it as a type IIL event. A detailed SYNOW modelling of spectra indicates the presence of some high velocity components in H$\alpha$ and H$\beta$ profiles, implying possible ejecta-CSM interaction. The nebular phase spectrum shows an unusual notch in the H$\alpha$ emission which may indicate bipolar distribution of $^{56}$Ni. Modelling of the bolometric light curve yields a progenitor mass of $ \sim14 $M$_{\odot}$ and a radius of $ \sim450 $R$_{\odot}$, with a total explosion energy of $ \sim2.3\times10^{51} $ erg.
Binary Systems are the most studied sources of gravitational waves. The mechanisms of emission and the behavior of the orbital parameters are well known and can be written in analytic form in several cases. Besides, the strongest indication of the existence of gravitational waves has arisen from the observation of binary systems. On the other hand, when the detection of gravitational radiation becomes a reality, one of the observed pattern of the signals will be probably of stochastic background nature, which are characterized by a superposition of signals emitted by many sources around the universe. Our aim here is to develop an alternative method of calculating such backgrounds emitted by cosmological compact binary systems during their periodic or quasiperiodic phases. We use an analogy with a problem of Statistical Mechanics in order to perform this sum as well as taking into account the temporal variation of the orbital parameters of the systems. Such a kind of background is of particular importance since it could well form an important foreground for the planned gravitational wave interferometers DECI-Hertz Interferometer Gravitational wave Observatory (DECIGO), Big Bang Observer (BBO), Laser Interferometer Space Antenna (LISA) or Evolved LISA (eLISA), Advanced Laser Interferometer Gravitational-Wave Observatory (ALIGO) and Einstein Telescope (ET).
The overtone and multi-mode RR Lyrae stars in the globular cluster M3 are studied using a 200-d long, $B,V$ and $I_{\mathrm C}$ time-series photometry obtained in 2012. 70\% of the 52 overtone variables observed show some kind of multi-periodicity (additional frequency at ${f_{0.61}}={f_{\mathrm {1O}}}/0.61$ frequency ratio, Blazhko effect, double/multi-mode pulsation, period doubling). A signal at 0.587 frequency ratio to the fundamental-mode frequency is detected in the double-mode star, V13, which may be identified as the second radial overtone mode. If this mode-identification is correct, than V13 is the first RR Lyrae star showing triple-mode pulsation of the first three radial modes. Either the Blazhko effect or the ${f_{0.61}}$ frequency (or both of these phenomena) appear in 7 double-mode stars. The $P_{\mathrm{1O}}/P_{\mathrm{F}}$ period ratio of RRd stars showing the Blazhko effect are anomalous. A displacement of the main frequency component at the fundamental-mode with the value of modulation frequency (or its half) is detected in three Blazhko RRd stars parallel with the appearance of the overtone-mode pulsation. The ${f_{0.61}}$ frequency appears in RRc stars that lie at the blue side of the double-mode region and in RRd stars, raising the suspicion that its occurrence may be connected to double-mode pulsation. The changes of the Blazhko and double-mode properties of the stars are also reviewed using the recent and archive photometric data.
We propose a new multiscale method to calculate the amplitude of the gradient of the linear polarisation vector using a wavelet-based formalism. We demonstrate this method using a field of the Canadian Galactic Plane Survey (CGPS) and show that the filamentary structure typically seen in gradients of linear polarisation maps depends strongly on the instrumental resolution. Our analysis reveals that different networks of filaments are present on different angular scales. The wavelet formalism allows us to calculate the power spectrum of the fluctuations seen in gradients of linear polarisation maps and to determine the scaling behaviour of this quantity. The power spectrum is found to follow a power law with gamma ~ 2.1. We identify a small drop in power between scales of 80 < l < 300 arcmin, which corresponds well to the overlap in the u-v plane between the Effelsberg 100-m telescope and the DRAO 26-m telescope data. We suggest that this drop is due to undersampling present in the 26-m telescope data. In addition, the wavelet coefficient distributions show higher skewness on smaller scales than at larger scales. The spatial distribution of the outliers in the tails of these distributions creates a coherent subset of filaments correlated across multiple scales, which trace the sharpest changes in the polarisation vector P within the field. We suggest that these structures may be associated with highly compressive shocks in the medium. The power spectrum of the field excluding these outliers shows a steeper power law with gamma ~ 2.5.
Establishing the relative role of internally and externally driven mechanisms responsible for disc and bulge growth is essential to understand the evolution of disc galaxies. In this context, we have studied the physical properties of disc galaxies without classical bulges in comparison to those with classical bulges since z~0.9. Using images from the Hubble Space Telescope and Sloan Digital Sky Survey, we have computed both parametric and non-parametric measures, and examined the evolution in size, concentration, stellar mass, effective stellar mass density and asymmetry. We find that both disc galaxies with and without classical bulges have gained more than 50% of their present stellar mass over the last ~8 Gyrs. Also, the increase in disc size is found to be peripheral. While the average total (Petrosian) radius almost doubles from z~0.9 to z~0, the average effective radius undergoes a marginal increase in comparison. Additionally, increase in the density of the inner region is evident through the evolution of both concentration and effective stellar mass density. We find that the asymmetry index falls from higher to lower redshifts, but this is more pronounced for the bulgeless disc sample. Also, asymmetry correlates with the global effective radius, and concentration correlates with the global Sersic index, but better so for higher redshifts only. The substantial increase in mass and size indicates that accretion of external material has been a dominant mode of galaxy growth, where the circumgalactic environment plays a significant role.
We present results from a pulse phase resolved spectroscopy of the complex emission lines around 1 keV in the unique accretion powered X-ray pulsar 4U 1626-67 using observation made with the XMM-Newton in 2003. In this source, the red and blue shifted emission lines and the line widths measured earlier with Chandra suggest their accretion disk origin. Another possible signature of lines produced in accretion disk can be a modulation of the line strength with pulse phase. We found the line fluxes to have pulse phase dependence, making 4U 1626-67 only the second pulsar after Her X-1 to show such variability. The O~VII line at 0.568~keV from 4U 1626-67 varied by a factor of $\sim$4, stronger than the continuum variability, that support their accretion disk origin. The line flux variability may appear due to variable illumination of the accretion disk by the pulsar or more likely, a warp like structure in the accretion disk. We also discuss some further possible diagnostics of the accretion disk in 4U 1626-67 with pulse phase resolved emission line spectroscopy.
The lateral density of a cosmic air shower with a non-zero zenith angle is azimuthally asymmetric. The azimuthal asymmetry consist of a stretching of the iso-density contours to ellipses and to a shift of the center of the elliptic contours with respect to the core of the shower. The aim of the paper is to investigate the shift of the center of the elliptic iso-density contours for different zenith angles . On the basis of a model a qualitative equation is derived for the iso-density contours of inclined showers including the shift. to obtain a quantitative equation MC densities are investigated. The shift can be incorporated in an analytic expression of the azimuthal asymmetry of the lateral density as a function of the polar coordinates and parameterized by the zenith angle. Its predictions for asymmetric lateral densities are compared with densities obtained with MC simulations.
The M 31 nova M31N 2008-12a was recently found to be a recurrent nova (RN) with a recurrence time of about 1 year. This is by far the fastest recurrence time scale of any known RNe. Our optical monitoring programme detected the predicted 2014 outburst of M31N 2008-12a in early October. We immediately initiated an X-ray/UV monitoring campaign with Swift to study the multiwavelength evolution of the outburst. We monitored M31N 2008-12a with daily Swift observations for 20 days after discovery, covering the entire supersoft X-ray source (SSS) phase. We detected SSS emission around day six after outburst. The SSS state lasted for approximately two weeks until about day 19. M31N 2008-12a was a bright X-ray source with a high blackbody temperature. The X-ray properties of this outburst were very similar to the 2013 eruption. Combined X-ray spectra show a fast rise and decline of the effective blackbody temperature. The short-term X-ray light curve showed strong, aperiodic variability which decreased significantly after about day 14. Overall, the X-ray properties of M31N 2008-12a are consistent with the average population properties of M 31 novae. The optical and X-ray light curves can be scaled uniformly to show similar time scales as those of the Galactic RNe U Sco or RS Oph. The SSS evolution time scales and effective temperatures are consistent with a high-mass WD. We predict the next outburst of M31N 2008-12a to occur in autumn 2015.
We use numerical simulations to investigate effect of turbulent velocity on the power spectrum of \HI intensity from external galaxies when (a) all emission is considered, (b) emission with velocity range smaller than the turbulent velocity dispersion is considered. We found that for case (a) the intensity fluctuation depends directly only on the power spectrum of the column density, whereas for case (b) it depends only on the turbulent velocity fluctuation. We discuss the implications of this result in real observations of \HI fluctuations.
We present a precise measurement of downward-going albedo proton fluxes for kinetic energy above $\sim$ 70 MeV performed by the PAMELA experiment at an altitude between 350 and 610 km. On the basis of a trajectory tracing simulation, the analyzed protons were classified into quasi-trapped, concentrating in the magnetic equatorial region, and un-trapped spreading over all latitudes, including both short-lived (precipitating) and long-lived (pseudo-trapped) components. In addition, features of the penumbra region around the geomagnetic cutoff were investigated in detail. PAMELA results significantly improve the characterization of the high energy albedo proton populations at low Earth orbits.
We investigate whether the fuelling of low excitation radio galaxies (LERGs) is linked to major galaxy interactions. Our study utilizes a sample of 10,800 spectroscopic galaxy pairs and 97 post-mergers selected from the Sloan Digital Sky Survey with matches to multi-wavelength datasets. The LERG fraction amongst interacting galaxies is a factor of 3.5 higher than that of a control sample matched in local galaxy density, redshift and stellar mass. However, the LERG excess in pairs does not depend on projected separation and remains elevated out to at least 500 kpc, suggesting that major mergers are not their main fuelling channel. In order to identify the primary fuelling mechanism of LERGs, we compile samples of control galaxies that are matched in various host galaxy and environmental properties. The LERG excess is reduced, but not completely removed, when halo mass or D4000 are included in the matching parameters. However, when BOTH M_halo and D4000 are matched, there is no LERG excess and the 1.4 GHz luminosities (which trace jet mechanical power) are consistent between the pairs and control. In contrast, the excess of optical and mid-IR selected AGN in galaxy pairs is unchanged when the additional matching parameters are implemented. Our results suggest that whilst major interactions may trigger optically and mid-IR selected AGN, the gas which fuels the LERGs has two secular origins: one associated with the large scale environment, such as accretion from the surrounding medium or minor mergers, plus an internal stellar mechanism, such as winds from evolved stars.
We introduce a new set of simulations of a Milky Way like galaxy using the AMR code ART + hydrodynamics in a $\Lambda$CDM cosmogony. The simulation series is named GARROTXA and follow the formation of a late type galaxy from z=60 with a final virial mass of \sim$7.4$\times$10$^{11}$M$_{\odot}$. This system has no major mergers since z=3 and at z=0 becomes a disk late-type spiral galaxy. Several of its large scale properties fall inside recent observational limits of our Galaxy, like the rotation curve shape, the presence of a stellar bar and flare, and a gaseous disk warp, as well as the stellar and baryonic mass. Here, as a first scientific exploitation of the model we study the total amount and spatial distribution of hot X-ray luminous gas. We do not observe in our models a significant presence of a hot gas thick disk as has been recently discussed in observational studies. The analysis of hot gas mock observations (column density and emission measure) revealed that commonly used hypothesis assumed to derive the total hot gas mass in the MW halo from observations lead to biases of about one order of magnitude. Our results suggest that such hot gas distribution is highly anisotropic and that its total mass can account for a non-negligible portion of the missing cosmic baryons leaving still open the contribution of circumgalactic gas. Finally, we have found a clear correlation between the total hot gas mass and the dark matter halo mass of galactic systems. If confirmed, this correlation can become a new method to constrain the total mass of galaxies, in particular the one of the Milky Way.
We aim at determining the spatial distribution of the gas and dust in star-forming regions and address their relative abundances in quantitative terms. We also examine the dust opacity exponent beta for spatial and/or temporal variations. Using mapping observations of the very dense rho Oph A core, we examined standard 1D and non-standard 3D methods to analyse data of far-infrared and submillimeter (submm) continuum radiation. The resulting dust surface density distribution can be compared to that of the gas. The latter was derived from the analysis of accompanying molecular line emission, observed with Herschel from space and with APEX from the ground. As a gas tracer we used N2H+, which is believed to be much less sensitive to freeze-out than CO and its isotopologues. Radiative transfer modelling of the N2H+(J=3-2) and (J=6-5) lines with their hyperfine structure explicitly taken into account provides solutions for the spatial distribution of the column density N(H2), hence the surface density distribution of the gas. The gas-to-dust mass ratio is varying across the map, with very low values in the central regions around the core SM 1. The global average, =88, is not far from the canonical value of 100, however. In rho Oph A, the exponent beta of the power-law description for the dust opacity exhibits a clear dependence on time, with high values of 2 for the envelope-dominated emission in starless Class -1 sources to low values close to 0 for the disk-dominated emission in Class III objects. beta assumes intermediate values for evolutionary classes in between. Since beta is primarily controlled by grain size, grain growth mostly occurs in circumstellar disks. The spatial segregation of gas and dust, seen in projection toward the core centre, probably implies that, like C18O, also N2H+ is frozen onto the grains.
We present a kinematical study of 29 spiral galaxies included in the Spitzer Survey of Stellar Structure in Galaxies, using Halpha Fabry-Perot data obtained with the Galaxy Halpha Fabry-Perot System instrument at the William Herschel Telescope in La Palma, complemented with images in the R-band and in Halpha. The primary goal is to study the evolution and properties of the main structural components of galaxies through the kinematical analysis of the FP data, complemented with studies of morphology, star formation and mass distribution. In this paper we describe how the FP data have been obtained, processed and analysed. We present the resulting moment maps, rotation curves, velocity model maps and residual maps. Images are available in FITS format through the NASA/IPAC Extragalactic Database and the Centre de Donn\'ees Stellaires. With these data products we study the non-circular motions, in particular those found along the bars and spiral arms. The data indicate that the amplitude of the non-circular motions created by the bar does not correlate with the bar strength indicators. The amplitude of those non-circular motions in the spiral arms does not correlate with either arm class or star formation rate along the spiral arms. This implies that the presence and the magnitude of the streaming motions in the arms is a local phenomenon.
We measure the clustering of X-ray, radio, and mid-IR-selected active galactic nuclei (AGN) at 0.2 < z < 1.2 using multi-wavelength imaging and spectroscopic redshifts from the PRIMUS and DEEP2 redshift surveys, covering 7 separate fields spanning ~10 square degrees. Using the cross-correlation of AGN with dense galaxy samples, we measure the clustering scale length and slope, as well as the bias, of AGN selected at different wavelengths. Similar to previous studies, we find that X-ray and radio AGN are more clustered than mid-IR-selected AGN. We further compare the clustering of each AGN sample with matched galaxy samples designed to have the same stellar mass, star formation rate, and redshift distributions as the AGN host galaxies and find no significant differences between their clustering properties. The observed differences in the clustering of AGN selected at different wavelengths can therefore be explained by the clustering differences of their host populations, which have different distributions in both stellar mass and star formation rate. Selection biases inherent in AGN selection, therefore, determine the clustering of observed AGN samples. We further find no significant difference between the clustering of obscured and unobscured AGN, using IRAC or WISE colors or X-ray hardness ratio.
Numerical simulations show the formation of self-gravitating primordial disks during the assembly of the first structures in the Universe, in particular during the formation of Pop. III and supermassive stars. Their subsequent evolution is expected to be crucial to determine the mass scale of the first cosmological objects, which depends on the temperature of the gas and the dominant cooling mechanism. Here, we derive a one-zone framework to explore the chemical evolution of such disks and show that viscous heating leads to the collisional dissociation of an initially molecular gas. The effect is relevant on scales of 10 AU (1000 AU) for a central mass of 10 M_solar (10^4 M_solar) at an accretion rate of 0.1 M_solar/yr, and provides a substantial heat input to stabilize the disk. If the gas is initially atomic, it remains atomic during the further evolution, and the effect of viscous heating is less significant. The additional thermal support is particularly relevant for the formation of very massive objects, such as the progenitors of the first supermassive black holes. The stabilizing impact of viscous heating thus alleviates the need for a strong radiation background as a means of keeping the gas atomic.
We present a detailed analysis of the local evolution of 206 Lagrangian
Volumes (LVs) selected at high redshift around galaxy seeds, identified in a
large-volume $\Lambda$CDM hydrodynamical simulation. The LVs have a mass range
of $1 - 1500 \times 10^{10} M_\odot$. We follow the dynamical evolution of the
density field inside these initially spherical LVs from $z=10$ up to $z_{\rm
low}= 0.05$, witnessing highly non-linear, anisotropic mass rearrangements
within them, leading to the emergence of the local cosmic web (CW). These mass
arrangements have been analysed in terms of the reduced inertia tensor
$I_{ij}^r$, focusing on the evolution of the principal axes of inertia and
their corresponding eigen directions, and paying particular attention to the
times when the evolution of these two structural elements declines. In
addition, mass and component effects along this process have also been
investigated.
We have found that deformations are led by DM dynamics and they transform
most of the initially spherical LVs into prolate shapes, i.e. filamentary
structures. An analysis of the individual freezing-out time distributions for
shapes and eigen directions shows that first most of the LVs fix their three
axes of symmetry (like a skeleton), while accretion flows towards them still
continue. Very remarkably, we have found that more massive LVs fix their
skeleton earlier on than less massive ones. We briefly discuss the
astrophysical implications our findings could have, including the galaxy
mass-morphology relation and the effects on the galaxy-galaxy merger parameter
space, among others.
We present newly processed archival Herschel images of molecular cloud MCLD 123.5+24.9 in the Polaris Flare. This cloud contains five starless cores. Using the spectral synthesis code Cloudy, we explore uncertainties in the derivation of column densities, hence, masses of molecular cores from Herschel data. We first consider several detailed grain models that predict far-IR grain opacities. Opacities predicted by the models differ by more than a factor of two, leading to uncertainties in derived column densities by the same factor. Then we consider uncertainties associated with the modified blackbody fitting process used by observers to estimate column densities. For high column density clouds (N(H) $\gg$ 10$^{22}$ cm$^{-2}$), this fitting technique can underestimate column densities by about a factor of three. Finally, we consider the virial stability of the five starless cores in MCLD 123.5+24.9. All of these cores appear to have strongly sub-virial masses, assuming, as we argue, that $^{13}$CO line data provide reliable estimates of velocity dispersions. Evidently, they are not self-gravitating, so it is no surprise that they are starless.
It is shown that holographic cosmology implies an evolving Hubble radius $c^{-1}\dot{R}_H = -1 + 3\Omega_m$ in the presence of a dimensionless matter density $\Omega_m$ scaled to the closure density $3H^2/8\pi G$, where $c$ denotes the velocity of light and $H$ and $G$ denote the Hubble parameter and Newton's constant. It reveals a dynamical dark energy and a sixfold increase in gravitational attraction to matter on the scale of the Hubble acceleration. It reproduces the transition redshift $z_t\simeq 0.4$ to the present epoch of accelerated expansion and is consistent with $(q_0,(dq/dz)_0)$ of the deceleration parameter $q(z)=q_0+(dq/dz)_0z$ observed in Type Ia supernovae.
We consider the cross-correlation search for periodic GWs and its potential application to the LMXB Sco X-1. This method coherently combines data from different detectors at the same time, as well as different times from the same or different detectors. By adjusting the maximum time offset between a pair of data segments to be coherently combined, one can tune the method to trade off sensitivity and computing costs. In particular, the detectable signal amplitude scales as the inverse fourth root of this coherence time. The improvement in amplitude sensitivity for a search with a coherence time of 1hr, compared with a directed stochastic background search with 0.25Hz wide bins is about a factor of 5.4. We show that a search of 1yr of data from Advanced LIGO and Advanced Virgo with a coherence time of 1hr would be able to detect GWs from Sco X-1 at the level predicted by torque balance over a range of signal frequencies from 30-300Hz; if the coherence time could be increased to 10hr, the range would be 20-500Hz. In addition, we consider several technical aspects of the cross-correlation method: We quantify the effects of spectral leakage and show that nearly rectangular windows still lead to the most sensitive search. We produce an explicit parameter-space metric for the cross-correlation search in general and as applied to a neutron star in a circular binary system. We consider the effects of using a signal template averaged over unknown amplitude parameters: the search is sensitive to a combination of the intrinsic signal amplitude and the inclination of the neutron star rotation axis, and the peak of the expected detection statistic is systematically offset from the true signal parameters. Finally, we describe the potential loss of SNR due to unmodelled effects such as signal phase acceleration within the Fourier transform timescale and gradual evolution of the spin frequency.
This is the summary of two lectures that aim to give an overview of cosmology. I will not try to be too rigorous in derivations, nor to give a full historical overview. The idea is to provide a "taste" of cosmology and some of the interesting topics it covers. The standard cosmological model is presented and I highlight the successes of cosmology over the past decade or so. Keys to the development of the standard cosmological model are observations of the cosmic microwave background and of large-scale structure, which are introduced. Inflation and dark energy and the outlook for the future are also discussed. Slides from the lectures are available from the school website: physicschool.web.cern.ch/PhysicSchool/CLASHEP/CLASHEP2011/.
We construct the minimal effective field theory (EFT) of supersymmetric inflation, whose field content is a real scalar, the goldstone for time-translation breaking, and a Weyl fermion, the goldstino for supersymmetry (SUSY) breaking. The inflating background can be viewed as a single SUSY-breaking sector, and the degrees of freedom can be efficiently parameterized using constrained superfields. Our EFT is comprised of a chiral superfield X_NL containing the goldstino and satisfying X_NL^2 = 0, and a real superfield B_NL containing both the goldstino and the goldstone, satisfying X_NL B_NL = B_NL^3 = 0. We match results from our EFT formalism to existing results for SUSY broken by a fluid background, showing that the goldstino propagates with subluminal velocities. The same effect can also be derived from the unitary gauge gravitino action after embedding our EFT in supergravity. If the gravitino mass is comparable to the Hubble scale during inflation, we identify a new parameter in the EFT related to a time-dependent phase of the gravitino mass parameter. We briefly comment on the leading contributions of goldstino loops to inflationary observables.
Due to their superfluid properties some compact astrophysical objects, like neutron or quark stars, may contain a significant part of their matter in the form of a Bose-Einstein Condensate. Observationally distinguishing between neutron/quark stars and Bose-Einstein Condensate stars is a major challenge for this latter theoretical model. An observational possibility of indirectly distinguishing Bose-Einstein Condensate stars from neutron/quark stars is through the study of the thin accretion disks around compact general relativistic objects. In the present paper, we perform a detailed comparative study of the electromagnetic and thermodynamic properties of the thin accretion disks around rapidly rotating Bose-Einstein Condensate stars, neutron stars and quark stars, respectively. Due to the differences in the exterior geometry, the thermodynamic and electromagnetic properties of the disks (energy flux, temperature distribution, equilibrium radiation spectrum and efficiency of energy conversion) are different for these classes of compact objects. Hence in this preliminary study we have pointed out some astrophysical signatures that may allow to observationally discriminate between Bose-Einstein Condensate stars and neutron/quark stars, respectively.
The Tayler instability is a kink-type, current driven instability that plays an important role in plasma physics but might also be relevant in liquid metal applications with high electrical currents. In the framework of the Tayler-Spruit dynamo model of stellar magnetic field generation, the question of spontaneous helical (chiral) symmetry breaking during the saturation of the Tayler instability has received considerable interest. Focusing on fluids with low magnetic Prandtl numbers, for which the quasistatic approximation can be applied, we utilize an integro-differential equation approach in order to investigate the saturation mechanism of the Tayler instability. Both the exponential growth phase and the saturated phase are analyzed in terms of the action of the alpha and beta effects of mean-field magnetohydrodynamics. In the exponential growth phase we always find a spontaneous chiral symmetry breaking which, however, disappears in the saturated phase. For higher degrees of supercriticality, we observe helicity oscillations in the saturated regime. For Lundquist numbers in the order of one we also obtain chiral symmetry breaking of the saturated magnetic field.
We test and compare a number of existing models predicting the location of magnetic reconnection at Earth's dayside magnetopause for various solar wind conditions. We employ robust image processing techniques to determine the locations where each model predicts reconnection to occur. The predictions are then compared to the magnetic separators, the magnetic field lines separating different magnetic topologies. The predictions are tested in distinct high-resolution simulations with interplanetary magnetic field (IMF) clock angles ranging from 30 to 165 degrees in global magnetohydrodynamic simulations using the three-dimensional Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) code with a uniform resistivity, although the described techniques can be generally applied to any self-consistent magnetosphere code. Additional simulations are carried out to test location model dependence on IMF strength and dipole tilt. We find that most of the models match large portions of the magnetic separators when the IMF has a southward component, with the models saying reconnection occurs where the local reconnection rate and reconnection outflow speed are maximized performing best. When the IMF has a northward component, none of the models tested faithfully map the entire magnetic separator, but the maximum magnetic shear model is the best at mapping the separator in the cusp region where reconnection has been observed. Predictions for some models with northward IMF orientations improve after accounting for plasma flow shear parallel to the reconnecting components of the magnetic fields. Implications for observations are discussed.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
