Understanding the surface and atmospheric conditions of Earth-size, rocky planets in the habitable zones (HZs) of low-mass stars is currently one of the greatest astronomical endeavors. Knowledge of the planetary effective surface temperature alone is insufficient to accurately interpret biosignature gases when they are observed in the coming decades. The UV stellar spectrum drives and regulates the upper atmospheric heating and chemistry on Earth-like planets, is critical to the definition and interpretation of biosignature gases, and may even produce false-positives in our search for biologic activity. This white paper briefly describes the scientific motivation for panchromatic observations of exoplanetary systems as a whole (star and planet), argues that a future NASA UV/Vis/near-IR space observatory is well-suited to carry out this work, and describes technology development goals that can be achieved in the next decade to support the development of a UV/Vis/near-IR flagship mission in the 2020s.
We interpret recent ALMA observations of z > 6 normal star forming galaxies by means of a semi-numerical method, which couples the output of a cosmological hydrodynamical simulation with a chemical evolution model which accounts for the contribution to dust enrichment from supernovae, asymptotic giant branch stars and grain growth in the interstellar medium. We find that while stellar sources dominate the dust mass of small galaxies, the higher level of metal enrichment experienced by galaxies with Mstar > 10^9 Msun allows efficient grain growth, which provides the dominant contribution to the dust mass. Even assuming maximally efficient supernova dust production, the observed dust mass of the z = 7.5 galaxy A1689-zD1 requires very efficient grain growth. This, in turn, implies that in this galaxy the average density of the cold and dense gas, where grain growth occurs, is comparable to that inferred from observations of QSO host galaxies at similar redshifts. Although plausible, the upper limits on the dust continuum emission of galaxies at 6.5 < z < 7.5 show that these conditions must not apply to the bulk of the high redshift galaxy population
One possible channel for black hole formation is the collapse of a rigidly rotating massive neutron star as it loses its angular momentum or gains excessive mass through accretion. It was proposed that part of the neutron star may form a debris disk around the black hole. Such short-lived massive disks could be the sources of powerful jets emitting cosmological gamma-ray bursts. Whether the collapse creates a disk depends on the equation of state of the neutron star. We survey a wide range of equations of states allowed by observations and find that disk formation is unfeasible. We conclude that this channel of black hole formation is incapable of producing powerful jets, and discuss implications for models of gamma-ray bursts.
Kepler-296 is a binary star system with two M-dwarf components separated by 0.2 arcsec. Five transiting planets have been confirmed to be associated with the Kepler-296 system; given the evidence to date, however, the planets could in principle orbit either star. This ambiguity has made it difficult to constrain both the orbital and physical properties of the planets. Using both statistical and analytical arguments, this paper shows that all five planets are highly likely to orbit the primary star in this system. We performed a Markov-Chain Monte Carlo simulation using a five transiting planet model, leaving the stellar density and dilution with uniform priors. Using importance sampling, we compared the model probabilities under the priors of the planets orbiting either the brighter or the fainter component of the binary. A model where the planets orbit the brighter component, Kepler-296A, is strongly preferred by the data. Combined with our assertion that all five planets orbit the same star, the two outer planets in the system, Kepler-296 Ae and Kepler-296 Af, have radii of 1.53 +/- 0.26 and 1.80 +/- 0.31 R_earth, respectively, and receive incident stellar fluxes of 1.40 +/- 0.23 and 0.62 +/- 0.10 times the incident flux the Earth receives from the Sun. This level of irradiation places both planets within or close to the circumstellar habitable zone of their parent star.
We present a new analysis of the ionizing emissivity ($\dot{N}_{\rm{ion}}$, s$^{-1}$ Mpc$^{-3}$) for galaxies during the epoch of reionization and their potential for completing and maintaining reionization. We use extensive SED modelling -- incorporating two plausible mechanisms for the escape of Lyman continuum photon -- to explore the range and evolution of ionizing efficiencies consistent with new results on galaxy colours ($\beta$) during this epoch. We estimate $\dot{N}_{\rm{ion}}$ for the latest observations of the luminosity and star-formation rate density at $z<10$, outlining the range of emissivity histories consistent with our new model. Given the growing observational evidence for a UV colour-magnitude relation in high-redshift galaxies, we find that for any plausible evolution in galaxy properties, red (brighter) galaxies are less efficient at producing ionizing photons than their blue (fainter) counterparts. The assumption of a redshift and luminosity evolution in $\beta$ leads to two important conclusions. Firstly, the ionizing efficiency of galaxies naturally increases with redshift. Secondly, for a luminosity dependent ionizing efficiency, we find that galaxies down to a rest-frame magnitude of $M_{\rm{UV}} \approx -15$ alone can potentially produce sufficient numbers of ionizing photons to maintain reionization as early as $z\sim8$ for a clumping factor of $C_{\rm{H {\small II}}} \leq 3$.
We present an analysis of the binary and physical parameters of a unique pulsating white dwarf with a main-sequence companion, SDSS J1136+0409, observed for more than 77 d during the first pointing of the extended Kepler mission: K2 Campaign 1. Using new ground-based spectroscopy, we show that this post-common-envelope binary has an orbital period of 6.89760103(60) hr, which is also seen in the photometry as a result of Doppler beaming and ellipsoidal variations of the secondary. We spectroscopically refine the temperature of the white dwarf to 12330(260) K and its mass to 0.601(36) Msun. We detect seven independent pulsation modes in the K2 light curve. A preliminary asteroseismic solution is in reasonable agreement with the spectroscopic atmospheric parameters. Three of the pulsation modes are clearly rotationally split multiplets, which we use to demonstrate that the white dwarf is not synchronously rotating with the orbital period but has a rotation period of 2.49(53) hr. This is faster than any known isolated white dwarf, but slower than almost all white dwarfs measured in non-magnetic cataclysmic variables, the likely future state of this binary.
We present spatially and spectrally resolved Atacama Large Millimeter/submillimeter Array (ALMA) observations of gas and dust in the disk orbiting the pre-main sequence binary AK Sco. By forward-modeling the disk velocity field traced by CO J=2-1 line emission, we infer the mass of the central binary, $M_\ast = 2.49 \pm 0.10~M_\odot$, a new dynamical measurement that is independent of stellar evolutionary models. Assuming the disk and binary are co-planar within $\sim$2{\deg}, this disk-based binary mass measurement is in excellent agreement with constraints from radial velocity monitoring of the combined stellar spectra. These ALMA results are also compared with the standard approach of estimating masses from the location of the binary in the Hertzsprung-Russell diagram, using several common pre-main sequence model grids. These models predict stellar masses that are marginally consistent with our dynamical measurement (at $\sim 2\,\sigma$), but are systematically high (by $\sim$10%). These same models consistently predict an age of $18\pm1$ Myr for AK Sco, in line with its membership in the Upper Centaurus-Lupus association but surprisingly old for it to still host a gas-rich disk. As ALMA accumulates comparable data for large samples of pre-main sequence stars, the methodology employed here to extract a dynamical mass from the disk rotation curve should prove extraordinarily useful for efforts to characterize the fundamental parameters of early stellar evolution.
The discovery of over 50 planets around evolved stars and more than 35 debris discs orbiting white dwarfs highlight the increasing need to understand small body evolution around both early and asymptotic giant branch (GB) stars. Pebbles and asteroids are susceptible to strong accelerations from the intense luminosity and winds of GB stars. Here, we establish equations that can model time-varying GB stellar radiation, wind drag and mass loss. We derive the complete three-dimensional equations of motion in orbital elements due to (1) the Epstein and Stokes regimes of stellar wind drag, (2) Poynting-Robertson drag, and (3) the Yarkovsky drift with seasonal and diurnal components. We prove through averaging that the potential secular eccentricity and inclination excitation due to Yarkovsky drift can exceed that from Poynting-Robertson drag and radiation pressure by at least three orders of magnitude, possibly flinging asteroids which survive YORP spin-up into a widely dispersed cloud around the resulting white dwarf. The GB Yarkovsky effect alone may change an asteroid's orbital eccentricity by ten per cent in just one Myr. Damping perturbations from stellar wind drag can be just as extreme, but are strongly dependent on the highly uncertain local gas density and mean free path length. We conclude that GB radiative and wind effects must be considered when modelling the post-main-sequence evolution of bodies smaller than about 1000 km.
We present optical photometry and spectroscopy of the optical transient SN 2011A. Our data spans 140 days after discovery including $BVRIu'g'r'i'z'$ photometry and 11 epochs of optical spectroscopy. Originally classified as a type IIn supernova (SN IIn) due to the presence of narrow H$\alpha$ emission, this object shows exceptional characteristics. Firstly, the light curve shows a double plateau; a property only observed before in the impostor SN 1997bs. Secondly, SN 2011A has a very low luminosity ($M_{V}=-15.72$), placing it between normal luminous SNe IIn and SN impostors. Thirdly, SN 2011A shows low velocity and high equivalent width absorption close to the sodium doublet, which increases with time and is most likely of circumstellar origin. This evolution is also accompanied by a change of line profile; when the absorption becomes stronger, a P-Cygni profile appears. We discuss SN 2011A in the context of interacting SNe IIn and SN impostors, which appears to confirm the uniqueness of this transient. While we favour an impostor origin for SN 2011A, we highlight the difficulty in differentiating between terminal and non-terminal interacting transients.
We calibrate here cosmological radiative transfer simulation with ATON/RAMSES with a range of measurements of the Lyman alpha opacity from QSO absorption spectra. We find the Lyman alpha opacity to be very sensitive to the exact timing of hydrogen reionisation. Models reproducing the measured evolution of the mean photoionisation rate and average mean free path reach overlap at z ~ 7 and predict an accelerated evolution of the Lyman alpha opacity at z > 6 consistent with the rapidly evolving luminosity function of Lyman alpha emitters in this redshift range. Similar to "optically thin" simulations our full radiative transfer simulations fail, however, to reproduce the high-opacity tail of the Lyman alpha opacity PDF at z > 5. We argue that this is due to spatial UV fluctuations in the post-overlap phase of reionisation on substantially larger scales than predicted by our source model, where the ionising emissivity is dominated by large numbers of sub-L* galaxies. We further argue that this suggests a significant contribution to the ionising UV background by much rarer bright sources at high redshift.
Extremely metal-poor (EMP) stars ([Fe/H] < -3.0 dex) provide a unique window into understanding the first generation of stars and early chemical enrichment of the Universe. EMP stars are exceptionally rare, however, and the relatively small number of confirmed discoveries limits our ability to exploit these near-field probes of the first ~500 Myr after the Big Bang. Here, a new method to photometrically estimate [Fe/H] from only broadband photometric colors is presented. I show that the method, which utilizes machine-learning algorithms and a training set of ~170,000 stars with spectroscopically measured [Fe/H], produces a typical scatter of ~0.29 dex. This performance is similar to what is achievable via low-resolution spectroscopy, and outperforms other photometric techniques, while also being more general. I further show that a slight alteration to the model, wherein synthetic EMP stars are added to the training set, yields the robust identification of EMP candidates. In particular, this synthetic-oversampling method recovers ~20% of the EMP stars in the training set, at a precision of ~0.05. Furthermore, ~65% of the false positives from the model are very metal-poor stars ([Fe/H] < -2.0 dex). The synthetic-oversampling method is biased towards the discovery of warm (~F-type) stars, a consequence of the targeting bias from the SDSS/SEGUE survey. This EMP selection method represents a significant improvement over alternative broadband optical selection techniques. The models are applied to >12 million stars, with an expected yield of ~600 new EMP stars, which promises to open new avenues for exploring the early universe.
We present deep 10h VLT/XSHOOTER spectroscopy for an extraordinarily luminous and extended Lya emitter at z=6.595 referred to as Himiko and first discussed by Ouchi et al. (2009), with the purpose of constraining the mechanisms powering its strong emission. Complementary to the spectrum, we discuss NIR imaging data from the CANDELS survey. We find neither for HeII nor any metal line a significant excess, with 3 sigma upper limits of 6.8, 3.1, and 5.8x10^{-18} erg/s/cm^2 for CIV $\lambda$1549, HeII $\lambda$1640, CIII] $\lambda$1909, respectively, assuming apertures with 200 km/s widths and offset by -250 km/s w.r.t to the peak Lya redshift. These limits provide strong evidence that an AGN is not a major contribution to Himiko's Lya flux. Strong conclusions about the presence of PopIII star-formation or gravitational cooling radiation are not possible based on the obtained HeII upper limit. Our Lya spectrum confirms both spatial extent and flux (8.8+/-0.5x10^{-17} erg/s/cm^2) of previous measurements. In addition, we can unambiguously exclude any remaining chance of it being a lower redshift interloper by significantly detecting a continuum redwards of Lya, while being undetected bluewards.
In this work we are reporting on the measurement of the proton-air inelastic cross section $\sigma^{\rm inel}_{\rm p-air}$ using the Telescope Array (TA) detector. Based on the measurement of the $\sigma^{\rm inel}_{\rm p-air}$ the proton-proton cross section $\sigma_{\rm p-p}$ value is also determined at $\sqrt{s} = 95$ TeV. Detecting cosmic ray events at ultra high energies with Telescope Array enables us to study this fundamental parameter that we are otherwise unable to access with particle accelerators. The data used in this report is collected over five years using hybrid events observed by the Middle Drum fluorescence detector together with the surface array detector. The value of the $\sigma^{\rm inel}_{\rm p-air}$ is found to be equal to $ 567.0 \pm 70.5 [{\rm Stat.}] ^{+25}_{-29} [{\rm Sys.}]$ mb. The total proton-proton cross section is subsequently inferred from Glauber Formalism and Block, Halzen and Stanev QCD inspired fit and is found to be equal to $170_{-44}^{+48} [{\rm Stat.}] \pm _{-19}^{+17} [{\rm Sys.}] $mb.
When neutron stars reside in transient X-ray binaries, their crustal layers become heated during accretion outbursts and subsequently cool in quiescence. Observing and modeling this thermal response has yielded valuable insight into the physics of neutron star crusts. However, one unresolved problem is the evidence in several neutron stars for an extra energy source, located at shallow depth in the crust, that is not accounted for by standard heating models. Its origin remains puzzling, and it is currently unclear whether this additional heating occurs in all neutron stars, and if the magnitude is always the same. Here, we report on Chandra observations that cover two years after the 2012 outburst of the transient neutron star X-ray binary Swift J174805.3-244637 in the globular cluster Terzan 5. The temperature of the neutron star was elevated during the first two months following its ~8 week accretion episode, but had decayed to the pre-outburst level within ~100 days. Interpreting this as rapid cooling of the accretion-heated crust, we model the observed temperature curve with a thermal evolution code. We find that there is no need to invoke shallow heating for this neutron star, although an extra energy release up to ~1.4 MeV/nucleon is allowed by the current data (2-sigma confidence). We also present two new data points on the crust cooling curve of the 11-Hz X-ray pulsar IGR J17480-2446 in Terzan 5, which was active in 2010. The temperature of this neutron star remained significantly above its pre-outburst level, but we detect no temperature change since the previous measurements of 2013 February. This is consistent with the slower cooling expected several years post-outburst.
We consider evolution properties of galaxies and quasars with steep radio spectrum at the decametre band from the UTR-2 catalogue. The ratios of source's monochromatic luminosities at the decametre and high-frequency bands display the dependence on the redshift, linear size, characteristic age of examined objects. At that, the mean values of corresponding ratios for considered galaxies and quasars have enough close quantities, testifying on the unified model of sources. We analyse obtained relations for two types of steep-spectrum sources (with linear steep spectrum (S) and low-frequency steepness after a break (C+)) from the UTR-2 catalogue.
We present a mass map reconstructed from weak gravitational lensing shear measurements over 139 sq. deg from the Dark Energy Survey (DES) Science Verification data. The mass map probes both luminous and dark matter, thus providing a tool for studying cosmology. We find good agreement between the mass map and the distribution of massive galaxy clusters identified using a red-sequence cluster finder. Potential candidates for super-clusters and voids are identified using these maps. We measure the cross-correlation between the mass map and a magnitude-limited foreground galaxy sample and find a detection at the 5-7 sigma level on a large range of scales. These measurements are consistent with simulated galaxy catalogs based on LCDM N-body simulations, suggesting low systematics uncertainties in the map. We summarize our key findings in this letter; the detailed methodology and tests for systematics are presented in a companion paper.
The spiral structure of the Milky Way is highly uncertain and is the subject of much discussion nowadays. We present the first result from a program that determines the properties of the Local spiral arm (LOA), together with a full description of the program. In this context we have made a comprehensive study of the young LOA open cluster NGC 2302, which includes a UBVRI photometric analysis and determination of its kinematic properties - proper motion and radial velocity - and of its orbital parameters. We determined the mean PM of NGC 2302 relative to the local field of disk stars, and, through a comparison with the UCAC4 catalog, we transformed this relative PM into an absolute one. Using medium-resolution spectroscopy of 26 stars in the field of NGC 2302, we derived its mean RV. Isochrone fits to the photometric diagrams allowed us to determine the fundamental parameters of NGC 2302, including reddening, distance, and age. The kinematic data and derived distance allowed us to determine the space motion of NGC 2302. This was done by adopting a time-independent, axisymmetric, and fully analytic gravitational potential for the MW. We obtained an absolute PM for NGC 2302 of ($\mu_{\alpha} \cos\delta,\mu_{\delta}) = (-2.09,-2.11)$ mas/yr, with standard errors of 0.410 and 0.400 mas/yr. The mean RV of NGC 2302 turned out to be 31.2 km/sec with a standard error of 0.7 km/sec. Isochrone fits displaced for this reddening and for a distance modulus of (m-M)o = 10.69 indicate an age of log(t) = 7.90-8.00 with a slight tendency toward the younger age. Inspection of the shape of the orbit of NGC 2302 and the resulting orbital parameters indicate that it is a typical population I object.
Distance measurements are currently the most powerful tool to study the expansion history of the universe without specifying its matter content nor any theory of gravitation. Assuming only an isotropic, homogeneous and flat universe, in this work we introduce a model-independent method to reconstruct directly the deceleration function via a piecewise function. Including a penalty factor, we are able to vary continuously the complexity of the deceleration function from a linear case to an arbitrary $(n+1)$-knots spline interpolation. We carry out a Monte Carlo analysis to determine the best penalty factor, evaluating the bias-variance trade-off, given the uncertainties of the SDSS-II and SNLS supernova combined sample (JLA), compilations of baryon acoustic oscillation (BAO) and $H(z)$ data. We show that, evaluating a single fiducial model, the conclusions about the bias-variance ratio are misleading. We determine the reconstruction method in which the bias represents at most $10\%$ of the total uncertainty. In all statistical analyses, we fit the coefficients of the deceleration function along with four nuisance parameters of the supernova astrophysical model. For the full sample, we also fit $H_0$ and the sound horizon $r_s(z_d)$ at the drag redshift. The bias-variance trade-off analysis shows that, apart from the deceleration function, all other estimators are unbiased. Finally, we apply the Ensemble Sampler Markov Chain Monte Carlo (ESMCMC) method to explore the posterior of the deceleration function up to redshift $1.3$ (using only JLA) and $2.3$ (JLA+BAO+$H(z)$). We obtain that the standard cosmological model agrees within $3\sigma$ level with the reconstructed results in the whole studied redshift intervals. Since our method is calibrated to minimize the bias, the error bars of the reconstructed functions are a good approximation for the total uncertainty.
The Vertical Shear Instability is one of two known mechanisms potentially active in the so-called dead zones of protoplanetary accretion disks. A recent analysis indicates that a subset of unstable modes shows unbounded growth - both as resolution is increased and when the nominal lid of the atmosphere is extended, possibly indicating ill-posedness in previous attempts of linear analysis. The reduced equations governing the instability are revisited and the generated solutions are examined using both the previously assumed separable forms and an improved non-separable solution form that is herewith introduced. Analyzing the reduced equations using the separable form shows that, while the low-order body modes have converged eigenvalues and eigenfunctions (as both the vertical boundaries of the atmosphere are extended and with increased radial resolution), it is also confirmed that the corresponding high-order body modes and the surface modes do indeed show unbounded growth rates. However, the energy contained in both the higher-order body modes and surface modes diminishes precipitously due to the disk's Gaussian density profile. Most of the energy of the instability is contained in the low-order modes. An inseparable solution form is introduced which filters out the inconsequential surface modes leaving only body modes (both low and high-order ones). The analysis predicts a fastest growing mode with a specific radial length scale. The growth rates associated with the fundamental corrugation and breathing modes matches the growth and length scales observed in previous nonlinear studies of the instability.
We studied the distribution of projected ellipticity $\epsilon$ for the galaxies in a sample of 20 rich ($\mathcal{R} > 2$) nearby ($z < 0.1$) clusters of galaxies. We find no evidence of differences in the distribution of $\epsilon$ between different clusters, except possibly for the Coma Cluster (the probability that Coma has the same distribution as the other clusters is $P < 10^{-3}$). We then study the distribution of $\epsilon$ within the clusters, and find that $\epsilon$ increases with the projected cluster-centric radius R. The trend is preserved at fixed magnitude, showing that our result exists over and above the trend of more luminous galaxies to be both rounder and more common in the centre of clusters. The trend of $\epsilon$ with R persists even for the subsample of intrinsically flat galaxies ($\epsilon$ > 0.4), therefore it is not a consequence of the increasing fraction of intrinsically round Slow Rotator galaxies nearer the centre of clusters. The trend is also observed for flat, smooth galaxies and dividing the galaxies according to the shape of their light profile. This suggests that the systematic variation of the fraction of spiral galaxies with R does not explain alone the observed trend. We interpret our findings in light of the classification of Early Type Galaxies (ETGs) as Fast and Slow Rotators. We conclude that the observed trend of decreasing $\epsilon$ nearer to the centre of clusters is evidence for physical effects in clusters causing Fast Rotator ETGs to have a lower average intrinsic ellipticity near the centre of rich clusters.
During the "WISE at 5: Legacy and Prospects" conference in Pasadena, CA -- which ran from February 10 - 12, 2015 -- attendees were invited to engage in an interactive session exploring the future uses of the Wide-field Infrared Survey Explorer (WISE) data. The 65 participants -- many of whom are extensive users of the data -- brainstormed the top questions still to be answered by the mission, as well as the complementary current and future datasets and additional processing of WISE/NEOWISE data that would aid in addressing these most important scientific questions. The results were mainly bifurcated between topics related to extragalactic studies (e.g. AGN, QSOs) and substellar mass objects. In summary, participants found that complementing WISE/NEOWISE data with cross-correlated multiwavelength surveys (e.g. SDSS, Pan-STARRS, LSST, Gaia, Euclid, etc.) would be highly beneficial for all future mission goals. Moreover, developing or implementing machine-learning tools to comb through and understand cross-correlated data was often mentioned for future uses. Finally, attendees agreed that additional processing of the data such as co-adding WISE and NEOWISE and extracting a multi-epoch photometric database and parallax and proper motion catalog would greatly improve the scientific results of the most important projects identified. In that respect, a project such as MaxWISE which would execute the most important additional processing and extraction as well as make the data and catalogs easily accessible via a public portal was deemed extremely important.
Guided by non-relativistic Effective Field Theory (EFT) we classify the most general spin-dependent interactions between a fermionic Weakly Interacting Massive Particle (WIMP) and nuclei, and within this class of models we discuss the viability of an interpretation of the DAMA modulation result in terms of a signal from WIMP elastic scatterings using a halo-independent approach. We find that, although several relativistic EFT's can lead to a spin-dependent cross section, in some cases with an explicit, non-negligible dependence on the WIMP incoming velocity, three main scenarios can be singled out in the non-relativistic limit which approximately encompass them all, and that only differ by their dependence on the transferred momentum. For two of them compatibility between DAMA and other constraints is possible for a WIMP mass below 30 GeV, but only for a WIMP velocity distribution in the halo of our Galaxy which departs from a Maxwellian. This is achieved by combining a suppression of the WIMP effective coupling to neutrons (to evade constraints from xenon and germanium detectors) to an explicit quadratic or quartic dependence of the cross section on the transferred momentum (that leads to a relative enhancement of the expected rate off sodium in DAMA compared to that off fluorine in droplet detectors and bubble chambers). For larger WIMP masses the same scenarios are excluded by scatterings off iodine in COUPP.
Long expected transition states between the rotation powered and accretion powered non-thermal emission in the millisecond pulsar binary systems have been recently observed in the case of three objects PSR J1023+0038, PSR J1824-2452, and PSR J1227-4859. Surprisingly, the transition is related to the significant change in the $\gamma$-ray flux being a factor of a few higher with the presence of an accretion disk. The origin of this enhanced emission seems to be related to the penetration of the inner pulsar magnetosphere by the accretion disk. We propose that the radiation processes, characteristic for the rotation powered pulsar, can co-exist with the presence of an accretion disk in the inner pulsar magnetosphere. In our scenario additional $\gamma$-ray emission is produced by secondary leptons, originated close to the acceleration gap, which Compton up-scatter thermal radiation from the accretion disk to GeV energies. The accretion disk penetrates deep into the pulsar magnetosphere allowing the matter to fall onto the NS surface producing pulsed X-ray emission. We show that the sum of the rotation powered pulsar $\gamma$-ray emission, produced by the primary electrons in the curvature process, and the $\gamma$-ray emission, produced by secondary leptons, can explain the observed high energy radiation from the redback binary pulsar PSR J1227-4853 in the state with evidences of the accretion disk.
ALMA observations of a small sample of transitional disks with large dust cavities observed in Cycle 0 and 1 are summarized. The gas and dust surface density structures are inferred from the continuum and 12CO, 13CO and C18O line data using the DALI physical-chemical code. Thanks to its ability to self-shield, CO can survive inside dust cavities in spite of being exposed to intense UV radiation and can thus be used as a probe of the gas structure. Modeling of the existing data shows that gas is present inside the dust cavities in all cases, but at a reduced level compared with the gas surface density profile of the outer disk. The gas density decrease inside the dust cavity radius by factors of up to 10^4 suggests clearing by one or more planetary-mass companions. The accompanying pressure bumps naturally lead to trapping of the mm-sized dust grains observed in the ALMA images.
Using 2D numerical hydrodynamical simulations of type Ia supernova remnants (SNR Ia) we show that iron clumps few times denser than the rest of the SN ejecta might form protrusions in an otherwise spherical SNR. Such protrusions exist in some SNR Ia, e.g., SNR 1885 and Tycho. Iron clumps are expected to form in the deflagration to detonation explosion model. In SNR Ia where there are two opposite protrusions, termed ears, such as Kepler's SNR and SNR G1.9+0.3, our scenario implies that the dense clumps, or iron bullets, were formed along an axis. Such a preferred axis can result from a rotating white dwarf progenitor. If our claim holds, this offers an important clue to the SN Ia explosion scenario.
The variation of solar irradiance is one of the natural forcing mechanisms of the terrestrial climate. Hence, the time-dependent solar irradiance is an important input parameter for climate modelling. The solar surface magnetic field is a powerful proxy for solar irradiance reconstruction. The analyses of data obtained with the Michelson Doppler Imager (MDI) on board the SOHO mission are therefore useful for the identification of solar surface magnetic features to be used in solar irradiance reconstruction models. However, there is still a need for automated technologies that would enable the identification of solar activity features from large databases. To achieve this we present a series of enhanced segmentation algorithms developed to detect and calculate the area coverages of specific magnetic features from MDI intensitygrams and magnetograms. These algorithms are part of the Automated Solar Activity Prediction (ASAP) tool. The segmentation algorithms allow us to identify the areas on the solar disk covered by magnetic elements inside and outside boundaries of active regions. Depending on their contrast properties, magnetic features within an active region boundary are classified as sunspot umbra and penumbra, or faculae. Outside an active region boundary magnetic elements are identified as network. We present the detailed steps involved in the segmentation process and provide the area coverages of the segmented MDI intensitygrams and magnetograms. The feature segmentation has been carried out on daily intensitygrams and magnetograms from April 21, 1996 to April 11, 2011. This offers an exciting opportunity to undertake further investigations that benefit from solar features segmentations, such as solar irradiance reconstruction, which we plan to investigate in the future.
We estimate the rate of dark matter scattering in collapsed structures throughout the history of the Universe. If the scattering cross-section is velocity-independent, then the canonical picture is correct that scatterings occur mainly at late times. The scattering rate peaks slightly at redshift z~6, and remains significant today. Half the scatterings occur after z~1, in structures more massive than 10^12 M_sun. Within a factor of two, these numbers are robust to changes in the assumed astrophysics, and the scatterings would be captured in cosmological simulations. However, for particle physics models with a velocity-dependent cross-section (as for Yukawa potential interactions via a massive mediator), the scattering rate peaks before z~20, in objects with mass less than 10^4 M_sun. These precise values are sensitive to the redshift-dependent mass-concentration relation and the small-scale cutoff in the matter power spectrum. In extreme cases, the qualitative effect of early interactions may be reminiscent of warm dark matter and strongly affect the subsequent growth of structure. However, these scatterings are being missed in existing cosmological simulations with limited mass resolution.
Time is a parameter playing a central role in our most fundamental modeling of natural laws. Relativity theory shows that the comparison of times measured by different clocks depends on their relative motions and on the strength of the gravitational field in which they are embedded. In standard cosmology, the time parameter is the one measured by fundamental clocks, i.e. clocks at rest with respect to the expanding space. This proper time is assumed to flow at a constant rate throughout the whole history of the Universe. We make the alternative hypothesis that the rate at which cosmological time flows depends on the dynamical state of the Universe. In thermodynamics, the arrow of time is strongly related to the second law, which states that the entropy of an isolated system will always increase with time or, at best, stay constant. Hence, we assume that time measured by fundamental clocks is proportional to the entropy of the region of the Universe that is causally connected to them. Under that simple assumption, we build a cosmological model that explains the Type Ia Supernovae data (the best cosmological standard candles) without the need for exotic dark matter nor dark energy.
Chemical tagging of stellar debris from disrupted open clusters and associations underpins the science cases for next-generation multi-object spectroscopic surveys. As part of the Galactic Archaeology project TraCD (Tracking Cluster Debris), a preliminary attempt at reconstructing the birth clouds of now phase-mixed thin disk debris is undertaken using a parametric minimum spanning tree (MST) approach. Empirically-motivated chemical abundance pattern uncertainties (for a 10-dimensional chemistry-space) are applied to NBODY6-realised stellar associations dissolved into a background sea of field stars, all evolving in a Milky Way potential. We demonstrate that significant population reconstruction degeneracies appear when the abundance uncertainties approach 0.1 dex and the parameterised MST approach is employed; more sophisticated methodologies will be required to ameliorate these degeneracies.
The N-independence observed in the evolution of massive black hole binaries (MBHBs) in recent simulation of merging stellar bulges suggests a simple interpretation beyond complex time-dependent relaxation processes. We conjecture that the MBHB hardening rate is equivalent to that of a binary immersed in a field of unbound stars with density $\rho$ and typical velocity $\sigma$, provided that $\rho$ and $\sigma$ are the stellar density and the velocity dispersion at the influence radius of the MBHB. By comparing direct N-body simulations to an hybrid model based on 3-body scattering experiments, we verify this hypothesis: when normalized to the stellar density and velocity dispersion at the binary influence radius, the N-body MBHB hardening rate approximately matches that predicted by 3-body scatterings in the investigated cases. The eccentricity evolution obtained with the two techniques is also in reasonable agreement. This result is particularly practical because it allows to estimate the lifetime of MBHBs forming in dry mergers based solely on the stellar density profile of the host galaxy. We briefly discuss some implications of our finding for the gravitational wave signal observable by pulsar timing arrays and for the expected population of MBHBs lurking in massive ellipticals.
Lincoln Near-Earth Asteroid Research asteroid survey (LINEAR) observed approximately 10,000 deg$^2$ of the northern sky in the period roughly from 1998 to 2013. Long baseline of observations combined with good cadence and depth ($14.5 < r_{SDSS}< 17.5$) provides excellent basis for investigation of variable and transient objects in this relatively faint and underexplored part of the sky. Details covering the repurposing of this survey for use in time domain astronomy, creation of a highly reliable catalogue of approximately 7,200 periodically variable stars (RR Lyrae, eclipsing binaries, SX Phe stars and LPVs) as well as search for optical signatures of exotic transient events (such as tidal disruption event candidates), are presented.
We present a statistical study of prominence and filament eruptions observed by the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO). Several properties are recorded for 904 events that were culled from the Heliophysics Event Knowledgebase (HEK) and incorporated into an online catalog for general use. These characteristics include the filament and eruption type, eruption symmetry and direction, apparent twisting and writhing motions, and the presence of vertical threads and coronal cavities. Associated flares and white-light coronal mass ejections (CME) are also recorded. Total rates are given for each property along with how they differ among filament types. We also examine the kinematics of 106 limb events to characterize the distinct slow- and fast-rise phases often exhibited by filament eruptions. The average fast-rise onset height, slow-rise duration, slow-rise velocity, maximum field-of-view (FOV) velocity, and maximum FOV acceleration are 83 Mm, 4.4 hours, 2.1 km/s, 106 km/s, and 111 m/s^2, respectively. All parameters exhibit lognormal probability distributions similar to that of CME speeds. A positive correlation between latitude and fast-rise onset height is found, which we attribute to a corresponding negative correlation in the average vertical magnetic field gradient, or decay index, estimated from potential field source surface (PFSS) extrapolations. We also find the decay index at the fast-rise onset point to be 1.1 on average, consistent with the critical instability threshold theorized for straight current channels. Finally, we explore relationships between the derived kinematics properties and apparent twisting motions. We find that events with evident twist have significantly faster CME speeds and significantly lower fast-rise onset heights, suggesting relationships between these values and flux rope helicity.
We report on NuSTAR, XMM-Newton and Swift observations of the gamma-ray binary 1FGL J1018.6-5856. We measure the orbital period to be 16.544+/-0.008 days using Swift data spanning 1900 days. The orbital period is different from the 2011 gamma-ray measurement which was used in the previous X-ray study of An et al. (2013) using ~400 days of Swift data, but is consistent with a new gamma-ray solution reported in 2014. The light curve folded on the new period is qualitatively similar to that reported previously, having a spike at phase 0 and broad sinusoidal modulation. The X-ray flux enhancement at phase 0 occurs more regularly in time than was previously suggested. A spiky structure at this phase seems to be a persistent feature, although there is some variability. Furthermore, we find that the source flux clearly correlates with the spectral hardness throughout all orbital phases, and that the broadband X-ray spectra measured with NuSTAR, XMM-Newton, and Swift are well fit with an unbroken power-law model. This spectrum suggests that the system may not be accretion-powered.
A very large lobe overflow event is suggested to explain the $0.^m4$ brightening observed in K band at pericenter passage of the star known as S2 that orbits the Galaxy's supermassive black hole (SMBH). Known observed properties of S2 that contribute to lobe filling are 1) the enormous mass ratio, $M_{SMBH}/M_{S2}$, 2) S2's fast rotation, and 3) S2's large orbital eccentricity. Published estimates have given limiting lobe sizes of order 100 to 300 $R_\odot$ but, with S2's fast rotation taken into account, the computed lobe size is much smaller, being compatible with either a main sequence OB star or a stripped evolved star. An important evolutionary consideration that predicts very large pericenter overflows is envelope expansion following mass loss that is characteristic of highly evolved stars. Material removed by lobe overflow at pericenter is replenished by envelope expansion as an evolved star awaits its next pericenter passage. An observational signature of lobe overflow for upcoming pericenter passages would be appearance of emission lines as the ejected gas expands and becomes optically thin.
We compare predictions of large-scale cosmological hydrodynamical simulations for neutral hydrogen absorption signatures in the vicinity of 1e11 - 1e12.5 MSun haloes with observational measurements. Two different hydrodynamical techniques and a variety of prescriptions for gas removal in high density regions are examined. Star formation and wind feedback play only secondary roles in the HI absorption signatures outside the virial radius, but play important roles within. Accordingly, we identify three distinct gaseous regions around a halo: the virialized region, the mesogalactic medium outside the virial radius arising from the extended haloes of galaxies out to about two turnaround radii, and the intergalactic medium beyond. Predictions for the amount of absorption from the mesogalactic and intergalactic media are robust across different methodologies, and the predictions agree with the amount of absorption observed around star-forming galaxies and QSO host galaxies. Recovering the measured amount of absorption within the virialized region, however, requires either a higher dynamic range in the simulations, additional physics, or both.
Microvariability consists in small time scale variations of low amplitude in the photometric light curves of quasars, and represents an important tool to investigate their inner core. Detection of quasar microvariations is challenging for their non-periodicity, as well as the need for high monitoring frequency and high signal-to-noise ratio. Statistical tests developed for the analysis of quasar differential light curves usually show either low power or low reliability, or both. In this paper we compare two statistical procedures that include several stars to perform tests with enhanced power and high reliability. We perform light curve simulations of variable quasars and non-variable stars, and analyze them with statistical procedures developed from the F-test and the analysis of variance. The results show a large improvement in the power of both statistical probes, and a larger reliability, when several stars are included in the analysis. The results from the simulations agree with those obtained from observations of real quasars. The high power and high reliability of the tests discussed in this paper improve the results that can be obtained from short and long time scale variability studies. These techniques are not limited to quasar variability; on the contrary, they can be easily implemented to other sources such as variable stars. Their applications to future research and to the analysis of large field photometric monitoring archives can reveal new variable sources.
We have carried out a study of radio emission from a small sample of magnetic O- and B-type stars using the Giant Metrewave Radio Telescope, with the goal of investigating their magnetospheres at low frequencies. These are the lowest frequency radio measurements ever obtained of hot magnetic stars. The observations were taken at random rotational phases in the 1390 and the 610 MHz bands. Out of the 8 stars, we detect five B-type stars in both the 1390 and the 610 MHz bands. The O-type stars were observed only in the 1390 MHz band, and no detections were obtained. We explain this result as a consequence of free-free absorption by the free-flowing stellar wind exterior to the closed magnetosphere. We also study the variability of individual stars. One star - HD 133880 - exhibits remarkably strong and rapid variability of its low frequency flux density. We discuss the possibility of this emission being coherent emission as reported for CU Vir by Trigilio et al. (2000).
The sharp change in slope of the ultra-high energy cosmic ray (UHECR) spectrum around 10^18.6 eV (the ankle), combined with evidence of a light but extragalactic component near and below the ankle and intermediate composition above, has proved exceedingly challenging to understand theoretically. We propose a mechanism whereby photo-disintegration of ultra-high energy nuclei in the region surrounding a UHECR accelerator {\it naturally} accounts for the observed spectrum and inferred composition at Earth. Under a range of reasonable source conditions, the model reproduces the spectrum and the composition over the entire extragalactic CR range, i.e. above 10^18.6 eV. Predictions for the spectrum and flavors of neutrinos resulting from this process are also presented.
Evidence for dark matter self-interactions has recently been reported based on the observation of a spatial offset between the dark matter halo and the stars in a galaxy in the cluster Abell 3827. Interpreting the offset as due to dark matter self-interactions leads to a cross section measurement of sigma_DM/m ~ (1-1.5) cm^2/g, where m is the mass of the dark matter particle. We use this observation to constrain singlet scalar dark matter coupled to the Standard Model and to two-Higgs-doublet models. We show that the most natural scenario in this class of models is very light dark matter, below about 0.1 GeV, whose relic abundance is set by freeze-in, i.e., by slow production of dark matter in the early universe via extremely tiny interactions with the Higgs boson, never reaching thermal equilibrium. We also show that the dark matter abundance can be established through the usual thermal freeze-out mechanism in the singlet scalar extension of the Yukawa-aligned two-Higgs-doublet model, but that it requires rather severe fine tuning of the singlet scalar mass.
We study the effects on the spectrum and distribution of high-energy neutrinos due to scattering with dark matter both outside and within our galaxy, focusing on the neutrinos observed by the IceCube experiment with energies up to several PeV. If these neutrinos originate from extra-galactic astrophysical sources, then scattering in transit with dark matter particles will delay their arrival to Earth. This results in a cut-off in their spectrum at an energy set by the scattering cross section, allowing us to place an upper limit on cross sections $\sigma$ which increase with energy E at the level of $\sigma$ < 10^{-17} x (m / GeV) x (E / PeV)^2 cm^2, for dark matter particles of mass m. Once these neutrinos enter our galaxy, the large dark matter densities result in further scattering, especially towards the Galactic Centre. Intriguingly, we find that for $\sigma$ ~ 10^{-22} x (m / GeV) x (E / PeV)^2 cm^2, the distribution of the neutrinos on the sky has a small cluster of events towards the centre of the galaxy, potentially explaining the ~2 sigma excess seen by IceCube in this region without needing a galactic source.
We study collisions between nearly planar domain walls including the effects of small initial nonplanar fluctuations. These perturbations represent the small fluctuations that must exist in a quantum treatment of the problem. In a previous paper, we demonstrated that at the linear level a subset of these fluctuations experience parametric amplification as a result of their coupling to the planar symmetric background. Here we study the full three-dimensional nonlinear dynamics using lattice simulations, including both the early time regime when the fluctuations are well described by linear perturbation theory as well as the subsequent stage of fully nonlinear evolution. We find that the nonplanar fluctuations have a dramatic effect on the overall evolution of the system. Specifically, once these fluctuations begin to interact nonlinearly the split into a planar symmetric part of the field and the nonplanar fluctuations loses its utility. At this point the colliding domain walls dissolve, with the endpoint of this being the creation of a population of oscillons in the collision region. The original (nearly) planar symmetry has been completely destroyed at this point and an accurate study of the system requires the full three-dimensional simulation.
We study a universal structure of two and three dimensional self-gravitating systems in the quasi-equilibrium state. It is shown numerically that the two dimensional self-gravitating system in the quasi-equilibrium state has the same kind of density profile as the three dimensional one. We develop a phenomenological model to describe this universal structure by using a special Langevin equation with a distinctive random noise to self-gravitating systems. We find that the density profile derived theoretically is consistent well with results of observations and simulations.
In the Bona-Masso formulation, Einstein equations are written as a set of flux conservative first order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for Numerical Relativity. Our model is validated with well-established tests, showing good agreement with analytical solutions. Furthermore, we show that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improves. Finally, in order to show the potential of our approach a linear scaling law for parallelisation with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.
A GPS carrier-phase frequency transfer link along a baseline of 450 km has been established and is characterized by comparing it to a phase-stabilized optical fiber link of 920 km length, established between the two endpoints, the Max-Planck-Institut f\"ur Quantenoptik in Garching and the Physikalisch-Technische Bundesanstalt in Braunschweig. The characterization is accomplished by comparing two active hydrogen masers operated at both institutes. The masers serve as local oscillators and cancel out when the double differences are calculated, such that they do not constitute a limitation for the GPS link characterization. We achieve a frequency instability of 3 x 10^(-13) in 30 s and 5 x 10^(-16) for long averaging times. Frequency comparison results obtained via both links show no deviation larger than the statistical uncertainty of 6 x 10^(-16). These results can be interpreted as a successful cross-check of the measurement uncertainty of a truly remote end fiber link.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
We employ the all-sky map of the anomalous microwave emission (AME) produced by component separation of the microwave sky to study correlations between the AME and Galactic dust properties. We find that while the AME is highly correlated with all tracers of dust emission, fluctuations in the AME intensity per dust optical depth are uncorrelated with fluctuations in the emission from polycyclic aromatic hydrocarbons (PAHs), casting doubt on the association between AME and PAHs. Further, we find that the best predictor of the AME strength is the dust radiance and that the AME intensity increases with increasing radiation field strength, at variance with predictions from the spinning dust hypothesis. A reconsideration of other emission mechanisms, such as magnetic dipole emission, is warranted.
Giant protoplanets evacuate a gap in their host protoplanetary disc, which gas must cross before it can be accreted. A magnetic field is likely carried into the gap, potentially influencing the flow. Gap crossing has been simulated with varying degrees of attention to field evolution (pure hydrodynamical, ideal, and resistive MHD), but as yet there has been no detailed assessment of the role of the field accounting for all three key non-ideal MHD effects: Ohmic resistivity, ambipolar diffusion, and Hall drift. We present a detailed investigation of gap magnetic field structure as determined by non-ideal effects. We assess susceptibility to turbulence induced by the magnetorotational instability, and angular momentum loss from large-scale fields. As full non-ideal simulations are computationally expensive, we take an a posteriori approach, estimating MHD quantities from the pure hydrodynamical gap crossing simulation by Tanigawa et al. (2012). We calculate the ionisation fraction and estimate field strength and geometry to determine the strength of non-ideal effects. We find that the protoplanetary disc field would be easily drawn into the gap and circumplanetary disc. Hall drift dominates, so that much of the gap is conditionally MRI unstable depending on the alignment of the field and disc rotation axes. Field alignment also influences the strong toroidal field component permeating the gap. Large-scale magnetic forces are small in the circumplanetary disc, indicating they cannot drive accretion there. However, turbulence will be key during satellite growth as it affects critical disc features, such as the location of the ice line.
ALMA Cycle 2 observations of the long wavelength dust emission in 180 star-forming (SF) galaxies are used to investigate the evolution of ISM masses at z = 1 to 6.4. The ISM masses exhibit strong increases from z = 0 to $\rm <z>$ = 1.15 and further to $\rm <z>$ = 2.2 and 4.8, particularly amongst galaxies above the SF galaxy main sequence (MS). The galaxies with highest SFRs at $\rm <z>$ = 2.2 and 4.8 have gas masses 100 times that of the Milky Way and gas mass fractions reaching 50 to 80\%, i.e. gas masses 1 - 4$\times$ their stellar masses. For the full sample of galaxies, we find a single, very simple SF law: $\rm SFR \propto M_{\rm ISM}^{0.9}$, i.e. a `linear' dependence on the ISM mass -- on and above the MS. Thus, the galaxies above the MS are converting their larger ISM masses into stars on a timescale similar to those on the MS. At z $> 1$, the entire population of star-forming galaxies has $\sim$5 - 10$\times$ shorter gas depletion times ($\sim0.2$ Gyr) than galaxies at low redshift. These {\bf shorter depletion times are due to a different, dominant mode of SF in the early universe} -- dynamically driven by compressive, high dispersion gas motions and/or galaxy interactions. The dispersive gas motions are a natural consequence of the extraordinarily high gas accretion rates which must occur to maintain the prodigious SF.
We present the results of a deep study of the isolated dwarf galaxies Andromeda XXVIII and Andromeda XXIX with Gemini/GMOS and Keck/DEIMOS. Both galaxies are shown to host old, metal-poor stellar populations with no detectable recent star formation, conclusively identifying both of them as dwarf spheroidal galaxies (dSphs). And XXVIII exhibits a complex horizontal branch morphology, which is suggestive of metallicity enrichment and thus an extended period of star formation in the past. Decomposing the horizontal branch into blue (metal poor, assumed to be older) and red (relatively more metal rich, assumed to be younger) populations shows that the metal rich are also more spatially concentrated in the center of the galaxy. We use spectroscopic measurements of the Calcium triplet, combined with the improved precision of the Gemini photometry, to measure the metallicity of the galaxies, confirming the metallicity spread and showing that they both lie on the luminosity-metallicity relation for dwarf satellites. Taken together, the galaxies exhibit largely typical properties for dSphs despite their significant distances from M31. These dwarfs thus place particularly significant constraints on models of dSph formation involving environmental processes such as tidal or ram pressure stripping. Such models must be able to completely transform the two galaxies into dSphs in no more than two pericentric passages around M31, while maintaining a significant stellar populations gradient. Reproducing these features is a prime requirement for models of dSph formation to demonstrate not just the plausibility of environmental transformation but the capability of accurately recreating real dSphs.
We study collisions between pairs of bubbles nucleated in an ambient false vacuum. For the first time, we include the effects of small initial (quantum) fluctuations around the instanton profiles describing the most likely initial bubble profile. Past studies of this problem neglect these fluctuations and work under the assumption that the collisions posess an exact SO(2,1) symmetry. We use three-dimensional lattice simulations to demonstrate that for double-well potentials, small initial perturbations to this symmetry can be amplified as the system evolves. Initially the amplification is well-described by linear perturbation theory around the SO(2,1) background, but the onset of strong nonlinearities amongst the fluctuations quickly leads to a drastic breaking of the original SO(2,1) symmetry and the production of oscillons in the collision region. We explore several single-field models, and we find it is hard to both realize inflation inside of a bubble and produce oscillons in a collision. Finally, we extend our results to a simple two-field model. The additional freedom allowed by the second field allows us to construct viable inflationary models that allow oscillon production in collisions. The breaking of the SO(2,1) symmetry allows for a new class of observational signatures from bubble collisions that do not posess azimuthal symmetry, including the production of gravitational waves which cannot be supported by an SO(2,1) spacetime.
It is difficult to understand how cold circumstellar disks accrete onto their central stars. A hydrodynamic mechanism, the vertical shear instability (VSI), offers a means to drive angular momentum transport in cold accretion disks such as protoplanetary disks (PPDs). The VSI is driven by a weak vertical gradient in the disk's orbital motion. In order to grow, the VSI must overcome vertical buoyancy, a strongly stabilizing influence in cold disks, where heating is dominated by external irradiation. Rapid cooling, via radiative losses, reduces the effective buoyancy and allows the VSI to operate. In this paper, we quantify the cooling timescale, $t_c$, needed for growth of the VSI. We perform a linear analysis of the VSI with cooling in vertically global and radially local disk models. For irradiated disks, we find that the VSI is most vigorous for rapid cooling with $t_c < \Omega_\mathrm{K}^{-1} h |q| / (\gamma -1)$ in terms of the Keplerian orbital frequency, $\Omega_\mathrm{K}$, the disk's aspect ratio, $ h \ll 1$, the radial power-law temperature gradient, $q$, and the adiabatic index, $\gamma$. For longer cooling times, the VSI is much less effective because growth slows and shifts to smaller length scales, which are more prone to viscous or turbulent decay. We apply our results to PPD models where $t_c$ is determined by the opacity of dust grains. We find that the VSI is most effective at intermediate radii, from $\sim 5$AU to $\sim 50$AU with a characteristic growth time of $\sim 30$ local orbital periods. Growth is suppressed by long cooling times both in the opaque inner disk and the optically thin outer disk. A reduction in the dust opacity by a factor of 10 increases cooling times enough to quench the VSI at all disk radii. Thus the formation of solid protoplanets, a sink for dust grains, can impede the VSI.
We study the characteristics of the galaxy cluster samples expected from the European Space Agency's Euclid satellite and forecast constraints on cosmological parameters describing a variety of cosmological models. The method used in this paper, based on the Fisher Matrix approach, is the same one used to provide the constraints presented in the Euclid Red Book (Laureijs et al.2011). We describe the analytical approach to compute the selection function of the photometric and spectroscopic cluster surveys. Based on the photometric selection function, we forecast the constraints on a number of cosmological parameter sets corresponding to different extensions of the standard LambdaCDM model. The dynamical evolution of dark energy will be constrained to Delta w_0=0.03 and Delta w_a=0.2 with free curvature Omega_k, resulting in a (w_0,w_a) Figure of Merit (FoM) of 291. Including the Planck CMB covariance matrix improves the constraints to Delta w_0=0.02, Delta w_a=0.07 and a FoM=802. The amplitude of primordial non-Gaussianity, parametrised by f_NL, will be constrained to \Delta f_NL ~ 6.6 for the local shape scenario, from Euclid clusters alone. Using only Euclid clusters, the growth factor parameter \gamma, which signals deviations from GR, will be constrained to Delta \gamma=0.02, and the neutrino density parameter to Delta Omega_\nu=0.0013 (or Delta \sum m_\nu=0.01). We emphasise that knowledge of the observable--mass scaling relation will be crucial to constrain cosmological parameters from a cluster catalogue. The Euclid mission will have a clear advantage in this respect, thanks to its imaging and spectroscopic capabilities that will enable internal mass calibration from weak lensing and the dynamics of cluster galaxies. This information will be further complemented by wide-area multi-wavelength external cluster surveys that will already be available when Euclid flies. [Abridged]
By means of fully kinetic simulations, we show that magnetic reconnection is highly efficient at converting magnetic energy into nonthermal electrons in the nonrelativistic low-beta regimes. This leads to a nonthermally dominated distribution with a power-law tail that contains more than half of the electrons in the system. The main acceleration mechanism is a Fermi-type acceleration accomplished by the curvature drift motion along the electric field induced by Alfv\'enic plasma flows. A guiding-center description is used to reveal the role of drift motions during the bulk nonthermal energization. The nonthermally dominated acceleration resulting from magnetic reconnection may help explain the highly efficient electron acceleration in solar flares and other astrophysical systems.
This is an introduction to models of accretion discs around black holes. After a presentation of the non-relativistic equations describing the structure and evolution of geometrically thin accretion discs we discuss their steady-state solutions and compare them to observation. Next we describe in detail the thermal-viscous disc instability model and its application to dwarf novae for which it was designed and its X-ray irradiated-disc version which explains the soft X--ray transients, i.e. outbursting black-hole low-mass X-ray binaries. We then turn to the role of advection in accretion flow onto black holes illustrating its action and importance with a toy model describing both ADAFs and slim discs. We conclude with a presentation of the general-relativistic formalism describing accretion discs in the Kerr space-time.
The Milankovitch theory states that the orbital eccentricity, precession, and obliquity of the Earth influence our climate by modulating the summer insolation at high latitudes in the northern hemisphere. Despite considerable success of this theory in explaining climate change over the Pleistocene epoch (2.6 to 0.01 Myr ago), it is inconclusive with regard to which combination of orbital elements paced the 100 kyr glacial-interglacial cycles over the late Pleistocene. Here we explore the role of the orbital elements in pacing the Pleistocene deglaciations by modeling ice-volume variations in a Bayesian approach. When comparing models, this approach takes into account the uncertainties in the data as well as the different degrees of model complexity. We find that the Earth's obliquity (axial tilt) plays a dominant role in pacing the glacial cycles over the whole Pleistocene, while precession only becomes important in pacing major deglaciations after the transition of the dominant period from 41 kyr to 100 kyr (the mid-Pleistocene transition). We also find that geomagnetic field and orbital inclination variations are unlikely to have paced the Pleistocene deglaciations. We estimate that the mid-Pleistocene transition took place over a 220 kyr interval centered on a time 715 kyr ago, although the data permit a range of 600--1000 kyr. This transition, occurring within just two 100\,kyr cycles, indicates a relatively rapid change in the climate response to insolation.
We examine the cosmological correlators induced by the simultaneous breaking of parity and of statistical isotropy, e.g., in presence of the coupling ${\cal L} = f(\phi) ( - \frac{1}{4} F^2 + \frac{\gamma}{4} F \tilde{F} )$ between the inflaton $\phi$ and a vector field with vacuum expectation value ${\bf A}$. For a suitably chosen function $f$, the energy in the vector field $\rho_{\rm A}$ does not decay during inflation. This results in nearly scale-invariant signatures of broken statistical isotropy and parity. Specifically, we find that the scalar-scalar correlator of primordial curvature perturbations includes a quadrupolar anisotropy, $P_\zeta ( {\bf k}) = P(k)[ 1 + g_* ( \hat{\bf k} \cdot \hat{\bf A})^2]$, and a (angle-averaged) scalar bispectrum that is a linear combination of the first $3$ Legendre polynomials, $B_\zeta(k_1, k_2, k_3) = \sum_L c_L P_L (\hat{\bf k}_1 \cdot \hat{\bf k}_2) P(k_1) P(k_2) + 2~{\rm perms} $, with $c_0 : c_1 : c_2 = 2 : -3 : 1$ ($c_1 \neq 0$ is a consequence of parity violation, corresponding to the constant $\gamma \neq 0$). The latter is one of the main results of this paper, which provides for the first time a clear example of an inflationary model where a non-negligible $c_1$ contribution to the bispectrum is generated. The scalar-tensor and tensor-tensor correlators induce characteristic signatures in the Cosmic Microwave Background temperature anisotropies (T) and polarization (E/B modes); namely, non-diagonal contributions to $\langle a_{\ell_1 m_1} a_{\ell_2 m_2}^* \rangle$, with $|\ell_1 - \ell_2| = 1$ in TT, TE, EE and BB, and $|\ell_1 - \ell_2| = 2$ in TB and EB. The latest CMB bounds on the scalar observables ($g_*$, $c_0$, $c_1$ and $c_2$), translate into the upper limit $\rho_{\rm A} / \rho_\phi \lesssim 10^{-9}$ at $\gamma=0$. We find that the upper limit on the vector energy density becomes much more stringent as $\gamma$ grows.
We use astrophysical and atomic clock tests of the stability of the fine-structure constant $\alpha$, together with Type Ia supernova and Hubble parameter data, to constrain the simplest class of dynamical dark energy models where the same degree of freedom is assumed to provide both the dark energy and (through a dimensionless coupling, $\zeta$, to the electromagnetic sector) the $\alpha$ variation. We show how current data tightly constrains a combination of $\zeta$ and the dark energy equation of state $w_0$. At the $95\%$ confidence level and marginalizing over $w_0$ we find $|\zeta|<5\times10^{-6}$, with the atomic clock tests dominating the constraints. The forthcoming generation of high-resolution ultra-stable spectrographs will enable significantly tighter constraints.
There exist a variety of star-galaxy classification techniques, each with their own strengths and weaknesses. In this paper, we present a novel meta-classification framework that combines and fully exploits different techniques to produce a more robust star-galaxy classification. To demonstrate this hybrid, ensemble approach, we combine a purely morphological classifier, a supervised machine learning method based on random forest, an unsupervised machine learning method based on self-organizing maps, and a hierarchical Bayesian template fitting method. Using data from the CFHTLenS survey, we consider different scenarios: when a high-quality training set is available with spectroscopic labels from DEEP2, SDSS, VIPERS, and VVDS, and when the demographics of sources in a low-quality training set do not match the demographics of objects in the test data set. We demonstrate that our Bayesian combination technique improves the overall performance over any individual classification method in these scenarios. Thus, strategies that combine the predictions of different classifiers may prove to be optimal in currently ongoing and forthcoming photometric surveys, such as the Dark Energy Survey and the Large Synoptic Survey Telescope.
We model the multi-wavelength emission in the southern hotspot of the radio quasar 4C74.26. The synchrotron radio emission is resolved near the shock with the MERLIN radio-interferometer, and the rapid decay of this emission behind the shock is interpreted as the decay of the amplified downstream magnetic field as expected for small scale turbulence. Electrons are accelerated to only 0.3 TeV, consistent with a diffusion coefficient many orders of magnitude greater than in the Bohm regime. If the same diffusion coefficient applies to the protons, their maximum energy is only ~100 TeV.
Warm Neptune- and sub-Neptune-sized exoplanets in orbits smaller than Mercury's are thought to have experienced extensive atmospheric evolution. Here we propose that a potential outcome of this atmospheric evolution is the formation of helium-dominated atmospheres. The hydrodynamic escape rates of Neptune- and sub-Neptune-sized exoplanets are comparable to the diffusion-limited escape rate of hydrogen, and therefore the escape is heavily affected by diffusive separation between hydrogen and helium. A helium atmosphere can thus be formed -- from a primordial hydrogen-helium atmosphere -- via atmospheric hydrodynamic escape from the planet. The helium atmosphere has very different abundances of major carbon and oxygen species from those of a hydrogen atmosphere, leading to distinctive transmission and thermal emission spectral features. In particular, the hypothesis of a helium-dominated atmosphere can explain the thermal emission spectrum of GJ 436 b, a warm Neptune-sized exoplanet, while also consistent with the transmission spectrum. This model atmosphere contains trace amounts of hydrogen, carbon, and oxygen, with the predominance of CO over CH4 as the main form of carbon. With our atmospheric evolution model, we find that if the mass of the initial atmosphere envelope is 1E-3 planetary mass, hydrodynamic escape can reduce the hydrogen abundance in the atmosphere by several orders of magnitude in ~10 billion years. Observations of exoplanet transits may thus detect signatures of helium atmospheres and probe the evolutionary history of small exoplanets.
We present a result of a blind search for [CII] 158 $\mu$m emitters at $z\sim 4.5$ using ALMA Cycle~0 archival data. We collected extra-galactic data covering at 330-360 GHz (band~7) from 8 Cycle~0 projects from which initial results have been already published. The total number of fields is 243 and the total on-source exposure time is 19.2 hours. We searched for line emitters in continuum-subtracted data cubes with spectral resolutions of 50, 100, 300 and 500 km/s. We could not detect any new line emitters above a 6-$\sigma$ significance level. This result provides upper limits to the [CII] luminosity function at $z\sim 4.5$ over $L_{\rm [CII]} \sim 10^8 - 10^{10} L_{\odot}$ or star formation rate, SFR $\sim$ 10-1000 M$_{^\odot}$/yr. These limits are at least 2 orders of magnitude larger than the [CII] luminosity functions expected from the $z \sim 4$ UV luminosity function or from numerical simulation. However, this study demonstrates that we would be able to better constrain the [CII] luminosity function and to investigate possible contributions from dusty galaxies to the cosmic star-formation rate density by collecting Cycle~1+2 archival data as the ALMA Patchy Deep Survey.
There has been a persistent conundrum in attempts to model the nucleosynthesis of heavy elements by rapid neutron capture (the $r$-process). Although the location of the abundance peaks near nuclear mass numbers 130 and 195 identify an environment of rapid neutron capture near closed nuclear shells, the abundances of elements just above and below those peaks are often underproduced by more than an order of magnitude in model calculations. At the same time there is a debate in the literature as to what degree the $r$-process elements are produced in supernovae or the mergers of binary neutron stars. In this paper we propose a novel solution to both problems. We demonstrate that the underproduction of elements above and below the $r$-process peaks characteristic in the main or weak $r$-process events (like magnetohydrodynamic jets or neutrino-driven winds in core-collapse supernovae) can be supplemented via fission fragment distributions from the recycling of material in a neutron-rich environment such as that encountered in neutron star mergers. In this paradigm, the abundance peaks themselves are well reproduced by a moderately neutron rich, main $r$-process environment such as that encountered in the magnetohydrodynamical jets in supernovae supplemented with a high-entropy, weakly neutron rich environment such as that encountered in the neutrino-driven-wind model to produce the lighter $r$-process isotopes. Moreover, we show that the relative contributions to the $r$-process abundances in both the solar-system and metal-poor stars from the weak, main, and fission-recycling environments required by this proposal are consistent with estimates of the relative Galactic event rates of core-collapse supernovae for the weak and main $r$-process and neutron star mergers for the fission-recycling $r$-process.
The extent of nucleosynthesis in models of type I X-ray bursts and the associated impact on the energy released in these explosive events are sensitive to nuclear masses and reaction rates around the $^{64}$Ge waiting point. Using a recent high precision mass measurement of $^{65}$As along with large-scale shell model calculations, we have determined new thermonuclear rates of the $^{64}$Ge($p$,$\gamma$)$^{65}$As and $^{65}$As($p$,$\gamma$)$^{66}$Se reactions. We examine the impact of available rates for these two reactions through a representative one-zone X-ray burst model. We find that our recommended rates may strongly suppress the flow of abundances toward $A\approx100$, in sharp contrast to recent work claiming that $^{64}$Ge is not a significant $rp$-process waiting point. Indeed, the summed mass fractions for species with $A > 70$ varies by about factors of 3 or 2 depending upon the adopted $^{64}$Ge($p$,$\gamma$)$^{65}$As or $^{65}$As($p$,$\gamma$)$^{66}$Se rates, respectively. Furthermore, the predictions for nuclear energy generation rate E$_\mathrm{gen}$ at late times during the burst varies rather significantly between the models using the different rates, with differences as large as about a factor of 2.
Scientific exploitation of the ever increasing volumes of astronomical data requires efficient and practical methods for data access, visualisation, and analysis. Hierarchical sky tessellation techniques enable a multi-resolution approach to organising data on angular scales from the full sky down to the individual image pixels. Aims. We aim to show that the Hierarchical progressive survey (HiPS) scheme for describing astronomical images, source catalogues, and three-dimensional data cubes is a practical solution to managing large volumes of heterogeneous data and that it enables a new level of scientific interoperability across large collections of data of these different data types. Methods. HiPS uses the HEALPix tessellation of the sphere to define a hierarchical tile and pixel structure to describe and organise astronomical data. HiPS is designed to conserve the scientific properties of the data alongside both visualisation considerations and emphasis on the ease of implementation. We describe the development of HiPS to manage a large number of diverse image surveys, as well as the extension of hierarchical image systems to cube and catalogue data. We demonstrate the interoperability of HiPS and Multi-Order Coverage (MOC) maps and highlight the HiPS mechanism to provide links to the original data. Results. Hierarchical progressive surveys have been generated by various data centres and groups for ~200 data collections including many wide area sky surveys, and archives of pointed observations. These can be accessed and visualised in Aladin, Aladin Lite, and other applications. HiPS provides a basis for further innovations in the use of hierarchical data structures to facilitate the description and statistical analysis of large astronomical data sets.
We present a detailed description of the physical properties of our current census of white dwarfs within 40 pc of the Sun, based on an exhaustive spectroscopic survey of northern hemisphere candidates from the SUPERBLINK proper motion database. Our method for selecting white dwarf candidates is based on a combination of theoretical color-magnitude relations and reduced proper motion diagrams. We reported in an earlier publication the discovery of nearly 200 new white dwarfs, and we present here the discovery of an additional 133 new white dwarfs, among which we identify 96 DA, 3 DB, 24 DC, 3 DQ, and 7 DZ stars. We further identify 178 white dwarfs that lie within 40 pc of the Sun, representing a 40% increase of the current census, which now includes 492 objects. We estimate the completeness of our survey at between 66 and 78%, allowing for uncertainties in the distance estimates. We also perform a homogeneous model atmosphere analysis of this 40 pc sample and find a large fraction of massive white dwarfs, indicating that we are successfully recovering the more massive, and less luminous objects often missed in other surveys. We also show that the 40 pc sample is dominated by cool and old white dwarfs, which populate the faint end of the luminosity function, although trigonometric parallaxes will be needed to shape this part of the luminosity function more accurately. Finally, we identify 4 probable members of the 20 pc sample, 4 suspected double degenerate binaries, and we also report the discovery of two new ZZ Ceti pulsators.
Wavelength-dependent point spread functions (PSFs) violate an implicit assumption in current galaxy shape measurement algorithms that deconvolve the PSF measured from stars (which have stellar spectral energy distributions (SEDs)) from images of galaxies (which have galactic SEDs). Since the absorption length of silicon depends on wavelength, CCDs are a potential source of PSF chromaticity. Here we develop two toy models to estimate the sensitivity of the cosmic shear survey from the Large Synoptic Survey Telescope to chromatic effects in CCDs. We then compare these toy models to simulated estimates of PSF chromaticity derived from the LSST photon simulator PhoSim. We find that even though sensor contributions to PSF chromaticity are subdominant to atmospheric contributions, they can still significantly bias cosmic shear results if left uncorrected, particularly in the redder filter bands and for objects that are off-axis in the field of view.
EK Draconis (HD 129333: G1.5 V) is a well-known young (50 Myr) solar analog. In 2012, Hubble Space Telescope returned to EK Dra to follow up a far-ultraviolet (FUV) SNAPshot visit by Cosmic Origins Spectrograph (COS) two years earlier. The brief SNAP pointing had found surprisingly redshifted, impulsively variable subcoronal "hot-line" emission of Si IV 140 nm (T~ 80,000 K). Serendipitously, the 2012 follow-on program witnessed one of the largest FUV flares ever recorded on a sunlike star, which again displayed strong redshifts (downflows) of 30-40 km/s, even after compensating for small systematics in the COS velocity scales, uncovered through a cross-calibration by Space Telescope Imaging Spectrograph (STIS). The (now reduced, but still substantial) ~10 km/s hot-line redshifts outside the flaring interval did not vary with rotational phase, so cannot be caused by "Doppler Imaging" (bright surface patches near a receding limb). Density diagnostic O IV] 140 nm multiplet line ratios of EK Dra suggest log(Ne)~ 11, an order of magnitude larger than in low-activity solar twin Alpha Centauri A, but typical of densities inferred in large stellar soft X-ray events. The self-similar FUV hot-line profiles between the flare decay and the subsequent more quiet periods, and the unchanging but high densities, reinforce a long-standing idea that the coronae of hyperactive dwarfs are flaring all the time, in a scale-free way; a "flare-ona" if you will. In this picture, the subsonic hot-line downflows probably are a byproduct of the post-flare cooling process, something like "coronal rain" on the Sun. All in all, the new STIS/COS program documents a complex, energetic, dynamic outer atmosphere of the young sunlike star.
The Sun provides the energy necessary to sustain our existence. While the Sun provides for us, it is also capable of taking away. The weather and climatic scales of solar evolution and the Sun-Earth connection are not well understood. There has been tremendous progress in the century since the discovery of solar magnetism - magnetism that ultimately drives the electromagnetic, particulate and eruptive forcing of our planetary system. There is contemporary evidence of a decrease in solar magnetism, perhaps even indicators of a significant downward trend, over recent decades. Are we entering a minimum in solar activity that is deeper and longer than a typical solar minimum, a "grand minimum"? How could we tell if we are? What is a grand minimum and how does the Sun recover? These are very pertinent questions for modern civilization. In this paper we present a hypothetical demonstration of entry and exit from grand minimum conditions based on a recent analysis of solar features over the past 20 years and their possible connection to the origins of the 11(-ish) year solar activity cycle.
The collisionless dissipation of anisotropic Alfv\'enic turbulence is a promising candidate to solve the solar wind heating problem. Extensive studies examined the kinetic properties of Alfv\'en waves in simple Maxwellian or bi-Maxwellian plasmas. However, the observed electron velocity distribution functions in the solar wind are more complex. In this study, we analyze the properties of kinetic Alfv\'en waves in a plasma with two drifting electron populations. We numerically solve the linearized Maxwell-Vlasov equations and find that the damping rate and the proton-electron energy partition for kinetic Alfv\'en waves are significantly modified in such plasmas, compared to plasmas without electron drifts. We suggest that electron drift is an important factor to take into account when considering the dissipation of Alfv\'enic turbulence in the solar wind or other $\beta \sim 1$ astrophysical plasmas.
We present the results of a Hubble Space Telescope Wide Field Camera 3 imaging survey of 11 of the lowest mass brown dwarfs in the Pleiades known (25-40 Mjup). These objects represent the predecessors to T dwarfs in the field. Using a semi-empirical binary PSF-fitting technique, we are able to probe to 0.03" (0.75 pixel), better than 2x the WFC3/UVIS diffraction limit. We did not find any companions to our targets. From extensive testing of our PSF-fitting method on simulated binaries, we compute detection limits which rule out companions to our targets with mass ratios of $q\gtrsim0.7$ and separations $a\gtrsim4$ AU. Thus, our survey is the first to attain the high angular resolution needed to resolve brown dwarf binaries in the Pleiades at separations that are most common in the field population. We constrain the binary frequency over this range of separation and mass ratio of 24-40 Mjup Pleiades brown dwarfs to be <11% for 1$\sigma$ (<26% at 2$\sigma$). This binary frequency is consistent with both younger and older brown dwarfs in this mass range.
In the article we present a hardware meteor detector. The detection principle is based on the electromagnetic wave reflection from the ionized meteor trail in the atmosphere. The detector uses the ANADIGM field programmable analogue array (FPAA), which is an attractive alternative for a typically used detecting equipment - a PC computer with dedicated software. We implement an analog signal path using most of available FPAA resources to obtain precise audio signal detection. Our new detector was verified in collaboration with the Polish Fireball Network - the organization which monitors meteor activity in Poland. When compared with currently used signal processing PC software employing real radio meteor scatter signals, our low-cost detector proved to be more precise and reliable. Due to its cost and efficiency superiority over the current solution, the presented module is going to be implemented in the planned distributed detectors system.
We present the physical properties of [OIII] emission line galaxies at z>3 as the tracers of active galaxies at 1Gyr before the peak epoch at z~2. We have performed deep narrow-band imaging surveys in the Subaru/XMM-Newton Deep Survey Field with MOIRCS on the Subaru Telescope and have constructed coherent samples of 34 [OIII] emitters at z=3.2 and 3.6, as well as 107 H$\alpha$ emitters at z=2.2 and 2.5. We investigate their basic physical quantities, such as stellar masses, star formation rates (SFRs), and sizes using the publicly available multi-wavelength data and high resolution images by the Hubble Space Telescope. The stellar masses and SFRs show a clear correlation known as the "main sequence" of star-forming galaxies. It is found that the location of the main sequence of the [OIII] emitters at z=3.2 and 3.6 is almost identical to that of the H$\alpha$ emitters at z=2.2 and 2.5. Also, we investigate their mass-size relation and find that the relation does not change between the two epochs. When we assume that the star-forming galaxies at z=3.2 grow simply along the same main sequence down to z=2.2, galaxies with $M_* = 10^{9}$-$10^{11} M_{\odot}$ increase their stellar masses significantly by a factor of 10-2. They climb up the main sequence, and their star formation rates also increase a lot as their stellar masses grow. This indicates that star formation activities of galaxies are accelerated from z>3 towards the peak epoch of galaxy formation at z~2.
We present the results of multi wavelength, co-spatial and near co-temporal observations of jets above a sunspot light bridge. The data were obtained with the Solar Optical Telescope (SOT) on board Hinode, the Interface Region Spectrograph (IRIS) and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO). Most of the jets in the Ca II H images show decreasing brightness with height while in the IRIS slit jaw images at 1330 \AA~ jets show a bright leading edge. These jets show rising and falling motion as evident from the parabolic profile obtained from the time-distance diagram. The rising and falling speeds of the jets are similar. These jets show a coordinated behaviour between neighbouring jets moving jointly up and down. Some of the jets show a plasma ejection from the leading edge which is also hotter at the transition region (TR) and coronal temperatures . A Similar behaviour is seen in the AIA wave bands that suggests that jets above the LB reach up to the lower corona and the leading edges are heated up to coronal temperatures. Such jets are important means of transfer mass and energy to the transition region and corona above sunspots.
In this paper we report 23 magnetic field measurements of the B3IV star HD 23478: 12 obtained from high resolution Stokes $V$ spectra using the ESPaDOnS (CFHT) and Narval (TBL) spectropolarimeters, and 11 from medium resolution Stokes $V$ spectra obtained with the DimaPol spectropolarimeter (DAO). HD 23478 was one of two rapidly rotating stars identified as potential "centrifugal magnetosphere" hosts based on IR observations from the Apache Point Observatory Galactic Evolution Experiment survey. We derive basic physical properties of this star including its mass ($M=6.1^{+0.8}_{-0.7}\,M_\odot$), effective temperature ($T_{\rm eff}=20\pm2\,$kK), radius ($R=2.7^{+1.6}_{-0.9}\,R_\odot$), and age ($\tau_{\rm age}=3^{+37}_{-1}\,$Myr). We repeatedly detect weakly-variable Zeeman signatures in metal, He and H lines in all our observations corresponding to a longitudinal magnetic field of $\langle B_z\rangle\approx-2.0\,$kG. The rotational period is inferred from Hipparcos photometry ($P_{\rm rot}=1.0498(4)\,$d). Under the assumption of the Oblique Rotator Model, our obsevations yield a surface dipole magnetic field of strength $B_d\geq9.5\,$kG that is approximately aligned with the stellar rotation axis. We confirm the presence of strong and broad H$\alpha$ emission and gauge the volume of this star's centrifugal magnetosphere to be consistent with those of other H$\alpha$ emitting centrifugal magnetosphere stars based on the large inferred Alfv\'en to Kepler radius ratio.
Evolutionary models of pre-main sequence stars remain largely uncalibrated, especially for masses below that of the Sun, making each new dynamical mass and radius measurement a valuable test of theoretical models. Stellar mass dependent features of star formation (such as disk evolution, planet formation, and even the IMF) are fundamentally tied to these models, which implies a systematic uncertainty that can only be improved with precise measurements of calibrator stars. We present the discovery that UScoCTIO 5, a known spectroscopic binary (P = 34 days, Mtot sin(i) = 0.64 Msun), is an eclipsing system with both primary and secondary eclipses apparent in K2 light curves obtained during Campaign 2. We have simultaneously fit the eclipse profiles from the K2 light curves and the existing RV data to demonstrate that UScoCTIO 5 consists of a pair of nearly identical M4.5 stars with M_A = 0.329 +/- 0.002 Msun, R_A = 0.834 +/- 0.006 Rsun, M_B = 0.317 +/- 0.002 Msun, and R_B = 0.810 +/- 0.006 Rsun. The radii are broadly consistent with pre-main sequence ages predicted by stellar evolutionary models, but none agree to within the uncertainties. All models predict systematically incorrect masses at the 25--50% level for the HR diagram position of these mid-M dwarfs, suggesting significant modifications for any trend of other properties that vary with stellar mass. The form of the discrepancy for most model sets is not that they predict luminosities that are too low, but rather that they predict temperatures that are too high, suggesting that the models do not fully encompass the physics of energy transport (via convection and/or missing opacities) and/or a miscalibration of the SpT-Teff scale. The simplest modification to the models (changing Teff to match observations) would yield an older age for this system, in line with the recently proposed older age of Upper Scorpius (~11 Myr).
Massive B-type stars with strong magnetic fields and fast rotation are very
rare and provide a mystery for theories of both star formation and magnetic
field evolution. Only two such stars, called sigma Ori E analogs, were
previously known. Recently, a team involved in APOGEE, one of the Sloan Digital
Sky Survey III programs, announced the discovery of two additional rigidly
rotating magnetosphere stars, HD23478 and HD345439. The presence of magnetic
fields in these newly discovered sigma Ori E analogs was not investigated in
the past.
In the framework of our ESO Large Programme, and one normal ESO programme, we
carried out low-resolution FORS2 spectropolarimetric observations of HD23478
and HD345439.
From the measurements using hydrogen lines, we discover a rather strong
longitudinal magnetic field of the order of up to 1.5kG in HD23478, and up to
1.3kG using the entire spectrum. The analysis of HD345439 using four subsequent
spectropolarimetric subexposures does not reveal the presence of a magnetic
field at a significance level of 3sigma. On the other hand, the inspection of
individual subexposures indicates that HD345439 may host a strong magnetic
field, rapidly varying over 88 minutes. A hint at the fast rotation of HD345439
is also given by the behaviour of several metallic and He I lines in the
low-resolution FORS2 spectra, showing profile variations already on such a
short time scale.
We study the possibility of realizing the growth rate of matter density perturbations lower than that in General Relativity. Using the approach of the effective field theory of modified gravity encompassing theories beyond Horndeski, we derive the effective gravitational coupling $G_{\rm eff}$ and the gravitational slip parameter $\eta$ for the perturbations deep inside the Hubble radius. In Horndeski theories the necessary condition for achieving weak gravity is that the tensor propagation speed $c_{\rm t}$ is smaller than the speed of light, but this is not a sufficient condition due to the presence of a scalar-matter interaction that always enhances $G_{\rm eff}$. Beyond the Horndeski domain it is possible to realize $G_{\rm eff}$ smaller than the Newton gravitational constant $G$, while the scalar and tensor perturbations satisfy no-ghost and stability conditions. We present a concrete dark energy scenario with varying $c_{\rm t}$ and numerically study the evolution of perturbations to confront the model with the observations of red-shift space distortions and weak lensing.
The credibility of an eclipse timing variation (ETV) diagram analysis is investigated for various manifestations of the mass transfer and gravitational radiation processes in binary systems. The monotonicity of the period variations and the morphology of the respective ETV diagrams are thoroughly explored in both the direct impact and the accretion disk mode of mass transfer, accompanied by different types of mass and angular momentum losses (through a hot-spot emission from the gainer and via the L2/L3 points). Mass transfer rates comparable to or greater than 10^{-8} M_sun/yr are measurable for typical noise levels of the ETV diagrams, regardless of whether the process is conservative. However, the presence of a transient disk around the more massive component defines a critical mass ratio q_cr ~ 0.83 above which the period turns out to decrease when still in the conservative regime, rendering the measurability of the anticipated variations a much more complicated task. The effects of gravitational radiation proved to be rather undetectable, except for systems with physical characteristics that only refer to cataclysmic variables. Unlike the hot-spot effects, the Lagrangian points L2 and L3 support very efficient routes of strong angular momentum loss. It is further shown that escape of mass via the L3 point - when the donor is the less massive component - safely provides critical mass ratios above which the period is expected to decrease, no matter how intense the process is.
Usually the effects of isotropic inhomogeneities are not seriously taken into account in the determination of the cosmological parameters because of Copernican principle whose statement is that we do not live in the privileged domain in the universe. But Copernican principle has not been observationally confirmed yet in sufficient accuracy, and there is the possibility that there are non-negligible large-scale isotropic inhomogeneities in our universe. In this paper, we study the effects of the isotropic inhomogeneities on the determination of the cosmological parameters and show the probability that non-Copernican isotropic inhomogeneities mislead us into believing, for example, the phantom energy of the equation of state, $p=w\rho$ with $w<-1$, even in case that $w=-1$ is the true value.
We have studied the impact of dust feedback on the survival and structure of vortices in protoplanetary discs using 2-D shearing box simulations with Lagrangian dust particles. We consider dust with a variety of sizes (stopping time $t_s = 10^{-2}\Omega^{-1} - 10^{2}\Omega^{-1}$, from fully coupled with the gas to the decoupling limit. We find that a vortex is destroyed by dust feedback when the total dust-to-gas mass ratio within the vortex is larger than 30-50%, independent of the dust size. The dust distribution can still be asymmetric in some cases after the vortex has been destroyed. With smaller amounts of dust, a vortex can survive for at least 100 orbits, and the maximum dust surface density within the vortex can be more than 100 times larger than the gas surface density, potentially facilitating planetesimal formation. On the other hand, in these stable vortices, small ($t_s < \Omega^{-1}$) and large ($t_s > \Omega^{-1}$) dust grains concentrate differently and affect the gas dynamics in different ways. The distribution of large dust is more elongated than that of small dust. Large dust ($t_s > \Omega^{-1}$) concentrates in the centre of the vortex and feedback leads to turn-over in vorticity towards the centre, forming a quiescent region within an anticyclonic vortex. Such a turn-over is absent if the vortex is loaded with small grains. We demonstrate that, in protoplanetary discs where both large and small dust grains are present and under the right condition, the concentration of large dust towards the vortex centre can lead to a quiescent centre, repelling the small dust and forming a small dust ring around the vortex centre. Such anticorrelations between small and large dust within vortices may explain the discrepancy between ALMA and near-IR scattered light observations in the asymmetric region of transitional discs.
We report results of the CO(J=3-2) and HCO+(J=4-3) observations of the W40 HII region with the ASTE 10 m telescope (HPBW~22 arcsec) to search for molecular outflows and dense clumps.We found that the velocity field in the region is highly complex, consisting of at least four distinct velocity components at V LSR ~ 3, 5, 7, and 10 km/s. The ~7 km/s component represents the systemic velocity of cold gas surroundingthe entire region, and causes heavy absorption in the CO spectra over the velocity range 6 <V LSR< 9 km/s. The ~5 and ~10 km/s components exhibit high CO temperature (>40 K) and are found mostly around the HII region, suggesting that these components are likely to be tracing dense gas interacting with the expanding shell around the HII region. Based on the CO data, we identified 13 regions of high velocity gas which we interpret as candidate outflow lobes. Using the HCO+ data, we also identified six clumps and estimated their physical parameters. On the basis of the ASTE data and near-infrared images from 2MASS, we present an updated three-dimensional model of this region. In order to investigate molecular outflows in W40, the SiO (J=1-0, v=0) emission line and some other emission lines at 40 GHz were also observed with the 45 m telescope at the Nobeyama Radio Observatory, but they were not detected at the present sensitivity.
Detection of the global redshifted 21~cm signal is an excellent means of deciphering the physical processes during the Dark Ages and subsequent Epoch of Reionization (EoR). However, the detection of this faint signal is challenging due to the high precision required in instrumental calibration and modeling of substantially brighter foregrounds and instrumental systematics. In particular, modeling and removal of receiver noise with mK accuracy remains a formidable task in experiments aiming to detect the global signal using single-element spectral radiometers. Interferometers do not respond to receiver noise; therefore, we explore here the theory of the response of interferometers to global signals. We first derive the response to uniform sky of interferometers made of element antennas, then extend the analysis to interferometers made of 1-D arrays and finally consider 2-D aperture antennas. The analysis suggests that short-spacing interferometers made of omnidirectional antennas have the best sensitivity for the detection of the global signal. These antennas may be wideband 1-D arrays of dipoles or 1-D aperture antennas; we argue that the interferometer is best configured to be EW with the elements oriented NS. However, the performance of such interferometers would be limited by crosstalk, mode-coupling of foreground continuum sources into spectral confusion, and uncertainty in its telescope filter function. We conclude that the only useful interferometer for the global EoR signal is an interferometer between two collinear wideband short dipoles, with sensitivity to the global signal realized by placing a space beam splitter between the elements to form a Zero-Spacing Interferometer.
We revisit the evidence for a "dust cloud" observed by the Cassini spacecraft at Saturn in 2006. The simultaneous data of 3 instruments are compared to interpret the signatures of a coherent swarm of dust that could have remained floating near the equatorial plane. The conspicuous pattern, as seen in the dust counters of the Cosmic Dust Analyser (CDA) and in the magnetic field (MAG), clearly repeats on three consecutive revolutions of the spacecraft. The data of the Radio Plasma and Wave Science (RPWS) appear less decisive but do back our conclusions. The results support the idea of a "magnetic bubble" as reported from both Voyager flybys in the early 1980ies. That particular cloud, which we firstly discovered in the CDA data, is estimated to about 1.36 Saturnian radii in size, and probably broadening. Both the bulk of dust particles and the peak of the magnetic depression seem to drift apart, but this can also be an effect of hitting the cloud at different parts during the traverse.
Providing an early warning of galactic supernova explosions from neutrino signals is important in studying supernova dynamics and neutrino physics. A dedicated supernova trigger system has been designed and installed in the data acquisition system at Daya Bay and integrated into the worldwide Supernova Early Warning System (SNEWS). Daya Bay's unique feature of eight identically-designed detectors deployed in three separate experimental halls makes the trigger system naturally robust against cosmogenic backgrounds, enabling a prompt analysis of online triggers and a tight control of the false-alert rate. The trigger system is estimated to be fully sensitive to 1987A-type supernova bursts throughout most of the Milky Way. The significant gain in sensitivity of the eight-detector configuration over a mass-equivalent single detector is also estimated. The experience of this online trigger system is applicable to future projects with spatially distributed detectors.
We report the detection of 19 new high-energy gamma-ray sources closely associated with BL Lac objects at energies higher than 10 GeV in the 6.3 years Fermi-Large Area Telescope sky, selected using the Minimum Spanning Tree (MST) clustering method. Photon clusters with good selection parameters were found matching the positions of known blazars in the fifth Roma-BZCAT catalogue. A brief summary of the properties of these sources is presented.
The neutrino mechanism of core-collapse supernova is investigated via non-relativistic, two-dimensional (2D), neutrino radiation-hydrodynamic simulations. For the transport of electron flavor neutrinos, we use the interaction rates defined by Bruenn (1985) and the isotropic diffusion source approximation (IDSA) scheme, which decomposes the transported particles into trapped particle and streaming particle components. Heavy neutrinos are described by a leakage scheme. Unlike the "ray-by-ray" approach in other multi-dimensional IDSA implementations in spherical coordinates, we use cylindrical coordinates and solve the trapped particle component in multiple dimensions, improving the proto-neutron star resolution and the neutrino transport in angular and temporal directions. We perform Newtonian 1D and 2D ab initio simulations from prebounce core collapse to several hundred milliseconds postbounce with 11, 15, 21, and 27 $M_\odot$ progenitors from Woosley et al.~(2002) with the HS(DD2) equation of state. We obtain robust explosions with diagnostic energies $E_{\rm dig} \gtrsim 0.1- 0.5$~B for all considered 2D models within approximately $100-300$ milliseconds after bounce and find that explosions are mostly dominated by the neutrino-driven convection, although standing accretion shock instabilities are observed as well. We also find that the level of electron deleptonization during collapse dramatically affects the postbounce evolution, e.g.~the ignorance of neutrino-electron scattering during collapse will lead to a stronger explosion.
We present T-PHOT, a publicly available software aimed at extracting accurate photometry from low resolution images of deep extragalactic fields, where the blending of sources can be a serious problem for the accurate and unbiased measurement of fluxes and colours. T-PHOT has been developed within the ASTRODEEP project and it can be considered as the next generation to TFIT, providing significant improvements above it and other similar codes. T-PHOT gathers data from a high resolution image of a region of the sky, and uses it to obtain priors for the photometric analysis of a lower resolution image of the same field. It can handle different types of datasets as input priors: i) a list of objects that will be used to obtain cutouts from the real high resolution image; ii) a set of analytical models; iii) a list of unresolved, point-like sources, useful e.g. for far infrared wavelength domains. We show that T-PHOT yields accurate estimations of fluxes within the intrinsic uncertainties of the method, when systematic errors are taken into account (which can be done thanks to a flagging code given in the output). T-PHOT is many times faster than similar codes like TFIT and CONVPHOT (up to hundreds, depending on the problem and the method adopted), whilst at the same time being more robust and more versatile. This makes it an optimal choice for the analysis of large datasets. In addition we show how the use of different settings and methods significantly enhances the performance. Given its versatility and robustness, T-PHOT can be considered the preferred choice for combined photometric analysis of current and forthcoming extragalactic UV to FIR imaging surveys. [abridged]
Spicules are ubiquitous, fast moving jets observed off-limb in chromospheric spectral lines. Combining the recently-launched Interface Region Imaging Spectrograph with the Solar Dynamics Observatory and Hinode, we have a unique opportunity to study spicules simultaneously in multiple passbands and from a seeing free environment. This makes it possible to study their thermal evolution over a large range of temperatures. A recent study showed that spicules appear in several chromospheric and transition region spectral lines, suggesting that spicules continue their evolution in hotter passbands after they fade from Ca II H. In this follow-up paper we answer some of the questions that were raised in the introductory study. In addition, we study spicules off-limb in C II 1330 {\AA} for the first time. We find that Ca II H spicules are more similar to Mg II 2796 {\AA} spicules than initially reported. For a sample of 54 spicules, we find that 44% of Si IV 1400 {\AA} spicules are brighter toward the top; 56% of the spicules show an increase in Si IV emission when the Ca II H component fades. We find several examples of spicules that fade from passbands other than Ca II H, and we observe that if a spicule fades from a passband, it also generally fades from the passbands with lower formation temperatures. We discuss what these new, multi-spectral results mean for the classification of type I and II spicules.
The 2011 outburst of the black hole candidate IGR J17091-3624 followed the canonical track of state transitions along with the evolution of Quasi-Periodic Oscillation (QPO) frequencies before it began exhibiting various variability classes similar to GRS 1915+105. We use this canonical evolution of spectral and temporal properties to determine the mass of IGR J17091-3624, using three different methods, viz : Photon Index ($\Gamma$) - QPO frequency ($\nu$) correlation, QPO frequency ($\nu$) - Time (day) evolution and broadband spectral modelling based on Two Component Advective Flow. We provide a combined mass estimate for the source using a Naive Bayes based joint likelihood approach. This gives a probable mass range of 11.8 M$_{\odot}$ - 13.7 M$_{\odot}$. Considering each individual estimate and taking the lowermost and uppermost bounds among all three methods, we get a mass range of 8.7 M$_{\odot}$ - 15.6 M$_{\odot}$ with 90% confidence. We discuss the probable implications of our findings in the context of two component accretion flow.
We present new nonlinear, time-dependent convective hydrodynamical models of RR Lyrae stars computed assuming a constant helium-to-metal enrichment ratio and a broad range in metal abundances (Z=0.0001--0.02). The stellar masses and luminosities adopted to construct the pulsation models were fixed according to detailed central He burning Horizontal Branch evolutionary models. The pulsation models cover a broad range in stellar luminosity and effective temperatures and the modal stability is investigated for both fundamental and first overtones. We predict the topology of the instability strip as a function of the metal content and new analytical relations for the edges of the instability strip in the observational plane. Moreover, a new analytical relation to constrain the pulsation mass of double pulsators as a function of the period ratio and the metal content is provided. We derive new Period-Radius-Metallicity relations for fundamental and first-overtone pulsators. They agree quite well with similar empirical and theoretical relations in the literature. From the predicted bolometric light curves, transformed into optical (UBVRI) and near-infrared (JHK) bands, we compute the intensity-averaged mean magnitudes along the entire pulsation cycle and, in turn, new and homogenous metal-dependent (RIJHK) Period-Luminosity relations. Moreover, we compute new dual and triple band optical, optical--NIR and NIR Period-Wesenheit-Metallicity relations. Interestingly, we find that the optical Period-W(V,B-V) is independent of the metal content and that the accuracy of individual distances is a balance between the adopted diagnostics and the precision of photometric and spectroscopic datasets.
During an extensive multiwavelength campaign that we performed in 2013-14 the prototypical Seyfert 1 galaxy NGC 5548 has been found in an unusual condition of heavy and persistent obscuration. The newly discovered "obscurer" absorbs most of the soft X-ray continuum along our line of sight and lowers the ionizing luminosity received by the classical warm absorber. Here we present the analysis of the high resolution X-ray spectra collected with XMM-Newton and Chandra throughout the campaign, which are suitable to investigate the variability of both the obscurer and the classical warm absorber. The time separation between these X-ray observations range from 2 days to 8 months. On these timescales the obscurer is variable both in column density and in covering fraction. This is consistent with the picture of a patchy wind. The most significant variation occurred in September 2013 when the source brightened for two weeks. A higher and steeper intrinsic continuum and a lower obscurer covering fraction are both required to explain the spectral shape during the flare. We suggest that a geometrical change of the soft X-ray source behind the obscurer cause the observed drop in the covering fraction. Due to the higher soft X-ray continuum level the September 2013 Chandra spectrum is the only X ray spectrum of the campaign where individual features of the warm absorber could be detected. The spectrum shows absorption from Fe-UTA, \ion{O}{iv}, and \ion{O}{v}, consistent to belong to the lower-ionization counterpart of the historical NGC 5548 warm absorber. Hence, we confirm that the warm absorber has responded to the drop in the ionizing luminosity caused by the obscurer.
This paper presents SunPy (version 0.5), a community-developed Python package for solar physics. Python, a free, cross-platform, general-purpose, high-level programming language, has seen widespread adoption among the scientific community, resulting in the availability of a large number of software packages, from numerical computation (NumPy, SciPy) and machine learning (scikit-learn) to visualisation and plotting (matplotlib). SunPy is a data-analysis environment specialising in providing the software necessary to analyse solar and heliospheric data in Python. SunPy is open-source software (BSD licence) and has an open and transparent development workflow that anyone can contribute to. SunPy provides access to solar data through integration with the Virtual Solar Observatory (VSO), the Heliophysics Event Knowledgebase (HEK), and the HELiophysics Integrated Observatory (HELIO) webservices. It currently supports image data from major solar missions (e.g., SDO, SOHO, STEREO, and IRIS), time-series data from missions such as GOES, SDO/EVE, and PROBA2/LYRA, and radio spectra from e-Callisto and STEREO/SWAVES. We describe SunPy's functionality, provide examples of solar data analysis in SunPy, and show how Python-based solar data-analysis can leverage the many existing tools already available in Python. We discuss the future goals of the project and encourage interested users to become involved in the planning and development of SunPy.
One of the key ingredients of the Unified Model of Active Galactic Nuclei (AGN) is the presence of a torus-like optically thick medium composed by dust and gas around the putative supermassive black hole. However, the structure, size and composition of this circumnuclear medium are still matter of debate. To this end, the search for column density variations through X-ray monitoring on different timescales (months, weeks and few days) is fundamental to constrain size, kinematics and location of the X-ray absorber(s). Here we describe our project of mining the Swift-XRT archive to assemble a sample of AGN with extreme column density variability and determining the physical properties of the X-ray absorber(s). We also present the results obtained from a daily-weekly Swift-XRT follow-up monitoring recently performed on one of the most interesting new candidates for variability discovered so far, Mrk 915.
The structural and dynamical properties of star clusters are generally derived by means of the comparison between steady-state analytic models and the available observables. With the aim of studying the biases of this approach, we fitted different analytic models to simulated observations obtained from a suite of direct N-body simulations of star clusters in different stages of their evolution and under different levels of tidal stress to derive mass, mass function and degree of anisotropy. We find that masses can be under/over-estimated up to 50% depending on the degree of relaxation reached by the cluster, the available range of observed masses and distances of radial velocity measures from the cluster center and the strength of the tidal field. The mass function slope appears to be better constrainable and less sensitive to model inadequacies unless strongly dynamically evolved clusters and a non-optimal location of the measured luminosity function are considered. The degree and the characteristics of the anisotropy developed in the N-body simulations are not adequately reproduced by popular analytic models and can be detected only if accurate proper motions are available. We show how to reduce the uncertainties in the mass, mass-function and anisotropy estimation and provide predictions for the improvements expected when Gaia proper motions will be available in the near future.
Rich in HII regions, giant molecular clouds are natural laboratories to study massive stars and sequential star formation. The Galactic star forming complex W33 is located at l=~12.8deg and at a distance of 2.4 kpc, has a size of ~10 pc and a total mass of (~0.8 - ~8.0) X 10^5 Msun. The integrated radio and IR luminosity of W33 - when combined with the direct detection of methanol masers, the protostellar object W33A, and protocluster embedded within the radio source W33 main - mark the region out as a site of vigorous ongoing star formation. In order to assess the long term star formation history, we performed an infrared spectroscopic search for massive stars, detecting for the first time fourteen early-type stars, including one WN6 star and four O4-7 stars. The distribution of spectral types suggests that this population formed during the last ~2-4 Myr, while the absence of red supergiants precludes extensive star formation at ages 6-30 Myr. This activity appears distributed throughout the region and does not appear to have yielded the dense stellar clusters that characterize other star forming complexes such as Carina and G305. Instead, we anticipate that W33 will eventually evolve into a loose stellar aggregate, with Cyg OB2 serving as a useful, albeit richer and more massive, comparator. Given recent distance estimates, and despite a remarkably similar stellar population, the rich cluster Cl 1813-178 located on the north-west edge of W33 does not appear to be physically associated with W33.
Using observations of Spitzer/IRAC, we report the serendipitous discovery of excess infrared emission from a single white dwarf PG 0010+280. At a temperature of 27,220 K and a cooling age of 16 Myr, it is the hottest and youngest white dwarf to display an excess at 3-8 $\mu$m. The infrared excess can be fit by either an opaque dust disk within the tidal radius of the white dwarf or a 1300 K blackbody, possibly from an irradiated substellar object or a re-heated giant planet. PG 0010+280 has two unique properties that are different from white dwarfs with a dust disk: (i) relatively low emission at 8 $\mu$m and (ii) non-detection of heavy elements in its atmosphere from high-resolution spectroscopic observations with Keck/HIRES. The origin of the infrared excess remains unclear.
We describe here a focal plane array of Cold-Electron Bolometer (CEB) detectors integrated in a Frequency Selective Surface (FSS) for the 350 GHz detection band of the OLIMPO balloon-borne telescope. In our architecture, the two terminal CEB has been integrated in the periodic unit cell of the FSS structure and is impedance matched to the embedding impedance seen by it and provides a resonant interaction with the incident sub-mm radiation. The detector array has been designed to operate in background noise limited condition for incident powers of 20 pW to 80 pW, making it possible to use the same pixel in both photometric and spectrometric configurations. We present high frequency and dc simulations of our system, together with fabrication details. The frequency response of the FSS array, optical response measurements with hot/cold load in front of optical window and with variable temperature black body source inside cryostat are presented. A comparison of the optical response to the CEB model and estimations of Noise Equivalent power (NEP) is also presented.
We have observed regions of three galaxy clusters at z$\sim$ [0.06, 0.09]
(Abell85, Abell1205, Abell2440), as well as calibration sources with the Nancay
radiotelescope (NRT) to search for 21 cm emission and fully characterize the
FPGA based BAORadio digital backend. The total observation time of few hours
per source have been distributed over few months, from March 2011 to January
2012, due to scheduling constraints of the NRT, which is a transit telescope.
Data have been acquired in parallel with the NRT standard correlator (ACRT)
back-end, as well as with the BAORadio data acquisition system. The latter
enables wide band instantaneous observation of the [1250, 1500]MHz frequency
range, as well as the use of powerful RFI mitigation methods thanks to its fine
time sampling. A number of questions related to instrument stability, data
processing and calibration are discussed. We have obtained the radiometer
curves over the integration time range [0.01,10 000] seconds and we show that
sensitivities of few mJy over most of the wide frequency band can be reached
with the NRT.
It is clearly shown that in blind line search, which is the context of HI
intensity mapping for Baryon Acoustic Oscillations, the new acquisition system
and processing pipeline outperforms the standard one. We report a positive
detection of 21 cm emission at 3 sigma-level from galaxies in the outer region
of Abell85 at 1352 MHz (14 400 km/s) corresponding to a line strength of 0.8 Jy
km/s. We observe also an excess power around 1318 MHz (21 600 km/s), although
at lower statistical significance, compatible with emission from Abell1205
galaxies. Detected radio line emissions have been cross matched with optical
catalogs and we have derived hydrogen mass estimates.
Long-term records of sunspot number and concentrations of cosmogenic radionuclides (10Be and 14C) on the Earth reveal the variation of the Sun's magnetic activity over hundreds and thousands of years. We identify several clear periods in sunspot, 10Be, and 14C data as 1000, 500, 350, 200 and 100 years. We found that the periods of the first five spherical harmonics of the slow magnetic Rossby mode in the presence of a steady toroidal magnetic field of 1200-1300 G in the lower tachocline are in perfect agreement with the time scales of observed variations. The steady toroidal magnetic field can be generated in the lower tachocline either due to the steady dynamo magnetic field for low magnetic diffusivity or due to the action of the latitudinal differential rotation on the weak poloidal primordial magnetic field, which penetrates from the radiative interior. The slow magnetic Rossby waves lead to variations of the steady toroidal magnetic field in the lower tachocline, which modulate the dynamo magnetic field and consequently the solar cycle strength. This result constitutes a key point for long-term prediction of the cycle strength. According to our model, the next deep minimum in solar activity is expected during the first half of this century.
The statistical validation of transiting exoplanets proved to be an efficient technique to secure the nature of small exoplanet signals which cannot be established by purely spectroscopic means. However, the spectroscopic diagnoses are providing us with useful constraints on the presence of blended stellar contaminants. In this paper, we present how a contaminating star affects the measurements of the various spectroscopic diagnoses as function of the parameters of the target and contaminating stars using the model implemented into the PASTIS planet-validation software. We find particular cases for which a blend might produce a large radial velocity signal but no bisector variation. It might also produce a bisector variation anti-correlated with the radial velocity one, as in the case of stellar spots. In those cases, the full width half maximum variation provides complementary constraints. These results can be used to constrain blend scenarios for transiting planet candidates or radial velocity planets. We review all the spectroscopic diagnoses reported in the literature so far, especially the ones to monitor the line asymmetry. We estimate their uncertainty and compare their sensitivity to blends. Based on that, we recommend the use of BiGauss which is the most sensitive diagnosis to monitor line-profile asymmetry. In this paper, we also investigate the sensitivity of the radial velocities to constrain blend scenarios and develop a formalism to estimate the level of dilution of a blended signal. Finally, we apply our blend model to re-analyse the spectroscopic diagnoses of HD16702, an unresolved face-on binary which exhibits bisector variations.
One of the most peculiar characteristics of Active Galactic Nuclei (AGN) is
their variability over all wavelengths. This property has been used in the past
to select AGN samples and is foreseen to be one of the detection techniques
applied in future multi-epoch surveys, complementing photometric and
spectroscopic methods.
In this paper, we aim to construct and characterise an AGN sample using a
multi-epoch dataset in the r band from the SUDARE-VOICE survey.
Our work makes use of the VST monitoring program of an area surrounding the
Chandra Deep Field South to select variable sources. We use data spanning a six
month period over an area of 2 square degrees, to identify AGN based on their
photometric variability.
The selected sample includes 175 AGN candidates with magnitude r < 23 mag. We
distinguish different classes of variable sources through their lightcurves, as
well as X-ray, spectroscopic, SED, optical and IR information overlapping with
our survey.
We find that 12% of the sample (21/175) is represented by SN. Of the
remaining sources, 4% (6/154) are stars, while 66% (102/154) are likely AGNs
based on the available diagnostics. We estimate an upper limit to the
contamination of the variability selected AGN sample of about 34%, but we point
out that restricting the analysis to the sources with available
multi-wavelength ancillary information, the purity of our sample is close to
80% (102 AGN out of 128 non-SN sources with multi-wavelength diagnostics). Our
work thus confirms the efficiency of the variability selection method in
agreement with our previous work on the COSMOS field; in addition we show that
the variability approach is roughly consistent with the infrared selection.
While Lorentz invariance, the fundamental symmetry of Einstein's theory of General Relativity, has been tested to a great level of detail, Grand Unified Theories that combine gravity with the other three fundamental forces may result in a violation of Lorentz symmetry at the Planck scale. These energies are unattainable experimentally. However, minute deviations from Lorentz invariance may still be present at much lower energies. These deviations can accumulate over large distances, making astrophysical measurements the most sensitive tests of Lorentz symmetry. One effect of Lorentz invariance violation is an energy dependent photon dispersion of the vacuum resulting in differences of the light travel time from distant objects. The Standard-Model Extension (SME) is an effective theory to describe the low-energy behaviour of a more fundamental Grand Unified Theory, including Lorentz and CPT violating terms. In the SME the Lorentz violating operators can in part be classified by their mass-dimension d, with the lowest order being d = 5. However, measurements of photon polarization have constrained operators with d = 5 setting lower limits on the energy at which they become dominant well beyond the Planck scale. On the other hand, these operators also violate CPT, and thus d = 6 could be the leading order. In this paper we present constraints on all 25 real coefficients describing anisotropic non-birefringent Lorentz invariance violation at mass dimension d = 6 in the SME. We used Fermi-LAT observations of 25 active galactic nuclei to constrain photon dispersion and combined our results with previously published limits in order to simultaneously constrain all 25 coefficients. This represents the first set of constraints on these coefficients of a mass-dimension d = 6, whereas previous measurements were only able to constrain linear combinations of all 25 coefficients.
Estimations of black hole spin in the three Galactic microquasars GRS 1915+105, GRO J1655-40, and XTE J1550-564 have been carried out based on spectral and timing X-ray measurements and various theoretical concepts. Among others, a non-linear resonance between axisymmetric epicyclic oscillation modes of an accretion disc around a Kerr black hole has been considered as a model for the observed high-frequency quasi-periodic oscillations (HF QPOs). Estimates of spin predicted by this model have been derived based on the geodesic approximation of the accreted fluid motion. Here we assume accretion flow described by the model of a pressure-supported torus and carry out related corrections to the mass-spin estimates. We find that for dimensionless black hole spin a<0.9, the resonant eigenfrequencies are very close to those calculated for the geodesic motion. Their values slightly grow with increasing torus thickness. These findings agree well with results of a previous study carried out in the pseudo-Newtonian approximation. The situation becomes different for a>0.9, in which case the resonant eigenfrequencies rapidly decrease as the torus thickness increases. We conclude that the assumed non-geodesic effects shift the lower limit of the spin, implied for the three microquasars by the epicyclic model and independently measured masses, from a~0.7 to a~0.6. Their consideration furthermore confirms compatibility of the model with the rapid spin of GRS 1915+105 and provides highly testable predictions of the QPO frequencies. Individual sources with a moderate spin (a<0.9) should exhibit a smaller spread of the measured 3:2 QPO frequencies than sources with a near-extreme spin (a~1). This should be further examined using the large amount of high-resolution data expected to become available with the next generation of X-ray instruments, such as the proposed Large Observatory for X-ray Timing (LOFT).
Active galactic nuclei (AGNs) and quasars are important astrophysical objects to understand. Recently, microlensing observations have constrained the size of the quasar X-ray emission region to be of the order of 10 gravitational radii of the central supermassive black hole. For distances within a few gravitational radii, light paths are strongly bent by the strong gravity field of the central black hole. If the central black hole has nonzero angular momentum (spin), a photon's polarization plane will be rotated by the gravitational Faraday effect. The observed X-ray flux and polarization will then be influenced significantly by the strong gravity field near the source. Consequently, linear gravitational lensing theory is inadequate for such extreme circumstances. We present simple algorithms computing strong lensing effects of Kerr black holes, including effects on polarization. Our algorithms are realized in a program "KERTAP" in two versions: MATLAB and Python. The key ingredients of KERTAP are: a graphic user interface, a {\it backward} ray-tracing algorithm, a polarization propagator dealing with gravitational Faraday rotation, and algorithms computing observables such as flux magnification and polarization angles. Our algorithms can be easily realized in other programming languages such as FORTRAN, C, and C++. The MATLAB version of KERTAP is parallelized using the MATLAB Parallel Computing Toolbox and the Distributed Computing Server. The Python code was sped up using Cython and supports full implementation of MPI using 'mpi4py' package. As an example, we investigate the inclination angle dependence of the observed polarization and the strong lensing magnification of AGN X-ray emission. We conclude that it is possible to perform complex numerical-relativity-related computations using interpreted languages such as MATLAB and Python.
Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars. Aims. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit. The search pro- cedure is made on a N x N grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the companion in the number of sigmas. The code CANDID also incorporates a robust method to set a 3{\sigma} detection limit on the flux ratio, which is based on an analytical injection of a fake companion at each point in the grid. We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting V1334 Cyg and AX Cir were detected at > 13 sigmas. The companion around AW Per is detected at more than 15 sigmas with a flux ratio of f = 1.22 +/- 0.30 %. We made a possible detection of the companion orbiting RT Aur with f = 0.22 +/- 0.11 %. It was detected at 3.8{\sigma} using the closure phases only, and so more observations are needed to confirm the detection. We also set the detection limit for possible undetected companions. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively.
We study numerically the anisotropy of the cosmic ray (CR) flux emitted by a single source calculating the trajectories of individual CRs. We show that the contribution of a single source to the observed anisotropy is instead determined solely by the fraction the source contributes to the total CR intensity, its age and its distance,and does not depend on the CR energy at late times. Therefore the observation of a constant dipole anisotropy indicates that a single source dominates the CR flux in the corresponding energy range. A natural explanation for the plateau between 2--20 TeV observed in the CR anisotropy is thus the presence of a single, nearby source. For the source age of 2 Myr, as suggested by the explanation of the antiproton and positron data from PAMELA and AMS-02 through a local source [arXiv:astro-ph/1504.06472], we determine the source distance as $\sim 200$ pc. Combined with the contribution of the global CR sea calculated in the escape model, we can explain qualitatively the data for the dipole anisotropy. Our results suggest that the assumption of a smooth CR source distribution should be abandoned between 200 GeV and 1 PeV.
We examine the capability of a coronal flare ignited blast wave scenario to reproduce the chromospheric phenomenon. We numerically simulate the Moreton event of December 06, 2006 considering both the corona and the chromosphere. To obtain a sufficiently strong coronal shock -able to generate a detectable chromospheric Moreton wave- a relatively low magnetic field intensity is required, in comparison with the active region values. Employing reasonable coronal constraints, we show that the flare ignited blast wave scenario is capable to reproduce the observations.
Observations from IRIS allow us to study the chromospheric heating and evaporation processes during solar flares with unprecedented high resolution and high cadence. We track the complete evolution of ~11 MK evaporation flows in the M1.1 flare on 2014 September 6 and the X1.6 flare on 2014 September 10. These hot flows, as indicated by the entirely blueshifted Fe xxi 1354.08 emission line, evolve smoothly with a velocity decreasing from ~200 km/s to almost stationary within a few minutes. The velocity decrease appears to be exponential in time, especially for the X1.6 flare. There is a good correlation between the flow velocity and the energy deposition rate as represented by the RHESSI hard X-Ray flux, or the time derivative of the soft X-Ray flux observed by GOES and the HINODE/XRT, which is in general agreement with models of nonthermal electron heating. The maximum blue shift of Fe xxi appears approximately at the same time as or slightly after the impulsive enhancement of the UV continuum and the Mg ii 2798.8 line emission, demonstrating that the hot evaporation flow is closely related to heating of the lower chromosphere. In the M1.1 flare, the maximum Fe XXI blue shift, the peaks of RHESSI hard X-Ray flux and the time derivative of GOES flux occur at roughly the same time. In the X1.6 flare, the maximum blue shift is found slightly before the peak time of the soft X-Ray derivative, a result that is qualitatively consistent with a recent radiative hydrodynamic simulation (RADYN) of chromospheric evaporation. Finally, while the hot Fe xxi 1354.08 line is entirely blueshifted with no obvious rest component, cool chromospheric and transition region lines like Si iv 1402.77 are often not entirely redshifted but just reveal an obvious enhancement in the red wing of the line profiles at the flare ribbons, suggesting a faster speed of chromospheric condensation than previously thought.
The steep power-law dependence of element abundance enhancements on the mass-to-charge ratios [A/Q] of the ions in impulsive solar energetic-particle (SEP) events causes these enhancements to reflect the temperature-dependent pattern of Q of the ions in the source plasma. We searched for SEP events from coronal plasma that is hotter or cooler than the limited region of 2.5 - 3.2 MK previously found to dominate 111 impulsive SEP events. Fifteen new events were found, four (three) originated in 2-MK (4-MK) plasma, but none from outside this temperature range. Although the impulsive SEP events are strongly associated with flares, this result indicates that these ions are not accelerated from flare-heated plasma, which can often exceed 10 MK. Evidently the ions of 2 - 20 MeV/amu that we observe in space are accelerated from active-region plasma on open magnetic-field lines near the flare, but not from the closed loops of the flare. The power-law dependence of the abundance enhancements on A/Q of the ions is expected from theoretical models of acceleration from regions of magnetic reconnection.
VISTA variables in the Via Lactea is an ESO Public survey dedicated to scan the bulge and an adjacent portion of the Galactic disk in the fourth quadrant using the VISTA telescope and the near-infrared camera VIRCAM. One of the leading goals of the VVV survey is to contribute to the knowledge of the star cluster population of the Milky Way. To improve the census of the Galactic star clusters, we performed a systematic scan of the JHKs images of the Galactic plane section of the VVV survey. Our detection procedure is based on a combination of superficial density maps and visual inspection of promising features in the NIR images. The material examined are color-composite images corresponding to the DR1 of VVV. We report the discovery of 493 new star cluster candidates. The analysis of the spatial distribution show that the clusters are very concentrated in he Galactic plane, presenting some local maxima around the position of large star-forming complexes, such as G305, RCW 95, and RCW 106. The vast majority of the cluster candidates are quite compact and generally surrounded by bright and/or dark nebulosities. IRAS point sources are associated with 59% of the sample, while 88% are associated with MSX point sources. GLIMPSE 8 mum images of the cluster candidates show a variety of morphologies, with 292 clusters dominated by knotty sources, while 361 clusters show some kind of nebulosity. Spatial cross-correlation with young stellar objects, masers, and extended green-object catalogs suggest that a large sample of the new cluster candidates are extremely young. In particular, 104 star clusters associated to methanol masers are excellent candidates for ongoing massive star formation. Also, there is a special set of sixteen cluster candidates that present clear signspot of star-forming activity having associated simultaneosly dark nebulae, young stellar objects, EGOs, and masers.
We present a manifestly local, diffeomorphism invariant and locally Poincare invariant formulation of vacuum energy sequestering. In this theory, quantum vacuum energy generated by matter loops is cancelled by auxiliary fields. The auxiliary fields decouple from gravity almost completely. Their only residual effect is an {\it a priori} arbitrary, finite contribution to the curvature of the background geometry, which is radiatively stable. Its value is to be determined by a measurement, like the finite part of any radiatively stable UV-sensitive quantity in quantum field theory
Percolation theory is applied to the phase transition dynamics of domain pattern formation in segregating quasi-two-dimensional binary Bose--Einstein condensates. Our numerical experiments revealed that the percolation threshold is close to 0.5. A long-range open domain wall appears with a fractal dimension between two percolating domains. Such a wall can survive for a long time as a relic of the phase transition according to the dynamic finite-size-scaling hypothesis, which seems to be in contrast to the current understanding in cosmology that an infinite defect violates a scale invariance.
We study the proposal by Mersini et al. that the observed dark energy might be explained by the back-reaction of the set of tail modes in a theory with a dispersion relation in which the mode frequency decays exponentially in the trans-Planckian regime. The matter tail modes are frozen out, however they induce metric fluctuations. The energy-momentum tensor with which the tail modes effect the background geometry obtains contributions from both metric and matter fluctuations. We calculate the equation of state induced by the tail modes taking into account the gravitational contribution. We find that, in contrast to the case of frozen super-Hubble cosmological fluctuations, in this case the matter perturbations dominate, and they yield an equation of state which to leading order takes the form of a positive cosmological constant.
We review the status of String Gas Cosmology after the 2015 Planck data release. String gas cosmology predicts an almost scale-invariant spectrum of cosmological perturbations with a slight red tilt, like the simplest inflationary models. It also predicts a scale-invariant spectrum of gravitational waves with a slight blue tilt, unlike inflationary models which predict a red tilt of the gravitational wave spectrum. String gas cosmology yields two consistency relations which determine the tensor to scalar ratio and the slope of the gravitational wave spectrum given the amplitude and tilt of the scalar spectrum. We show that these consistency relations are in good agreement with the Planck data. We discuss future observations which will be able to differentiate between the predictions of inflation and those of string gas cosmology.
ArDM-1t is the first operating ton-scale liquid argon detector for direct search of Dark Matter particles. Developed at CERN as Recognized Experiment RE18, the experiment has been approved in 2010 to be installed in the Spanish underground site LSC (Laboratorio Subterraneo de Canfranc). Under the label of LSC EXP-08-2010 the ArDM detector underwent an intensive period of technical completion and safety approval until the recent filling of the target vessel with almost 2 ton of liquid argon. This report describes the experimental achievements during commissioning of ArDM and the transition into a stage of first physics data taking in single phase operational mode. We present preliminary observations from this run. A first indication for the background discrimination power of LAr detectors at the ton-scale is shown. We present an outlook for completing the detector with the electric drift field and upgrade of the scintillation light readout system with novel detector modules based on SiPMs in order to improve the light yield.
During inflation, scalar fields, including the Higgs boson, may acquire a nonzero vacuum expectation value, which must later relax to the equilibrium value during reheating. In the presence of the time-dependent condensate, the vacuum state can evolve into a state with a nonzero particle number. We show that, in the presence of lepton number violation in the neutrino sector, the particle production can explain the observed matter-antimatter asymmetry of the universe. We find that this form of leptogenesis is particularly effective when the Higgs condensate decays rapidly and at low reheat temperature. As part of the calculation, we present some exact results for the Bogoliubov transformations for Majorana fermions with a nonzero time-dependent chemical potential, in addition to a time-dependent mass.
It is known that after Affleck-Dine baryogenesis, spatial inhomogeneities of Affleck-Dine field grow into non-topological solitons called Q-balls. In gauge mediated SUSY breaking models, sufficiently large Q-balls with baryon charge are stable while Q-balls with lepton charge can always decay into leptons. For a Q-ball that carries nonzero $B$ and $L$ charges, the difference between the baryonic component and the leptonic component in decay rate may induce nonzero electric charge on the Q-ball. This implies that charged Q-ball, also called gauged Q-ball, may emerge in our universe. In this paper, we investigate two complex scalar fields, a baryonic scalar field and a leptonic one, in an Abelian gauge theory. We find stable solutions of gauged Q-balls for different baryon and lepton charges. Those solutions shows that a Coulomb potential arises and the Q-ball becomes electrically charged as expected. It is energetically favored that some amount of leptonic component decays, but there is an upper bound on its amount due to the Coulomb force. The baryonic decay also becomes possible by virtue of electrical repulsion and we find the condition to suppress it so that the charged Q-balls can survive in the universe.
We study the effect of gravitational focusing of the earth on dark matter. We find that the effect can produce a detectable diurnal modulation in the dark matter signal for part of the parameter space which for high dark matter masses is larger than the diurnal modulation induced by the fluctuations in the flux of dark matter particles due to the rotation of the earth around its own axis. The two sources of diurnal modulation have different phases and can be distinguished from each other. We demonstrate that the diurnal modulation can potentially check the self-consistency of experiments that observe annual modulated signals that can be attributed to dark matter. Failing to discover a daily varying signal can result conclusively to the falsification of the hypothesis that the annual modulation is due to dark matter. We also suggest that null result experiments should check for a daily modulation of their rejected background signal with specific phases. A potential discovery could mean that dark matter collisions have been vetoed out.
In de Rham-Gabadadze-Tolley (dRGT) massive gravity and bigravity, a non-minimal matter coupling involving both metrics generically re-introduces the Boulware--Deser (BD) ghost. A non-minimal matter coupling via a simple, yet specific composite metric has been proposed, which eliminates the BD ghost below the strong coupling scale. Working explicitly in the metric formulation and for arbitrary spacetime dimensions, we show that this composite metric is the unique consistent non-minimal matter coupling below the strong coupling scale, which emerges out of two diagnostics, namely, absence of Ostrogradski ghosts in the decoupling limit and absence of the BD ghost from matter quantum loop corrections.
We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations [Schekochihin et al., Astrophys. J. Suppl. 182:310 (2009)]. Two main applications of these equations are magnetised (Alfvenic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge-Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor-corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag-Tang-type decaying turbulence, both in fluid and kinetic regimes.
In the majority of the literature on plasma shock waves until now, electrons have played the role of "ghost particles," since they contribute to mass- and momentum flows only negligibly and have been treated as taking care of the electric plasma neutrality. In some more recent papers, however, electrons play a new important role in the shock dynamics and thermodynamics, especially at the solar-wind termination shock. They react on the shock electric field in a very specific way, leading to suprathermal non-equilibrium distributions of the downstream electrons that can be represented by a kappa distribution function. In this article, we discuss why these anticipated hot electron population has not been seen by the plasma detectors of the Voyager spacecraft downstream of the solar-wind termination shock. We show that hot non-equilibrium electrons induce a strong negative electric charge-up of any spacecraft cruising through this downstream plasma environment. This charge reduces electron fluxes at the spacecraft detectors to non-detectable intensities. Furthermore, we show that the Debye length $\lambda_{\mathrm D}^{\kappa}$ grows to values of about $\lambda_{\mathrm D}^{\kappa}/\lambda_{\mathrm D}\simeq 10^{6}$ compared to the classical value $\lambda_{\mathrm D}$ in this hot-electron environment. This unusual condition allows for the propagation of a certain type of electrostatic plasma waves that, at very large wavelengths, allow us to determine the effective temperature of the suprathermal electrons directly by means of the phase velocity of these waves. At moderate wavelengths, the electron-acoustic dispersion relation leads to non-propagating oscillations with the ion-plasma frequency $\omega_{\mathrm p}$ instead of the traditional electron plasma frequency.
A large set of independent astronomical observations have provided a strong evidence for nonbaryonic dark matter in the Universe. One of the most investigated candidates is an unknown long-lived Weakly Interacting Massive Particle (WIMP) which was in thermal equilibrium with the primeval plasma. Here we investigate the WIMP abundance based on the relativistic kinetic treatment for gravitationally induced particle production recently proposed in the literature (Lima \& Baranov, Phys. Rev. D {\bf 90}, 043515, 2014). The new evolution equation is deduced and solved both numerically and also through a semi-analytical approach. The predictions of the WIMP observables are discussed and compared with the ones obtained in the standard approach.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
X-ray emission associated with the west lobe of the giant radio galaxy, 3C 236, was investigated with the Suzaku observatory, to evaluate the energetics in the lobe. After removing contamination from X-ray sources detected with Chandra and subtracting the X-ray and non-X-ray backgrounds, the Suzaku spectrum from the lobe was reproduced by a power-low model with a photon index of $\Gamma = 2.23_{-0.38-0.12}^{+0.44+0.14}$ where the first and second errors represent the statistical and systematic ones, respectively. Within the errors, the X-ray index was consistent with the radio synchrotron one, $\Gamma_{\rm R} = 1.74 \pm 0.07$, estimated in the 326 -- 2695 MHz range. This agreement supports that the X-ray emission is attributed to the inverse-Compton (IC) radiation from the synchrotron electrons filling the lobe, where the cosmic microwave background photons are up-scattered. This result made 3C 236 the largest radio galaxy, of which the lobe has ever been probed through the IC X-ray photons. When the photon index was fixed at $\Gamma_{\rm R}$, the X-ray flux density at 1 keV was measured as $S_{\rm X} = 12.3 \pm 2.0 \pm 1.9$ nJy. A comparison of the X-ray flux to the radio one ($S_{\rm R} = 1.11 \pm 0.02$ Jy at 608.5 MHz) yields the energy densities of the electrons and magnetic field in the west lobe as $u_{\rm e} = 3.9_{-0.7 -0.9}^{+0.6 +1.0} \times 10^{-14} $ ergs cm$^{-3}$ and $u_{\rm m} = 0.92_{-0.15 -0.35}^{+0.21 +0.52}\times 10^{-14} $ ergs cm$^{-3}$, respectively, indicating a mild electron dominance of $u_{\rm e}/u_{\rm m} = 4.2_{-1.3 -2.3}^{+1.6 +4.1}$. The latter corresponds to the magnetic field strength of $B = 0.48_{-0.04 -0.10}^{+0.05 +0.12}$ $\mu$G.These are typical among the lobes of giant radio galaxies. A compilation of the $u_{\rm e}$-size relation for the IC-detected radio galaxies implies that the west lobe of 3C 236 is still actively energized by its jet.
Among the possible alternatives to the standard cosmological model ($\Lambda$CDM), coupled Dark Energy models postulate that Dark Energy (DE), seen as a dynamical scalar field, may interact with Dark Matter (DM), giving rise to a "fifth-force", felt by DM particles only. In this paper, we study the impact of these cosmologies on the statistical properties of galaxy populations by combining high-resolution numerical simulations with semi-analytic models (SAM) of galaxy formation and evolution. New features have been implemented in the reference SAM in order to have it run self-consistently and calibrated on these cosmological simulations. They include an appropriate modification of the mass temperature relation and of the baryon fraction in DM haloes, due to the different virial scalings and to the gravitational bias, respectively. Our results show that the predictions of our coupled-DE SAM do not differ significantly from theoretical predictions obtained with standard SAMs applied to a reference $\Lambda$CDM simulation, implying that the statistical properties of galaxies provide only a weak probe for these alternative cosmological models. On the other hand, we show that both galaxy bias and the galaxy pairwise velocity distribution are sensitive to coupled DE models: this implies that these probes might be successfully applied to disentangle among quintessence, $f(R)$-Gravity and coupled DE models.
Looking at the region connecting two clusters is a promising way to identify and study the Warm-Hot Intergalactic Medium. Observations show that the spectrum of the bridge between A3556 and A3558 has a stronger soft X-ray emission than the nearby region. Suzaku observations could not discriminate the origin of the extra emission. In this work we analyze a dedicated Chandra observation of the same target to identify point sources and characterize the background emission in the bridge. We find that the count number of the point sources is much higher than average field population (using CDFS~4~Ms as a reference). Moreover, the shape of the cumulative distribution resembles that of galaxy distribution suggesting that the point sources are galaxies in a filament. The Suzaku extra emission is well explained by the high abundance of point sources identified by Chandra. Furthermore, we used optical/IR observations of point sources in the same field to estimate the density of the putative filament as rho~150 rho_b$, below Suzaku sensitivity.
Vector mode of cosmological perturbation theory imprints characteristic signals on the weak lensing signals such as curl- and B-modes which are never imprinted by the scalar mode. However, the vector mode is neglected in the standard first-order cosmological perturbation theory since it only has a decaying mode. This situation changes if the cosmological perturbation theory is expanded up to second order. The second-order vector and tensor modes are inevitably induced by the product of the first-order scalar modes. We study the effect of the second-order vector mode on the weak lensing curl- and B-modes. The curl-mode induced by the second-order vector mode dominates instead of the primordial gravitational waves when the tensor-to-scalar ratio is $r = 0.1$ and the second-order tensor mode at $\ell \geq 200$. Furthermore, the B-mode cosmic shear induced by the second-order vector mode dominates on almost all scales. However, we find that the observational signatures of the second-order vector and tensor modes cannot exceed the expected noise of ongoing and upcoming weak lensing measurements. We conclude that the curl- and B-modes induced by the second-order vector and tensor modes are unlikely to be detected in future experiments.
The luminosity of [C II] is used to probe the star formation rate in galaxies, but the correlation breaks down in some active galactic nuclei (AGNs). Models of the [C II] emission from galactic nuclei do not include the influence of X-rays on the carbon ionization balance, which may be a factor in reducing the [C II] luminosity. We calculate the [C II] luminosity in galactic nuclei under the influence of bright sources of X-rays. We solve the balance equation of the ionization states of carbon as a function of X-ray flux, electron, atomic hydrogen, and molecular hydrogen density. These are input to models of [CII] emission from the interstellar medium (ISM) in galactic nuclei. We also solve the distribution of the ionization states of oxygen and nitrogen in highly ionized regions. We find that the dense warm ionized medium (WIM) and dense photon dominated regions (PDRs) dominate the [C II] emission when no X-rays are present. The X-rays in galactic nuclei can affect strongly the C$^+$ abundance in the WIM converting some fraction to C$^{2+}$ and higher ionization states and thus reducing its [C II] luminosity. For an X-ray luminosity > 10$^{43}$ erg/s the [C II] luminosity can be suppressed by a factor of a few, and for very strong sources, >10$^{44}$ erg/s, such as found for many AGNs by an order of magnitude. Comparison of the model with extragalactic sources shows that the [C II] to far-infrared ratio declines for an X-ray luminosity >10$^{43}$ erg/s, in reasonable agreement with our model.
Kapteyn's star is an old M subdwarf believed to be a member of the Galactic halo population of stars. A recent study has claimed the existence of two super-Earth planets around the star based on radial velocity (RV) observations. The innermost of these candidate planets--Kapteyn b (P = 48 days)--resides within the circumstellar habitable zone. Given recent progress in understanding the impact of stellar activity in detecting planetary signals, we have analyzed the observed HARPS data for signatures of stellar activity. We find that while Kapteyn's star is photometrically very stable, a suite of spectral activity indices reveals a large-amplitude rotation signal, and we determine the stellar rotation period to be 143 days. The spectral activity tracers are strongly correlated with the purported RV signal of "planet b," and the 48-day period is an integer fraction (1/3) of the stellar rotation period. We conclude that Kapteyn b is not a planet in the Habitable Zone, but an artifact of stellar activity.
The quest for a $B$-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the $353\,$GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We filter the intensity and polarization maps to isolate filaments in the range of angular scales where the power asymmetry between $E$-modes and $B$-modes is observed. Using the Smoothed Hessian Major Axis Filament Finder, we identify 259 filaments at high Galactic latitude, with lengths larger or equal to $2$\deg\ (corresponding to $3.5\,$pc in length for a typical distance of $100\,$pc). These filaments show a preferred orientation parallel to the magnetic field projected onto the plane of the sky, derived from their polarization angles. We present mean maps of the filaments in Stokes $I$, $Q$, $U$, $E$, and $B$, computed by stacking individual images rotated to align the orientations of the filaments. Combining the stacked images and the histogram of relative orientations, we estimate the mean polarization fraction of the filaments to be $11\,$%. Furthermore, we show that the correlation between the filaments and the magnetic field orientations may account for the $E$ and $B$ asymmetry and the $C_{\ell}^{TE}/C_{\ell}^{EE}$ ratio, reported in the power spectra analysis of the Planck $353\,$GHz polarization maps. Future models of the dust foreground for CMB polarization studies will need to take into account the observed correlation between the dust polarization and the structure of interstellar matter.
The mapping between the distributions of the observed galaxy stellar mass and the underlying dark matter halos provides the crucial link from theories of large-scale structure formation to interpreting the complex phenomena of galaxy formation and evolution. We develop a novel statistical method, based on the Halo Occupation Distribution model (HOD), to solve for this mapping by jointly fitting the galaxy clustering and the galaxy-galaxy lensing measured from the Sloan Digital Sky Survey (SDSS). The method, called the iHOD model, extracts maximum information from the survey by including ~80% more galaxies than the traditional HOD methods, and takes into account the incompleteness of the stellar mass samples in a statistically consistent manner. The derived stellar-to-halo mass relation not only explains the clustering and lensing of SDSS galaxies over almost four decades in stellar mass, but also successfully predicts the stellar mass functions observed in SDSS. Due to its capability of modelling significantly more galaxies, the iHOD is able to break the degeneracy between the logarithmic scatter in the stellar mass at fixed halo mass and the slope of the stellar-to-halo mass relation at high mass end, without the need to assume a strong prior on the scatter and/or use the stellar mass function as an input. We detect a decline of the scatter with halo mass, from 0.22 dex at below 10^{12} Msun to 0.18 dex at 10^{14} Msun. The model also enables stringent constraints on the satellite stellar mass functions at fixed halo mass, predicting a departure from the Schechter functional form in high mass halos. The iHOD model can be easily applied to existing and future spectroscopic datasets, greatly improving the statistical constraint on the stellar-to-halo mass relation compared to the traditional HOD methods within the same survey.
Radio relics in galaxy clusters are associated with powerful shocks that (re)accelerate relativistic electrons. It is widely believed that the acceleration proceeds via diffusive shock acceleration. In the framework of thermal leakage, the ratio of the energy in relativistic electrons to the energy in relativistic protons should should be smaller than $K_{\rm e/p} \sim 10^{-2}$. The relativistic protons interact with the thermal gas to produce $\gamma$-rays in hadronic interactions. Combining observations of radio relics with upper limits from $\gamma$-ray observatories can constrain the ratio $K_{\rm e/p}$. In this work we selected 10 galaxy clusters that contain double radio relics, and derive new upper limits from the stacking of $\gamma$-ray observations by FERMI. We modelled the propagation of shocks using a semi-analytical model, where we assumed a simple geometry for shocks and that cosmic ray protons are trapped in the intracluster medium. Our analysis shows that diffusive shock acceleration has difficulties in matching simultaneously the observed radio emission and the constraints imposed by FERMI, unless the magnetic field in relics is unrealistically large ($\gg 10 ~\rm \mu G$). In all investigated cases (also including realistic variations of our basic model and the effect of re-acceleration) the mean emission of the sample is of the order of the stacking limit by FERMI, or larger. These findings put tension on the commonly adopted model for the powering of radio relics, and imply that the relative acceleration efficiency of electrons and protons is at odds with predictions of diffusive shock acceleration, requiring $K_{\rm e/p} \geq 10-10^{-2}$.
Recently, we developed a pair of meshless finite-volume Lagrangian methods for hydrodynamics: the 'meshless finite mass' (MFM) and 'meshless finite volume' (MFV) methods. These capture advantages of both smoothed-particle hydrodynamics (SPH) and adaptive mesh-refinement (AMR) schemes. Here, we extend these to include ideal magneto-hydrodynamics (MHD). The MHD equations are second-order consistent and conservative. We augment these with a divergence-cleaning scheme, which maintains div*B~0 to high accuracy. We implement these in the code GIZMO, together with a state-of-the-art implementation of SPH MHD. In every one of a large suite of test problems, the new methods are competitive with moving-mesh and AMR schemes using constrained transport (CT) to ensure div*B=0. They are able to correctly capture the growth and structure of the magneto-rotational instability (MRI), MHD turbulence, and the launching of magnetic jets, in some cases converging more rapidly than AMR codes. Compared to SPH, the MFM/MFV methods exhibit proper convergence at fixed neighbor number, sharper shock capturing, and dramatically reduced noise, div*B errors, and diffusion. Still, 'modern' SPH is able to handle most of our tests, at the cost of much larger kernels and 'by hand' adjustment of artificial diffusion parameters. Compared to AMR, the new meshless methods exhibit enhanced 'grid noise' but reduced advection errors and numerical diffusion, velocity-independent errors, and superior angular momentum conservation and coupling to N-body gravity solvers. As a result they converge more slowly on some problems (involving smooth, slowly-moving flows) but more rapidly on others (involving advection or rotation). In all cases, divergence-control beyond the popular Powell 8-wave approach is necessary, or else all methods we consider will systematically converge to unphysical solutions.
We construct models of the structural evolution of super-Earth- and mini-Neptune-type exoplanets with hydrogen-helium envelopes, incorporating radiative cooling and XUV-driven mass loss. We conduct a parameter study of these models, focusing on initial mass, radius, and envelope mass fractions, as well as orbital distance, metallicity, and the specific prescription for mass loss. From these calculations, we investigate how the observed masses and radii of exoplanets today relate to the distribution of their initial conditions. Orbital distance and initial envelope mass fraction are the most important factors determining planetary evolution, particular radius evolution. Initial mass also becomes important below a "turnoff mass," which varies with orbital distance, with mass-radius curves being approximately flat for higher masses. Initial radius is the least important parameter we study, with very little difference between the hot start and cold start limits after an age of 100 Myr. Model sets with no mass loss fail to produce results consistent with observations, but a plausible range of mass loss scenarios is allowed. In addition, we present scenarios for the formation of the Kepler-11 planets. Our best fit to observations Kepler-11b and Kepler-11c involves formation beyond the snow line, after which they moved inward, circularized, and underwent a reduced degree mass loss.
Polarization of the >~400 keV hard tail of the microquasar Cygnus X-1 has been independently reported by INTEGRAL/IBIS, and INTEGRAL/SPI and interpreted as emission from a compact jet. These conclusions were, however, based on the accumulation of all INTEGRAL data regardless of the spectral state. We utilize additional INTEGRAL exposure accumulated until December 2012, and include the AMI/Ryle (15 GHz) radio data in our study. We separate the observations into hard, soft, and intermediate/transitional states and detect radio emission from a compact jet in hard and intermediate states, but not in the soft. The 10-400 keV INTEGRAL (JEM-X and IBIS) state resolved spectra are well modeled with thermal Comptonization and reflection components. We detect a hard tail in the 0.4-2 MeV range for the hard state only. We extract the state dependent polarigrams of Cyg X-1, which all are compatible to no or undetectable level of polarization except in 400-2000 keV range in the hard state where the polarization fraction is 75$\pm$32 % and the polarization angle 40.0 +-14 deg. An upper limit on the 0.4-2 MeV soft state polarization fraction is 70%. Due to the short exposure, we obtain no meaningful constraint for the intermediate state. The likely detection of a >400 keV polarized tail in the hard state, together with the simultaneous presence of a radio jet, reinforce the notion of a compact jet origin of the 400 keV emission.
The Ly$\alpha$ line of high-redshift galaxies has emerged as a powerful probe of both early galaxy evolution and the epoch of reionization (EoR). Motivated by the upcoming wide-field survey with the Subaru Hyper Supreme-Cam (HSC), we study the angular correlation function (ACF) of narrow-band selected, $z\approx7$ LAEs. The clustering of LAEs is determined by both: (i) their typical host halo masses, $\bar{M}_{\rm h}$; (ii) the absorption due to a patchy EoR, characterized by an average neutral fraction of the IGM, $\bar{x}_{\rm HI}$. We bracket the allowed LAE ACF by exploring extreme scenarios for both the intrinsic Ly$\alpha$ emission and the morphology of cosmic ionized patches in physical EoR models. Current LAE ACF measurements imply that the Universe is mostly ionized at $z\approx7$, with $\bar{x}_{\rm HI}\lesssim0.5$ (1-$\sigma$) even for an extremely conservative model of intrinsic emission. The upcoming Ultra Deep campaign with the HSC could improve on these constraints by tens of percent, or $\bar{x}_{\rm HI}\lesssim0.3$ if the mean value of the ACF remains unchanged. The ACF at a fixed observed LAE number density and $\bar{x}_{\rm HI}$ is extremely insensitive to the EoR morphology; distinguishing between different EoR models would therefore require more accurate redshift determinations with spectroscopic follow-up observations. We also find that the low values of the currently-observed ACF implies that LAEs are hosted by relatively small dark matter halos, with $\bar{M}_{\rm h}\lesssim10^{10}M_\odot$. Combined with their observed number densities, this implies a duty cycle $\lesssim$ few per cent. These values are over an order of magnitude lower than the analogous ones for color-selected, Lyman break galaxies. This discrepancy could be due to the narrow-band LAEs searches preferentially selecting a population of young, star-burst galaxies, residing in less massive halos.
We present spectral analysis of NuSTAR and Swift observations of Cep X-4 during its outburst in 2014. We observed the source once during the peak of the outburst and once during the decay, finding good agreement in the spectral shape between the observations. We describe the continuum using a powerlaw with a Fermi-Dirac cutoff at high energies. Cep X-4 has a very strong cyclotron resonant scattering feature (CRSF) around 30 keV. A simple absorption-like line with a Gaussian optical depth or a pseudo-Lorentzian profile both fail to describe the shape of the CRSF accurately, leaving significant deviations at the red side of the line. We characterize this asymmetry with a second absorption feature around 19 keV. The line energy of the CRSF, which is not influenced by the addition of this feature, shows a small but significant positive luminosity dependence. With luminosities between (1-6)e36 erg/s, Cep X-4 is below the theoretical limit where such a correlation is expected. This behavior is similar to Vela X-1 and we discuss parallels between the two systems.
A robust classification of Cepheids into their different sub-classes and, in particular, between classical and Type II Cepheids, is necessary to properly calibrate the period-luminosity relations and for populations studies in the Galactic disc. Type II Cepheids are, however, very diverse, and classifications based either on intrinsic (period, light curve) or external parameters (e.g., [Fe/H], |z|) do not provide a unique classification. We want to ascertain the classification of two Cepheids, HQ Car and DD Vel, that are sometimes classified as classical Cepheids and sometimes as Type II Cepheids. To achieve this goal, we examine both their chemical composition and the presence of specific features in their spectra. We find emission features in the H{\alpha} and in the 5875.64 {\AA} He I lines that are typical of W Vir stars. The [Na/Fe] (or [Na/Zn]) abundances are typical of thick-disc stars, while BL Her stars are Na-overabundant ([Na/Fe]>+0.5 dex). Finally, the two Cepheids show a possible (HQ Car) or probable (DD Vel) signature of mild dust-gas separation that is usually observed only in long-period type II Cepheids and RV Tau stars. These findings clearly indicate that HQ Car and DD Vel are both Type II Cepheids from the W Vir sub-class. Several studies have reported an increase in the Cepheids' abundance dispersion towards the outer (thin) disc. A detailed inspection of the Cepheid classification, in particular for those located in the outer disc, will indicate whether this feature is real or simply an artefact of the inclusion of type II Cepheids belonging to the thick disc in the current samples.
In this work we present high resolution spectroscopic data of the giant star-forming region of N11, obtained with the GIRAFFE instrument at the Very Large Telescope. By using this data set, we find that most of the H$\alpha$ emission lines profiles in this complex can be fitted by a single Gaussian, however, multiple emission line profiles can be observed in the central region of N11. By adding all the spectra, we derive the integrated H$\alpha$ profile of this complex, which displays a width ($\sigma$) of about 12 km s$^{-1}$ (corrected by instrumental and thermal width). We find that a single Gaussian fit on the integrated H$\alpha$ profile leaves remaining wings, which can be fitted by a secondary broad Gaussian component. In addition, we find high velocity features, which spatially correlate with soft diffuse X-ray emission.
As part of the Bluedisk survey we analyse the radial gas-phase metallicity profiles of 50 late-type galaxies We compare the metallicity profiles of a sample of HI-rich galaxies against a control sample of HI-'normal' galaxies. We find the metallicity gradient of a galaxy to be strongly correlated with its HI mass fraction (M(HI) / Mstar). We note that some galaxies exhibit a steeper metallicity profile in the outer disc than in the inner disc. These galaxies are found in both the HI-rich and control samples. This contradicts a previous indication that these outer drops are exclusive to HI-rich galaxies. These effects are not driven by bars, although we do find some indication that barred galaxies have flatter metallicity profiles. By applying a simple analytical model we are able to account for the variety of metallicity profiles that the two samples present. The success of this model implies that the metallicity in these isolated galaxies may be in a local equilibrium, regulated by star formation. This insight could provide an explanation of the observed local mass-metallicity relation.
Selecting centrally quiescent galaxies from the Sloan Digital Sky Survey (SDSS) to create high signal-to-noise (>100) stacked spectra with minimal emission line contamination, we accurately and precisely model the central stellar populations of barred and unbarred quiescent disk galaxies. By splitting our sample by redshift, we can use the fixed size of the SDSS fiber to model the stellar populations at different radii within galaxies. At 0.02<z<0.04, the SDSS fiber radius corresponds to ~1 kpc, which is the typical half-light radii of both classical bulges and disky pseudobulges. Assuming that the SDSS fiber primarily covers the bulges at these redshifts, our analysis shows that there are no significant differences in the stellar populations, i.e., stellar age, [Fe/H], [Mg/Fe], and [N/Fe], of the bulges of barred vs. unbarred quiescent disk galaxies. Modeling the stellar populations at different redshift intervals from z=0.020 to z=0.085 at fixed stellar masses produces an estimate of the stellar population gradients out to about half the typical effective radius of our sample, assuming null evolution over this ~1 Gyr epoch. We find that there are no noticeable differences in the slopes of the azimuthally averaged gradients of barred vs. unbarred quiescent disk galaxies. These results suggest that bars are not a strong influence on the chemical evolution of quiescent disk galaxies.
Binary-driven hypernovae (BdHNe) following the induced gravitational collapse (IGC) paradigm have been introduced to explain the concomitance of energetic long gamma-ray bursts (GRBs) with type Ic supernovae. The progenitor system is a tight binary system composed of a carbon-oxygen (CO) core and a neutron star (NS) companion. The supernova ejecta of the exploding CO core triggers a hypercritical accretion process onto the NS, which in a few seconds reach the NS critical mass, and gravitationally collapses to a black hole (BH) emitting a GRB. These tight binary systems evolve through the supernova explosion very differently than compact binary progenitors studied in population synthesis calculations. First, the hypercritical accretion onto the NS companion alters both the mass and momentum of the binary. Second, because the explosion timescale is on par with the orbital period, the mass ejection can not be assumed to be instantaneous. Finally, the bow shock created as the accreting NS plows through the supernova ejecta can transfer angular momentum, acting as a viscosity on the orbit. These systems remain bound even if a large fraction of the binary system's mass is lost in the explosion (well above the canonical 50% limit). Because orbital velocity is high in these tight binaries, even large kicks are unlikely to unbind the system. Indeed, the outcome of BdHNe is a new family of NS-BH binaries unaccounted for in current population synthesis analyses and, although such binaries may be rare, the fact that nearly 100\% of the systems remain bound argues that they may play an important role in the compact merger rate, important for both gravitational waves, and define a new class of ultra-short GRBs.
Solar system planets move on almost circular orbits. In strong contrast, many massive gas giant exoplanets travel on highly elliptical orbits, whereas the shape of the orbits of smaller, more terrestrial, exoplanets remained largely elusive. Knowing the eccentricity distribution in systems of small planets would be important as it holds information about the planet's formation and evolution, and influences its habitability. We make these measurements using photometry from the Kepler satellite and utilizing a method relying on Kepler's second law, which relates the duration of a planetary transit to its orbital eccentricity, if the stellar density is known. Our sample consists of 28 bright stars with precise asteroseismic density measurements. These stars host 74 planets with an average radius of 2.6 $R_\oplus$. We find that the eccentricity of planets in Kepler multi-planet systems is low and can be described by a Rayleigh distribution with $\sigma$ = 0.049 $\pm$ 0.013. This is in full agreement with solar system eccentricities, but in contrast to the eccentricity distributions previously derived for exoplanets from radial velocity studies. Our findings are helpful in identifying which planets are habitable because the location of the habitable zone depends on eccentricity, and to determine occurrence rates inferred for these planets because planets on circular orbits are less likely to transit. For measuring eccentricity it is crucial to detect and remove Transit Timing Variations (TTVs), and we present some previously unreported TTVs. Finally transit durations help distinguish between false positives and true planets and we use our measurements to confirm six new exoplanets.
LCDM is remarkably successful in predicting the cosmic microwave background and large-scale structure, and LCDM parameters have been determined with only mild tensions between different types of observations. Hydrodynamical simulations starting from cosmological initial conditions are increasingly able to capture the complex interactions between dark matter and baryonic matter in galaxy formation. Simulations with relatively low resolution now succeed in describing the overall galaxy population. For example, the EAGLE simulation in volumes up to 100 cubic Mpc reproduces the observed local galaxy mass function nearly as well as semi-analytic models. It once seemed that galaxies are pretty smooth, that they generally grow in size as they evolve, and that they are a combination of disks and spheroids. But recent HST observations combined with high-resolution hydrodynamic simulations are showing that most star-forming galaxies are very clumpy; that galaxies often undergo compaction which reduces their radius and increases their central density; and that most lower-mass star-forming galaxies are not spheroids or disks but are instead elongated when their centers are dominated by dark matter. We also review LCDM challenges on smaller scales: cusp-core, "too big to fail," and substructure issues. Although starbursts can rapidly drive gas out of galaxy centers and thereby reduce the dark matter density, it remains to be seen whether this or other baryonic physics can explain the observed rotation curves of the entire population of dwarf and low surface brightness galaxies. If not, perhaps more complicated physics such as self-interacting dark matter may be needed. But standard LCDM appears to be successful in predicting the dark matter halo substructure that is now observed via gravitational lensing and breaks in cold stellar streams, and any alternative theory must do at least as well.
The basic physical parameters of a poorly studied open cluster NGC 110 and an unstudied open cluster DOLIDZE 14 are estimated in the present study using the archival PPMXL and WISE catalogues. The radius of both the clusters are estimated by fitting the modified King's empirical model on their stellar density profiles. The other basic parameters of the clusters such as distance, reddening, and age are obtained by visual fitting of the Marigo's solar metallicity isochrone on their IR colour-magnitude diagrams (CMDs). The mean-proper motion of the clusters are estimated through the individual proper motion of probable members identified through the dynamical and statistical methods. The archival catalogues (JHKW1W2) are constructed for both the clusters by compiling the extracted data from the PPMXL and WISE catalogues. The various colour-excesses, such as E(J-H), E(H-K) and E(W1-W2), are estimated using the best fit theoretical isochrone on the (J-H)-H, (H-K)-H and (W1-W2)-H CMDs, respectively. The ratios of various infrared colours of the clusters are obtained through their two-colour diagrams. We also identify the most probable members in these clusters by estimating spatial, kinematic and spatio-kinematic probabilities of stars within the cluster. A correlation between the E(H-K) and E(W1-W2) is also established.
Rotation Powered-Pulsars are subjected to long-term changes in their period of rotation, which are measured by timing observations of their rotation frequency and its derivatives ($\Omega$, $\dot{\Omega}$, $\ddot{\Omega}$). If the spin-down is solely due to dipolar radiation, the braking index should be $n=3$. To date, only a handful of braking indices have been estimated for young pulsars, and in all cases one observes that $n<3$. These observations suggest that there are complex spin-down processes taking place in the pulsar that are not fully well understood. In the present work we revisit the spin-down of young pulsars by considering a possible magnetic field growth due ohmic diffusion. In order to perform such study we perform calculations with phenomenological growth functions for the magnetic field. With that we are able to calculate the spin evolution of the neutron star with all relevant quantities. We show that such approach could explain the low values of $n$ in very young neutron stars and may be relevant to explain why a small group of neutron stars (including central compact objects) exhibit no evidence of a standard magnetic field. We find that the possibility of magnetic field growth are relevant to the spin evolution of young neutron stars, and therefore should not be neglected.
The origin of dust in a galaxy is poorly understood. Recently, the surveys of the Large Magellanic Cloud (LMC) provide astrophysical laboratories for the dust studies. By a method of population synthesis, we investigate the contributions of dust produced by asymptotic giant branch (AGB) stars, common envelope (CE) ejecta and type II supernovae (SNe II) to the total dust budget in the LMC. Based on our models, the dust production rates (DPRs) of AGB stars in the LMC are between about $2.5\times10^{-5}$ and $4.0\times10^{-6}M_\odot{\rm yr^{-1}}$. The uncertainty mainly results from different models for the dust yields of AGB stars. The DPRs of CE ejecta are about $6.3\times10^{-6}$(The initial binary fraction is 50\%). These results are within the large scatter of several observational estimates. AGB stars mainly produce carbon grains, which is consistent with the observations. Most of dust grains manufactured by CE ejecta are silicate and iron grains. The contributions of SNe II are very uncertain. Compared with SNe II without reverse shock, the DPRs of AGB stars and CE ejecta are negligible. However, if only 2 \% of dust grains produced by SNe II can survive after reverse shock, the contributions of SNe II are very small. The total dust masses produced by AGB stars in the LMC are between $2.8\times10^4$ and $3.2\times10^5M_\odot$, and those produced by CE ejecta are about $6.3\times10^4$. They are much lower than the values estimated by observations. Therefore, there should be other dust sources in the LMC.
In this study we combine the multiwavelength ultraviolet -- optical (Solar Dynamics Observatory, SDO) and radio (Nobeyama Radioheliograph, NoRH) observations to get further insight into space-frequency distribution of oscillations at different atmospheric levels of the Sun. We processed the observational data on NOAA 11711 active region and found oscillations propagating from the photospheric level through the transition region upward into the corona. The power maps of low-frequency (1--2 mHz) oscillations reproduce well the fan-like coronal structures visible in the Fe ix 171A line. High frequency oscillations (5--7 mHz) propagate along the vertical magnetic field lines and concentrate inside small-scale elements in the umbra and at the umbra-penumbra boundary. We investigated the dependence of the dominant oscillation frequency upon the distance from the sunspot barycentre to estimate inclination of magnetic tubes in higher levels of sunspots where it cannot be measured directly, and found that this angle is close to 40 degrees above the umbra boundaries in the transition region.
We report the discovery and characterisation of a new M-dwarf binary, with component masses and radii of M1 = 0.244 -0.003/+0.003 Msun, R1 = 0.261 -0.009/+0.006 Rsun, M2 = 0.179 -0.001/+0.002 Msun, R2 = 0.218 -0.011/+0.007 Rsun, and orbital period of ~4.1 days. The M-dwarf binary HATS551-027 (LP 837-20) was identified as an eclipsing binary by the HATSouth survey, and characterised by a series of high precision photometric observations of the eclipse events, and spectroscopic determinations of the atmospheric parameters and radial velocity orbits. HATS551-027 is one of few systems with both stellar components lying in the fully-convective regime of very low mass stars, and can serve as a test for stellar interior models. The radius of HATS551-027A is consistent with models to 1 sigma, whilst HATS551-027B is inflated by 9% at 2 sigma significance. We measure the effective temperatures for the two stellar components to be Teff,1 = 3190 +/- 100 K and Teff,2 = 2990+/-110 K, both are slightly cooler than theoretical models predict, but consistent with other M-dwarfs of similar masses that have previously been studied. We also measure significant Halpha emission from both components of the binary system, and discuss this in the context of the correlation between stellar activity and the discrepancies between the observed and model temperatures.
We present a complete spatial and dynamical study of the poorly populated stellar system ESO65SC03. The radial distribution of the system gives a core and cluster radii of 1.10+/-0.63 arcmin and 5.36+/-0.24 arcmin, respectively. The surface number density profile (SNDP) does not show any clear enhancement of the surface stellar number density between the stars of the system and the field regions. We derive the optimum isochrone solution for a particular grid size in the colour-magnitude diagram (CMD) using the statistical cleaning procedure. Using the statistically cleaned CMDs, we find the distance modulus, (m-M)_0, and reddening, E({B-V}), of the system to be 11.8+/-0.2 mag and 0.45 mag, respectively. The mean proper motion of this system is -5.37+/-0.81 mas/yr and 0.31+/-0.40 in RA and DEC directions, respectively. The mean proper motion of this system is found to be almost similar to the field region. The mass function for the brighter stars is found to be too high for the system to be an open cluster. These combined results place constraints on whether stellar system ESO65SC03 is a possible open star cluster remnant (POCR) or an Asterism. Our understanding is that the ESO65SC03 is in a stage of POCR by loosing their main sequence stars in the dynamic evolution processes.
We investigate the evolution of hydromagnetic perturbations in a small section of accretion disks. It is known that molecular viscosity is negligible in accretion disks. Hence, it has been argued that a mechanism, known as Magnetorotational Instability (MRI), is responsible for transporting matter in the presence of weak magnetic field. However, there are some shortcomings, which question effectiveness of MRI. Now the question arises, whether other hydromagnetic effects, e.g. transient growth (TG), can play important role to bring nonlinearity in the system, even at weak magnetic fields. Otherwise, whether MRI or TG, which is primarily responsible to reveal nonlinearity to make the flow turbulent? Our results prove explicitly that the flows with high Reynolds number (Re ), which is the case of realistic astrophysical accretion disks, exhibit nonlinearity by TG of perturbation modes faster than that by modes producing MRI. For a fixed wave vector, MRI dominates over transient effects, only at low Re , lower than its value expected to be in astrophysical accretion disks, and low magnetic fields. This seriously questions (overall) suasiveness of MRI in astrophysical accretion disks.
Theories of stellar spiral arms in disk galaxies can be grouped into two classes based on the longevity of a spiral arm. Although the quasi-stationary density wave theory supposes that spirals are rigidly-rotating, long-lived patterns, the dynamic spiral theory predicts that spirals are differentially-rotating, transient, recurrent patterns. In order to distinguish between the two spiral models from observations, we performed hydrodynamic simulations with steady and dynamic spiral models. Hydrodynamics simulations in steady spiral models demonstrated that the dust lane locations relative to the stellar spiral arms (hereafter, arm-gas offsets) depend on radius, regardless of the strength and pitch angle of the spiral and the model of the inter-stellar medium (ISM). In contrast, we found that the dynamic spiral models show no systematic radial dependence of the arm-gas offsets. The arm-gas offset radial profile method, together with the other test methods, will help us to distinguish between the two spiral models in observed spiral galaxies.
We present a mass estimate of the Planck-discovered cluster PLCK G100.2-30.4, derived from a weak lensing analysis of deep SUBARU griz images. We perform a careful selection of the background galaxies using the multi-band imaging data, and undertake the weak lensing analysis on the deep (1hr) r-band image. The shape measurement is based on the KSB algorithm; we adopt the PSFex software to model the Point Spread Function (PSF) across the field and correct for this in the shape measurement. The weak lensing analysis is validated through extensive image simulations. We compare the resulting weak lensing mass profile and total mass estimate to those obtained from our re-analysis of XMM-Newton observations, derived under the hypothesis of hydrostatic equilibrium. The total integrated mass profiles are in remarkably good agreement, agreeing within 1$\sigma$ across their common radial range. A mass $M_{500} \sim 7 x 10^{14} M_\odot$ is derived for the cluster from our weak lensing analysis. Comparing this value to that obtained from our reanalysis of XMM-Newton data, we obtain a bias factor of (1-b) = 0.8 $\pm$ 0.1. This is compatible within 1$\sigma$ with the value of (1-b) obtained by Planck Collaboration XXIV from their calibration of the bias factor using newly-available weak lensing reconstructed masses.
Filaments are considered to be basic structures and molecular clouds consist of filaments. Filaments are often observed as extending in the direction perpendicular to the interstellar magnetic field. The structure of filaments has been studied based on a magnetohydrostatic equilibrium model (Tomisaka 2014). Here, we simulate the expected polarization pattern for isothermal magnetohydrostatic filaments. The filament exhibits a polarization pattern in which the magnetic field is apparently perpendicular to the filament when observed from the direction perpendicular to the magnetic field. When the line-of-sight is parallel to the global magnetic field, the observed polarization pattern is dependent on the center-to-surface density ratio for the filament and the concentration of the gas mass toward the central magnetic flux tube. Filaments with low center-to-surface density ratios have an insignificant degree of polarization when observed from the direction parallel to the global magnetic field. However, models with a large center-to-surface density ratio have polarization patterns that indicate the filament is perpendicularly threaded by the magnetic field. When mass is heavily concentrated at the central magnetic flux tube, which can be realized by the ambipolar diffusion process, the polarization pattern is similar to that expected for a low center-to-surface density contrast.
In this paper, we investigate the influences of the radiation pressures of the central engines on the black hole virial masses $M_{\rm{RM}}$ for 40 active galactic nuclei (AGNs) with high accretion rates. The radiation pressures on the clouds in broad-line regions are equivalent to decrease the gravitational forces of the central black holes, i.e., the clouds undergo the decreased gravitational forces of the black holes due to the central radiation. For the AGNs with the high accretion rates, a part of the gravitational forces of the black holes is counteracted by the central radiation, and the counteracted fraction depends on the percent of ionized hydrogen in clouds (or the ionized depth ratio of clouds). The black hole masses counteracted by the radiation pressures $M_{\rm{RP}}$ are not negligible compared to, or are comparable to $M_{\rm{RM}}$ at least for a part of the AGNs. The black hole masses $M_{\rm{\bullet}}$ are underestimated at least by a factor of 30--50 percent for the AGNs with the close- and super-Eddington limit, regardless of redshifts of sources. It is more appropriate to ignore the radiation pressures of the central sources for the AGNs with lower accretion rates, but the influences of the radiation pressures shall be taken into account to estimate the black hole masses of the extremely high accretion rate AGNs based on the reverberation mapping method.
Circular ribbon flares are usually related to spine-fan type magnetic topology containing null-points. In this paper, we investigate an X-class circular ribbon flare on 2012 October 23, using the multi-wavelength data from the \textit{Solar Dynamics Observatory}, \textit{Hinode}, and the \textit{Ramaty High Energy Solar Spectroscopic Imager}. In \ion{Ca}{2} H emission, the flare showed three ribbons with two highly elongated ones inside and outside a quasi-circular one, respectively. A hot channel was displayed in the extreme ultraviolet (EUV) emissions that infers the existence of a magnetic flux rope. Two hard X-ray (HXR) sources in the 12--25 keV energy band were located at the footpoints of this hot channel. Using a nonlinear force-free magnetic field extrapolation, we identify three topological structures: (1) a 3D null-point, (2) a flux rope below the fan of the null-point, and (3) a large-scale quasi-separatrix layers (QSL) induced by the quadrupolar-like magnetic field of the active region. We find that the null-point is embedded within the large-scale QSL. In our case, all three identified topological structures must be considered to explain all the emission features associated with the observed flare. Besides, the HXR sources are regarded as the consequence of the reconnection within or near the border of the flux rope.
The $m$-$z$ relation for type Ia supernovae is one of the key pieces of evidence supporting the cosmological `concordance model' with $\lambda_0 \approx 0.7$ and $\Omega_0 \approx 0.3$. However, it is well known that the $m$-$z$ relation depends not only on $\lambda_0$ and $\Omega_0$ (with $H_0$ as a scale factor) but also on the density of matter along the line of sight, which is not necessarily the same as the large-scale density. I investigate to what extent the measurement of $\lambda_0$ and $\Omega_0$ depends on this density when it is characterized by the parameter $\eta$ ($0 \le \eta \le 1$), which describes the ratio of density along the line of sight to the overall density. I also discuss what constraints can be placed on $\eta$, both with and without constraints on $\lambda_0$ and $\Omega_0$ in addition to those from the $m$-$z$ relation for type~Ia supernovae.
All atmosphere-less planetary bodies are covered with a dust layer, the so-called regolith, which determines the optical, mechanical and thermal properties of their surface. These properties depend on the regolith material, the size distribution of the particles it consists of, and the porosity to which these particles are packed. We performed experiments in parabolic flights to determine the gravity dependency of the packing density of regolith for solid-particle sizes of 60 $\mu$m and 1 mm as well as for 100-250 $\mu$m-sized agglomerates of 1.5 $\mu$m-sized solid grains. We utilized g-levels between 0.7 m s$^{-2}$ and 18 m s$^{-2}$ and completed our measurements with experiments under normal gravity conditions. Based on previous experimental and theoretical literature and supported by our new experiments, we developed an analytical model to calculate the regolith stratification of celestial rocky and icy bodies and estimated the mechanical yields of the regolith under the weight of an astronaut and a spacecraft resting on these objects.
Modified Newtonian dynamics (MOND) provides a paradigm alternative to dark matter that has been successful in fitting and predicting the rich phenomenology of rotating disc galaxies. There have also been attempts to test MOND in dispersion-supported early-type galaxies, but it remains unclear whether MOND can fit the various empirical properties of early-type galaxies. As a way of rigorously testing MOND in elliptical galaxies we calculate the MOND-predicted velocity dispersion profiles (VDPs) in the inner regions of $\sim 2000$ nearly round SDSS elliptical galaxies under a variety of assumptions on VD anisotropy, and then compare the predicted distribution of VDP slopes with the observed distribution in 11 ATLAS3d galaxies selected with essentially the same criteria. We find that the MOND model parameterised with an interpolating function that works well for rotating galaxies can also reproduce the observed distribution of VDP slopes based only on the observed stellar mass distribution without DM or any other galaxy-to-galaxy varying factor. This is remarkable in view that Newtonian dynamics with DM requires a specific amount and/or profile of DM for each galaxy in order to reproduce the observed distribution of VDP slopes. When we analyse non-round galaxy samples using the MOND-based spherical Jeans equation, we do not find any systematic difference in the mean property of the VDP slope distribution compared with the nearly round sample. However, in line with previous studies of MOND through individual analyses of elliptical galaxies, varying MOND interpolating function or VD anisotropy can lead to systematic change in the VDP slope distribution, indicating that a statistical analysis of VDPs can be used to constrain specific MOND models with an accurate measurement of VDP slopes or a prior constraint on VD anisotropy.
This document describes the Multi-Order Coverage map method (MOC) to specify arbitrary sky regions. The goal is to be able to provide a very fast comparison mechanism between coverage maps. The mechanism is based on the HEALPix sky tessellation algorithm. It is essentially a simple way to map regions of the sky into hierarchically grouped predefined cells.
We study the effects of planetary late migration on the gas giants obliquities. We consider the planetary instability models from Nesvorny & Morbidelli (2012), in which the obliquities of Jupiter and Saturn can be excited when the spin-orbit resonances occur. The most notable resonances occur when the $s_7$ and $s_8$ frequencies, changing as a result of planetary migration, become commensurate with the precession frequencies of Jupiter's and Saturn's spin vectors. We show that Jupiter may have obtained its present obliquity by crossing of the $s_8$ resonance. This would set strict constrains on the character of migration during the early stage. Additional effects on Jupiter's obliquity are expected during the last gasp of migration when the $s_7$ resonance was approached. The magnitude of these effects depends on the precise value of the Jupiter's precession constant. Saturn's large obliquity was likely excited by capture into the $s_8$ resonance. This probably happened during the late stage of planetary migration when the evolution of the $s_8$ frequency was very slow, and the conditions for capture into the spin-orbit resonance with $s_8$ were satisfied. However, whether or not Saturn is in the spin-orbit resonance with $s_8$ at the present time is not clear, because the existing observations of Saturn's spin precession and internal structure models have significant uncertainties.
Asteroids and meteorites provide key evidence on the formation of planetesimals in the Solar System. Asteroids are traditionally thought to form in a bottom-up process by coagulation within a population of initially km-scale planetesimals. However, new models challenge this idea by demonstrating that asteroids of sizes from 100 to 1000 km can form directly from the gravitational collapse of small particles which have organised themselves in dense filaments and clusters in the turbulent gas. Particles concentrate passively between eddies down to the smallest scales of the turbulent gas flow and inside large-scale pressure bumps and vortices. The streaming instability causes particles to take an active role in the concentration, by piling up in dense filaments whose friction on the gas reduces the radial drift compared to that of isolated particles. In this chapter we review new paradigms for asteroid formation and compare critically against the observed properties of asteroids as well as constraints from meteorites. Chondrules of typical sizes from 0.1 to 1 mm are ubiquitous in primitive meteorites and likely represent the primary building blocks of asteroids. Chondrule-sized particles are nevertheless tightly coupled to the gas via friction and are therefore hard to concentrate in large amounts in the turbulent gas. We review recent progress on understanding the incorporation of chondrules into the asteroids, including layered accretion models where chondrules are accreted onto asteroids over millions of years. We highlight in the end ten unsolved questions in asteroid formation where we expect that progress will be made over the next decade.
Terrestrial gamma-ray flashes (TGFs) are short intense flashes of gamma rays associated with lightning activity in thunderstorms. Using Monte Carlo simulations of the relativistic runaway electron avalanche (RREA) process, theoretical predictions for the temporal and spectral evolution of TGFs are compared to observations made with the Gamma-ray Burst Monitor (GBM) on board the Fermi Gamma-ray Space Telescope. Assuming a single source altitude of 15 km, a comparison of simulations to data is performed for a range of empirically chosen source electron variation time scales. The data exhibit a clear softening with increased source distance, in qualitative agreement with theoretical predictions. The simulated spectra follow this trend in the data, but tend to underestimate the observed hardness. Such a discrepancy may imply that the basic RREA model is not sufficient. Alternatively, a TGF beam that is tilted with respect to the zenith could produce an evolution with source distance that is compatible with the data. Based on these results, we propose that the source electron distributions of TGFs observed by GBM vary on time scales of at least tens of microseconds, with an upper limit of approx. 100 microseconds.
The short duration of mass ejection events during the Ascending Giant Branch (AGB) phase of stellar evolution can take different forms ranging thermal pulses to regular stellar pulsations. In most cases the shells ejected by the star are assumed to reach the escape velocity and expand into circumstellar space until they dissipate or merge with the surrounding gas. In this paper we investigate the case of an AGB star that emits a pulse of material below the escape velocity as may occur during a Common Envelope event. We explore the evolution of the shell created by this short mass loss event. We seek to determine when the shell falls back onto the star as opposed to being driven to escape velocity by the action of winds which occur after shell ejection. The problem is solved via 2.5D AMR AstroBEAR hydrodynamic simulations and a simplified one dimensional analytic model. We find that for given set of initial wind characteristics there is a critical shell velocity that distinguishes between shell fallback and shell escape. We provide an estimate of this critical velocity based on an analytic model. We use AstroBEAR to simulate models in which the shell material is ejected with a mild density perturbation. We find that the perturbed shell evolves into clumpy structures. We discuss the relevance of our results for both single and binary AGB stars and for binary stars our results provide a step towards understanding how fallback material may contribute to forming the substantial population of observed of post-AGB with dusty disks.
The Gamma-ray Burst Monitor (GBM) on board the Fermi Gamma-ray Space Telescope has triggered on over 300 terrestrial gamma-ray flashes (TGFs) since its launch in June 2008. With 14 detectors, GBM collects on average ~100 counts per triggered TGF, enabling unprecedented studies of the time profiles of TGFs. Here we present the first rigorous analysis of the temporal properties of a large sample of TGFs (278), including the distributions of the rise and fall times of the individual pulses and their durations. A variety of time profiles are observed with 19 of TGFs having multiple pulses separated in time and 31 clear cases of partially overlapping pulses. The effect of instrumental dead time and pulse pileup on the temporal properties are also presented. As the observed gamma ray pulse structure is representative of the electron flux at the source, TGF pulse parameters are critical to distinguish between relativistic feedback discharge and lightning leader models. We show that at least 67% of TGFs at satellite altitudes are significantly asymmetric. For the asymmetric pulses, the rise times are almost always shorter than the fall times. Those which are not are consistent with statistical fluctuations. The median rise time for asymmetric pulses is ~3 times shorter than for symmetric pulses while their fall times are comparable. The asymmetric shapes observed are consistent with the relativistic feedback discharge model when Compton scattering of photons between the source and Fermi is included, and instrumental effects are taken into account.
We investigate the effects of homogeneous general dark energy on the non-linear matter perturbation in fully general relativistic context. Taking into account the next-to-leading corrections, the total power spectrum with general dark energy deviates from the LambdaCDM spectrum, which is nearly identical to that in the Einstein-de Sitter universe, as large as a few percent at scales comparable to that for the baryon acoustic oscillations and increases on smaller scales. The contribution from the curvature perturbation, which is absent in the Newtonian theory, exhibits even more drastic difference larger than 100%, while the amplitude is heavily suppressed on all scales.
In this work, a gradual solar energetic particle (SEP) event observed by multi-spacecraft has been investigated with simulations. Based on a numerical solution of the Fokker-Planck focused transport equation, we obtain the intensity time profiles of SEPs accelerated by an interplanetary shock in the three-dimensional interplanetary space. The shock is treated as a moving source of energetic particles with a distribution function. By fitting the 1979/03/01 SEP event observed by $Helios$ 1, $Helios$ 2, and $IMP$ 8 with our simulations simultaneously, we obtain the best parameters for the shock acceleration strength model. And we also find that the particle perpendicular diffusion coefficient with the level of $\sim 1\%-3\%$ of parallel diffusion coefficient at $1$ AU should be included. In addition, the gradient of SEP fluxes in the decay phase is more sensitive to the shock acceleration strength parameters than that is to the perpendicular diffusion coefficient.
NGC 4203 is a nearby early-type galaxy surrounded by a very large, low-column-density HI disc. In this paper we study the star formation efficiency in the gas disc of NGC 4203 by using the UV, deep optical imaging and infrared data. We confirm that the HI disc consists of two distinct components: an inner star forming ring with radius from $\sim$ 1 to $\sim$ 3 R$_{eff}$, and an outer disc. The outer HI disc is 9 times more massive than the inner HI ring. At the location of the inner HI ring we detect spiral-like structure both in the deep $g'-r'$ image and in the 8 $\mu$m $Spitzer$-IRAC image, extending in radius up to $\sim$ 3 R$_{eff}$. These two gas components have a different star formation efficiency likely due to the different metallicity and dust content. The inner component has a star formation efficiency very similar to the inner regions of late-type galaxies. Although the outer component has a very low star formation efficiency, it is similar to that of the outer regions of spiral galaxies and dwarfs. We suggest that these differences can be explained with different gas origins for the two components such as stellar mass loss for the inner HI ring and accretion from the inter galactic medium (IGM) for the outer HI disc. The low level star formation efficiency in the outer HI disc is not enough to change the morphology of NGC 4203, making the depletion time of the HI gas much too long.
We investigate the early-time light-curves of a large sample of 223 type II supernovae (SNe) from the Sloan Digital Sky Survey and the Supernova Legacy Survey. Having a cadence of a few days and sufficient non-detections prior to explosion, we constrain rise-times, i.e. the durations from estimated first to maximum light, as a function of effective wavelength. At restframe g-band (4722A), we find a distribution of fast rise-times with median of (7.5+/-0.3) days. Comparing these durations with analytical shock models of Rabinak and Waxman (2013); Nakar and Sari (2010) and hydrodynamical models of Tominaga et al. (2009), which are mostly sensitive to progenitor radius at these epochs, we find a median characteristic radius of less than 400 solar radii. The inferred radii are on average much smaller than the radii obtained for observed red supergiants (RSG). Investigating the post-maximum slopes as a function of effective wavelength in the light of theoretical models, we find that massive hydrogen envelopes are still needed to explain the plateaus of SNe II. We therefore argue that the SN II rise-times we observe are either a) the shock cooling resulting from the core collapse of RSG with small and dense envelopes, or b) the delayed and prolonged shock breakout of the collapse of a RSG with an extended atmosphere or embedded within pre-SN circumstellar material.
In a recent study star-forming galaxies with HeII emission at moderate redshifts have been found to occur in two modes, distinguished by the width of their HeII emission lines. Broad HeII emission has been attributed to stellar emission from a population of evolved Wolf-Rayet (WR) stars while narrow HeII emission has been attributed to nebular emission excited by a population of very hot PopIII stars formed in pockets of pristine gas at moderate redshifts. In this work we propose an alternative scenario for the origin of the narrow HeII emission, namely very massive stars (VMS) at low metallicity (Z) which form strong but slow WR-type stellar winds due to their proximity to the Eddington limit. We estimate the expected HeII line fluxes and equivalent widths based on wind models for VMS and population synthesis models, and compare the results with recent observations of star-forming galaxies at moderate redshifts. The observed HeII line strengths and equivalent widths are in line with what is expected for a population of VMS in one or more young super-clusters located within these galaxies. In our scenario the two observed modes of HeII emission originate from massive stellar populations in distinct evolutionary stages at low Z. If this interpretation is correct there is no need to postulate the existence of PopIII stars at moderate redshifts to explain the observed narrow HeII emission. An interesting possibility is the existence of self-enriched VMS with similar WR-type spectra at extremely low Z. Stellar HeII emission from such very early generations of VMS may be detectable in future studies of star-forming galaxies at high redshifts, as they are planned for the James Webb Space Telescope. The fact that the HeII emission of VMS is largely neglected in current population synthesis models will generally affect the interpretation of the integrated spectra of young stellar populations.
Shift-and-add is an approach employed to mitigate the phenomenon of resolution degradation in images acquired through a turbulent medium. Using this technique, a large number of consecutive short exposures is registered below the coherence time of the atmosphere or other blurring medium. The acquired images are shifted to the position of the brightest speckle and stacked together to obtain high-resolution and high signal-to-noise frame. In this paper we present a highly efficient method for determination of frames shifts, even if in a single frame the object cannot be distinguished from the background noise. The technique utilizes our custom genetic algorithm, which iteratively evolves a set of image shifts. We used the maximal energy of stacked images as an objective function for shifts estimation and validate the efficiency of the method on simulated and real images of simple and complex sources. Obtained results confirmed, that our proposed method allows for the recovery of spatial distribution of objects even only 2% brighter than their background. The presented approach extends significantly current limits of image reconstruction with the use of shift-and-add method. The applications of our algorithm include both the optical and the infrared imaging. Our method may be also employed as a digital image stabilizer in extremely low light level conditions in professional and consumer applications.
In order to explain light curve (LC) for Supernova (SN) we derive a classical formula for the conversion of the flux of kinetic energy into radiation. We then introduce a correction for the absorption adopting an optical depth as function of the time. The developed framework allows to fit the LC of type Ia SN 2005cf ( B and V ) and type IIp SN 2004A (B,V,I and R ). A relativistic formula for the flux of kinetic energy is also derived in terms of a Taylor expansion and the application is done to the LC of GRB 050814. The decay of the radioactive isotopes as a driver the LC for SNs is also reviewed and a new formulation is introduced. The Arnett's formula for bolometric luminosity is corrected for the optical depth and applied to SN 2001ay.
Abell 1767 is a dynamically relaxed, cD cluster of galaxies with a redshift of 0.0703. Among 250 spectroscopically confirmed member galaxies within a projected radius of 2.5r_{200}, 243 galaxies (~ 97%) are spectroscopically covered by the Sloan Digital Sky Survey (SDSS). Based on this homogeneous spectral sample, the stellar evolutionary synthesis code, STARLIGHT, is applied to investigate the stellar populations and star formation histories (SFHs) of cluster galaxies. The star formation properties of galaxies, such as mean stellar ages, metallicities, stellar masses, and star formation rates (SFRs), are presented as the functions of local galaxy density. Strong environmental effect is found in the manner that massive galaxies in the high-density core region of cluster tend to have higher metallicities, longer mean stellar ages, and lower specific star formation rates (SSFRs), and their recent star formation activities have been remarkably suppressed. In addition, the correlations of the metallicity and SSFR with stellar mass are confirmed.
We present an analysis of Kitt Peak National Observatory and Lowell Observatory observations of comet 103P/Hartley 2 obtained from August through December 2010. The results are then compared with contemporaneous observations made by the EPOXI spacecraft. Each ground-based dataset has previously been investigated individually; the combined dataset has complementary coverage that reduces the time between observing runs and allows us to determine additional apparent periods at intermediate times. We compare CN coma morphology between ground-based datasets, making nine new measurements of apparent periods. The first five are consistent with the roughly linearly increasing apparent period during the apparition found by previous authors. The final four suggest that the change in apparent period slowed or stopped by late November. We also measure an inner coma lightcurve in both CN and R-band ground-based images, finding a single-peaked lightcurve which repeats in phase with the coma morphology. The apparent period from the lightcurve had significantly larger uncertainties than from the coma morphology, but varied over the apparition in a similar manner. Our ground-based lightcurve aligns with the published EPOXI lightcurve, indicating that the lightcurve represents changing activity rather than viewing geometry of structures in the coma. The EPOXI lightcurve can best be phased by a triple-peaked period near 54-55 hr that increases from October to November. This phasing reveals that the spacing between maxima is not constant, and that the overall lightcurve shape evolves from one triple-peaked cycle to the next. These behaviors suggest that much of the scatter in apparent periods derived from ground-based datasets acquired at similar epochs are likely due to limited sampling of the data.
Asymmetric X-ray emission and powerful cluster-scale radio halo indicate that A2319 is a merging cluster of galaxies. This paper presents our multicolor photometry for A2319 with 15 optical intermediate filters in the Beijing-Arizona-Taiwan-Connecticut (BATC) system. There are 142 galaxies with known spectroscopic redshifts within the viewing field, including 128 member galaxies (called sample I).A large velocity dispersion in the rest frame suggests a merger dynamics in A2319. The contour map of projected density and localized velocity structure confirm the so-called A2319B substructure, at ~ 10' NW to the main concentration A2319A. The spectral energy distributions (SEDs) of more than 30,000 sources are obtained in our BATC photometry down to V ~ 20 mag. With color-color diagrams and photometric redshift technique, 233 galaxies brighter than h=19.0 are newly selected as member candidates. The early-type galaxies are found to follow a tight color-magnitude correlation. Based on sample I and the enlarged sample of member galaxies (called sample II), subcluster A2319B is confirmed. A strong environmental effect on star formation histories is found in the manner that galaxies in the sparse regions have various star formation histories, while galaxies in the dense regions are found to have shorter SFR time scales, older stellar ages, and higher ISM metallicities. For the merging cluster A2319, local surface density is a better environmental indicator rather than the clustercentric distance. Compared with the well-relaxed cluster A2589, a higher fraction of star-forming galaxies is found in A2319, indicating that the galaxy-scale turbulence stimulated by the subcluster merger might have played a role in triggering the star formation activity.
We investigate the existence of magnetohydrostatic equilibria for topologically complex magnetic fields. The approach employed is to perform ideal numerical relaxation experiments. We use a newly-developed Lagrangian relaxation scheme that exactly preserves the magnetic field topology during the relaxation. Our configurations include both twisted and sheared fields, of which some fall into the category for which Parker (1972) predicted no force-free equilibrium. The first class of field considered contains no magnetic null points, and field lines connect between two perfectly conducting plates. In these cases we observe only resolved current layers of finite thickness. In further numerical experiments we confirm that magnetic null points are loci of singular currents.
PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.
A model of cyclical (sinusoidal) motion of the solar system, intercepting event lines distributed at fixed intervals, explains the pattern of timings of mass extinctions, earlier glaciations, largest impact craters and the largest known extrusions of magma in the history of the Earth. The model reveals links between several sets of key events, including the end-Cretaceous and end-Ordovician extinctions with the Marinoan glaciation, and the end-Permian with the end-Serpukhovian extinctions. The model is supported by significant clusters of events and a significant reduction of impact crater size with position (sine value). The pattern of event lines is sustained to the earliest-dated impact craters (2023 and 1849 Ma) and to the origin of the solar system, close to 4567.4 Ma. The implication is that, for the entirety of its existence, the solar system has passed in a consistent manner through a predictably structured galaxy. Dark matter is a possible contender for the structure determining the event lines.
We compare estimates of the speed and width of coronal mass ejections (CMEs) in several catalogs for the CMEs associated with ~200 solar energetic particle (SEP) events in 2006-2013 that included 25 MeV protons. The catalogs used are: CDAW, CACTUS, SEEDS and CORIMP, all derived from observations by the LASCO coronagraphs on the SOHO spacecraft, the CACTUS catalog derived from the COR2 coronagraphs on the STEREO-A and -B spacecraft, and the DONKI catalog, which uses observations from SOHO and the STEREO spacecraft. We illustrate how, for this set of events, CME parameters can differ considerably in each catalog. The well-known correlation between CME speed and proton event intensity is shown to be similar for most catalogs, but this is largely because it is determined by a few large particle events associated with fast CMEs, and small events associated with slow CMEs. Intermediate particle events "shuffle" in position when speeds from different catalogs are used. Quadrature spacecraft CME speeds do not improve the correlation. CME widths also vary widely between catalogs, and are influenced by plane of the sky projection and how the width is inferred from the coronagraph images. The high degree of association (~50%) between the 25 MeV proton events and "full halo" (360 deg.-width) CMEs as defined in the CDAW catalog is removed when other catalogs are considered. Using CME parameters from the quadrature spacecraft, the SEP intensity is correlated with CME width, which is also correlated with CME speed.
Polarimetry is one of the keys to enhanced direct imaging of exoplanets. Not only does it deliver a differential observable providing extra contrast, but when coupled with spectroscopy, it also reveals valuable information on the exoplanetary atmospheric composition. Nevertheless, angular separation and contrast ratio to the host-star make for extremely challenging observation. Producing detailed predictions for exactly how the expected signals should appear is of critical importance for the designs and observational strategies of tomorrow's telescopes. We aim at accurately determining the magnitudes and evolution of the main observational signatures for imaging an exoplanet: separation, contrast ratio to the host-star and polarization as a function of the orbital geometry and the reflectance parameters of the exoplanet. These parameters were used to construct polarized-reflectance model based on the input of orbital parameters and two albedo values. The model is able to calculate a variety of observational predictions for exoplanets at any orbital time. The inter-dependency of the three main observational criteria -separation, contrast ratio, polarization- result in a complex time-evolution of the system. They greatly affect the viability of planet observation by direct imaging. We introduce a new generic display of the main observational criteria, which enables an observer to determine whether an exoplanet is within detection limits: the Separation-POlarization-Contrast diagrams (SPOC). We explore the complex effect of orbital and albedo parameters on the visibility of an exoplanet. The code we developed is available for public use and collaborative improvement on the python package index, together with its documentation. It is another step towards a full comprehensive simulation tool for predicting and interpreting the results of future observational exoplanetary discovery campaigns.
We present the first large-scale study of the photometric and structural relations followed by early-type galaxies (ETGs) in the Antlia cluster. Antlia is the third nearest populous galaxy cluster after Fornax and Virgo (d $\sim 35$ Mpc). A photographic catalog of its galaxy content was built by Ferguson & Sandage in 1990 (FS90). Afterwards, we performed further analysis of the ETG population located at the cluster centre. Now, we extend our study covering an area four times larger, calculating new total magnitudes and colours, instead of isophotal photometry, as well as structural parameters obtained through S\'ersic model fits extrapolated to infinity. Our present work involves a total of 177 ETGs, out of them 56 per cent have been cataloged by FS90 while the rest (77 galaxies) are newly discovered ones. Medium-resolution GEMINI and VLT spectra are used to confirm membership when available. Including radial velocities from the literature, 59 ETGs are confirmed as Antlia members. Antlia scaling relations mainly support the existence of unique functions (linear and curved) that join bright and dwarf ETGs, excluding compact ellipticals (cEs). Lenticular galaxies are outliers only with respect to the curved relation derived for effective surface brightness versus absolute magnitude. The small number of bright ellipticals and cEs present in Antlia, prevents us from testing if the same data can be fitted with two different linear sequences, for bright and dwarf ETGs. However, adding data from other clusters and groups, the existence of such sequences is also noticeable in the same scaling relations.
We compare the shapes and intrinsic alignments of galaxies in the MassiveBlack-II cosmological hydrodynamic simulation (MBII) to those in a dark matter-only (DMO) simulation performed with the same volume (100$h^{-1}$Mpc)$^{3}$, cosmological parameters, and initial conditions. Understanding the impact of baryonic physics on galaxy shapes and alignments and their relation to the dark matter distribution should prove useful to map the intrinsic alignments of galaxies from hydrodynamic to dark matter-only simulations. We find that dark matter subhalos are typically rounder in MBII, and the shapes of stellar matter in low mass galaxies are more misaligned with the shapes of the dark matter of the corresponding subhalos in the DMO simulation. At $z=0.06$, the fractional difference in the mean misalignment angle between MBII and DMO simulations varies from $\sim 28 \% - 12 \%$ in the mass range $10^{10.8} - 6.0 \times 10^{14} h^{-1}M_{\odot}$. We study the dark matter halo shapes and alignments as a function of radius, and find that while galaxies in MBII are more aligned with the inner parts of their dark matter subhalos, there is no radial trend in their alignments with the corresponding subhalo in the DMO simulation. This result highlights the importance of baryonic physics in determining the alignment of the galaxy with respect to the inner parts of the halo. Finally, we compare the ellipticity-direction (ED) correlation for galaxies to that for dark matter halos, finding that it is suppressed on all scales by stellar-dark matter misalignment. In the projected shape-density correlation ($w_{\delta+}$), which includes ellipticity weighting, this effect is partially canceled by the higher mean ellipticities of the stellar component, but differences of order $30-40\%$ remain on scales $> 1$ Mpc over a range of subhalo masses, with scale-dependent effects below $1$ Mpc.
With the goal of deriving dissipative hydrodynamics from an action, we study classical actions for open systems, which follow from the generic structure of effective actions in the Schwinger-Keldysh Closed-Time-Path formalism with two time axes and a doubling of degrees of freedom. The central structural feature of such effective actions is the coupling between degrees of freedom on the two time axes. This reflects the fact that from an effective field theory point of view, dissipation is the loss of energy of the low-energy hydrodynamical degrees of freedom to the integrated-out, UV degrees of freedom of the environment. The dynamics of only the hydrodynamical modes may therefore not posses a conserved stress-energy tensor. After a general discussion of the CTP effective actions, we use the variational principle to derive the energy-momentum balance equation for a dissipative fluid from an effective Goldstone action of the long-range hydrodynamical modes. Despite the absence of conserved energy and momentum, we show that we can construct the first-order dissipative stress-energy tensor and derive the Navier-Stokes equations near hydrodynamical equilibrium. The shear viscosity is shown to vanish in the classical theory under consideration, while the bulk viscosity is determined by the form of the effective action. We also discuss the thermodynamics of the system and analyse the entropy production.
The description of nuclear reactions induced by supernova neutrinos has
witnessed significant progress during the recent years. At the energies and
momentum transfers relevant for supernova neutrinos neutrino-nucleus cross
sections are dominated by allowed transitions, however, often with
non-negligible contributions from (first) forbidden transitions. For several
nuclei allowed Gamow-Teller strength distributions could be derived from
charge-exchange reactions and from inelastic electron scattering data.
Importantly the diagonalization shell model has been proven to accurately
describe these data and hence became the appropriate tool to calculate the
allowed contributions to neutrino-nucleus cross sections for supernova
neutrinos. Higher multipole contributions are usually calculated within the
framework of the Quasiparticle Random Phase Approximation, which describes the
total strength and the position of the giant resonances quite well.
This manuscript reviews the recent progress achieved in calculating
supernova-relevant neutrino-nucleus cross sections and discusses its
verification by data. Moreover, the review summarizes also the impact which
neutrino-nucleus reactions have on the dynamics of supernovae and on the
associated nucleosynthesis. These include the absorption of neutrinos by nuclei
(the inverse of nuclear electron capture which is the dominating
weak-interaction process during collapse), inelastic neutrino-nucleus
scattering and nuclear de-excitation by neutrino-pair emission. We also discuss
the role of neutrino-induced reactions for the recently discovered $\nu p$
process, for the r-process and for the neutrino process, for which
neutrino-nucleus reactions have the largest impact. Finally, we briefly review
neutrino-nucleus reactions important for the observation of supernova neutrinos
by earthbound detectors. (Abridged)
Cosmology is making impressive progress and it is producing stringent bounds on the sum of the neutrino masses {\Sigma}, a parameter of great importance for the current laboratory experiments. In this letter, we exploit the potential relevance of the analysis of Palanque-Delabrouille et al. [JCAP 1502, 045 (2015)] to the neutrinoless double beta decay (0\nu\beta\beta) search. In fact, this result seems to favor the normal hierarchy spectrum for the light neutrino masses. If this is confirmed, the impact on the 0\nu\beta{\beta} experiments will be tremendous, bringing the possibility of detecting a signal out of the reach of the next generation of experiments.
It has been suggested that several small-scale structure anomalies in $\Lambda$CDM cosmology can be solved by strong self-interaction between dark matter particles. It was shown by Braaten and Hammer that the presence of a near threshold S-wave resonance can make the scattering cross section at nonrelativistic speeds come close to saturating the unitarity bound. This can result in the formation of a stable bound state of two asymmetric dark matter particles (which we call darkonium). Laha and Braaten studied the nuclear recoil energy spectrum in dark matter direct detection experiments due to this incident bound state. Here we study the angular recoil spectrum, and show that it is uniquely determined up to normalization by the S-wave scattering length. Observing this angular recoil spectrum in a dark matter directional detection experiment will uniquely determine many of the low-energy properties of dark matter independent of the underlying dark matter microphysics.
Sterile neutrinos with mass ~1 eV and order 10% mixing with active neutrinos have been proposed as a solution to anomalies in neutrino oscillation data, but are tightly constrained by cosmological limits. It was recently shown that these constraints are avoided if sterile neutrinos couple to a new MeV-scale gauge boson A'. However, even this scenario is restricted by structure formation constraints when A'-mediated collisional processes lead to efficient active-to-sterile neutrino conversion after neutrinos have decoupled. In view of this, we reevaluate in this paper the viability of sterile neutrinos with such "secret" interactions. We carefully dissect their evolution in the early Universe, including the various production channels and the expected modifications to large scale structure formation. We argue that there are two regions in parameter space - one at very small A' coupling, one at relatively large A' coupling - where all constraints from big bang nucleosynthesis (BBN), cosmic microwave background (CMB), and large scale structure (LSS) data are satisfied. Interestingly, the large A' coupling region is precisely the region that was previously shown to have potentially important consequences for the small scale structure of dark matter halos if the A' boson couples also to the dark matter in the Universe.
The formalism for describing a metric and the corresponding scalar in terms of multipole moments has recently been developed for scalar-tensor theories. We take advantage of this formalism in order to obtain expressions for the observables that characterise geodesics in terms of the moments. These expressions provide some insight into how the structure of a scalarized compact object affects observables. They can also be used to understand how deviations from general relativity are imprinted on the observables.
We describe a method for removing the effect of confounders in order to reconstruct a latent quantity of interest. The method, referred to as half-sibling regression, is inspired by recent work in causal inference using additive noise models. We provide a theoretical justification and illustrate the potential of the method in a challenging astronomy application.
In the scenario of single-field inflation, this field is done in terms of Jacobian elliptic functions. This approach provides, when constrained to particular cases, analytic solutions already known in the past, generalizing them to a bigger family of analytical solutions. The emergent cosmology is analysed using the Hamilton-Jacobi approach and then, the main results are contrasted with the recent measurements obtained from the Planck 2015 data.
The present paper is devoted to find a new generalization of the generalized Chaplygin gas. Therefore, starting from the Hubble parameter associated to the Chaplygin scalar field and using some elliptic identities, the elliptic generalization is straightforward. Thus, all relevant quantities that drive inflation are calculated exactly. Finally, using the measurement on inflation from the Planck 2015 results, observational constraints on the parameters are given.
Beginning with a set of simplified models for spin-0, spin-$\half$, and spin-1 dark matter candidates using completely general Lorentz invariant and renormalizable Lagrangians, we derive the full set of non-relativistic operators and nuclear matrix elements relevant for direct detection of dark matter, and use these to calculate rates and recoil spectra for scattering on various target nuclei. This allows us to explore what high energy physics constraints might be obtainable from direct detection experiments, what degeneracies exist, which operators are ubiquitous and which are unlikely or sub-dominant. We find that there are operators which are common to all spins as well operators which are unique to spin-$\half$ and spin-1 and elucidate two new operators which have not been previously considered. In addition we demonstrate how recoil energy spectra can distinguish fundamental microphysics if multiple target nuclei are used. Our work provides a complete roadmap for taking generic fundamental dark matter theories and calculating rates in direct detection experiments. This provides a useful guide for experimentalists designing experiments and theorists developing new dark matter models.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
Priorities in exo-planet research are rapidly moving from finding planets to characterizing their physical properties. Of key importance is their chemical composition, which feeds back into our understanding of planet formation. For the foreseeable future, far-ultraviolet spectroscopy of white dwarfs accreting planetary debris remains the only way to directly and accurately measure the bulk abundances of exo-planetary bodies. The exploitation of this method is limited by the sensitivity of HST, and significant progress will require a large-aperture space telescope with a high-throughput ultraviolet spectrograph.
Local galaxies are broadly divided into two main classes, star-forming (gas-rich) and quiescent (passive and gas-poor). The primary mechanism responsible for quenching star formation in galaxies and transforming them into quiescent and passive systems is still unclear. Sudden removal of gas through outflows or stripping is one of the mechanisms often proposed. An alternative mechanism is so-called "strangulation", in which the supply of cold gas to the galaxy is halted. Here we report that the difference between quiescent and star forming galaxies in terms of stellar metallicity (i.e. the fraction of metals heavier than helium in stellar atmospheres) can be used to discriminate efficiently between the two mechanisms. The analysis of the stellar metallicity in local galaxies, from 26,000 spectra, clearly reveals that strangulation is the primary mechanism responsible for quenching star formation, with a typical timescale of 4 billion years, at least for local galaxies with a stellar mass less than 10^11 solar masses. This result is further supported independently by the stellar age difference between quiescent and star-forming galaxies, which indicates that quiescent galaxies of less than 10^11 solar masses are on average observed four billion years after quenching due to strangulation.
Whilst the star formation rate (SFR) of molecular clouds and galaxies is key in understanding galaxy evolution, the physical processes which determine the SFR remain unclear. This uncertainty about the underlying physics has resulted in various different star formation laws, all having substantial intrinsic scatter. Extending upon previous works that define the column density of star formation (Sigma_SFR) by the gas column density (Sigma_gas), we develop a new universal star formation (SF) law based on the multi-freefall prescription of gas. This new SF law relies predominantly on the probability density function (PDF) and on the sonic Mach number of the turbulence in the star-forming clouds. By doing so we derive a relation where the star formation rate (SFR) correlates with the molecular gas mass per multi-freefall time, whereas previous models had used the average, single-freefall time. We define a new quantity called maximum (multi-freefall) gas consumption rate (MGCR) and show that the actual SFR is only about 0.4% of this maximum possible SFR, confirming the observed low efficiency of star formation. We show that placing observations in this new framework (Sigma_SFR vs. MGCR) yields a significantly improved correlation with 3-4 times reduced scatter compared to previous SF laws and a goodness-of-fit parameter R^2=0.97. By inverting our new relationship, we provide sonic Mach number predictions for kpc-scale observations of Local Group galaxies as well as unresolved observations of local and high-redshift disk and starburst galaxies that do not have independent, reliable estimates for the turbulent cloud Mach number.
We present an implementation of the Schwarzschild orbit superposition method which can be used for constructing self-consistent equilibrium models of barred or non-barred disc galaxies, or of elliptical galaxies with figure rotation. This is a further development of the publicly available code SMILE; its main improvements include a new efficient representation of an arbitrary gravitational potential using two-dimensional spline interpolation of Fourier coefficients in the meridional plane, as well as the ability to deal with rotation of the density profile and with multicomponent mass models. We compare several published methods for constructing composite axisymmetric disc--bulge--halo models and demonstrate that our code produces the models that are closest to equilibrium. We also apply it to create models of triaxial elliptical galaxies with cuspy density profiles and figure rotation, and find that such models can be found and are stable over many dynamical times in a wide range of pattern speeds and angular momenta, covering both slow- and fast-rotator classes. We then attempt to create models of strongly barred disc galaxies, using an analytic three-component potential, and find that it is not possible to make a stable dynamically self-consistent model for this density profile. Finally, we take snapshots of two N-body simulations of barred disc galaxies embedded in nearly-spherical haloes, and construct equilibrium models using only information on the density profile of the snapshots. We demonstrate that such reconstructed models are in near-stationary state, in contrast with the original N-body simulations, one of which displayed significant secular evolution.
A fraction of the early-type galaxy population hosts a prominent dust lane. Methods to quantify the dust content of these systems based on optical imaging data usually yield dust masses which are an order of magnitude lower than dust masses derived from the observed FIR emission. High-quality optical data from the Next Generation Virgo cluster Survey (NGVS) and FIR/submm observations from the Herschel Virgo Cluster Survey (HeViCS) allow us to revisit previous methods to determine the dust content in galaxies and explore new ones. We aim to derive the dust mass in NGC 4370 from both optical and FIR data, and investigate the need to invoke a putative diffuse dust component. We create color and attenuation maps, which are converted to approximate dust mass maps based on simple dust geometries. Dust masses are also derived from SED fits to FIR/submm observations. Finally, inverse radiative transfer fitting is performed to investigate more complex dust geometries. The empirical methods applied to the optical data yield lower limits of 3.4e5 solar masses, an order of magnitude below the total dust masses derived from SED fitting. In contrast, radiative transfer models yield dust masses which are slightly lower, but fully consistent with the FIR-derived mass. Dust is more likely to be distributed in a ring around the centre of NGC 4370 as opposed to an exponential disc or a simple foreground screen. Moreover, using inverse radiative transfer fitting, we are able to constrain most of the parameters describing these geometries. The resulting dust masses are high enough to account for the dust observed at FIR/submm wavelengths, so that no diffuse dust component needs to be invoked. We furthermore caution for the interpretation of dust masses and optical depths based on optical data alone, using overly simplistic star-dust geometries and the neglect of scattering effects. [ABRIDGED]
We perform numerical simulations on the merger of multiple black holes (BHs) in primordial gas at early cosmic epochs. We consider two cases of BH mass: $M_{BH} = 30 M_{\odot}$ and $M_{BH} = 10^4 M_{\odot}$. Attention is concentrated on the effect of the dynamical friction by gas in a host object. The simulations incorporate such general relativistic effects as the pericentre shift and gravitational wave emission. As a result, we find that multiple BHs are able to merge into one BH within 100 Myr in a wide range of BH density. The merger mechanism is revealed to be categorized into three types: gas-drag-driven merger (type A), interplay-driven merger (type B), and three-body-driven merger (type C). We find the relation between the merger mechanism and the ratio of the gas mass within the initial BH orbit ($M_{gas}$) to the total BH mass (${\Sigma}M_{BH}$). Type A merger occurs if $M_{gas} \gtrsim 10^5 {\Sigma}M_{BH}$, type B if $M_{gas} \lesssim 10^5 {\Sigma}M_{BH}$, and type C if $M_{gas} \ll 10^5 {\Sigma}M_{BH}$. Supposing the gas and BH density based on the recent numerical simulations on first stars, all the BH remnants from first stars are likely to merge into one BH through the type B or C mechanism. Also, we find that multiple massive BHs ($M_{BH} = 10^4 M_{\odot}$) distributed over several parsec can merge into one BH through the type B mechanism, if the gas density is higher than $5\times 10^6$ cm$^{-3}$. The present results imply that the BH merger may contribute significantly to the formation of supermassive BHs at high redshift epochs.
Understanding the 3D flow topology around a planet embedded in its natal disk is crucial to the study of planet formation. 3D modifications to the well-studied 2D flow topology have the potential to resolve longstanding problems in both planet migration and accretion. We present a detailed analysis of the 3D isothermal flow field around a 5 Earth-mass planet on a fixed circular orbit, simulated using our high-resolution multi-GPU hydrodynamics code PEnGUIn. We show that, overall, the horseshoe region has a columnar structure extending vertically much beyond the Hill sphere of the planet. This columnar structure is only broken for some of the widest horseshoe streamlines, along which high altitude fluid descends and converges rapidly toward the planet, enters its Bondi sphere, performs one horseshoe turn, and exits radially in the midplane. A portion of this flow gathers enough speed to exit the horseshoe region altogether. We call this newly identified feature the "transient" horseshoe flow. As the flow continues close to the disk midplane, it splits into the up-down symmetric parts, and rolls up into a pair of counter-rotating, horizontal vortex lines shed downstream into the wake of the planet. This flow, unique to 3D, affects both planet migration and accretion. It prevents the planet from sustaining a hydrostatic atmosphere due to its intrusion into the Bondi sphere, and it also leads to a significant corotation torque on the planet, unanticipated by 2D analysis. In the reported simulation, starting with a $\Sigma\sim r^{-3/2}$ radial surface density profile, this torque is positive and largely cancels with the negative differential Lindblad torque, resulting in an almost 2 orders of magnitude reduction in the inferred rate of planet migration. Finally, we report that 3D effects can be suppressed by a sufficiently large disk viscosity, leading to results similar to those in 2D.
We use optical integral-field spectroscopic (IFS) data from 103 nearby galaxies at different stages of the merging event, from close pairs to merger remnants provided by the CALIFA survey, to study the impact of the interaction in the specific star formation and oxygen abundance on different galactic scales. To disentangle the effect of the interaction and merger from internal processes, we compared our results with a control sample of 80 non-interacting galaxies. We confirm the moderate enhancement (2-3 times) of specific star formation for interacting galaxies in central regions as reported by previous studies; however, the specific star formation is comparable when observed in extended regions. We find that control and interacting star-forming galaxies have similar oxygen abundances in their central regions, when normalized to their stellar masses. Oxygen abundances of these interacting galaxies seem to decrease compared to the control objects at the large aperture sizes measured in effective radius. Although the enhancement in central star formation and lower metallicities for interacting galaxies have been attributed to tidally induced inflows, our results suggest that other processes such as stellar feedback can contribute to the metal enrichment in interacting galaxies.
We explore a chameleon type of interacting dark matter-dark energy scenario in which a scalar field adiabatically traces the minimum of an effective potential sourced by the dark matter density. We discuss extensively the effect of this coupling on cosmological observables, especially the parameter degeneracies expected to arise between the model parameters and other cosmological parameters, and then test the model against observations of the cosmic microwave background (CMB) anisotropies and other cosmological probes. We find that the chameleon parameters $\alpha$ and $\beta$, which determine respectively the slope of the scalar field potential and the dark matter-dark energy coupling strength, can be constrained to $\alpha < 0.17$ and $\beta < 0.19$ using CMB data alone. The latter parameter in particular is constrained only by the late Integrated Sachs-Wolfe effect. Adding measurements of the local Hubble expansion rate $H_0$ tightens the bound on $\alpha$ by a factor of two, although this apparent improvement is arguably an artefact of the tension between the local measurement and the $H_0$ value inferred from Planck data in the minimal $\Lambda$CDM model. The same argument also precludes chameleon models from mimicking a dark radiation component, despite a passing similarity between the two scenarios in that they both delay the epoch of matter-radiation equality. Based on the derived parameter constraints, we discuss possible signatures of the model for ongoing and future large-scale structure surveys.
In Galactic halo stars, sulphur has been shown to behave like other $\alpha$-elements, but until now, no comprehensive studies have been done on this element in stars of other galaxies. Here, we use high-resolution ESO VLT/FLAMES/GIRAFFE spectra to determine sulphur abundances for 85 stars in the Sculptor dwarf spheroidal galaxy, covering the metallicity range $-2.5\leq \text{[Fe/H]} \leq-0.8$. The abundances are derived from the S~I triplet at 9213, 9228, and 9238~\AA. These lines have been shown to be sensitive to departure from local thermodynamic equilibrium, i.e. NLTE effects. Therefore, we present new NLTE corrections for a grid of stellar parameters covering those of the target stars. The NLTE-corrected sulphur abundances in Sculptor show the same behaviour as other $\alpha$-elements in that galaxy (such as Mg, Si, and Ca). At lower metallicities ($\text{[Fe/H]}\lesssim-2$) the abundances are consistent with a plateau at $\text{[S/Fe]}\approx+0.16$, similar to what is observed in the Galactic halo, $\text{[S/Fe]}\approx+0.2$. With increasing [Fe/H], the [S/Fe] ratio declines, reaching negative values at $\text{[Fe/H]}\gtrsim-1.5$. The sample also shows an increase in [S/Mg] with [Fe/H], most probably because of enrichment from Type Ia supernovae.
We present full-orbit phase curve observations of the eccentric ($e$~0.08) transiting hot Jupiter WASP-14b obtained in the 3.6 and 4.5 $\mu$m bands using the $\textit{Spitzer Space Telescope}$. We use two different methods for removing the intrapixel sensitivity effect and compare their efficacy in decoupling the instrumental noise. Our measured secondary eclipse depths of 0.1857%$\pm$0.0104% and 0.2241%$\pm$0.0087% at 3.6 and 4.5 $\mu$m, respectively, are both consistent with a blackbody temperature of 2379$\pm$55 K. We place a $2\sigma$ upper limit on the nightside flux at 3.6 $\mu$m and find it to be 10%$\pm$1% of the dayside flux, corresponding to a 1322$\pm$212 K difference in brightness temperature. At 4.5 $\mu$m, the minimum planet flux is 30%$\pm$3% of the maximum flux, corresponding to a 1016$\pm$99 K difference in brightness temperature. We compare our measured phase curves to the predictions of one-dimensional radiative transfer and three-dimensional general circulation models. We find that WASP-14b's measured dayside emission is consistent with a model atmosphere with equilibrium chemistry and a moderate temperature inversion. These same models provide a poor match to the planet's nightside emission, which is lower than predicted at 3.6 $\mu$m and higher at 4.5 $\mu$m. We propose that this discrepancy might be explained by an enhanced global C/O ratio. In addition, we find that the phase curves of WASP-14b (7.3 $M_{Jup}$) are consistent with a much lower albedo than those of other Jovian mass planets with thermal phase curve measurements, suggesting that it may be emitting detectable heat from the deep atmosphere or interior processes.
Markarian 421 (Mrk 421) is one of the brightest, fastest and closest BL Lac object known. Its very high energy (VHE) spectrum has been successfully modeled with both leptonic and hadronic models and not conclusive results have been achieved yet about the origin of its VHE emission. Here we investigate the possibility that a fraction of the VHE flares of Mrk 421 are due to hadronic processes and calculate the expected neutrino flux associated. We introduce the obtained neutrino flux in a Monte Carlo simulation to see the expectation for a Km$^{3}$ Cherenkov neutrino telescope.
The nearby active galaxy NGC 1275, has widely been detected from radio to gamma rays. Its spectral energy distribution (SED) shows a double-peak feature, which is well explained by synchrotron self-Compton (SSC) model. However, recent TeV detections might suggest that very-high-energy $\gamma$-rays (E$\geq$100 GeV) may not have a leptonic origin. We test a lepto-hadronic model to describe the whole SED through SSC emission and neutral pion decay resulting from p$\gamma$ interactions. Also, we estimate the neutrino events expected in a Km$^3$ Cherenkov telescope.
In this paper, the integrated continuum radio-spectrum of supernova remnant (SNR) W44 was analyzed up to 70 GHz, testing the different emission models that can be responsible for its particular shape. {\it Planck's} observations made possible to analyze the high frequency part of radio-emission from SNRs. Although the quality of radio-continuum spectrum (a high scatter of data points at same frequencies) prevents us to make definite conclusions, we emphasize the possibility of spinning-dust emission detection towards this remnant. In addition, a concave-down feature, due to synchrotron losses, can not be definitely dismissed by the present knowledge of the integrated radio continuum spectrum of this SNR.
In this paper, four sets of evolutionary models are computed with different values of the mixing length parameter $\alpha_{\rm p}$ and the overshooting parameter $\delta_{\rm ov}$. The properties of the convective cores and the convective envelopes are studied in the massive stars. We get three conclusions: First, the larger $\alpha_{\rm p}$ leads to enhancing the convective mixing, removing the chemical gradient, and increasing the convective heat transfer efficiency. Second, core potential $\phi_{\rm c} = M_{\rm c} / R_{\rm c}$ describes sufficiently the evolution of a star, whether it is a red or blue supergiant at central helium ignition. Third, the discontinuity of hydrogen profile above the hydrogen burning shell seriously affect the occurrence of blue loops in the Hertzsprung--Russell diagram.
This paper consists of two related parts: In the first part we derive an expression of the moment of inertia (MOI) of a neutron star as a function of observables from a hypothetical r-mode gravitational wave detection. For a given r-mode detection we show how the value of the MOI of a neutron star constrains the equation of state (EOS) of the matter in the core of the neutron star. Subsequently, for each candidate EOS, we derive a possible value of the saturation amplitude, \alpha, of the r-mode oscillations on the neutron star. Additionally, we argue that a r-mode detection will provide clues about the cooling rate mechanism of the neutron star. The above physics that can be derived from a hypothetical r-mode detection constitute our motivation for the second part of the paper. In that part we present a detection strategy to efficiently search for r-modes in gravitational-wave data. R-mode signals were injected into simulated noise colored with the advanced LIGO (aLIGO) and Einstein Telescope (ET) sensitivity curves. The r-mode waveforms used are those predicted by early theories based on a polytropic equation of state (EOS) neutron star matter. In our best case scenario \alpha of order 10^{-1}, the maximum detection distance when using the aLIGO sensitivity curve is 1 Mpc (supernova event rate of 3-4 per century) while the maximum detection distance when using the ET sensitivity curve is 10 Mpc (supernova event rate of 1-2 per year).
In many ways, the HR8799 system resembles our Solar system more closely than
any other discovered to date - albeit on a larger, younger, and more dramatic
scale - featuring four giant planets and two debris belts. The first belt lies
beyond the orbit of the outer planet, and mirrors our Solar system's
Edgeworth-Kuiper belt. The second belt lies interior to the orbit of the inner
planet, HR8799e, and is analogous to our Asteroid Belt. With such a similar
architecture, the system is a valuable laboratory for examining exoplanet
dynamics, and the interaction between debris disks and planets.
In recent years, HR8799's outer disk has been relatively well characterised,
primarily using the Herschel Space Observatory. In contrast, the inner disk,
too close to HR8799 to be spatially resolved by Herschel, remains poorly
understood. This leaves significant questions over both the location of the
planetesimals responsible for producing the observed dust, and the physical
properties of those grains.
We have performed extensive simulations of HR8799's inner, unresolved debris
belt, using UNSW Australia's supercomputing facility, Katana. Here, we present
the results of integrations following the evolution of a belt of dynamically
hot debris interior to the orbit of HR8799e, for a period of 60 Myr, using an
initial population of 500,000 massless test particles. These simulations have
enable the characterisation of the extent and structure of the inner belt,
revealing that its outer edge must lie interior to the 3:1 mean-motion
resonance with HR8799, at approximately 7.5au, and highlighting the presence of
fine structure analogous to the Solar system's Kirkwood gaps. In the future,
out results will allow us to calculate a first estimate of the small-body
impact rate and water delivery prospects for any potential terrestrial
planet(s) that might lurk, undetected, in the inner system.
We present the results of high resolution (R$\ge$30,000) optical and near-IR spectroscopic monitoring observations of HBC 722, a recent FU Orionis object that underwent an accretion burst in 2010. We observed HBC 722 in optical/near-IR with the BOES, HET-HRS, and IGRINS spectrographs, at various points in the outburst. We found atomic lines with strongly blueshifted absorption features or P Cygni profiles, both evidence of a wind driven by the accretion. Some lines show a broad double-peaked absorption feature, evidence of disk rotation. However, the wind-driven and disk-driven spectroscopic features are anti-correlated in time; the disk features became strong as the wind features disappeared. This anti-correlation might indicate that the rebuilding of the inner disk was interrupted by the wind pressure during the first two years. The Half-Width at Half-Depth (HWHD) of the double-peaked profiles decreases with wavelength, indicative of the Keplerian rotation; the optical spectra with the disk feature are fitted by a G5 template stellar spectrum convolved with a rotation velocity of 70 km s$^{-1}$ while the near-IR disk features are fitted by a K5 template stellar spectrum convolved with a rotation velocity of 50 km s$^{-1}$. Therefore, the optical and near-IR spectra seem to trace the disk at 39 and 76 $\textit{R}_{\odot}$, respectively. We fit a power-law temperature distribution in the disk, finding an index of 0.8, comparable to optically thick accretion disk models.
The nature of the supernova leading to the Crab Nebula has long been controversial because of the low energy that is present in the observed nebula. One possibility is that there is significant energy in extended fast material around the Crab but searches for such material have not led to detections. An electron capture supernova model can plausibly account for the low energy and the observed abundances in the Crab. Here, we examine the evolution of the Crab pulsar wind nebula inside a freely expanding supernova and find that the observed properties are most consistent with a low energy event. Both the velocity and radius of the shell material, and the amount of gas swept up by the pulsar wind point to a low explosion energy ($\sim 10^{50}$ ergs). We do not favor a model in which circumstellar interaction powers the supernova luminosity near maximum light because the required mass would limit the freely expanding ejecta.
Observations of the Sun from the Earth are always limited by the presence of the atmosphere, which strongly disturbs the images. A solution to this problem is to place the telescopes in space satellites, which produce observations without any (or limited) atmospheric aberrations. However, even though the images from space are not affected by atmospheric seeing, the optical properties of the instruments still limit the observations. In the case of diffraction limited observations, the PSF establishes the maximum allowed spatial resolution, defined as the distance between two nearby structures that can be properly distinguished. In addition, the shape of the PSF induce a dispersion of the light from different parts of the image, leading to what is commonly termed as stray light or dispersed light. This effect produces that light observed in a spatial location at the focal plane is a combination of the light emitted in the object at relatively distant spatial locations. We aim to correct the effect produced by the telescope's PSF using a deconvolution method, and we decided to apply the code on Hinode/SP quiet Sun observations. We analyze the validity of the deconvolution process with noisy data and we infer the physical properties of quiet Sun magnetic elements after the deconvolution process.
We provide a general analysis on the properties of emitting material of some rapidly evolving and luminous transients discovered recently with the Pan-STARRS1 Medium Deep Survey. It is found that these transients are probably produced by a low-mass non-relativistic outflow that is continuously powered by a newly born, rapidly spinning, and highly magnetized neutron star. Such a system could originate from an accretion-induced collapse of a white dwarf or a merger of a neutron star-neutron star binary. Therefore, observations to these transients would be helpful for constraining white dwarf and neutron star physics and/or for searching and identifying gravitational wave signals from the mergers.
We explore the Quasi-Periodic Pulsations (QPPs) in a solar flare observed by Fermi Gamma-ray Burst Monitor (GBM), Solar Dynamics Observatory (SDO), Solar Terrestrial Relations Observatory (STEREO), and Interface Region Imaging Spectrograph (IRIS) on 2014 September 10. QPPs are identified as the regular and periodic peaks on the rapidly-varying components, which are the light curves after removing the slowly-varying components. The QPPs display only three peaks at the beginning on the hard X-ray (HXR) emissions, but ten peaks on the chromospheric and coronal line emissions, and more than seven peaks (each peak is corresponding to a type III burst on the dynamic spectra) at the radio emissions. An uniform quasi-period about 4 minutes are detected among them. AIA imaging observations exhibit that the 4-min QPPs originate from the flare ribbon, and tend to appear on the ribbon front. IRIS spectral observations show that each peak of the QPPs tends to a broad line width and a red Doppler velocity at C I, O IV, Si IV, and Fe XXI lines. Our findings indicate that the QPPs are produced by the non-thermal electrons which are accelerated by the induced quasi-periodic magnetic reconnections in this flare.
Rezzolla et al. [ApJ 531 (2000), L139; Phys. Rev. D 64 (2001), 104013; Phys. Rev. D 64 (2001), 104014] draw attention to the second order secular drift associated with r-modes and claimed that it should lead to magnetic field enhancement and suppression of r-mode instability in magnetized neutron stars. We critically revise these results. We present a particular second order r-mode solution with vanishing secular drift, thus refuting a widely believed statement that secular drift is an unavoidable feature of r-modes. This non-drifting solution is not affected by magnetic field $B$, if $B\ll B_{\mathrm{crit}}\approx 10^{17}\,(\nu/600\,\mathrm{Hz})$ G ($\nu$ is a spin frequency) and does not lead to secular evolution of magnetic field. For general second order r-mode solution the drift does not necessarily vanish, but the solution can be presented as a superposition of two solutions: one describes evolution of differential rotation in nonoscillating star (which describes secular drift; for nonmagnetized star it is arbitrary stationary rotation stratified on cylinders; for magnetized star differential rotation evolves on the Alfv\'{e}n timescale and may lead to magnetic energy enhancement), and another one is non-drifting r-mode solution mentioned above. This representation allows us to conclude that enhancement of magnetic field energy is limited by initial energy of differential rotation, which is much less (for a factor $\propto \alpha^2$, where $\alpha$ is mode amplitude) than the total energy of r-mode. Hence, magnetic field enhancement by drift cannot suppress r-mode instability. Results can be generalized for any oscillation mode in any medium, if this mode has non-drifting solution for $B=0$.
The rotational spectral lines of c-C$_3$H$_2$ and two kinds of the $^{13}$C isotopic species, c-$^{13}$CCCH$_2$ ($C_{2v}$ symmetry) and c-CC$^{13}$CH$_2$ ($C_s$ symmetry) have been observed in the 1-3 mm band toward the low-mass star-forming region L1527. We have detected 7, 3, and 6 lines of c-C$_3$H$_2$, c-$^{13}$CCCH$_2$ , and c-CC$^{13}$CH$_2$, respectively, with the Nobeyama 45 m telescope, and 34, 6, and 13 lines, respectively, with the IRAM 30 m telescope, where 7, 2, and 2 transitions, respectively, are observed with the both telescopes. With these data, we have evaluated the column densities of the normal and $^{13}$C isotopic species. The [c-C$_3$H$_2$]/[c-$^{13}$CCCH$_2$] ratio is determined to be $310\pm80$, while the [c-C$_3$H$_2$]/[c-CC$^{13}$CH$_2$] ratio is determined to be $61\pm11$. The [c-C$_3$H$_2$]/[c-$^{13}$CCCH$_2$] and [c-C$_3$H$_2$]/[c-CC$^{13}$CH$_2$] ratios expected from the elemental $^{12}$C/$^{13}$C ratio are 60-70 and 30-35, respectively, where the latter takes into account the statistical factor of 2 for the two equivalent carbon atoms in c-C$_3$H$_2$. Hence, this observation further confirms the dilution of the $^{13}$C species in carbon-chain molecules and their related molecules, which are thought to originate from the dilution of $^{13}$C$^+$ in the gas-phase C$^+$ due to the isotope exchange reaction: $\mathrm{^{13}C^++CO\rightarrow{}^{13}CO+C^+}$. Moreover, the abundances of the two $^{13}$C isotopic species are different from each other. The ratio of c-$\mathrm{^{13}CCCH_2}$ species relative to c-$\mathrm{CC^{13}CH_2}$ is determined to be $0.20\pm0.05$. If $^{13}$C were randomly substituted for the three carbon atoms, the [c-$\mathrm{^{13}CCCH_2}$]/[c-$\mathrm{CC^{13}CH_2}$] ratio would be 0.5. Hence, the observed ratio indicates that c-$\mathrm{CC^{13}CH_2}$ exists more favorably. Possible origins of the different abundances are discussed.
The quest for the particle nature of dark matter is one of the big open questions of modern physics. The CRESST-II experiment, located at the Gran Sasso laboratory in Italy, is optimised for the detection of the elastic scattering of dark matter particles with ordinary matter. We present the result obtained with an improved detector setup with increased radiopurity and enhanced background rejection. The limit obtained in the so-called low mass region between one and three GeV/c2 is at the present among the best limits obtained for direct dark matter experiments. In addition we give an outlook of the future potential for direct dark matter detection using further improved CRESST CaWO4 cryogenic detectors.
Hot subdwarf-B stars in long-period binaries are found to be on eccentric orbits, even though current binary-evolution theory predicts those objects to be circularised before the onset of Roche-lobe overflow (RLOF). We aim to find binary-evolution mechanisms that can explain these eccentric long-period orbits, and reproduce the currently observed period-eccentricity diagram. Three different processes are considered; tidally-enhanced wind mass-loss, phase-dependent RLOF on eccentric orbits and the interaction between a circumbinary disk and the binary. The binary module of the stellar-evolution code MESA (Modules for Experiments in Stellar Astrophysics) is extended to include the eccentricity-pumping processes. The effects of different input parameters on the final period and eccentricity of a binary-evolution model are tested with MESA. The end products of models with only tidally-enhanced wind mass-loss can indeed be eccentric, but these models need to lose too much mass, and invariably end up with a helium white dwarf that is too light to ignite helium. Within the tested parameter space, no sdBs in eccentric systems are formed. Phase-dependent RLOF can reintroduce eccentricity during RLOF, and could help to populate the short-period part of the period-eccentricity diagram. When phase-dependent RLOF is combined with eccentricity pumping via a circumbinary disk, the higher eccentricities can be reached as well. A remaining problem is that these models favour a distribution of higher eccentricities at lower periods, while the observed systems show the opposite. The models presented here are potentially capable of explaining the period-eccentricity distribution of long-period sdB binaries, but further theoretical work on the physical mechanisms is necessary.
GRB 130925A is an ultra-long GRB, and it shows clear evidences for a thermal emission in the soft X-ray data of \emph{Swift}/XRT ($\sim0.5$\,keV), lasting till the X-ray afterglow phase. Due to the long duration of the GRB, the burst could be studied in hard X-rays with high-resolution focusing detectors (\emph{NuSTAR}). The blackbody temperature, as measured by the \emph{Swift}/XRT, shows a decreasing trend till the late phase (Piro et al. 2014) whereas the high-energy data reveals a significant blackbody component during the late epochs at an order of magnitude higher temperature ($\sim5$\,keV), as compared to the contemporaneous low energy data (Bellm et al. 2014). We resolve this apparent contradiction by demonstrating that a model with two black bodies and a power-law (2BBPL) is consistent with the data right from the late prompt emission to the afterglow phase. Both the blackbodies show a similar cooling behaviour upto the late time. We invoke a structured jet, having a fast spine and a slower sheath layer, to identify the location of these blackbodies. Independent of the physical interpretation, we propose that the 2BBPL model is a generic feature of the prompt emission of all long GRBs, and the thermal emission found in the afterglow phase of different GRBs reflects the lingering thermal component of the prompt emission with diverse time-scales. We strengthen this proposal by pointing out a close similarity between the spectral evolutions of this GRB and GRB~090618, a source with significant wide band data during the early afterglow phase.
We present a set of 144 galactic chemical evolution models applied to a Milky Way analogue, computed using four sets of low and intermediate star nucleosynthetic yields, six massive star yield compilations, and six functional forms for the initial mass function. The integrated or true yields for each combination are derived. A comparison is made between a grid of multiphase chemical evolution models computed with these yield combinations and empirical data drawn from the Milky Way's disc, including the solar neighbourhood. By means of a chi2 methodology, applied to the results of these multiphase models, the best combination of stellar yields and initial mass function capable of reproducing these observations is identified.
Coronal waves are an important aspect of the dynamics of the plasma in the corona. Wavefront dislocations are topological features of most waves in nature and also of magnetohydrodynamic waves. Are there dislocations in coronal waves? The finding and explanation of dislocations may shed light on the nature and characteristics of the propagating waves, their interaction in the corona and in general on the plasma dynamics. We positively identify dislocations in coronal waves observed by the Coronal Multi-channel Polarimeter (CoMP) as singularities in the Doppler shifts of emission coronal lines. We study the possible singularities that can be expected in coronal waves and try to reproduce the observed dislocations in terms of localization and frequency of appearance. The observed dislocations can only be explained by the interference of a kink and a sausage wave modes propagating with different frequencies along the coronal magnetic field. In the plane transverse to the propagation, the cross-section of the oscillating plasma must be smaller than the spatial resolution, and the two waves result in net longitudinal and transverse velocity components that are mixed through projection onto the line of sight. Alfv\'en waves can be responsible of the kink mode, but a magnetoacoustic sausage mode is necessary in all cases. Higher (flute) modes are excluded. The kink mode has a pressure amplitude that is smaller than the pressure amplitude of the sausage mode, though its observed velocity is larger. This concentrates dislocations on the top of the loop. To explain dislocations, any model of coronal waves must include the simultaneous propagation and interference of kink and sausage wave modes of comparable but different frequencies, with a sausage wave amplitude much smaller than the kink one.
We present an in-depth analysis of the flavour and spectral composition of the 36 high-energy neutrino events observed after three years of observation by the IceCube neutrino telescope. While known astrophysical sources of HE neutrinos are expected to produce a nearly $(1:1:1)$ flavour ratio (electron : muon : tau) of neutrinos at earth, we show that the best fits based on the events detected above $E_\nu \ge 28$ TeV do not necessarily support this hypothesis. Crucially, the energy range that is considered when analysing the HE neutrino data can have a profound impact on the conclusions. We highlight two intriguing puzzles: an apparent deficit of muon neutrinos, seen via a deficit of track-like events; and an absence of $\bar \nu_e$'s at high energy, seen as an absence of events near the Glashow resonance. We discuss possible explanations, including the misidentification of tracks as showers, and a broken power law, in analogy to the observed HE cosmic ray spectrum.
During the violent relaxation of a self-gravitating system a significant fraction of its mass may be ejected. If the time varying gravitational field also breaks spherical symmetry this mass can potentially carry angular momentum. Thus starting initial configurations with zero angular momentum can in principle lead to a bound virialized system with non-zero angular momentum. We explore here, using numerical simulations, how much angular momentum can be generated in a virialized structure in this way, starting from configurations of cold particles which are very close to spherically symmetric. For initial configurations in which spherical symmetry is broken only by the Poissonian fluctuations associated with the finite particle number $N$, with $N$ in range $10^3$ to $10^5$, we find that the relaxed structures have standard "spin" parameters $\lambda \sim 10^{-3}$, and decreasing slowly with $N$. For slightly ellipsoidal initial conditions, in which the finite-$N$ fluctuations break the residual reflection symmetries, we observe values $\lambda \sim 10^{-2}$ i.e. of the same order of magnitude as those reported for elliptical galaxies and dark matter halos in cosmological simulations. The net angular momentum vector is typically aligned close to normal to the major semi-axis of the triaxial relaxed structure, and also with that of the ejected mass. This simple mechanism may provide an alternative, or complement, to "tidal torque theory" for understanding the origin of angular momentum in astrophysical structures.
We analyze the dependence of the stellar disc flatness on the galaxy
morphological type using 2D decomposition of galaxies from the reliable
subsample of the Edge-on Galaxies in SDSS (EGIS) catalogue. Combining these
data with the retrieved models of the edge-on galaxies from the Two Micron All
Sky Survey (2MASS) and the Spitzer Survey of Stellar Structure in Galaxies
(S$^4$G) catalogue, we make the following conclusions:
(1) The disc relative thickness $z_0/h$ in the near- and mid-infrared
passbands correlates weakly with morphological type and does not correlate with
the bulge-to-total luminosity ratio $B/T$ in all studied bands.
(2) Applying an 1D photometric profile analysis overestimates the disc
thickness in galaxies with large bulges making an illusion of the relationship
between the disc flattening and the ratio $B/T$.
(3) In our sample the early-type disc galaxies (S0/a) have both flat and
"puffed" discs. The early spirals and intermediate-type galaxies have a large
scatter of the disc flatness, which can be caused by the presence of a bar:
barred galaxies have thicker stellar discs, on average. On the other hand, the
late-type spirals are mostly thin galaxies, whereas irregular galaxies have
puffed stellar discs.
The chromospheric response to the input of flare energy is marked by extended extreme ultraviolet (EUV) ribbons and hard X-ray (HXR) footpoints. These are usually explained as the result of heating and bremsstrahlung emission from accelerated electrons colliding in the dense chromospheric plasma. We present evidence of impulsive heating of flare ribbons above 10 MK in a two-ribbon flare. We analyse the impulsive phase of SOL2013-11-09T06:38, a C2.6 class event using data from Atmospheric Imaging Assembly (AIA) on board of Solar Dynamics Observatory (SDO) and the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) to derive the temperature, emission measure and differential emission measure of the flaring regions and investigate the evolution of the plasma in the flaring ribbons. The ribbons were visible at all SDO/AIA EUV/UV wavelengths, in particular, at 94 and 131 \AA\ filters, sensitive to temperatures of 8 MK and 12 MK. Time evolution of the emission measure of the plasma above 10 MK at the ribbons has a peak near the HXR peak time. The presence of hot plasma in the lower atmosphere is further confirmed by RHESSI imaging spectroscopy analysis, which shows resolved sources at 11-13 MK associated with at least one ribbon. We found that collisional beam heating can only marginally explain the necessary power to heat the 10 MK plasma at the ribbons.
While recent observational progress is converging on the detection of compact regions of thermal emission due to embedded protoplanets, further theoretical predictions are needed to understand the response of a protoplanetary disk to the planet formation radiative feedback. This is particularly important to make predictions for the observability of circumplanetary regions. In this work we use 2D hydrodynamical simulations to examine the evolution of a viscous protoplanetary disk in which a luminous Jupiter-mass planet is embedded. We use an energy equation which includes the radiative heating of the planet as an additional mechanism for planet formation feedback. Several models are computed for planet luminosities ranging from $10^{-5}$ to $10^{-3}$ Solar luminosities. We find that the planet radiative feedback enhances the disk's accretion rate at the planet's orbital radius, producing a hotter and more luminous environement around the planet, independently of the prescription used to model the disk's turbulent viscosity. We also estimate the thermal signature of the planet feedback for our range of planet luminosities, finding that the emitted spectrum of a purely active disk, without passive heating, is appreciably modified in the infrared. We simulate the protoplanetary disk around HD 100546 where a planet companion is located at about 68 au from the star. Assuming the planet mass is 5 Jupiter masses and its luminosity is $\sim 2.5 \times 10^{-4} \, L_\odot$, we find that the radiative feedback of the planet increases the luminosity of its $\sim 5$ au circumplanetary disk from $10^{-5} \, \rm L_\odot$ (without feedback) to $10^{-3} \, \rm L_\odot$, corresponding to an emission of $\sim 1 \, \rm mJy$ in $L^\prime$ band after radiative transfer calculations, a value that is in good agreement with HD 100546b observations.
We present high-resolution observations, using both the European VLBI Network (EVN) at 1.7-GHz, and the Very Long Baseline Array (VLBA) at 5 and 8.4-GHz to image radio structures of 14 compact sources classified as broad absorption line (BAL) quasars based on the absorption index (AI). All source but one were resolved, with the majority showing core-jet morphology typical for radio-loud quasars. We discuss in details the most interesting cases. The high radio luminosities and small linear sizes of the observed objects indicate they are strong young AGNs. Nevertheless, the distribution of the radio-loudness parameter, log(Ri), of a larger sample of AI quasars shows that the objects observed by us constitute the most luminous, small subgroup of AI population. Additionally we report that for the radio-loudness parameter, the distribution of AI quasars and those selected by using the traditional balnicity index (BI), BI quasars differ significantly. Strong absorption is connected with the lower log(Ri), and thus probably with larger viewing angles. Since, the AI quasars have on average larger log(Ri), the orientation can cause that we see them less absorbed. However, we suggest that the orientation is not the only parameter that affects the detected absorption. The fact that the strong absorption is associated with the weak radio emission is equally important and worth exploring.
We present deep Hubble Space Telescope imaging at the locations of four, potentially hostless, long-faded Type Ia supernovae (SNe Ia) in low-redshift, rich galaxy clusters that were identified in the Multi-Epoch Nearby Cluster Survey. Assuming a steep faint-end slope for the galaxy cluster luminosity function ($\alpha_d=-1.5$), our data includes all but $\lesssim0.2\%$ percent of the stellar mass in cluster galaxies ($\lesssim0.005\%$ with $\alpha_d=-1.0$), a factor of 10 better than our ground-based imaging. Two of the four SNe Ia still have no possible host galaxy associated with them ($M_R>-9.2$), confirming that their progenitors belong to the intracluster stellar population. The third SNe Ia appears near a faint disk galaxy ($M_V=-12.2$) which has a relatively high probability of being a chance alignment. A faint, red, point source coincident with the fourth SN Ia's explosion position ($M_V=-8.4$) may be either a globular cluster (GC) or faint dwarf galaxy. We estimate the local surface densities of GCs and dwarfs to show that a GC is more likely, due to the proximity of an elliptical galaxy, but neither can be ruled out. This faint host implies that the SN Ia rate in dwarfs or GCs may be enhanced, but remains within previous observational constraints. We demonstrate that our results do not preclude the use of SNe Ia as bright tracers of intracluster light at higher redshifts, but that it will be necessary to first refine the constraints on their rate in dwarfs and GCs with deep imaging for a larger sample of low-redshift, apparently hostless SNe Ia.
Binaries with hot massive components are strong X-ray sources. Besides the intrinsic X-ray emission of individual binary members originating in their winds, X-ray emission stems from the accretion on the compact companion or from wind collision. Since hot star winds are driven by the light absorption in the lines of heavier elements, wind acceleration is sensitive to the ionization state. Therefore, the over-ionization induced by external X-ray source strongly influences the winds of individual components. We studied the effect of external X-ray irradiation on hot star winds. We used our kinetic equilibrium (NLTE) wind models to estimate the influence of external X-ray ionization for different X-ray luminosities and source distances. The models are calculated for parameters typical of O stars. The influence of X-rays is given by the X-ray luminosity, by the optical depth between a given point and the X-ray source, and by a distance to the X-ray source. Therefore, the results can be interpreted in the diagrams of X-ray luminosity vs.~the optical depth parameter. X-rays are negligible in binaries with low X-ray luminosities or at large distances from the X-ray source. The influence of X-rays is stronger for higher X-ray luminosities and in closer proximity of the X-ray source. There is a forbidden area with high X-ray luminosities and low optical depth parameters, where the X-ray ionization leads to wind inhibition. There is excellent agreement between the positions of observed stars in these diagrams and our predictions. All wind-powered high-mass X-ray binary (HMXB) primaries lie outside the forbidden area. Many of them lie close to the border of the forbidden area, indicating that their X-ray luminosities are self-regulated. We discuss the implications of our work for other binary types.
Though in-orbit calibration is adopted to reduce position error of individual star spot down to 0.02pixel on star tracker, little study has been conducted on the accuracy to what extent for some significant error sources which often leads to in-orbit correction inefficiency. This study presents the general theory and estimates of the minimum error constraints, including not only on position but also on intensity and scale of Gaussian shaped profile based on Cramer Rao Lower Bound(CRLB) theory. By imposing those constraints on motion, drift in focal length and so on, margins of in-flight error sources and the final accuracy of star tracker can be analytically determined before launch.
The presence of a white dwarf in a resolved binary system, such as Sirius, provides an opportunity to combine dynamical information about the masses, from astrometry and spectroscopy, with a gravitational red-shift measurement and spectrophotometry of the white dwarf atmosphere to provide a test of theoretical mass-radius relations of unprecedented accuracy. We demonstrated this with the first Balmer line spectrum of Sirius B to be obtained free of contamination from the primary, with STIS on HST. However, we also found an unexplained discrepancy between the spectroscopic and gravitational red-shift mass determinations. With the recovery of STIS, we have been able to revisit our observations of Sirius B with an improved observation strategy designed to reduce systematic errors on the gravitational red-shift measurement. We provide a preliminary report on the refined precision of the Sirius B mass-radius measurements and the extension of this technique to a larger sample of white dwarfs in resolved binaries. Together these data can provide accurate mass and radius determinations capable of testing the theoretical mass-radius relation and distinguishing between possible structural models.
Quiescent solar prominence fine structures are typically modelled as density enhancements, called threads, which occupy a fraction of a longer magnetic flux tube. The profile of the mass density along the magnetic field is however unknown and several arbitrary alternatives are employed in prominence wave studies. We present a comparison of theoretical models for the field-aligned density along prominence fine structures. We consider Lorentzian, Gaussian, and parabolic profiles. We compare their theoretical predictions for the period ratio between the fundamental transverse kink mode and the first overtone to obtain estimates for the ratio of densities between the central part of the tube and its foot-points and to assess which one would better explain observed period ratio data. Bayesian parameter inference and model comparison techniques are developed and applied. Parameter inference requires the computation of the posterior distribution for the density gradient parameter conditional on the observable period ratio. Model comparison involves the computation of the marginal likelihood as a function of the period ratio to obtain the plausibility of each density model and the computation of Bayes Factors to quantify the relative evidence for each model, given a period ratio observation. A Lorentzian density profile, with plasma density concentrated around the centre of the tube seems to offer the most plausible inversion result. A Gaussian profile would require unrealistically large values of the density gradient parameter and a parabolic density distribution does not enable us to obtain well constrained posterior probability distributions. However, our model comparison results indicate that the evidence points to the Gaussian and parabolic profiles for period ratios in between 2 and 3, while the Lorentzian profile is preferred for larger period ratio values.
We investigate the nature of the long-period radial velocity variations in Alpha Tau first reported over 20 years ago. We analyzed precise stellar radial velocity measurements for Alpha Tau spanning over 30 years. An examination of the Halpha and Ca II 8662 spectral lines, and Hipparcos photometry was also done to help discern the nature of the long-period radial velocity variations. Our radial velocity data show that the long-period, low amplitude radial velocity variations are long-lived and coherent. Furthermore, Halpha equivalent width measurements and Hipparcos photometry show no significant variations with this period. Another investigation of this star established that there was no variability in the spectral line shapes with the radial velocity period. An orbital solution results in a period of P = 628.96 +/- 0.90 d, eccentricity, e = 0.10 +/- 0.05, and a radial velocity amplitude, K = 142.1 +/- 7.2 m/s. Evolutionary tracks yield a stellar mass of 1.13 +/- 0.11 M_sun, which corresponds to a minimum companion mass of 6.47 +/- 0.53 M_Jup with an orbital semi-major axis of a = 1.46 +/- 0.27 AU. After removing the orbital motion of the companion, an additional period of ~ 520 d is found in the radial velocity data, but only in some time spans. A similar period is found in the variations in the equivalent width of Halpha and Ca II. Variations at one-third of this period are also found in the spectral line bisector measurements. The 520 d period is interpreted as the rotation modulation by stellar surface structure. Its presence, however, may not be long-lived, and it only appears in epochs of the radial velocity data separated by $\sim$ 10 years. This might be due to an activity cycle. The data presented here provide further evidence of a planetary companion to Alpha Tau, as well as activity-related radial velocity variations.
Magnetar flares generate Alfv\'en waves bouncing in the closed magnetosphere with energy up to $\sim 10^{46}$ erg. We show that on a 10-ms timescale the waves are transmitted into the star and form a compressed packet of high energy density. This packet strongly shears the stellar crust and initiates a plastic flow, heating the crust and melting it hundreds of meters below the surface. A fraction of the deposited plastic heat is eventually conducted to the stellar surface, contributing to the surface afterglow months to years after the flare. A large fraction of heat is lost to neutrino emission or conducted into the core of the neutron star.
We determine the concentration-mass relation of 19 X-ray selected galaxy clusters from the CLASH survey in theories of gravity that directly modify the lensing potential. We model the clusters as NFW haloes and fit their lensing signal, in the Cubic Galileon and Nonlocal gravity models, to the lensing convergence profiles of the clusters. We discuss a number of important issues that need to be taken into account, associated with the use of nonparametric and parametric lensing methods, as well as assumptions about the background cosmology. Our results show that the concentration and mass estimates in the modified gravity models are, within the errorbars, the same as in $\Lambda$CDM. This result demonstrates that, for the Nonlocal model, the modifications to gravity are too weak at the cluster redshifts, and for the Galileon model, the screening mechanism is very efficient inside the cluster radius. However, at distances $\sim \left[2-20\right] {\rm Mpc}/h$ from the cluster center, we find that the surrounding force profiles are enhanced by $\sim20-40\%$ in the Cubic Galileon model. This has an impact on dynamical mass estimates, which means that tests of gravity based on comparisons between lensing and dynamical masses can also be applied to the Cubic Galileon model.
Supernova (SN) 2002hh was unusual among core-collapse SNe because it was highly reddened, and displayed a bright infrared (IR) excess due to radiatively heated dust in its circumstellar medium (CSM). Estimates for the mass of dust responsible for the IR echo suggested the presence of a massive shell within 0.26 pc of the star. For a velocity of 5000 - 10000 km/s, this material should be hit by the SN blast wave at late times, starting at roughly 12 years post-explosion. We have obtained deep late-time spectra with the MMT Blue Channel spectrograph to search for any spectral signatures of ongoing shock interaction. Interaction with a strength comparable to SN 1987A's collision with the equatorial ring would be detected in our data. However, in the spectra reported here, we do not detect clear signs of strong CSM interaction, contrary to expectations based on the reported radii of the dust shell. We do, however, detect emission associated with the old SN, and we find that the broad lines in the spectrum indicate a continuation of an ongoing reflected light echo, which appears similar to the spectrum at peak luminosity for this Type II-P event.
We report the first detection of >100 MeV gamma rays associated with a behind-the-limb solar flare, which presents a unique opportunity to probe the underlying physics of high-energy flare emission and particle acceleration. On 2013 October 11 a GOES M1.5 class solar flare occurred ~ 9.9 degrees behind the solar limb as observed by STEREO-B. RHESSI observed hard X-ray emission above the limb, most likely from the flare loop-top, as the footpoints were occulted. Surprisingly, the Fermi Large Area Telescope (LAT) detected >100 MeV gamma-rays for ~30 minutes with energies up to GeV. The LAT emission centroid is consistent with the RHESSI hard X-ray source, but its uncertainty does not constrain the source to be located there. The gamma-ray spectra can be adequately described by bremsstrahlung radiation from relativistic electrons having a relatively hard power-law spectrum with a high-energy exponential cutoff, or by the decay of pions produced by accelerated protons and ions with an isotropic pitch-angle distribution and a power-law spectrum with a number index of ~3.8. We show that high optical depths rule out the gamma rays originating from the flare site and a high-corona trap model requires very unusual conditions, so a scenario in which some of the particles accelerated by the CME shock travel to the visible side of the Sun to produce the observed gamma rays may be at work.
In this paper, we present a detailed high-resolution spectroscopic study of post main sequence stars in the Globular Cluster M68. Our sample, which covers a range of 4000 K in $T_{eff}$, and 3.5 dex in $log(g)$, is comprised of members from the red giant, red horizontal, and blue horizontal branch, making this the first high-resolution globular cluster study covering such a large evolutionary and parameter space. Initially, atmospheric parameters were determined using photometric as well as spectroscopic methods, both of which resulted in unphysical and unexpected $T_{eff}$, $log(g)$, $\xi_{t}$, and [Fe/H] combinations. We therefore developed a hybrid approach that addresses most of these problems, and yields atmospheric parameters that agree well with other measurements in the literature. Furthermore, our derived stellar metallicities are consistent across all evolutionary stages, with $\langle$[Fe/H]$\rangle$ = $-$2.42 ($\sigma$ = 0.14) from 25 stars. Chemical abundances obtained using our methodology also agree with previous studies and bear all the hallmarks of globular clusters, such as a Na-O anti-correlation, constant Ca abundances, and mild $r$-process enrichment.
A model equation for the Reynolds number dependence of the dimensionless dissipation rate in freely decaying homogeneous magnetohydrodynamic turbulence in the absence of a mean magnetic field is derived from the real-space energy balance equation, leading to $C_{\varepsilon}=C_{\varepsilon, \infty}+C/R_- +O(1/R_-^2))$, where $R_-$ is a generalized Reynolds number. The constant $C_{\varepsilon, \infty}$ describes the total energy transfer flux. This flux depends on magnetic and cross helicities, because these affect the nonlinear transfer of energy, suggesting that the value of $C_{\varepsilon,\infty}$ is not universal. Direct numerical simulations were conducted on up to $2048^3$ grid points, showing good agreement between data and the model. The model suggests that the magnitude of cosmological-scale magnetic fields is controlled by the values of the vector field correlations. The ideas introduced here can be used to derive similar model equations for other turbulent systems.
Many modern cosmological scenarios feature large volumes of spacetime in a de Sitter vacuum phase. Such models are said to be faced with a "Boltzmann Brain problem" - the overwhelming majority of observers with fixed local conditions are random fluctuations in the de Sitter vacuum, rather than arising via thermodynamically sensible evolution from a low-entropy past. We argue that this worry can be straightforwardly avoided in the Many-Worlds (Everett) approach to quantum mechanics, as long as the underlying Hilbert space is infinite-dimensional. In that case, de Sitter settles into a truly stationary quantum vacuum state. While there would be a nonzero probability for observing Boltzmann-Brain-like fluctuations in such a state, "observation" refers to a specific kind of dynamical process that does not occur in the vacuum (which is, after all, time-independent). Observers are necessarily out-of-equilibrium physical systems, which are absent in the vacuum. Hence, the fact that projection operators corresponding to states with observers in them do not annihilate the vacuum does not imply that such observers actually come into existence. The Boltzmann Brain problem is therefore much less generic than has been supposed.
We study the contribution of Galactic sources to the flux of astrophysical neutrinos recently observed by the IceCube Collaboration. We show that the Galactic diffuse neutrino emission consistent with $\gamma$-ray (Fermi-LAT) and cosmic ray data (KASCADE, KASCADE-Grande and CREAM) is expected to account for only 4%$-$8% of the IceCube flux above 60 TeV. Direct neutrino emission from cosmic ray-gas ($pp$) interactions in the sources would require an unusually large average opacity above 0.01. On the other hand, we find that the IceCube events already probe Galactic neutrino scenarios via the distribution of event arrival directions. We show that most Galactic scenarios can only have a limited contribution to the astrophysical signal: diffuse Galactic emission ($\lesssim50$%), quasi-diffuse emission of neutrino sources ($\lesssim65$%), extended diffuse emission from the Fermi Bubbles ($\lesssim25$%) or unidentified TeV $\gamma$-ray sources ($\lesssim25$%). Presently, dark matter decay remains unconstrained.
In this chapter, we present some reflections about the learning process -and its implications in the teaching- of notions related to the phases of the Moon.
In this paper we present a dynamical analysis for a coupled tachyonic dark energy with dark matter. The tachyonic field $\phi$ is considered in the presence of barothropic fluids (matter and radiation) and the autonomous system due to the evolution equations is studied. The three cosmological eras (radiation, matter and dark energy) are described through the critical points, for a generic potential $V(\phi)$.
In the last 10 years, we have observed an important increase of interest in the application of time-dependent energy density functional theory (TD-EDF). This approach allows to treat nuclear structure and nuclear reaction from small to large amplitude dynamics in a unified framework. The possibility to perform unrestricted three-dimensional simulations using state of the art effective interactions has opened new perspectives. In the present article, an overview of applications where the predictive power of TD-EDF has been benchmarked is given. A special emphasize is made on processes that are of astrophysical interest. Illustrations discussed here include giant resonances, fission, binary and ternary collisions leading to fusion, transfer and deep inelastic processes.
The statistical response of a Kerr black hole to incoming quantum radiation has heretofore been studied by the methods of maximum entropy or quantum field theory in curved spacetime. Neither approach pretends to take into account the quantum structure of the black hole itself. To address this last issue we calculate here the conditional probability distribution associated with the hole's response by assuming that the horizon area has a discrete quantum spectrum, and that its quantum evolution corresponds to jumps between adjacent area eigenvalues, possibly occurring in series, with consequent emission or absorption of quanta, possibly in the same mode. This "atomic" model of the black hole is implemented in two different ways and recovers the previously calculated radiation statistics in both cases. The corresponding conditional probably distribution is here expressed in closed form in terms of an hypergeometric function.
$f(R)$ gravity is one of the simplest theories of modified gravity to explain the accelerated cosmic expansion. Although it is usually assumed that the quasi-Newtonian approach for cosmic perturbations is good enough to describe the evolution of large scale structure in $f(R)$ models, some studies have suggested that this method is not valid for all $f(R)$ models. Here, we show that in the matter-dominated era, the pressure and shear equations alone, which can be recast into four first-order equations to solve for cosmological perturbations exactly, are sufficient to solve for the Newtonian potential, $\Psi$, and the curvature potential, $\Phi$. Based on these two equations, we are able to clarify how the exact linear perturbations fit into different limits. We find that in the subhorizon limit, the so called quasi-static assumption plays no role in reducing the exact linear perturbations in any viable $f(R)$ gravity. Our findings also disagree with previous studies where we find little difference between our exact solutions and that of the quasi-Newtonian approach even up to $k=10 c^{-1} H_0$.
In a series of papers Kallosh, Linde, and collaborators have provided a unified description of single-field inflation with several types of potentials, ranging from power law to supergravity, in terms of just one parameter $\alpha$. These so-called $\alpha$-attractors predict a spectral index $n_{s}$ and a tensor-to-scalar ratio $r$, which are fully compatible with the latest Planck data. The only common feature of all $\alpha$-attractors is the analyticity of the scalar potential in the non-canonical Einstein frame. In this paper we explore the case of non-analytic potentials and we find that they lead to a class of attractors characterized by quasi-scale invariance in the Jordan frame. In the canonical Einstein frame they all converge to a model with a linear potential and a universal relation between $r$ and $n_{s}$ that can fit the observational data. We show that the breaking of exact, classical, scale invariance in the Jordan frame can be attributed to one-loop corrections, in line with previous results.
Undergraduate animation students at the Maryland Institute College of Art teamed up with scientists from the Fermi Gamma-ray Space Telescope to produce a set of animations on several astronomy topics. We describe the process and discuss the results, including educational benefits and the cross-cultural experience. These animations are freely available online.
Thusfar, there does not appear to be an agreed definition of homogeneous dark energy (DE). In this work, we argue that a correct definition of homogeneous DE is one whose density perturbation in comoving gauge vanishes. Using different DE models, we then investigate the consequence of this approach in the power spectrum -- with all the power spectra being normalized to match each other on small scales, at z = 0. We find that on super-Hubble scales, relativistic corrections in the observed galaxy power spectrum are able to distinguish a homogeneous DE from the concordance model and from a clustering DE, at low z and for high magnification bias. However, the matter power spectrum: is incapable of distinguishing a homogeneous DE from the concordance model (on all scales), at z = 0; but is able to differentiate it from a clustering DE, particularly at low z. Moreover, we found that relativistic effects become enhanced with decreasing magnification bias, and with increasing z.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
An analytical model is presented for calculating the surface density as a function of radius $\Sigma(r)$ in protoplanetary disks in which a planet has opened a gap. This model is also applicable to circumbinary disks with extreme binary mass ratios. The gap profile can be solved for algebraically, without performing any numerical integrals. In contrast with previous one-dimensional gap models, this model correctly predicts that low-mass (sub-Jupiter) planets can open gaps in sufficiently low-viscosity disks, and it correctly recovers the power-law dependence of gap depth on planet-to-star mass ratio $q$, disk aspect ratio $h/r$, and dimensionless viscosity $\alpha$ found in previous numerical studies. Analytical gap profiles are compared with numerical calculations over a range of parameter space in $q$, $h/r$, and $\alpha$, demonstrating accurate reproduction of the "partial gap" regime, and general agreement over a wide range of parameter space.
The next generation of large sky photometric surveys will finally be able to use arc statistics as a cosmological probe. Here we present the first of a series of papers on this topic. In particular, we study how arc counts are sensitive to the variation of two cosmological parameters: the (total) matter density parameter, $\Omega_m$, and the normalisation of the primordial power spectrum, expressed in terms of $\sigma_8$. Both these parameters influence the abundances of collapsed structures and their internal structure. We compute the expected number of gravitational arcs with various length-to-width ratios in mock light cones, by varying these cosmological parameters in the ranges $0.1\leq\Omega_m\leq0.5$ and $0.6\leq\sigma_8\leq 1$. We find that the arc counts dependence on $\Omega_m$ and $\sigma_8$ is similar, but not identical, to that of the halo counts. We investigate how the precision of the constraints on the cosmological parameters based on arc counts depends on the survey area. We find that the constraining power of arc statistics degrades critically only for surveys covering an area smaller than $10\%$ of the whole sky. Finally, we consider the case in which the search for arcs is done only in frames where galaxy clusters have been previously identified. Adopting the selection function for galaxy clusters expected to be detected from photometric data in future wide surveys, we find that less than $10\%$ of the arcs will be missed, with only a small degradation of the corresponding cosmological constraints.
The water ice or snow line is one of the key properties of protoplanetary disks that determines the water content of terrestrial planets in the habitable zone. Its location is determined by the properties of the star, the mass accretion rate through the disk, and the size distribution of dust suspended in the disk. We calculate the snow line location from recent observations of mass accretion rates and as a function of stellar mass. By taking the observed dispersion in mass accretion rates as a measure of the dispersion in initial disk mass, we find that stars of a given mass will exhibit a range of snow line locations. At a given age and stellar mass, the observed dispersion in mass accretion rates of 0.4 dex naturally leads to a dispersion in snow line locations of 0.2 dex. For ISM-like dust sizes, the one-sigma snow line location among solar mass stars of the same age ranges from 2 to 5 au. For more realistic dust opacities that include larger grains, the snow line is located up to two times closer to the star. We use these locations and the outcome of N-body simulations to predict the amount of water delivered to terrestrial planets that formed in situ in the habitable zone. We find that the dispersion in snow line locations leads to a large range in water content. For ISM-like dust sizes, a significant fraction of habitable-zone terrestrial planets around sun-like stars remain dry, and no water is delivered to the habitable zones of low-mass M stars (less than half a solar mass) as in previous works. The closer-in snow line in disks with larger grains enables water delivery to the habitable zone for a significant fraction of M stars and all FGK stars. Considering their larger numbers and higher planet occurrence, M stars may host most of the water-rich terrestrial planets in the galaxy if these planets are able to hold on to their water in their subsequent evolution.
We present optical and near-infrared photometry of GRB 140606B ($z=0.384$), and optical photometry and spectroscopy of its associated supernova (SN). The bolometric properties of the SN are: a nickel mass of M$_{\rm Ni}$=0.4$\pm$0.2 M$_{\odot}$, an ejecta mass of M$_{\rm ej}$=5$\pm$2 M$_{\odot}$, and a kinetic energy of E$_{\rm K}$=2$\pm1\times10^{52}$ erg. The uncertain value of M$_{\rm Ni}$ is primarily due to the poorly constrained rest-frame extinction ($E(B-V)_{\rm rest}$=0.16$\pm$0.14 mag). The photospheric velocity of the SN near maximum light is $v_{\rm ph}\approx$20,000 km/s. The photospheric velocity and bolometric properties are fully consistent with the statistical averages determined for other GRB-SNe. However, in terms of its $\gamma$-ray emission, GRB 140606B is an outlier of the Amati relation, and occupies the same region as low-luminosity ($ll$) and short GRBs. The $\gamma$-ray emission in $ll$GRBs is thought to arise, at least in some events, from a shock-breakout (SBO), rather than from a jet. The measured peak photon energy is $E_{\rm p}\approx800$ keV, which is close to the value expected for gamma-rays created by a SBO ($\approx 1$ MeV). Moreover, based on its position in the $M_{V,\rm p}$--$L_{\rm iso,\gamma}$ plane and the $E_{\rm K}$--$\Gamma\beta$ plane, GRB 140606B has properties similar to both SBO-GRBs and jetted-GRBs. Finally, we searched for correlations between the isotropic $\gamma$-ray emission and the bolometric properties of a sample of GRB-SNe, finding that no statistically significant correlation is present. The average kinetic energy of the sample is <$E_{\rm K}$>=2.1$\times10^{52}$ erg. All of the GRB-SNe in our sample, with the exception of SN 2006aj, are within this range, which has implications for the total energy budget available to power both the relativistic and non-relativistic components in a GRB-SN event. [abridged]
We present statistics of 89 faint 1.2-mm continuum sources with a flux density of ~0.01-1 mJy detected by about 100 deep ALMA pointing data that include the complete deep datasets archived by 2015 March. These faint sources are identified in 50 blank fields and behind one cluster, Abell 1689, that magnifies the background sources by gravitational lensing. Evaluating various important effects including the false detection, detection completeness, and flux boosting as well as the lensing magnification by modeling and simulations, we derive number counts of 1.2 mm continuum sources. We find that the number counts are well represented by the Schechter function down to ~0.01 mJy, and that the total integrated 1.2 mm flux of the securely identified sources is 22.8^(+6.1)_(-6.4) Jy deg^(-2) that corresponds to 104^(+27)_(-30)% of the extragalactic background light (EBL) measured by COBE observations. These results suggest that the major 1.2 mm EBL contributors are sources with >~0.01 mJy, and that very faint 1.2 mm sources with <~ 0.01 mJy contribute negligibly to the EBL with the possible flattening and/or truncation of number counts in this very faint flux regime. To understand the physical origin of our faint ALMA sources, we measure the galaxy bias bg by the counts-in-cells technique under the assumption that the sources reside at z~2.5, and place a stringent upper limit of bg<4.1 that is not similar to bg values of massive DRGs and SMGs but comparable to those of UV-bright sBzKs and LBGs. Moreover, in optical and near-infrared (NIR) deep fields, we identify optical-NIR counterparts for 54% of our faint ALMA sources, majority of which have luminosities and colors same as sBzKs and LBGs. We thus conclude that about a half of our faint ALMA sources are dust-poor high-z galaxies as known as sBzKs and LBGs in optical studies.
We present new NuSTAR and Chandra observations of NGC 3393, a galaxy reported to host the smallest separation dual AGN resolved in the X-rays. While past results suggested a 150 pc separation dual AGN, three times deeper Chandra imaging, combined with adaptive optics and radio imaging suggest a single, heavily obscured, radio-bright AGN. Using VLA and VLBA data, we find an AGN with a two-sided jet rather than a dual AGN and that the hard X-ray, UV, optical, NIR, and radio emission are all from a single point source with a radius <0.2". We find that the previously reported dual AGN is most likely a spurious detection resulting from the low number of X-ray counts (<160) at 6-7 keV and Gaussian smoothing of the data on scales much smaller than the PSF (0.25" vs. 0.80" FWHM). We show that statistical noise in a single Chandra PSF generates spurious dual peaks of the same separation (0.55$\pm$0.07" vs. 0.6") and flux ratio (39$\pm$9% vs. 32% of counts) as the purported dual AGN. With NuSTAR, we measure a Compton-thick source (NH=$2.2\pm0.4\times10^{24}$ cm$^{-2}$) with a large torus half-opening angle, {\theta}=79 which we postulate results from feedback from strong radio jets. This AGN shows a 2-10 keV intrinsic to observed flux ratio of 150. Using simulations, we find that even the deepest Chandra observations would severely underestimate the intrinsic luminosity of NGC 3393 above z>0.2, but would detect an unobscured AGN of this luminosity out to high redshift (z=5).
In this review we discuss the population of stellar-mass black holes in our galaxy and beyond, which are the extreme endpoints of massive star evolution. In particular we focus on how we can attempt to balance the available accretion energy with feedback to the environment via radiation, jets and winds, considering also possible contributions to the energy balance from black hole spin and advection. We review quantitatively the methods which are used to estimate these quantities, regardless of the details of the astrophysics close to the black hole. Once these methods have been outlined, we work through an outburst of a black hole X-ray binary system, estimating the flow of mass and energy through the different accretion rates and states. While we focus on feedback from stellar mass black holes in X-ray binary systems, we also consider the applicability of what we have learned to supermassive black holes in active galactic nuclei. As an important control sample we also review the coupling between accretion and feedback in neutron stars, and show that it is very similar to that observed in black holes, which strongly constrains how much of the astrophysics of feedback can be unique to black holes.
We have used FMOS on Subaru to obtain near-infrared spectroscopy of 123 far-infrared selected galaxies in COSMOS and obtain the key rest-frame optical emission lines. This is the largest sample of infrared galaxies with near-infrared spectroscopy at these redshifts. The far-infrared selection results in a sample of galaxies that are massive systems that span a range of metallicities in comparison with previous optically selected surveys, and thus has a higher AGN fraction and better samples the AGN branch. We establish the presence of AGN and starbursts in this sample of (U)LIRGs selected as Herschel-PACS and Spitzer-MIPS detections in two redshift bins (z~0.7 and z~1.5) and test the redshift dependence of diagnostics used to separate AGN from star-formation dominated galaxies. In addition, we construct a low redshift (z~0.1) comparison sample of infrared selected galaxies and find that the evolution from z~1.5 to today is consistent with an evolving AGN selection line and a range of ISM conditions and metallicities from the models of Kewley et al. (2013b). We find that a large fraction of (U)LIRGs are BPT-selected AGN using their new, redshift-dependent classification line. We compare the position of known X-ray detected AGN (67 in total) with the BPT selection and find that the new classification line accurately selects most of these objects (> 70%). Furthermore, we identify 35 new (likely obscured) AGN not selected as such by their X-ray emission. Our results have direct implications for AGN selection at higher redshift with either current (MOSFIRE,
We describe a novel and fully Bayesian approach to detect halos of different masses in cosmological observations and to quantify corresponding uncertainties. To demonstrate the capability of this approach we perform a Bayesian analysis of a realistic spectroscopic galaxy mock catalogue with the previously developed HADES inference algorithm. This procedure provides us with accurately inferred three-dimensional density fields and corresponding quantification of uncertainties inherent to any cosmological observation. Based upon these results we develop a novel Bayesian methodology to detect halos of different masses in cosmological observations subject to noise and systematic uncertainties. Specifically, we use a Bayesian chain rule to connect properties of halos found in simulations with actual observations. In an entirely Bayesian manner, this approach returns detection probabilities for halos above specified mass thresholds throughout the entire observed domain. We present maps of such detection probabilities and demonstrate the validity of this approach within our mock scenario. The proposed methodology can easily be extended to account for more complex scientific questions and is a promising novel tool to analyse the cosmic large-scale structure in observations.
We present and make publicly available the first data release (DR1) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. The KODIAQ survey is aimed at studying galactic and circumgalactic gas in absorption at high-redshift, with a focus on highly-ionized gas traced by OVI, using the HIRES spectrograph on the Keck-I telescope. KODIAQ DR1 consists of a fully-reduced sample of 170 quasars at 0.29 < z_em < 5.29 observed with HIRES at high resolution (36,000 <= R <= 103,000) between 2004 and 2012. DR1 contains 247 spectra available in continuum normalized form, representing a sum total exposure time of ~1.6 megaseconds. These co-added spectra arise from a total of 567 individual exposures of quasars taken from the Keck Observatory Archive (KOA) in raw form and uniformly processed using a HIRES data reduction package made available through the XIDL distribution. DR1 is publicly available to the community, housed as a higher level science product at the KOA. We will provide future data releases that make further QSOs, including those with pre-2004 observations taken with the previous-generation HIRES detectors.
We report on the determination of the astrometric, spin and orbital parameters for PSR J1953+1846A, a "black widow" binary millisecond pulsar in the globular cluster M71. By using the accurate position and orbital parameters obtained from radio timing, we identified the optical companion in ACS/Hubble Space Telescope images. It turns out to be a faint (m_F606W>=24, m_F814W>=23) and variable star located at only ~0.06" from the pulsar timing position. The light curve shows a maximum at the pulsar inferior conjunction and a minimum at the pulsar superior conjunction, thus confirming the association with the system. The shape of the optical modulation suggests that the companion star is heated, likely by the pulsar wind. The comparison with the X-ray light curve possibly suggests the presence of an intra-binary shock due to the interaction between the pulsar wind and the material released by the companion. This is the second identification (after COM-M5C) of an optical companion to a black widow pulsar in a globular cluster. Interestingly, the two companions show a similar light curve and share the same position in the color magnitude diagram.
We present a multiwavelength morphological analysis of star forming clouds and filaments in the central ($< 50$ kpc) regions of 16 low redshift ($z<0.3$) cool core brightest cluster galaxies (BCGs). New Hubble Space Telescope (HST) imaging of far ultraviolet continuum emission from young ($\sim 10$ Myr), massive ($> 5$ \Msol) stars reveals filamentary and clumpy morphologies, which we quantify by means of structural indices. The FUV data are compared with X-ray, Ly$\alpha$, narrowband H$\alpha$, broadband optical/IR, and radio maps, providing a high spatial resolution atlas of star formation locales relative to the ambient hot ($\sim10^{7-8}$ K) and warm ionised ($\sim 10^4$ K) gas phases, as well as the old stellar population and radio-bright AGN outflows. Nearly half of the sample possesses kpc-scale filaments that, in projection, extend toward and around radio lobes and/or X-ray cavities. These filaments may have been uplifted by the propagating jet or buoyant X-ray bubble, or may have formed {\it in situ} by cloud collapse at the interface of a radio lobe or rapid cooling in a cavity's compressed shell. The morphological diversity of nearly the entire FUV sample is reproduced by recent hydrodynamical simulations in which the AGN powers a self-regulating rain of thermally unstable star forming clouds that precipitate from the hot atmosphere. In this model, precipitation triggers where the cooling-to- freefall time ratio is $t_{\mathrm{cool}}/t_{\mathrm{ff}}\sim 10$. This condition is roughly met at the maxmial projected FUV radius for more than half of our sample, and clustering about this ratio is stronger for sources with higher star formation rates.
The Spitzer Survey of Stellar Structure in Galaxies (S4G) is a volume, magnitude, and size-limited survey of 2352 nearby galaxies with deep imaging at 3.6 and 4.5um. In this paper we describe our surface photometry pipeline and showcase the associated data products that we have released to the community. We also identify the physical mechanisms leading to different levels of central stellar mass concentration for galaxies with the same total stellar mass. Finally, we derive the local stellar mass-size relation at 3.6um for galaxies of different morphologies. Our radial profiles reach stellar mass surface densities below 1 Msun pc-2. Given the negligible impact of dust and the almost constant mass-to-light ratio at these wavelengths, these profiles constitute an accurate inventory of the radial distribution of stellar mass in nearby galaxies. From these profiles we have also derived global properties such as asymptotic magnitudes (and the corresponding stellar masses), isophotal sizes and shapes, and concentration indices. These and other data products from our various pipelines (science-ready mosaics, object masks, 2D image decompositions, and stellar mass maps), can be publicly accessed at IRSA (this http URL).
We investigate the occurrence rate of rapidly rotating ($v\sin i$$>$10 km s$^{-1}$), low-mass giant stars in the APOGEE-Kepler (APOKASC) fields with asteroseismic mass and surface gravity measurements. Such stars are likely merger products and their frequency places interesting constraints on stellar population models. We also identify anomalous rotators, i.e. stars with 5 km s$^{-1}$$<$$v\sin i$$<$10 km s$^{-1}$ that are rotating significantly faster than both angular momentum evolution predictions and the measured rates of similar stars. Our data set contains fewer rapid rotators than one would expect given measurements of the Galactic field star population, which likely indicates that asteroseismic detections are less common in rapidly rotating red giants. The number of low-mass moderate (5-10 km s$^{-1}$) rotators in our sample gives a lower limit of 7% for the rate at which low-mass stars interact on the upper red giant branch because single stars in this mass range are expected to rotate slowly. Finally, we classify the likely origin of the rapid or anomalous rotation where possible. KIC 10293335 is identified as a merger product and KIC 6501237 is a possible binary system of two oscillating red giants.
B and V time-series photometry of the M31 dwarf spheroidal satellite Andromeda XXI (And XXI) was obtained with the Large Binocular Cameras at the Large Binocular Telescope. We have identified 50 variables in And XXI, of which 41 are RR Lyrae stars (37 fundamental-mode RRab, and 4 first-overtone RRc, pulsators) and 9 are Anomalous Cepheids (ACs). The average period of the RRab stars (<Pab> = 0.64 days) and the period-amplitude diagram place And~XXI in the class of Oosterhoff II - Oosterhoff-Intermediate objects. From the average luminosity of the RR Lyrae stars we derived the galaxy distance modulus of (m-M)$_0$=$24.40\pm0.17$ mag, which is smaller than previous literature estimates, although still consistent with them within 1 $\sigma$. The galaxy color-magnitude diagram shows evidence for the presence of three different stellar generations in And~XXI: 1) an old ($\sim$ 12 Gyr) and metal poor ([Fe/H]=$-$1.7 dex) component traced by the RR Lyrae stars; 2) a slightly younger (10-6 Gyr) and more metal rich ([Fe/H]=$-$1.5 dex) component populating the red horizontal branch, and 3) a young age ($\sim$ 1 Gyr) component with same metallicity, that produced the ACs. Finally, we provide hints that And~XXI could be the result of a minor merging event between two dwarf galaxies.
We report the trigonometric parallax of IRAS 07427-2400 with VERA to be 0.185 $\pm$ 0.027 mas, corresponding to a distance of 5.41$^{+0.92}_{-0.69}$ kpc. The result is consistent with the previous result of 5.32$^{+0.49}_{-0.42}$ kpc obtained by Choi et al. (2014) within error. To remove the effect of internal maser motions (e.g., random motions), we observed six maser features associated with IRAS 07427-2400 and determined systematic proper motions of the source by averaging proper motions of the six maser features. The obtained proper motions are ($\mu_{\alpha}$cos$\delta$, $\mu_{\delta}$) = ($-$1.79 $\pm$ 0.32, 2.60 $\pm$ 0.17) mas yr$^{-1}$ in equatorial coordinates, while Choi et al. (2014) showed ($\mu_{\alpha}$cos$\delta$, $\mu_{\delta}$) = ($-$2.43 $\pm$ 0.02, 2.49 $\pm$ 0.09) mas yr$^{-1}$ with one maser feature. Our astrometry results place the source in the Perseus arm, the nearest main arm in the Milky Way. Using our result with previous astrometry results obtained from observations of the Perseus arm, we conducted direct (quantitative) comparisons between 27 astrometry results and an analytic gas dynamics model based on the density-wave theory and obtained two results. First is the pitch angle of the Perseus arm determined by VLBI astrometry, 11.1 $\pm$ 1.4 deg, differing from what is determined by the spiral potential model (probably traced by stars), $\sim$ 20 deg. The second is an offset between a dense gas region and the bottom of the spiral potential model. The dense gas region traced by VLBI astrometry is located downstream of the spiral potential model, which was previously confirmed in the nearby grand-design spiral galaxy M51 in Egusa et al. (2011).
We present a new and simple method to measure the instantaneous mass and radial growth rates of the stellar discs of spiral galaxies, based on their star formation rate surface density (SFRD) profiles. Under the hypothesis that discs are exponential with time-varying scalelengths, we derive a universal theoretical profile for the SFRD, with a linear dependence on two parameters: the specific mass growth rate $\nu_\textrm{M} \equiv \dot{M_\star}/M_\star$ and the specific radial growth rate $\nu_\textrm{R} \equiv \dot{R}_\star/R_\star$ of the disc. We test our theory on a sample of 35 nearby spiral galaxies, for which we derive a measurement of $\nu_\textrm{M}$ and $\nu_\textrm{R}$. 32/35 galaxies show the signature of ongoing inside-out growth ($\nu_\textrm{R} > 0$). The typical derived e-folding timescales for mass and radial growth in our sample are ~ 10 Gyr and ~ 30 Gyr, respectively, with some systematic uncertainties. More massive discs have a larger scatter in $\nu_\textrm{M}$ and $\nu_\textrm{R}$, biased towards a slower growth, both in mass and size. We find a linear relation between the two growth rates, indicating that our galaxy discs grow in size at ~ 0.35 times the rate at which they grow in mass; this ratio is largely unaffected by systematics. Our results are in very good agreement with theoretical expectations if known scaling relations of disc galaxies are not evolving with time.
We report on 230 GHz (1.3 mm) VLBI observations of M87 with the Event Horizon Telescope using antennas on Mauna Kea in Hawaii, Mt. Graham in Arizona and Cedar Flat in California. For the first time, we have acquired 230 GHz VLBI interferometric phase information on M87 through measurement of closure phase on the triangle of long baselines. Most of the measured closure phases are consistent with 0$^{\circ}$ as expected by physically-motivated models for 230 GHz structure such as jet models and accretion disk models. The brightness temperature of the event-horizon-scale structure is $\sim 1 \times 10^{10}$ K derived from the compact flux density of $\sim 1$ Jy and the angular size of $\sim 40 $ $\rm \mu$as $\sim$ 5.5 $R_{{\rm s}}$, which is broadly consistent with the peak brightness of the radio cores at 1-86 GHz located within $\sim 10^2$ $R_{{\rm s}}$. Our observations occurred in the middle of an enhancement in very-high-energy (VHE) $\rm \gamma$-ray flux, presumably originating in the vicinity of the central black hole. Our measurements, combined with results of multi-wavelength observations, favor a scenario in which the VHE region has an extended size of $\sim$20-60 $R_{{\rm s}}$.
We present the thermal evolution and emergent spectra of solidifying terrestrial planets along with the formation of steam atmospheres. The lifetime of a magma ocean and its spectra through a steam atmosphere depends on the orbital distance of the planet from the host star. For a type-I planet, which is formed beyond a certain critical distance from the host star, the thermal emission declines on a timescale shorter than approximately $10^6$ years. Therefore, young stars should be targets when searching for molten planets in this orbital region. In contrast, a type-II planet, which is formed inside the critical distance, will emit significant thermal radiation from near-infrared atmospheric windows during the entire lifetime of the magma ocean. The Ks and L bands will be favorable for future direct imaging because the planet-to-star contrasts of these bands are higher than approximately 10$^{-7}$-10$^{-8}$. Our model predicts that, in the type-II orbital region, molten planets would be present over the main sequence of the G-type host star if the initial bulk content of water exceeds approximately 1 wt%. In visible atmospheric windows, the contrasts of the thermal emission drop below $10^{-10}$ in less than $10^5$ years, whereas those of the reflected light remain $10^{-10}$ for both types of planets. Since the contrast level is comparable to those of reflected light from Earth-sized planets in the habitable zone, the visible reflected light from molten planets also provides a promising target for direct imaging with future ground- and space-based telescopes.
The unification model for powerful radio galaxies and radio-loud quasars postulates that these objects are intrinsically the same but viewed along different angles. Herschel Space Observatory data permit the assessment of that model in the far-infrared spectral window. We analyze photometry from Spitzer and Herschel for the distant 3CR hosts, and find that radio galaxies and quasars have different mid-infrared, but indistinguishable far-infrared colors. Both these properties, the former being orientation dependent and the latter orientation invariant, are in line with expectations from the unification model. Adding powerful radio-quiet active galaxies and typical massive star-forming galaxies to the analysis, we demonstrate that infrared colors not only provide an orientation indicator, but can also distinguish active from star-forming galaxies.
We report on a 350-ks NuSTAR observation of the magnetar 1E 1841-045 taken in 2013 September. During the observation, NuSTAR detected six bursts of short duration, with $T_{90}<1$ s. An elevated level of emission tail is detected after the brightest burst, persisting for $\sim$1 ks. The emission showed a power-law decay with a temporal index of 0.5 before returning to the persistent emission level. The long observation also provided detailed phase-resolved spectra of the persistent X-ray emission of the source. By comparing the persistent spectrum with that previously reported, we find that the source hard-band emission has been stable over approximately 10 years. The persistent hard X-ray emission is well fitted by a coronal outflow model, where $e^{+/-}$ pairs in the magnetosphere upscatter thermal X-rays. Our fit of phase-resolved spectra allowed us to estimate the angle between the rotational and magnetic dipole axes of the magnetar, $\alpha_{mag}=0.25$, the twisted magnetic flux, $2.5\times10^{26}\rm \ G\ cm^2$, and the power released in the twisted magnetosphere, $L_j=6\times10^{36}\rm \ erg\ s^{-1}$. Assuming this model for the hard X-ray spectrum, the soft X-ray component is well fit by a two-blackbody model, with the hotter blackbody consistent with the footprint of the twisted magnetic field lines on the star. We also report on the 3-year Swift monitoring observations obtained since 2011 July. The soft X-ray spectrum remained stable during this period, and the timing behavior was noisy, with large timing residuals.
The Faint Young Sun Paradox comes from the fact that solar luminosity (2-4)x10^9 years ago was insufficient to support the Earth's temperature necessary for the efficient development of geological and biological evolution (particularly, for the existence of considerable volumes of liquid water). It remains unclear by now if the so-called greenhouse effect on the Earth can resolve this problem. An interesting alternative explanation was put forward recently by M.Krizek (New Ast. 2012, 17, 1), who suggested that planetary orbits expand with time due to the local Hubble effect, caused by the uniformly-distributed Dark Energy. Then, under a reasonable value of the local Hubble constant, it is easy to explain why the Earth was receiving an approximately constant amount of solar irradiation for a long period in the past and will continue to do so for a quite long time in future.
We propose a new and robust method to test the consistency of the cosmic evolution given by a cosmological model. It is realized by comparing the combined quantity r_d^CMB/D_V^SN, which is derived from the comoving sound horizon r_d from cosmic microwave background (CMB) measurements and the effective distance D_V derived from low-redshift Type-Ia supernovae (SNe Ia) data, with direct and independent r_d/D_V obtained by baryon acoustic oscillation (BAO) measurements at median redshifts. We apply this test method for the LCDM and wCDM models, and investigate the consistency of the derived value of r_d/D_V from Planck 2015 and the SN Ia data sets of Union2.1 and JLA (z<1.5), and the r_d/D_V directly given by BAO data from six-degree-field galaxy survey (6dFGS), Sloan Digital Sky Survey Data Release 7 Main Galaxy Survey (SDSS-DR7 MGS), DR11 of SDSS-III, WiggleZ and Ly-alpha forecast surveys from Baryon Oscillation Spectroscopic Data (BOSS) DR-11 over 0.1<z<2.36. We find that the non-flat LCDM model is well consistent with the BAO measurements within 1-sigma CL across the whole redshift range (0<z<2.4), while the flat wCDM model is in tension with BAO data at 1.5 sigma CL at z=2.34 and z=2.36. Future surveys will further tight up the constraints significantly, and provide stronger test on the consistency.
Galaxy growth depends critically on the interplay between radiative cooling of cosmic gas and the resulting energetic feedback that cooling triggers. This interplay has proven exceedingly difficult to model, even with large supercomputer simulations, because of its complexity. Nevertheless, real galaxies are observed to obey simple scaling relations among their primary observable characteristics. Here we show that a generic emergent property of the interplay between cooling and feedback can explain the observed scaling relationships between a galaxy's stellar mass, its total mass, and its chemical enrichment level, as well as the relationship between the average orbital velocity of its stars and the mass of its central black hole. These relationships naturally result from any feedback mechanism that strongly heats a galaxy's circumgalactic gas in response to precipitation of colder clouds out of that gas, because feedback then suspends the gas in a marginally precipitating state.
We analyze simultaneously six composite RXTE spectra of GX~339-4 in the hard state comprising 77 million counts collected over 196 ks.The source spectra are ordered by luminosity and span the range 1.6% to 17% of the Eddington luminosity. Crucially, using our new tool PCACORR, we re-calibrate the data to a precision of 0.1%, an order of magnitude improvement over all earlier work. Using our advanced reflection model RELXILL, we target the strong features in the component of emission reflected from the disk, namely, the relativistically-broadened Fe K emission line, the Fe K edge and the Compton hump. We report results for two joint fits to the six spectra: For the first fit, we fix the spin parameter to its maximal value ($a_*=0.998$) and allow the inner disk radius $R_{in}$ to vary. Results include (i) precise measurements of $R_{in}$, with evidence that the disk becomes slightly truncated at a few percent of Eddington; and (ii) an order-of-magnitude swing with luminosity in the high energy cutoff, which reaches >890 keV at our lowest luminosity. For the second fit, implementing the standard assumption of $R_{in} = R_{ISCO}$, we estimate the spin parameter to be $a_* = 0.95^{+0.03}_{-0.05}$ (90% confidence, statistical). For both fits, and at the same level of statistical confidence, we estimate that the disk inclination is $i = 48 \pm 1$ deg and that the Fe abundance is super-solar, $A_{Fe} = 5 \pm 1$.
We present new high-resolution ($\sim$0\farcs09) $H$-band imaging observations of the circumstellar disk around the T Tauri star SU Aur. Our observations with Subaru-HiCIAO have revealed the presence of scattered light as close as 0\farcs15 ($\sim$20 AU) to the star. Within our image, we identify bright emission associated with a disk with a minimum radius of $\sim$90 AU, an inclination of $\sim$35$\degr$ from the plane of the sky, and an approximate P.A. of 15$\degr$ for the major axis. We find a brightness asymmetry between the northern and southern sides of the disk due to a non-axisymmetric disk structure. We also identify a pair of asymmetric tail structures extending east and west from the disk. The western tail extends at least 2\farcs5 (350 AU) from the star, and is probably associated with a reflection nebula previously observed at optical and near-IR wavelengths. The eastern tail extends at least 1\arcsec (140 AU) at the present signal-to-noise. These tails are likely due to an encounter with an unseen brown dwarf, but our results do not exclude the explanation that these tails are outflow cavities or jets.
The fundamental parameters describing the coronal spectrum of an accreting black hole are the slope $\Gamma$ of the power-law continuum and the energy $E_{cut}$ at which it rolls over. Remarkably, this parameter can be accurately measured for values as high as 1 MeV by modeling the spectrum of X-rays reflected from a black hole accretion disk at energies below 100 keV. This is possible because the details in the reflection spectrum, rich in fluorescent lines and other atomic features, are very sensitive to the spectral shape of the hardest coronal radiation illuminating the disk. We show that fitting simultaneous NuSTAR (3-79 keV) and low-energy (e.g., Suzaku) data with the most recent version of our reflection model RELXILL, one can obtain reasonable constraints on $E_{cut}$ at energies from tens of keV up to 1 MeV, for a source as faint as 1 mCrab in a 100 ks observation.
We study the atomic neutral hydrogen (HI) content of $\sim$1600 galaxies up
to $z \sim 0.1$ using stacking techniques. The observations were carried out
with the Westerbork Synthesis Radio Telescope (WSRT) in the area of the SDSS
South Galactic Cap (SSGC), where we selected a galaxy sample from the SDSS
spectroscopic catalog. Multi-wavelength information is provided by SDSS, NVSS,
GALEX, and WISE. We use the collected information to study HI trends with
color, star-forming, and active galactic nuclei (AGN) properties.
Using NUV-r colors, galaxies are divided into blue cloud, green valley and
red sequence galaxies. As expected based on previous observations, we detect HI
in green valley objects with lower amounts of HI than blue galaxies, while
stacking only produces a 3-$\sigma$ upper limit for red galaxies with M$_{\rm
HI}$ $<$ (5 $\pm$ 1.5) $\times$ 10$^{8}$ M$_{\odot}$ and M$_{\rm HI}/\rm{L}_r$
$<$ 0.02 $\pm$ 0.006 $\rm M_{\odot} / \rm L_{\odot} $. We find that the HI
content is more dependent on NUV-r color, and less on ionization properties, in
the sense that regardless of the presence of an optical AGN (based on optical
ionization line diagnostics), green-valley galaxies always show HI, whereas red
galaxies only produce an upper limit. This suggests that feedback from optical
AGN is not the (main) reason for depleting large-scale gas reservoirs.
Low-level radio continuum emission in our galaxies can stem either from star
formation, or from AGN. We use the WISE color-color plot to separate these
phenomena by dividing the sample into IR late-type and IR early-type galaxies.
We find that the radio emission in IR late-type galaxies stems from enhanced
star formation, and this group is detected in HI. However, IR early-type
galaxies lack any sign of HI gas and star formation activity, suggesting that
radio AGN are likely to be the source of radio emission in this population.
Aims. We model the present-day population of 'classical' low-mass X-ray binaries (LMXBs) with neutron star accretors, which have hydrogen-rich donor stars. Their population is compared with that of hydrogen-deficient LMXBs, known as ultracompact X-ray binaries (UCXBs). We model the observable LMXB population and compare it to observations. Methods. We combine the binary population synthesis code SeBa with detailed LMXB evolutionary tracks to model the size and properties of the present-day LMXB population in the Galactic Bulge. Whether sources are persistent or transient, and what their instantaneous X-ray luminosities are, is predicted using the thermal-viscous disk instability model. Results. We find a population of ~2.1 x 10^3 LMXBs with neutron star accretors. Of these about 15 - 40 are expected to be persistent (depending on model assumptions), with luminosities higher than 10^35 erg s^-1. About 7 - 20 transient sources are expected to be in outburst at any given time. Within a factor of two these numbers are consistent with the observed population of bright LMXBs in the Bulge. This gives credence to our prediction of the existence of a population of ~1.6 x 10^3 LMXBs with low donor masses that have gone through the period minimum, and have present-day mass transfer rates below 10^-11 Msun yr^-1. Conclusions. Even though the observed population of hydrogen-rich LMXBs in the Bulge is larger than the observed population of (hydrogen-deficient) UCXBs, the latter have a higher formation rate. While UCXBs may dominate the total LMXB population at the present, the majority would be very faint, or may have become detached and produced millisecond radio pulsars. In that case UCXBs would contribute significantly more to the formation of millisecond radio pulsars than hydrogen-rich LMXBs. [abridged]
The problem of origin and age of asteroid families studied very intensively. First of all youngest families are interesting due to possibility of the reconstruction collisional history. But in oldest families present objects with very close orbits some of them are products of recent collisions. The search for and studying of dynamics of such pairs of orbits is a main aim of the present paper. In result, two new young groups of minor planets are detected as well as some new members of known families.
Using hydrodynamical simulations, we show that an episode of star formation in the center of the Milky Way, with a star-formation-rate (SFR) $\sim 0.5$ M$_\odot$ yr$^{-1}$ for $\sim 30$ Myr, can produce bubbles that resemble the Fermi Bubbles (FBs), when viewed from the solar position. The morphology, extent and multi-wavelength observations of FBs, especially X-rays, constrain various physical parameters such as SFR, age, and the circum-galactic medium (CGM) density. We show that the interaction of the CGM with the Galactic wind driven by a star formation in the central region can explain the observed surface brightness and morphological features of X-rays associated with the Fermi Bubbles. Furthermore, assuming that cosmic ray electrons are accelerated {\it in situ} by shocks and/or turbulence, the brightness and morphology of gamma-ray emission and the microwave haze can be explained. The kinematics of the cold and warm clumps in our model also matches with recent observations of absorption lines through the bubbles.
Upcoming and current large astronomical survey experiments often seek to constrain cosmological parameters via measurements of subtle effects such as weak lensing, which can only be measured statistically. In these cases, instrumental effects in the image plane CCDs need to be accounted and/or corrected for in measurement algorithms. Otherwise, the systematic errors induced in the measurements might overwhelm the size of the desired effects. Lateral electric fields in the bulk of the CCDs caused by field shaping potentials or space charge build up as the electrons in the image are acquired can cause lateral deflections of the electrons drifting in the CCD bulk. Here, I report on the LSST effort to model these effects on a photon-by-photon basis by the use of a Monte Carlo technique. The eventual goal of this work is to produce a CCD model validated by laboratory data which can then be used to evaluate its effects on weak lensing science.
We show that emerging tension between the direct astronomical measurements at low redshifts and cosmological parameters deduced from the Planck measurements of the CMB anisotropies can be alleviated if the dark matter consists of two fractions, stable part being dominant and a smaller unstable fraction constituting of about 5 - 10 per cent at the recombination epoch. The latter decays after cosmological recombination but earlier than the present epoch.
Using the data from the public database of the All Sky Automated Survey ({\tt ASAS}) we performed the very first light curve analyses of the three eclipsing binary systems \astrobj{AO~Aqr}, \astrobj{CW~Aqr} and \astrobj{ASAS~012206-4924.7}. The physical parameters of the systems were determined by the {\tt PHOEBE} software. From an analysis of the ASAS data it was concluded that AO~Aqr was found to be a contact binary system while CW~Aqr and ASAS~012206-4924.7 were found to be near--contact and detached binaries, respectively. Finally, the locations of the components, corresponding to the estimated physical parameters, in the HR diagram were also discussed.
High-mass X-ray binaries are fundamental in the study of stellar evolution, nucleosynthesis, structure and evolution of galaxies and accretion processes. Hard X-rays observations by INTEGRAL and Swift have broad- ened significantly our understanding in particular for the super-giant systems in the Milky Way, which number has increased by almost a factor of three. INTEGRAL played a crucial role in the discovery, study and understanding of heavily obscured systems and of fast X-ray transients. Most super-giant systems can now be classified in three categories: classical/obscured, eccentric and fast transient. The classical systems feature low eccentricity and variability factor of about 1000, mostly driven by hydrodynamic phenomena occurring on scales larger than the accretion radius. Among them, systems with short orbital periods and close to Roche-Lobe overflow or with slow winds, appear highly obscured. In eccentric systems, the variability amplitude can reach even higher factors, because of the contrast of the wind density along the orbit. Four super-giant systems, featuring fast outbursts, very short orbital periods and anomalously low accretion rates, are not yet understood. Simulations of the accretion processes on relatively large scales have pro- gressed and reproduce parts of the observations. The combined effects of wind clumps, magnetic fields, neutron star rotation and eccentricity ought to be included in future modelling work. Observations with INTEGRAL in combination with other observatories were also important for detecting cyclotron resonant scattering features in spectra of X-ray pulsars, probing their variations and the geometry of the accretion column and emission regions. Finally, the unique characteristics of INTEGRAL and its long life time played a fundamental role for building a complete catalogue of HXMBs, to study the different populations of these systems ...
We compared the total mass density profiles of three different types of galaxies using weak gravitational lensing: (i) 29 galaxies that host quasars at z~0.32 that are in a post-starburst (PSQ) phase with high star formation indicating recent merger activity, (ii) 22 large elliptical galaxies from the SLACS sample that do not host a quasar at z~0.23, and (iii) 17 galaxies that host moderately luminous quasars at z~0.36 powered by disk instabilities, but with no intense star formation. On an initial test we found no evidence for a connection between the merger state of a galaxy and the profile of the halo, with the PSQ profile comparable to that of the other two samples and consistent with the Leauthaud et al. (2014) study of moderately luminous quasars in COSMOS. Given the compatibility of the two quasar samples, we combined these and found no evidence for any connection between black hole activity and the dark matter halo. All three mass profiles remained compatible with isothermality given the present data.
Using N-body simulations, we measure the power spectrum of the effective dark matter density field, which is defined through the modified Poisson equation in $f(R)$ cosmologies. We find that when compared to the conventional dark matter power spectrum, the effective power spectrum deviates more significantly from the $\Lambda$CDM model. For models with $f_{R0}=-10^{-4}$, the deviation can exceed 150\% while the deviation of the conventional matter power spectrum is less than 50\%. Even for models with $f_{R0}=-10^{-6}$, for which the conventional matter power spectrum is very close to the $\Lambda$CDM prediction, the effective power spectrum shows sizeable deviations. Our results indicate that traditional analyses based on the dark matter density field may seriously underestimate the impact of $f(R)$ gravity on galaxy clustering. We therefore suggest the use of the effective density field in such studies.
We study the cross-correlation of distribution of galaxies, the Sunyaev-Zel'dovich (SZ) and X-ray power spectra of galaxies from current and upcoming surveys and show these to be excellent probes of the nature, i.e. extent, evolution and energetics, of the circumgalactic medium (CGM). For a flat pressure profile, the SZ cross power spectrum shows oscillations at $l$-values corresponding to the length scales smaller than $\sim \frac{2}{3}$ times the virial radius of the galaxy. These oscillations are sensitive to the steepness of the pressure profile of the CGM and vanish for a sufficiently steep profile. Similar oscillations are also present in the X-ray cross power spectrum which is, however, more sensitive to the density profile. We forecast the detectability of the cross-correlated galaxy distribution, SZ and X-ray signals by combining SPT-DES and eROSITA-DES/eROSITA-LSST surveys, respectively. We find that, for the SPT-DES survey, the signal-to-noise ratio (SNR) peaks at high mass and redshift with SNR $\sim 9$ around $M_h\sim 10^{13} h^{-1} M_{\odot}$ and $z\sim 1.5\hbox{--} 2$ for flat density and temperature profiles. The SNR peaks at $\sim 6 (12 )$ for the eROSITA-DES (eROSITA-LSST) surveys. We also perform a Fisher matrix analysis and find that the gas fraction in the CGM can be constrained to a precision of $\sim 34\% (23 \%)$ by the SPT-DES and $\sim 23\% (14 \%)$ by the eROSITA-DES surveys in the presence (absence) of an unknown redshift evolution of the gas fraction. Finally, we demonstrate that the cross-correlated SZ-galaxy and X-ray-galaxy power spectrum can be used as powerful probes of the CGM energetics and potentially discriminate between different feedback models recently proposed in the literature; for example, one can distinguish a `no AGN feedback' scenario from a CGM energized by `fixed-velocity hot winds' at greater than $3\sigma$.
We select a sample of 34 gamma-ray bursts (GRBs) whose $\Gamma_0$ values are derived with the onset peak observed in the afterglow lightcurves (except for GRB 060218 whose $\Gamma_0$ is estimated with its radio data), and investigate the correlations among $\Gamma_0$, the isotropic peak luminosity ($L_{\rm p, iso}$), and the peak energy of the $\nu f_\nu$ spectrum in the cosmological rest frame ($E^{'}_{\rm p}$). An analysis of pair correlations among these observables well confirms the results reported by previous papers. In addition, a tight correlation among $L_{\rm p, iso}$, $E^{'}_{p}$, and $\Gamma_0$ is found from a multiple regression analysis, which takes the form of $L_{\rm p,iso} \propto {E'_{\rm p}}^{1.34\pm 0.14} \Gamma_0^{1.32\pm 0.19}$ or $E'_{\rm p} \propto L_{\rm p,iso}^{0.55\pm 0.06}\Gamma_0^{-0.50\pm 0.17}$. Nine other GRBs whose $\Gamma_0$ are derived via the pair production opacity constraint also follow such a correlation. We argue that this tighter $L_{\rm p, iso} - E^{'}_{p} - \Gamma_0$ correlation may be more physical than the $L_{\rm p,iso} - E'_{\rm p}$ correlation, since the relationship between the observed $L_{\rm p,iso}$ and $E'_{\rm p}$ not only depends on radiation physics, but also depends on the bulk motion of the jet. We explore possible origins of this correlation and discuss its physical implications for understanding GRB jet composition and radiation mechanism.
We use three domain wall simulations from the radiation era to the late time dark energy domination era based on the PRS algorithm to calculate the energy-momentum tensor components of domain wall networks in an expanding universe. Unequal time correlators in the radiation, matter and cosmological constant epochs are calculated using the scaling regime of each of the simulations. The CMB power spectrum of a network of domain walls is determined. The first ever quantitative constraint for the domain wall surface tension is obtained using a Markov chain Monte Carlo method; an energy scale of domain walls of 0.93 MeV, which is close but below the Zel'dovich bound, is determined.
The variability of accretion rate, which is believed to induce the aperiodic variability of X-ray emission from disk, may affect the energy injection into the jet. In this spirit, a correlation between disk emission and jet emission can be formed even if the mean luminosity of disk emission keeps constant. In this work, these correlations are found in the situation that the luminosity of disk emission is variable and kept with a constant mean value. The obtained correlations may be shallower than that of the fundamental plane of black hole activity. In addition, the slope of correlation may increase with increasing observed frequency of jet emission. For the luminosities spacing with three days, the slope of correlation decreases with increasing black hole mass. The deviation of our found correlations from that of the fundamental plane is related to the suppression of variability in the jet emission in comparison with that in the disk emission. This mechanism may work in some sources in which shallower correlations have been reported. Moreover, it implies that luminosities used to estimate the relation of fundamental plane should cover an appropriate timescale, in which the variability of jet emission is not significantly suppressed.
We estimate the degree of circular polarization for the gravitational waves generated during the electroweak and QCD phase transitions from the kinetic and magnetic helicity generated by bubble collisions during those cosmological phase transitions.
In this paper we report the discovery and detailed radio/X-ray analysis of a peculiar giant radio galaxy (GRG) detected by INTEGRAL, IGR J14488-4008. The source has been recently classified as a Seyfert 1.2 galaxy at redshift 0.123; the radio data denote the source to be a type II Fanaroff-Riley radio galaxy, with a linear projected size exceeding 1.5 Mpc, clearly assigning IGR J14488-4008 to the class of GRG. In the X-rays, the source shows a remarkable spectrum, characterised by absorption by ionised elements, a characteristic so far found in only other four broad line radio galaxies.
FUor outbursts in young stellar objects (YSOs) are the most dramatic events among episodic accretion phenomena. The origin of these bursts is not clear: disk instabilities and/or disk perturbations by an external body being the most viable hypotheses. Here, we report our VLT/SINFONI high angular resolution AO-assisted observations of 2MASS J06593158-0405277, which is undergoing a recently discovered FUor outburst. Our observations reveal the presence of an extended disc-like structure around the FUor, a very low-mass companion (2MASS J06593158-0405277B) at ~100 au in projection, and, possibly, a third closer companion at ~11 au. These sources appear to be young, displaying accretion signatures. Assuming the components are physically linked, 2MASS J06593158-0405277 would then be one of the very few triple systems observed in FUors.
We obtain a constraint on the parameters of a generalized cubic Galileon gravity model exhibiting the Vainshtein mechanism by using multi-wavelength observations of the Coma Cluster. The generalized cubic Galileon model is characterized by three parameters of the turning scale associated with the Vainshtein mechanism, and the amplitude of modifying a gravitational potential and a lensing potential. X-ray and Sunyaev-Zel'dovich (SZ) observations of the intra-cluster medium are sensitive to the gravitational potential, while the weak-lensing (WL) measurement is specified by the lensing potential. A joint fit of a complementary multi-wavelength dataset of X-ray, SZ and WL measurements enables us to simultaneously constrain these three parameters of the generalized cubic Galileon model for the first time.
We present a spectroscopic study of eight extremely low-metallicity candidate
emission-line galaxies with oxygen abundances possibly below 12 +log O/H =
7.35. These galaxies were selected from Data Release 10 of the Sloan Digital
Sky Survey/Baryon Oscillation Spectroscopic Survey (SDSS/BOSS DR10). We will
call these extremely metal-deficient galaxies XMD galaxies. The electron
temperature-sensitive emission line [O~{\sc iii}] $\lambda$4363 is detected in
three galaxies and marginally detected in two galaxies, allowing for abundance
determination by a "direct" method. Because of large uncertainties in the [O
{\sc iii}]$\lambda$4363\AA\ line fluxes, we also calculated oxygen abundance in
these galaxies together with the remaining three galaxies using a strong-line
semi-empirical method. This method gives oxygen abundances higher than 7.35 for
three galaxies with detected [O {\sc iii}]$\lambda$4363\AA\ line and lower than
7.35 for the remaining five objects of the sample. The newly-discovered
galaxies represent excellent targets for follow-up spectroscopic observations
with the largest telescopes to improve the oxygen abundance determination and
to increase the number of these very rare low-metallicity objects. The extreme
location of the most massive and luminous XMD galaxies and XMD candidates in
the stellar mass-metallicity diagram implies that these galaxies may be genuine
young objects.
With stellar masses of up to $\sim$ 10$^7$ - 10$^8$$M_{\odot}$, the galaxies
are not chemically enriched and strongly deviate to lower metallicity as
compared to the relation obtained for a large sample of low-redshift,
star-forming galaxies.
We investigate the dynamics of dust grains with various sizes in protoplanetary disk winds driven by magnetorotational turbulence, by simulating the time evolution of the dust grain distribution in the vertical direction. Small dust grains, which are well coupled to the gas, are dragged upward with the upflowing gas, while large grains remain near the midplane of a disk. Intermediate--size grains float at several scale heights from the midplane in time-averated force balance between the downward gravity and the upward gas drag. For the minimum mass solar nebula at 1 AU, dust grains with size of 20 -- 40 $\mu m$ float at 5-10 scale heights from the midplane. Considering the dependence on the distance from the central star, smaller-size grains remain only in an outer region of the disk, while larger-size grains are distributed in a broader region. This implies that the dust depletion is expected to take place in small-to-large and inside-out manners. We also discuss the implication of our result to the observation of dusty material around young stellar objects.
Extending previous studies of nonthermal electron transport in solar flares which include the effects of collisional energy diffusion and thermalization of fast electrons, we present an analytic method to infer more accurate estimates of the accelerated electron spectrum in solar flares from observations of the hard X-ray spectrum. Unlike for the standard cold-target model, the spatial characteristics of the flaring region, especially the necessity to consider a finite volume of hot plasma in the source, need to be taken into account in order to correctly obtain the injected electron spectrum from the source-integrated electron flux spectrum (a quantity straightforwardly obtained from hard X-ray observations). We show that, for a given source-integrated electron flux spectrum, the overall power in the injected electrons could be reduced by an order of magnitude or more relative to its cold-target value. Indeed, the extent of electron thermalization can be significant enough to nullify the need to introduce an {\it ad hoc} low-energy cutoff to the injected electron spectrum.
The Cybele region, located between the 2J:-1A and 5J:-3A mean-motion
resonances, is adjacent and exterior to the asteroid main belt. An increasing
density of three-body resonances makes the region between the Cybele and Hilda
populations dynamically unstable, so that the Cybele zone could be considered
the last outpost of an extended main belt. The presence of binary asteroids
with large primaries and small secondaries suggested that asteroid families
should be found in this region, but only relatively recently the first
dynamical groups were identified in this area. Among these, the Sylvia group
has been proposed to be one of the oldest families in the extended main belt.
In this work we identify families in the Cybele region in the context of the
local dynamics and non-gravitational forces such as the Yarkovsky and
stochastic YORP effects. We confirm the detection of the new Helga group at
$\simeq$3.65~AU, that could extend the outer boundary of the Cybele region up
to the 5J:-3A mean-motion resonance. We obtain age estimates for the four
families, Sylvia, Huberta, Ulla and Helga, currently detectable in the Cybele
region, using Monte Carlo methods that include the effects of stochastic YORP
and variability of the Solar luminosity. The Sylvia family should be $T = 1220
\pm 40$ Myr old, with a possible older secondary solution. Any collisional
Cybele group formed prior to the late heavy bombardment would have been most
likely completely dispersed in the jumping Jupiter scenario of planetary
migration.
The active galaxy NGC1275 lies at the centre of the Perseus cluster of galaxies, which is the X-ray brightest cluster in the Sky. The nucleus shows large variability over the past few decades. We compile a lightcurve of its X-ray emission covering about 40 years and show that the bright phase around 1980 explains why the inner X-ray bubbles were not seen in the images taken with the Einstein Observatory. The flux had dropped considerably by 1992 when images with the ROSAT HRI led to their discovery. The nucleus is showing a slow X-ray rise since the first Chandra images in 2000. If it brightens back to the pre-1990 level, then X-ray absorption spectroscopy by ASTRO-H can reveal the velocity structure of the shocked gas surrounding the inner bubbles.
SMSS J031300.36-670839.3 (hereafter SM0313-6708) is a sub-giant halo star, with no detectable Fe lines and large overabundances of C and Mg relative to Ca. We obtained VLT-UVES spectra extending to 3060 Angstroms showing strong OH A-X band lines enabling an oxygen abundance to be derived. The OH A-X band lines in SM0313-6708 are much stronger than the CH C-X band lines. Spectrum synthesis fits indicate an [O/C] ratio of 0.02 +- 0.175. Our high S/N UVES data also enabled us to lower the Fe abundance limit to [Fe/H]{3D},NLTE < -7.52 (3 sigma). These data support our previous suggestion that the star formed from the iron-poor ejecta of a single massive star Population III supernova.
All galaxies once passed through a hyperluminous quasar phase powered by accretion onto a supermassive black hole. But because these episodes are brief, quasars are rare objects typically separated by cosmological distances. In a survey for Lyman-alpha emission at redshift z ~ 2, we discovered a physical association of four quasars embedded in a giant nebula. Located within a substantial overdensity of galaxies, this system is probably the progenitor of a massive galaxy cluster. The chance probability of finding a quadruple quasar is estimated to be ~10^-7, implying a physical connection between Lyman-alpha nebulae and the locations of rare protoclusters. Our findings imply that the most massive structures in the distant universe have a tremendous supply (~ 10^11 solar masses) of cool dense (volume density ~1 cm^-3) gas, which is in conflict with current cosmological simulations.
We perform numerical simulations of impulsively generated Alfv\'en waves in an isolated photospheric flux tube, and explore the propagation of these waves along such magnetic structure that extends from the photosphere, where these waves are triggered, to the solar corona, and analyze resulting magnetic shells. Our model of the solar atmosphere is constructed by adopting the temperature distribution based on the semi-empirical model and specifying the curved magnetic field lines that constitute the magnetic flux tube which is rooted in the solar photosphere. The evolution of the solar atmosphere is described by 3D, ideal magnetohydrodynamic equations that are numerically solved by the FLASH code. Our numerical simulations reveal, based on the physical properties of the multi-shell magnetic twisters and the amount of energy and momentum associated with them, that these multi-shell magnetic twisters may be responsible for the observed heating of the lower solar corona and for the formation of solar wind. Moreover, it is likely that the existence of these twisters can be verified by high-resolution observations.
We discuss a new IRAS Faint Source Catalog galaxy redshift catalogue (RIFSCz) which incorporates data from Galex, SDSS, 2MASS, WISE, Akari and Planck. Akari fluxes are consistent with photometry from other far infrared and submillimetre missions provided an aperture correction is applied. Results from the Hermes-SWIRE survey in Lockman are also discussed briefly, and the strong contrast between the galaxy populations selected at 60 and 500 mu is summarized.
The Asteroid Redirect Mission (ARM) proposes to retrieve a near-Earth asteroid and position it in a lunar distant retrograde orbit (DRO) for later study, crewed exploration, and ultimately resource exploitation. During the Caltech Space Challenge, a recent workshop to design a crewed mission to a captured asteroid in a DRO, it became apparent that the asteroid's low escape velocity (<1 cm s$^{-1}$) would permit the escape of asteroid particles during any meaningful interaction with astronauts or robotic probes. This Note finds that up to 5% of escaped asteroid fragments will cross Earth-geosynchronous orbits and estimates the risk to satellites from particle escapes or complete disruption of a loosely bound rubble pile.
We spectroscopically measure multiple hydrogen Balmer line profiles from laboratory plasmas to investigate the theoretical line profiles used in white dwarf atmosphere models. X-ray radiation produced at the Z Pulsed Power Facility at Sandia National Laboratories initiates plasma formation in a hydrogen-filled gas cell, replicating white dwarf photospheric conditions. Here we present time-resolved measurements of H$\beta$ and fit this line using different theoretical line profiles to diagnose electron density, $n_{\rm e}$, and $n=2$ level population, $n_2$. Aided by synthetic tests, we characterize the validity of our diagnostic method for this experimental platform. During a single experiment, we infer a continuous range of electron densities increasing from $n_{\rm e}\sim4$ to $\sim30\times10^{16}\,$cm$^{-3}$ throughout a 120-ns evolution of our plasma. Also, we observe $n_2$ to be initially elevated with respect to local thermodynamic equilibrium (LTE); it then equilibrates within $\sim55\,$ns to become consistent with LTE. This supports our electron-temperature determination of $T_{\rm e}\sim1.3\,$eV ($\sim15,000\,$K) after this time. At $n_{\rm e}\gtrsim10^{17}\,$cm$^{-3}$, we find that computer-simulation-based line-profile calculations provide better fits (lower reduced $\chi^2$) than the line profiles currently used in the white dwarf astronomy community. The inferred conditions, however, are in good quantitative agreement. This work establishes an experimental foundation for the future investigation of relative shapes and strengths between different hydrogen Balmer lines.
We present the Gaia-Groundbased Observational Service for Asteroids (GOSA). Gaia-GOSA is an interactive tool which supports observers in planning photometric observations of asteroids. Each user is able to personalise the observation plan taking into account the equipment used and the observation site. The list of targets has been previously selected among the most relevant and scientifically remarkable objects, while the prediction of the transits in the Gaia's field of view have been calculated at the Observatoire de la C\^ote d'Azur. The data collected by the GOSA community will be exploited to enhance the reliability of the Gaia's Solar system science. The service is publicly available at www.gaiagosa.eu.
We present an analysis of the energy and frequency dependence of the Fourier time lags of the hectoHertz quasi-periodic oscillations (QPOs) and of the QPOs at the frequency at which the power density spectrum shows a break in the neutron-star low-mass X-ray binary 4U 1636-53, using a large data set obtained with the Rossi X-ray Timing Explorer. We found that: (i) For the break frequency QPO: for low frequencies, in general the time lag is positive, but it is decreasing with increasing frequency, reaching zero lag at 20 Hz. Between 20 and 35 Hz there is a small fluctuation around zero, from where the time lags become positive again and increase slightly above zero up to 65 Hz. (ii) For the hHz QPO: we see that when the frequency is 100 Hz the time lag is negative, but it increases to zero already at 110 Hz, being consistent with this value up to 130 Hz from where it increases to 0.5 msec at around 140 Hz. From 140 Hz the time lag decreases sharply, being strongly negative for hHz greater than 220 Hz. We compare these lags with our previous results for the lags of the kiloHertz QPOs in this same source and discuss possible scenarios for producing the lags in this system in the context of reflection off the accretion disc or (up-)down-scattering in a hot medium close to the neutron star.
We comprehensively compile and review N content in geologic materials to
calculate a new N budget for Earth. Using analyses of rocks and minerals in
conjunction with N-Ar geochemistry demonstrates that the Bulk Silicate Earth
(BSE) contains \sim7\pm4 times present atmospheric N (4\times10^18 kg N, PAN),
with 27\pm16\times10^18 kg N. Comparison to chondritic composition, after
subtracting N sequestered into the core, yields a consistent result, with BSE N
between 17\pm13\times10^18 kg to 31\pm24\times10^18 kg N. In the chondritic
comparison we calculate a N mass in Earth's core (180\pm110 to
300\pm180\times10^18 kg) and discuss the Moon as a proxy for the early mantle.
Significantly, we find the majority of the planetary budget of N is in the
solid Earth. The N estimate herein precludes the need for a "missing N"
reservoir. Nitrogen-Ar systematics in mantle rocks and basalts identify two
mantle reservoirs: MORB-source like (MSL) and high-N. High-N mantle is composed
of young, N-rich material subducted from the surface and is identified in OIB
and some xenoliths. In contrast, MSL appears to be made of old material, though
a component of subducted material is evident in this reservoir as well.
Using our new budget, we calculate a {\delta}15N value for BSE plus
atmosphere of \sim2\permil. This value should be used when discussing bulk
Earth N isotope evolution. Additionally, our work indicates that all surface N
could pass through the mantle over Earth history, and the mantle may act as a
long-term sink for N. Since N acts as a tracer of exchange between the
atmosphere, oceans, and mantle over time, clarifying its distribution in the
Earth is critical for evolutionary models concerned with Earth system
evolution. We suggest that N be viewed in the same vein as carbon: it has a
fast, biologically mediated cycle which connects it to a slow,
tectonically-controlled geologic cycle.
We revisit the kink-like parametrization of the deceleration parameter ($q(z)$) \cite{ishida08}, which considers a transition, at redshift $z_t$, from cosmic deceleration to acceleration. In this parametrization the initial ($z \gg z_t$) value of the q-parameter is $q_i$, its final ($z=-1$) value is $q_f$ and the duration of the transition is parametrized by $\tau$. We obtain constraints on the free parameters of the model using recent data from type Ia supernovae (SN Ia), baryon acoustic oscillations (BAO), cosmic microwave background (CMB) and the Hubble parameter (H(z)). The use of H(z) data introduces an explicit dependence of the combined likelihood on the present value of the Hubble parameter ($H_0$), allowing us to explore the influence of different priors when marginalizing over this parameter. We also study the importance of the CMB information in the results by considering data from WMAP7, WMAP9 (Wilkinson Microwave Anisotropy Probe - 7 and 9 years) and the Planck satellite. Assuming a flat space geometry, $q_i=1/2$ and expressing the present value of the deceleration parameter ($q_0$) as a function of the others three free parameters, we obtain $z_t=0.68^{+0.10}_{-0.09}$, $\tau=0.24^{+0.14}_{-0.10}$ and $q_0=-0.48^{+0.11}_{-0.14}$, at 68\% of confidence level, with flat prior on $H_0$. If in addition we fix $q_f=-1$, we get $z_t=0.67^{+0.03}_{-0.04}$, $\tau=0.32^{+0.04}_{-0.03}$ and $q_0=-0.55^{+0.05}_{-0.06}$.
The scale of Baryon Acoustic Oscillations (BAO) imprinted in the matter power spectrum provides an almost-perfect standard ruler: it only suffers sub-percent deviations from fixed comoving length due to non-linear effects. We study the BAO shift in the large Horndeski class of gravitational theories and compute its magnitude in momentum space using second-order perturbation theory and a peak-background split. The standard prediction is affected by the modified linear growth, as well as by non-linear gravitational effects that alter the mode-coupling kernel. For covariant Galileon models, we find a $14-45\%$ enhancement of the BAO shift with respect to standard gravity and a distinct time evolution depending on the parameters. Despite the larger values, the shift remains well below the forecasted precision of next-generation galaxy surveys. Models that produce significant BAO shift would cause large redshift-space distortions or affect the bispectrum considerably. Our computation therefore validates the use of the BAO scale as a comoving standard ruler for tests of general dark energy models.
We study the observable signatures of self-gravitating MHD turbulence by applying the probability density functions (PDFs) and the spatial density power spectrum to synthetic column density maps. We find that there exists three characterizable stages of the evolution of the collapsing cloud which we term "early," "intermediate," and "advanced." At early times, i.e. $t<0.15t_{ff}$, the column density has a power spectral slope similar to nongravitating supersonic turbulence and a lognormal distribution. At an intermediate stage, i.e. $0.15t_{ff}< t \leq 0.35t_{ff}$, there exists signatures of the prestellar cores in the shallower PDF and power spectrum power law slopes. The column density PDF power law tails at these times have line of sight averaged slopes ranging from -2.5 to -1.5 with shallower values belonging to simulations with lower magnetic field strength. The density power spectrum slope becomes shallow and can be characterized by $P(k)=A_1k^{\beta_2}e^{-k/k_c}$, where $A_1$ describes the amplitude, $k^{\beta_2}$ describes the classical power law behavior and the scale $k_c$ characterizes the turn over from turbulence dominated to self-gravity dominated. At advanced stages of collapse, i.e. $\approx t>0.35t_{ff}$, the power spectral slope is positive valued, and a dramatic increase is observed in the PDF moments and the Tsallis incremental PDF parameters, which gives rise to deviations between PDF-sonic Mach number relations. Finally, we show that the imprint of gravity on the density power spectrum can be replicated in non-gravitating turbulence by introducing a delta-function with amplitude equivalent to the maximum valued point in a given self-gravitating map. We find that the turbulence power spectrum restored through spatial filtering of the high density material.
The data analysis of the gravitational wave signals emitted by coalescing neutron star binaries requires the availability of an accurate analytical representation of the dynamics and waveforms of these systems. We propose an effective-one-body (EOB) model that describes the general relativistic dynamics of neutron star binaries from the early inspiral up to merger. Our EOB model incorporates an enhanced attractive tidal potential motivated by recent analytical advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. No fitting parameters are introduced for the description of tidal interaction in the late, strong-field dynamics. We compare the model energetics and the gravitational wave phasing with new high-resolution multi-orbit numerical relativity simulations of equal-mass configurations with different equations of state. We find agreement within the uncertainty of the numerical data for all configurations. Our model is the first semi-analytical model which captures the tidal amplification effects close to merger. It thereby provides the most accurate analytical representation of binary neutron star dynamics and waveforms currently available.
In this work we discuss the mechanism behind the large electroluminescence signals observed at relatively low electric fields in the holes of a Thick Gas Electron Multiplier (THGEM) electrode immersed in liquid xenon. We present strong evidence that the scintillation light is generated in xenon bubbles trapped below the THGEM holes. The process is shown to be remarkably stable over months of operation, providing - under specific thermodynamic conditions - energy resolution similar to that of present dual-phase liquid xenon experiments. The observed mechanism may serve as the basis for the development of Liquid Hole Multipliers (LHMs), capable of producing local charge-induced electroluminescence signals in large-volume single-phase noble-liquid detectors for dark matter and neutrino physics experiments.
Self-accelerating solutions in massive gravity provide explicit, calculable examples that exhibit the general interplay between superluminality, the well-posedness of the Cauchy problem, and strong coupling. For three particular classes of vacuum solutions, one of which is new to this work, we construct the conformal diagram for the characteristic surfaces on which isotropic stress-energy perturbations propagate. With one exception, all solutions necessarily possess spacelike characteristics, indicating perturbative superluminality. Foliating the spacetime with these surfaces gives a pathological frame where kinetic terms of the perturbations vanish, confusing the Hamiltonian counting of degrees of freedom. This frame dependence distinguishes the vanishing of kinetic terms from strong coupling of perturbations or an ill-posed Cauchy problem. We give examples where spacelike characteristics do and do not originate from a point where perturbation theory breaks down and where spacelike surfaces do or do not intersect all characteristics in the past light cone of a given observer. The global structure of spacetime also reveals issues that are unique to theories with two metrics: in all three classes of solutions, the Minkowski fiducial space fails to cover the entire de Sitter spacetime allowing worldlines of observers to end in finite proper time at determinant singularities. Characteristics run tangent to these surfaces requiring {\it ad hoc} rules to establish continuity across singularities.
Theories of modified gravity, in both the linear and fully non-linear regime, are often studied under the assumption that the evolution of the new (often scalar) degree of freedom present in the theory is quasi-static. This approximation significantly simplifies the study of the theory, and one often has good reason to believe that it should hold. Nevertheless it is a crucial assumption that should be explicitly checked whenever possible. In this paper we do so for the Vainshtein mechanism. By solving for the full spatial and time evolution of the Dvali-Gabadadze-Porrati and the Cubic Galileon model, in a spherical symmetric spacetime, we are able to demonstrate that the Vainshtein solution is a stable attractor and forms no matter what initial conditions we take for the scalar field. Furthermore,the quasi-static approximation is also found to be a very good approximation whenever it exists. For the best-fit Cubic Galileon model, however, we find that for deep voids at late times, the numerical solution blows up at the same time as the quasi-static solution ceases to exist. We argue that this phenomenon is a true instability of the model.
In this work we consider an inflationary potential in fractional form coupled non-minimally to gravity. This potential has a dominant constant energy density at early times which can realize successful inflation. We show that this potential predicts small tensor-to-scalar ratio of the order of $r\approx 0.01$ which is fully consistent with Planck constraints. Using the lower and upper bounds on reheating temperature, we provide additional constraints on the non-minimally coupling parameter $\xi$ of the model. We also study the preheating stage predicted by this kind of potentials using numerical calculations.
Dark energy is investigated from the perspective of quantum cosmology. By treating the existence of a classical universe as a constraint, it is found that the normal ordering ambiguity factor q in Wheeler-DeWitt equation tends to take its value on domain (-1, 3). Furthermore, to ensure the existence of a classical universe, there must be dark energy in the universe. It is in this sense we propose that dark energy is the reason for the existence of a classical universe.
It is customary to employ a semi-spherical scale model to describe the apparent path of the Sun across the sky, whether it be its diurnal motion or its variation throughout the year. A flat surface and three bent semi-rigid wires (representing the three solar arcs during solstices and equinoxes) will do the job. On the other hand, since very early times, there have been famous armillary spheres built and employed by the most outstanding astronomers for the description of the celestial movements. In those instruments, many of them now considered true works of art, Earth lies in the center of the cosmos and the observer looks at the whole "from the outside." Of course, both devices, the scale model of the sky and the armillary sphere, serve to represent the movement of the Sun, and in this paper we propose to show their equivalence by a simple construction. Knowing the basics underlying the operation of the armillary sphere will give us confidence to use it as a teaching resource in school.
The fast rates of magnetic reconnection found in both nature and experiments are important to understand theoretically. Recently, it was demonstrated that two-fluid magnetic reconnection remains fast in the strong guide field regime, regardless of the presence of fast-dispersive waves. This conclusion is in agreement with recent results from kinetic simulations, and is in contradiction to the findings in an earlier two-fluid study, where it was suggested that fast-dispersive waves are necessary for fast reconnection. In this paper, we give a more detailed derivation of the analytic model presented in a recent letter, and present additional simulation results to support the conclusions that the magnetic reconnection rate in this regime is independent of both collisional dissipation and system-size. In particular, we present a detailed comparison between fluid and kinetic simulations, finding good agreement in both the reconnection rate and overall length of the current layer. Finally, we revisit the earlier two-fluid study, which arrived at different conclusions, and suggest an alternative interpretation for the numerical results presented therein.
In a previous paper, we developed tools for studying the horizon geometry of a Kerr black hole that is tidally distorted by a binary compansion, using techniques that require large mass ratios but can be applied to any bound orbit and allow for arbitrary black hole spin. We now apply these techniques to generic Kerr black hole orbits. This allows us to investigate horizon dynamics: the tidal field perturbing the horizon's geometry varies over a generic orbit, with significant variations for eccentric orbits. We find that many of the features of the horizon's behavior found in our previous analysis carry over to the dynamical case in a natural way. In particular, we find significant offsets between the applied tide and the horizon's response, which leads to bulging in the horizon's geometry which can lag or lead the orbit, depending upon the hole's rotation and the orbit's geometry. An interesting new feature we find are small amplitude coherent wiggles in the horizon's response to the applied tide, which we explain using a teleological Green's function relating the shear to the tide. These wiggles only appear when the black hole's spin is very large ($a \ge 0.99M$). Future work which extends our abilities to embed the horizon to large spins will look for this feature in the distorted horizon's geometry.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)