We employ the Planck 2013 CMB temperature anisotropy and lensing data, and baryon acoustic oscillation (BAO) data to constrain a phenomenological $w$CDM model, where dark matter and dark energy interact. We assume time-dependent equation of state parameter for dark energy, and treat dark matter and dark energy as fluids whose energy-exchange rate is proportional to the dark-matter density. The CMB data alone leave a strong degeneracy between the interaction rate and the physical CDM density parameter today, $\omega_c$, allowing a large interaction rate $|\Gamma| \sim H_0$. However, as has been known for a while, the BAO data break this degeneracy. Moreover, we exploit the CMB lensing potential likelihood, which probes the matter perturbations at redshift $z \sim 2$ and is very sensitive to the growth of structure, and hence one of the tools for discerning between the $\Lambda$CDM model and its alternatives. However, we find that in the non-phantom models ($w_{\mathrm{de}}>-1$), the constraints remain unchanged by the inclusion of the lensing data and consistent with zero interaction, $-0.14 < \Gamma/H_0 < 0.02$ at 95\% CL. On the contrary, in the phantom models ($w_{\mathrm{de}}<-1$), energy transfer from dark energy to dark matter is moderately favoured over the non-interacting model; $-0.57 < \Gamma/H_0 < -0.10$ at 95\% CL with CMB+BAO, while addition of the lensing data shifts this to $-0.46 < \Gamma/H_0 < -0.01$.
We present FIRE/Gizmo hydrodynamic zoom-in simulations of isolated dark
matter halos, two each at the mass of classical dwarf galaxies ($M_{\rm vir}
\simeq 10^{10} M_{\odot}$) and ultra-faint galaxies ($M_{\rm vir} \simeq 10^9
M_{\odot}$), and with two feedback implementations. The resultant central
galaxies lie on an extrapolated abundance matching relation from $M_{\star}
\simeq 10^6$ to $10^4 M_{\odot}$ without a break. Every host is filled with
subhalos, many of which form stars. Our dwarfs with $M_{\star} \simeq 10^6
M_{\odot}$ each have 1-2 well-resolved satellites with $M_{\star} = 3-200
\times 10^3 M_{\odot}$. Even our isolated ultra-faint galaxies have
star-forming subhalos. If this is representative, dwarf galaxies throughout the
universe should commonly host tiny satellite galaxies of their own. We combine
our results with the ELVIS simulations to show that targeting $\sim 50~ \rm
kpc$ regions around nearby isolated dwarfs could increase the chances of
discovering ultra-faint galaxies by $\sim 35\%$ compared to random halo
pointings, and specifically identify the region around the Phoenix dwarf galaxy
as a good potential target.
The well-resolved ultra-faint galaxies in our simulations ($M_{\star} \simeq
3 - 30 \times 10^3 M_{\odot}$) form within $M_{\rm peak} \simeq 0.5 - 3 \times
10^9 M_{\odot}$ halos. Each has a uniformly ancient stellar population ($ > 10~
\rm Gyr$) owing to reionization-related quenching. More massive systems, in
contrast, all have late-time star formation. Our results suggest that $M_{\rm
halo} \simeq 5 \times 10^9 M_{\odot}$ is a probable dividing line between halos
hosting reionization "fossils" and those hosting dwarfs that can continue to
form stars in isolation after reionization.
Using astrometric measurements obtained with the FORS2/VLT camera, we are searching for low-mass companions around 20 nearby ultracool dwarfs. With a single-measurement precision of 0.1 milli-arcseconds, our survey is sensitive to a wide range of companion masses from planetary companions to binary systems. Here, we report the discovery and orbit characterisation of a new ultracool binary at a distance of 19.5 pc from Earth that is composed of the M8.5-dwarf primary DE0630-18 and a substellar companion. The nearly edge-on orbit is moderately eccentric (e=0.23) with an orbital period of 1120 d, which corresponds to a relative separation in semimajor axis of approximately 1.1 AU. We obtained a high-resolution optical spectrum with UVES/VLT and measured the system's heliocentric radial velocity. The spectrum does not exhibit lithium absorption at 670.8 nm, indicating that the system is not extremely young. A preliminary estimate of the binary's physical parameters tells us that it is composed of a primary at the stellar-substellar limit and a massive brown-dwarf companion. DE0630-18 is a new very low-mass binary system with a well-characterised orbit.
Dwarf galaxies generally follow a mass-metallicity (MZ) relation, where more massive objects retain a larger fraction of heavy elements. Young tidal dwarf galaxies (TDGs), born in the tidal tails produced by interacting gas-rich galaxies, have been thought to not follow the MZ relation, because they inherit the metallicity of the more massive parent galaxies. We present chemical evolution models to investigate if TDGs that formed at very high redshifts, where the metallicity of their parent galaxy was very low, can produce the observed MZ relation. Assuming that galaxy interactions were more frequent in the denser high-redshift universe, TDGs could constitute an important contribution to the dwarf galaxy population. The survey of chemical evolution models of TDGs presented here captures for the first time an initial mass function (IMF) of stars that is dependent on both the star formation rate and the gas metallicity via the integrated galactic IMF (IGIMF) theory. As TDGs form in the tidal debris of interacting galaxies, the pre-enrichment of the gas, an underlying pre-existing stellar population, infall, and mass dependent outflows are considered. The models of young TDGs that are created in strongly pre-enriched tidal arms with a pre-existing stellar population can explain the measured abundance ratios of observed TDGs. The same chemical evolution models for TDGs, that form out of gas with initially very low metallicity, naturally build up the observed MZ relation. The modelled chemical composition of ancient TDGs is therefore consistent with the observed MZ relation of satellite galaxies.
We examine the effects of gas-expulsion on initially substructured distributions of stars. We perform N-body simulations of the evolution of these distributions in a static background potential to mimic the gas. We remove the static potential instantaneously to model gas-expulsion. We find that the exact dynamical state of the cluster plays a very strong role in affecting a cluster's survival, especially at early times: they may be entirely destroyed or only weakly affected. We show that knowing both detailed dynamics and relative star-gas distributions can provide a good estimate of the post-gas expulsion state of the cluster, but even knowing these is not an absolute way of determining the survival or otherwise of the cluster.
Studies of Fermi data indicate an excess of GeV gamma rays around the Galactic center (GC), possibly due to dark matter. We show that young gamma-ray pulsars can yield a similar signal. First, a high concentration of GC supernovae naturally leads to a population of kicked pulsars symmetric about the GC. Second, while very-young pulsars with soft spectra reside near the Galactic plane, pulsars with spectra that have hardened with age accumulate at larger angles. This combination, including unresolved foreground pulsars, traces the morphology and spectrum of the Excess.
We present rest-frame NIR luminosities and stellar masses for a large and uniformly-selected population of GRB host galaxies using deep Spitzer Space Telescope imaging of 117 targets from the Swift GRB Host Galaxy Legacy Survey spanning 0.03 < z < 6.3, and determine the effects of galaxy evolution and chemical enrichment on the mass distribution of the GRB host population across cosmic history. We find strong evolution in the host luminosity distribution between z~0.5 (median absolute NIR AB magnitude ~ -18.5, corresponding to M* ~ 3x10^8 M_sun and z~1.5), but negligible variation between z~1.5 and z~5 (median magnitude ~ -21.2, corresponding to M* ~ 5x10^9 M_sun). Dust-obscured GRBs dominate the massive host population but are only rarely seen associated with low-mass hosts, indicating that massive star-forming galaxies are universally and (to some extent) homogeneously dusty at high-redshift while low-mass star-forming galaxies retain little dust in their ISM. Comparing our luminosity distributions to field surveys and measurements of the high-z mass-metallicity relation, our results have good consistency with a model in which the GRB rate per unit star-formation is constant in galaxies with gas-phase metallicity below approximately the Solar value but heavily suppressed in more metal-rich environments. This model also naturally explains the previously-reported "excess" in the GRB rate beyond z>2; metals stifle GRB production in most galaxies at z<1.5 but have only minor impact at higher redshifts. The metallicity threshold we infer is much higher than predicted by single-star models and favors a binary progenitor. Our observations also constrain the fraction of cosmic star-formation in low-mass galaxies undetectable to Spitzer to be a small minority at most redshifts (~10% at z~2, ~25% at z~3, and ~50% at z=3.5-6.0).
In earlier work we showed that a dark matter halo with a virial mass of $10^7$ M$_\odot$ can survive feedback from its own massive stars and form stars for $\gtrsim100$ Myr. We also found that our modelled systems were consistent with observations of ultrafaint dwarfs (UFDs), the least massive known galaxies. Very metal-poor damped Lyman-$\alpha$ systems (DLAs) recently identified at $z\sim2$ may represent the gas that formed at least some of the observed stars in UFDs. We compare projected sightlines from our simulations to the observed metal-poor DLAs and find that our models can reach the densities of the observed sightlines; however the metallicities are inconsistent with the single supernova simulations, suggesting enrichment by multiple supernovae. We model two scenarios for the history of these systems. The first explains the gas abundances in DLAs by a single burst of star formation. This model can produce the observed DLA abundances, but does not provide an explanation as to why the DLAs show suppressed [$\alpha$/Fe] compared to the stellar population of UFDs. The second scenario splits the DLAs into a population which is enriched by a single burst, and a population that is enriched by a second burst after the accretion of metal-poor gas. In this scenario, the suppressed average [$\alpha$/Fe] in DLAs compared to UFDs results from enrichment of second-burst systems by Type Ia supernovae.
According to our current cosmological model, galaxies like the Milky Way are expected to experience many mergers over their lifetimes. The most massive of the merging galaxies will be dragged towards the disc-plane, depositing stars and dark matter into an accreted disc structure. In this work, we utilize the chemo-dynamical template developed in Ruchti et al. to hunt for accreted stars. We apply the template to a sample of 4,675 stars in the third internal data release from the Gaia-ESO Spectroscopic Survey. We find a significant component of accreted halo stars, but find no evidence of an accreted disc component. This suggests that the Milky Way has had a rather quiescent merger history since its disc formed some 8-10 billion years ago and therefore possesses no significant dark matter disc.
We introduce the Swift Gamma-Ray Burst Host Galaxy Legacy Survey ("SHOALS"), a multi-observatory high-redshift galaxy survey targeting the largest unbiased sample of long-duration gamma-ray burst hosts yet assembled (119 in total). We describe the motivations of the survey and the development of our selection criteria, including an assessment of the impact of various observability metrics on the success rate of afterglow-based redshift measurement. We briefly outline our host-galaxy observational program, consisting of deep Spitzer/IRAC imaging of every field supplemented by similarly-deep, multi-color optical/NIR photometry, plus spectroscopy of events without pre-existing redshifts. Our optimized selection cuts combined with host-galaxy follow-up have so far enabled redshift measurements for 110 targets (92%) and placed upper limits on all but one of the remainder. About 20% of GRBs in the sample are heavily dust-obscured, and at most 2% originate from z>5.5. Using this sample we estimate the redshift-dependent GRB rate density, showing it to peak at z~2.5 and fall by about an order of magnitude towards low (z=0) redshift, while declining more gradually towards high (z~7) redshift. This behavior is consistent with a progenitor whose formation efficiency varies modestly over cosmic history. Our survey will permit the most detailed examination to date of the connection between the GRB host population and general star-forming galaxies, directly measure evolution in the host population over cosmic time and discern its causes, and provide new constraints on the fraction of cosmic star-formation occurring in undetectable galaxies at all redshifts.
We examine stellar population gradients in ~100 massive early type galaxies spanning 180 < sigma* < 370 km/s and M_K of -22.5 to -26.5 mag, observed as part of the MASSIVE survey (Ma et al. 2014). Using integral-field spectroscopy from the Mitchell Spectrograph on the 2.7m telescope at McDonald Observatory, we create stacked spectra as a function of radius for galaxies binned by their stellar velocity dispersion, stellar mass, and group richness. With excellent sampling at the highest stellar mass, we examine radial trends in stellar population properties extending to beyond twice the effective radius (~2.5 R_e). We examine radial trends in age, metallicity, and in the abundance ratios of Mg, C, N, and Ca, and discuss the implications for star formation histories and elemental yields. When weighted to a fixed physical radius of 3-6 kpc (the likely size of the galaxy cores formed at high redshift) stellar age and [Mg/Fe] increase with increasing sigma* and depend only weakly on stellar mass, as we might expect if denser galaxies form their central cores earlier and faster. If we instead focus on trends weighted towards R_e, the trends in abundance and abundance ratio weaken, as might be expected if the stars at large radius were accreted by smaller galaxies. Finally, we show that when controlling for sigma*, there are only very subtle differences in stellar population properties or gradients as a function of group richness; even at large radius internal properties matter more than environment in determining star formation history.
The Gaia-ESO Survey has recently unveiled the complex kinematic signature of the Gamma Velorum cluster: this cluster is composed of two kinematically distinct populations (hereafter, population A and B), showing two different velocity dispersions and a relative ~2 km s^-1 radial velocity (RV) shift. In this paper, we propose that the two populations of the Gamma Velorum cluster originate from two different sub-clusters, born from the same parent molecular cloud. We investigate this possibility by means of direct-summation N-body simulations. Our scenario is able to reproduce not only the RV shift and the different velocity dispersions, but also the different centroid (~0.5 pc), the different spatial concentration and the different line-of-sight distance (~5 pc) of the two populations. The observed 1-2 Myr age difference between the two populations is also naturally explained by our scenario, in which the two sub-clusters formed in two slightly different star formation episodes. Our simulations suggest that population B is strongly supervirial, while population A is close to virial equilibrium. We discuss the implications of our models for the formation of young star clusters and OB associations in the Milky Way.
The lack of detected pulsars at the Galactic Center (GC) region is a long-standing mystery. We argue that the high stellar density in the central parsec around the GC is likely to result in a pulsar population dominated by millisecond pulsars (MSPs), similar to the situation in globular cluster environments. Earlier GC pulsar searches have been largely insensitive to such an MSP population, accounting for the lack of pulsar detections. We estimate the best search frequency for such an MSP population with present and upcoming broad-band radio telescopes for two possible scattering scenarios, the "weak-scattering" case suggested by the recent detection of a magnetar close to the GC, and the "strong-scattering" case, with the scattering screen located close to the GC. The optimal search frequencies are $\approx 8$ GHz (weak-scattering) and $\approx 25$ GHz (strong-scattering), for pulsars with periods 1-20 ms, assuming that GC pulsars have a luminosity distribution similar to that those in the rest of the Milky Way. We find that 10-30 hour integrations with the Very Large Array and the Green Bank Telescope would be sufficient to detect MSPs at the GC distance in the weak-scattering case. However, if the strong-scattering case is indeed applicable to the GC, observations with the full Square Kilometre Array would be needed to detect the putative MSP population.
We develop an analytic model for galaxy intrinsic alignments (IA) based on the theory of tidal alignment. We calculate all relevant nonlinear corrections at one-loop order, including effects from nonlinear density evolution, galaxy biasing, and source density weighting. Contributions from density weighting are found to be particularly important and lead to bias dependence of the IA amplitude, even on large scales. This effect may be responsible for much of the luminosity dependence in IA observations. The increase in IA amplitude for more highly biased galaxies reflects their locations in regions with large tidal fields. We also consider the impact of smoothing the tidal field on halo scales. We compare the performance of this consistent nonlinear model in describing the observed alignment of luminous red galaxies with the linear model as well as the frequently used "nonlinear alignment model," finding a significant improvement on small and intermediate scales. We also show that the cross-correlation between density and IA (the "GI" term) can be effectively separated into source alignment and source clustering, and we accurately model the observed alignment down to the one-halo regime using the tidal field from the fully nonlinear halo-matter cross correlation. Inside the one-halo regime, the average alignment of galaxies with density tracers no longer follows the tidal alignment prediction, likely reflecting nonlinear processes that must be considered when modeling IA on these scales. Finally, we discuss tidal alignment in the context of cosmic shear measurements.
Ellerman bombs (EBs) are tiny brightenings often observed near sunspots. The most impressive characteristic of the EB spectra is the two emission bumps in both wings of the H$\alpha$ and \ion{Ca}{II} 8542 {\AA} lines. High-resolution spectral data of three small EBs were obtained on 2013 June 6 with the largest solar telescope, the 1.6 meter New Solar Telescope (NST), at the Big Bear Solar Observatory. The characteristics of these EBs are analyzed. The sizes of the EBs are in the range of 0.3\arcsec\--0.8\arcsec\ and their durations are only 3--5 minutes. Our semi-empirical atmospheric models indicate that the heating occurs around the temperature minimum region with a temperature increase of 2700--3000 K, which is surprisingly higher than previously thought. The radiative and kinetic energies are estimated to be as high as 5$\times$10$^{25}$--3.0$\times$10$^{26}$ ergs despite the small size of these EBs. Observations of the magnetic field show that the EBs appeared just in a parasitic region with mixed polarities and accompanied by mass motions. Nonlinear force-free field extrapolation reveals that the three EBs are connected with a series of magnetic field lines associated with bald patches, which strongly implies that these EBs should be produced by magnetic reconnection in the solar lower atmosphere. According to the lightcurves and the estimated magnetic reconnection rate, we propose that there is a three phase process in EBs: pre-heating, flaring and cooling phases.
Simulations of the formation of large-scale structure predict that dark matter, low density highly ionized gas, and galaxies form 10 40 Mpc scale filaments. These structure are easily recognized in the distribution of galaxies, but have not been directly observed in the distribution of the gas. We use Ly-alpha absorption lines in the spectra of 24 AGN to present a new way to probe these filaments. We use a new catalogue of nearby (cz<10,000 km/s) galaxies, complete down to a luminosity of about 0.05 L* for the region of space analyzed here. Using HST spectra of 24 AGN we sample the gas associated with a 30x5 Mpc galaxy filament at cz~3500 km/s. All of our sightlines pass outside the virial radius of any known filament galaxy. Within 500 kpc of the filament axis the detection rate is ~80%, while no detections are seen more than 2.1 Mpc from the filament. The width of the Lya lines correlates with filament impact parameter and the four BLAs in our sample all occur within 400 kpc of the filament axis, indicating increased temperature and/or turbulence. Comparing to simulations, we find that the recent Haardt & Madau (2012) extragalactic ionizing background predicts a factor 3-5 too few ionizing photons. Using a more intense radiation field matches the hydrogen density profile within 2.1 Mpc of the filament axis, but the simulations still overpredict the detection rate between 2.1 and 5 Mpc from the axis.
To test the theoretical understanding that finding bright CO emission depends primarily on dust shielding, we investigate the relationship between CO emission ($I_{\rm CO}$) and the amount of dust (estimated from IR emission and expressed as "$A_V$") across the Large Magellanic Cloud, the Small Magellanic Cloud, and the Milky Way. We show that at our common resolution of 10 pc scales, $I_{\rm CO}$ given a fixed line-of-sight $A_V$ is similar across all three systems despite the difference in metallicity. We find some evidence for a secondary dependence of $I_{\rm CO}$ on radiation field; in the LMC, $I_{\rm CO}$ at a given $A_V$ is smaller in regions of high $T_{\rm dust}$, perhaps because of an increased photodissociating radiation field. We suggest a simple but useful picture in which the CO-to-H$_2$ conversion factor (\xco) depends on two separable factors: (1) the distribution of gas column densities, which maps to an extinction distribution via a dust-to-gas ratio; and (2) the dependence of $I_{\rm CO}$ on $A_V$. Assuming that the probability distribution function (PDF) of local Milky Way clouds is universal, this approach predicts a dependence of \xco\ on $Z$ between $Z^{-1}$ and $Z^{-2}$ above about a third solar metallicity. Below this metallicity, CO emerges from only the high column density parts of the cloud and so depends very sensitively on the adopted PDF and the H$_2$/{\sc Hi} prescription. The PDF of low metallicity clouds is thus of considerable interest and the uncertainty associated with even an ideal prescription for \xco\ at very low metallicity will be large.
Large-scale and deep sky survey missions are rapidly collecting a large amount of stellar spectra, which necessitate the estimation of atmospheric parameters directly from spectra and makes it feasible to statistically investigate latent principles in a large dataset. We present a technique for estimating parameters $T_{eff}$, log$~g$ and [Fe/H] from stellar spectra. With this technique, we first extract features from stellar spectra using the LASSO algorithm; then, the parameters are estimated from the extracted features using the SVR. On a subsample of 20~000 stellar spectra from SDSS with reference parameters provided by SDSS/SEGUE Pipeline SSPP, estimation consistency are 0.007458 dex for log$~T_{eff}$ (101.609921 K for $T_{eff}$), 0.189557 dex for log$~g$ and 0.182060 for [Fe/H], where the consistency is evaluated by mean absolute error. Prominent characteristics of the proposed scheme are sparseness, locality, and physical interpretability. In this work, every spectrum consists of 3821 fluxes, and 10, 19, and 14 typical wavelength positions are detected respectively for estimating $T_{eff}$, log$~g$ and [Fe/H]. It is shown that the positions are related to typical lines of stellar spectra. This characteristic is important in investigating physical indications from analysis results. Then, stellar spectra can be described by the individual fluxes on the detected positions (PD) or local integration of fluxes near them (LI). The abovementioned consistency is the result based on features described by LI. If features are described by PD, consistency are 0.009092 dex for log$~T_{eff}$ (124.545075 K for $T_{eff}$), 0.198928 dex for log$~g$, and 0.206814 dex for [Fe/H].
The alignment of Saturn's magnetic pole with its rotation axis precludes the use of magnetic field measurements to determine its rotation period. The period was previously determined from radio measurements by the Voyager spacecraft to be 10h 39m 22.4s. When the Cassini spacecraft measured a period of 10h 47m 6s, which was additionally found to change between sequential measurements, it became clear that the radio period could not be used to determine the bulk planetary rotation period. Estimates based upon Saturn's measured wind fields have increased the uncertainty even more, giving numbers smaller than the Voyager rotation period, and at present Saturn's rotation period is thought to be between 10h 32m and 10h 47m, which is unsatisfactory for such a fundamental property. Here we report a period of 10h 32m 45s +- 46s, based upon an optimization approach using Saturn's measured gravitational field and limits on the observed shape and possible internal density profiles. Moreover, even when solely using the constraints from its gravitational field, the rotation period can be inferred with a precision of several minutes. To validate our method, we applied the same procedure to Jupiter and correctly recovered its well-known rotation period.
The color-magnitude diagram (CMD) of globular cluster NGC1651 has special structures including a broad main sequence, an extended main sequence turn-off and an extended red giant clump. The reason for such special CMDs remains unclear. In order to test how different the results from various stellar population assumptions are, we study a high-quality CMD of NGC1651 from the Hubble Space Telescope archive via eight kinds of models. Distance modulus, extinction, age ranges, star formation mode, fraction of binaries, and fraction of rotational stars are determined and then compared. The results show that stellar populations both with and without age spread can reproduce the special structure of the observed CMD. A composite population with extended star formation from 1.8\,Gyrs ago to 1.4\,Gyrs ago, which contains 50 per cent binaries and 70 per cent rotational stars, fits the observed CMD best. Meanwhile, a 1.5\,Gyr-old simple population that consists of rotational stars can also fit the observed CMD well. The results of CMD fitting are shown to depend strongly on stellar population type (simple or composite), and fraction of rotators. If the member stars of NGC1651 formed in a single star burst, the effect of stellar rotation should be very important for the explanation of observed CMDs. Otherwise, the effect may be small. It is also possible that the special observed CMD is a result of the combined effects of stellar binarity, rotation and age spread. Therefore, further work on stellar population type and fraction of rotational stars of intermediate-age clusters are necessary to understand their observed CMDs.
We study the high-quality CMDs of three star clusters, NGC 1831, NGC 1868 and NGC 2249 in detail, via the two most likely causes (stellar rotation and age spread) for CMDs with extended main-sequence turn-offs. The results show evident failure of stellar rotation (including resolved and unresolved binary stars) to interpret the CMDs of three star clusters, and the unexpected success of age spread. In particular, the special structures of turn-off and red clump parts cannot be generated by stellar rotation, but age spread perfectly reproduces the observed features. This suggests that these three clusters contain multiple populations of stars, rather than a single population of rotating stars and binaries. The thick subgiant branch of NGC 1831 gives the strongest support to this. The results demonstrate that stellar rotation cannot save the widely accepted view (simple population) of star clusters, and an extended star formation history is needed for explaining the observed CMDs. In addition, this work shows that a narrow subgiant branch does not correspond to a simple population. The judgement of simple population in NGC 1651 by a previous work is not necessarily reliable. We anticipate our assay to be a starting point for more precise study of Hertzsprung-Russell diagrams of star clusters. This involves many factors, such as binaries, rotating stars, and star formation history (SFH).
We present a fast iterative FFT-based reconstruction algorithm that allows for non- parallel redshift-space distortions (RSD). We test our algorithm on both N-body dark matter simulations and mock distributions of galaxies designed to replicate galaxy survey conditions. We compare solenoidal and irrotational components of the redshift distortion and show that an approximation of this distortion leads to a better estimate of the real-space potential (and therefore faster convergence) than ignoring the RSD when estimating the displacement field. Our iterative reconstruction scheme converges in two iterations for the mock samples corresponding to BOSS CMASS DR11 when we start with an approximation of the RSD. The scheme takes six iterations when the initial estimate, measured from the redshift-space overdensity, has no RSD correction. Slower convergence would be expected for surveys covering a larger angle on the sky. We show that this FFT based method provides a better estimate of the real space displacement field than a configuration space method that uses finite difference routines to compute the potential for the same grid resolution. Finally we show that a lognormal transform of the overdensity, used as a proxy for the linear overdensity, is beneficial in estimating the full displacement field from a dense sample of tracers. However the lognormal transform of the overdensity does not perform well when estimating the displacements from sparser simulations with a more realistic galaxy density.
The JEM-EUSO mission aims to explore the origin of the extreme energy cosmic rays (EECRs) through the observation of air-shower fluorescence light from space. The superwide-field telescope looks down from the International Space Station onto the night sky to detect UV photons (fluorescence and Cherenkov photons) emitted from air showers. Such a space detector offers the remarkable opportunity to observe a huge volume of atmosphere at once and will achieve an unprecedented statistics within a few years of operation. Several test experiments are currently in operation: e.g., one to observe the fluorescence background from the edge of the Atmosphere (EUSO-Balloon), or another to demonstrate on ground the capability of detecting air showers with a EUSO-type telescope (EUSO-TA). In this contribution a short review on the scientific objectives of the mission and an update of the instrument definition, performances and status, as well as status of the test experiments will be given.
In recent years, our understanding of gamma-ray bursts (GRB) prompt emission has been revolutionized, due to a combination of new instruments, new analysis methods and novel ideas. In this review, I describe the most recent observational results and the current theoretical interpretation. Observationally, a major development is the rise of time-resolved spectral analysis. These led to (I) identification of a distinguished high energy component, with GeV photons often seen at a delay; and (II) firm evidence for the existence of a photospheric (thermal) component in a large number of bursts. These results triggered many theoretical efforts aimed at understanding the physical conditions in the inner jet regions from which the prompt photons are emitted, as well as the spectral diversity observed. I highlight some areas of active theoretical research. These include: (I) understanding the role played by magnetic fields in shaping the dynamics of GRB outflow and spectra; (II) understanding the microphysics of kinetic and magnetic energy transfer, namely accelerating particle to high energies in both shock waves and magnetic reconnection layers; (III) understanding how sub-photospheric energy dissipation broadens the "Planck" spectrum; and (IV) geometrical light aberration effects. I highlight some of these efforts, and point towards gaps that still exist in our knowledge as well as promising directions for the future.
Over the last 15 years, the supernova community has endeavoured to identify progenitor stars of core-collapse supernovae in high resolution archival images of their galaxies.This review compiles results (from 1999 - 2013) in a distance limited sample and discusses the implications. The vast majority of the detections of progenitor stars are of type II-P, II-L or IIb with one type Ib progenitor system detected and many more upper limits for progenitors of Ibc supernovae (14). The data for these 45 supernovae progenitors illustrate a remarkable deficit of high luminosity stars above an apparent limit of Log L ~= 5.1 dex. For a typical Salpeter IMF, one would expect to have found 13 high luminosity and high mass progenitors. There is, possibly, only one object in this time and volume limited sample that is unambiguously high mass (the progenitor of SN2009ip). The possible biases due to the influence of circumstellar dust and sample selection methods are reviewed. It does not appear likely that these can explain the missing high mass progenitor stars. This review concludes that the observed populations of supernovae in the local Universe are not, on the whole, produced by high mass (M > ~18Msun) stars. Theoretical explosions of model stars also predict that black hole formation and failed supernovae tend to occur above M > ~18Msun. The models also suggest there are islands of explodability for stars in the 8-120Msun range. The observational constraints are quite consistent with the bulk of stars above M > ~18Msun collapsing to form black holes with no visible supernovae. (Abridged).
Broad absorption lines seen in some quasars prove the existence of ionized plasma outflows from the accretion disk. Outflows together with powerful jets are important feedback processes. Understanding physics behind BAL outflows might be a key to comprehend Galaxy Evolution as a whole. First radio-loud BAL quasar was discovered in 1997 and this discovery has opened new possibilities for studies of the BAL phenomena, this time on the basis of radio emission. However, information about the radio structures, orientation and age of BAL quasars is still very limited due to weak radio emission and small sizes of these objects. Our high-resolution radio survey of a sample of BAL quasars aims to increase our knowledge about these objects. In this article, we present some conclusions arising from our research.
The fundamental plane (FP) is a widely used tool to investigate the properties of early-type galaxies, and the tight relation between its parameters has spawned several cosmological applications, including its use as a distance indicator for peculiar velocity surveys and as a means to suppress intrinsic noise in cosmic size magnification measurements. Systematic trends with the large-scale structure across the FP could cause serious biases for these cosmological probes, but may also yield new insights into the early-type population. Here we report the first detection of spatial correlations among offsets in galaxy size from an FP that explicitly accounts for redshift trends, using a sample of about $95,000$ elliptical galaxies from the Sloan Digital Sky Survey. We show that these offsets correlate with the density field out to at least $10h^{-1}$Mpc at $4\sigma$ significance in a way that cannot be explained by systematic errors in galaxy size estimates. We propose a physical explanation for the correlations by dividing the sample into central, satellite, and field galaxies, identifying trends for each galaxy type separately. Central (satellite) galaxies lie on average above (below) the FP, which we argue could be due to a higher (lower) than average mass-to-light ratio. We fit a simple model to the correlations of FP residuals and use it to predict the impact on peculiar velocity power spectra, finding a contamination larger than $10\,\%$ for $k>0.04\,h/$Mpc. Moreover, cosmic magnification measurements based on an FP could be severely contaminated over a wide range of scales by the intrinsic FP correlations.
The amount of dust estimated from infrared to sub-millimetre (submm) observations strongly depends on assumptions of different grain sizes, compositions and optical properties. Here we use a simple model of thermal emission from cold silicate/carbon dust at a range of dust grain temperatures and fit the spectral energy distribution (SED) of the Crab Nebula as a test. This can lower the derived dust mass for the Crab by ~50% and 30-40% for astronomical silicates and amorphous carbon grains compared to recently published values (0.25M_sun -> 0.12M_sun and 0.12M_sun -> 0.072M_sun, respectively), but the implied dust mass can also increase by as much as almost a factor of six (0.25M_sun -> 1.14M_sun and 0.12M_sun -> 0.71M_sun) depending on assumptions regarding the sizes/temperatures of the coldest grains. The latter values are clearly unrealistic due to the expected metal budget, though. Furthermore, we show by a simple numerical experiment that if a cold-dust component does have a grain-temperature distribution, it is almost unavoidable that a two-temperature fit will yield an incorrect dust mass estimate. But we conclude that grain temperatures is not a greater uncertainty than the often poorly constrained emissivities (i.e., material properties) of cosmic dust, although there is clearly a need for improved dust emission models. The greatest complication associated with deriving dust masses still arises in the uncertainty in the dust composition.
Phase angle-induced spectral effects are important to characterize since they affect spectral band parameters such as band depth and band center, and therefore skew mineralogical interpretations of planetary bodies via reflectance spectroscopy. Dwarf planet (1) Ceres is the next target of NASA's Dawn mission, which is expected to arrive in March 2015. The visible and near-infrared mapping spectrometer (VIR) onboard Dawn has the spatial and spectral range to characterize the surface between 0.25-5.0 microns. Ceres has an absorption feature at 3.0 microns due to hydroxyl- and/or water-bearing minerals (e.g. Lebofsky et al. 1981, Rivkin et al. 2003). We analyzed phase angle-induced spectral effects on the 3-micron absorption band on Ceres using spectra measured with the long-wavelength cross-dispersed (LXD: 1.9-4.2 microns) mode of the SpeX spectrograph/imager at the NASA Infrared Telescope Facility (IRTF). Ceres LXD spectra were measured at different phase angles ranging from 0.7o to 22o. We found that the band center slightly increases from 3.06 microns at lower phase angles (0.7o and 6o) to 3.07 microns at higher phase angles (11 o and 22o), the band depth decreases by ~20% from lower phase angles to higher phase angles, and the band area decreases by ~25% from lower phase angles to higher phase angles. Our results will have implications for constraining the abundance of OH on the surface of Ceres from VIR spectral data, which will be acquired by Dawn starting spring 2015.
The radio-loud BCG at the center of the cool core cluster RBS 797 is known to exhibit a misalignment of its 5 GHz radio emission observed at different VLA resolutions, with the innermost kpc-scale jets being almost orthogonal to the radio lobes which extends for tens of kpc filling the X-ray cavities seen by Chandra. The different radio directions may be caused by rapid jet reorientation due to interaction with a secondary supermassive black hole (SMBH), or to the presence of two AGN, probably in a merging phase, which are emitting radio jets in different directions. We present the results of new 5 GHz observations performed with the EVN in May 2013. In particular, we detected two compact radio components, with a projected separation of 77 pc. We discuss two possible scenarios for the origin and nature of the EVN double source, showing that both interpretations are consistent with the presence of a SMBH binary system in the BCG of RBS 797.
We performed a detailed analysis of the main theoretical uncertainties affecting the age at the lithium depletion boundary (LDB). To do that we computed almost 12000 pre-main sequence models with mass in the range [0.06, 0.4] M_sun by varying input physics (nuclear reaction cross-sections, plasma electron screening, outer boundary conditions, equation of state, and radiative opacity), initial chemical elements abundances (total metallicity, helium and deuterium abundances, and heavy elements mixture), and convection efficiency (mixing length parameter, alpha_ML). As a first step, we studied the effect of varying these quantities individually within their extreme values. Then, we analysed the impact of simultaneously perturbing the main input/parameters without an a priori assumption of independence. Such an approach allowed us to build for the first time the cumulative error stripe, which defines the edges of the maximum uncertainty region in the theoretical LDB age. We found that the cumulative error stripe is asymmetric and dependent on the adopted mixing length value. For alpha_ML = 1.00, the positive relative age error ranges from 5 to 15 per cent, while for solar-calibrated mixing length, the uncertainty reduces to 5-10 per cent. A large fraction of such an error (about 40 per cent) is due to the uncertainty in the adopted initial chemical elements abundances.
Binary systems emit gravitational waves in a well-known pattern; for binaries in circular orbits, the emitted radiation has a frequency that is twice the orbital frequency. Systems in eccentric orbits, however, emit gravitational radiation in the higher harmonics too. In this paper, we are concerned with the stochastic background of gravitational waves generated by double neutron star systems of cosmological origin in eccentric orbits. We consider in particular the long-lived systems, that is, those binaries for which the time to coalescence is longer than the Hubble time ($\sim 10$Gyr). Thus, we consider double neutron stars with orbital frequencies ranging from $10^{-8}$ to $2\times 10^{-6}$Hz. Although in the literature some papers consider the spectra generated by eccentric binaries, there is still space for alternative approaches for the calculation of the backgrounds. In this paper, we use a method that consists in summing the spectra that would be generated by each harmonic separately in order to obtain the total background. This method allows us to clearly obtain the influence of each harmonic on the spectra. In addition, we consider different distribution functions for the eccentricities in order to investigate their effects on the background of gravitational waves generated. At last, we briefly discuss the detectability of this background by space-based gravitational wave antennas and pulsar timing arrays.
By statistically analyzing a large sample which includes blazars of Fermi detection (FBs) and non-Fermi detection (NFBs), we find that there are significant differences between FBs and NFBs for redshift, black hole mass, jet kinetic power from "cavity" power, broad-line luminosity, and ratio of core luminosity to absolute V-band magnitude ($R_{\rm v}$), but not for ratio of radio core to extended flux ($R_{\rm c}$) and Eddington ratio. Compared with NFBs, FBs have larger mean jet power, $R_{\rm c}$ and $R_{\rm v}$ while smaller mean redshift, black hole mass, broad-line luminosity. These results support that the beaming effect is main reason for differences between FBs and NFBs, and that FBs are likely to have a more powerful jet. For both Fermi and non-Fermi blazars, there are significant correlations between jet power and the accretion rate (traced by the broad-emission-lines luminosity), between jet power and black hole mass; for Fermi blazars, the black hole mass does not have significant influence on jet power while for non-Fermi blazars, both accretion rate and black hole mass have contributions to the jet power. Our results support the "blazar sequence" and show that synchrotron peak frequency ($\nu_{\rm peak}$) is associated with accretion rate but not with black hole mass.
CIZA J1358.9-4750 is a nearby (z = 0.074) pair of clusters of galaxies located close to the Galactic plane. It consists of two X-Ray extended humps at north-west and south-east separated by 14 arcmin (~ 1.2 Mpc), and an X-Ray bright bridge-like structure in between. With Suzaku, the south-east hump was shown to have a temperature of 5.6 keV and the north-west one 4.6 keV. Neither humps exhibit significant central cool component. The bridge region has a temperature higher than 9 keV at the maximum, and this hot region is distributed almost orthogonal to the bridge axis in agreement with the shock heating seen in numerical simulations at an early phase of a head-on major merger. This resemblance is supported by good positional coincidence between the X-Ray peaks and cD galaxies associated with each cluster. In a short exposure XMM-Newton image, a significant intensity jump was found at a position where the Suzaku-measured temperature exhibits a steep gradient. These properties indicate the presence of a shock discontinuity. The Mach number is estimated to be 1.32 from the temperature difference across the identified shock front, which gives the colliding velocity of approximately 1800 km/s. From optical redshifts of the member galaxies, the two clusters are indicated to be merging nearly on the sky plane. Thus, CIZA J1358.9-4750 is considered as a valuable nearby example of early-phase merger with a clear shock feature.
The primary goal of this study is to observe how complex the behaviors of
Cepheids can be, and to show how the continued monitoring of Cepheids at
multiple wavelengths can begin to reveal their "secret lives."
We aim to achieve this through optical photometry, UV spectroscopy and X-ray
imaging. Through Villanova's guaranteed access to ground-based telescopes, we
have secured well-covered light curves as regularly as possible. Amplitudes and
times of max brightness were obtained and compared to previous literature
results. At UV wavelengths, we have secured hi-res spectra of 2 nearby Cepheids
- delta Cep and beta Dor - with HST-COS. Also, we have obtained X-ray images of
5 Cepheids with XMM-Newton and the Chandra X-ray Observatory, and further
observations with both satellites have been proposed for (XMM) and approved
(Chandra).
Optical photometry has shown that 8 of the 10 observed Cepheids have
amplitude variability, or hints thereof, and all 10 show period variability
(recent, long-term or possibly periodic). The UV spectra reveal emission lines
from heated atmospheric plasmas of 10^4 - 10^5 K that vary in phase with the
Cepheid pulsations. The X-ray images have detected the three nearest Cepheids
observed (Polaris, delta Cep and beta Dor), while the distances of the other
two place their fluxes likely at or below detector background levels. The X-ray
fluxes for delta Cep show possible phased variability, but anti-correlated with
the UV emission lines (i.e. high X-ray flux during low UV flux, and vice
versa).
Further data are required to ultimately confirm Blazhko-like cycles in
Cepheids, X-ray variability with phase and the particulars of the high-energy
variability such as phase-lags between atmospheric plasma emissions of
different temperature and the exact contributions of the possible heating
mechanism.
We investigate the 2D excitation structure of the ISM in a sample of LIRGs and Seyferts using near-IR IFS. This study extends to the near-IR the well-known optical and mid-IR emission line diagnostics used to classify activity in galaxies. Based on the spatially resolved spectroscopy of prototypes, we identify in the [FeII]1.64/Br$\gamma$ - H_2 1-0S(1)/Br$\gamma$ plane regions dominated by the different heating sources, i.e. AGNs, young MS massive stars, and evolved stars i.e. supernovae. The ISM in LIRGs occupy a wide region in the near-IR diagnostic plane from -0.6 to +1.5 and from -1.2 to +0.8 (in log units) for the [FeII]/Br$\gamma$ and H_2/Br$\gamma$ line ratios, respectively. The corresponding median(mode) ratios are +0.18(0.16) and +0.02(-0.04). Seyferts show on average larger values by factors ~2.5 and ~1.4 for the [FeII]/Br$\gamma$ and H_2/Br$\gamma$ ratios, respectively. New areas and relations in the near-IR diagnostic plane are defined for the compact, high surface brightness regions dominated by AGN, young ionizing stars, and SNe explosions, respectively. In addition, the diffuse regions affected by the AGN radiation field cover an area similar to that of Seyferts, but with high values in [FeII]/Br$\gamma$ that are not as extreme. The extended, non-AGN diffuse regions cover a wide area in the diagnostic diagram that overlaps that of individual excitation mechanisms (i.e. AGN, young stars, and SNe), but with its mode value to that of the young SF clumps. This indicates that the excitation conditions of the diffuse ISM are likely due to a mixture of the different ionization sources. The integrated line ratios in LIRGs show higher excitation conditions i.e. towards AGNs, than those measured by the spatially resolved spectroscopy. If this behaviour is representative, it would have clear consequences when classifying high-z, SF galaxies based on their near-IR integrated spectra.
Aims: We present a new upgraded version of ARES. The new version includes a series of interesting new features such as automatic radial velocity correction, a fully automatic continuum determination, and an estimation of the errors for the equivalent widths. Methods: The automatic correction of the radial velocity is achieved with a simple cross-correlation function, and the automatic continuum determination, as well as the estimation of the errors, relies on a new approach to evaluating the spectral noise at the continuum level. Results: ARES v2 is totally compatible with its predecessor. We show that the fully automatic continuum determination is consistent with the previous methods applied for this task. It also presents a significant improvement on its performance thanks to the implementation of a parallel computation using the OpenMP library.
Raman scattering by H$_2$ in Neptune's atmosphere has significant effects on its reflectivity for $\lambda <$ 0.5 $\mu$m, producing baseline decreases of $\sim$ 20% in a clear atmosphere and $\sim$ 10% in a hazy atmosphere. Here we present the first radiation transfer algorithm that includes both polarization and Raman scattering and facilitates computation of spatially resolved spectra. New calculations show that Cochran and Trafton's (1978, Astrophys. J. 219, 756-762) suggestion that light reflected in the deep CH$_4$ bands is mainly Raman scattered is not valid for current estimates of the CH$_4$vertical distribution, which implies only a 4% Raman contribution. Comparisons with IUE, HST, and groundbased observations confirm that high altitude haze absorption is reducing Neptune's geometric albedo by $\sim$6% in the 0.22-0.26 $\mu$m range and by $\sim$13% in the 0.35-0.45 $\mu$m range. We used accurate calculations to evaluate several approximations of Raman scattering. The Karkoschka (1994, Icarus 111, 174-192) method of removing Raman effects from observed spectra is shown to have limited applicability and to undercorrect the depths of weak CH$_4$ absorption bands. The Wallace (1972, Astrophys. J. 176, 249-257) approximation produces geometric albedo values $\sim$5% low as originally proposed, but can be much improved by adding scattering contributions from the vibrational transition. The Pollack et al. (1986, Icarus 65, 442-466) approximation is inaccurate and unstable, but can also be improved greatly by several simple modifications. A new approximation provides low errors for zenith angles below 70\deg in a clear atmosphere, although intermediate clouds present problems at longer wavelengths.
We have detected four new HgMn stars, while monitoring a sample of apparently slowly rotating superficially normal bright late B and early A stars in the northern hemisphere. Important classification lines of Hg II and Mn II are found as conspicuous features in the high resolution SOPHIE spectra of these stars (R = 75000). Several lines of Hg II, Mn II and Fe II have been synthesized using model atmospheres and the spectrum synthesis code SYNSPEC48 including hyperfine structure of various isotopes when relevant. These synthetic spectra have been compared to high resolution high signal-to-noise observations of these stars in order to derive abundances of these key elements. The four stars are found to have distinct enhancements of Hg and Mn which show that these stars are not superficially normal B and A stars, but actually are new HgMn stars and should reclassified as such.
Planets are thought to form in the gas and dust disks around young stars. In particular, it has been proposed that giant planets can form via gravitational instability of massive extended disks around intermediate mass stars. However, direct observations to constrain this mechanism lack. We have spatially resolved the 8.6 and 11.2 $\mu$m emission of a massive edge on protoplanetary disk around an A star, Gomez's Hamburger (GoHam), using VISIR at the Very Large Telescope. A compact region situated at a projected distance of $350\pm50$ AU South of the central star is found to have a reduced emission. This asymmetry is fully consistent with the presence of a cold density structure, or clump, identified in earlier CO observations, and we derive physical characteristics consistent with those observations: a mass of 0.8-11.4 Jupiter masses (for a dust to gas mass ratio of 0.01), a radius of the order of 10$^2$ astronomical units, a local density of the order of $10^{7}$ cm$^{-3}$. Based on this evidence, we argue that this clump, which we call GoHam b, is a promising candidate for a young protoplanet formed by gravitational instability, that could be representative of the precursors of massive planets observed around A stars, like HR 8799. Further studies at high angular resolution are needed to better constrain the physical properties of this object in order to confirm this proposal.
Ultraviolet and 21-cm observations suggest that the extremely low-metallicity galaxy, I Zw 18, is a stream-fed galaxy containing a "pocket" of pristine stars responsible for producing nebular He II recombination emission observed in I Zw18-NW. Far-UV spectra by Hubble/COS and the Far Ultraviolet Spectroscopic Explorer (FUSE) make this suggestion conclusive by demonstrating that the spectrum of I Zw 18-NW shows no metal lines like O VI 1032, 1038 of comparable ionization as the He II recombination emission.
The emission spectrum has been calculated of a homogeneous pure hydrogen
layer, which parameters are typical for a flare on a red dwarf. The ionization
and excitation states were determined by the solution of steady-state equations
taking into account the continuum and all discrete hydrogen levels. We consider
the following elementary processes: electron-impact transitions, spontaneous
and induced radiative transitions, and ionization by the bremsstrahlung and
recombination radiation of the layer itself. The Biberman--Holstein
approximation was used to calculate the scattering of line radiation.
Asymptotic formulae for the escape probability are obtained for a symmetric
line profile taking into account the Stark and Doppler effects. The
approximation for the core of the H$-\alpha$ line by a gaussian curve has been
substantiated.
The spectral intensity of the continuous spectrum, the intensity of the lines
of the Balmer series and the magnitude of the Balmer jump have been calculated.
The conditions have been determined for which the Balmer jump and the emission
line intensities above the continuum decrease to such low values that the
emission spectrum can be assumed to be continuum as well as the conditions at
which the emission spectrum becomes close to the blackbody.
Separating active regions that are quiet from potentially eruptive ones is a key issue in Space Weather applications. Traditional classification schemes such as Mount Wilson and McIntosh have been effective in relating an active region large scale magnetic configuration to its ability to produce eruptive events. However, their qualitative nature prevents systematic studies of an active region's evolution for example. We introduce a new clustering of active regions that is based on the local geometry observed in Line of Sight magnetogram and continuum images. We use a reduced-dimension representation of an active region that is obtained by factoring (i.e. applying dictionary learning to) the corresponding data matrix comprised of local image patches. Two factorizations can be compared via the definition of appropriate metrics on the resulting factors. The distances obtained from these metrics are then used to cluster the active regions. We find that these metrics result in natural clusterings of active regions. The clusterings are related to large scale descriptors of an active region such as its size, its local magnetic field distribution, and its complexity as measured by the Mount Wilson classification scheme. We also find that including data focused on the neutral line of an active region can result in an increased correspondence between the Mount Wilson classifications and our clustering results. We provide some recommendations for which metrics and matrix factorization techniques to use to study small, large, complex, or simple active regions.
Context. ABDoradus is the main system of the ABDoradus moving group. It is a quadruple system formed by two widely separated binaries of pre-main-sequence (PMS) stars: ABDor A/C and ABDor Ba/Bb. The pair ABDor A/C has been extensively studied and its dynamical masses have been determined with high precision, thus making of ABDor C a benchmark for calibrating PMS stellar models. If the orbit and dynamical masses of the pair ABDor Ba/Bb can be determined, they could not only play a similar role to that of ABDor C in calibrating PMS models, but would also help to better understand the dynamics of the whole ABDoradus system. Aims. We aim to determine the individual masses of the pair ABDor Ba/Bb using VLBI observations and archive infrared data, as part of a larger program directed to monitor binary systems in the ABDoradus moving group. Methods. We observed the system ABDor B between 2007 and 2013 with the Australian Long Baseline Array (LBA), at a frequency of 8.4 GHz in phase-reference mode. Results. We detected, for the first time, compact radio emission from both stars in the binary, ABDor Ba and ABDor Bb. This result allowed us to determine the orbital parameters of both the relative and absolute orbits and, consequently, their individual dynamical masses: 0.28+/-0.05 Msun and 0.25+/-0.05 Msun, respectively. Conclusions. Comparisons of the dynamical masses with the prediction of PMS evolutionary models show that the models underpredict the dynamical masses of the binary components Ba and Bb by 10-30% and 10-40%, respectively, although they all still agree at the 2-sigma level. Some of stellar models considered favour an age between 50 and 100 Myr for this system, meanwhile others predict older ages. We also discuss the evolutionary status of ABDor Ba/Bb in terms of an earlier double-double star scenario that might explain the strong radio emission detected in both components.
Flat-spectrum radio-loud Narrow-Line Seyfert 1 galaxies (NLS1s) are a recently discovered class of $\gamma$-ray emitting Active Galactic Nuclei (AGN), that exhibit some blazar-like properties which are explained with the presence of a relativistic jet viewed at small angles. When blazars are observed at larger angles they appear as radio-galaxies, and we expect to observe an analogue parent population for beamed NLS1s. However, the number of known NLS1s with the jet viewed at large angles is not enough. Therefore, we tried to understand the origin of this deficit. Current hypotheses about the nature of parent sources are steep-spectrum radio-loud NLS1s, radio-quiet NLS1s and disk-hosted radio-galaxies. To test these hypotheses we built three samples of candidate sources plus a control sample, and calculated their black hole mass and Eddington ratio using their optical spectra. We then performed a Kolmogorov-Smirnov statistical test to investigate the compatibility of our different samples with a beamed population. Our results indicate that, when the inclination angle increases, a beamed source appears as a steep-spectrum radio-loud NLS1, or possibly even as a disk-hosted radio-galaxy with low black hole mass and high Eddington ratio. Further investigations, involving larger complete samples and observations at radio frequency, are needed to understand the incidence of disk-hosted radio-galaxies in the parent population, and to assess whether radio-quiet NLS1s can play a role, as well.
I discuss an upper bound on the boost and the energy of elementary particles. The limit is derived utilizing the core principle of relativistic quantum mechanics stating that there is a lower limit for localization of an elementary quantum system and the suggestion that when the localization scale reaches the Planck length, elementary particles are removed from observables. The limit for the boost and energy, $M_{Planck}/m$ and $M_{Planck}c^{2}\approx\,8.6* 10^{27}$ eV, is defined in terms of fundamental constants and the mass of elementary particle and does not involve any dynamic scale. These bounds imply that the cosmic ray flux of any flavor may stretch up to energies of order $10^{18}$ GeV and will cut off at this value.
Ekpyrotic instantons describe the emergence of classical contracting universes out of the no-boundary quantum state. However, up to now these instantons ended in a big crunch singularity. We remedy this by adding a higher-derivative term, allowing a ghost condensate to form. This causes a smooth, non-singular bounce from the contracting phase into an expanding, kinetic-dominated phase. Remarkably, and although there is a non-trivial evolution during the bounce, the wavefunction of the universe is "classical" in a WKB sense just as much after the bounce as before. These new non-singular instantons can thus form the basis for a fully non-singular and calculable ekpyrotic history of the universe, from creation until now.
We have explored Natural Supersymmetry (NSUSY) scenarios with low values of the $\mu$ parameter which are characterised by higgsino-like Dark Matter (DM) and compressed spectra for the lightest MSSM particles, $\chi^0_1$, $\chi^0_2$ and $\chi^\pm_1$. This scenario could be probed via monojet signatures, but as the signal-to-background ratio (S/B) is low we demonstrate that the 8 TeV LHC cannot obtain limits on the DM mass beyond those of LEP2. On the other hand, we have found, for the 13 TeV run of the LHC, that by optimising kinematical cuts we can bring the S/B ratio up to the 5(3)% level which would allow the exclusion of the DM mass up to 200(250) GeV respectively, significantly extending LEP2 limits. Moreover, we have found that LUX/XENON1T and LHC do play very complementary roles in exploring the parameter space of NSUSY, as the LHC has the capability to access regions where DM is quasi-degenerate with other higgsinos, which are challenging for direct detection experiments.
No-hair like relations between the multipole moments of the exterior gravitational field of neutron stars have recently been found to be approximately independent of the star's internal structure. This approximate, equation-of-state universality arises after one adimensionalizes the multipole moments appropriately, which then begs the question of whether there are better ways to adimensionalize the moments to obtain stronger universality. We here investigate this question in detail by considering slowly-rotating neutron stars both in the non-relativistic limit and in full General Relativity. We find that there exist normalizations that lead to stronger equation-of-state universality in the relations among the moment of inertia and the quadrupole, octopole and hexadecapole moments of neutron stars. We determine the optimal normalization that minimizes the equation-of-state dependence in these relations. The results found here may have applications in the modeling of X-ray pulses and atomic line profiles from millisecond pulsars with NICER and LOFT.
Because of physical processes ranging from microscopic particle collisions to macroscopic hydrodynamic fluctuations, any plasma in thermal equilibrium emits gravitational waves. For the largest wavelengths the emission rate is proportional to the shear viscosity of the plasma. In the Standard Model at T > 160 GeV, the shear viscosity is dominated by the most weakly interacting particles, right-handed leptons, and is relatively large. We estimate the order of magnitude of the corresponding spectrum of gravitational waves. Even though at small frequencies (corresponding to the sub-Hz range relevant for planned observatories such as eLISA) this background is tiny compared with that from non-equilibrium sources, the total energy carried by the high-frequency part of the spectrum is non-negligible if the production continues for a long time. We suggest that this may constrain (weakly) the highest temperature of the radiation epoch. Observing the high-frequency part directly sets a very ambitious goal for future generations of GHz-range detectors.
We consider numerical models of Boussinesq convection in non-rotating spherical shells for a fluid with a unity Prandtl number and Rayleigh numbers up to $10^9$. In this geometry, curvature and radial variations of the gravitationnal acceleration yield asymmetric boundary layers. A systematic parameter study for various radius ratios (from $\eta=r_i/r_o=0.2$ to $\eta=0.95$) and gravity profiles allows us to explore the dependence of the asymmetry on these parameters. We find that the average plume spacing is conserved between the spherical inner and outer bounding surfaces. An estimate of the average plume separation allows us to accurately predict the boundary layer asymmetry for the various spherical shell configurations explored here. The mean temperature and horizontal velocity profiles are in good agreement with classical Prandtl-Blasius laminar boundary layer profiles, provided the boundary layers are analysed in a dynamical frame, that fluctuates with the local and instantaneous boundary layer thicknesses. The scaling properties of the Nusselt and Reynolds numbers are investigated by separating the different contributions to the thermal and viscous dissipation rates using numerical models with $\eta=0.6$ and a gravity proportional to $1/r^2$. We show that our spherical models are consistent with the predictions of Grossmann \& Lohse's (2000) theory and that $Nu(Ra)$ and $Re(Ra)$ scalings are in good agreement with plane layer results.
Based on the analogy with superconductor physics we consider a scalar-vector-tensor gravitational model, in which the dark energy action is described by a gauge invariant electromagnetic type functional. By assuming that the ground state of the dark energy is in a form of a condensate with the U(1) symmetry spontaneously broken, the gauge invariant electromagnetic dark energy can be described in terms of the combination of a vector and of a scalar field (corresponding to the Goldstone boson), respectively. The gravitational field equations are obtained by also assuming the possibility of a non-minimal coupling between the cosmological mass current and the superconducting dark energy. The cosmological implications of the dark energy model are investigated for a Friedmann-Robertson-Walker homogeneous and isotropic geometry for two particular choices of the electromagnetic type potential, corresponding to a pure electric type field, and to a pure magnetic field, respectively. The time evolution of the scale factor, matter energy density and deceleration parameter are obtained for both cases, and it is shown that in the presence of the superconducting dark energy the Universe ends its evolution in an exponentially accelerating vacuum de Sitter state. By using the formalism of the irreversible thermodynamic processes for open systems we interpret the generalized conservation equations in the superconducting dark energy model as describing matter creation. The particle production rates, the creation pressure and the entropy evolution are explicitly obtained.
Matrix determinants play an important role in data analysis, in particular when Gaussian processes are involved. Due to currently exploding data volumes linear operations - matrices - acting on the data are often not accessible directly, but are only represented indirectly in form of a computer routine. Such a routine implements the transformation a data vector undergoes under matrix multiplication. Meanwhile efficient probing routines to estimate a matrix's diagonal or trace, based solely on such computationally affordable matrix-vector multiplications, are well known and frequently used in signal inference, a stochastic estimate for its determinant is still lacking. In this work a probing method for the logarithm of a determinant of a linear operator is introduced. This method rests upon a reformulation of the log-determinant by an integral representation and the transformation of the involved terms into stochastic expressions. This stochastic determinant determination enables large-size applications in Bayesian inference, in particular evidence calculations, model comparison, and posterior determination.
In this communication, we consider a wide class of extensions to General Relativity that break explicitly the Einstein Equivalence Principle by introducing a multiplicative coupling between a scalar field and the electromagnetic Lagrangian. In these theories, we show that 4 cosmological observables are intimately related to each other: a temporal variation of the fine structure constant, a violation of the distance-duality relation, the evolution of the cosmic microwave background (CMB) temperature and CMB spectral distortions. This enables one to put very stringent constraints on possible violations of the distance-duality relation, on the evolution of the CMB temperature and on admissible CMB spectral distortions using current constraints on the fine structure constant. Alternatively, this offers interesting possibilities to test a wide range of theories of gravity by analyzing several data sets concurrently.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
Using the spectroscopic catalogue of the Sloan Digital Survey Data Release 10 (SDSS DR10), we have explored the abundance of satellites around a sample of 307 massive (M_star > 10^11 M_sun) local (z < 0.025) galaxies. We have divided our sample into 4 morphological groups (E, S0, Sa, Sb/c). We find that the number of satellites with M_star > 10^9 M_sun and R < 300 kpc depends drastically on the morphology of the central galaxy. The average number of satellites per galaxy host (N_Sat/N_Host) down to a mass ratio of 1:100 is: 5.5 +/- 1.0 for E hosts, 2.7 +/- 0.4 for S0, 1.4 +/- 0.3 for Sa and 1.2 +/- 0.3 for Sb/c. The amount of stellar mass enclosed by the satellites around massive E-type galaxies is a factor of 2, 4, and 6 larger than the mass in the satellites of S0, Sa and Sb/c-types, respectively. If these satellites would eventually infall into the host galaxies, for all the morphological types, the merger channel will be largely dominated by satellites with a mass ratio satellite-host $\mu$ < 0.1. The fact that massive elliptical galaxies have a significant larger number of satellites than massive spirals could point out that elliptical galaxies inhabit heavier dark matter halos than equally massive galaxies with later morphological types. If this hypothesis is correct, the dark matter halos of late-type spiral galaxies are a factor ~3 more efficient on producing galaxies with the same stellar mass than those dark matter halos of massive ellipticals.
Weak lensing by large-scale structure is a powerful technique to probe the dark components of the universe. To understand the measurement process of weak lensing and the associated systematic effects, image simulations are becoming increasingly important. For this purpose we present a first implementation of the $\textit{Monte Carlo Control Loops}$ ($\textit{MCCL}$; Refregier & Amara 2014), a coherent framework for studying systematic effects in weak lensing. It allows us to model and calibrate the shear measurement process using image simulations from the Ultra Fast Image Generator (UFig; Berge et al. 2013). We apply this framework to a subset of the data taken during the Science Verification period (SV) of the Dark Energy Survey (DES). We calibrate the UFig simulations to be statistically consistent with DES images. We then perform tolerance analyses by perturbing the simulation parameters and study their impact on the shear measurement at the one-point level. This allows us to determine the relative importance of different input parameters to the simulations. For spatially constant systematic errors and six simulation parameters, the calibration of the simulation reaches the weak lensing precision needed for the DES SV survey area. Furthermore, we find a sensitivity of the shear measurement to the intrinsic ellipticity distribution, and an interplay between the magnitude-size and the pixel value diagnostics in constraining the noise model. This work is the first application of the $\textit{MCCL}$ framework to data and shows how it can be used to methodically study the impact of systematics on the cosmic shear measurement.
The current standard model of cosmology, $\Lambda$CDM, requires dark matter to make up around 25% of the total energy budget of the universe. Yet, quite puzzlingly, there appears to be no candidate particle in the current Standard Model of particle physics. Assuming the validity of the CDM paradigm, dark matter has evaded detection thus far either because it is intrinsically a weakly-interacting substance or because its interactions are suppressed by its high constituent mass and low number density. Most approaches to explain dark matter to date assume the former and therefore require beyond-the-Standard-Model particles that have yet to be observed directly or indirectly. Given the dearth of evidence for this class of candidates it is timely to consider the latter possibility, which allows for candidates that may or may not arise from the Standard Model. In this work we extend a recent study of this general class of so-called macro dark matter -- candidates with characteristic masses of grams and geometric cross sections of cm$^2$. We consider new bounds that can be set using existing data from the resonant bar gravitational wave detectors NAUTILUS and EXPLORER.
We present results on the dust attenuation curve of z~2 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Our sample consists of 224 star-forming galaxies with nebular spectroscopic redshifts in the range z= 1.36-2.59 and high S/N measurements of, or upper limits on, the H-alpha and H-beta emission lines obtained with Keck/MOSFIRE. We construct composite SEDs of galaxies in bins of specific SFR and Balmer optical depth in order to directly constrain the dust attenuation curve from the UV through near-IR for typical star-forming galaxies at high redshift. Our results imply an attenuation curve that is very similar to the SMC extinction curve at wavelengths redward of 2500 Angstroms. At shorter wavelengths, the shape of the curve is identical to that of the Calzetti relation, but with a lower normalization (R_V). Hence, the new attenuation curve results in SFRs that are ~20% lower, and log stellar masses that are 0.16 dex lower, than those obtained with the Calzetti attenuation curve. Moreover, we find that the difference in the reddening---and the total attenuation---of the ionized gas and stellar continuum correlates strongly with SFR, such that for dust-corrected SFRs larger than 20 Msun/yr assuming a Chabrier IMF, the nebular emission lines suffer an increasing degree of obscuration relative to the continuum. A simple model that can account for these trends is one in which the UV through optical stellar continuum is dominated by a population of less reddened stars, while the nebular line and bolometric luminosities become increasingly dominated by dustier stellar populations for galaxies with large SFRs, as a result of the increased dust enrichment that accompanies such galaxies. Consequently, UV- and SED-based SFRs may underestimate the total SFR at even modest levels of ~20 Msun/yr. [Abridged]
In Aligned Natural Inflation, an alignment between different potential terms produces an inflaton excursion greater than the axion scales in the potential. We show that, starting from a general potential of two axions with two aligned potential terms, the effective theory for the resulting light direction is characterized by four parameters: an effective potential scale, an effective axion constant, and two extra parameters (related to ratios of the axion scales and the potential scales in the $2-$field theory). For all choices of these extra parameters, the model can support inflation along valleys (in the $2-$field space) that end in minima of the potential. This leads to a phenomenology similar to that of single field Natural Inflation. For a significant range of the extra two parameters, the model possesses also higher altitude inflationary trajectories passing through saddle points of the $2-$field potential, and disconnected from any minimum. These plateaus end when the heavier direction becomes unstable, and therefore all of inflation takes place close to the saddle point, where - due to the higher altitude - the potential is flatter (smaller $\epsilon$ parameter). As a consequence, a tensor-to-scalar ratio $r = {\rm O } \left( 10^{-4} - 10^{-2} \right)$ can be easily achieved in the allowed $n_s$ region, well within the latest $1 \sigma$ CMB contours.
A star on a nearly radial trajectory approaching a massive black hole (MBH) gets tidally disrupted if it comes sufficiently close to the MBH. Here we explore what happens to binary stars whose centers of mass approach the MBH on nearly radial orbits. The interaction with the MBH often leads to both stars being disrupted in sequence. We argue that such events could produce light curves that are substantially different from those of the single disruptions, with possible features such as two local maxima. Tidal forces from the MBH can also lead the binary components to collide; these merger products can form highly magnetized stars, whose subsequent tidal disruption may enable prompt jet formation.
In this paper we develop two models for the steady states and evolution of two dimensional isothermal self gravitating and rotating incompressible gas which are based on the hydrodynamic equations for stratified fluid. The first model is for the steady states of the gas while the second addresses the time evolution of the gas subject to some constraints. These models reduce the initial five partial differential equations that govern this system to two for the steady state model and to three for the time dependent model. Analytical and numerical solutions of the model equations are used to study the structure of the resulting steady and time dependent states of the fluid with some possible astrophysical applications.
Blazars are the dominant type of extragalactic sources at microwave and at $\gamma$-ray energies. In the most energetic part of the electromagnetic spectrum (E>100GeV) a large fraction of high Galactic latitude sources are blazars of the High Synchrotron Peaked (HSP) type, that is BL Lac objects with synchrotron power peaking in the UV or in the X-ray band. HSP blazars are remarkably rare, with only a few hundreds of them expected to be above the sensitivity limits of currently available surveys. To find these very uncommon objects, we have devised a method that combines ALLWISE survey data with multi-frequency selection criteria. The sample was defined starting from a primary list of infrared colour-colour selected sources from the ALLWISE all sky survey database, and applying further restrictions on IR-radio and IR-X-ray flux ratios. Using a polynomial fit to the multi-frequency data (radio to X-ray) we estimated synchrotron peak frequencies and fluxes of each object. We assembled a sample including 992 sources, which is currently the largest existing list of confirmed and candidates HSP blazars. All objects are expected to radiate up to the highest $\gamma$-ray photon energies. In fact, 299 of these are confirmed emitters of GeV $\gamma$-ray photons (based on Fermi-LAT catalogues), and 36 have already been detected in the TeV band. The majority of sources in the sample are within reach of the upcoming Cherenkov Telescope Array (CTA), and many may be detectable even by the current generation of Cherenkov telescopes during flaring episodes. The sample includes 425 previously known blazars, 151 new identifications, and 416 HSP candidates (mostly faint sources) for which no optical spectra is available yet. The full 1WHSP catalogue is on-line at this http URL providing a direct link to the SED building tool where multifrequency data can be easily visualised.
Accurate rotational-vibrational line lists for $^{31}$P$^{14}$N and $^{31}$P$^{15}$N in their ground electronic states are computed. The line lists are produced using an empirical potential energy curve obtained by fitting to the experimental transition frequencies available in the literature in conjunction with an accurate, high level \textit{ab initio} dipole moment curve. In these calculations the programs DPotFit and LEVEL~8.0 were used. The new line lists reproduce the experimental wavenumbers with a root-mean-square error of 0.004~cm$^{-1}$. The line lists cover the frequency range 0--51000 cm$^{-1}$, contain almost 700~000 lines each and extend up to a maximum vibrational level of $v$=66 and a maximum rotational level of $J$=357. They should be applicable for a large range of temperature up to, at least, 5000~K. These new line lists are used to simulate spectra for PN at a range of temperatures and are deposited in the Strasbourg data centre. This work is performed as part of the ExoMol project.
We investigated the reliability of the genetic algorithm which will be used to invert the photometric measurements of asteroids collected by the European Space Agency Gaia mission. To do that, we performed several sets of simulations for 10 000 asteroids having different spin axis orientations, rotational periods and shapes. The observational epochs used for each simulation were extracted from the Gaia mission simulator developed at the Observatoire de la C\^{o}te d'Azur, while the brightness was generated using a Z-buffer standard graphic method. We also explored the influence on the inversion results of contaminating the data set with Gaussian noise with different $\sigma$ values. The research enabled us to determine a correlation between the reliability of the inversion method and the asteroid's pole latitude. In particular, the results are biased for asteroids having quasi-spherical shapes and low pole latitudes. This effect is caused by the low lightcurve amplitude observed under such circumstances, as the periodic signal can be lost in the photometric random noise when both values are comparable, causing the inversion to fail. Such bias might be taken into account when analysing the inversion results, not to mislead it with physical effects such as non-gravitational forces. Finally, we studied what impact on the inversion results has combining a full lightcurve and Gaia photometry collected simultaneously. Using this procedure we have shown that it is possible to reduce the number of wrong solutions for asteroids having less than 50 data points. The latter will be of special importance for planning ground-based observations of asteroids aiming to enhance the scientific impact of Gaia on Solar system science.
Massive young star clusters contain dozens or hundreds of massive stars that inject mechanical energy in the form of winds and supernova explosions, producing an outflow which expands into their surrounding medium, shocking it and forming structures called superbubbles. The regions of shocked material can have temperatures in excess of 10$^6$ K, and emit mainly in thermal X-rays (soft and hard). This X-ray emission is strongly affected by the action of thermal conduction, as well as by the metallicity of the material injected by the massive stars. We present three-dimensional numerical simulations exploring these two effects, metallicity of the stellar winds and supernova explosions, as well as thermal conduction.
We performed the deepest search for an X-ray emission line between 0.5 and 7 keV from non-baryonic dark matter with the Suzaku XIS. Dark matter associated with the Milky Way galaxy was selected as the target to obtain the best signal-to-noise ratio. From the Suzaku archive, we selected 187 data sets of blank sky regions which were dominated by the X-ray diffuse background. The data sets were from 2005 to 2013. Instrumental responses were adjusted by multiple calibration data sets of the Crab Nebula. We also improved the technique of subtracting lines of instrumental origin. These energy spectra were well described by X-ray emission due to charge exchange around the Solar System, hot plasma in and around the Milky Way and superposition of extra-galactic point sources. A signal of a narrow emission line was searched for, and the significance of detection was evaluated in consideration of the blind search method (the Look-elsewhere Effect). Our results exhibited no significant detection of an emission line feature from dark matter. The 3$\sigma$ upper limit for the emission line intensity between 1 and 7 keV was $\sim10^{-2}$ photons cm$^{-2}$ s$^{-1}$ sr$^{-1}$, or $\sim 5\times10^{-4}$ photons cm$^{-2}$ s$^{-1}$ sr$^{-1}$ per $M_\odot$ pc$^{-2}$, assuming a dark matter distribution with the Galactic rotation curve. The parameters of sterile neutrinos as candidates of dark matter were also constrained.
We present Magellan/M2FS, VLT/GIRAFFE, and Gemini South/GMOS spectroscopy of the newly discovered Milky Way satellite Reticulum II. Based on the spectra of 25 Ret II member stars selected from Dark Energy Survey imaging, we measure a mean heliocentric velocity of 62.8 +/- 0.5 km/s and a velocity dispersion of 3.3 +/- 0.7 km/s. The mass-to-light ratio of Ret II within its half-light radius is 470 +/- 210 Msun/Lsun, demonstrating that it is a strongly dark matter-dominated system. Despite its spatial proximity to the Magellanic Clouds, the radial velocity of Ret II differs from that of the LMC and SMC by 199 and 83 km/s, respectively, suggesting that it is not gravitationally bound to the Magellanic system. The likely member stars of Ret II span 1.3 dex in metallicity, with a dispersion of 0.28 +/- 0.09 dex, and we identify several extremely metal-poor stars with [Fe/H] < -3. In combination with its luminosity, size, and ellipticity, these results confirm that Ret II is an ultra-faint dwarf galaxy. With a mean metallicity of [Fe/H] = -2.65 +/- 0.07, Ret II matches Segue~1 as the most metal-poor galaxy known. Although Ret II is the third-closest dwarf galaxy to the Milky Way, the line-of-sight integral of the dark matter density squared is log J = 18.8 +/- 0.6 Gev^2/cm^5 within 0.2 degrees, indicating that the predicted gamma-ray flux from dark matter annihilation in Ret II is lower than that of several other dwarf galaxies.
The Dark Energy Camera is a new imager with a 2.2-degree diameter field of view mounted at the prime focus of the Victor M. Blanco 4-meter telescope on Cerro Tololo near La Serena, Chile. The camera was designed and constructed by the Dark Energy Survey Collaboration, and meets or exceeds the stringent requirements designed for the wide-field and supernova surveys for which the collaboration uses it. The camera consists of a five element optical corrector, seven filters, a shutter with a 60 cm aperture, and a CCD focal plane of 250 micron thick fully-depleted CCDs cooled inside a vacuum Dewar. The 570 Mpixel focal plane comprises 62 2kx4k CCDs for imaging and 12 2kx2k CCDs for guiding and focus. The CCDs have 15 microns x15 microns pixels with a plate scale of 0.263 arc sec per pixel. A hexapod system provides state-of-the-art focus and alignment capability. The camera is read out in 20 seconds with 6-9 electrons readout noise. This paper provides a technical description of the camera's engineering, construction, installation, and current status.
The standard model of eruptive, two-ribbon flares involves reconnection of over- lying magnetic fields beneath a rising ejection. Numerous observers have reported evidence linking this reconnection, indicated by photospheric flux swept out by flare ribbons, to coronal mass ejection (CME) acceleration. This acceleration might be caused by reconnected fields that wrap around the ejection producing an increased outward hoop force. Other observations have linked stronger over- lying fields, measured by the power-law index of the fitted decay rate of field strength overlying eruption sites, to slower CME speeds. This might be caused by greater downward magnetic tension in stronger overlying fields. So overlying fields might both help and hinder the acceleration of CMEs: reconnection that converts overlying fields into flux winding about the ejection might help, but unreconnected overlying fields might hurt. Here, we investigate the roles of both ribbon fluxes and the decay rates of overlying fields in a set of 16 eruptive events. We confirm previous results that higher CME speeds are associated with both larger ribbon fluxes and more rapidly decaying overlying fields. We find the association with ribbon fluxes to be weaker than a previous report, but stronger than the dependence on the decay rate of overlying fields.
The equation of state of cold baryonic matter is studied within a relativistic mean-field model with hadron masses and coupling constants depending on the scalar field. All hadron masses undergo a universal scaling, whereas the coupling constants are scaled differently. The appearance of hyperons in dense neutron star interiors is accounted for, however the equation of state remains sufficiently stiff if a reduction of the $\phi$ meson mass is included. Our equation of state matches well the constraints known from analyses of the astrophysical data and the particle production in heavy-ion collisions.
Recent observations have shown that circumstellar and circumbinary discs in young stellar binaries are often misaligned with respect to the binary orbital plane. We analyze the tidal truncation of such misaligned discs due to torques applied to the disc at the Lindblad resonances from the tidal forcings of the binary. We consider eccentric binaries with arbitrary binary-disc inclination angles. We determine the dependence of the tidal forcing strengths on the binary parameters and show that they are complicated non-monotonic functions of eccentricity and inclination. We adopt a truncation criterion determined by the balance between resonant torque and viscous torque, and use it to calculate the outer radii of circumstellar discs and the inner radii of circumbinary discs. Misaligned circumstellar discs have systematically larger outer radii than aligned discs, and are likely to fill their Roche lobes if inclined by more than $45^\circ - 90^\circ$, depending on the binary mass ratio and disc viscosity parameter. Misaligned circumbinary discs generally have smaller inner radii than aligned discs, but the details depend sensitively on the binary and disc parameters.
We investigate whether Gaia can specify the binary fractions of massive stellar populations in the Galactic disk through astrometric microlensing. Furthermore, we study if some information about their mass distributions can be inferred via this method. In this regard, we simulate the binary astrometric microlensing events due to massive stellar populations according to the Gaia observing strategy by considering (a) stellar-mass black holes, (b) neutron stars, (c) white dwarfs and (d) main-sequence stars as microlenses. The Gaia efficiency for detecting the binary signatures in binary astrometric microlensing events is $\sim 10-20$ per cent. By calculating the optical depth due to the mentioned stellar populations, the number of the binary astrometric microlensing events being observed with Gaia with detectable binary signatures, for the binary fraction about 0.1, is estimated as 6, 11, 77 and 1316 respectively. Consequently, Gaia can potentially specify the binary fractions of these massive stellar populations. However, the binary fraction of black holes measured with this method has the large uncertainty owing to a low number of the estimated events. Knowing the binary fractions in massive stellar populations helps for studying the gravitational waves. Moreover, we investigate the number of massive microlenses which Gaia specifies their masses through astrometric microlensing of single lenses toward the Galactic bulge. The resulted efficiencies of measuring the mass of mentioned populations are 9.8, 2.9, 1.2 and 0.8 per cent respectively. The number of their astrometric microlensing events being observed in the Gaia era in which the lens mass can be inferred with the relative error less than 0.5 toward the Galactic bulge is estimated as 45, 34, 76 and 786 respectively.
In this second paper of a series, we present a treatment procedure for data cubes obtained with the Spectrograph for Integral Field Observations in the Near Infrared of the Very Large Telescope. We verified that the treatment procedure improves significantly the quality of the images of the data cubes, allowing a more detailed analysis. The images of the Br$\gamma$ and H$_2 \lambda 21218$ emission lines from the treated data cube of the nuclear region of NGC 5643 reveal the existence of ionized and molecular-gas clouds around the nucleus, which cannot be seen clearly in the images from the non-treated data cube of this galaxy. The ionized-gas clouds represent the narrow-line region, in the form of a bicone. We observe a good correspondence between the positions of the ionized-gas clouds in the Br$\gamma$ image and in an [O III] image, obtained with the Hubble Space Telescope, of the nuclear region of this galaxy convolved with an estimate of the point-spread function of the data cube of NGC 5643. The morphologies of the ionized and molecular gas seem to be compatible with the existence of a molecular torus/disc that collimates the active galactic nucleus (AGN) emission. The molecular gas may also flow along this torus/disc, feeding the AGN. This scenario is compatible with the unified model for AGNs.
We describe an algorithm for identifying point-source transients and moving objects on reference-subtracted optical images containing artifacts of processing and instrumentation. The algorithm makes use of the supervised machine learning technique known as Random Forest. We present results from its use in the Dark Energy Survey Supernova program (DES-SN), where it was trained using a sample of 898,963 signal and background events generated by the transient detection pipeline. After reprocessing the data collected during the first DES-SN observing season (Sep. 2013 through Feb. 2014) using the algorithm, the number of transient candidates eligible for human scanning decreased by a factor of 13.4, while only 1 percent of the artificial Type Ia supernovae (SNe) injected into search images to monitor survey efficiency were lost, most of which were very faint events. Here we characterize the algorithm's performance in detail, and we discuss how it can inform pipeline design decisions for future time-domain imaging surveys, such as the Large Synoptic Survey Telescope and the Zwicky Transient Facility.
Accurate line lists are calculated for aluminium monoxide covering the pure rotation, rotation-vibration and electronic (B -- X blue-green and A -- X infrared bands) spectrum. Line lists are presented for the main isotopologue, $^{27}$Al$^{16}$O, as well as for $^{27}$Al$^{17}$O, $^{27}$Al$^{18}$O and $^{26}$Al$^{16}$O. These line lists are suitable for high temperatures (up to 8000 K) including those relevant to exoplanetary atmospheres and cool stars. A combination of empirical and \textit{ab initio} methods is used: the potential energy curves were previously determined to high accuracy by fitting to extensive data from analysis of laboratory spectra; a high quality {\it ab initio} dipole moment curve is calculated using quadruple zeta basis set and the multi-reference configuration interaction (MRCI) method. Partition functions plus full line lists of transitions are made available in an electronic form as supplementary data to this article and at \url{www.exomol.com}.
We present results of numerical simulations of flux and linear polarization variations in transiting exoplanetary systems, caused by the host star disk symmetry breaking. We consider different configurations of planetary transits depending on orbital parameters. Starspot contribution to the polarized signal is also estimated. Applying the method to known systems and simulating observational conditions, a number of targets is selected where transit polarization effects could be detected. We investigate several principal benefits of the transit polarimetry, particularly, for determining orbital spatial orientation and distinguishing between grazing and near-grazing planets. Simulations show that polarization parameters are also sensitive to starspots, and they can be used to determine spot positions and sizes.
With the advent of the imaging atmospheric Cherenkov technique in late 1980's, ground-based observations of Very High-Energy gamma rays came into reality. Since the first source detected at TeV energies in 1989 by Whipple, the number of high energy gamma-ray sources has rapidly grown up to more than 150 thanks to the second generation experiments like MAGIC, H.E.S.S. and VERITAS. The Cherenkov Telescope Array observatory is the next generation of Imaging Atmospheric Cherenkov Telescopes, with at least 10 times higher sensitivity than current instruments. Cherenkov Telescopes have to be equipped with optical dishes of large diameter -- in general based on segmented mirrors -- with typical angular resolution of a few arc-minutes. To evaluate the mirror's quality specific metrological systems are required that possibly take into account the environmental conditions in which typically Cherenkov telescopes operate (in open air without dome protection). For this purpose a new facility for the characterization of mirrors has been developed at the labs of the Osservatorio Astronomico di Brera of the Italian National Institute of Astrophysics. The facility allows the precise measurement of the radius of curvature and the distribution of the concentred light in terms of focused and scattered components and it works in open air. In this paper we describe the facility and report some examples of its measuring capabilities from recent tests campaigns carried out performed on mirrors devoted to Cherenkov telescopes.
The accretion histories of embedded protostars are an integral part of descriptions of their physical and chemical evolution. In particular, are the accretion rates smoothly declining from the earlier toward later stages or in fact characterized by variations such as intermittent bursts? We aim to characterize the impact of possible accretion variations in a sample of embedded protostars by measuring the size of the inner regions of their envelopes where CO is sublimated and relate those to their temperature profiles dictated by their current luminosities. Using observations from the Submillimeter Array we measure the extents of the emission from the C18O isotopologue toward 16 deeply embedded protostars. We compare these measurements to the predicted extent of the emission given the current luminosities of the sources through dust and line radiative transfer calculations. Eight out of sixteen sources show more extended C18O emission than predicted by the models. The modeling shows that the likely culprit for these signatures is sublimation due to increases in luminosities of the sources by about a factor five or more during the recent 10,000 years - the time it takes for CO to freeze-out again on dust grains. For four of those sources the increase would have had to have been a factor 10 or more. The compact emission seen toward the other half of the sample suggests that C18O only sublimates when the temperature exceeds 30 K - as one would expect if CO is mixed with H2O in the grain ice-mantles. The small-number statistics from this survey suggest that protostars undergo significant bursts about once every 20,000 years. This also illustrates the importance of taking the physical evolutionary histories into account for descriptions of the chemical structures of embedded protostars.
Using the science verification data of the Dark Energy Survey (DES) for a new sample of 106 X-Ray selected clusters and groups, we study the stellar mass growth of Bright Central Galaxies (BCGs) since redshift 1.2. Compared with the expectation in a semi-analytical model applied to the Millennium Simulation, the observed BCGs become under-massive/under-luminous with decreasing redshift. We incorporate the uncertainties associated with cluster mass, redshift, and BCG stellar mass measurements into analysis of a redshift-dependent BCG-cluster mass relation, $m_{*}\propto(\frac{M_{200}}{1.5\times 10^{14}M_{\odot}})^{0.24\pm 0.08}(1+z)^{-0.19\pm0.34}$, and compare the observed relation to the simulation prediction. We estimate the average growth rate since z = 1.0 for BCGs hosted by clusters of $M_{200, z}=10^{13.8}M_{\odot}$, at $z=1.0$: $m_{*, BCG}$ appears to have grown by $0.13\pm0.11$ dex, in tension at $\sim 2.5 \sigma$ significance level with the 0.4 dex growth rate expected in the simulation. We show that the buildup of extended intra-cluster light after $z=1.0$ may alleviate this tension in BCG growth rates.
The Dark Energy Survey (DES) is currently undertaking an observational program imaging $1/4$ of the southern hemisphere sky with unprecedented photometric accuracy. In the process of observing millions of faint stars and galaxies to constrain the parameters of the dark energy equation of state, the DES will obtain pre-discovery images of the regions surrounding an estimated 100 gamma-ray bursts (GRBs) over five years. Once GRBs are detected by, e.g., the Swift satellite, the DES data will be extremely useful for follow-up observations by the transient astronomy community. We describe a recently-commissioned suite of software that listens continuously for automated notices of GRB activity, collates useful information from archival DES data, and promulgates relevant data products back to the community in near-real-time. Of particular importance are the opportunities that DES data provide for relative photometry of GRBs or their afterglows, as well as for identifying key characteristics (e.g., photometric redshifts) of potential GRB host galaxies. We provide the functional details of the DESAlert software as it presently operates, as well as the data products that it produces, and we show sample results from the application of DESAlert to several previously-detected GRBs.
We advance an exact, explicit form for the solutions to the fractional diffusion-advection equation. Numerical analysis of this equation shows that its solutions resemble power-laws.
Weak gravitational lensing allows one to reconstruct the spatial distribution of the projected mass density across the sky. These "mass maps" provide a powerful tool for studying cosmology as they probe both luminous and dark matter. In this paper, we present a weak lensing mass map reconstructed from shear measurements in a 139 deg^2 area from the Dark Energy Survey (DES) Science Verification (SV) data overlapping with the South Pole Telescope survey. We compare the distribution of mass with that of the foreground distribution of galaxies and clusters. The overdensities in the reconstructed map correlate well with the distribution of optically detected clusters. Cross-correlating the mass map with the foreground galaxies from the same DES SV data gives results consistent with mock catalogs that include the primary sources of statistical uncertainties in the galaxy, lensing, and photo-z catalogs. The statistical significance of the cross-correlation is at the 6.8 sigma level with 20 arcminute smoothing. A major goal of this study is to investigate systematic effects arising from a variety of sources, including PSF and photo-z uncertainties. We make maps derived from twenty variables that may characterize systematics and find the principal components. We find that the contribution of systematics to the lensing mass maps is generally within measurement uncertainties. We test and validate our results with mock catalogs from N-body simulations. In this work, we analyze less than 3% of the final area that will be mapped by the DES; the tools and analysis techniques developed in this paper can be applied to forthcoming larger datasets from the survey.
We have obtained deep NIR narrow and broad (J and Y) band imaging data of the GOODS-South field. The narrow band filter is centered at 1060 nm corresponding to redshifts $z = 0.62, 1.15, 1.85$ for the strong emission lines H$\alpha$, $[$OIII$]$/H$\beta$ and $[$OII$]$, respectively. From those data we extract a well defined sample ($M(AB)=24.8$ in the narrow band) of objects with large emission line equivalent widths in the narrow band. Via SED fits to published broad band data we identify which of the three lines we have detected and assign redshifts accordingly. This results in a well defined, strong emission line selected sample of galaxies down to lower masses than can easily be obtained with only continuum flux limited selection techniques. We compare the (SED fitting-derived) main sequence of star-formation (MS) of our sample to previous works and find that it has a steeper slope than that of samples of more massive galaxies. We conclude that the MS steepens at lower (below $M_{\star} = 10^{9.4} M_{\odot}$) galaxy masses. We also show that the SFR at any redshift is higher in our sample. We attribute this to the targeted selection of galaxies with large emission line equivalent widths, and conclude that our sample presumably forms the upper boundary of the MS. We briefly investigate and outline how samples with accurate redshifts down to those low stellar masses open a new window to study the formation of large scale structure in the early universe. In particular we report on the detection of a young galaxy cluster at $z=1.85$ which features a central massive galaxy which is the candidate of an early stage cD galaxy, and we identify a likely filament mapped out by $[$OIII$]$ and $H\beta$ emitting galaxies at $z=1.15$.
Non-rotating neutron stars are generally treated in theoretical studies as perfect spheres. Such a treatment, however, may not be correct if strong magnetic fields are present (such as for magnetars) and/or the pressure of the matter in the cores of neutron stars is non-isotropic (e.g., color superconducting). In this paper, we investigate the structure of non-spherical neutron stars in the framework of general relativity. Using a parameterized metric to model non-spherical mass distributions, we first derive a stellar structure equation for deformed neutron stars. Numerical investigations of this model equation show that the gravitational masses of deformed neutron stars depend rather strongly on the degree and type (oblate or prolate) of stellar deformation. In particular, we find that the mass of a neutron star increases with increasing oblateness but decreases with increasing prolateness. If this feature carries over to a full two-dimensional treatment of deformed neutron stars, this opens up the possibility that, depending on the type of stellar deformation, there may exist multiple maximum-mass neutron stars for one and for the same model for the nuclear equation of state.
Blazars are astrophysical sources whose emission is dominated by non-thermal processes, typically interpreted as synchrotron and inverse Compton emission. Although the general picture is rather robust and consistent with observations, many aspects are still unexplored. Polarimetric monitoring can offer a wealth of information about the physical processes in blazars. Models with largely different physical ingredients can often provide almost indistinguishable predictions for the total flux, but usually are characterized by markedly different polarization properties. We explore, with a pilot study, the possibility to derive structural information about the emitting regions of blazars by means of a joint analysis of rapid variability of the total and polarized flux at optical wavelengths. Short timescale (from tens of seconds to a couple of minutes) optical linear polarimetry and photometry for two blazars, BL Lacertae and PKS 1424+240, was carried out with the PAOLO polarimeter at the 3.6m Telescopio Nazionale Galileo. Several hours of almost continuous observations were obtained for both sources. Our intense monitoring allowed us to draw strongly different scenarios for BL Lacertae and PKS 1424+240, with the former characterized by intense variability on time-scales from hours to a few minutes and the latter practically constant in total flux. Essentially the same behavior is observed for the polarized flux and the position angle. The variability time-scales turned out to be as short as a few minutes, although involving only a few percent variation of the flux. The polarization variability time-scale is generally consistent with the total flux variability. Total and polarized flux appear to be essentially uncorrelated. However, even during our relatively short monitoring, different regimes can be singled out. (abridged)
Reconstruction of the local velocity field from the overdensity field and a gravitational acceleration that falls off from a point mass as r^-2 yields velocities in broad agreement with peculiar velocities measured with galaxy distance indicators. MONDian gravity does not. To quantify this, we introduce the velocity angular correlation function as a diagnostic of peculiar velocity field alignment and coherence as a function of scale. It is independent of the bias parameter of structure formation in the standard model of cosmology and the acceleration parameter of MOND. A modified gravity acceleration consistent with observed large scale structure would need to asymptote to zero at large distances more like r^-2, than r^-1.
As part of the OzDES spectroscopic survey we are carrying out a large scale reverberation mapping study of $\sim$500 quasars over five years in the 30 deg$^2$ area of the Dark Energy Survey (DES) supernova fields. These quasars have redshifts ranging up to 4 and have apparent AB magnitudes between $16.8<r<22.5$ mag. The aim of the survey is to measure time lags between fluctuations in the quasar continuum and broad emission line fluxes of individual objects in order to measure black hole masses for a broad range of AGN and constrain the radius-luminosity ($R-L$) relationship. Here we investigate the expected efficiency of the OzDES reverberation mapping campaign and its possible extensions. We expect to recover lags for $\sim$35-45\% of the quasars. AGN with shorter lags and greater variability are more likely to yield a lag, and objects with lags $\lesssim$6 month or $\sim$1 year are expected be recovered the most accurately. The baseline OzDES reverberation mapping campaign is predicted to produce an unbiased measurement of the $R-L$ relationship parameters for H$\beta$, MgII $\lambda$2798, and CIV $\lambda$1549. However, extending the baseline survey by either increasing the spectroscopic cadence, extending the survey season, or improving the emission line flux measurement accuracy will significantly improve the $R-L$ parameter constraints for all broad emission lines.
We carried out 2D-axisymmetric MHD simulations of core-collapse supernovae for rapidly-rotating magnetized progenitors. By changing both the strength of the magnetic field and the spatial resolution, the evolution of the magnetorotational instability (MRI) and its impacts upon the dynamics are investigated. We found that the MRI greatly amplifies the seed magnetic fields in the regime where not the Alfv\'en mode but the buoyant mode plays a primary role in the exponential growth phase. The MRI indeed has a powerful impact on the supernova dynamics. It makes the shock expansion faster and the explosion more energetic, with some models being accompanied by the collimated-jet formations. These effects, however, are not made by the magnetic pressure except for the collimated-jet formations. The angular momentum transfer induced by the MRI causes the expansion of the heating region, by which the accreting matter gain an additional time to be heated by neutrinos. The MRI also drifts low-$Y_p$ matter from the deep inside of the core to the heating region, which makes the net neutrino heating rate larger by the reduction of the cooling due to the electron capture. These two effects enhance the efficiency of the neutrino heating, which is found to be the key to boost the explosion. Indeed we found that our models explode far more weakly when the net neutrino heating is switched off. The contribution of the neutrino heating to the explosion energy is estimated to be more than 90\%.
OzDES is a five-year, 100-night, spectroscopic survey on the Anglo-Australian
Telescope, whose primary aim is to measure redshifts of approximately 2,500
Type Ia supernovae host galaxies over the redshift range 0.1 < z < 1.2, and
derive reverberation-mapped black hole masses for approximately 500 active
galactic nuclei and quasars over 0.3 < z < 4.5. This treasure trove of data
forms a major part of the spectroscopic follow-up for the Dark Energy Survey
for which we are also targeting cluster galaxies, radio galaxies, strong
lenses, and unidentified transients, as well as measuring luminous red galaxies
and emission line galaxies to help calibrate photometric redshifts.
Here we present an overview of the OzDES program and our first-year results.
Between Dec 2012 and Dec 2013, we observed over 10,000 objects and measured
more than 6,000 redshifts. Our strategy of retargeting faint objects across
many observing runs has allowed us to measure redshifts for galaxies as faint
as m_r=25 mag. We outline our target selection and observing strategy, quantify
the redshift success rate for different types of targets, and discuss the
implications for our main science goals. Finally, we highlight a few
interesting objects as examples of the fortuitous yet not totally unexpected
discoveries that can come from such a large spectroscopic survey.
We observed three massive subhalos in the Coma cluster with {\it Suzaku}. These subhalos, labeled "ID 1", "ID 2", and "ID 32", were detected with a weak-lensing survey using the Subaru/Suprime-Cam (Okabe et al. 2014a), and are located at the projected distances of 1.4 $r_{500}$, 1.2 $r_{500}$, and 1.6 $r_{500}$ from the center of the Coma cluster, respectively. The subhalo "ID 1" has a compact X-ray excess emission close to the center of the weak-lensing mass contour, and the gas mass to weak-lensing mass ratio is about 0.001. The temperature of the emission is about 3 keV, which is slightly lower than that of the surrounding intracluster medium (ICM) and that expected for the temperature vs. mass relation of clusters of galaxies. The subhalo "ID 32" shows an excess emission whose peak is shifted toward the opposite direction from the center of the Coma cluster. The gas mass to weak-lensing mass ratio is also about 0.001, which is significantly smaller than regular galaxy groups. The temperature of the excess is about 0.5 keV and significantly lower than that of the surrounding ICM and far from the temperature vs. mass relation of clusters. However, there is no significant excess X-ray emission in the "ID 2" subhalo. Assuming an infall velocity of about 2000 $\rm km~s^{-1}$, at the border of the excess X-ray emission, the ram pressures for "ID 1" and "ID 32" are comparable to the gravitational restoring force per area. We also studied the effect of the Kelvin-Helmholtz instability to strip the gas. Although we found X-ray clumps associated with the weak-lensing subhalos, their X-ray luminosities are much lower than the total ICM luminosity in the cluster outskirts.
Cosmological applications of the "redshift - angular size" test require knowledge of the linear size of the "standard rod" used. In this paper, we study the properties of a large sample of 140 milliarcsecond compact radio sources with flux densities measured at 6 cm and 20 cm, compiled by Gurvits et al.(1999). Using the best-fitted cosmological parameters given by Planck/WMAP9 observations, we investigate the characteristic length $l_m$ as well as its dependence on the source luminosity $L$ and redshift $l_m=l L^\beta (1+z)^n$. For the full sample, measurements of the angular size $\theta$ provide a tight constraint on the linear size parameters. We find that cosmological evolution of the linear size is small ($|n|\simeq 10^{-2}$) and consistent with previous analysis. However, a substantial evolution of linear sizes with luminosity is still required ($\beta\simeq 0.17$). Furthermore, similar analysis done on sub-samples defined by different source optical counterparts and different redshift ranges, seems to support the scheme of treating radio galaxies and quasars with distinct strategies. Finally, a cosmological-model-independent method is discussed to probe the properties of angular size of milliarcsecond radio quasars. Using the corrected redshift - angular size relation for quasar sample, we obtained a value of the matter density parameter, $\Omega_m=0.292^{+0.065}_{-0.090}$, in the spatially flat $\Lambda$CDM cosmology.
We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of 182 stellar targets along the line of sight to the newly-discovered `ultrafaint' object Reticulum 2 (Ret 2). For 38 of these targets, the spectra are sufficient to provide simultaneous estimates of line-of-sight velocity ($v_{\rm los}$, median random error $\delta_{v_{\rm los}}=1.3$ km s$^{-1}$), effective temperature ($T_{\rm eff}$, $\delta_{T_{\rm eff}}=464$ K), surface gravity ($\log g$, $\delta_{\rm logg}=0.54$ dex) and iron abundance ([Fe/H], $\delta_{\mathrm{[Fe/H]}}=0.45$ dex). We use these results to confirm 18 stars as members of Ret 2. From the member sample we estimate a velocity dispersion of $\sigma_{v_{\rm los}}=3.6_{-0.6}^{+0.9}$ km s$^{-1}$ about a mean of $\langle v_{\rm los}\rangle =64.8_{-1.0}^{+1.1}$ km s$^{-1}$ in the solar rest frame ($\sim -90.9$ km s$^{-1}$ in the Galactic rest frame), and a metallicity dispersion of $\sigma_{\rm [Fe/H]}=0.50_{-0.13}^{+0.17}$ dex about a mean of $\langle \mathrm{[Fe/H]} \rangle =-2.67_{-0.34}^{+0.34}$. These estimates marginalize over possible velocity and metallicity gradients, which are consistent with zero. Our results place Ret 2 on chemodynamical scaling relations followed by the Milky Way's dwarf-galactic satellites. Under assumptions of dynamic equilibrium and negligible contamination from binary stars---both of which must be checked with deeper imaging and repeat spectroscopic observations---the estimated velocity dispersion suggests a dynamical mass of $M(R_{\rm h})\approx 5R_{\rm h}\sigma_{v_{\rm los}}^2/(2G)=2.4_{-0.8}^{+1.3}\times 10^5$ $M_{\odot}$ enclosed within projected halflight radius $R_{\rm h}\sim 32$ pc, with mass-to-light ratio $\approx 2M(R_{\rm h})/L_{\rm V}=462_{-157}^{+264}$ in solar units.
We report on monitoring observations of the TeV gamma-ray binary HESS J0632+057, which were carried out to constrain the interaction between the Be circumstellar disk and the compact object of unknown nature, and provide for the first time high-dispersion (R > 50000) optical spectra in the second half of the orbital cycle, from apastron through periastron. The Halpha, Hbeta, and Hgamma line profiles are found to exhibit remarkable short-term variability for ~1 month after the apastron (phase 0.6--0.7), whereas they show little variation near the periastron. These emission lines show "S-shaped" variations with timescale of ~150 days, which is about twice that reported previously. In contrast to the Balmer lines, no profile variability is seen in any FeII emission line. We estimate the radii of emitting regions of the Halpha, Hbeta, Hgamma, and FeII emission lines to be ~30, 11, 7, and 2 stellar radii (R_*), respectively. The amplitudes of the line profile variations in different lines indicate that the interaction with the compact object affects the Be disk down to, at least, the radius of 7 R_* after the apastron. This fact, together with little profile variability near the periastron, rules out the tidal force as the major cause of disk variability. Although this leaves the pulsar wind as the most likely candidate mechanism for disk variations, understanding the details of the interaction, particularly the mechanism for causing a large-scale disk disturbance after the apastron, remains an open question.
Hard xray data from the RXTE observatory (HEXTE energy range 15 to 240 keV) have been analyzed to obtain a phase coherent timing solution for the Crab pulsar glitch of 15 July 2000. The results are: (1) step change in the rotation frequency $\nu_0$ of the Crab pulsar at the epoch of the glitch is $\Delta \nu_0 = (30 \pm 3) \times 10^{-9} \times \nu_0$, (2) step change in its time derivative is $\Delta \dot \nu_0 = (4.8 \pm 0.6) \times 10^{-3} \times \dot \nu_0$, and (3) the time scale of decay of the the step change is $\tau_d = 4.7 \pm 0.5$ days. The first two results are consistent with those obtained at radio frequencies by the Jodrell Bank observatory. The last result has not been quoted in the literature, but could be an underestimate due to lack of observations very close to the glitch epoch. By comparing with the monthly timing ephemeris published by the Jodrell group for the Crab pulsar, the time delay between the main peaks of the hard xray and radio pulse profiles is estimated to be $+411 \pm 167$ $\mu$sec. Although this number is not very significant, it is consistent with the number derived for the 2 to 16 keV energy range, using the PCA instrument of RXTE. The separation between the two peaks of the integrated pulse profile of the Crab pulsar, and the ratio of their intensities, both are statistically similar before and after the glitch. The dead time corrected integrated photon flux within the integrated pulse profile appears to decrease after the glitch, although this is not a statistically strong result. This work achieves what can be considered to be almost absolute timing analysis of the Crab pulsar hard xray data.
The eruption of the filament with the kink fashion is often regarded as a signature of the kink instability. However, the kink instability threshold for the filament magnetic structure has been not widely understood. Using the H-alpha observation from the New Vacuum Solar Telescope (NVST), we present a partial eruptive filament. In the eruption, a filament thread appeared to split from the middle portion of the filament and to break out in a kinklike fashion. During this period, the left filament material remained below, which erupted without the kinking motion later on. The coronal magnetic field lines associated with the filament are obtained from the nonlinear force-free field (NLFFF) extrapolations using the 12 minutes cadence vector magnetograms of the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamic Observatory (SDO). We studied the extrapolated field lines passing through the magnetic dips that are in good agreement with the observed filament. The field lines are non-uniformly twisted and appear to be made up by two twisted flux ropes winding about each other. One of them has higher twist than the other, and the highly twisted one has its dips aligned with the kinking eruptive thread at the beginning of its eruption. Before the eruption, moreover, the highly twisted flux rope was found to expand with the approximately constant field twist. In addition, the helicity flux maps deduced from the HMI magnetograms show that some helicity is injected into the overlying magnetic arcade, but no significant helicity is injected into the flux ropes. Accordingly, we suggest that the highly twisted flux rope became kink unstable when the instability threshold declined with the expansion of the flux rope.
Shock waves in clusters of galaxies have been observed as radio relics and also as discontinuities in X-ray. While it has been suggested that radio relics could a consequence of mergers, merger shocks are expected to be weak ($M_s \la 3$), so may be difficult to account for several relics with flat radio spectra. Moreover, for some radio relics, there are inconsistencies in properties such as shock Mach number and location derived from radio and X-ray observations. In this paper, using the data of cosmological hydrodynamic simulations for the large-scale structure formation of the Universe, we construct mock radio and X-ray maps of simulated clusters, projected onto the two-dimensional (2D) plane, and study various properties of shocks and synthetic radio relics in the maps that could be inferred from observations. We find that a substantial fraction of radio and X-ray shocks in 2D maps involve multiple shock surfaces along line of sights (LoSs), and the morphology of shock distributions in 2D maps depends on the projection direction. A fraction of synthetic radio relics include significant contributions from infall shocks with $M_s \ga 3$. Among the multiple shocks in a LoS, radio observations tend to pick up stronger shocks with flatter radio spectra, while X-ray observations preferentially select weaker shocks with larger kinetic energy flux. As a result, the shock Mach numbers and locations derived from radio and X-ray observations could be different from each other.
Aims. Having established the binary status of nineteen O-type stars located in four Cygnus OB associations, we now determine their fundamental parameters to constrain their properties and their evolutionary status. We also investigate their surface nitrogen abundances, which we compare with other results from the literature obtained for galactic O-type stars. Methods. Using optical spectra collected for each object in our sample and some UV data from the archives, we apply the CMFGEN atmosphere code to determine their main properties. For the binary systems, we have disentangled the components to obtain their individual spectra and investigate them as if they were single stars. Results. We find that the distances of several presumably single O-type stars seem poorly constrained because their luminosities are not in agreement with the "standard" luminosities of stars with similar spectral types. The ages of these O-type stars are all less than 7 Myrs. Therefore, the ages of these stars agree with those, quoted in the literature, of the four associations, except for CygOB8 for which the stars seem older than the association itself. However, we point out that the distance of certain stars is debatable relative to values found in the literature. The N content of these stars put in perspective with N contents of several other galactic O-type stars seems to draw the same five groups as found in the "Hunter" diagram for the O and B-type stars in the LMC even though their locations are obviously different. We determine mass-loss rates for several objects from the Halpha line and UV spectra. Finally, we confirm the "mass discrepancy" especially for O stars with masses smaller than 30 Msun. .
Relativistic jets are one of the most powerful manifestations of the release of energy related to the supermassive black holes at the centre of active galactic nuclei (AGN). Their emission is observed across the entire electromagnetic spectrum, from the radio band to gamma rays. Despite decades of efforts, many aspects of the physics of relativistic jets remain elusive. In particular, the location and the mechanisms responsible for the high-energy emission and the connection of the variability at different wavelengths are among the greatest challenges in the study of AGN. Recent high resolution radio observations of flaring objects locate the high-energy emitting region downstream the jet at parsec scale distance from the central engine, posing questions on the nature of the seed photons upscattered to gamma-rays. Furthermore, monitoring campaigns of the most active blazars indicate that not all the high energy flares have the same characteristics in the various energy bands, even from the same source, making the interpretation of the mechanism responsible for the high-energy emission not trivial. Although the variability of the most luminous blazars is well explained by the "shock-in-jet" scenario, the sub-class of TeV emitting objects suggests a more complex emission model with velocity gradients in a structured jet. This contribution presents results obtained by recent multiwavelength campaigns of blazars aimed at studying the radio and gamma-ray connection and the physical mechanisms at the basis of the emission in these low and high energy bands.
High energy resolution spectroscopy of the 1.8 MeV radioactive decay line of 26Al with the SPI instrument on board the INTEGRAL satellite has recently revealed that diffuse 26Al has large velocities in comparison to other components of the interstellar medium in the Milky Way. 26Al shows Galactic rotation in the same sense as the stars and other gas tracers, but reaches excess velocities up to 300 km/s. We investigate if this result can be understood in the context of superbubbles, taking into account the statistics of young star clusters and H I supershells, as well as the association of young star clusters with spiral arms. We derive energy output and 26Al mass of star clusters as a function of the cluster mass via population synthesis from stellar evolutionary tracks of massive stars. [...] We link this to the size distribution of HI supershells and assess the properties of likely 26Al-carrying superbubbles. 26Al is produced by star clusters of all masses above about 200 solar masses, roughly equally contributed over a logarithmic star cluster mass scale, and strongly linked to the injection of feedback energy. The observed superbubble size distribution cannot be related to the star cluster mass function in a straight forward manner. In order to avoid that the added volume of all superbubbles exceeds the volume of the Milky Way, individual superbubbles have to merge frequently. If any two superbubbles merge, or if 26Al is injected off-centre in a bigger HI supershell we expect the hot 26Al-carrying gas to obtain velocities of the order of the typical sound speed in superbubbles, about 300 km/s before decay. [...]
PAMELA and, more recently, AMS-02, are ushering us into a new era of greatly reduced statistical uncertainties in experimental measurements of cosmic-ray fluxes. In particular, new determinations of traditional diagnostic tools such as the boron-to-carbon ratio (B/C) are expected to significantly reduce errors on cosmic-ray diffusion parameters, with important implications for astroparticle physics, ranging from inferring primary source spectra to indirect dark matter searches. It is timely to stress, however, that the conclusions obtained crucially depend on the framework in which the data are interpreted as well as from some nuclear input parameters. We aim at assessing the theoretical uncertainties affecting the outcome, with models as simple as possible while still retaining the key dependencies. We compare different semi-analytical, two-zone model descriptions of cosmic-ray transport in the Galaxy. We test for the effect of a primary source contamination in the boron flux by parametrically altering its flux, as well as for nuclear cross section uncertainties. Our study on preliminary results from AMS-02 suggests that, differently for instance from the leptonic case, realistic modelling of the geometry of the Galaxy and of the source distribution are of minor importance to correctly reproduce B/C data at high energies and thus, to a large extent, for the extraction of diffusion parameters. The Ansatz on the lack of primary injection of boron represents the most serious bias, and requires multi-messenger studies to be addressed. If this uncertainty could be lifted, nuclear uncertainties would still represent a serious concern, which degrade the systematic error on the inferred parameters to the 20% level, or three times the estimated experimental sensitivity. In order to reduce this, a new nuclear cross section measurement campaign is probably required.
The Global mm-VLBI Array (GMVA) is a network of 14 3\,mm and 7\,mm capable telescopes spanning Europe and the United States, with planned extensions to Asia. The array is capable of sensitive maps with angular resolution often exceeding 50\,$\mu$as. Using the GMVA, a large sample of prominent $\gamma$-ray blazars have been observed approximately 6 monthly from later 2008 until now. Combining 3\,mm maps from the GMVA with near-in-time 7\,mm maps from the VLBA-BU-BLAZAR program and 2\,cm maps from the MOJAVE program, we determine the sub-pc morphology and high frequency spectral structure of $\gamma$-ray blazars. The magnetic field strength can be estimated at different locations along the jet under the assumption of equipartition between magnetic field and relativistic particle energies. Making assumptions on the jet magnetic field configuration (e.g. poloidal or toroidal), we can estimate the separation of the mm-wave "core" and the jet base, and estimate the strength of the magnetic field there. The results of this analysis show that on average, the magnetic field strength decreases with a power-law $B \propto r^{-n}$, $n=0.3 \pm 0.2$. This suggests that on average, the mm-wave "core" is $\sim 1-3$\,pc downstream of the de-projected jet apex and that the magnetic field strength is of the order $B_{\rm{apex}} \sim 5-20$\,kG, broadly consistent with the predictions of magnetic jet launching (e.g. via magnetically arrested disks (MAD)).
We find that the formation of MWC 656 (the first Be binary containing a black
hole) involves a common envelope phase and a supernova explosion. This result
supports the idea that a rapidly rotating Be star can emerge out of a common
envelope phase, which is very intriguing because this evolutionary stage is
thought to be too fast to lead to significant accretion and spin up of the B
star.
We predict $\sim 10-100$ of Be BH binaries to currently reside in the
Galactic disk, but there is only a small chance to observe a system with
parameters resembling MWC 656. If MWC 656 is representative of intrinsic
Galactic Be BH binary population, it may indicate that standard evolutionary
theory needs to be revised. This would pose another evolutionary problem in
understanding BH binaries, with BH X-ray Novae formation issue being the prime
example.
The future evolution of MWC 656 with a $\sim 5$ M$_{\odot}$ black hole and
with a $\sim 13$ M$_{\odot}$ main sequence companion on a $\sim 60$ day orbit
may lead to the formation of a coalescing BH-NS system. The estimated Advanced
LIGO/Virgo detection rate of such systems is up to $\sim 0.2$ yr$^{-1}$. This
empirical estimate is a lower limit as it is obtained with only one particular
evolutionary scenario, the MWC 656 binary. This is only a third such estimate
available (after Cyg X-1 and Cyg X-3), and it lends additional support to the
existence of so far undetected BH--NS binaries.
T and Y-dwarfs are among the coolest and least luminous objects detected, and they can help to understand the properties of giant planets. Their multiplicity properties can shed light on the formation process. We observed a sample six T dwarfs and one L9 dwarf with the Laser Guide Star (LGS) and NAOS-CONICA, the adaptive optics (AO) facility, and the near infrared camera at the ESO Very Large Telescope. From the seven observed objects, we have identified a subarcsecond binary system, WISE~J0612-3036, composed of two similar components with spectral types of T6. We measure a separation of $\rho$ = 350$\pm$5 mas and a position angle of $PA = 235\pm1^{\circ}$. Using the mean absolute magnitudes of T6 dwarfs in the 2MASS $JHK_s$ bands, we estimate a distance of $d$=31$\pm$6 pc and derive a projected separation of $\rho$ = 11$\pm$2 au. Another target, WISE J2255-3118, shows a very faint object at 1.3 arcsec in the $K_{\rm s}$ image. The object is marginally detected in $H$, and we derive a near infrared color of $H-K_{\rm s}$$>$ 0.1\,mag. $HST/WFC3$ public archival data reveals that the companion candidate is an extended source.Together with the derived color, this suggests that the source is most probably a background galaxy. The five other sources are apparently single, with 3-$\sigma$ sensitivity limits between $H$=19-21 for companions at separations $\geq$ 0.5 arcsec. WISE 0612-3036 is probably a new T-dwarf binary composed of two T6 dwarfs. As in the case of other late T-dwarf binaries, it shows a mass ratio close to 1, although its projected separation, $\sim$11 au, is larger than the average ($\sim$ 5 au). Additional observations are needed to confirm that the system is bound.
In the past years, several observations of AGN and X-ray binaries have suggested the existence of a warm T around 0.5-1 keV and optically thick, \tau ~ 10-20, corona covering the inner parts of the accretion disk. These properties are directly derived from spectral fitting in UV to soft-X-rays using Comptonization models. However, whether such a medium can be both in radiative and hydrostatic equilibrium with an accretion disk is still uncertain. We investigate the properties of such warm, optically thick coronae and put constraints on their existence. We solve the radiative transfer equation for grey atmosphere analytically in a pure scattering medium, including local dissipation as an additional heating term in the warm corona. The temperature profile of the warm corona is calculated assuming it is cooled by Compton scattering, with the underlying dissipative disk providing photons to the corona. Our analytic calculations show that a dissipative thick, (\tau_{cor} ~ 10-12) corona on the top of a standard accretion disk can reach temperatures of the order of 0.5-1 keV in its upper layers provided that the disk is passive. But, in absence of strong magnetic fields, the requirement of a Compton cooled corona in hydrostatic equilibrium in the vertical direction sets an upper limit on the Thomson optical depth \tau_{cor} < 5 . We show this value cannot be exceeded independently of the accretion disk parameters. However, magnetic pressure can extend this result to larger optical depths. Namely, a dissipative corona might have an optical depth up to ~ 20 when the magnetic pressure is 100 times higher that the gas pressure. The observation of warm coronae with Thomson depth larger than ~ 5 puts tights constraints on the physics of the accretion disk/corona systems and requires either strong magnetic fields or vertical outflows to stabilize the system.
Svestka (Solar Phys. 1989, 121, 399) on the basis of the Solar Maximum Mission observations introduced a new class of flares, the so-called flare hybrids. When they start, they look as typical compact flares (phase 1), but later on they look like flares with arcades of magnetic loops (phase 2). We summarize the features of flare hybrids in soft and hard X-rays as well as in extreme-ultraviolet; these allow us to distinguish them from other flares. Additional energy release or long plasma cooling timescales have been suggested as possible cause of phase 2. Estimations of frequency of flare hybrids have been given. Magnetic configurations supporting their origin have been presented. In our opinion, flare hybrids are quite frequent and a difference between lengths of two interacting systems of magnetic loops is a crucial parameter for recognizing their features.
The abundance of galaxy clusters is in principle a powerful tool to constrain cosmological parameters, especially $\Omega_\mathrm{m}$ and $\sigma_8$, due to the exponential dependence in the high-mass regime. While the best observables are the X-ray temperature and luminosity, the abundance of galaxy clusters, however, is conventionally predicted as a function of mass. Hence, the intrinsic scatter and the uncertainties in the scaling relations between mass and either temperature or luminosity lower the reliability of galaxy clusters to constrain cosmological parameters. In this article we further refine the X-ray temperature function for galaxy clusters by Angrick et al., which is based on the statistics of perturbations in the cosmic gravitational potential and proposed to replace the classical mass-based temperature function, by including a refined analytic merger model and compare the theoretical prediction to results from a cosmological hydrodynamical simulation. Although we find already a good agreement if we compare with a cluster temperature function based on the mass-weighted temperature, including a redshift-dependent scaling between mass-based and spectroscopic temperature yields even better agreement between theoretical model and numerical results. Incorporating this additional scaling in our model, we constrain the cosmological parameters $\Omega_\mathrm{m}$ and $\sigma_8$ from an X-ray sample of galaxy clusters and find agreement with the recent CMB-based results from the Planck mission at 1$\sigma$-level.
Stellar mixing length theory is modified to include the effects of a nongradient term that originates from the motion of convective elements with entropy perturbations of either sign. It is argued that such a term, first studied by Deardorff in the meteorological context, represents the effects of thin intense downdrafts caused by the rapid cooling in the granulation layer at the top of the convection zone. They transport heat nonlocally, as originally anticipated by Spruit in the 1990s, who describes the convection in the strongly stratified simulations of Stein & Nordlund as entropy rain. Although our model has ill-determined free parameters, it demonstrates that solutions can be found that look similar to the original ones, except that the deeper layers are now Schwarzschild stable, so no giant cells are produced and the typical convective scale is that of granules even at larger depth. Consequences for modeling solar differential, the global dynamo, and sunspots are briefly discussed.
An increasing body of data reveals a one-to-one linear correlation between galaxy halo mass and the total mass in its globular cluster (GC) population, M_{GCS} ~ M_h^{1.03 \pm 0.03}, valid over 5 orders of magnitude. We explore the nature of this correlation for galaxies of different morphological types, and for the subpopulations of metal-poor (blue) and metal-rich (red) GCs. For the subpopulations of different metallicity we find M_{GCS}(blue) ~ M_h^{0.96 \pm 0.03} and M_{GCS}(red) ~ M_h^{1.21 \pm 0.03} with similar scatter. The numerical values of these exponents can be derived from the detailed behavior of the red and blue GC fractions with galaxy mass and provide a self-consistent set of relations. In addition, all morphological types (E, S0, S/Irr) follow the same relation, but with a second-order trend for spiral galaxies to have a slightly higher fraction of metal-rich GCs for a given mass. These results suggest that the amount of gas available for GC formation at high redshift was in nearly direct proportion to the dark-matter halo potential, in strong contrast to the markedly nonlinear behavior of total stellar mass versus halo mass. Of the few available theoretical treatments that directly discuss the formation of GCs in a hierarchical merging framework,the model of Kravtsov & Gnedin (2005) best matches these observations. They find that the blue, metal-poor GCs formed in small halos at $z > 3$ and did so in nearly direct proportion to halo mass. Similar models addressing the formation rate of the red, metal-richer GCs in the same detail and continuing to lower redshift are still needed for a comprehensive picture.
Far-infrared Herschel PACS imaging and spectroscopic observations of the nebula around the luminous blue variable (LBV) star AG Car have been obtained along with optical imaging in the Halpha+[NII] filter. In the infrared light, the nebula appears as a clumpy ring shell that extends up to 1.2 pc with an inner radius of 0.4 pc. It coincides with the Halpha nebula, but extends further out. Dust modeling of the nebula was performed and indicates the presence of large grains. The dust mass is estimated to be ~ 0.2 Msun. The infrared spectrum of the nebula consists of forbidden emission lines over a dust continuum. Apart from ionized gas, these lines also indicate the existence of neutral gas in a photodissociation region that surrounds the ionized region. The abundance ratios point towards enrichment by processed material. The total mass of the nebula ejected from the central star amounts to ~ 15 Msun, assuming a dust-to-gas ratio typical of LBVs. The abundances and the mass-loss rate were used to constrain the evolutionary path of the central star and the epoch at which the nebula was ejected, with the help of available evolutionary models. This suggests an ejection during a cool LBV phase for a star of ~ 55 Msun with little rotation.
We develop a numerical formalism for calculating the distribution with energy
of the (internal) pairs formed in a relativistic source from unscattered
MeV--TeV photons.
For GRB afterglows, this formalism is more suitable if the relativistic
reverse-shock that energizes the ejecta is the source of the GeV photons.
The number of pairs formed is set by the source GeV output (calculated from
the Fermi-LAT fluence), the unknown source Lorentz factor, and the unmeasured
peak energy of the LAT spectral component.
We show synchrotron and inverse-Compton light-curves expected from pairs
formed in the shocked medium and identify some criteria for testing a pair
origin of GRB optical counterparts.
Pairs formed in bright LAT afterglows with a Lorentz factor in the few
hundreds may produce bright optical counterparts (R < 10) lasting for up to one
hundred seconds.
The number of internal pairs formed from unscattered seed photons decreases
very strongly with the source Lorentz factor, thus bright GRB optical
counterparts cannot arise from internal pairs if the afterglow Lorentz factor
is above several hundreds.
We performed time lag analysis on the X-ray light curves of Centaurus A (Cen A) obtained by the Gas Slit Camera (GSC) aboard the Monitor of All-sky X-ray Image (MAXI) in three energy bands (2--4 keV, 4--10 keV, and 10--20 keV). We discovered a soft X-ray lag relative to higher energies (soft lag) on a time scale of days by employing the discrete correlation function (DCF) and the z-transformed discrete correlation function (ZDCF) method in a flare episode. A peak in the DCF and the ZDCF was observed at a soft lag of $\sim 5$ days in 2--4 keV versus 4--10 keV and in 4--10 keV versus 10--20 keV, and $\sim 10$ days in 2--4 keV versus 10--20 keV. We found it difficult to explain the observed X-ray variation with the one-zone synchrotron self-Compton (SSC) model, in which the soft lags reflect the different cooling times of the relativistic electrons in these three energy bands. Alternatively, if the X-ray variation was produced in a corona surrounding or along the inner part of the accretion disk, we can explain this phenomenon due to different Compton cooling times of the electrons with different energies.
Fermi-LAT analyses show that the gamma-ray photon spectral indices Gamma_gamma of a large sample of blazars correlate with the vFv peak synchrotron frequency v_s according to the relation Gamma_gamma = d - k log v_s. The same function, with different constants d and k, also describes the relationship between Gamma_gamma and peak Compton frequency v_C. This behavior is derived analytically using an equipartition blazar model with a log-parabola description of the electron energy distribution (EED). In the Thomson regime, k = k_EC = 3b/4 for external Compton processes and k = k_SSC = 9b/16 for synchrotron self-Compton (SSC) processes, where b is the log-parabola width parameter of the EED. The BL Lac object Mrk 501 is fit with a synchrotron/SSC model given by the log-parabola EED, and is best fit away from equipartition. Departures from equipartition are limited by jet power and make small corrections to the spectral-index diagrams. Analytic expressions are compared with numerical values derived from self-Compton and external Compton scattered gamma-ray spectra calculated for Ly alpha broad-line region and IR target photons. The Gamma_gamma vs. nu_s behavior in the model depends strongly on b, with progressively and predictably weaker dependences on gamma-ray detection range, variability time, and isotropic gamma-ray luminosity. Implications of the results for blazar unification and blazars as ultra-high energy cosmic-ray sources are presented.
The statistical study of the parsec scale properties of radio sources is crucial to get information on the nature of the central engine and to provide the foundations of the current unified theories, suggesting that the appearance of active galactic nuclei depends strongly on orientation. We started a project to observe at sub-arcsec resolution a complete sample of 94 nearby (z<0.1) radio galaxies, the Bologna Complete Sample, which is not affected by any selection effect on the jet velocity and orientation with respect to the line of sight. Up to now, we published our parsec scale analysis of 77/94 sources. Here, we describe the last VLBA observations at 5 GHz and EVN data at 18 cm obtained for the 17 remaining faintest radio core (<5 mJy at 5 GHz in VLA images) BCS sources and we report our preliminary results on the whole complete sample.
We present VLA H I and 6 cm radio continuum observations of the spiral NGC 3145 and H I observations of its two companions, NGC 3143 and PGC 029578. In optical images NGC 3145 has stellar arms that appear to cross, forming "X"-features. Our radio continuum observations rule out shock fronts at 3 of the 4 "X"-features. In the middle-to-outer disk, the H I line-profiles of NGC 3145 are skewed. Relative to the disk, the gas in the skewed wing of the line-profiles has z-motions away from us on the approaching side of the galaxy and z-motions of about the same magnitude (about 40 km/s) towards us on the receding side. These warping motions imply that there has been a perturbation perpendicular to the disk over large spatial scales. Two features in NGC 3145 have velocities indicating that they are out-of-plane tidal arms. One is an apparent branch of a main spiral arm; the velocity of the branch is 150 km/s greater than the spiral arm where they appear to intersect in projection. The other is an arm that forms 3 of the "X"-features. It differs in velocity by 56 km/s from that of the disk at the same projected location. Based on its SFR and H I properties, NGC 3143 is the more likely of the two companions to have interacted with NGC 3145 recently. A simple analytic model demonstrates that an encounter between NGC 3143 and NGC 3145 is a plausible explanation for the observed warping motions in NGC 3145.
A large fraction of stars host one or multiple close-in super-Earth planets. There is an active debate about whether these planets formed in situ or at greater distances from the central star and migrated to their current position. It has been shown that part of their observed properties (e.g., eccentricity distribution) can be reproduced by N-body simulations of in situ formation starting with a population of protoplanets of high masses and neglecting the effects of the disk gas. We plan to reassess the in situ formation of close-in super-Earths through more complete simulations. We performed N-body simulations of a population of small planetary embryos and planetesimals that include the effects of disk-planet interactions (e.g., eccentricity damping, type I migration). In addition, we also consider the accretion of a primitive atmosphere from a protoplanetary disk. We find that planetary embryos grow very quickly well before the gas dispersal, and thus undergo rapid inward migration, which means that one cannot neglect the effects of a gas disk when considering the in-situ formation of close-in super-Earths. Owing to their rapid inward migration, super-Earths reach a compact configuration near the disk's inner edge whose distribution of orbital parameters matches the observed close-in super-Earths population poorly. On the other hand, simulations including eccentricity damping, but no type I migration, reproduce the observed distributions better. Including the accretion of an atmosphere does not help reproduce the bulk architecture of observations. Interestingly, we find that the massive embryos can migrate inside the disk edge while capturing only a moderately massive hydrogen/helium atmosphere. By this process they avoid becoming giant planets. The bulk of close-in super-Earths cannot form in situ, unless type I migration is suppressed in the entire disk inside 1 AU.
We have measured the orbital parameters of seven close binaries, including six new objects, in a radial velocity survey of 38 objects comprising a hot subdwarf star with orbital periods ranging from ~0.17 to 3 d. One new system, GALEX J2205-3141, shows reflection on a M dwarf companion. Three other objects show significant short-period variations, but their orbital parameters could not be constrained. Two systems comprising a hot subdwarf paired with a bright main-sequence/giant companion display short-period photometric variations possibly due to irradiation or stellar activity and are also short-period candidates. All except two candidates were drawn from a selection of subluminous stars in the Galaxy Evolution Explorer ultraviolet sky survey. Our new identifications also include a low-mass subdwarf B star and likely progenitor of a low mass white dwarf (GALEX J0805-1058) paired with an unseen, possibly substellar, companion. The mass functions of the newly identified binaries imply minimum secondary masses ranging from 0.03 to 0.39 M_solar. Photometric time series suggest that, apart from GALEX J0805-1058 and J2205-3141, the companions are most likely white dwarfs. We update the binary population statistics: Close to 40 per cent of hot subdwarfs have a companion. Also, we found that the secondary mass distribution shows a low-mass peak attributed to late-type dwarfs, and a higher-mass peak and tail distribution attributed to white dwarfs and a few spectroscopic composites. Also, we found that the population kinematics imply an old age and include a few likely halo population members.
Recent studies have shown that the extended main-sequence turn off (eMSTO) is
a common feature of intermediate-age star clusters in the Magellanic Clouds
(MCs). The most simple explanation is that these stellar systems harbor
multiple generations of stars with an age difference of a few hundred Myrs.
However, while an eMSTO has been detected in a large number of clusters with
ages between ~1-2 Gyrs, several studies of young clusters in both MCs and in
nearby galaxies do not find any evidence for a prolonged star-formation
history, i.e. for multiple stellar generations. These results have suggested
alternative interpretation of the eMSTOs observed in intermediate-age star
clusters. The eMSTO could be due to stellar rotation mimicking an age spread or
to interacting binaries. In these scenarios, intermediate-age MC clusters would
be simple stellar populations, in close analogy with younger clusters.
Here we provide the first evidence for an eMSTO in a young stellar cluster.
We exploit multi-band Hubble Space Telescope photometry to study the ~300-Myr
old star cluster NGC1856 in the Large Magellanic Cloud and detected a broadened
MSTO that is consistent with a prolonged star-formation which had a duration of
about 150 Myrs. Below the turn-off, the MS of NGC1856 is split into a red and
blue component, hosting 33+/-5% and 67+/-5% of the total number of MS stars,
respectively. We discuss these findings in the context of
multiple-stellar-generation, stellar-rotation, and interacting-binary
hypotheses.
We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include (1) idealized gaussian DEM distributions, (2) 3D models of NOAA Active Region 11158 comprising quasi-steady loop atmospheres in a non-linear force-free field, and (3) thermodynamic models from a fully-compressible, 3D MHD simulation of AR corona formation following magnetic flux emergence. We then present results from the application of the method to AIA observations of Active Region 11158, comparing the region's thermal structure on two successive solar rotations. Additionally, we show how the DEM inversion method can be adapted to simultaneously invert AIA and XRT data, and how supplementing AIA data with the latter improves the inversion result. The speed of the method allows for routine production of DEM maps, thus facilitating science studies that require tracking of the thermal structure of the solar corona in time and space.
Recently Blum, Katz and Waxman have claimed that the flux of high energy cosmic ray (CR) positrons near Earth that has been measured with the Alpha Magnetic Spectrometer (AMS) aboard the International Space Station can be produced in the collisions of Galactic CR protons and nuclei with the ambient matter in the Galactic interstellar medium (ISM). Their claim was based on an alleged "robust upper limit to the positron flux" which neglected the energy loss of e+'s in the ISM. Inclusion of this energy loss, however, yields a much smaller upper limit, which excludes secondary production in the ISM by the Galactic cosmic rays as the main origin of the CR e^+ flux above 10 GeV.
We present the first results of a survey for high redshift, z $\ge$ 6, quasars using izY multi-colour photometric observations from the Dark Energy Survey (DES). Here we report the discovery and spectroscopic confirmation of the $\rm z_{AB}, Y_{AB}$ = 20.2, 20.2 (M$_{1450}$ = $-$26.5) quasar DES J0454$-$4448 with an emission line redshift of z = 6.10$\pm$0.03 and a HI near zone size of 4.6 $\pm$ 1.7 Mpc.The quasar was selected as an i-band drop out with i$-$z = 2.46 and z$_{AB} < 21.5$ from an area of $\rm \sim$300 deg$^2$. It is the brightest of our 43 candidates and was identified for follow-up spectroscopically solely based on the DES i$-$z and z$-$Y colours. The quasar is detected by WISE and has $W1_{AB} = 19.68$. The discovery of one spectroscopically confirmed quasar with 5.7 $<$ z $<$ 6.5 and z$_{AB} \leq$ 20.2 is consistent with recent determinations of the luminosity function at z $\sim$ 6. DES when completed will have imaged $\rm \sim$5000 deg$^2$ to $Y_{AB}$ = 23.0 ($5\sigma$ point source) and we expect to discover $>$ 50-100 new quasars with z $>$ 6 including 3-10 with z $>$ 7 dramatically increasing the numbers of quasars currently known that are suitable for detailed studies including determination of the neutral HI fraction of the intergalactic medium (IGM) during the epoch of Hydrogen reionization.
We apply two methods to reconstruct the Hubble parameter $H(z)$ as a function of redshift from 15 measurements of the expansion rate obtained from age estimates of passively evolving galaxies. These reconstructions enable us to derive the luminosity distance to a certain redshift $z$, calibrate the light-curve fitting parameters accounting for the (unknown) intrinsic magnitude of type Ia supernova (SNe Ia) and construct cosmological model-independent Hubble diagrams of SNe Ia. In order to test the compatibility between the reconstructed functions of $H(z)$, we perform a statistical analysis considering the latest SNe Ia sample, the so-called JLA compilation. We find that, while one of the reconstructed functions leads to a value of the local Hubble parameter $H_0$ in excellent agreement with the one reported by the Planck collaboration, the other requires a higher value of $H_0$, which is consistent with recent measurements of this quantity from Cepheids and other local distance indicators.
We propose that observations of "hidden" magnetars in central compact objects can be used to probe crustal activity of neutron stars with large internal magnetic fields. Estimates based on calculations by Perna \& Pons (2011), Pons \& Rea (2012) and Kaminker et al. (2014) suggest that central compact objects, which are proposed to be "hidden" magnetars, must demonstrate flux variations on the time scale of months-years. However, the most prominent candidate for the "hidden" magnetars --- CXO J1852.6+0040 in Kes 79, -- shows constant (within error bars) flux. This can be interpreted by lower variable crustal activity than in typical magnetars. Alternatively, CXO J1852.6+0040 can be in a high state of variable activity during the whole period of observations. Then we consider the source 1E161348-5055 in RCW103 as another candidate. Employing a simple 2D-modeling we argue that properties of the source can be explained by the crustal activity of the magnetar type. Thus, this object may be supplemented for the three known candidates for the "hidden" magnetars among central compact objects discussed in literature.
New views on the diffuse interstellar bands are discussed. In particular results from DIB surveys and the study of near-infrared DIBs.
We report on GalevNB (Galev for N-Body simulations), an integrated software solution that provides N-body users direct access to the software package GALEV (GALaxy EVolutionary synthesis models). GalevNB is developed for the purpose of a direct comparison between N-body simulations and observations. It converts the fundamental stellar properties of N-body simulations, i.e., stellar mass, temperature, stellar luminosity and metallicity, into observational magnitudes for a variety of filters of widely used instruments/telescopes (HST, ESO, SDSS, 2MASS), and into spectra that span from far-UV (90 $\rm \AA$) to near-IR (160 $\rm \mu$m).
We analyze a high resolution spectropolarimetric dataset collected for the He-weak B3p IV star HR 2949. The Zeeman effect is visible in the circularly polarized component of numerous spectral lines. The longitudinal magnetic field varies between approximately $-650$ and $+150$ G. The polar strength of the surface magnetic dipole is calculated to be 2.4$^{+0.3}_{-0.2}$ kG. The star has strong overabundances of Fe-peak elements, along with extremely strong overabundances of rare-earth elements; however, He, Al, and S are underabundant. This implies that HR 2949 is a chemically peculiar star. Variability is seen in all photospheric lines, likely due to abundance patches as seen in many Ap/Bp stars. Longitudinal magnetic field variations measured from different spectral lines yield different results, likely a consequence of uneven sampling of the photospheric magnetic field by the abundance patches. Analysis of photometric and spectroscopic data for both HR 2949 and its companion star, HR 2948, suggests a revision of HR 2949's fundamental parameters: in particular, it is somewhat larger, hotter, and more luminous than previously believed. There is no evidence of optical or ultraviolet emission originating in HR 2949's magnetosphere, despite its moderately strong magnetic field and relatively rapid rotation; however, when calculated using theoretical and empirical boundaries on the initial rotational velocity, the spindown age is compatible with the stellar age. With the extensive phase coverage presented here, HR 2949 will make an excellent subject for Zeeman Doppler Imaging.
We present a new method of finding cosmic voids using tomographic maps of Ly{\alpha} forest flux. We identify cosmological voids with radii of 2 - 12 $h^{-1}$Mpc in a large N-body simulation at $z = 2.5$, and characterize the signal of the high-redshift voids in density and Ly{\alpha} forest flux. The void properties are similar to what has been found at lower redshifts, but they are smaller and have steeper radial density profiles. Similarly to what has been found for low-redshift voids, the radial velocity profiles have little scatter and agree very well with the linear theory prediction. We run the same void finder on an ideal Ly{\alpha} flux field and tomographic reconstructions at various spatial samplings. We compare the tomographic map void catalogs to the density void catalog and find good agreement even with modest-sized voids ($r > 6 \, h^{-1}$Mpc). Using our simple void-finding method, the configuration of the ongoing CLAMATO survey covering 1 deg$^2$ would provide a sample of about 100 high-redshift voids. We also provide void-finding forecasts for larger area surveys, and discuss how these void samples can be used to test modified gravity models, study high-redshift void galaxies, and to make an Alcock-Paczynski measurement. To aid future work in this area, we provide public access to our simulation products, catalogs, and sample tomographic flux maps.
In this proceedings contribution we briefly discuss about the consequences of the presence of Majorana dark matter in a dense neutron star environment focusing on a particularly interesting possible indirect effect, namely that of bubble nucleation. This is somewhat similar to current techniques developed for direct detection using bubble chamber or superheated droplet detectors.
We use the Wide Field Camera 3 onboard the Hubble Space Telescope to obtain deep, high-resolution photometry of the young (age ~ 300 Myr) star cluster NGC1856 in the Large Magellanic Cloud. We compare the observed colour-magnitude diagram (CMD), after having applied a correction for differential reddening, with Monte Carlo simulations of simple stellar populations (SSPs) of various ages. We find that the main sequence turn-off (MSTO) region is wider than that derived from the simulation of a single SSP. Using constraints based on the distribution of stars in the MSTO region and the red clump, we find that the CMD is best reproduced using a combination of two different SSPs with ages separated by 80 Myr (0.30 and 0.38 Gyr, respectively). However, we can not formally exclude that the width of the MSTO could be due to a range of stellar rotation velocities if the efficiency of rotational mixing is higher than typically assumed. Using a King-model fit to the surface number density profile in conjunction with dynamical evolution models, we determine the evolution of cluster mass and escape velocity from an age of 10 Myr to the present age, taking into account the possible effects of primordial mass segregation. We find that the cluster has an escape velocity Vesc ~ 17 km/s at an age of 10 Myr, and it remains high enough during a period of ~ 100 Myr to retain material ejected by slow winds of first-generation stars. Our results are consistent with the presence of an age spread in NGC1856, in contradiction to the results of Bastian & Silva-Villa (2013).
We calculated the time evolution of the free magnetic energy during the 2014-Mar-29 flare (SOL2014-03-29T17:48), the first X-class flare detected by IRIS. The free energy was calculated from the difference between the nonpotential field, constrained by the geometry of observed loop structures, and the potential field. We use AIA/SDO and IRIS images to delineate the geometry of coronal loops in EUV wavelengths, as well as to trace magnetic field directions in UV wavelengths in the chromosphere and transition region. We find an identical evolution of the free energy for both the coronal and chromospheric tracers, as well as agreement between AIA and IRIS results, with a peak free energy of $E_{free}(t_{peak}) \approx (45 \pm 2) \times 10^{30}$ erg, which decreases by an amount of $\Delta E_{free} \approx (29 \pm 3) \times 10^{30}$ erg during the flare decay phase. The consistency of free energies measured from different EUV and UV wavelengths for the first time here, demonstrates that vertical electric currents (manifested in form of helically twisted loops) can be detected and measured from both chromospheric and coronal tracers.
The dwarf spheroidal galaxies (dSph) of the Milky Way are among the most attractive targets for indirect searches of dark matter. In this work, we reconstruct the dark matter annihilation (J-factor) and decay profiles for the newly discovered dSph Reticulum~II. This is done using an optimized spherical Jeans analysis of kinematic data obtained from the Michigan/Magellan Fiber System (M2FS). We find Reticulum~II to have one of the highest J-factor when compared to the other Milky Way dSphs. We have also checked the robustness of this result against several ingredients of the analysis. Unless it suffers from tidal disruption or significant inflation of its velocity dispersion from binary stars, Reticulum~II may provide a unique window on dark matter particle properties.
We study a simple Standard Model (SM) extension, which includes three families of right-handed neutrinos with generic non-trivial flavor structure and an economic implementation of the invisible axion idea. We find that in some regions of the parameter space this model accounts for all experimentally confirmed pieces of evidence for physics beyond the SM: it explains neutrino masses (via the type-I see-saw mechanism), dark matter, baryon asymmetry (through leptogenesis), solve the strong CP problem and has a stable electroweak vacuum. The last property may allow us to identify the Higgs field with the inflaton.
We study the possibility of improving the constraints on the lifetime of gravitino dark matter in scenarios with bilinear R-parity violation by estimating the amount of cosmic-ray antideuterons that can be produced in gravitino decays. Taking into account all different sources of theoretical uncertainties, we find that the margin of improvement beyond the limits already set by cosmic-ray antiproton data are quite narrow and unachievable for the next generation of experiments. However, we also identify more promising energy ranges for future experiments.
The equation of state (EoS) of asymmetric nuclear matter at high densities is a key topic for the description of matter inside neutron stars. The determination of the properties of asymmetric nuclear matter, such as the symmetry energy ($a_{sym}$) and the slope of the symmetry energy ($L_0$) at saturation density, has been exaustively studied in order to better constrain the nuclear matter EoS. However, differently from symmetric matter properties that are reasonably constrained, the symmetry energy and its slope still large uncertainties in their experimental values. Regarding this subject, some studies point towards small values of the slope of the symmetry energy, while others suggest rather higher values. Such a lack of agreement raised a certain debate in the scientific community. In this paper, we aim to analyse the role of these properties on the behavior of asymmetric hyperonic matter. Using the formalism presented in Ref. (R.O. Gomes et al 2014}, which considers many-body forces contributions in the meson-baryon coupling, we calculate the EoS of asymmetric hyperonic matter and apply it to describe hyperonic matter and hyperon stars.
We will review the approach used for studying the conversion of a hadronic star into a quark star based on the assumption of a infinitely thin combustion zone and we will discuss why, in this scheme, the combustion stops before the whole hadronic star is converted.
A realistic mission scenario for the deflection of fictitious asteroid 2015PDC is investigated that makes use of the ion beam shepherd concept as primary deflection technique. The article deals with the design of a low thrust rendezvous trajectory to the asteroid, the estimation of the propagated covariance ellipsoid and the outcome of a slow-push deflection starting from three worst case scenarios (impacts in New Delhi, Dhaka and Teheran). Displacing the impact point towards very low populated areas, as opposed to full deflection, is found to be the simplest and most effective mitigation approach. Mission design, technical and political aspects are discussed.
It has recently been discovered that Einstein once attempted - and subsequently abandoned - a 'steady-state' model of the universe, i.e., a cosmic model in which the expanding universe remains essentially unchanged due to a continuous creation of matter from empty space. The discovery offers several new insights into Einstein's cosmology, from his view of the role of the cosmological constant to his attitude to the question of cosmic origins. More generally, Einstein's exploration of steady-state cosmology casts new light on his philosophical journey from a static, bounded cosmology to the dynamic, evolving universe, and is indicative of a pragmatic, empiricist approach to cosmology.
Intra Cluster Media (ICMs) located at galaxy clusters is in the state of hot, tenuous, magnetized, and highly ionized X-ray emitting plasmas. This overall collisionless, viscous, and conductive magnetohydrodynamic (MHD) turbulence in ICM is simulated using hyper magnetic diffusivity with weak background magnetic field. The result shows that fluctuating random plasma motion amplifies the magnetic field, which cascades toward the diffusivity scale passing through the viscous scale. The kinetic eddies in the subviscous scale are driven and constrained by the magnetic tension which is eventually balanced with the highly damping effect of the kinetic eddies. Simulation results show the saturated kinetic energy spectrum is $\sim k^{-3}$, deeper than that of the incompressible or compressible fluid. To explain this unusual field profile we set up two simultaneous differential equations for the kinetic and magnetic energy spectrum using an Eddy Damped Quasi Normal Markovianized (EDQNM) approximation. The analytic solution shows that the magnetic energy in addition to the viscous damping effect constrains the plasma motion leading to the power spectra: kinetic energy spectrum $E_V^k\sim k^{-3}$ and corresponding magnetic energy spectrum $E_M^k\sim k^{-1/2}$. Also the comparison of simulation results with different resolutions implies the role of small scale magnetic energy in the nonlocal energy transfer from kinetic to magnetic eddy.
We present details of numerical simulations of the gravitational radiation produced by a first order {thermal} phase transition in the early universe. We confirm that the dominant source of gravitational waves is sound waves generated by the expanding bubbles of the low-temperature phase. We demonstrate that the sound waves have a power spectrum with power-law form between the scales set by the average bubble separation (which sets the length scale of the fluid flow $L_\text{f}$) and the bubble wall width. The sound waves generate gravitational waves whose power spectrum also has a power-law form, at a rate proportional to $L_\text{f}$ and the square of the fluid kinetic energy density. We identify a dimensionless parameter $\tilde\Omega_\text{GW}$ characterising the efficiency of this "acoustic" gravitational wave production whose value is $8\pi\tilde\Omega_\text{GW} \simeq 0.8 \pm 0.1$ across all our simulations. We compare the acoustic gravitational waves with the standard prediction from the envelope approximation. Not only is the power spectrum steeper (apart from an initial transient) but the gravitational wave energy density is generically two orders of magnitude or more larger.
Using fully kinetic simulations, we study the x-line orientation of magnetic reconnection in an asymmetric configuration. A spatially localized perturbation is employed to induce a single x-line, that has sufficient freedom to choose its orientation in three-dimensional systems. The effect of ion to electron mass ratio is investigated, and the x-line appears to bisect the magnetic shear angle across the current sheet in the large mass ratio limit. The orientation can generally be deduced by scanning through corresponding 2D simulations to find the reconnection plane that maximizes the peak reconnection electric field. The deviation from the bisection angle in the lower mass ratio limit can be explained by the physics of tearing instability.
We consider a Lorentz-violating theory of gravity where the aether vector is taken to be nondynamical. This "ponderable aether theory" is almost the same as Einstein-aether theory (where the aether vector is dynamical), but involves additional integration constants arising due to the loss of initial value constraints. One of these produces an effective energy density for the aether fluid, similar to the appearance of dark matter in projectable Ho\v{r}ava gravity and the mimetic dark matter theory. Here we investigate the extent to which this energy density can reproduce the phenomenology of dark matter. Although it is indistinguishable from cold dark matter in homogeneous, isotropic cosmology, it encounters phenomenological problems in both spherically symmetric configurations and cosmological perturbations. Furthermore, inflationary considerations lead us to expect a tiny value for the ponderable aether energy density today unless a sourcing effect is added to the theory. The theory then effectively reduces to dynamical Einstein-aether theory, rendering moot the question of whether an aether must be dynamical in order to be consistent.
AMS-02 recently published its lepton spectrums measurement. The results show that the positron fraction no longer increases above $\sim$200 GeV. The aim of this work is to investigate the possibility that the excess of positron fraction is due to nearby pulsars. We consider nearby known pulsars from ATNF catalogue as the possible extra single sources of the high energy positrons. We find that the pulsars with age $T\simeq (0.45\sim4.5)\times10^{5}$ yr and distance $d<0.5$ kpc can explain the behavior of positron fraction of AMS-02 in the range of high energy. We show that each of the four pulsars --- Geminga, J1741-2054, Monogem and J0942-5552 --- is able to be a single source satisfying all considered physical requirements. We also discuss the possibility that these high energy $e^{\pm}$ are from multiple pulsars. The multiple pulsars contribution predicts a positron fraction with some structures at higher energies.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
We present the results of a new, non-parametric method to reconstruct the Galactic dark matter profile directly from observations. Using the latest kinematic data to track the total gravitational potential and the observed distribution of stars and gas to set the baryonic component, we infer the dark matter contribution to the circular velocity across the Galaxy. The radial derivative of this dynamical contribution is then estimated to extract the dark matter profile. The innovative feature of our approach is that it makes no assumption on the functional form nor shape of the profile, thus allowing for a clean determination with no theoretical bias. We illustrate the power of the method by constraining the spherical dark matter profile between 2.5 and 25 kpc away from the Galactic centre. The results show that the proposed method, free of widely used assumptions, can already be applied to pinpoint the dark matter distribution in the Milky Way with competitive accuracy, and paves the way for future developments.
We recently identified a population of low surface brightness objects in the field of the z=0.023 Coma cluster, using the Dragonfly Telephoto Array. Here we present Keck spectroscopy of one of the largest of these "ultra-diffuse galaxies" (UDGs), confirming that it is a member of the cluster. The galaxy has prominent absorption features, including the Ca II H+K lines and the G-band, and no detected emission lines. Its radial velocity of cz=6280 +- 120 km/s is within the 1 sigma velocity dispersion of the Coma cluster. The galaxy has an effective radius of 4.3 +- 0.3 kpc and a Sersic index of 0.89 +- 0.06, as measured from Keck imaging. We find no indications of tidal tails or other distortions, at least out to a radius of ~2 r_e. We show that UDGs are located in a previously sparsely populated region of the size - magnitude plane of quiescent stellar systems, as they are ~6 magnitudes fainter than normal early-type galaxies of the same size. It appears that the luminosity distribution of large quiescent galaxies is not continuous, although this could largely be due to selection effects. Dynamical measurements are needed to determine whether the dark matter halos of UDGs are similar to those of galaxies with the same luminosity or to those of galaxies with the same size.
We present high-redshift predictions of the star-formation-rate distribution function (SFR DF), UV luminosity function (UV LF), galactic stellar mass function (GSMF), and specific star-formation rates (sSFRs) of galaxies from the latest version of the Munich semi-analytic model L-Galaxies. We find a good fit to both the shape and normalisation of the SFR DF at $z=4-7$, apart from a slight under-prediction at the low SFR end at $z=4$. Likewise, we find a good fit to the faint number counts for the observed UV LF; at brighter magnitudes our predictions lie below the observations, increasingly so at higher redshifts. At all redshifts and magnitudes, the raw (unattenuated) number counts for the UV LF lie above the observations. Because of the good agreement with the SFR we interpret our under-prediction as an over-estimate of the amount of dust in the model for the brightest galaxies, especially at high-redshift. While the shape of our GSMF matches that of the observations, we lie between (conflicting) observations at $z=4-5$, and under-predict at $z=6-7$. The sSFRs of our model galaxies show the observed trend of increasing normalisation with redshift, but do not reproduce the observed mass dependence. Overall, we conclude that the latest version of L-Galaxies, which is tuned to match observations at $z\leq3$, does a fair job of reproducing the observed properties of galaxies at $z\geq4$. More work needs to be done on understanding observational bias at high-redshift, and upon the dust model, before strong conclusions can be drawn on how to interpret remaining discrepancies between the model and observations.
We report on the first measurement of the position-dependent correlation function from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 10 CMASS sample. This new observable measures the correlation between two-point functions of galaxy pairs within different subvolumes, $\hat{\xi}({\rm r},{\rm r}_L)$, where ${\rm r}_L$ is the location of a subvolume, and the corresponding mean overdensities, $\bar{\delta}({\rm r}_L)$. This correlation, which we call the "integrated three-point function", $i\zeta(r)=\langle\hat{\xi}({\rm r},{\rm r}_L)\bar{\delta}({\rm r}_L)\rangle$, measures a three-point function of two short- and one long-wavelength modes, and is generated by nonlinear gravitational evolution and possibly also by the physics of inflation. The $i\zeta(r)$ measured from the BOSS data lies within the scatter of those from the mock galaxy catalogs in redshift space, yielding a ten-percent-level determination of the amplitude of $i\zeta(r)$. The tree-level perturbation theory in redshift space predicts how this amplitude depends on the linear and quadratic nonlinear galaxy bias parameters ($b_1$ and $b_2$), as well as on the amplitude and linear growth rate of matter fluctuations ($\sigma_8$ and $f$). Combining $i\zeta(r)$ with the constraints on $b_1\sigma_8$ and $f\sigma_8$ from the global two-point correlation function and that on $\sigma_8$ from the weak lensing signal of BOSS galaxies, we measure $b_2 = 1.30 \pm 0.54$ (68% CL).
Integral field unit (IFU) data of the iconic Pillars of Creation in M16 are presented. The ionisation structure of the pillars was studied in great detail over almost the entire visible wavelength range, and maps of the relevant physical parameters, e.g. extinction, electron density, electron temperature, line-of-sight velocity of the ionised and neutral gas are shown. In agreement with previous authors, we find that the pillar tips are being ionised and photo-evaporated by the massive members of the nearby cluster NGC 6611. They display a stratified ionisation structure where the emission lines peak in a descending order according to their ionisation energies. The IFU data allowed us to analyse the kinematics of the photo-evaporative flow in terms of the stratified ionisation structure, and we find that, in agreement with simulations, the photo-evaporative flow is traced by a blueshift in the position-velocity profile. The gas kinematics and ionisation structure have allowed us to produce a sketch of the 3D geometry of the Pillars, positioning the pillars with respect to the ionising cluster stars. We use a novel method to detect a previously unknown bipolar outflow at the tip of the middle pillar and suggest that it has an embedded protostar as its driving source. Furthermore we identify a candidate outflow in the leftmost pillar. With the derived physical parameters and ionic abundances, we estimate a mass loss rate due to the photo-evaporative flow of 70 M$_{\odot}$ Myr$^{-1}$ which yields an expected lifetime of approximately 3 Myr.
We perform in-depth dynamical modelling of the luminous and dark matter (DM) content of the elliptical galaxy NGC 1407. Our strategy consists of solving the spherical Jeans equations for three independent dynamical tracers: stars, blue GCs and red GCs in a self-consistent manner. We adopt a maximum-likelihood Markov-Chain Monte Carlo fitting technique in the attempt to constrain the inner slope of the DM density profile (the cusp/core problem), and the stellar initial mass function (IMF) of the galaxy. We find the inner logarithmic slope of the DM density profiles to be $\gamma = 0.6\pm0.4$, which is consistent with either a DM cusp ($\gamma = 1$) or with a DM core $(\gamma = 0)$. Our findings are consistent with a Salpeter IMF, and marginally consistent with a Kroupa IMF. We infer tangential orbits for the blue GCs, and radial anisotropy for red GCs and stars. The modelling results are consistent with the virial mass--concentration relation predicted by $\Lambda$CDM simulations. The virial mass of NGC 1407 is $\log$ $M_{\rm vir} = 13.3 \pm 0.2 M_{\odot}$, whereas the stellar mass is $\log M_* = 11.8 \pm 0.1 M_{\odot}$. The overall uncertainties on the mass of NGC 1407 are only 5 per cent at the projected stellar effective radius. We attribute the disagreement between our results and previous X-ray results to the gas not being in hydrostatic equilibrium in the central regions of the galaxy. The halo of NGC 1407 is found be DM dominated, with a dynamical mass-to-light ratio of $M/L=260_{-100} ^{+174} M_{\odot}/L_{\odot, B}$. However, this value can be larger up to a factor of 3 depending on the assumed prior on the DM scale radius.
There is ongoing debate regarding the extent that environment affects galaxy size growth beyond z>1. To investigate the differences in star-forming and quiescent galaxy properties as a function of environment at z=2.1, we create a mass-complete sample of 59 cluster galaxies Spitler et al. (2012) and 478 field galaxies with log(M)>9 using photometric redshifts from the ZFOURGE survey. We compare the mass-size relation of field and cluster galaxies using measured galaxy semi-major axis half-light radii ($r_{1/2,maj}$) from CANDELS HST/F160W imaging. We find consistent mass normalized (log(M)=10.7) sizes for quiescent field galaxies ($r_{1/2,maj}=1.81\pm0.29$ kpc) and quiescent cluster galaxies ($r_{1/2,maj}=2.17\pm0.63$ kpc). The mass normalized size of star-forming cluster galaxies ($r_{1/2,maj}=4.00\pm0.26$ kpc ) is 12% larger (KS test $2.1\sigma$) than star-forming field galaxies ($r_{1/2,maj}=3.57\pm0.10$ kpc). From the mass-color relation we find that quiescent field galaxies with 9.7<log(M)<10.4 are slightly redder (KS test $3.6\sigma$) than quiescent cluster galaxies, while cluster and field quiescent galaxies with log(M)>10.4 have consistent colors. We find that star-forming cluster galaxies are on average 20% redder than star-forming field galaxies at all masses. Furthermore, we stack galaxy images to measure average radial color profiles as a function of mass. Negative color gradients are only present for massive star-forming field and cluster galaxies with log(M)>10.4, the remaining galaxy masses and types have flat profiles. Our results suggest given the observed differences in size and color of star-forming field and cluster galaxies, that the environment has begun to influence/accelerate their evolution. However, the lack of differences between field and cluster quiescent galaxies indicates that the environment has not begun to significantly influence their evolution at z~2.
It is now well-established that the elemental abundance patterns of stars holds key clues not only to their formation but also to the assembly histories of galaxies. One of the most exciting possibilities is the use of stellar abundance patterns as "chemical tags" to identify stars that were born in the same molecular cloud. In this paper we assess the prospects of chemical tagging as a function of several key underlying parameters. We build an observationally-motivated model for the star formation history (SFH), the gas and stellar mass distributions, and the radial size growth of the Milky Way through cosmic time. The multidimensional grid of parameters includes the fraction of stars that were born in-situ in the Solar annulus, the evolution and slope of the zero age cluster mass function (CMF), the survey geometry, number of stars in the survey, and the dimensionality of the chemical space. We show that in the fiducial case of $10^4$ distinct cells in chemical space and $10^5-10^6$ stars in the survey, one can expect to detect $\sim 10^2-10^3$ groups that are $\ge 5 \sigma$ overdensities in the chemical space. However, we find that even very large overdensities in chemical space do not guarantee that the overdensity is due to a single set of stars from a common birth cloud. In fact, for our fiducial model parameters, the typical $5 \sigma$ overdensity is comprised of stars from a wide range of clusters with the most dominant cluster contributing only 25% of the stars. [...] While recovering individual clusters through chemical tagging may prove challenging, we show, in agreement with previous work, that different CMFs imprint different degrees of clumpiness in chemical space. These differences provide the opportunity to statistically reconstruct the slope and high mass cutoff of CMF and its evolution through cosmic time.
X-ray and UV line emission in X-ray binaries can be accounted for by a hot corona. Such a corona forms through irradiation of the outer disk by radiation produced in the inner accretion flow. The same irradiation can produce a strong outflow from the disk at sufficiently large radii. Outflowing gas has been recently detected in several X-ray binaries via blue-shifted absorption lines. However, the causal connection between winds produced by irradiation and the blue-shifted absorption lines is problematic, particularly in the case of GRO J1655-40. Observations of this source imply wind densities about two orders of magnitude higher than theoretically predicted. This discrepancy does not mean that these `thermal disk-winds' cannot explain blue-shifted absorption in other systems, nor that they are unimportant as a sink of matter. Motivated by the inevitability of thermal disk-winds and wealth of data taken with current observatories such as Chandra, XMM-Newton and Suzaku, as well as the future AstroH mission, we decided to investigate the requirements to produce very dense winds. Using physical arguments, hydrodynamical simulations and absorption line calculations, we found that modification of the heating and cooling rates by a factor of a few results in an increase of the wind density of up to an order of magnitude and the wind velocity by a factor of about two. Therefore, the mass loss rate from the disk can be one, if not even two orders of magnitude higher than the accretion rate onto the central object. Such a high mass loss rate is expected to destabilize the disk and perhaps provides a mechanism for state change.
We present the infrared (IR) and X-ray properties of a sample of 33 mid-IR luminous quasars ($\nu$L(6 micron)>6x10$^{44}$ erg/s) at redshift z~1-3, identified through detailed spectral energy distribution analyses of distant star-forming galaxies, using the deepest IR data from Spitzer and Herschel in the GOODS-Herschel fields. The aim is to constrain the fraction of obscured, and Compton-thick (CT, N$_H$>1.5x10$^{24}$ cm$^{-2}$) quasars at the peak era of nuclear and star-formation activities. Despite being very bright in the mid-IR band, ~30% of these quasars are not detected in the extremely deep 2 Ms and 4 Ms Chandra X-ray data available in these fields. X-ray spectral analysis of the detected sources reveals that the majority (~67%) are obscured by column densities N$_H$>10$^{22}$ cm$^{-2}$; this fraction reaches ~80% when including the X-ray undetected sources (9 out of 33), which are likely to be the most heavily-obscured, CT quasars. We constrain the fraction of CT quasars in our sample to be ~21-45%, and their space density to be $\Phi$=(6.2$\pm$2.2)x10$^{-6}$ Mpc$^{-3}$. From the investigation of the quasar host galaxies in terms of star-formation rates (SFRs) and morphological distortions, as a sign of galaxy mergers/interactions, we do not find any direct relation between SFRs and quasar luminosity or X-ray obscuration. On the other hand, there is tentative evidence that the most heavily-obscured quasars have, on average, more disturbed morphologies than the unobscured/moderately-obscured quasar hosts, which preferentially live in undisturbed systems. However, the fraction of quasars with disturbed morphology amongst the whole sample is ~40%, suggesting that galaxy mergers are not the main fuelling mechanism of quasars at z~2.
This is the summary chapter of a review book on galaxy bulges. Bulge properties and formation histories are more varied than those of ellipticals. I emphasize two advances: 1 - "Classical bulges" are observationally indistinguishable from ellipticals, and like them, are thought to form by major galaxy mergers. "Disky pseudobulges" are diskier and more actively star-forming (except in S0s) than are ellipticals. Theys are products of the slow ("secular") evolution of galaxy disks: bars and other nonaxisymmetries move disk gas toward the center, where it starbursts and builds relatively flat, rapidly rotating components. This secular evolution is a new area of galaxy evolution work that complements hierarchical clustering. 2 - Disks of high-redshift galaxies are unstable to the formation of mass clumps that sink to the center and merge - an alternative channel for the formation of classical bulges. I review successes and unsolved problems in the formation of bulges+ellipticals and their coevolution (or not) with supermassive black holes. I present an observer's perspective on simulations of dark matter galaxy formation including baryons. I review how our picture of the quenching of star formation is becoming general and secure at redshifts z < 1. The biggest challenge is to produce realistic bulges+ellipticals and disks that overlap over a factor of 10**3 in mass but that differ from each other as observed over that whole range. Second, how does hierarchical clustering make so many giant, bulgeless galaxies in field but not cluster environments? I argue that we rely too much on AGN and star-formation feedback to solve these challenges.
We present an analysis of the X-ray spectrum and long-term variability of the nearby dwarf starburst galaxy Henize 2-10. Recent observations suggest that this galaxy hosts an actively accreting black hole with mass ~10^6 M_sun. The presence of an AGN in a low-mass starburst galaxy marks a new environment for active galactic nuclei (AGNs), with implications for the processes by which "seed" black holes may form in the early Universe. In this paper, we analyze four epochs of X-ray observations of Henize 2-10, to characterize the long-term behavior of its hard nuclear emission. We analyze observations with Chandra from 2001 and XMM-Newton from 2004 and 2011, as well as an earlier, less sensitive observation with ASCA from 1997. Based on detailed analysis of the source and background, we find that the hard (2-10 keV) flux of the putative AGN has decreased by approximately an order of magnitude between the 2001 Chandra observation and exposures with XMM-Newton in 2004 and 2011. The observed variability confirms that the emission is due to a single source. It is unlikely that the variable flux is due to a supernova or ultraluminous X-ray source, based on the observed long-term behavior of the X-ray and radio emission, while the observed X-ray variability is consistent with the behavior of well-studied AGNs.
This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold ones with spectral type T and Y, whose modeling is a current challenge. Silicate and iron clouds are believed to disappear at the photosphere at the L/T transition, but cloudless models fail to reproduce correctly the spectra of T dwarfs, advocating for the addition of more physics, e.g. other types of clouds or internal energy transport mechanisms. We use a one-dimensional (1D) radiative/convective equilibrium code ATMO to investigate this issue. This code includes both equilibrium and out-of-equilibrium chemistry and solves consistently the PT structure. Included opacity sources are H2-H2, H2-He, H2O, CO, CO2, CH4, NH3, K, Na, and TiO, VO if they are present in the atmosphere. We show that the spectra of Y dwarfs can be accurately reproduced with a cloudless model if vertical mixing and NH3 quenching are taken into account. T dwarf spectra still have some reddening in e.g. J - H compared to cloudless models. This reddening can be reproduced by slightly reducing the temperature gradient in the atmosphere. We propose that this reduction of the stabilizing temperature gradient in these layers, leading to cooler structures, is due to the onset of fingering convection, triggered by the destabilizing impact of condensation of very thin dust.
The evolution of the Universe between inflation and the onset of Big Bang Nucleosynthesis is difficult to probe and largely unconstrained. This ignorance profoundly limits our understanding of dark matter: we cannot calculate its thermal relic abundance without knowing when the Universe became radiation dominated. Fortunately, small-scale density perturbations provide a probe of the early Universe that could break this degeneracy. If dark matter is a thermal relic, density perturbations that enter the horizon during an early matter-dominated era grow linearly with the scale factor prior to reheating. The resulting abundance of substructure boosts the annihilation rate by several orders of magnitude, which can compensate for the smaller annihilation cross sections that are required to generate the observed dark matter density in these scenarios. In particular, thermal relics with masses less than a TeV that thermally and kinetically decouple prior to reheating may already be ruled out by Fermi-LAT observations of dwarf spheroidal galaxies. Although these constraints are subject to uncertainties regarding the internal structure of the microhalos that form from the enhanced perturbations, they open up the possibility of using gamma-ray observations to learn about the reheating of the Universe.
Recent results from the AMS-02 data have confirmed that the cosmic ray positron fraction increases with energy between 10 and 200GeV. This quantity should not exceed 50%, and it is hence expected that it will either converge towards 50% or fall. We study the possibility that future data may show the positron fraction dropping down abruptly to the level expected with only secondary production, and forecast the implications of such a feature in term of possible injection mechanisms that include both Dark Matter and pulsars.
Galaxies in Hickson Compact Group 91 (HCG 91) were observed with the WiFeS integral field spectrograph as part of our ongoing campaign targeting the ionized gas physics and kinematics inside star forming members of compact groups. Here, we report the discovery of HII regions with abundance and kinematic offsets in the otherwise unremarkable star forming spiral HCG 91c. The optical emission line analysis of this galaxy reveals that at least three HII regions harbor an oxygen abundance ~0.15 dex lower than expected from their immediate surroundings and from the abundance gradient present in the inner regions of HCG 91c. The same star forming regions are also associated with a small kinematic offset in the form of a lag of 5-10 km/s with respect to the local circular rotation of the gas. HI observations of HCG 91 from the Very Large Array and broadband optical images from Pan-STARRS suggest that HCG 91c is caught early in its interaction with the other members of HCG 91. We discuss different scenarios to explain the origin of the peculiar star forming regions detected with WiFeS, and show that evidence point towards infalling and collapsing extra-planar gas clouds at the disk-halo interface, possibly as a consequence of long-range gravitational perturbations of HCG 91c from the other group members. As such, HCG 91c provides evidence that some of the perturbations possibly associated with the early phase of galaxy evolution in compact groups impact the star forming disk locally, and on sub-kpc scales.
We measure the recent star formation history (SFH) across M31 using optical images taken with the \texit{Hubble Space Telescope} as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We fit the color-magnitude diagrams in ~9000 regions that are ~100 pc $\times$ 100 pc in projected size, covering a 0.5 square degree area (~380 kpc$^2$, deprojected) in the NE quadrant of M31. We show that the SFHs vary significantly on these small spatial scales but that there are also coherent galaxy-wide fluctuations in the SFH back to ~500 Myr, most notably in M31's 10-kpc star-forming ring. We find that the 10-kpc ring is at least 400 Myr old, showing ongoing star formation over the past ~500 Myr. This indicates the presence of molecular gas in the ring over at least 2 dynamical times at this radius. We also find that the ring's position is constant throughout this time, and is stationary at the level of 1 km/s, although there is evidence for broadening of the ring due to diffusion of stars into the disk. Based on existing models of M31's ring features, the lack of evolution in the ring's position makes a purely collisional ring origin highly unlikely. We find that the global SFR has been fairly constant over the last ~500 Myr, though it does show a small increase at 50 Myr that is 1.3 times the average SFR over the past 100 Myr. During the last ~500 Myr, ~60% of all SF occurs in the 10-kpc ring. Finally, we find that in the past 100 Myr, the average SFR over the PHAT survey area is $0.28\pm0.03$ M$_\odot$ yr$^{-1}$ with an average deprojected intensity of $7.3 \times 10^{-4}$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$, which yields a total SFR of ~0.7 M$_\odot$ yr$^{-1}$ when extrapolated to the entire area of M31's disk. This SFR is consistent with measurements from broadband estimates. [abridged]
By looking at the kinetic Sunyaev-Zeldovich effect (kSZ) in Planck nominal mission data, we present a significant detection of baryons participating in large-scale bulk flows around central galaxies (CGs) at redshift $z\approx 0.1$. We estimate the pairwise momentum of the kSZ temperature fluctuations at the positions of the CGC (Central Galaxy Catalogue) samples extracted from Sloan Digital Sky Survey (DR7) data. For the foreground-cleaned maps, we find $1.8$-$2.5\sigma$ detections of the kSZ signal, which are consistent with the kSZ evidence found in individual Planck raw frequency maps, although lower than found in the WMAP-9yr W band ($3.3\sigma$). We further reconstruct the peculiar velocity field from the CG density field, and compute for the first time the cross-correlation function between kSZ temperature fluctuations and estimates of CG radial peculiar velocities. This correlation function yields a $3.0$-$3.7$$\sigma$ detection of the peculiar motion of extended gas on Mpc scales, in flows correlated up to distances of 80-100 $h^{-1}$ Mpc. Both the pairwise momentum estimates and kSZ temperature-velocity field correlation find evidence for kSZ signatures out to apertures of 8 arcmin and beyond, corresponding to a physical radius of $> 1$ Mpc, more than twice the mean virial radius of halos. This is consistent with the predictions from hydro simulations that most of the baryons are outside the virialized halos. We fit a simple model, in which the temperature-velocity cross-correlation is proportional to the signal seen in a semi-analytic model built upon N-body simulations, and interpret the proportionality constant as an "effective" optical depth to Thomson scattering. We find $\tau_T=(1.4\pm0.5)\times 10^{-4}$; the simplest interpretation of this measurement is that much of the gas is in a diffuse phase, which contributes little signal to X-ray or thermal SZ observations.
The SkyMapper Southern Sky Survey is carrying out a search for the most metal-poor stars in the Galaxy. It identifies candidates by way of its unique filter set that allows for estimation of stellar atmospheric parameters. The set includes a narrow filter centered on the Ca II K 3933A line, enabling a robust estimate of stellar metallicity. Promising candidates are then confirmed with spectroscopy. We present the analysis of Magellan-MIKE high-resolution spectroscopy of 122 metal-poor stars found by SkyMapper in the first two years of commissioning observations. 41 stars have [Fe/H] <= -3.0. Nine have [Fe/H] <= -3.5, with three at [Fe/H] ~ -4. A 1D LTE abundance analysis of the elements Li, C, Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Co, Ni, Zn, Sr, Ba and Eu shows these stars have [X/Fe] ratios typical of other halo stars. One star with low [X/Fe] values appears to be "Fe-enhanced," while another star has an extremely large [Sr/Ba] ratio: >2. Only one other star is known to have a comparable value. Seven stars are "CEMP-no" stars ([C/Fe] > 0.7, [Ba/Fe] < 0). 21 stars exhibit mild r-process element enhancements (0.3 <=[Eu/Fe] < 1.0), while four stars have [Eu/Fe] >= 1.0. These results demonstrate the ability to identify extremely metal-poor stars from SkyMapper photometry, pointing to increased sample sizes and a better characterization of the metal-poor tail of the halo metallicity distribution function in the future.
Intermediate-mass young stellar objects (YSOs) provide a link to understand how feedback from shocks and UV radiation scales from low to high-mass star forming regions. Aims: Our aim is to analyze excitation of CO and H$_2$O in deeply-embedded intermediate-mass YSOs and compare with low-mass and high-mass YSOs. Methods: Herschel/PACS spectral maps are analyzed for 6 YSOs with bolometric luminosities of $L_\mathrm{bol}\sim10^2 - 10^3$ $L_\odot$. The maps cover spatial scales of $\sim 10^4$ AU in several CO and H$_2$O lines located in the $\sim55-210$ $\mu$m range. Results: Rotational diagrams of CO show two temperature components at $T_\mathrm{rot}\sim320$ K and $T_\mathrm{rot}\sim700-800$ K, comparable to low- and high-mass protostars probed at similar spatial scales. The diagrams for H$_2$O show a single component at $T_\mathrm{rot}\sim130$ K, as seen in low-mass protostars, and about $100$ K lower than in high-mass protostars. Since the uncertainties in $T_\mathrm{rot}$ are of the same order as the difference between the intermediate and high-mass protostars, we cannot conclude whether the change in rotational temperature occurs at a specific luminosity, or whether the change is more gradual from low- to high-mass YSOs. Conclusions: Molecular excitation in intermediate-mass protostars is comparable to the central $10^{3}$ AU of low-mass protostars and consistent within the uncertainties with the high-mass protostars probed at $3\cdot10^{3}$ AU scales, suggesting similar shock conditions in all those sources.
During the late stages of their evolution, Sun-like stars bring the products of nuclear burning to the surface. Most of the carbon in the Universe is believed to originate from stars with masses up to a few solar masses. Although there is a chemical dichotomy between oxygen-rich and carbon-rich evolved stars, the dredge-up itself has never been directly observed. In the last three decades, however, a few stars have been shown to display both carbon- and oxygen-rich material in their circumstellar envelopes. Two models have been proposed to explain this dual chemistry: one postulates that a recent dredge-up of carbon produced by nucleosynthesis inside the star during the Asymptotic Giant Branch changed the surface chemistry of the star. The other model postulates that oxygen-rich material exists in stable keplerian rotation around the central star. The two models make contradictory, testable, predictions on the location of the oxygen-rich material, either located further from the star than the carbon-rich gas, or very close to the star in a stable disk. Using the Faint Object InfraRed CAmera (FORCAST) instrument on board the Stratospheric Observatory for Infrared Astronomy (SOFIA) Telescope, we obtained images of the carbon-rich planetary nebula (PN) BD+30 3639 which trace both carbon-rich polycyclic aromatic hydrocarbons (PAHs) and oxygen-rich silicate dust. With the superior spectral coverage of SOFIA, and using a 3D photoionisation and dust radiative transfer model we prove that the O-rich material is distributed in a shell in the outer parts of the nebula, while the C-rich material is located in the inner parts of the nebula. These observations combined with the model, suggest a recent change in stellar surface composition for the double chemistry in this object. This is evidence for dredge-up occurring ~1000yr ago.
Chemical abundance of r-process elements in nearby dwarf spheroidal (dSph) galaxies is a powerful tool to probe the site of r-process since their small-mass scale can sort out individual events producing r-process elements. A merger of binary neutron stars is a promising candidate of this site. In faint, or less massive dSph galaxies such as the Draco, a few binary neutron star mergers are expected to have occurred at most over the whole past. We have measured chemical abundances including Eu and Ba of three red giants in the Draco dSph by Subaru/HDS observation. The Eu detection for one star with [Fe/H]=-1.45 confirms a broadly constant [Eu/H] of ~-1.3 for stars with [Fe/H]>-2. This feature is shared by other dSphs with similar masses, i.e., the Sculptor and the Carina, and suggests that neutron star merger is the origin of r-process elements in terms of its rarity. In addition, two very metal-poor stars with [Fe/H]=-2.12 and -2.51 are found to exhibit very low Eu abundances such as [Eu/H]<-2 with an implication of a sudden increase of Eu abundance by more than 0.7 dex at [Fe/H] ~ -2.2 in the Draco dSph. In addition, the detection of Ba abundances for these stars suggests that the r-process enrichment initiated no later than the time when only a few % of stars in the present-day Draco dSph was formed. Though an identification of the origin of an early Eu production inside the Draco dSph should be awaited until more abundance data of stars with [Fe/H]<-2 in the Draco as well as other faint dSphs become available, the implied early emergence of Eu production event might be reconciled with the presence of extremely metal-poor stars enriched by r-process elements in the Galactic halo.
This is the third of a series of papers of low X-ray luminosity galaxy clusters. In this work we present the weak lensing analysis of eight clusters, based on observations obtained with the Gemini Multi-Object Spectrograph in the $g'$, $r'$ and $i'$ passbands. For this purpose, we have developed a pipeline for the lensing analysis of ground-based images and we have performed tests applied to simulated data. We have determined the masses of seven galaxy clusters, six of them measured for the first time. For the four clusters with availably spectroscopic data, we find a general agreement between the velocity dispersions obtained via weak lensing assuming a Singular Isothermal Sphere profile, and those obtained from the redshift distribution of member galaxies. The correlation between our weak lensing mass determinations and the X-ray luminosities are suitably fitted by other observations of the $M-L_{X}$ relation and models.
We assess the fragmentation of strange quark matter in astrophysical events, showing that the application of statistical multifragmentation models suggests that most of the fragments (strangelets) should belong to an unstable baryon number $A$. While there are some caveats to be addressed, the flux of strangelets onto the Earth could be orders of magnitude lower than previous estimates, with negative prospects for the ongoing search experiments.
We present a search for laser emission coming from point sources in the vicinity of 2796 stars, including 1368 Kepler Objects of Interest (KOIs) that host one or more exoplanets. We search for extremely narrow emission lines in the wavelength region between 3640 and 7890 Angstroms using the Keck 10-meter telescope and spectroscopy with high resolution ($\lambda/\Delta \lambda$ = 60,000). Laser emission lines coming from non-natural sources are distinguished from natural astrophysical sources by being monochromatic and coming from an unresolved point in space. We search for laser emission located 2-7 arcsec from the 2796 target stars. The detectability of laser emission is limited by Poisson statistics of the photons and scattered light, yielding a detection threshold flux of approximately $10^{-2}$ photons $m^{-2} s^{-1}$ for typical Kepler stars and 1 photon $m^{-2} s^{-1}$ for solar-type stars within 100 light-years. Diffraction-limited lasers having a 10-meter aperture can be detected from 100 light-years away if their power exceeds 90 W, and from 1000 light-years away (Kepler planets), if their power exceeds 1 kW (from lasers located 60-200 AU, and 2000-7000 AU from the nearby and Kepler stars, respectively). We did not find any such laser emission coming from any of the 2796 target stars. We discuss the implications for the search for extraterrestrial intelligence (SETI).
The spiral arm tangencies are ideal lines of sight in which to determine the distribution of interstellar gas components in the spiral arms and study the influence of spiral density waves on the interarm gas in the Milky Way. We present a large scale (~15deg) position-velocity map of the Galactic plane in [CII] from l = 326.6 to 341.4deg observed with Herschel HIFI. We use [CII] l-v maps along with those for Hi and 12CO to derive the average spectral line intensity profiles over the longitudinal range of each tangency. Using the VLSR of the emission features, we locate the [CII], HI, and 12CO emissions along a cross cut of the spiral arm. In the spectral line profiles at the tangencies [CII] has two emission peaks, one associated with the compressed WIM and the other the molecular gas PDRs. When represented as a cut across the inner to outer edge of the spiral arm, the [CII]-WIM peak appears closest to the inner edge while 12CO and [CII] associated with molecular gas are at the outermost edge. HI has broader emission with an intermediate peak located nearer to that of 12CO. The velocity resolved spectral line data of the spiral arm tangencies unravel the internal structure in the arms locating the emission lanes within them. We interpret the excess [CII] near the tangent velocities as shock compression of the WIM induced by the spiral density waves and as the innermost edge of spiral arms. For the Norma and Perseus arms, we estimate widths of ~250 pc in [CII]-WIM and ~400 pc in 12CO and overall spiral arm widths of ~500 pc in [CII] and 12CO emissions. The electron densities in the WIM are ~ 0.5 cm^-3, about an order of magnitude higher than the average for the disk. The enhanced electron density in the WIM is a result of compression of the WIM by the spiral density wave potential.
Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright elliptical galaxies within a core of radius 10kpc. Such corrugation of the stellar distribution is very rare, and suggests recent formation by several simultaneous mergers. We map the distribution of associated dark matter, using new Hubble Space Telescope imaging and VLT/MUSE integral field spectroscopy of a gravitationally lensed system threaded through the cluster core. We find that each of the central galaxies retains a dark matter halo, but that (at least) one of these is spatially offset from its stars. The best-constrained offset is 1.62+/-0.48kpc, where the 68% confidence limit includes both statistical error and systematic biases in mass modelling. Such offsets are not seen in field galaxies, but are predicted during the long infall to a cluster, if dark matter self-interactions generate an extra drag force. With such a small physical separation, it is difficult to definitively rule out astrophysical effects operating exclusively in dense cluster core environments - but if interpreted solely as evidence for self-interacting dark matter, this offset implies a cross-section sigma/m=(1.7+/-0.7)x10^{-4}cm^2/g x (t/10^9yrs)^{-2}, where t is the infall duration.
Observation of the solar atmosphere reveals a wide range of motions, from small scale jets and spicules to global-scale coronal mass ejections. Identifying and characterizing these motions are essential to advancing our understanding the drivers of space weather. Both automated and visual identifications are currently used in identifying CMEs. To date, eruptions near the solar surface (which may be precursors to CMEs) have been identified primarily by visual inspection. Here we report on EruptionPatrol (EP): a software module that is designed to automatically identify eruptions from data collected by SDO/AIA. We describe the method underlying the module and compare its results to previous identifications found in the Heliophysics Event Knowledgebase. EP identifies eruptions events that are consistent with those found by human annotations, but in a significantly more consistent and quantitative manner. Eruptions are found to be distributed within 15Mm of the solar surface. They possess peak speeds ranging from 4 to 100 km/sec and display a power-law probability distribution over that range. These characteristics are consistent with previous observations of prominences.
We present photometry of the large scale environments of a sample of twelve broad line AGN with $0.06 < z < 0.37$ from deep images in the SDSS $u$, $g$, $r$, and $i$ filters taken with the 90Prime prime focus camera on the Steward Observatory Bok Telescope. We measure galaxy clustering around these AGN using two standard techniques: correlation amplitude (B$_{gq}$) and the two point correlation function. We find average correlation amplitudes for the 10 radio quiet objects in the sample equal to (9$\pm$18, 144$\pm$114, -39$\pm$56, 295$\pm$260) Mpc$^{1.77}$ in ($u$, $g$, $r$, $i$), all consistent with the expectation from galaxy clustering. Using a ratio of the galaxy-quasar cross-correlation function to the galaxy autocorrelation function, we calculate the relative bias of galaxies and AGN, $b_{gq}$. The bias in the $u$ band, $b_{gq}=3.08\pm0.51$ is larger compared to that calculated in the other bands, but it does not correlate with AGN luminosity, black hole mass, or AGN activity via the luminosity of the [OIII] emission line. Thus ongoing nuclear accretion activity is not reflected in the large scale environments from $\sim$10 h$^{-1}$ kpc to $\sim$0.5 h$^{-1}$ Mpc and may indicate a non-merger mode of AGN activity and/or a significant delay between galaxy mergers and nuclear activity in this sample of mostly radio quiet quasars.
Titan's equatorial regions are covered by eastward propagating linear dunes. This direction is opposite to mean surface winds simulated by Global Climate Models (GCMs), which are oriented westward at these latitudes, similar to trade winds on Earth. Different hypotheses have been proposed to address this apparent contradiction, involving Saturn's gravitational tides, large scale topography or wind statistics, but none of them can explain a global eastward dune propagation in the equatorial band. Here we analyse the impact of equinoctial tropical methane storms developing in the superrotating atmosphere (i.e. the eastward winds at high altitude) on Titan's dune orientation. Using mesoscale simulations of convective methane clouds with a GCM wind profile featuring superrotation, we show that Titan's storms should produce fast eastward gust fronts above the surface. Such gusts dominate the aeolian transport, allowing dunes to extend eastward. This analysis therefore suggests a coupling between superrotation, tropical methane storms and dune formation on Titan. Furthermore, together with GCM predictions and analogies to some terrestrial dune fields, this work provides a general framework explaining several major features of Titan's dunes: linear shape, eastward propagation and poleward divergence, and implies an equatorial origin of Titan's dune sand.
Re-acceleration of $\pi$'s and $\mu$'s modifies the flavor ratio at Earth (at astrophysical sources) of neutrinos produced by $\pi$ decay, $\nu_e:\nu_{\mu}:\nu_{\tau}$, from $1:1:1$ ($1:2:0$) to $1:1.8:1.8$ ($0:1:0$) at high energy, because $\pi$'s decay more than $\mu$'s during re-acceleration. The neutrino spectrum accompanies a flat excess, differently from the case of energy losses. With the flavor spectra, we can probe timescales of cosmic-ray acceleration and shock dynamics. We obtain general solutions of convection-diffusion equations and apply to gamma-ray bursts, which may have the flavor modification at around PeV -- EeV detectable by IceCube and next-generation experiments.
Nearby Type II (galaxy-spanning) Kardashev supercivilizations would have high mid-infrared (MIR) luminosities. We have used the Wide-field Infrared Survey Explorer (WISE) to survey ~$1 \times 10^5$ galaxies for extreme MIR emission, $10^3$ times more galaxies than the only previous such search. We have calibrated the WISE All-sky Catalog pipeline products to improve its photometry for extended sources. We present 563 extended sources with $|b| \ge 10$ and red MIR colors, having visually vetted them to remove artifacts. No galaxies in our sample host an alien civilization reprocessing more than 85% of its starlight into the MIR, and only 50 galaxies, including Arp 220, have MIR luminosities consistent with >50% reprocessing. Ninety of these (likely) extragalactic sources have little literature presence; in most cases they are likely barely resolved galaxies or pairs of galaxies undergoing large amounts of star formation. Five are new to science and deserve further study. The Be star 48 Librae sits within a MIR nebula, and we suggest that it may be creating dust. WISE, 2MASS, and Spitzer imagery shows that IRAS 04287+6444 is consistent with a previously unnoticed, heavily extinguished cluster of young stellar objects. We identify five "passive" (i.e. red) spiral galaxies with unusually high MIR and low NUV luminosity. We search a set of optically "dark" HI galaxies for MIR emission, and find none. These 90 poorly understood sources and five anomalous passive spirals deserve follow-up via both SETI and conventional astrophysics.
We present optical and infrared spectroscopy of V1309 Sco, an object which erupted in 2008 in a stellar-merger event. During the outburst, V1309 Sco displayed characteristics typical for red transients, a class of objects similar to V838 Mon. Our observations were obtained in 2009 and 2012, i.e. months and years after the eruption of V1309 Sco, and illustrate severe changes in the remnant, mainly in its circumstellar surroundings. In addition to atomic gas observed in earlier epochs, we identified molecular bands of TiO, VO, H$_2$O, ScO, AlO, and CrO. Infrared bands of CrO we analyze are first astronomical identification of those features. Over the whole period covered by our data, the remnant was associated with a cool ($\lesssim$1000 K) outflow with a terminal velocity of about 200 km\s. Signatures of warmer atomic gas, likely to be still dissipating the energy of the 2008 outburst, dramatically decreased their brightness between 2009 and 2012. Also, sometime before 2012 the source of optical continuum disappeared, likely owing to formation of new dust. The final stage captured by our spectra show an object remarkably similar to an older red transient, V4332 Sgr. In addition to providing a detailed view on the settling of the eruptive object, the observations presented here reinforce the conclusion that all the Galactic red transients are a manifestation of the same phenomenon, i.e. a stellar merger. Also, the late spectra of V1309 Sco suggest peculiarities in the chemical composition of the remnant, which yet need to be explored.
The masses of the black holes powering quasars represent a fundamental parameter of active galaxies. Estimates of quasar black hole masses using single-epoch spectra are quite uncertain, and require quantitative improvement. We recently identified a correction for C IV $\lambda$1549-based scaling relationships used to estimate quasar black hole masses that relies on the continuum-subtracted peak flux ratio of the ultraviolet emission-line blend Si IV + OIV] (the $\lambda$1400 feature) to that of C IV. This parameter correlates with the suite of associated quasar spectral properties collectively known as "Eigenvector 1" (EV1). Here we use a sample of 85 quasars with quasi-simultaneous optical-ultraviolet spectrophotometry to demonstrate how biases in the average EV1 properties can create systematic biases in C IV-based black hole mass scaling relationships. This effect results in nearly an order of magnitude moving from objects with small $<$ peak $\lambda$1400/C IV $>$, which have overestimated black hole masses, to objects with large $<$ peak $\lambda1400$/C IV $>$, which have underestimated values. We show that existing reverberation-mapped samples of quasars with ultraviolet spectra -- used to calibrate C IV-based scaling relationships -- have significant EV1 biases that result in predictions of black hole masses nearly 50\% too high for the average quasar. We offer corrections and suggestions to account for this bias.
We probe the possible anisotropy in the accelerated expanding Universe by using the JLA compilation of type-Ia supernovae. We constrain the amplitude and direction of anisotropy in the anisotropic cosmological models. For the dipole-modulated $\Lambda$CDM model, the anisotropic amplitude has an upper bound $D<1.04\times10^{-3}$ at the $68\%$ confidence level. Similar results are found in the dipole-modulated $w$CDM and CPL models. Our studies show that there are no significant evidence for the anisotropic expansion of the Universe. Thus the Universe is still well compatible with the isotropy.
We investigate the influence of photoevaporation of protoplanetary discs on the final distribution of exoplanets semi-major axis distances. We model giant planet migration in viscous discs affected by photoevaporation driven by either pure EUV or soft X-ray radiation (XEUV). We show that the final exoplanet distributions are strongly dependant on the choice of the photoevaporation model. In particular, we find that XEUV is more efficient than pure EUV radiation at parking planets at approximately 1-2 AU distance from their central star, hence roughly reproducing the observed peak in the exoplanets semi-major axis distributions. We note however that a more quantitative comparison with the observations is hindered by the oversimplified treatment of planetary accretion, which severely affects migration rates. For this reason, caution should be used when using these models to constrain details of disc clearing and/or migration from the observations. Nevertheless our results indicate that disc dispersal by photoevaporation may be the main driver of the features in the exoplanets semi-major axis distribution observed by recent surveys.
In this paper, the fractional projective Riccati expansion method is proposed to solve fractional differential equations. To illustrate the effectiveness of the method, we discuss the space-time fractional Burgers equation, the space-time fractional mKdV equation and time fractional biological population model. The solutions are expressed in terms of fractional hyperbolic functions. These solutions are useful to understand the mechanisms of the complicated nonlinear physical phenomena and fractional differential equations. Among these solutions, some are found for the first time. The fractal index for the obtained results is equal to one. Counter examples to compute the fractal index are introduced in appendix.
Rubano and Barrow have discussed the emergence of a dark energy, with late-time cosmic acceleration arising from a self-interacting homogeneous scalar field with a potential of hyperbolic power type. Here, we study the evolution of this scalar field potential back in the inflationary era. Using the hyperbolic power potential in the framework of inflation, we find that the main slow-roll parameters, like the scalar spectral index, the running of the spectral index and the tensor-to-scalar fluctuation ratio can be computed analytically. Finally, in order to test the viability of this hyperbolic scalar field model at the early stages of the Universe, we compare the predictions of that model against the latest observational data, namely Planck 2015.
The formation and dynamics of coronal rain are currently not fully understood. Coronal rain is the fall of cool and dense blobs formed by thermal instability in the solar corona towards the solar surface with acceleration smaller than gravitational free fall. We aim to study the observational evidence of the formation of coronal rain and to trace the detailed dynamics of individual blobs. We used time series of the 171 \AA\, and 304 \AA\, spectral lines obtained by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) above active region AR 11420 on February 22, 2012. Observations show that a coronal loop disappeared in the 171 \AA\ channel and appeared in the 304 \AA\ line$\text{}\text{}$ more than one hour later, which indicates a rapid cooling of the coronal loop from 1 MK to 0.05 MK. An energy estimation shows that the radiation is higher than the heat input, which indicates so-called catastrophic cooling. The cooling was accompanied by the formation of coronal rain in the form of falling cold plasma. We studied two different sequences of falling blobs. The first sequence includes three different blobs. The mean velocities of the blobs were estimated to be 50 km s$^{-1}$, 60 km s$^{-1}$ and 40 km s$^{-1}$. A polynomial fit shows the different values of the acceleration for different blobs, which are lower than free-fall in the solar corona. The first and second blob move along the same path, but with and without acceleration, respectively. We performed simple numerical simulations for two consecutive blobs, which show that the second blob moves in a medium that is modified by the passage of the first blob. Therefore, the second blob has a relatively high speed and no acceleration, as is shown by observations. The second sequence includes two different blobs with mean velocities of 100 km s$^{-1}$ and 90 km s$^{-1}$, respectively.
The question of whether ultra-long GRBs form a population different from that of "regular" long GRBs has been much debated recently and during the conference. We discuss here the data and the evidence that lead to the conclusion that indeed ultra-long GRBs form a different class of high energy transients. The sample of ultra-long GRBs is still poor and the discussion on their origin remain opens, though they might be the signature of PopIII stars. We urge that the design of new instrumentation, such as the SVOM satellite, takes into account the need for the detection of distant ultra-long GRBs.
The Hard X-ray Modulation Telescope (HXMT) will perform an all-sky survey in hard X-ray band as well as deep imaging of a series of small sky regions. We expect various compact objects to be detected in these imaging observations. Point source detection performance of HXMT imaging observation depends not only on the instrument but also on its data analysis since images are reconstructed from HXMT observed data with numeric methods. Denoising technique plays an import part in HXMT imaging data analysis pipeline alongside with demodulation and source detection. In this paper we have implemented several methods for denoising HXMT data and evaluated the point source detection performances in terms of sensitivities and location accuracies. The results show that direct demodulation with 1-fold cross correlation should be the default reconstruction and regularization methods, although both sensitivity and location accuracy could be further imporved by selecting and tuning numerical methods in data analysis of HXMT imaging observations.
We present a new way of describing the flares from Sgr A* with a self-consistent calculation of the particle distribution. All relevant radiative processes are taken into account in the evolution of the electron distribution and resulting spectrum. We present spectral modelling for new X-ray flares observed by NuSTAR, together with older observations in different wavelengths, and discuss the changes in plasma parameters to produce a flare. We show that under certain conditions, the real particle distribution can differ significantly from standard distributions assumed in most studies. We conclude that the flares are likely generated by magnetized plasma consistent with our understanding of the accretion flow. Including non-thermal acceleration, injection, escape, and cooling losses produces a spectrum with a break between the infrared and the X-ray, allowing a better simultaneous description of the different wavelengths. We favour the non-thermal synchrotron interpretation, assuming the infrared flare spectrum used is representative. We also consider the effects on Sgr A*s quiescent spectrum in the case of a density increase due to the G2 encounter with Sgr A*.
We present the abundance analysis for a sample of 7 red giant branch stars in the metal-poor globular cluster NGC 4372 based on UVES spectra acquired as part of the Gaia-ESO Survey. This is the first extensive study of this cluster from high resolution spectroscopy. We derive abundances of O, Na, Mg, Al, Si, Ca, Sc, Ti, Fe, Cr, Ni, Y, Ba, and La. We find a metallicity of [Fe/H] = -2.19 $\pm$ 0.03 and find no evidence for a metallicity spread. This metallicity makes NGC 4372 one of the most metal-poor galactic globular clusters. We also find an {\alpha}-enhancement typical of halo globular clusters at this metallicity. Significant spreads are observed in the abundances of light elements. In particular we find a Na-O anti-correlation. Abundances of O are relatively high compared with other globular clusters. This could indicate that NGC 4372 was formed in an environment with high O for its metallicity. A Mg-Al spread is also present which spans a range of more than 0.5 dex in Al abundances. Na is correlated with Al and Mg abundances at a lower significance level. This pattern suggests that the Mg-Al burning cycle is active. This behavior can also be seen in giant stars of other massive, metal-poor clusters. A relation between light and heavy s-process elements has been identified.
We consider a two-parametric family of radially anisotropic models with
non-singular density distribution in the centre. If highly eccentric orbits are
locked near the centre, the characteristic growth rate of the instability is
much less than the Jeans and dynamic frequencies of the stars (slow modes). The
instability occurs only for even spherical harmonics and the perturbations are
purely growing (aperiodic). On the contrary, if all orbits nearly reach the
outer radius of the sphere, both even and odd harmonics are unstable. Unstable
odd modes oscillate having characteristic frequencies of the order of the
dynamical frequencies (fast modes). Unstable even harmonics contain a single
aperiodic mode and several oscillatory modes, the aperiodic mode being the most
unstable.
The question of the nature of the radial orbit instability (ROI) is
revisited. Two main interpretations of ROI were suggested in the literature.
The first one refers to the classical Jeans instability associated with the
lack of velocity dispersion of stars in the transverse direction. The second
one refers to Lynden-Bell's orbital approach to bar formation in disc galaxies,
which implies slowness and bi-symmetry of the perturbation. Oscillatory modes,
odd spherical harmonics modes, and non-slow modes found in one of the models
show that the orbital interpretation is not the only possible.
QUIJOTE (Q-U-I JOint TEnerife) is a new polarimeter aimed to characterize the polarization of the Cosmic Microwave Background and other Galactic and extragalactic signals at medium and large angular scales in the frequency range 10-40 GHz. The multi-frequency (10-20~GHz) instrument, mounted on the first QUIJOTE telescope, saw first light on November 2012 from the Teide Observatory (2400~m a.s.l). During 2014 the second telescope has been installed at this observatory. A second instrument at 30~GHz will be ready for commissioning at this telescope during summer 2015, and a third additional instrument at 40~GHz is now being developed. These instruments will have nominal sensitivities to detect the B-mode polarization due to the primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r=0.05.
We present a simple approach for obtaining robust values of astrophysical parameters from the observed colour-magnitude diagrams (CMDs) of star clusters. The basic inputs are the Hess diagram built with the photometric measurements of a star cluster and a set of isochrones covering wide ranges of age and metallicity. In short, each isochrone is shifted in apparent distance modulus and colour excess until it crosses over the maximum possible Hess density. Repeating this step for all available isochrones leads to the construction of the solution map, in which the optimum values of age and metallicity - as well as foreground/background reddening and distance from the Sun - can be searched for. Controlled tests with simulated CMDs show that the approach is efficient in recovering the input values. We apply the approach to the open clusters M\,67, NGC\,6791, and NGC\,2635, which are characterised by different ages, metallicities and distances from the Sun.
In this work we selected one particular fibril from a high resolution solar chromosphere observation from the Dutch Open Telescope, and tried to obtain a broad picture of the intricate mechanism that might be incorporated in the multiple layer of the Solar atmosphere in high cadence multi-wavelength observation. We analyzed the changingvfibril patter using multi-wavelength tomography, which consists of both H$\alpha$ line center \& the blue wing, Doppler-signal, Ca II H, and the G-band. We have found that the intermittent ejected material through fibril from Doppler images has clearly shown oscillation mode, as seen in the H$\alpha$ blue wing. The oscillations in the umbrae and penumbrae magnetic field lines that are above the sunspot cause a broadening and forms the area like a ring shape from 3 to 15-minute oscillations as function of height. These made a distinct boundary of umbrae and penumbrae which suggest the comb structure, and indicate that the oscillations could propagate along the inclined magnetic flux tubes from below. The 3-minute strongly appeared in the broadly inclined penumbrae magnetic filed lines and gave the clear light-bridge. The well known 5-minute was dominated in the umbrae-penumbrae region boundary, the long 7-minute one was transparent in the H$\alpha$ blue wing, but this was the same with 10 and 15-minute, it was concentrated in the inner-penumbrae, as seen in the H$\alpha$ line center. From these findings we propose a picture on the role of fibril as the fabric of interaction between the layers, also the related activites around the active region under investigation.
With the high tempo-spatial \emph{Interface Region Imaging Spectrograph} 1330 {\AA} images, we find that many bright structures are rooted in the light bridge of NOAA 12192, forming a \emph{light wall}. The light wall is brighter than the surrounding areas, and the wall top is much brighter than the wall body. The New Vacuum Solar Telescope H$\alpha$ and the \emph{Solar Dynamics Observatory} 171 {\AA} and 131 {\AA} images are also used to study the light wall properties. In 1330 {\AA}, 171 {\AA}, and 131 {\AA}, the top of the wall has a higher emission, while in the H$\alpha$ line, the wall top emission is very low. The wall body corresponds to bright areas in 1330 {\AA} and dark areas in the other lines. The top of the light wall moves upward and downward successively, performing oscillations in height. The deprojected mean height, amplitude, oscillation velocity, and the dominant period are determined to be 3.6 Mm, 0.9 Mm, 15.4 km s$^{-1}$, and 3.9 min, respectively. We interpret the oscillations of the light wall as the leakage of \emph{p}-modes from below the photosphere. The constant brightness enhancement of the wall top implies the existence of some kind of atmospheric heating, e.g., via the persistent small-scale reconnection or the magneto-acoustic waves. In another series of 1330 {\AA} images, we find that the wall top in the upward motion phase is significantly brighter than in the downward phase. This kind of oscillations may be powered by the energy released due to intermittent impulsive magnetic reconnection.
Over the next decade, observations conducted with ALMA and the SKA will reveal the process of mass assembly and accretion onto young stars and will be revolutionary for studies of star formation. Here we summarise the capabilities of ALMA and discuss recent results from its early science observations. We then review infrared and radio variability observations of both young low-mass and high-mass stars. A time domain SKA radio continuum survey of star forming regions is then outlined. This survey will produce radio light-curves for hundreds of young sources, providing for the first time a systematic survey of radio variability across the full range of stellar masses. These light-curves will probe the magnetospheric interactions of young binary systems, the origins of outflows, trace episodic accretion on the central sources and potentially constrain the rotation rates of embedded sources.
ISO-Oph-50 is a young low-mass object in the ~Myr old Ophiuchus star forming region undergoing dramatic changes in its optical/near/mid-infrared brightness by 2-4 mag. We present new multi-band photometry and near-infrared spectra, combined with a synopsis of the existing literature data. Based on the spectroscopy, the source is confirmed as a mid M dwarf, with evidence for ongoing accretion. The near-infrared lightcurves show large-scale variations, with 2-4 mag amplitude in the bands IJHK, with the object generally being bluer when faint. Near its brightest state, the object shows colour changes consistent with variable extinction of dAV~7 mag. High-cadence monitoring at 3.6mu reveals quasi-periodic variations with a typical timescale of 1-2 weeks. The best explanation for these characteristics is a low-mass star seen through circumstellar matter, whose complex variability is caused by changing inhomogeneities in the inner parts of the disk. When faint, the direct stellar emission is blocked; the near-infrared radiation is dominated by scattered light. When bright, the emission is consistent with a photosphere strongly reddened by circumstellar dust. Based on the available constraints, the inhomogeneities have to be located at or beyond ~0.1 AU distance from the star. If this scenario turns out to be correct, a major portion of the inner disk has to be clumpy, structured, and/or in turmoil. In its observational characteristics, this object resembles other types of YSOs with variability caused in the inner disk. Compared to other objects, however, ISO-Oph-50 is clearly an extreme case, given the large amplitude of the brightness and colour changes combined with the erratic behaviour. ISO-Oph-50 has been near its brightest state since 2013; further monitoring is highly encouraged.
In these lectures we first concentrate on the cosmological problems which, hopefully, have to do with the new physics to be probed at the LHC: the nature and origin of dark matter and generation of matter-antimatter asymmetry. We give several examples showing the LHC cosmological potential. These are WIMPs as cold dark matter, gravitinos as warm dark matter, and electroweak baryogenesis as a mechanism for generating matter-antimatter asymmetry. In the remaining part of the lectures we discuss the cosmological perturbations as a tool for studying the epoch preceeding the conventional hot stage of the cosmological evolution.
Burn-UD is a hydrodynamic combustion code used to model the phase transition of hadronic to quark matter with particular application to the interior of neutron stars. Burn-UD models the flame micro-physics for different equations of state (EoS) on both sides of the interface, i.e. for both the ash (up-down-strange quark phase) and the fuel (up-down quark phase). It also allows the user to explore strange quark seeding produced by different processes including DM annihilation inside neutron stars. The simulations provide a physical window to diagnose whether the combustion process will simmer quietly and slowly, lead to a transition from deflagration to detonation or a (quark) core-collapse explosion. Such an energetic phase transition (a Quark-Nova) would have consequences in high-energy astrophysics and could aid in our understanding of many still enigmatic astrophysical transients. Furthermore, having a precise understanding of the phase transition dynamics for different EoSs could aid further in constraining the nature of the non-perturbative regimes of QCD in general. We hope that Burn-UD will evolve into a platform/software to be used and shared by the QCD community exploring the phases of Quark Matter and astrophysicists working on Compact Stars.
The photometry profile of the integrated Sachs-Wolfe (ISW) effect recently obtained by the Planck consortium by stacking patches of Cosmic Microwave Background (CMB) sky maps around a large number of cosmic voids, contains a cold ring at about half the void's effective radius surrounded by a hot ring near the void's boundary. The source of the temperature structure is assumed to be the ISW effect but the exact cause of the ringed structure is not currently well understood, particularly the outer hot ring. Numerical simulations have suggested that hot/cold ring structures can be produced by motions associated with nonlinear growths of cosmic structures whose gravitational potentials produce the ISW effect. We have recently developed the embedded lens theory and the Fermat potential formalism which can be used to model the ISW effect caused by intervening individual lens inhomogeneities evolving arbitrarily. This theory only requires knowledge of the void's projected mass profile as a function of the passing CMB photons' impact radius and the rate of change of that mass distribution at passage. We present two simple embedded void lens models with evolving mass densities and investigate the ISW effect caused by these lenses. Both models posses expanding mass shells which produce hot rings around central cold regions, consistent with the recent observations. By adding a small over-density at the void's center we can produce the slight positive temperature excess hinted at in Planck's photometric results. We conclude that the embedded lens theory and the Fermat potential formalism is well suited for modeling the ISW effect.
We present the first photometric observations of trans-Neptunian objects (TNOs) taken with the Kepler space telescope, obtained in the course of the K2 ecliptic survey. Two faint objects have been monitored in specifically designed pixel masks that were centered on the stationary points of the objects, when their daily motion was the slowest. In the design of the experiment, only the apparent path of these objects were retrieved from the detectors, i.e. the costs in terms of Kepler pixels were minimized. Because of the faintness of the targets we employ specific reduction techniques and co-added images. We measure rotational periods and amplitudes in the unfiltered Kepler band as follows: for (278331) 2007 JJ43 and 2002 GV31 we get P_rot=12.097 h and P_rot=29.2 h while 0.10 and 0.35 mag for the total amplitudes, respectively. Future space missions, like TESS and PLATO are not well suited to this kind of observations. Therefore, we encourage to include the brightest TNOs around their stationary points in each observing campaign to exploit this unique capability of the K2 Mission -- and therefore to provide unbiased rotational, shape and albedo characteristics of many objects.
We study the spot distribution on the surface of LQ~Hya during the observing seasons October 1998 -- November 2002. We look for persistent active longitudes, trends in the level of spot activity and compare to photometric data, specifically to the derived time epochs of the lightcurve minima. We apply the Doppler Imaging technique on photospheric spectral lines using an inversion code to retrieve images of the surface temperature. We present new temperature maps using multiple spectral lines for a total of 7 seasons. We calculate spot coverage fractions from each map, and as a result we find a general trend that is interpreted to be an indication of a spot cycle. There is a minimum during the observing season of March 1999. After this the activity increases until November 2000, followed by a general decrease in activity again. We find no evidence for active longitudes persisting over multiple observing seasons. The spot activity appears to be concentrated to two latitude regions. The high latitude spots are particularly strong when the spot coverage is at a maximum. Using the currently accepted rotation period, we find spot structures to show a trend in the phase-time plot, indicative of a need for a longer period. We conclude that the long-term activity of LQ~Hya is more chaotic than that of some magnetically active binary stars analyzed with similar methods, but still with clear indications of an activity cycle.
The Zeldovich approximation, 1st order Lagrangian perturbation theory, provides a good description of the clustering of matter and galaxies on large scales. The acoustic feature in the large-scale correlation function of galaxies imprinted by sound waves in the early Universe has been successfully used as a `standard ruler' to constrain the expansion history of the Universe. The standard ruler can be improved if a process known as density field reconstruction is employed. In this paper we develop the Zeldovich formalism to compute the correlation function of biased tracers in both real- and redshift-space using the simplest reconstruction algorithm with a Gaussian kernel and compare to N-body simulations. The model qualitatively describes the effects of reconstruction on the simulations, though its quantitative success depends upon how redshift-space distortions are handled in the reconstruction algorithm.
Simulations of binary neutron stars have seen great advances in terms of physical detail and numerical quality. However, the spin of the neutron stars, one of the simplest global parameters of binaries, remains mostly unstudied. We present the first, fully nonlinear general relativistic dynamical evolutions of the last three orbits for constraint satisfying initial data of spinning neutron star binaries, with astrophysically realistic spins aligned and anti-aligned to the orbital angular momentum. The initial data is computed with the constant rotational velocity approach. The dynamics of the systems is analyzed in terms of gauge-invariant binding energy vs. orbital angular momentum curves. By comparing to a binary black hole configuration we can estimate the different tidal and spin contributions to the binding energy for the first time. First results on the gravitational wave forms are presented. The phase evolution during the orbital motion is significantly affected by spin-orbit interactions, leading to delayed or early mergers. Furthermore, a frequency shift in the main emission mode of the hyper massive neutron star is observed. Our results suggest that a detailed modeling of merger waveforms requires the inclusion of spin, even for the moderate magnitudes observed in binary neutron star systems.
Binary neutron star mergers are studied using nonlinear 3+1 numerical relativity simulations and the analytical effective-one-body (EOB) model. The EOB model predicts quasiuniversal relations between the mass-rescaled gravitational wave frequency and the binding energy at the moment of merger, and certain dimensionless binary tidal coupling constants depending on the stars Love numbers, compactnesses and the binary mass ratio. These relations are quasiuniversal in the sense that, for a given value of the tidal coupling constant, they depend significantly neither on the equation of state nor on the mass ratio, though they do depend on stars spins. The spin dependence is approximately linear for small spins aligned with the orbital angular momentum. The quasiuniversality is a property of the conservative dynamics; nontrivial relations emerge as the binary interaction becomes tidally dominated. This analytical prediction is qualitatively consistent with new, multi-orbit numerical relativity results for the relevant case of equal-mass irrotational binaries. Universal relations are thus expected to characterize neutron star mergers dynamics. In the context of gravitational wave astronomy, these universal relations may be used to constrain the neutron star equation of state using waveforms that model the merger accurately.
The concept of magnetic connections is extended to non-ideal relativistic magnetohydrodynamical plasmas. Adopting a general set of equations for relativistic magnetohydrodynamics including thermal-inertial, thermal electromotive, Hall and current-inertia effects, we derive a new covariant connection equation showing the existence of generalized magnetofluid connections that are preserved during the dissipationless plasma dynamics. These connections are intimately linked to a general antisymmetric tensor that unifies the electromagnetic and fluid fields, allowing the extension of the magnetic connection notion to a much broader concept.
We present a frequentist analysis of the parameter space of the pMSSM10, in which the following 10 soft SUSY-breaking parameters are specified independently at the mean scalar top mass scale Msusy = Sqrt[M_stop1 M_stop2]: the gaugino masses M_{1,2,3}, the 1st-and 2nd-generation squark masses M_squ1 = M_squ2, the third-generation squark mass M_squ3, a common slepton mass M_slep and a common trilinear mixing parameter A, the Higgs mixing parameter mu, the pseudoscalar Higgs mass M_A and tan beta. We use the MultiNest sampling algorithm with 1.2 x 10^9 points to sample the pMSSM10 parameter space. A dedicated study shows that the sensitivities to strongly-interacting SUSY masses of ATLAS and CMS searches for jets, leptons + MET signals depend only weakly on many of the other pMSSM10 parameters. With the aid of the Atom and Scorpion codes, we also implement the LHC searches for EW-interacting sparticles and light stops, so as to confront the pMSSM10 parameter space with all relevant SUSY searches. In addition, our analysis includes Higgs mass and rate measurements using the HiggsSignals code, SUSY Higgs exclusion bounds, the measurements B-physics observables, EW precision observables, the CDM density and searches for spin-independent DM scattering. We show that the pMSSM10 is able to provide a SUSY interpretation of (g-2)_mu, unlike the CMSSM, NUHM1 and NUHM2. As a result, we find (omitting Higgs rates) that the minimum chi^2/dof = 20.5/18 in the pMSSM10, corresponding to a chi^2 probability of 30.8 %, to be compared with chi^2/dof = 32.8/24 (31.1/23) (30.3/22) in the CMSSM (NUHM1) (NUHM2). We display 1-dimensional likelihood functions for SUSY masses, and show that they may be significantly lighter in the pMSSM10 than in the CMSSM, NUHM1 and NUHM2. We discuss the discovery potential of future LHC runs, e+e- colliders and direct detection experiments.
We explore some particle physics implications of the growing evidence for a helical primordial magnetic field (PMF). From the interactions of magnetic monopoles and the PMF, we derive an upper bound on the monopole number density, $n(t_0) < 1 \times 10^{-20}~{\rm cm}^{-3}$, which is a "primordial" analog of the Parker bound for the survival of galactic magnetic fields. Our bound is weaker than existing constraints, but it is derived under independent assumptions. We also show how improved measurements of the PMF at different redshifts can lead to further constraints on magnetic monopoles. Axions interact with the PMF due to the $\varphi {\bf E}\cdot {\bf B}$ interaction. Including the effects of the cosmological plasma, we find that the helicity of the PMF is a source for the axion field. Although the magnitude of the source is small for the PMF, it could potentially be of interest in astrophysical environments. Finally we apply constraints on the neutrino magnetic dipole moment that arise from requiring successful big bang nucleosynthesis in the presence of a PMF, and using the suggested strength $\sim 10^{-14}~{\rm G}$ and coherence length $\sim 10~{\rm Mpc}$ we find $\mu_\nu \lesssim 10^{-16} \mu_B$.
We investigate how long wavelength inflationary fluctuations can cause the background field to deviate from classical dynamics. For generic potentials, we show that, in the Hartree approximation, the long wavelength dynamics can be encapsulated by a two-field model operating in an effective potential. The latter is given by a simple Gaussian integral transformation of the original inflationary potential. We use this new expression to study backreaction effects in quadratic, hilltop, flattened, and axion monodromy potentials. We find that the net result of the altered dynamics is to slightly modify the spectral tilt, drastically decrease the tensor-to-scalar ratio, and to effectively smooth over any features of the potential, with the size of these deviations set by the initial value of power in large scale modes and the shape of the potential during the entire evolution.
In this paper, we considered an inflationary model that effectively behaves as a modified Chaplygin gas in the context of quintessence cosmology. We reconstructed the inflaton potential bottom-up and using the recent observational data we fixed the free parameters of the model. We showed that the modified Chaplygin gas inspired model is suitable for both the early and the late time acceleration but has shortcomings between the two periods.
We test the physical relevance of the full and truncated versions of the Israel-Stewart theory of irreversible thermodynamics in a cosmological setting. Using a dynamical systems method, we determine the asymptotic future of plane symmetric Bianchi type I spacetimes filled with a viscous {\gamma}-fluid, keeping track of the magnitude of relative dissipative fluxes, which determines the applicability of the Israel-Stewart theory. We consider the situations when the dissipative mechanisms of shear and bulk viscosity are involved separately and simultaneously. Also, we apply two different temperature models in the full version of the theory in order to compare the results. We demonstrate that the only case when the fluid asymptotically approaches local equilibrium, and the underlying assumptions of the IS theory are therefore not violated, is that of a dissipative fluid with vanishing bulk viscosity. The truncated Israel-Stewart equations for shear viscosity are found to produce solutions which manifest pathological dynamical features and are in addition strongly sensitive with respect to the choice of initial conditions. The possible role of bulk and shear viscosity in cosmological evolution is also discussed.
We perform unbiased tests for the clumpiness of universe by confronting the Zel'dovich-Kantowski-Dyer-Roeder luminosity distance which describes the effect of local inhomogeneities on the propagation of light with the observational one estimated from measurements of standard candles, i.e., type Ia supernovae (SNe Ia) and gamma-ray bursts (GRBs). Methodologically, we first determine the light-curve fitting parameters which account for distance estimation in SNe Ia observations and luminosity/energy relations which are responsible for distance estimation of GRBs in the global fit to reconstruct the Hubble diagrams in the context of a clumpy universe. Subsequently, these Hubble diagrams allow us to achieve unbiased constraints on the matter density parameter $\Omega_m$ as well as clumpiness parameter $\eta$ which quantifies the fraction of homogeneously distributed matter within a given light cone. At 1$\sigma$ confidence level, the constraints are $\Omega_m=0.34\pm0.02$ and $\eta=1.00^{+0.00}_{-0.02}$ from the joint analysis. The results suggest that the universe full of Friedman-Lema\^{i}tre-Robertson-Walker fluid is favored by observations of standard candles with very high statistical significance. On the other hand, they may also indicate that the Zel'dovich-Kantowski-Dyer-Roeder approximation is a sufficient accurate form to describe the effects of local homogeneity in the expanding universe.
Hubble's announcement of the magnitude-redshift relation \cite{Hub29} brought about a major change in our understanding of the Universe. After tracing the pre-history of Hubble's work, and the hiatus in our understanding which his underestimate of distances led to, this review focuses on the development and success of our understanding of the expanding universe up to the present day, and the part which General Relativity plays in that success.
We propose a dark matter model in which a fermionic dark matter (DM) candidate communicates with standard model particles through two distinct portals: Higgs and vector portals. The dark sector is charged under a $U(1)'$ gauge symmetry while the standard model has a leptophobic interaction with the dark vector boson. The leading contribution of DM-nucleon elastic scattering cross section begins at one-loop level. The model meets all the constraints imposed by direct detection experiments provided by LUX and XENON100, observed relic abundance according to WMAP and Planck, and the invisible Higgs decay width measured at the LHC. It turns out that the dark matter mass in the viable parameter space can take values from a few GeV up to 1 TeV. In addition, we can find in the constrained regions of the parameter space a DM mass of $\sim 34$ GeV annihilating into $b$ quark pair, which explains the Fermi-LAT gamma-ray excess.
What is the Higgs boson telling us? What else is there, maybe supersymmetry and/or dark matter? How do we find it? These are now the big questions in collider physics that I discuss in this talk, from a personal point of view.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
A short inflationary phase may not erase all traces of the primordial universe. Associated observables include both spatial curvature and "anomalies" in the microwave background or large scale structure. The present curvature $\Omega_{K,0}$ reflects the initial curvature, $\Omega_{K,\mathrm{start}}$, and the angular size of anomalies depends on $k_\mathrm{start}$, the comoving horizon size at the onset of inflation. We estimate posteriors for $\Omega_{K,\mathrm{start}}$ and $k_\mathrm{start}$ using current data and simulations, and show that if either quantity is measured to have a non-zero value, both are likely to be observable. Mappings from $\Omega_{K,\mathrm{start}}$ and $k_\mathrm{start}$ to present-day observables depend strongly on the primordial equation of state; $\Omega_{K,0}$ spans ten orders of magnitude for a given $\Omega_{K,\mathrm{start}}$ while a simple and general relationship connects $\Omega_{K,0}$ and $k_\mathrm{start}$. We show that current bounds on $\Omega_{K,0}$ imply that if $k_\mathrm{start}$ is measurable, the curvature was already small when inflation began. Finally, since the energy density changes slowly during inflation, primordial gravitational wave constraints require that a short inflationary phase is preceded by a nontrivial pre-inflationary phase with critical implications for the expected value of $\Omega_{K,\mathrm{start}}$.
We report the discovery of the millisecond pulsar (MSP) PSR J1950+2414 ($P=4.3$ ms) in a binary system with an eccentric ($e=0.08$) orbit in Pulsar ALFA survey observations with the Arecibo telescope. Its companion star has a median mass of 0.3 $M_\odot$ and is most likely a white dwarf. Fully recycled MSPs like this one are thought to be old neutron stars spun-up by mass transfer from a companion star. This process should circularize the orbit, as is observed for the vast majority of binary MSPs, which predominantly have orbital eccentricities $e < 0.001$. However, four recently discovered binary MSPs have orbits with larger eccentricities ($0.03 < e < 0.4$); PSR J1950+2414 is only the fifth such system to be discovered. The upper limits for the the intrinsic spin period derivative and inferred surface magnetic field strength are comparable to those of the general MSP population. The large eccentricities of these systems are not compatible with the predictions of the standard recycling scenario: something unusual happened during their formation or evolution. Some of the proposed scenarios are a) the initial evolution of the pulsar in a triple system which became dynamically unstable, b) origin in an exchange encounter in an environment with high stellar density, like that of the core of a globular cluster, c) rotationally delayed accretion-induced collapse of a super-Chandrasekhar white dwarf and d) dynamical interaction of the binary with a circumbinary disk. We compare the properties of all five known eccentric MSPs with the predictions of these formation channels. We also outline how future measurements of the mass and proper motion of PSR J1950+2414 might allow us to firmly exclude some of the proposed formation scenarios.
Massive early-type galaxies commonly have gas discs which are kinematically misaligned with the stellar component. These discs feel a torque from the stars, however, and the angular momentum vectors are naively expected to align within a few dynamical times. We present results on the evolution of a misaligned gas disc in a cosmological 'zoom-in' simulation of a massive early-type galaxy from the Feedback In Realistic Environments (FIRE) project. This galaxy experiences a merger at z=0.3, which, together with a strong galactic wind, removes most of the gas disc that was in place. The galaxy subsequently reforms a gas disc through accretion of cold gas, but it is initially 120 degrees misaligned with the stellar rotation axis. This misalignment persists for about 2 Gyr before the gas-star misalignment angle drops below 20 degrees. This is about 150 times longer than the dynamical time in the central kpc and varies with galactocentric radius. The time it takes for the gaseous and stellar components to align is much longer than previously thought, because the gas disc is accreting a significant amount of mass for about 1.5 Gyr after the merger, during which the angular momentum change induced by accreted gas dominates over that induced by stellar torques. Once the gas accretion rate has decreased sufficiently, the gas disc decouples from the surrounding halo gas (which remains misaligned) and realigns with the stellar component in about 6 dynamical times, independent of radius. When stellar torques dominate the evolution of the misaligned gas disc, the centre aligns faster than the outskirts, temporarily resulting in a warped disc. We discuss the observational consequences of the long survival of our misaligned gas disc and how our results can be used to calibrate merger rate estimates from observed gas misalignments.
We present the laboratory verification of a method for re- moving the effects of frozen-flow atmospheric turbulence using a Linear Quadratic Gaussian (LQG) controller, also known as a Kalman Filter. This method, which we term "Predictive Fourier Control," can identify correlated atmospheric motions due to layers of frozen flow turbulence, and can predictively remove the effects of these correlated motions in real-time. Our laboratory verification suggests a factor of 3 improvement in the RMS residual wavefront error and a 10% improvement in measured Strehl of the system. We found that the RMS residual wavefront error was suppressed from 35.0 nm to 11.2 nm due to the use of Predictive Fourier Control, and that the far field Strehl improved from 0.479 to 0.520.
Accurate direct $N$-body simulations help to obtain detailed information about the dynamical evolution of star clusters. They also enable comparisons with analytical models and Fokker-Planck or Monte-Carlo methods. NBODY6 is a well-known direct $N$-body code for star clusters, and NBODY6++ is the extended version designed for large particle number simulations by supercomputers. We present NBODY6++GPU, an optimized version of NBODY6++ with hybrid parallelization methods (MPI, GPU, OpenMP, and AVX/SSE) to accelerate large direct $N$-body simulations, and in particular to solve the million-body problem. We discuss the new features of the NBODY6++GPU code, benchmarks, as well as the first results from a simulation of a realistic globular cluster initially containing a million particles. For million-body simulations, NBODY6++GPU is $400-2000$ times faster than NBODY6 with 320 CPU cores and 32 NVIDIA K20X GPUs. With this computing cluster specification, the simulations of million-body globular clusters including $5\%$ primordial binaries require about an hour per half-mass crossing time.
The recent discovery by Cantalupo et al. (2014) of the largest (~500 kpc) and luminous Ly-alpha nebula associated with the quasar UM287 (z=2.279) poses a great challenge to our current understanding of the astrophysics of the halos hosting massive z~2 galaxies. Either an enormous reservoir of cool gas is required $M\simeq10^{12}$ $M_{\odot}$, exceeding the expected baryonic mass available, or one must invoke extreme gas clumping factors not present in high-resolution cosmological simulations. However, observations of Ly-alpha emission alone cannot distinguish between these two scenarios. We have obtained the deepest ever spectroscopic integrations in the HeII and CIV lines with the goal of detecting extended line emission, but detect neither line to a 3$\sigma$ limiting SB $\simeq10^{-18}$ erg/s/cm$^2$/arcsec$^2$. We construct models of the expected emission spectrum in the highly probable scenario that the nebula is powered by photoionization from the central hyper-luminous quasar. The non-detection of HeII implies that the nebular emission arises from a mass $M_{\rm c}\lesssim6.4\times10^{10}$ $M_{\odot}$ of cool gas on ~200 kpc scales, distributed in a population of remarkably dense ($n_{\rm H}\gtrsim3$ cm$^{-3}$) and compact ($R\lesssim20$ pc) clouds, which would clearly be unresolved by current cosmological simulations. Given the large gas motions suggested by the Ly-alpha line ($v\simeq$ 500 km/s), it is unclear how these clouds survive without being disrupted by hydrodynamic instabilities. Our study serves as a benchmark for future deep integrations with current and planned wide-field IFU such as MUSE, KCWI, and KMOS. Our work suggest that a $\simeq$ 10 hr exposure would likely detect ~10 rest-frame UV/optical emission lines, opening up the possibility of conducting detailed photoionization modeling to infer the physical state of gas in the CGM.
Star formation is inefficient. Only a few percent of the available gas in molecular clouds forms stars, leading to the observed low star formation rate (SFR). The same holds when averaged over many molecular clouds, such that the SFR of whole galaxies is again surprisingly low. Indeed, considering the low temperatures, molecular clouds should be highly gravitationally unstable and collapse on their global mean freefall timescale. And yet, they are observed to live about 10-100 times longer, i.e., the SFR per freefall time (SFR_ff) is only a few percent. Thus, other physical mechanisms must provide support against quick global collapse. Magnetic fields, turbulence and stellar feedback have been proposed as stabilising agents, but it is still unclear which of these processes is the most important and what their relative contributions are. Here we run high-resolution simulations including gravity, turbulence, magnetic fields, and jet/outflow feedback. We confirm that clouds collapse on a mean freefall time, if only gravity is considered, producing stars at an unrealistic rate. In contrast, if turbulence, magnetic fields, and feedback are included step-by-step, the SFR is reduced by a factor of 2-3 with each additional physical ingredient. When they all act in concert, we find a constant SFR_ff = 0.04, currently the closest match to observations, but still about a factor of 2-4 higher than the average. A detailed comparison with other simulations and with observations leads us to conclude that only models with turbulence producing large virial parameters, and including magnetic fields and feedback can produce realistic SFRs.
NGC3718 is a LINER $L1.9$ galaxy, lying at a distance of about $\sim 17.4$ Mpc away from earth and its similarities with NGC5128 often award it the name "northern Centaurus A". We use high angular resolution ($\sim100$ mas) e-Merlin radio and SUBARU NIR ($\sim170$ mas) data, to take a detailed view of the processes taking place in its central region. In order to preserve some objectivity in our interpretation, we combine our results with literature values and findings from previous studies. Our NIR maps suggest, on one hand, that towards the stellar bulge there are no large scale absorption phenomena caused by the apparent dust lane and, on the other, that there is a significant (local) contribution from hot ($\sim1000$ K) dust to the nuclear NIR emission. The position where this takes place appears to be closer to the offset compact radio emission from our e-Merlin $6$ cm map, lying offset by $\sim4.25$ pc from the center of the underlying stellar bulge. The shape of the radio map suggests the presence of one (or possibly two, forming an X-shape) bipolar structure(s) $\sim1$ ($\sim0.6$) arcsec across, which combined with the balance between the gas and the stellar velocity dispersions and the presence of hard X-ray emission, point towards effects expected by AGN feedback. We also argue that NGC3718 has a "core" in its surface brightness profile, despite the fact that it is a gas-rich galaxy and we discuss its mixed photometric and spectroscopic characteristics. The latter combined with the observed spatial and radio offsets, the relative redshift between the broad and the narrow $H{\mathrm{\alpha}}$ line, the limited star formation activity and AGN feedback, strongly imply the existence of an SMBH recoil. Finally, we discuss a possible interpretation, that could naturally incorporate all these findings into one physically consistent picture.
We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release (Desvignes et al. in prep.). Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar in this dataset, along with common correlated signals including clock, and Solar System ephemeris errors to obtain a robust 95$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $\Omega_\mathrm{gw}(f)h^2 < 3.6\times10^{-10}$ at 2.8 nHz. We show that performing such an analysis when fixing the intrinsic noise parameters for the individual pulsars leads to an erroneously stringent upper limit, by a factor $\sim 1.7$. We obtain a difference in the logarithm of the Bayesian evidence between models that include either a correlated background, or uncorrelated common red noise of $-1.0 \pm 0.5$, indicating no support for the presence of a correlated GWB in this dataset. We discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $G\mu/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit $G\mu/c^2<1.3\times10^{-7}$, identical to that set by the Planck Collaboration, combining Planck and high-$\ell$ Cosmic Microwave Background data from other experiments. (Abridged)
We present the results of a Lya profile analysis of 12 Lya emitters (LAEs) at z = 2.2 with high-resolution Lya spectra. We find that all 12 objects have a Lya profile with the main peak redward of the systemic redshift defined by nebular lines, and five have a weak, secondary peak blueward of the systemic redshift (blue bump). The average velocity offset of the red main peak (the blue bump, if any) with respect to the systemic redshift is Delta_v_Lya,r = 174+/- 19 km s-1 (Delta_v_Lya,b = -316+/-45 km s-1), which is smaller than (comparable to) that of Lyman-break galaxies (LBGs). The outflow velocities inferred from metal absorption lines in three individual and one stacked spectra are comparable to those of LBGs. The expanding shell model constructed by Verhamme et al. (2006) reproduces not only the Lya profiles but also other observed quantities including the outflow velocity and the FWHM of nebular lines for the non-blue bump objects. On the other hand, the model predicts too high FWHMs of nebular lines for the blue bump objects, although this discrepancy may disappear if we introduce additional Lya photons produced by gravitational cooling. We show that the small Delta_v_Lya,r values of our sample can be explained by low neutral-hydrogen column densities of log(NHI) = 18.9 cm-2 on average. This value is more than one order of magnitude lower than those of LBGs but is consistent with recent findings that LAEs have high ionization parameters and low Hi gas masses. This result suggests that low NHI values, giving reduced numbers of resonant scattering of Lya photons, are the key to the strong Lya emission of LAEs.
Red Supergiants (RSGs) are cool (~4000K), highly luminous stars (L - 10^5 Lsun), and are among the brightest near-infrared (NIR) sources in star-forming galaxies. This makes them powerful probes of the properties of their host galaxies, such as kinematics and chemical abundances. We have developed a technique whereby metallicities of RSGs may be extracted from a narrow spectral window around 1{\mu}m from only moderate resolution data. The method is therefore extremely efficient, allowing stars at large distances to be studied, and so has tremendous potential for extragalactic abundance work. Here, we present an abundance study of the Large and Small Magellanic Clouds (LMC and SMC respectively) using samples of 9-10 RSGs in each. We find average abundances for the two galaxies of [Z]LMC = -0.37 +/- 0.14 and [Z]SMC = -0.53 +/- 0.16 (with respect to a Solar metallicity of Zsun=0.012). These values are consistent with other studies of young stars in these galaxies, and though our result for the SMC may appear high it is consistent with recent studies of hot stars which find 0.5-0.8dex below Solar. Our best-fit temperatures are on the whole consistent with those from fits to the optical-infrared spectral energy distributions, which is remarkable considering the narrow spectral range being studied. Combined with our recent study of RSGs in the Galactic cluster Per OB1, these results indicate that this technique performs well over a range of metallicities, paving the way for forthcoming studies of more distant galaxies beyond the Local Group.
(abridged) We present a dynamical analysis of the extended stellar stream encircling NGC 1097. Within a statistical framework, we model its surface brightness using mock streams as in Amorisco (2015) and deep imaging data from the CHART32 telescope (Stellar Tidal Stream Survey). We reconstruct the post-infall evolution of the progenitor, which has experienced 3 pericentric passages and lost more than 2 orders of magnitude in mass. At infall, $5.4\pm0.6$ Gyr ago, the progenitor was a disky dwarf with mass of $\log_{10}[m(<3.4\pm1 {\rm kpc})/ M_\odot]=10.35\pm0.25$. We illustrate how the 90$^\circ$ turn in the stream, identifying the `dog leg', is the signature of the progenitor's prograde rotation. Today, the remnant is a nucleated dwarf, with a LOS velocity of $v_{\rm p, los}^{\rm obs}=-30\pm 30$ kms$^{-1}$, and a luminosity of $3.3\times 10^7 L_{V,\odot}$ (Galianni et al. 2010). Our independent analysis predicts $v_{\rm p, los}=-51^{-17}_{+14}$ kms$^{-1}$, and measures $\log_{10}(m/ M_\odot)=7.4^{+0.6}_{-0.8}$, so that the compact nucleus is soon becoming a low-luminosity UCD. We find that NGC 1097 has a mass of $M_{200}=1.8^{+0.5}_{-0.4} \times 10^{12}\; M_{\odot}$, and its concentration $c_{200}=6.7^{+2.4}_{-1.3}$ is in agreement with LCDM. The stream is described almost down to the noise in a spherical host potential, we find this would not be possible if the halo was substantially triaxial at large radii. Its morphology shows that the slope of the total density profile bends from an inner $\gamma(r_{\rm peri})=1.5\pm0.15$. The progenitor's orbit reaches $r_{\rm apo}=150\pm 15$ kpc, more than a half of the virial radius of the host, so that, for the first time on an individual extragalactic halo, we measure the outer density slope, $\gamma(0.6r_{200,c})=3.9\pm0.5$. This demonstrates the promise of the newborn field of detailed, statistical modelling of extragalactic tidal streams.
The mean absolute extinction towards the central parsec of the Milky Way is A_K~3 mag, including both foreground and Galactic center dust. Here we present a measurement of dust extinction within the Galactic old nuclear star cluster (NSC), based on combining differential extinctions of NSC stars with their u_l proper motions along Galactic longitude. Extinction within the NSC preferentially affects stars at its far side, and because the NSC rotates, this causes higher extinctions for NSC stars with negative u_l, as well as an asymmetry in the u_l-histograms. We model these effects using an axisymmetric dynamical model of the NSC in combination with simple models for the dust distribution. Comparing the predicted asymmetry to data for ~7100 stars in several NSC fields, we find that dust associated with the Galactic center mini-spiral with extinction A_K~=0.15-0.8 mag explains most of the data. The largest extinction A_K~=0.8 mag is found in the region of the Western arm of the mini-spiral. Comparing with total A_K determined from stellar colors, we determine the extinction in front of the NSC. Finally, we estimate that for a typical extinction of A_K~=0.4 the statistical parallax of the NSC changes by ~0.4%.
We study 14 large solar jets observed in polar coronal holes. In EUV movies from SDO/AIA, each jet appears similar to most X-ray jets and EUV jets that erupt in coronal holes, but each is exceptional in that it goes higher than most, so high that it is observed in the outer corona beyond 2.2 RSun in images from the SOHO/LASCO/C2 coronagraph. From AIA He II 304 {\AA} movies and LASCO/C2 running-difference images of these high-reaching jets, we find: (1) the front of the jet transits the corona below 2.2 RSun at a speed typically several times the sound speed; (2) each jet displays an exceptionally large amount of spin as it erupts; (3) in the outer corona, most of the jets display measureable swaying and bending of a few degrees in amplitude; in three jets the swaying is discernibly oscillatory with a period of order 1 hour. These characteristics suggest that the driver in these jets is a magnetic-untwisting wave that is basically a large-amplitude (i.e., non-linear) torsional Alfven wave that is put into the reconnected open field in the jet by interchange reconnection as the jet erupts. From the measured spinning and swaying we estimate that the magnetic-untwisting wave loses most of its energy in the inner corona below 2.2 RSun. We point out that the torsional waves observed in Type-II spicules might dissipate in the corona in the same way as the magnetic-untwisting waves in our big jets and thereby power much of the coronal heating in coronal holes.
We present multi-wavelength observations and modeling of Gamma-ray Bursts (GRBs) that exhibit a simultaneous re-brightening in their X-ray and optical light curves, and are also detected at radio wavelengths. We show that the re-brightening episodes can be modeled by injection of energy into the blastwave and that in all cases the energy injection rate falls within the theoretical bounds expected for a distribution of energy with ejecta Lorentz factor. Our measured values of the circumburst density, jet opening angle, and beaming corrected kinetic energy are consistent with the distribution of these parameters for long-duration GRBs at both z~1 and z>6, suggesting that the jet launching mechanism and environment of these events are similar to that of GRBs that do not have bumps in their light curves. However, events exhibiting re-brightening episodes have lower radiative efficiencies than average, suggesting that a majority of the kinetic energy of the outflow is carried by slow-moving ejecta, which is further supported by steep measured distributions of the ejecta energy as a function of Lorentz factor. We do not find evidence for reverse shocks over the energy injection period, implying that the onset of energy injection is a gentle process. We further show that GRBs exhibiting simultaneous X-ray and optical re-brightenings are likely the tail of a distribution of events with varying rates of energy injection, forming the most extreme events in their class. Future X-ray observations of GRB afterglows with Swift and its successors will thus likely discover several more such events, while radio follow-up and multi-wavelength modeling of similar events will unveil the role of energy injection in GRB afterglows.
We present Gemini GMOS-IFU data of eight compact low-mass early-type galaxies (ETGs) in the Virgo cluster. We analyse their stellar kinematics, stellar population, and present two-dimensional maps of these properties covering the central 5"x 7" region. We find a large variety of kinematics: from non- to highly-rotating objects, often associated with underlying disky isophotes revealed by deep images from the Next Generation Virgo Cluster Survey. In half of our objects, we find a centrally-concentrated younger and more metal-rich stellar population. We analyze the specific stellar angular momentum through the lambdaR parameter and find six fast-rotators and two slow-rotators, one having a thin counter-rotating disk. We compare the local galaxy density and stellar populations of our objects with those of 39 more extended low-mass Virgo ETGs from the SMAKCED survey and 260 massive ($M>10^{10}$\Msun) ETGs from the A3D sample. The compact low-mass ETGs in our sample are located in high density regions, often close to a massive galaxy and have, on average, older and more metal-rich stellar populations than less compact low-mass galaxies. We find that the stellar population parameters follow lines of constant velocity dispersion in the mass-size plane, smoothly extending the comparable trends found for massive ETGs. Our study supports a scenario where low-mass compact ETGs have experienced long-lived interactions with their environment, including ram-pressure stripping and gravitational tidal forces, that may be responsible for their compact nature.
We present the discovery of two new X-ray transients in archival Chandra data. The first transient, XRT 110103, occurred in January 2011 and shows a sharp rise of at least three orders of magnitude in count rate in less than 10 s, a flat peak for about 20 s and decays by two orders of magnitude in the next 60 s. We find no optical or infrared counterpart to this event in preexisting survey data or in an observation taken by the SIRIUS instrument at the Infrared Survey Facility 2.1 yr after the transient, providing limiting magnitudes of J>18.1, H>17.6 and Ks>16.3. This event shows similarities to the transient previously reported in Jonker et al. which was interpreted as the possible tidal disruption of a white dwarf by an intermediate mass black hole. We discuss the possibility that these transients originate from the same type of event. If we assume these events are related a rough estimate of the rates gives 1.4*10^5 per year over the whole sky with a peak 0.3-7 keV X-ray flux greater than 2*10^-10 erg cm^-2 s^-1 . The second transient, XRT 120830, occurred in August 2012 and shows a rise of at least three orders of magnitude in count rate and a subsequent decay of around one order of magnitude all within 10 s, followed by a slower quasi-exponential decay over the remaining 30 ks of the observation. We detect a likely infrared counterpart with magnitudes J=16.70+/-0.06, H=15.92+/-0.04 and Ks=15.37+/-0.06 which shows an average proper motion of 74+/-19 milliarcsec per year compared to archival 2MASS observations. The JHKs magnitudes, proper motion and X-ray flux of XRT 120830 are consistent with a bright flare from a nearby late M or early L dwarf.
We compare the performance of mass estimators for elliptical galaxies that rely on the directly observable surface brightness and velocity dispersion profiles, without invoking computationally expensive detailed modeling. These methods recover the mass at a specific radius where the mass estimate is expected to be least sensitive to the anisotropy of stellar orbits. One method (Wolf et al. 2010) uses the total luminosity-weighted velocity dispersion and evaluates the mass at a 3D half-light radius $r_{1/2}$, i.e., it depends on the GLOBAL galaxy properties. Another approach (Churazov et al. 2010) estimates the mass from the velocity dispersion at a radius $R_2$ where the surface brightness declines as $R^{-2}$, i.e., it depends on the LOCAL properties. We evaluate the accuracy of the two methods for analytical models, simulated galaxies and real elliptical galaxies that have already been modeled by the Schwarzschild's orbit-superposition technique. Both estimators recover an almost unbiased circular speed estimate with a modest RMS scatter ($\lesssim 10 \%$). Tests on analytical models and simulated galaxies indicate that the local estimator has a smaller RMS scatter than the global one. We show by examination of simulated galaxies that the projected velocity dispersion at $R_2$ could serve as a good proxy for the virial galaxy mass. For simulated galaxies the total halo mass scales with $\sigma_p(R_2)$ as $M_{vir} \left[M_{\odot}h^{-1}\right] \approx 6\cdot 10^{12} \left( \frac{\sigma_p(R_2)}{200\, \rm km\, s^{-1}} \right)^{4}$ with RMS scatter $\approx 40 \%$.
A0535+26 is a slowly rotating pulsar accreting from the wind of a massive Be star, and that exhibits two cyclotron absorption lines in its X-ray spectrum, at about 45 and 100 keV, respectively. Unlike similar sources, no significant variations of the energy of its cyclotron lines with flux were observed to date. The bright outburst of February 2011 thus offers a unique occasion to probe this peculiar behavior at flux levels not yet observed with present-day instruments. Here we report on the spectral and timing analysis of the data from the spectrometer SPI on-board INTEGRAL collected during the outburst. At the peak of the outburst the estimated luminosity is ~4.9x10^37 erg/s. The fundamental cyclotron feature is detected at all flux levels, and its centroid energy is positively-correlated with the flux of the source, confirming that A0535+26 is accreting at a sub-critical regime. The correlation seems to fall off at ~10^37 erg/s, suggesting the transition from a Coulomb-stopping regime to a gas-mediated shock regime. From the timing analysis we found that the pulsar was spinning up during most of the outburst, and that the spin-up rate correlates with the flux of the source, albeit the correlation is steeper than the one expected from the standard disk accretion theory. Finally, we show that the pulse profile of the source changes dramatically as the flux increases. At high luminosity the profile is highly asymmetric, implying an asymmetry in the geometry of the accretion flow.
We used the IRAM 30m telescope to observe the frequency range [86-116]GHz towards the central regions of the starburst galaxies M83, M82, and NGC253, the AGNs M51, NGC1068, and NGC7469, and the ULIRGs Arp220 and Mrk231. Assuming LTE conditions, we calculated the column densities of 27 molecules and 10 isotopologues. Among others, we report the first tentative detections of CH3CHO, HNCO, and NS in M82 and, for the first time in the extragalactic medium, HC5N in NGC253. Halpha recombination lines were only found in M82 and NGC253. Vibrationally excited lines of HC3N were only detected in Arp220. CH3CCH emission is only seen in the starburst-dominated galaxies. By comparison of the fractional abundances among the galaxies, we looked for the molecules that are best suited to characterise the chemistry of starbursts, AGNs and ULIRGs, as well as the differences among galaxies within the same group.
We have mapped the NGC 2023 reflection nebula in [CII] and CO(11--10) with the heterodyne receiver GREAT on SOFIA and obtained slightly smaller maps in 13CO(3--2), CO(3--2), CO(4--3), CO(6--5), and CO(7--6) with APEX in Chile. We use these data to probe the morphology, kinematics, and physical conditions of the C II region, which is ionized by FUV radiation from the B2 star HD37903. The [CII] emission traces an ellipsoidal shell-like region at a position angle of ~ -50 deg, and is surrounded by a hot molecular shell. In the southeast, where the C II region expands into a dense, clumpy molecular cloud ridge, we see narrow and strong line emission from high-J CO lines, which comes from a thin, hot molecular shell surrounding the [CII] emission. The [CII] lines are broader and show photo evaporating gas flowing into the C II region. Based on the strength of the [13CII] F=2--1 line, the [CII] line appears to be somewhat optically thick over most of the nebula with an optical depth of a few. We model the physical conditions of the surrounding molecular cloud and the PDR emission using both RADEX and simple PDR models. The temperature of the CO emitting PDR shell is ~ 90 -- 120 K, with densities of 10^5 -- 10^6 cm^-3, as deduced from RADEX modeling. Our PDR modeling indicates that the PDR layer where [CII] emission dominates has somewhat lower densities, 10^4 to a few times 10^5 cm^-3
New CCD photometry is presented for the hot overcontact binary DK Cyg, together with reasonable explanations for the light and period variations. Historical light and velocity curves from 1962 to 2012 were simultaneously analyzed with the Wilson-Devinney (W-D) synthesis code. The brightness disturbances were satisfactorily modeled by applying a magnetic cool spot on the primary star. Based on 261 times of minimum light including 116 new timings and spanning more than 87 yrs, a period study reveals that the orbital period has varied due to a periodic oscillation superposed on an upward parabola. The period and semi-amplitude of the modulation are about 78.1 yrs and 0.0037 d, respectively. This detail is interpreted as a light-travel-time effect due to a circumbinary companion with a minimum mass of $M_3 =0.065 $M_\odot$, within the theoretical limit of $\sim$0.07 M$_\odot$ for a brown dwarf star. The observed period increase at a fractional rate of $+$2.74 $\times $10$^{-10}$ is in excellent agreement with that calculated from our W-D synthesis. Mass transfer from the secondary to the primary component is mainly responsible for the secular period change. We examined the evolutionary status of the DK Cyg system from the absolute dimensions.
We present a measurement of the atmospheric $\nu_e$ spectrum at energies between 0.1 TeV and 100 TeV using data from the first year of the complete IceCube detector. Atmospheric $\nu_e$ originate mainly from the decays of kaons produced in cosmic-ray air showers. This analysis selects 1078 fully contained events in 332 days of livetime, then identifies those consistent with particle showers. A likelihood analysis with improved event selection extends our previous measurement of the conventional $\nu_e$ fluxes to higher energies. The data constrain the conventional $\nu_e$ flux to be $1.3^{+0.4}_{-0.3}$ times a baseline prediction from a Honda's calculation, including the knee of the cosmic-ray spectrum. A fit to the kaon contribution ($\xi$) to the neutrino flux finds a kaon component that is $\xi =1.3^{+0.5}_{-0.4}$ times the baseline value. The fitted/measured prompt neutrino flux from charmed hadron decays strongly depends on the assumed astrophysical flux and shape. If the astrophysical component follows a power law, the result for the prompt flux is $0.0^{+3.0}_{-0.0}$ times a calculated flux based on the work by Enberg, Reno and Sarcevic.
We report twelve variables not previously detected in the globular cluster M5 (NGC 5904); one SX Phe and eleven semi-regular variables (SR). Their identifications, equatorial coordinates, ephemerides, and light curves are given. Furthermore, we have explored the light curves of a group of stars whose variability has not been confirmed and that are marked as probable non- variables in the CVSGC. Finally, we offer detailed identifications for some of the known variables in crowded regions that were misidentified in previous studies. We shall also address the cases of the cataclysmic variable or U Gem type V101 and of the variable blue straggler V159.
Infrared-faint radio sources (IFRS) form a new class of galaxies characterised by radio flux densities between tenths and tens of mJy and faint or absent infrared counterparts. It has been suggested that these objects are radio-loud active galactic nuclei (AGNs) at significant redshifts (z >~ 2). Whereas the high redshifts of IFRS have been recently confirmed based on spectroscopic data, the evidence for the presence of AGNs in IFRS is mainly indirect. So far, only two AGNs have been unquestionably confirmed in IFRS based on very long baseline interferometry (VLBI) observations. In this work, we test the hypothesis that IFRS contain AGNs in a large sample of sources using VLBI. We observed 57 IFRS with the Very Long Baseline Array (VLBA) down to a detection sensitivity in the sub-mJy regime and detected compact cores in 35 sources. Our VLBA detections increase the number of VLBI-detected IFRS from 2 to 37 and provide strong evidence that most - if not all - IFRS contain AGNs. We find that IFRS have a marginally higher VLBI detection fraction than randomly selected sources with mJy flux densities at arcsec-scales. Moreover, our data provide a positive correlation between compactness - defined as the ratio of milliarcsec- to arcsec-scale flux density - and redshift for IFRS, but suggest a decreasing mean compactness with increasing arcsec-scale radio flux density. Based on these findings, we suggest that IFRS tend to contain young AGNs whose jets have not formed yet or have not expanded, equivalent to very compact objects. We found two IFRS that are resolved into two components. The two components are spatially separated by a few hundred milliarcseconds in both cases. They might be components of one AGN, a binary black hole, or the result of gravitational lensing.
If vector type perturbations are present in the primordial plasma before recombination, the generation of magnetic fields is known to be inevitable through the Harrison mechanism. In the context of the standard cosmological perturbation theory, non-linear couplings of first-order scalar perturbations create second-order vector perturbations, which generate magnetic fields. Here we reinvestigate the generation of magnetic fields at second-order in cosmological perturbations on the basis of our previous study, and extend it by newly taking into account the time evolution of purely second-order vector perturbations with a newly developed second-order Boltzmann code. We confirm that the amplitude of magnetic fields from the product-terms of the first-order scalar modes is consistent with the result in our previous study. However, we find, both numerically and analytically, that the magnetic fields from the purely second-order vector perturbations partially cancel out the magnetic fields from one of the product-terms of the first-order scalar modes, in the tight coupling regime in the radiation dominated era. Therefore, the amplitude of the magnetic fields on small scales, $k \gtrsim 10~h{\rm Mpc}^{-1}$, is smaller than the previous estimates. The amplitude of the generated magnetic fields at cosmological recombination is about $B_{\rm rec} =5.0\times 10^{-24}~{\rm Gauss}$ on $k = 5.0 \times 10^{-1}~h{\rm Mpc}^{-1}$. Finally, we discuss the reason of the discrepancies that exist in estimates of the amplitude of magnetic fields among other authors.
HESS J1943+213 is a TeV source close to the Galactic plane proposed to be a BL Lac object. Our high resolution EVN observation failed to recover two thirds of the source flux density detected simultaneously by the WSRT. Our recent e-MERLIN observations in L and C bands show only a point source with flux density comparable to the EVN detection. Thus the structure responsible for the missing flux density has to be larger than 2". It may be related to the presumed extragalactic source (thus would have a kpc-scale size), or to the Galactic foreground material close to the line of sight to the source.
We devise and explore an iterative optimization procedure for controlling
particle populations in particle-in-cell (PIC) codes via merging and splitting
of computational macro-particles.
Our approach, is to compute an optimal representation of the global particle
phase space structure while decreasing or increasing the entire particle
population, based on k-means clustering of the data. In essence the procedure
amounts to merging or splitting particles by statistical means, throughout the
entire simulation volume in question, while minimizing a 6-dimensional total
distance measure to preserve the physics.
Particle merging is by far the most demanding procedure when considering
conservation laws of physics; it amounts to lossy compression of particle phase
space data. We demonstrate that our k-means approach conserves energy and
momentum to high accuracy, even for high compression ratios of about four ---
i.e., N_f/N_i = 0.25.
Interestingly, we find that the most intuitive naiive approach to particle
splitting does not yield accurate results over hundreds of simulation time
steps. Rather, a more complex particle splitting procedure is superior.
Implementation and testing is done using an electromagnetic PIC code, the
Photon-Plasma code. Nonetheless, the k-means framework is general; it is not
limited to Vlasov-Maxwell type PIC codes. We discuss advantages and drawbacks
of this optimal phase space reconstruction.
Together with their fingerprint modes, molecules carry coherent vibrations of all their atoms (phonons). Phonon spectra extend from $\sim$20 to more than $10^{4}\,\mu$m, depending on molecular size. These spectra are discrete but large assemblies of molecules of the same family, differing only by minor structural details, will produce continua. As such assemblies are expected to exist in regions where dust accumulates, they are bound to contribute to the observed continua underlying the Unidentified Infrared Bands and the 21-mum band of planetary nebulae as well as to the diffuse galactic emission surveyed by the Planck astronomical satellite and other means. The purpose of this work is to determine, for carbon-rich molecules, the intensity of such continua and their extent into the millimetric range, and to evaluate their detectability in this range. The rules governing the spectral distributions of phonons are derived and shown to differ from those which obtain in the solid state. Their application allow the extinction cross-section per H atom, and its maximum wavelength, to be determined as a function of molecular size and dimensionality. Chemical modeling of more than 15 large molecules illustrate these results. It is found that the maximum phonon wavelength of a 2D structure increases roughly as the square of its larger dimension. Spectral energy distributions were computed as far as 4000 mum, for molecules up to 50 A{\deg} in length.
A revolution in radio receiving technology is underway with the development of densely packed phased arrays for radio astronomy. This technology can provide an exceptionally large field of view, while at the same time sampling the sky with high angular resolution. Such an instrument, with a field of view of over 100 square degrees, is ideal for performing fast, all-sky, surveys, such as the "intensity mapping" experiment to measure the signature of Baryonic Acoustic Oscillations in the HI mass distribution at cosmological redshifts. The SKA, built with this technology, will be able to do a billion galaxy survey. I will present a very brief introduction to radio interferometry, as well as an overview of the Square Kilometre Array project. This will be followed by a description of the EMBRACE prototype and a discussion of results and future plans.
We discuss whether modern machine learning methods can be used to characterize the physical nature of the large number of objects sampled by the modern multi-band digital surveys. In particular, we applied the MLPQNA (Multi Layer Perceptron with Quasi Newton Algorithm) method to the optical data of the Sloan Digital Sky Survey - Data Release 10, investigating whether photometric data alone suffice to disentangle different classes of objects as they are defined in the SDSS spectroscopic classification. We discuss three groups of classification problems: (i) the simultaneous classification of galaxies, quasars and stars; (ii) the separation of stars from quasars; (iii) the separation of galaxies with normal spectral energy distribution from those with peculiar spectra, such as starburst or starforming galaxies and AGN. While confirming the difficulty of disentangling AGN from normal galaxies on a photometric basis only, MLPQNA proved to be quite effective in the three-class separation. In disentangling quasars from stars and galaxies, our method achieved an overall efficiency of 91.31% and a QSO class purity of ~95%. The resulting catalogue of candidate quasars/AGNs consists of ~3.6 million objects, of which about half a million are also flagged as robust candidates, and will be made available on CDS VizieR facility.
We present spectroscopy of nine planetary nebulae (PNe) in the outskirts of M31, all but one obtained with the 10.4m GTC telescope. These sources extend our previous study of the oxygen abundance gradient of M31 to galactocentric radii as large as 100 kpc. None of the targets are bona fide members of a classical, metal-poor and ancient halo. Two of the outermost PNe have solar oxygen abundances, as well as radial velocities consistent with the kinematics of the extended disk of M31. The other PNe have a slightly lower oxygen content ([O/H] ~ -0.4) and in some cases large deviations from the disk kinematics. These PNe support the current view that the external regions of M31 are the result of a complex interaction and merger process, with evidence for a widespread population of solar-metallicity stars produced in a starburst that occurred ~2 Gyr ago.
The age calibration of the Washington deltaT1 index is mainly used to estimate ages of star clusters older than 1 Gyr, no age-metallicity degeneracy effect is considered. We have profusely exploited synthetic T1 versus C-T1 colour magnitude diagrams aiming at exploring the intrinsic behaviour of the deltaT1 index. The analysis shows that deltaT1 varies with age and metal content as well. In general, the dependence on age weakens for ages greater than ~ 6 Gyr, and results even less sensitive to age as the metallicity decreases. For ages younger than ~ 5 Gyr deltaT1 shows a strong correlation with both age and metallicity. The deltaC index -defined as deltaT1 for the C passband- is also a combined measurement of age and metallicity. We introduce a new age-metallicity diagnostic diagram, deltaT1 versus deltaC - deltaT1, which has shown the ability of unambiguously providing age and metallicity estimates, simultaneously. The new procedure allows to derive ages from 1 up to 13 Gyr and metallicities [Fe/H] from -2.0 up to +0.5 dex, and is independent of the cluster reddening and distance modulus. It does solve the constraints found in the deltaT1 index and surpasses the performance of the standard giant branch metallicity method. All these features make the diagnostic diagram a powerful tool for estimating accurate ages as well as metallicities.
Compressed sensing theory is slowly making its way to solve more and more astronomical inverse problems. We address here the application of sparse representations, convex optimization and proximal theory to radio interferometric imaging. First, we expose the theory behind interferometric imaging, sparse representations and convex optimization, and second, we illustrate their application with numerical tests with SASIR, an implementation of the FISTA, a Forward-Backward splitting algorithm hosted in a LOFAR imager. Various tests have been conducted in Garsden et al., 2015. The main results are: i) an improved angular resolution (super resolution of a factor ~2) with point sources as compared to CLEAN on the same data, ii) correct photometry measurements on a field of point sources at high dynamic range and iii) the imaging of extended sources with improved fidelity. SASIR provides better reconstructions (five time less residuals) of the extended emissions as compared to CLEAN. With the advent of large radiotelescopes, there is scope of improving classical imaging methods with convex optimization methods combined with sparse representations.
The extreme physical conditions of Gamma Ray Bursts can constitute a useful observational laboratory to test theories of gravity where very high curvature regimes are involved. Here we propose a sort of curvature engine capable, in principle, of explaining the huge energy emission of Gamma Ray Bursts. Specifically, we investigate the emission of radiation by charged particles non-minimally coupled to the gravitational background where higher order curvature invariants are present. The coupling gives rise to an additional force inducing a non-geodesics motion of particles. This fact allows a strong emission of radiation by gravitationally accelerated particles. As we will show with some specific model, the energy emission is of the same order of magnitude of that characterizing the Gamma Ray Burst physics. Alternatively, strong curvature regimes can be considered as a natural mechanism for the generation of highly energetic astrophysical events.
Gamma-ray bursts (GRBs) are violent explosions, coming from cosmological distances. They are detected in gamma-rays (also X-rays, UV, optical, radio) almost every day, and have typical durations of a few seconds to a few minutes. Some GRBs have been reported with extraordinary duration of 10^4 sec. These are called Ultra Long GRBs. It has been debated whether these form a new distinct class of events or whether they are similar to long GRBs. According to Blandford & Znajek (1977), the spin energy of a rotating black hole can be extracted electromagnetically, should the hole be endowed with a magnetic field supported by electric currents in a surrounding disk. We argue that this can be the case for the central engines of GRBs and we show that the duration of the burst depends on the magnetic flux accumulated on the event horizon of the black hole. We thus estimate the surface magnetic field of a possible progenitor star, and we conclude that an Ultra Long GRB may originate from a progenitor star with a relatively low magnetic field.
Numerical algorithms to load relativistic Maxwell distributions in particle-in-cell (PIC) and Monte-Carlo simulations are presented. For stationary relativistic Maxwellian, the inverse transform method and the Sobol algorithm are reviewed. To boost particles to obtain relativistic shifted-Maxwellian, two rejection methods are proposed in a physically transparent manner. Their acceptance efficiencies are ${\approx}50\%$ for generic cases and $100\%$ for symmetric distributions. They can be combined with arbitrary base algorithms.
Cosmic ray propagation is diffusive because of pitch angle scattering by waves. We demonstrate that if the high-amplitude magnetohydrodynamic turbulence with $\tilde B/\langle B\rangle \sim 1$ is present on top of the mean field gradient, the diffusion becomes asymmetric. As an example, we consider the vertical transport of cosmic rays in our Galaxy propagating away from a point-like source. We solve this diffusion problem analytically using a one-dimensional Markov chain analysis. We obtained that the cosmic ray density markedly differs from the standard diffusion prediction and has a sizable effect on their distribution throughout the galaxy. The equation for the continuous limit is also derived, which shows limitations of the convection-diffusion equation.
We study the possibility of probing dark energy behaviour using gravitational wave experiments like LISA and Advanced LIGO. Using two popular parameterizations for dark energy equation of state, we show that with current sensitivities of LISA and Advanced LIGO to detect the stochastic gravitational waves, it is possible to probe a large section of parameter space for the dark energy equation of state which is allowed by present cosmological observations.
Transitional disks with large dust cavities are important laboratories to study planet formation and disk evolution. Cold gas may still be present inside these cavities, but the quantification of this gas is challenging. The gas content is important to constrain the origin of the dust cavity. We use Atacama Large Millimeter/submillimeter Array (ALMA) observations of 12CO 6--5 and 690 GHz (Band 9) continuum of five well-studied transitional disks. In addition, we analyze previously published Band 7 observations of a disk in 12CO 3--2 line and 345 GHz continuum. The observations are used to set constraints on the gas and dust surface density profiles, in particular the drop delta-gas of the gas density inside the dust cavity. The physical-chemical modeling code DALI is used to analyze the gas and dust images simultaneously. We model SR21, HD135344B, LkCa15, SR24S and RXJ1615-3255 (Band 9) and J1604-2130 (Band 7). The SED and continuum visibility curve constrain the dust surface density. Subsequently, the same model is used to calculate the 12CO emission, which is compared with the observations through spectra and intensity cuts. The amount of gas inside the cavity is quantified by varying the delta-gas parameter. Model fits to the dust and gas indicate that gas is still present inside the dust cavity for all disks but at a reduced level. The gas surface density drops inside the cavity by at least a factor 10, whereas the dust density drops by at least a factor 1000. Disk masses are comparable with previous estimates from the literature, cavity radii are found to be smaller than in the 345 GHz SubMillimeter Array (SMA) data. The derived gas surface density profiles suggest clearing of the cavity by one or more companions in all cases, trapping the millimeter-sized dust at the edge of the cavity.
HII regions are particularly interesting because they can generate dense layers of gas and dust, elongated columns or pillars of gas pointing towards the ionizing sources, and cometary globules of dense gas, where triggered star formation can occur. Understanding the interplay between the ionizing radiation and the dense surrounding gas is very important to explain the origin of these peculiar structures, and hence to characterize triggered star formation. G46.5-0.2 (G46), a poorly studied galactic HII region located at about 4 kpc, is an excellent target to perform this kind of studies. Using public molecular data extracted from the Galactic Ring Survey (13CO J=1-0) and from the James Clerk Maxwell Telescope data archive (12CO, 13CO, C18O J=3-2, HCO+ and HCN J=4-3), and infrared data from the GLIMPSE and MIPSGAL surveys, we perform a complete study of G46, its molecular environment and the young stellar objects placed around it. We found that G46, probably excited by an O7V star, is located close to the edge of the GRSMC G046.34-00.21 molecular cloud. It presents a horse-shoe morphology opening in direction of the cloud. We observed a filamentary structure in the molecular gas likely related to G46 and not considerable molecular emission towards its open border. We found that about 10' towards the southwest of G46 there are some pillar-like features, shining at 8 um and pointing towards the HII region open border. We propose that the pillar-like features were carved and sculpted by the ionizing flux from G46. We found several young stellar objects likely embedded in the molecular cloud grouped in two main concentrations: one, closer to the G46 open border consisting of Class II type sources, and other one mostly composed by Class I type YSOs located just ahead the pillars-like features, strongly suggesting an age gradient in the YSOs distribution.
The chameleon gravity model postulates the existence of a scalar field that couples with matter to mediate a fifth force. If it exists, this fifth force would influence the hot X-ray emitting gas that fills the potential wells of galaxy clusters. However, it would not influence the weak lensing signal from clusters. Therefore, by comparing X-ray and weak lensing profiles, one can place upper limits on the strength of a fifth force. This technique has been attempted before using a single, nearby cluster (Coma, $z=0.02$, Terukina et al. 2014). In this paper we apply the technique to the stacked profiles of 58 clusters at higher redshifts ($0.1<z<1.2$ ), including 12 new to the literature. X-ray data are taken from the XMM Cluster Survey (XCS) and weak lensing data are taken from the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). Using a simultaneous multi-parameter MCMC analysis, we are able to put constraints on the two chameleon gravity parameters ($\beta$ and $\phi_{\infty}$). Like the Terukina et al. (2014) study, our fits are consistent with general relativity, i.e. they do not require a fifth force. In the special case of $f(R)$ gravity (when the value of $\beta$ is fixed to $\sqrt{1/6}$), we can set an upper limit on the background field amplitude today of $|f_{\rm{R0}}| < 6 \times 10^{-5}$ (95% CL). This is the same level of constraint as Terukina et al. (2014) and demonstrates the validity of the stacking technique. It is one of the strongest constraints to date on $|f_{\rm{R0}}|$ on cosmological scales. We hope to improve this constraint in future by extending the study to hundreds of clusters using using weak lensing data from the Dark Energy Survey.
We present a theoretical study on the metallicity dependence of the initial$-$to$-$final mass relation and its influence on white dwarf age determinations. We compute a grid of evolutionary sequences from the main sequence to $\sim 3\, 000$ K on the white dwarf cooling curve, passing through all intermediate stages. During the thermally-pulsing asymptotic giant branch no third dredge-up episodes are considered and thus the photospheric C/O ratio is below unity for sequences with metallicities larger than $Z=0.0001$. We consider initial metallicities from $Z=0.0001$ to $Z=0.04$, accounting for stellar populations in the galactic disk and halo, with initial masses below $\sim 3M_{\odot}$. We found a clear dependence of the shape of the initial$-$to$-$final mass relation with the progenitor metallicity, where metal rich progenitors result in less massive white dwarf remnants, due to an enhancement of the mass loss rates associated to high metallicity values. By comparing our theoretical computations with semi empirical data from globular and old open clusters, we found that the observed intrinsic mass spread can be accounted for by a set of initial$-$to$-$final mass relations characterized by different metallicity values. Also, we confirm that the lifetime spent before the white dwarf stage increases with metallicity. Finally, we estimate the mean mass at the top of the white dwarf cooling curve for three globular clusters NGC 6397, M4 and 47 Tuc, around $0.53 M_{\odot}$, characteristic of old stellar populations. However, we found different values for the progenitor mass, lower for the metal poor cluster, NGC 6397, and larger for the younger and metal rich cluster 47 Tuc, as expected from the metallicity dependence of the initial$-$to$-$final mass relation.
We consider the shape of the posterior distribution to be used when fitting cosmological models to power spectra measured from galaxy surveys. At very large scales, Gaussian posterior distributions in the power do not approximate the posterior distribution $\mathcal P_R$ we expect for a Gaussian density field $\delta_\mathbf{k}$, even if we vary the covariance matrix according to the model to be tested. We compare alternative posterior distributions with $\mathcal P_R$, both mode-by-mode and in terms of expected $f_\mathrm{NL}$-measurements. Marginalising over a Gaussian posterior distribution $\mathcal P_f$ with fixed covariance matrix yields a posterior mean value of $f_\mathrm{NL}$ which, for a data set with the characteristics of Euclid, will be underestimated by $\triangle f_\mathrm{NL}=0.4$, while for the data release 9 (DR9) of the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) it will be underestimated by $\triangle f_\mathrm{NL}=19.1$. The inverse cubic normal distribution ($\mathcal P_\mathrm{ICN}$) agrees very well with $\mathcal P_R$ at all scales and for all data sets, hence providing the same marginalised value. Adopting this likelihood function means that we do not require a different covariance matrix for each model to be tested: this dependence is absorbed into the functional form of the posterior. Thus, the computational burden of analysis is significantly reduced.
We explore the thermal properties of hot and dense matter using a model that reproduces the empirical properties of isospin symmetric and asymmetric bulk nuclear matter, optical model fits to nucleon-nucleus scattering data, heavy-ion flow data in the energy range 0.5-2 GeV/A, and the largest well-measured neutron star mass of 2 $\rm{M}_\odot$. Results of this model which incorporates finite range interactions through Yukawa type forces are contrasted with those of a zero-range Skyrme model that yields nearly identical zero-temperature properties at all densities for symmetric and asymmetric nucleonic matter and the maximum neutron star mass, but fails to account for heavy-ion flow data due to the lack of an appropriate momentum dependence in its mean field. Similarities and differences in the thermal state variables and the specific heats between the two models are highlighted. Checks of our exact numerical calculations are performed from formulas derived in the strongly degenerate and non-degenerate limits. Our studies of the thermal and adiabatic indices, and the speed of sound in hot and dense matter for conditions of relevance to core-collapse supernovae, the thermal evolution of neutron stars from their birth and mergers of compact binary stars reveal that substantial variations begin to occur at sub-saturation densities before asymptotic values are reached at supra-nuclear densities.
We have investigated the critical conditions required for an efficient steady propeller mechanism in the spin-down phases of magnetized neutron stars with optically thick accretion disks. We have shown through simple analytical calculations that: (1) the strength of the dipole field at the Alfven radius is not sufficient to sustain an efficient mass-outflow even when the magnetic dipole field lines rotate much faster than the escape speed, (2) in the spin-down phase, mass accretion onto the star could persist above a minimum disk mass-flow rate that is orders of magnitude lower than the rate corresponding to the transition between the spin-up and the spin-down states, (3) below this critical mass-flow rate, a steady propeller state could be established with a maximum inner disk radius about 25 times smaller than the Alfven radius. Our results indicate that only for spherical accretion, the inner disk radius is likely to approach the Alfven radius, and for all realistic cases, the accretion-propeller transition could take place at a mass-flow rate lower than the rate equating the Alfven and the co-rotation radii. Our results are consistent with the properties of the transitional millisecond pulsars which show transitions between the accretion powered X-ray pulsar and the rotational powered radio pulsar states, and emit X-ray pulses in the sub-luminous X-ray phases.
The excitation of a non-linear ion-wake by a train of ultra-relativistic plasmons is modeled and its use for a novel regime of positron acceleration is explored. Its channel-like structure is independent of the energy-source driving the bubble-shaped slowly-propagating high phase-velocity electron density waves. The back of the bubble electron compression sucks-in the ions and the space-charge within the bubble expels them, forming a near-void channel with on-axis and bubble-edge density-spikes. The channel-edge density-spike is driven radially outwards as a non-linear ion acoustic-wave by the wake electron thermal pressure. OSIRIS PIC simulations are used to study the ion-wake structure, its evolution and its use for positron acceleration.
We investigate the allowed range of reheating temperature values in light of the Planck 2015 results and the recent joint analysis of Cosmic Microwave Background (CMB) data from the BICEP2/Keck Array and Planck experiments, using monomial and binomial inflationary potentials. While the well studied $\phi^2$ inflationary potential is no longer favored by current CMB data, as well as $\phi^p$ with $p>2$, a $\phi^1$ potential and canonical reheating ($w_{re}=0$) provide a good fit to the CMB measurements. In this last case, we find that the Planck 2015 $68\%$ confidence limit upper bound on the spectral index, $n_s$, implies an upper bound on the reheating temperature of $T_{re}\lesssim 6\times 10^{10}\,{\rm GeV}$, and excludes instantaneous reheating. The low reheating temperatures allowed by this model open the possiblity that dark matter could be produced during the reheating period instead of when the Universe is radiation dominated, which could lead to very different predictions for the relic density and momentum distribution of WIMPs, sterile neutrinos, and axions. We also study binomial inflationary potentials and show the effects of a small departure from a $\phi^1$ potential. We find that as a subdominant $\phi^2$ term in the potential increases, first instantaneous reheating becomes allowed, and then the lowest possible reheating temperature of $T_{re}=4\,{\rm MeV}$ is excluded by the Planck 2015 $68\%$ confidence limit.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
We estimate the amount of the missing baryons detected by the Planck measurements of the kinetic Sunyaev-Zeldovich effect (kSZ) around members of the Central Galaxy Catalogue (CGC) from the seventh release of the Sloan Survey. We use two statistics yielding evidence for kSZ signal, namely the pairwise peculiar momentum and the correlation function of the kSZ temperature estimates and predicted line-of-sight peculiar velocities. We find that both statistics yield consistent measurements of the Thomson optical depth $\tau_{\rm T}$ in the range of $0.5$-$1.4\times 10^{-4}$ for angular apertures that, on average, correspond to a range of distances of $>1$ to almost $3$ virial radii from the centres of the CG host halos. We find that, for the larger apertures for which we still have significant ($2$-$2.5\,\sigma$) kSZ detection, the regions probed around CGs contain roughly half the total amount of baryons present in the cosmological volume sampled by the Sloan footprint at $z\simeq 0.12$. Furthermore, under the assumption that baryons trace the dark matter distribution, our $\tau_{\rm T}$ measurements are compatible with having detected all the missing baryons around the CGs. Finally, our kSZ measurements yield no evidence for a kSZ dipole on the positions the CGs, providing the strongest constraints on the local bulk flow at a distance of $350$ $h^{-1}$ Mpc (below $290$ km s$^{-1}$ at 95% confidence level), and adding further evidence for the Copernican principle of homogeneity.
We present an analysis of the optical line emission from nearby barred galaxies, and in particular look at the radial range occupied by the bar. In many cases this region is marked by what we term a 'star formation desert', with a marked deficit of HII regions in optical narrow-band H-alpha imaging. Here we present long-slit spectroscopy revealing that such regions do have line emission, but that it is low-level, spatially smooth and almost ubiquitous. The relative strengths of the H-alpha and the spectrally adjacent [NII] lines in the regions are completely discrepant from those associated with star formation regions, and more closely match expectations for 'LINER' regions. We quantify the total line emission from these extended, kpc-scale regions, and determine the spurious contribution it would make to the determined star formation rate of these galaxies if interpreted as normal H-alpha emission. We concur with previous studies that link this LINER emission to old stellar populations, e.g. post-asymptotic giant branch stars, and propose these strongly-barred early-type spirals as a prime location for further tests of such emission.
We propose that luminous transients, including novae and supernovae, can be used to detect the faintest galaxies in the universe. Beyond a few Mpc, dwarf galaxies with stellar masses $<10^6 M_{\odot}$ will likely be too faint and/or too low in surface brightness to be directly detected in upcoming large area ground-based photometric surveys. However, single epoch LSST photometry will be able to detect novae to distances of $\sim30$ Mpc and SNe to Gpc-scale distances. Depending on the form of the stellar mass-halo mass relation and the underlying star formation histories of low mass dwarfs, the expected nova rates will be a few to $\sim100$ yr$^{-1}$ and the expected SN rates (including both type Ia and core-collapse) will be $\sim10^2-10^4$ within the observable ($4\pi$ sr) volume. The transient rate associated with intrahalo stars will be comparably large, but these transients will be located close to bright galaxies, in contrast to the dwarfs, which should trace the underlying large scale structure of the cosmic web. Aggressive follow-up of hostless transients has the potential to uncover the predicted enormous population of low mass field dwarf galaxies.
What strange new worlds will our next-generation telescopes find?
Most present-day galaxies with stellar masses $\geq10^{11}$ solar masses show no ongoing star formation and are dense spheroids. Ten billion years ago, similarly massive galaxies were typically forming stars at rates of hundreds solar masses per year. It is debated how star formation ceased, on which timescales, and how this "quenching" relates to the emergence of dense spheroids. We measured stellar mass and star-formation rate surface density distributions in star-forming galaxies at redshift 2.2 with $\sim1$ kiloparsec resolution. We find that, in the most massive galaxies, star formation is quenched from the inside out, on timescales less than 1 billion years in the inner regions, up to a few billion years in the outer disks. These galaxies sustain high star-formation activity at large radii, while hosting fully grown and already quenched bulges in their cores.
The astrophysics community is considering plans for a variety of gamma-ray telescopes (including ACT and GRIPS) in the energy range 1--100 MeV, which can fill in the so-called "MeV gap" in current sensitivity. We investigate the utility of such detectors for the study of low-mass dark matter annihilation or decay. For annihilating (decaying) dark matter with a mass below about 140 MeV (280 MeV) and couplings to first generation quarks, the final states will be dominated by photons or neutral pions, producing striking signals in gamma-ray telescopes. We determine the sensitivity of future detectors to the kinematically allowed final states. In particular, we find that planned detectors can improve on current sensitivity to this class of models by up to a few orders of magnitude.
We report the alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations. The combination of these state-of-the-art hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($11 < log_{10}(M_{200}/ [h^{-1}M_\odot]) < 15$), a wide radial range ($-2.3 < log_{10}(r/[h^{-1}Mpc ]) < 1.3$) and redshifts $0 < z < 1$. The shape parameters of the dark matter, stellar and hot gas distributions follow qualitatively similar trends: they become more aspherical (and triaxial) with increasing halo mass, radius and redshift. We measure the misalignment of the baryonic components (hot gas and stars) of galaxies with their host halo as a function of halo mass, radius, redshift, and galaxy type (centrals vs satellites and early- vs late-type). Overall, galaxies align well with local distribution of the total (mostly dark) matter. However, the stellar distributions on galactic scales exhibit a median misalignment of about 45-50 degrees with respect to their host haloes. This misalignment is reduced to 25-30 degrees in the most massive haloes ($13 < log_{10}(M_{200}/ [h^{-1}M_\odot ]) < 15$). Half of the disc galaxies in the EAGLE simulations have a misalignment angle with respect to their host haloes larger than 40 degrees. We present fitting functions and tabulated values for the probability distribution of galaxy-halo misalignment to enable a straightforward inclusion of our results into models of galaxy formations based on purely collisionless N-body simulations.
A sterile neutrino of ~keV mass is a well motivated dark matter candidate. Its decay generates a X-ray line which offers a unique target for X-ray telescopes. For the first time, we use the Gamma-ray Burst Monitor (GBM) onboard the Fermi Gamma-Ray Space Telescope to search for sterile neutrino decay lines; our analysis covers the energy range 10-25 keV (sterile neutrino mass 20-50 keV), which is inaccessible to X-ray and gamma-ray satellites such as Chandra, Suzaku, XMM-Newton, and INTEGRAL. The extremely wide field of view of the GBM enables a large fraction of the Milky Way dark matter halo to be probed. After implementing careful data cuts, we obtain ~53 days of full sky observational data. We search for sterile neutrino decay lines in the energy spectrum, and find no significant signal. From this, we obtain upper limits on the sterile neutrino mixing angle as a function of mass. In the sterile neutrino mass range 25-40 keV, we improve upon previous upper limits by approximately an order of magnitude. Better understanding of detector and astrophysical backgrounds, as well as detector response, will further improve the sensitivity of a search with the GBM.
Mazeh, Holczer, and Shporer (2015) have presented an approach that can, in principle, use the derived transit timing variation (TTV) of some transiting planets observed by the Kepler mission to distinguish between prograde and retrograde motion of their orbits with respect to the rotation of their parent stars. The approach utilizes TTVs induced by spot-crossing events that occur when the transiting planet moves across a spot on the stellar surface, by looking for a correlation between the derived TTVs and the stellar brightness derivatives at the corresponding transits, even in data that can not resolve the spot-crossing events themselves. We present here the application of this approach to the Kepler KOIs, identifying nine systems where the photometric spot modulation is large enough and the transit timing accurate enough to allow detection of a TTV-brightness-slope correlation. Excluding KOI-1546, which has been found recently to be a stellar binary, we are left with eight hot-Jupiter systems with high sensitivity to the correlation detection. Five of those eight systems show highly significant prograde motion, including two confirmed planets (KOI-203.01 = Kepler-17b and KOI-217.01 = Kepler-71b) and three planetary candidates (KOI-883.01, KOI-895.01, and KOI-1074.01), while no system displays retrograde motion, consistent with the suggestion that planets orbiting cool stars have prograde motion. All five systems have derived impact parameter $\lesssim$0.5, and all systems with an impact parameter in that range show significant correlation, except KOI-3.01 (= Kepler-3b = HAT-P-11b) where the lack of a correlation is explained by its large stellar obliquity. Although our sample is small, these findings hint that stellar spots, or at least the larger ones, have a tendency to be located at a low latitude on the stellar disc, similar to the Sun.
We present the results of a blind exercise to test the recoverability of
stellar rotation and differential rotation in Kepler light curves. The
simulated light curves lasted 1000 days and included activity cycles, Sun-like
butterfly patterns, differential rotation and spot evolution. The range of
rotation periods, activity levels and spot lifetime were chosen to be
representative of the Kepler data of solar like stars. Of the 1000 simulated
light curves, 770 were injected into actual quiescent Kepler light curves to
simulate Kepler noise. The test also included five 1000-day segments of the
Sun's total irradiance variations at different points in the Sun's activity
cycle.
Five teams took part in the blind exercise, plus two teams who participated
after the content of the light curves had been released. The methods used
included Lomb-Scargle periodograms and variants thereof, auto-correlation
function, and wavelet-based analyses, plus spot modelling to search for
differential rotation. The results show that the `overall' period is well
recovered for stars exhibiting low and moderate activity levels. Most teams
reported values within 10% of the true value in 70% of the cases. There was,
however, little correlation between the reported and simulated values of the
differential rotation shear, suggesting that differential rotation studies
based on full-disk light curves alone need to be treated with caution, at least
for solar-type stars.
The simulated light curves and associated parameters are available online for
the community to test their own methods.
Fairall 51 is a polar-scattered Seyfert 1 galaxy, a type of active galaxies believed to represent a bridge between unobscured type-1 and obscured type-2 objects. Fairall 51 has shown complex and variable X-ray absorption but only little is known about its origin. In our research, we observed Fairall 51 with the X-ray satellite Suzaku in order to constrain a characteristic time-scale of its variability. We performed timing and spectral analysis of four observations separated by 1.5, 2 and 5.5 day intervals. We found that the 0.5-50 keV broadband X-ray spectra are dominated by a primary power-law emission (with the photon index ~ 2). This emission is affected by at least three absorbers with different ionisations (log(xi) ~ 1-4). The spectrum is further shaped by a reprocessed emission, possibly coming from two regions -- the accretion disc and a more distant scattering region. The accretion disc emission is smeared by the relativistic effects, from which we measured the spin of the black hole as a ~ 0.8 (+-0.2). We found that most of the spectral variability can be attributed to the least ionised absorber whose column density changed by a factor of two between the first (highest-flux) and the last (lowest-flux) observation. A week-long scale of the variability indicates that the absorber is located at the distance ~ 0.05 pc from the centre, i.e., in the Broad Line Region.
A literature compilation of nuclear cluster (NSC) masses is used to study the correlation between global and NSC properties. A comparison of observational data to the predictions of semi-analytical galaxy formation models places constraints on the co-evolution of NSCs, massive black holes (MBHs) and host galaxies. Both data and theoretical predictions show an increased scatter in the NSC scaling correlations at high galaxy masses, and we show that this is due to the progressively more efficient ejection of stars from NSCs caused by MBH binaries in more massive stellar spheroids. Our results provide a natural explanation of why in nucleated galaxies hosting a MBH, the ratio (M_NSC+M_MBH)/M_bulge (with M_bulge the host spheroid's mass) shows significantly less scatter than M_NSC/M_bulge, and suggest that the formation of MBHs and NSCs are not mutually exclusive, as also supported by observations of co-existing systems. Both MBHs and NSCs represent generic products of galaxy formation, with NSCs being destroyed or modified by the merger evolution of their companion MBHs.
We propose an anisotropic generalisation of the line correlation function (ALCF) to separate and quantify phase information in the large-scale structure of galaxies. The line correlation function probes the strictly non-linear regime of structure formation and since phase information drops out of the power spectrum, the line correlation function provides a complementary tool to commonly used techniques based on two-point statistics. Furthermore, it is independent of linear bias as well as the Gaussian variance on the modulus of the density field and thus may also prove to be advantageous compared to the bispectrum or similar higher-order statistics for certain cases. For future applications it is vital, though, to be able to account for observational effects that cause anisotropies in the distribution of galaxies. Based on a number of numerical studies, we find that our ALCF is well suited to accomplish this task and we demonstrate how the Alcock-Paczynski effect and kinematical redshift-space distortions can in principle be measured via the ALCF.
We fit 54,296 sparsely-sampled asteroid lightcurves in the Palomar Transient Factory to a combined rotation plus phase-function model. Each lightcurve consists of 20+ observations acquired in a single opposition. Using 805 asteroids in our sample that have reference periods in the literature, we find the reliability of our fitted periods is a complicated function of the period, amplitude, apparent magnitude and other attributes. Using the 805-asteroid ground-truth sample, we train an automated classifier to estimate (along with manual inspection) the validity of the remaining 53,000 fitted periods. By this method we find 9,033 of our lightcurves (of 8,300 unique asteroids) have reliable periods. Subsequent consideration of asteroids with multiple lightcurve fits indicate 4% contamination in these reliable periods. For 3,902 lightcurves with sufficient phase-angle coverage and either a reliably-fit period or low amplitude, we examine the distribution of several phase-function parameters, none of which are bimodal though all correlate with the bond albedo and with visible-band colors. Comparing the theoretical maximal spin rate of a fluid body with our amplitude versus spin-rate distribution suggests that, if held together only by self-gravity, most asteroids are in general less dense than 2 g/cm$^3$, while C types have a lower limit of between 1 and 2 g/cm$^3$, in agreement with previous density estimates. For 5-20km diameters, S types rotate faster and have lower amplitudes than C types. If both populations share the same angular momentum, this may indicate the two types' differing ability to deform under rotational stress. Lastly, we compare our absolute magnitudes and apparent-magnitude residuals to those of the Minor Planet Center's nominal $G=0.15$, rotation-neglecting model; our phase-function plus Fourier-series fitting reduces asteroid photometric RMS scatter by a factor of 3.
Direct collapse black holes (DCBH) have been proposed as a solution to the challenge of assembling supermassive black holes by $z>6$ to explain the bright quasars observed at this epoch. The formation of a DCBH seed with $\rm M_{BH}\sim10^{4-5}\rm M_{\odot}$, requires a pristine atomic-cooling halo to be illuminated by an external radiation field that is sufficiently strong to entirely suppress H$_{2}$ cooling in the halo. Many previous studies have attempted to constrain the critical specific intensity that is likely required to suppress H$_{2}$ cooling, denoted as $J_{\rm crit}$. However, these studies have typically assumed that the incident external radiation field can be modeled with a black-body spectrum. Under this assumption, it is possible to derive a {unique} value for $J_{\rm crit}$ that depends only on the temperature of the black-body. In this study we consider a more realistic spectral energy distribution (SED) for the external source of radiation that depends entirely on its star formation history and age. The rate of destruction of the species responsible for suppressing molecular hydrogen cooling depends on the detailed shape of the SED. Therefore the value of $J_{\rm crit}$ is tied to the shape of the incident SED of the nearest star-forming protogalaxy. We fit a parametric form to the rates of destruction of H$_2$ and H$^-$ that permit direct collapse. Owing to this, we find that $J_{\rm crit}$ is not a fixed threshold but can lie anywhere in the range $J_{\rm crit} \sim 0.5$--$10^{3}$, depending on the details of the source stellar population, and its distance from the atomic cooling halo.
Investigations of the environments of SNe allow statistical constraints to be made on progenitor properties. We review progress that has been made in this field. Pixel statistics using tracers of e.g. star formation within galaxies show differences in the explosion sites of, in particular SNe types II and Ibc (SNe II and SNe Ibc), suggesting differences in population ages. Of particular interest is that SNe Ic are significantly more associated with H-alpha emission than SNe Ib, implying shorter lifetimes for the former. In addition, such studies have shown that the interacting SNe IIn do not explode in regions containing the most massive stars, which suggests that at least a significant fraction of their progenitors arise from the lower end of the core-collapse SN mass range. Host HII region spectroscopy has been obtained for a significant number of core-collapse events, however definitive conclusions have to-date been elusive. Single stellar evolution models predict that the fraction of SNe Ibc to SNe II should increase with increasing metallicity, due to the dependence of mass-loss rates on progenitor metallicity. We present a meta-analysis of host HII region oxygen abundances for CC SNe. It is concluded that the SN II to SN Ibc ratio shows little variation with oxygen abundance, with only a suggestion that the ratio increases in the lowest bin. Radial distributions of different SNe are discussed, where a central excess of SNe Ibc has been observed within disturbed galaxy systems, which is difficult to ascribe to metallicity or selection effects. Environment studies are evolving to enable studies at higher spatial resolutions than previously possible, while in addition the advent of wide-field integral field unit instruments allows galaxy-wide spectral analyses which will provide fruitful results to this field. Some example contemporary results are shown in that direction.
The electronic structure of six low-lying electronic states of scandium hydride, $X\,{}^{1}\Sigma^+$, $a\,{}^{3}\Delta$, $b\,{}^{3}\Pi$, $A\,{}^{1}\Delta$ $c\,{}^{3}\Sigma^+$, and $B\,{}^{1}\Pi$, is studied using multi-reference configuration interaction as a function of bond length. Diagonal and off-diagonal dipole moment, spin-orbit coupling and electronic angular momentum curves are also computed. The results are benchmarked against experimental measurements and calculations on atomic scandium. The resulting curves are used to compute a line list of molecular ro-vibronic transitions for $^{45}$ScH.
An on-going effort in the characterization of exoplanetary systems is the accurate determination of host star properties. This effort extends to the relatively bright host stars of planets discovered with the radial velocity method. The Transit Ephemeris Refinement and Monitoring Survey (TERMS) is aiding in these efforts as part of its observational campaign for exoplanet host stars. One of the first known systems is that of 70 Virginis, which harbors a jovian planet in an eccentric orbit. Here we present a complete characterization of this system with a compilation of TERMS photometry, spectroscopy, and interferometry. We provide fundamental properties of the host star through direct interferometric measurements of the radius (1.5\% uncertainty) and through spectroscopic analysis. We combined 59 new Keck HIRES radial velocity measurements with the 169 previously published from the ELODIE, Hamilton, and HIRES spectrographs, to calculate a refined orbital solution and construct a transit ephemeris for the planet. These newly determined system characteristics are used to describe the Habitable Zone of the system with a discussion of possible additional planets and related stability simulations. Finally, we present 19 years of precision robotic photometry that constrain stellar activity and rule out central planetary transits for a Jupiter-radius planet at the 5$\sigma$ level, with reduced significance down to an impact parameter of $b = 0.95$.
A neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH) is one of the plausible candidates for the central engine of gamma-ray bursts (GRBs). Two mechanisms, i.e., Blandford-Znajek (BZ) mechanism and neutrino annihilation process, are generally considered to power GRBs. Using the analytic solutions from Xue et al. (2013) and ignoring the effects of the magnetic field configuration, we estimate the BZ and neutrino annihilation luminosities as the functions of the disk masses and BH spin parameters to contrast the observational jet luminosities of GRBs. The results show that, although the neutrino annihilation processes could account for most of GRBs, the BZ mechanism is more effective, especially for long-duration GRBs. Actually, if the energy of afterglows and flares of GRBs is included, the distinction between these two mechanisms is more significant. Furthermore, massive disk mass and high BH spin are beneficial to power high luminosities of GRBs. Finally, we discuss possible physical mechanisms to enhance the disk mass or the neutrino emission rate of NDAFs and relevant difference between these two mechanisms.
We present observations of the first ten degrees of longitude in the Mopra carbon monoxide (CO) survey of the southern Galactic plane (Burton et al. 2013), covering Galactic longitude l = 320-330{\deg} and latitude b = $\pm$0.5{\deg}, and l = 327-330{\deg}, b = +0.5-1.0{\deg}. These data have been taken at 35 arc sec spatial resolution and 0.1 km/s spectral resolution, providing an unprecedented view of the molecular clouds and gas of the southern Galactic plane in the 109-115 GHz J = 1-0 transitions of 12CO, 13CO, C18O and C17O. Together with information about the noise statistics from the Mopra telescope, these data can be retrieved from the Mopra CO website and the CSIRO-ATNF data archive.
We present HI observations of the Sculptor Group starburst spiral galaxy NGC 253, obtained with the Karoo Array Telescope (KAT-7). KAT-7 is a pathfinder for the SKA precursor MeerKAT, under construction. The short baselines and low system temperature of the telescope make it very sensitive to large scale, low surface brightness emission. The KAT-7 observations detected 33% more flux than previous VLA observations, mainly in the outer parts and in the halo for a total HI mass of $2.1 \pm 0.1$ $\times 10^{9}$ M$_{\odot}$. HI can be found at large distances perpendicular to the plane out to projected distances of ~9-10 kpc away from the nucleus and ~13-14 kpc at the edge of the disk. A novel technique, based on interactive profile fitting, was used to separate the main disk gas from the anomalous (halo) gas. The rotation curve (RC) derived for the HI disk confirms that it is declining in the outer parts, as seen in previous optical Fabry-Perot measurements. As for the anomalous component, its RC has a very shallow gradient in the inner parts and turns over at the same radius as the disk, kinematically lagging by ~100 km/sec. The kinematics of the observed extra planar gas is compatible with an outflow due to the central starburst and galactic fountains in the outer parts. However, the gas kinematics shows no evidence for inflow. Analysis of the near-IR WISE data, shows clearly that the star formation rate (SFR) is compatible with the starburst nature of NGC 253.
(Abridged) We detect the large-scale structure of Lya emission in the Universe at redshifts z=2-3.5 by measuring the cross-correlation of Lya surface brightness with quasars in SDSS/BOSS. We use a million spectra targeting Luminous Red Galaxies at z<0.8, after subtracting a best fit model galaxy spectrum from each one, as an estimate of the high-redshift Lya surface brightness. The quasar-Lya emission cross-correlation we detect has a shape consistent with a LambdaCDM model with Omega_M =0.30^+0.10-0.07. The predicted amplitude of this cross-correlation is proportional to the product of the mean Lya surface brightness, <mu_alpha>, the amplitude of mass fluctuations, and the quasar and Lya emission bias factors. Using known values, we infer <mu_alpha>(b_alpha/3) = (3.9 +/- 0.9) x 10^-21 erg/s cm^-2 A^-1 arcsec^-2, where b_alpha is the Lya emission bias factor. If the dominant sources of Lya emission are star forming galaxies, we infer rho_SFR = (0.28 +/- 0.07) (3/b_alpha) /yr/Mpc^3 at z=2-3.5. For b_alpha=3, this value is a factor of 21-35 above previous estimates from individually detected Lya emitters, although consistent with the total rho_SFR derived from dust-corrected, continuum UV surveys. 97% of the Lya emission in the Universe at these redshifts is therefore undetected in previous surveys of Lya emitters. Our measurement is much greater than seen from stacking analyses of faint halos surrounding previously detected Lya emitters, but we speculate that it arises from similar Lya halos surrounding all luminous star-forming galaxies. We also detect redshift space anisotropy of the quasar-Lya emission cross-correlation, finding evidence at the 3.0 sigma level that it is radially elongated, consistent with distortions caused by radiative-transfer effects (Zheng et al. (2011)). Our measurements represent the first application of the intensity mapping technique to optical observations.
We present the detection and characterization of the transiting warm Jupiter KOI-12b, first identified with Kepler with an orbital period of 17.86 days. We combine the analysis of Kepler photometry with Doppler spectroscopy and line-profile tomography of time-series spectra obtained with the SOPHIE spectrograph to establish its planetary nature and derive its properties. To derive reliable estimates for the uncertainties on the tomographic model parameters, we devised an empirical method to calculate statistically independent error bars on the time-series spectra. KOI-12b has a radius of 1.43$\pm$0.13$ R_\mathrm{Jup}$ and a 3$\sigma$ upper mass limit of 10$M_\mathrm{Jup}$. It orbits a fast-rotating star ($v$sin$i_{\star}$ = 60.0$\pm$0.9 km s$^{-1}$) with mass and radius of 1.45$\pm$0.09 $M_\mathrm{Sun}$ and 1.63$\pm$0.15 $R_\mathrm{Sun}$, located at 426$\pm$40 pc from the Earth. Doppler tomography allowed a higher precision on the obliquity to be reached by comparison with the analysis of the Rossiter-McLaughlin radial velocity anomaly, and we found that KOI-12b lies on a prograde, slightly misaligned orbit with a low sky-projected obliquity $\lambda$ = 12.6$\stackrel{+3.0}{_{-2.9}}^\circ$. The properties of this planetary system, with a 11.4 magnitude host-star, make of KOI-12b a precious target for future atmospheric characterization.
We use our Galactic Globular Cluster Catalog (G2C2) photometry for 111
Galactic globular clusters (GC) in g and z, as well as r and i photometry for a
subset of 60 GCs and u photometry for 22 GCs, to determine the structural
parameters assuming King (1962) models.
In general, the resulting core radii are in good comparison with the current
literature values. However, our half-light radii are slightly lower than the
literature. The concentrations (and therefore also the tidal radii) are poorly
constrained mostly because of the limited radial extent of our imaging.
Therefore, we extensively discuss the effects of a limited field-of-view on the
derived parameters using mosaicked SDSS data, which do not suffer from this
restriction. We also illustrate how red giant branch (RGB) stars in cluster
cores can stochastically induce artificial peaks in the surface brightness
profiles. The issues related to these bright stars are scrutinised based on
both our photometry and simulated clusters. We also examine colour gradients
and find that the strongest central colour gradients are caused by central RGB
stars and thus not representative for the cluster light or colour distribution.
We recover the known relation between the half-light radius and the
Galactocentric distance in the g-band, but find a lower slope for redder
filters. We did not find a correlation between the scatter on this relation and
other cluster properties. We find tentative evidence for a correlation between
the half-light radii and the [Fe/H], with metal-poor GCs being larger than
metal-rich GCs. However, we conclude that this trend is caused by the position
of the clusters in the Galaxy, with metal-rich clusters being more centrally
located.
We constrain radio source clustering towards $Planck$-selected galaxy clusters using the NVSS point source catalogue. The constraint can be utilised for generating realistic Sunyaev-Zeldovich effect (SZE) mocks, and for predicting detectable clusters count and quantifying source confusion in radio surveys.
The IDV observations of S4 0917+624 were carried out monthly from August 2005 to January 2010, with the Urumqi 25m radio telescope at 4.8 GHz. The quasar S4 0917+624 exhibits either no or only very weak IDV during our 4.5 year observing interval. Prior to year 2000, the source S4 0917+624 was one of the most prominent IDV sources. Our new data indicate that the previous strong IDV has ceased and is not recovered. We analysed the long term VLBI structural variability using Gaussian model-fitting. From this we obtained the flux densities and the deconvolved sizes of core and inner jet components of the source. We studied the properties such as core fraction, angular size, spectral index, and brightness temperature of VLBI core for S4 0917+624, as well as the time delay between 5 and 15 GHz variations, and compared with the IDV properties of S4 0917+624. The source shows ejection of several jet components that are suspected to have partially reduced the IDV amplitude of S4 0917+624. However, during 2005-2006, the VLBI core size was comparable to the size before year 2000, but no strong IDV was detected in the period, suggesting that the quenching effect due to source size changes may not be responsible for the lack of strong IDV after year 2000. The refractive scattering properties for the strong IDV phase of S4 0917+624 before the year 2000 are discussed. The vanishing of strong IDV in S4 0917+624 after year 2000 is a mystery and cannot be explained via the quenching effect by changes in the observable VLBI structure. It however may be caused by changes in the interstellar medium, i.e. by interstellar weather, which induces changes in the scintillation pattern on timescales of several years. Further coordinated multi-frequency observations will be required to distinguish between the effect of source intrinsic variability and changing properties of the interstellar medium.
We argue that magnetic fields amplified within a very high accretion-rate disk around main sequence stars can lead to the formation of massive bipolar outflows that can remove most of the disk's mass and energy. This efficient directional removal of energy and mass allows the high accretion-rate disk to be built. We construct thick disks where the magnetic fields are amplified by an Alpha-Omega dynamo in the disk, bringing the fluctuating components of the magnetic field to be much stronger than the large-scale component. By examining the possible activity of the magnetic fields we conclude that main sequence stars can accrete mass at very high rates, up to 0.01Mo/yr for solar type stars, and up to 1Mo/yr for very massive stars. Such energetic outflows can account for the powering of some eruptive objects, such as merging main sequence stars, major eruptions of luminous blue variables, such as the Great Eruption of Eta Carinae, and other intermediate luminosity optical transients (ILOTs; Red Novae; Red Transients). Such powerful outflows can also supply the extra energy required in the common envelope process and in the grazing envelope evolution of binary systems where the companion is a main sequence star.
The current methods to determine the primary energy of ultra-high energy
cosmic rays (UHECRs) are different when dealing with hadron or photon
primaries. The current experiments combine two different techniques, an array
of surface detectors and fluorescence telescopes. The latter allow an almost
calorimetric measurement of the primary energy. Thus, hadron-initiated showers
detected by both type of detectors are used to calibrate the energy estimator
from the surface array (usually the interpolated signal at a certain distance
from the shower core S(r0)) with the primary energy. On the other hand, this
calibration is not feasible when searching for photon primaries since no high
energy photon has been unambiguously detected so far. Therefore, pure Monte
Carlo parametrizations are used instead.
In this work, we present a new method to determine the primary energy of
hadron-induced showers in a hybrid experiment based on a technique previously
developed for photon primaries. It consists on a set of calibration curves that
relate the surface energy estimator, S(r0), and the depth of maximum
development of the shower, Xmax, obtained from the fluorescence telescopes.
Then, the primary energy can be determined from pure surface information since
S(r0) and the zenith angle of the incoming shower are only needed. Considering
a mixed sample of ultra-high energy proton and iron primaries and taking into
account the reconstruction uncertainties and shower to shower fluctuations, we
demonstrate that the primary energy may be determined with a systematic
uncertainty below 1% and resolution around 16% in the energy range from
10^{18.5} to 10^{19.6} eV. Several array geometries, the shape of the energy
error distributions and the uncertainties due to the unknown composition of the
primary flux have been analyzed as well.
For the sake of complete theoretical research of atmospheric refraction, the atmospheric refraction under the condition of lower angles of elevation is still worthy to be analyzed and explored. In some engineering applications, the objects with larger zenith distance must be observed sometimes. Carrying out observational research of the atmospheric refraction at lower angles of elevation has an important significance. It has been considered difficult to measure the atmospheric refraction at lower angles of elevation. A new idea for determining atmospheric refraction by utilizing differential measurement with double fields of view is proposed. Taking the observational principle of HIPPARCOS satellite as a reference, a schematic prototype with double fields of view was developed. In August of 2013, experimental observations were carried out and the atmospheric refractions at lower angles of elevation can be obtained by the schematic prototype. The measured value of the atmospheric refraction at the zenith distance of 78.8 degree is $240.23"\pm0.27"$, and the feasibility of differential measurement of atmospheric refraction with double fields of view was justified. The limitations of the schematic prototype such as inadequate ability of gathering light, lack of accurate meteorological data recording and lower automatic level of observation and data processing were also pointed out, which need to be improved in subsequent work.
The main objective of the present project is to explore the viability of an adaptive optics control system based exclusively on Field Programmable Gate Arrays (FPGAs), making strong use of their parallel processing capability. In an Adaptive Optics (AO) system, the generation of the Deformable Mirror (DM) control voltages from the Wavefront Sensor (WFS) measurements is usually through the multiplication of the wavefront slopes with a predetermined reconstructor matrix. The ability to access several hundred hard multipliers and memories concurrently in an FPGA allows performance far beyond that of a modern CPU or GPU for tasks with a well defined structure such as Adaptive Optics control. The target of the current project is to generate a signal for a real time wavefront correction, from the signals coming from a Wavefront Sensor, wherein the system would be flexible to accommodate all the current Wavefront Sensing techniques and also the different methods which are used for wavefront compensation. The system should also accommodate for different data transmission protocols (like Ethernet, USB, IEEE 1394 etc.) for transmitting data to and from the FPGA device, thus providing a more flexible platform for Adaptive Optics control. Preliminary simulation results for the formulation of the platform, and a design of a fully scalable slope computer is presented.
In the analysis of thermal infrared data of asteroids by means of thermophysical models (TPMs) it is a common practice to neglect the uncertainty of the shape model and the rotational state, which are taken as an input for the model. Here, we present a novel method of investigating the importance of the shape model and the pole orientation uncertainties in the thermophysical modeling - the varied shape TPM (VS-TPM). Our method uses optical photometric data to generate various shape models that map the uncertainty in the shape and the rotational state. The TPM procedure is then run for all these shape models. We apply the implementation of the classical TPM as well as our VS-TPM to the convex shape models of several asteroids together with their thermal infrared data acquired by the NASA's Wide-field Infrared Survey Explorer (WISE) and compare the results. These show that the uncertainties of the shape model and the pole orientation can be very important (e.g., for the determination of the thermal inertia) and should be considered in the thermophysical analyses. We present thermophysical properties for six asteroids - (624) Hektor, (771) Libera, (1036) Ganymed, (1472) Muonio, (1627) Ivar, and (2606) Odessa.
We analyze the stellar absorption features in high signal-to-noise ratio near-infrared (NIR) spectra of the nuclear region of 12 nearby galaxies, mostly spirals. The features detected in some or all of the galaxies in this sample are the TiO (0.843 $\mu$m\ and 0.886 $\mu$m), VO (1.048 $\mu$m), CN (1.1 $\mu$m\ and 1.4 $\mu$m), H$\rm _2$O (1.4 $\mu$m\ and 1.9 $\mu$m) and CO (1.6 $\mu$m\ and 2.3 $\mu$m) bands. The C$\rm _2$ (1.17 $\mu$m\ and 1.76 $\mu$m) bands are generally weak or absent, although C$\rm _2$ (1.76 $\mu$m) may be weakly present in the mean galaxy spectrum. A deep feature near 0.93 $\mu$m, likely caused by CN, TiO and/or ZrO, is also detected in all objects. Fitting a combination of stellar spectra to the mean spectrum shows that the absorption features are produced by evolved stars: cool giants and supergiant stars in the early- or thermally-pulsing asymptotic giant branch (E-AGB or TP-AGB) phases. The high luminosity of TP-AGB stars, and the appearance of VO and ZrO features in the data, suggest that TP-AGB stars dominate these spectral features. However, a contribution from other evolved stars is also likely. Comparison with evolutionary population synthesis models shows that models based on empirical libraries that predict relatively strong NIR features provide a more accurate description of the data. However, none of the models tested accurately reproduces all of the features observed in the spectra. To do so, the models will need to not only improve the treatment of TP-AGB stars, but also include good quality spectra of red giant and E-AGB stars. The uninterrupted wavelength coverage, high S/N, and quantity of features we present here will provide a benchmark for the next generation of models aiming to explain and predict the NIR properties of galaxies.
The results of a search for peculiar astronomical objects using very low resolution spectra obtained with the NASA Orbital Debris Observatory (NODO) 3 meter diameter liquid mirror telescope (LMT) are compared with results of spectra obtained with the Sloan Digital Sky Survey (SDSS). The main purpose of this comparison is to verify whether observations taken with this novel type of telescope are reliable. This comparison is important because LMTs are a novel type of inexpensive telescope that are very useful for astronomical surveys, particularly surveys in the time domain, and a validation of data taken with an LMT, by comparison with data from a classical telescope, will validate their reliability. We start from a published data analysis that classified only 206 of the 18,000 astronomical objects observed with the NODO liquid mirror telescope as peculiar. A total of 29 of these 206 objects were found in the SDSS. The reliability of the NODO data can be seen by the results of the detailed analysis that, in practice, less than 0.3% of the 18,000 spectra were incorrectly identified as peculiar objects, most probably because they are variable stars. We conclude that the liquid mirror telescope gave reliable observations, comparable to those that would have been obtained with a telescope using a glass mirror.
There are many Slowly Pulsating B (SPB) stars and gamma Dor stars in the Kepler Mission data set. The light curves of these pulsating stars have been classified phenomenologically into stars with symmetric light curves and with asymmetric light curves. In the same effective temperature ranges as the gamma Dor and SPB stars, there are variable stars with downward light curves that have been conjectured to be caused by spots. Among these phenomenological classes of stars, some show `frequency groups' in their amplitude spectra that have not previously been understood. While it has been recognised that nonlinear pulsation gives rise to combination frequencies in a Fourier description of the light curves of these stars, such combination frequencies have been considered to be a only a minor constituent of the amplitude spectra. In this paper we unify the Fourier description of the light curves of these groups of stars, showing that many of them can be understood in terms of only a few base frequencies, which we attribute to g mode pulsations, and combination frequencies, where sometimes a very large number of combination frequencies dominate the amplitude spectra. The frequency groups seen in these stars are thus tremendously simplified. We show observationally that the combination frequencies can have amplitudes greater than the base frequency amplitudes, and we show theoretically how this arises. Thus for some gamma Dor and SPB stars combination frequencies can have the highest observed amplitudes. Among the B stars are pulsating Be stars that show emission lines in their spectra from occasional ejection of material into a circumstellar disk. Our analysis gives strong support to the understanding of these pulsating Be stars as rapidly rotating SPB stars, explained entirely by g mode pulsations.
The IACOB project is an ambitious long-term project which is contributing to step forward in our knowledge about the physical properties and evolution of Galactic massive stars. The project aims at building a large database of high-resolution, multi-epoch, spectra of Galactic OB stars, and the scientific exploitation of the database using state-of-the-art models and techniques. In this proceeding, we summarize the latest updates of the IACOB spectroscopic database and highlight some of the first scientific results from the IACOB project; we also announce the first data release and the first public version of the iacob-broad tool for the line-broadening characterization of OB-type spectra.
Context. Stellar bow shocks have been studied not only observationally, but also theoretically since the late 1980s. Only a few catalogues of them exist. The bow shocks show emission along all the electromagnetic spectrum, but they are detected more easily in infrared wavelengths. The release of new and high-quality infrared data eases the discovery and subsequent study of new objects. Aims. We search stellar bow-shock candidates associated with nearby runaway stars, and gather them together with those found elsewhere, to enlarge the list of the E-BOSS first release. We aim to characterize the bow-shock candidates and provide a database suitable for statistical studies. We investigate the low-frequency radio emission at the position of the bow-shock features, that can contribute to further studies of high-energy emission from these objects. Methods. We considered samples from different literature sources and searched for bow-shaped structures associated with stars in the Wide-field Infrared Survey Explorer (WISE) images. We looked for each bow-shock candidate on centimeter radio surveys. Results. We reunited 45 bow-shock candidates and generated composed WISE images to show the emission in different infrared bands. Among them there are new sources, previously studied objects, and bow shocks found serendipitously. Five bow shocks show evidence of radio emission. Conclusions. Stellar bow shocks constitute an active field with open questions and enormous amounts of data to be analyzed. Future research at all wavelengths databases, and use of instruments like Gaia, will provide a more complete picture of these objects. For instance, infrared spectral energy distributions can give information about physical parameters of the bow shock matter. In addition, dedicated high-sensitivity radio observations can help to understand the radio-$\gamma$ connection.
The optical luminous quasar PG0043+039 has not been detected before in deep X-ray observations indicating the most extreme optical-to-X-ray slope index ${\alpha}_{ox}$ of all quasars. This study aims to detect PG0043+039 in a deep X-ray exposure. Furthermore, we wanted to check out whether this object shows specific spectral properties in other frequency bands. We took deep X-ray (XMM-Newton), far-ultraviolet (HST), and optical (HET, SALT telescopes) spectra of PG0043+039 simultaneously in July 2013. We just detected PG0043+039 in our deep X-ray exposure. The steep ${\alpha}_{ox} = -2.37 {\pm} 0.05$ gradient is consistent with an unusual steep gradient $F_{\nu} {\sim} {\nu}^{\alpha}$ with ${\alpha} = -2.67 {\pm} 0.02$ seen in the UV/far-UV continuum. The optical/UV continuum flux has a clear maximum near 2500 {\AA}. The UV spectrum is very peculiar because it shows broad humps in addition to known emission lines. A modeling of these observed humps with cyclotron lines can explain their wavelength positions, their relative distances, and their relative intensities. We derive plasma temperatures of T ${\sim}$ 3keV and magnetic field strengths of B ${\sim}$ 2 ${\times} 10^8$ G for the line-emitting regions close to the black hole.
Polarization of radio emission in extragalactic jets at a sub-milliarcsecond angular resolution holds important clues for understanding the structure of the magnetic field in the inner regions of the jets and in close vicinity of the supermassive black holes in the centers of active galaxies. Space VLBI observations provide a unique tool for polarimetric imaging at a sub-milliarcsecond angular resolution and studying the properties of magnetic field in active galactic nuclei on scales of less than 10^4 gravitational radii. A space VLBI observation of high-redshift quasar TXS 0642+449 (OH 471), made at a wavelength of 18 cm (frequency of 1.6 GHz) as part of the Early Science Programme (ESP) of the RadioAstron} mission, is used here to test the polarimetric performance of the orbiting Space Radio Telescope (SRT) employed by the mission, to establish a methodology for making full Stokes polarimetry with space VLBI at 1.6 GHz, and to study the polarized emission in the target object on sub-milliarcsecond scales. Polarization leakage of the SRT at 18 cm is found to be within 9 percents in amplitude, demonstrating the feasibility of high fidelity polarization imaging with RadioAstron at this wavelength. A polarimetric image of 0642+449 with a resolution of 0.8 mas (signifying an ~4 times improvement over ground VLBI observations at the same wavelength) is obtained. The image shows a compact core-jet structure with low (~2%) polarization and predominantly transverse magnetic field in the nuclear region. The VLBI data also uncover a complex structure of the nuclear region, with two prominent features possibly corresponding to the jet base and a strong recollimation shock. The maximum brightness temperature at the jet base can be as high as 4*10^13 K.
Using the updated proton and helium fluxes just released by the AMS-02 experiment we reevaluate the secondary astrophysical antiproton to proton ratio and its uncertainties, and compare it with the ratio preliminarly reported by AMS-02. We find no unambiguous evidence for a significant excess with respect to expectations. Yet, some preference for thicker halos and a flatter energy dependence of the diffusion coefficient starts to emerge. Also, we provide an assessment of the room left for exotic components such as Galactic Dark Matter annihilation or decay, deriving new stringent constraints.
The beginning of photoionization marks the transition between the post-Asymptotic Giant Branch (post-AGB) and planetary nebula (PN) phases of stars with masses < 8 M_sun. This critical phase is difficult to observe, as it lasts only a few decades. The combination of jets and magnetic fields, the key agents of PNe shaping, could give rise to synchrotron emission, but this has never been observed before in any PNe, since free-free emission from the ionized gas is expected to dominate its radio spectrum. In this paper we report radio continuum observations taken with the Australia Telescope Compact Array between 1 and 46 GHz of the young PN IRAS 15103-5754. Our observations in 2010-2011 show non-thermal emission compatible with synchrotron emission from electrons accelerated at a shock with spectral index $\alpha \simeq -0.54$. However, in 2012, the spectral index $\alpha \simeq -0.28$ is no longer compatible with synchrotron emission in these types of processes. Several hypothesis are discussed to explain this change. The more plausible ones are related to the presence of the newly photoionized region in this young PN: either energy loss of electrons due to Coulomb collisions with the plasma, or selective suppression of synchrotron radiation due to the Razin effect. We postulate that the observed flattening of non-thermal radio spectra could be a hallmark identifying the beginning of the PN phase.
In the study of gravitational waves (GWs), the stochastic background generated by compact binary systems are among the most important kinds of signals. The reason for such an importance has to do with their probable detection by the interferometric detectors [such as the Advanced LIGO (ALIGO) and Einstein Telescope (ET)] in the near future. In this paper we are concerned with, in particular, the stochastic background of GWs generated by double neutron star (DNS) systems in circular orbits during their periodic and quasi--periodic phases. Our aim here is to describe a new method to calculate such spectra, which is based on an analogy with a problem of Statistical Mechanics. Besides, an important characteristic of our method is to consider the time evolution of the orbital parameters.
We present and test TesseRACt, a non-parametric technique for recovering the concentration of simulated dark matter halos using Voronoi tessellation. TesseRACt is tested on idealized N-body halos that are axisymmetric, triaxial, and contain substructure and compared to traditional least-squares fitting as well as two non-parametric techniques that assume spherical symmetry. TesseRACt recovers halo concentrations within 0.3% of the true value regardless of whether the halo is spherical, axisymmetric, or triaxial. Traditional fitting and non-parametric techniques that assume spherical symmetry can return concentrations that are systematically off by as much as 10% from the true value for non-spherical halos. TesseRACt also performs significantly better when there is substructure present outside $0.5R_{200}$. Given that cosmological halos are rarely spherical and often contain substructure, we discuss implications for studies of halo concentration in cosmological N-body simulations including how choice of technique for measuring concentration might bias scaling relations.
We investigate the physical conditions of the sources of two metric Type-II bursts associated with CME expansions with the aim of verifying the relationship between the shocks and the CMEs, comparing the heights of the radio sources and the heights of the EUV waves associated with the CMEs. The heights of the EUV waves associated with the events were determined in relation to the wave fronts. The heights of the shocks were estimated by applying two different density models to the frequencies of the Type-II emissions and compared with the heights of the EUV waves. For the 13 June 2010 event, with band-splitting, the shock speed was estimated from the frequency drifts of the upper and lower branches of the harmonic lane, taking into account the H/F frequency ratio fH/fF = 2. Exponential fits on the intensity maxima of the branches revealed to be more consistent with the morphology of the spectrum of this event. For the 6 June 2012 event, with no band-splitting and with a clear fundamental lane on the spectrum, the shock speed was estimated directly from the frequency drift of the fundamental emission, determined by linear fit on the intensity maxima of the lane. For each event, the most appropriate density model was adopted to estimate the physical parameters of the radio source. The 13 June 2010 event presented a shock speed of 664-719 km/s, consistent with the average speed of the EUV wave fronts of 609 km/s. The 6 June 2012 event was related to a shock of speed of 211-461 km/s, also consistent with the average speed of the EUV wave fronts of 418 km/s. For both events, the heights of the EUV wave revealed to be compatible with the heights of the radio source, assuming a radial propagation of the shock.
We present a detailed study of second-order matter perturbations for the general Horn- deski class of models. Being the most general scalar-tensor theory having second-order equations of motion, it includes many known gravity and dark energy theories and General Relativity with a cosmological constant as a specific case. This enables us to estimate the leading order dark matter bispectrum generated at late-times by gravitational instability. We parametrize the evolution of the first and second-order equations of motion as proposed by Bellini and Sawicki (2014), where the free functions of the theory are assumed to be proportional to the dark energy density. We show that it is unnatural to have large 10% ( 1%) deviations of the bispectrum introducing even larger ~ 30% (~ 5%) deviations in the linear growth rate. Considering that measurements of the linear growth rate have much higher signal-to-noise than bispectrum measurements, this indicates that for Horndeski models which reproduce the expansion history and the linear growth rate as predicted by GR the dark matter bispectrum kernel can be effectively modelled as the standard GR one. On the other hand, an observation of a large bispectrum deviation that can not be explained in terms of bias would imply either that the evolution of perturbations is strongly different than the evolution predicted by GR or that the theory of gravity is exotic (e.g., breaks the weak equivalence principle) and/or fine-tuned.
O and early B stars are at the apex of galactic ecology, but in the Milky Way, only a minority of them may yet have been identified. We present the results of a pilot study to select and parametrise OB star candidates in the Southern Galactic plane, down to a limiting magnitude of $g=20$. A 2 square-degree field capturing the Carina Arm around the young massive star cluster, Westerlund 2, is examined. The confirmed OB stars in this cluster are used to validate our identification method, based on selection from the $(u-g, g-r)$ diagram for the region. Our Markov Chain Monte Carlo fitting method combines VPHAS+ $u, g, r, i$ with published $J, H, K$ photometry in order to derive posterior probability distributions of the stellar parameters $\log(\rm T_{\rm eff})$ and distance modulus, together with the reddening parameters $A_0$ and $R_V$. The stellar parameters are sufficient to confirm OB status while the reddening parameters are determined to a precision of $\sigma(A_0)\sim0.09$ and $\sigma(R_V)\sim0.08$. There are 489 objects that fit well as new OB candidates, earlier than $\sim$B2. This total includes 74 probable massive O stars, 5 likely blue supergiants and 32 reddened subdwarfs. This increases the number of previously known and candidate OB stars in the region by nearly a factor of 10. Most of the new objects are likely to be at distances between 3 and 6 kpc. We have confirmed the results of previous studies that, at these longer distances, these sight lines require non-standard reddening laws with $3.5<R_V<4$.
We study the evolution of planetesimals in evolved gaseous disks, which orbit a solar-mass star and harbor a Jupiter-mass planet at a_p~5AU. The gas dynamics is modeled with a three-dimensional hydrodynamics code that employes nested-grids and achieves a resolution of one Jupiter's radius in the circumplanetary disk. The code models solids as individual particles. Planetesimals are subjected to gravitational forces by the star and the planet, drag force by the gas, disruption via ram pressure, and mass loss through ablation. The mass evolution of solids is calculated self-consistently with their temperature, velocity, and position. We consider icy and icy/rocky bodies of radius 0.1-100km, initially deployed on orbits around the star within a few Hill radii (Rhill) of the planet's orbit. Planetesimals are scattered inward, outward, and toward disk regions of radius r>>a_p. Scattering can relocate significant amounts of solids, provided that regions |r-a_p|~ 3Rhill are replenished with planetesimals. Scattered bodies can be temporarily captured on planetocentric orbits. Ablation consumes nearly all solids at gas temperatures > ~220K. Super-keplerian rotation around and beyond the outer edge of the gas gap can segregate < ~0.1km bodies, producing solid gap edges at size-dependent radial locations. Capture, break-up, and ablation of solids result in a dust-laden circumplanetary disk with low surface densities of km-size planetesimals, implying relatively long timescales for satellite formation. After a giant planet acquires most of its mass, accretion of solids is unlikely to alter significantly its heavy-element content. The luminosity generated by solids' accretion can be of a similar order of magnitude to the contraction luminosity.
We study the dominant effect of a long wavelength density perturbation $\delta(\lambda_L)$ on short distance physics. In the non-relativistic limit, the result is a uniform acceleration, fixed by the equivalence principle, and typically has no effect on statistical averages due to translational invariance. This same reasoning has been formalized to obtain a "consistency condition" on the cosmological correlation functions. In the presence of a feature, such as the acoustic peak at $l_{\rm BAO}$, this naive expectation breaks down for $\lambda_L<l_{\rm BAO}$. We calculate a universal piece of the three-point correlation function in this regime. The same effect is shown to underlie the spread of the acoustic peak, and is calculable to all orders in the long modes. This can be used to improve the result of perturbative calculations - a technique known as "infra-red resummation"- and is explicitly applied to the one-loop calculation of power spectrum. Finally, the success of BAO reconstruction schemes is argued to be another empirical evidence for the validity of the results.
It is well-known that the SU(2) Reissner-Nordstr\"om black-hole solutions of the Einstein-Yang-Mills theory are characterized by an infinite set of unstable (imaginary) eigenvalues $\{\omega_n(T_{\text{BH}})\}_{n=0}^{n=\infty}$ (here $T_{\text{BH}}$ is the black-hole temperature). In this paper we analyze the excited instability spectrum of these magnetically charged black holes. The numerical results suggest the existence of a universal behavior for these black-hole excited eigenvalues. In particular, we show that unstable eigenvalues in the regime $\omega_n\ll T_{\text{BH}}$ are characterized, to a very good degree of accuracy, by the simple universal relation $\omega_n(r_+-r_-)={\text{constant}}$, where $r_{\pm}$ are the horizon radii of the black hole.
Pulsars apparently missing from the galactic center could have been destroyed by asymmetric fermionic dark matter ($m_X = 1-100$ GeV) coupled to a light scalar ($m_{\phi}= 5-20$ MeV), which mixes with the Higgs boson. We point out that this pulsar-collapsing dark sector can resolve the core-cusp problem and will either be excluded or discovered by upcoming direct detection experiments. Another implication is a maximum pulsar age curve that increases with distance from the galactic center, with a normalization that depends on the couplings and masses of dark sector particles. In addition, we use old pulsars outside the galactic center to place bounds on asymmetric Higgs portal models.
In this paper we study a real scalar field as a possible candidate to explain the dark matter in the universe. In the context of a free scalar field with quadratic potential, we have used observational $H(z)$ data to constrain the dark matter mass to $m=\left(3.46^{+0.38+0.75+1.1}_{-0.43-0.92-1.5}\right)\times10^{-33}$ eV. This value is much below some previous estimates of $m\sim 10^{-22}$ eV found in some models, which we explain as being due to a slightly different formulation, but in complete agreement with a recent model based on a cosmological scalar field harmonic oscillator, for which $m\sim 10^{-32}$ eV. Although scalar field dark matter (SFDM) is much disfavored, as it gives rise to ultra hot dark matter and could halt structure formation, different scalar field potentials could alleviate this issue.
A novel model of spontaneous Leptogenesis is investigated, where it takes place in the thermal equilibrium due to a background Nambu-Goldstone field in motion. In particular, we identify the Nambu-Goldstone field to be the Majoron which associates with spontaneous breakdown of (discrete) $B-L$ symmetry. In this scenario sufficient lepton number asymmetry is generated in primordial thermal bath without having $CP$-violating out-of-equilibrium decay of the heavy right-handed Majorana neutrinos. To obtain the observed baryon asymmetry, the neutrino masses are predicted in certain ranges, which can be translated into the effective mass of the neutrinoless double beta decay.
We propose a QCD axion model that avoids the cosmological domain wall problem, introducing a global SU(3)_f family symmetry to which we embed the unwanted PQ discrete symmetry. The spontaneous breaking of SU(3)_f and PQ symmetry predicts eight NG bosons as well as axion, all of which contribute to dark radiation in the Universe. The derivation from the standard model prediction of dark radiation can be observed by future observations of CMB fluctuations. Our model also predicts a sizable exotic kaon decay rate, which is marginally consistent with the present collider data and would be tested by future collider experiments.
Quantum gravity corrections to accretion onto a Schwarzschild black hole are considered in the context of asymptotically safe scenario. The possible positions of the critical points are discussed and the general conditions for critical points are obtained. The explicit expressions for matter density compression and temperature profile both below the critical radius and at the event horizon are derived. For polytropic matter, we determine the corrected temperature and the integrated flux resulting from quantum gravity effects at the event horizon, which may can be as a test of asymptotically safe scenario.
The dynamics of spinning bodies in General Relativity is studied in the test-mass limit. Equations of motion are obtained both by a hamiltonian construction and by energy-momentum conservation using generalized one-particle energy-momentum tensors. The latter approach also allows the computation of gravitational perturbations created by rotating compact objects.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)