In the cold dark matter paradigm, structures form hierarchically, implying that large structures contain smaller substructures. These subhalos will enhance signatures of dark matter annihilation such as gamma rays. In the literature, typical estimates of this boost factor are based on field-halo modelling, where halos are assumed to be virialized without any mass loss.However, since subhalos accreted in the gravitational potential of their host lose mass through tidal stripping and dynamical friction, they have a quite characteristic density profile, different from that of the field halos of the same mass. In this work, we quantify the effect of tidal stripping on the boost factor, by developing a semi-analytic model that combines mass-accretion history of both the host and subhalos as well as subhalo accretion rates. We find that compared with the field-halo models, the boost factor increases by a factor 2-3 for host halos ranging from sub-galaxy to cluster masses. The results are particularly relevant for indirect dark matter searches in the extragalactic gamma-ray sky.
We present Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) observations of the source and lens stars for planetary microlensing event OGLE-2005-BLG-169, which confirm the relative proper motion prediction due to the planetary light curve signal observed for this event. This (and the companion Keck result) provide the first confirmation of a planetary microlensing signal, for which the deviation was only 2%. The follow-up observations determine the flux of the planetary host star in multiple passbands and remove light curve model ambiguity caused by sparse sampling of part of the light curve. This leads to a precise determination of the properties of the OGLE-2005-BLG-169Lb planetary system. Combining the constraints from the microlensing light curve with the photometry and astrometry of the HST/WFC3 data, we find star and planet masses of M_* = 0.69+- 0.02 M_solar and m_p = 14.1 +- 0.9 M_earth. The planetary microlens system is located toward the Galactic bulge at a distance of D_L = 4.1 +- 0.4 kpc, and the projected star-planet separation is a_perp = 3.5 +- 0.3 AU, corresponding to a semi-major axis of a = 4.0 (+2.2 -0.6) AU.
Several long, dynamically cold stellar streams have been observed around the Milky Way Galaxy, presumably formed from the tidal disruption of globular clusters. In integrable potentials---where all orbits are dynamically regular---tidal debris phase-mixes close to the orbit of the progenitor system. However, cosmological simulations of structure formation suggest that the Milky Way's dark matter halo is expected not to be fully integrable; an appreciable fraction of orbits will be chaotic. This paper examines the influence of chaos on the phase-space morphology of cold tidal streams. We find very stark results: Streams in chaotic regions look very different from those in regular regions. We find that streams (simulated using test particle ensembles of nearby orbits) can be sensitive to chaos on a much shorter time-scale than any standard prediction (from the Lyapunov or frequency-diffusion times). For example, on a weakly chaotic orbit with a chaotic timescale predicted to be >1000 orbital periods (>1000 Gyr), the resulting stellar stream is, after just a few 10's of orbits, substantially more diffuse than any formed on a nearby but regular orbit. We find that the enhanced diffusion of the stream stars can be understood by looking at the variance in orbital frequencies of orbit ensembles centered around the parent (progenitor) orbit. Our results suggest that long, cold streams around our Galaxy must exist only on regular (or very nearly regular) orbits; they potentially provide a map of the regular regions of the Milky Way potential. This suggests a promising new direction for the use of tidal streams to constrain the distribution of dark matter around our Galaxy.
This document is an addendum to "One-point remapping of Lagrangian perturbation theory in the mildly non-linear regime of cosmic structure formation" (arXiv:1305.4642).
Axions currently provide the most compelling solution to the strong CP problem. These particles may be copiously produced in the early universe, including via thermal processes. Therefore, relic axions constitute a hot dark matter component and their masses are strongly degenerate with those of the three active neutrinos, as they leave identical signatures in the different cosmological observables. In addition, thermal axions, while still relativistic states, also contribute to the relativistic degrees of freedom, parameterised via $N_{eff}$. We present the cosmological bounds on the relic axion and neutrino masses, exploiting the full Planck mission data, which include polarization measurements. In the mixed hot dark matter scenario explored here, we find the tightest and more robust constraint to date on the sum of the three active neutrino masses, $\sum m_\nu <0.136$ eV at $95\%$ CL, obtained in the well-known linear perturbation regime. The Planck Sunyaev-Zeldovich cluster number count data further tightens this bound, providing a $95\%$ CL upper limit of $\sum m_\nu <0.126$ eV in this very same mixed hot dark matter model, a value which is very close to the expectations in the inverted hierarchical neutrino mass scenario. Using this same combination of data sets we find the most stringent bound to date on the thermal axion mass, $m_a<0.529$ eV at $95\%$ CL.
Numerical hydrodynamics calculations are performed to determine conditions under which giant planet eccentricities can be excited by parent gas disks. Unlike in other studies, Jupiter-mass planets are found to have their eccentricities amplified --- provided their orbits start eccentric. We disentangle the web of co-rotation, co-orbital, and external resonances to show that this finite-amplitude instability is consistent with that predicted analytically. Ellipticities can grow until they reach of order the disk's aspect ratio, beyond which the external Lindblad resonances that excite eccentricity are weakened by the planet's increasingly supersonic epicyclic motion. Forcing the planet to still larger eccentricities causes catastrophic eccentricity damping as the planet collides into gap walls. If the eccentricity driving documented here survives in 3D, it may explain the low-to-moderate eccentricities $\lesssim 0.1$ exhibited by many giant planets (including Jupiter and Saturn), especially those without planetary or stellar companions.
We present custom-processed UV, optical, and near-IR photometry for the RESOLVE survey, a volume-limited census of stellar, gas, and dynamical mass within two subvolumes of the nearby universe (RESOLVE-A and -B), complete down to baryonic mass ~10^9.1-9.3 Msun. In contrast to standard pipeline photometry (e.g., SDSS), our photometry uses optimal background subtraction, avoids suppressing color gradients, and includes systematic errors. With these improvements, we measure brighter magnitudes, larger radii, bluer colors, and a real increase in scatter around the red sequence. Combining stellar masses from our photometry with the RESOLVE-A HI mass census, we create volume-limited calibrations of the photometric gas fractions (PGF) technique, which predicts gas-to-stellar mass ratios (G/S) from galaxy colors and optional additional parameters. We analyze G/S-color residuals vs. potential third parameters, finding that axial ratio is the best independent and physically meaningful third parameter. We define a "modified color" from planar fits to G/S as a function of both color and axial ratio. In the complete galaxy population, upper limits on G/S bias linear and planar fits. We therefore model the entire PGF probability density field, enabling iterative statistical modeling of upper limits and prediction of full G/S probability distributions for individual galaxies, with two-component structure for red colors. We use the RESOLVE-B 21cm census to test several PGF calibrations, finding that most systematically under- or overestimate gas masses, but the full probability density method performs well.
Gamma-Ray Bursts (GRBs) are the strongest explosions in the Universe, which due to their extreme character likely involve some of the strongest magnetic fields in nature. This review discusses the possible roles of magnetic fields in GRBs, from their central engines, through the launching, acceleration and collimation of their ultra-relativistic jets, to the dissipation and particle acceleration that power their $\gamma$-ray emission, and the powerful blast wave they drive into the surrounding medium that generates their long-lived afterglow emission. An emphasis is put on particular areas in which there have been interesting developments in recent years.
We present the characterization and initial results from the QUEST-La Silla AGN variability survey. This is an effort to obtain well sampled optical light curves in extragalactic fields with unique multi-wavelength observations. We present photometry obtained from 2010 to 2012 in the XMM-COSMOS field, which was observed over 150 nights using the QUEST camera on the ESO-Schmidt telescope. The survey uses a broadband filter, the $Q$-band, similar to the union of the $g$ and the $r$ filters, achieving an intrinsic photometric dispersion of $0.05$ mag, and a systematic error of $0.05$ mag in the zero-point. Since some detectors of the camera show significant non-linearity, we use a linear correlation to fit the zero-points as a function of the instrumental magnitudes, thus obtaining a good correction to the non-linear behavior of these detectors. We obtain good photometry to an equivalent limiting magnitude of $r\sim 20.5$. Studying the optical variability of X-ray detected sources in the XMM-COSMOS field, we find that the survey is $\sim75-80$% complete to magnitudes $r\sim20$, and $\sim67$% complete to a magnitude $r\sim21$. The determination and parameterization of the structure function (${SF}_{norm}(\tau) = A \tau^{\gamma}$) of the variable sources shows that most BL AGN are characterized by $A > 0.1$ and $\gamma > 0.025$. It is further shown that variable NL AGN and GAL sources occupying the same parameter space in $A$ and $\gamma$ are very likely to correspond to obscured or low luminosity AGN. Our samples are, however, small, and we expect to revisit these results using larger samples with longer light curves obtained as part of our ongoing survey.
We present here a quantitative analysis of the recent AMS-02 data with the purpose of investigating the interplay between astrophysical sources and Dark Matter in their interpretation. First, we show that AMS-02 leptonic measurements are in a remarkably good agreement with the hypothesis that all electrons and positrons are the outcome of primary or secondary astrophysical processes. Then, we add Dark Matter to the picture, in order to establish which are the informations on its annihilation cross section (or lifetime) that can be inferred by fitting AMS-02 data within a scenario in which Dark Matter and astrophysical sources jointly contribute to the different leptonic observables. In particular, by performing a Markov Chain Monte Carlo sampling of the parameters space of the theory, we attempt at characterizing the significance of a possible Dark Matter contribution to the observed data and we derive robust upper limits on the Dark Matter annihilation/decay rate.
Molecular absorption lines of OH (99 lines) and CH (105 lines) are measured for the carbon-enhanced metal-poor star BD+44 493 with [Fe/H]=-3.8. The abundances of oxygen and carbon determined from individual lines based on an 1D-LTE analysis exhibit significant dependence on excitation potentials of the lines; d log e/d chi ~ -0.15 - -0.2 dex/eV, where e and chi are elemental abundances from individual spectral lines and their excitation potentials, respectively. The dependence is not explained by the uncertainties of stellar parameters, but suggests that the atmosphere of this object possesses a cool layer that is not reproduced by the 1D model atmosphere. This result agrees with the predictions by 3D model calculations. Although absorption lines of neutral iron exhibit similar trend, it is much weaker than found in molecular lines and that predicted by 3D LTE models.
We model the generation of a magnetic field in a protostellar disc using an \alpha-dynamo and perform axisymmetric magnetohydrodynamics (MHD) simulations of a T Tauri star. We find that for small values of the dimensionless dynamo parameter $\alpha_d$ the poloidal field grows exponentially at a rate ${\sigma} \propto {\Omega}_K \sqrt{\alpha_d}$ , before saturating to a value $\propto \sqrt{\alpha_d}$ . The dynamo excites dipole and octupole modes, but quadrupole modes are suppressed, because of the symmetries of the seed field. Initial seed fields too weak to launch MHD outflows are found to grow sufficiently to launch winds with observationally relevant mass fluxes of order $10^{-9} M_{\odot}/\rm{yr}$ for T Tauri stars. For large values of $\alpha_d$ magnetic loops are generated over the entire disc. These quickly come to dominate the disc dynamics and cause the disc to break up due to the magnetic pressure.
We present a Bayesian reconstruction method which maps a galaxy distribution from redshift-space to real-space inferring the distances of the individual galaxies. The method is based on sampling density fields assuming a lognormal prior with a likelihood given by the negative binomial distribution function modelling stochastic bias. We assume a deterministic bias given by a power law relating the dark matter density field to the expected halo or galaxy field. Coherent redshift-space distortions are corrected in a Gibbs-sampling procedure by moving the galaxies from redshift-space to real-space according to the peculiar motions derived from the recovered density field using linear theory with the option to include tidal field corrections from second order Lagrangian perturbation theory. The virialised distortions are corrected by sampling candidate real-space positions (being in the neighbourhood of the observations along the line of sight), which are compatible with the bulk flow corrected redshift-space position adding a random dispersion term in high density collapsed regions. The latter are defined according to the eigenvalues of the Hessian. This approach presents an alternative method to estimate the distances to galaxies using the three dimensional spatial information, and assuming isotropy. Hence the number of applications is very broad. In this work we show the potential of this method to constrain the growth rate up to $k$ ~ 0.3 $h$ Mpc$^{-1}$. Furthermore it could be useful to correct for photo-metric redshift errors, and to obtain improved BAO reconstructions.
We analyze high resolution simulations of compressible, MHD turbulence with properties resembling conditions in galaxy clusters. The flow is driven to turbulence Mach number $\mathcal{M}_t \sim 1/2$ in an isothermal medium with an initially very weak, uniform seed magnetic field ($\beta = P_g/P_B = 10^6$). Since cluster turbulence is likely to result from a mix of sheared (solenoidal) and compressive forcing processes, we examine the distinct turbulence properties for both cases. In one set of simulations velocity forcing is entirely solenoidal ($\nabla\cdot \delta {\vec u} = 0$), while in the other it is entirely compressive ($\nabla\times \delta {\vec u} = 0$). Both cases develop a mixture of solenoidal and compressive turbulent motions, since each generates the other. The development of compressive turbulent motions leads to shocks, even when the turbulence is solenoidally forced and subsonic. Shocks, in turn, produce and amplify vorticity, which is especially important in compressively forced turbulence. To clarify those processes we include a pair of appendices that look in detail at vorticity evolution in association with shocks. From our simulation analyses we find that magnetic fields amplified to near saturation levels in predominantly solenoidal turbulence can actually enhance vorticity on small scales by concentrating and stabilizing shear. The properties, evolution rates and relative contributions of the kinetic and magnetic turbulent elements depend strongly on the character of the forcing. Specifically, shocks are stronger, but vorticity evolution and magnetic field amplification are slower and weaker when the turbulence is compressively forced. We identify a simple relation to estimate characteristic shock strengths in terms of the turbulence Mach number and the character of the forcing. Our results will be helpful in understanding flow motions in galaxy clusters.
Aurorae are detected from all the magnetized planets in our Solar System, including Earth. They are powered by magnetospheric current systems that lead to the precipitation of energetic electrons into the high-latitude regions of the upper atmosphere. In the case of the gas-giant planets, these aurorae include highly polarized radio emission at kilohertz and megahertz frequencies produced by the precipitating electrons, as well as continuum and line emission in the infrared, optical, ultraviolet and X-ray parts of the spectrum, associated with the collisional excitation and heating of the hydrogen-dominated atmosphere. Here we report simultaneous radio and optical spectroscopic observations of an object at the end of the stellar main sequence, located right at the boundary between stars and brown dwarfs, from which we have detected radio and optical auroral emissions both powered by magnetospheric currents. Whereas the magnetic activity of stars like our Sun is powered by processes that occur in their lower atmospheres, these aurorae are powered by processes originating much further out in the magnetosphere of the dwarf star that couple energy into the lower atmosphere. The dissipated power is at least four orders of magnitude larger than what is produced in the Jovian magnetosphere, revealing aurorae to be a potentially ubiquitous signature of large-scale magnetospheres that can scale to luminosities far greater than those observed in our Solar System. These magnetospheric current systems may also play a part in powering some of the weather phenomena reported on brown dwarfs.
Galactic Center (GC) molecular cloud Sgr B2 is the best manifestation of an X-ray reflection nebula (XRN) reprocessing a past giant outburst from the supermassive black hole Sgr A*. Alternatively, Sgr B2 could be illuminated by low-energy cosmic ray electrons (LECRe) or protons (LECRp). In 2013, NuSTAR for the first time resolved Sgr B2 hard X-ray emission on sub-arcminute scales. Two prominent features are detected above 10 keV - a newly emerging cloud G0.66-0.13 and the central 90" radius region containing two compact cores Sgr B2(M) and Sgr B2(N) surrounded by diffuse emission. It is inconclusive whether the remaining level of Sgr B2 emission is still decreasing or has reached a constant background level. A decreasing Fe K$\alpha$ emission can be best explained by XRN while a constant background emission can be best explained by LECRp. In the XRN scenario, the 3-79 keV Sgr B2 spectrum can well constrain the past Sgr A* outburst, resulting in an outburst spectrum with a peak luminosity of $L_{3-79\rm~keV} \sim 5\times10^{38} \rm~erg~s^{-1}$ derived from the maximum Compton-scattered continuum and the Fe K$\alpha$ emission consistently. The XRN scenario is preferred by the fast variability of G0.66-0.13, which could be a molecular clump located in the Sgr B2 envelope reflecting the same Sgr A* outburst. In the LECRp scenario, we derived the required CR ion power $dW/dt=(1-4)\times10^{39}\rm~erg~s^{-1}$ and the CR ionization rate $\zeta_{H}=(6-10)\times 10^{-15}\rm~H^{-1}~s^{-1}$. The Sgr B2 background level X-ray emission will be a powerful tool to constrain GC CR population.
Decades-long repeat observations of supernova SN1987A offer us unique, real-time insights into the violent death of a massive star and its long-term environmental effects, until its eventual switch-off.
We present the first search for spinning dust emission from a sample of 34 Galactic cold cores, performed using the CARMA interferometer. For each of our cores we use photometric data from the Herschel Space Observatory to constrain N_{H}, T_{d}, n_{H}, and G_{0}. By computing the mass of the cores and comparing it to the Bonnor-Ebert mass, we determined that 29 of the 34 cores are gravitationally unstable and undergoing collapse. In fact, we found that 6 cores are associated with at least one young stellar object, suggestive of their proto-stellar nature. By investigating the physical conditions within each core, we can shed light on the cm emission revealed (or not) by our CARMA observations. Indeed, we find that only 3 of our cores have any significant detectable cm emission. Using a spinning dust model, we predict the expected level of spinning dust emission in each core and find that for all 34 cores, the predicted level of emission is larger than the observed cm emission constrained by the CARMA observations. Moreover, even in the cores for which we do detect cm emission, we cannot, at this stage, discriminate between free-free emission from young stellar objects and spinning dust emission. We emphasise that, although the CARMA observations described in this analysis place important constraints on the presence of spinning dust in cold, dense environments, the source sample targeted by these observations is not statistically representative of the entire population of Galactic cores.
We use 350 mu angular diameter estimates from Planck to test the idea that some galaxies contain exceptionally cold (10-13 K) dust, since colder dust implies a lower surface brightness radiation field illuminating the dust, and hence a greater physical extent for a given luminosity. The galaxies identified from their spectral energy distributions as containing cold dust do indeed show the expected larger 350 mu diameters. For a few cold dust galaxies where Herschel data are available we are able to use submillimetre maps or surface brightness profiles to locate the cold dust, which as expected generally lies outside the optical galaxy.
We have studied the filaments extracted from the column density maps of the nearby Lupus 1, 3, and 4 molecular clouds, derived from photometric maps observed with the Herschel satellite. Filaments in the Lupus clouds have quite low column densities, with a median value of $\sim$1.5$\times$10$^{21}$ cm$^{-2}$ and most have masses per unit length lower than the maximum critical value for radial gravitational collapse. Indeed, no evidence of filament contraction has been seen in the gas kinematics. We find that some filaments, that on average are thermally subcritical, contain dense cores that may eventually form stars. This is an indication that in the low column density regime, the critical condition for the formation of stars may be reached only locally and this condition is not a global property of the filament. Finally, in Lupus we find multiple observational evidences of the key role that the magnetic field plays in forming filaments, and determining their confinement and dynamical evolution.
We present infrared multi-epoch observations of the dust forming nova V1280 Sco over $\sim$2000 days from the outburst. The temporal evolution of the infrared spectral energy distributions at 1272, 1616 and 1947 days can be explained by the emissions produced by amorphous carbon dust of mass (6.6--8.7)$\times$10$^{-8}$M$_{\odot}$ with a representative grain size of 0.01$~\mu$m and astronomical silicate dust of mass (3.4--4.3)$\times$10$^{-7}$M$_{\odot}$ with a representative grain size of 0.3--0.5$~\mu$m. Both of these dust species travel farther away from the white dwarf without an apparent mass evolution throughout those later epochs. The dust formation scenario around V1280 Sco suggested from our analyses is that the amorphous carbon dust is formed in the nova ejecta followed by the formation of silicate dust in the expanding nova ejecta or as a result of the interaction between the nova wind and the circumstellar medium.
NGC 6397 is one of the most interesting, well observed and theoretically studied globular clusters. The existing wealth of observations allows us to study the reliability of the theoretical white dwarf cooling sequences of low metallicity progenitors,to determine its age and the percentage of unresolved binaries, and to assess other important characteristics of the cluster, like the slope of the initial mass function, or the fraction of white dwarfs with hydrogen deficient atmospheres. We present a population synthesis study of the white dwarf population of NGC 6397. In particular, we study the shape of the color-magnitude diagram, and the corresponding magnitude and color distributions. We do this using an up-to-date Monte Carlo code that incorporates the most recent and reliable cooling sequences and an accurate modeling of the observational biases. We find a good agreement between our theoretical models and the observed data. In particular, we find that this agreement is best for those cooling sequences that take into account residual hydrogen burning. This result has important consequences for the evolution of progenitor stars during the thermally-pulsing asymptotic giant branch phase, since it implies that appreciable third dredge-up in low-mass, low-metallicity progenitors is not expected to occur. Using a standard burst duration of 1.0 Gyr, we obtain that the age of the cluster is 12.8+0.50-0.75 Gyr. Larger ages are also compatible with the observed data, but then realistic longer durations of the initial burst of star formation are needed to fit the luminosity function. We conclude that a correct modeling of the white dwarf opulation of globular clusters, used in combination with the number counts of main sequence stars provides an unique tool to model the properties of globular clusters.
Most of super-Earths detected by the radial velocity (RV) method have significantly smaller eccentricities than the eccentricities corresponding to velocity dispersion equal to their surface escape velocity ("escape eccentricities"). If orbital instability followed by giant impacts among protoplanets that have migrated from outer region is considered, it is usually considered that eccentricities of the merged bodies become comparable to those of orbital crossing bodies, which are excited up to their escape eccentricities by close scattering. However, the eccentricity evolution in the {\it in situ} accretion model has not been studied in detail. Here, we investigate the eccentricity evolution through {\it N}-body simulations. We have found that the merged planets tend to have much smaller eccentricities than the escape eccentricities due to very efficient collision damping. If the protoplanet orbits are initially well separated and their eccentricities are securely increased, an inner protoplanet collides at its apocenter with an outer protoplanet at its pericenter. The eccentricity of the merged body is the smallest for such configuration. Orbital inclinations are also damped by this mechanism and planets tend to share a same orbital plane, which is consistent with {\it Kepler} data. Such efficient collision damping is not found when we start calculations from densely packed orbits of the protoplanets. If the protoplanets are initially in the mean-motion resonances, which corresponds to well separated orbits, the {\it in situ} accretion model well reproduces the features of eccentricities and inclinations of multiple super-Earths/Earth systems discovered by RV and {\it Kepler} surveys.
We study the host galaxy properties of the tidal disruption object, Swift J164449.3+573451 using long-term optical to near-infrared (NIR) data. First, we decompose the galaxy surface brightness distribution and analyze the morphology of the host galaxy using high resolution \emph{HST} WFC3 images. We conclude that the host galaxy is a bulge-dominant galaxy that is well described by a single S\'{e}rsic model with S\'{e}rsic index $n=3.43\pm0.05$. Adding a disk component, the bulge to total host galaxy flux ratio (B/T) is $0.83\pm0.03$, which still indicates a bulge-dominant galaxy. Second, we estimate multi-band fluxes of the host galaxy through long-term light curves. Our long-term NIR light curves reveal the pure host galaxy fluxes $\sim500$ days after the burst. We fit spectral energy distribution (SED) models to the multi-band fluxes from the optical to NIR of the host galaxy and determine its properties. The stellar mass, the star formation rate, and the age of stellar population are $\log(M_{\star}/M_{\odot}) = 9.14^{+0.13}_{-0.10}$, $0.03^{+0.28}_{-0.03}\, M_{\odot}$/yr, and $0.63^{+0.95}_{-0.43}$ Gyr. Finally, we estimate the mass of the central super massive black hole which is responsible for the tidal disruption event. The black hole mass is estimated to be $10^{6.7\pm0.4}\, M_{\odot}$ from $M_{\mathrm{BH}}$ - $M_{\star,\mathrm{bul}}$ and $M_{\mathrm{BH}}$ - $L_{\mathrm{bul}}$ relations for the $K$ band, although a smaller value of $\sim10^5\, M_{\odot}$ cannot be excluded convincingly if the host galaxy harbors a pseudobulge.
We present deep near-infrared (NIR) J, Ks photometry of the old, metal-poor Galactic globular cluster M\,15 obtained with images collected with the LUCI1 and PISCES cameras available at the Large Binocular Telescope (LBT). We show how the use of First Light Adaptive Optics system coupled with the (FLAO) PISCES camera allows us to improve the limiting magnitude by ~2 mag in Ks. By analyzing archival HST data, we demonstrate that the quality of the LBT/PISCES color magnitude diagram is fully comparable with analogous space-based data. The smaller field of view is balanced by the shorter exposure time required to reach a similar photometric limit. We investigated the absolute age of M\,15 by means of two methods: i) by determining the age from the position of the main sequence turn-off; and ii) by the magnitude difference between the MSTO and the well-defined knee detected along the faint portion of the MS. We derive consistent values of the absolute age of M15, that is 12.9+-2.6 Gyr and 13.3+-1.1 Gyr, respectively.
We present the current accounting of systematic effect uncertainties for the Low Frequency Instrument (LFI) that are relevant to the 2015 release of the Planck cosmological results, showing the robustness and consistency of our data set, especially for polarization analysis. We use two complementary approaches: (i) simulations based on measured data and physical models of the known systematic effects; and (ii) analysis of difference maps containing the same sky signal ("null-maps"). The LFI temperature data are limited by instrumental noise. At large angular scales the systematic effects are below the cosmic microwave background (CMB) temperature power spectrum by several orders of magnitude. In polarization the systematic uncertainties are dominated by calibration uncertainties and compete with the CMB $E$-modes in the multipole range 10-20. Based on our model of all known systematic effects, we show that these effects introduce a slight bias of around $0.2\,\sigma$ on the reionization optical depth derived from the 70 GHz $EE$ spectrum using the 30 and 353\,GHz channels as foreground templates. At 30 GHz the systematic effects are smaller than the Galactic foreground at all scales in temperature and polarization, which allows us to consider this channel as a reliable template of synchrotron emission. We assess the residual uncertainties due to LFI effects on CMB maps and power spectra after component separation and show that these effects are smaller than the CMB at all scales. We also assess the impact on non-Gaussianity studies and find it to be negligible. Some residuals still appear in null maps from particular sky survey pairs, particularly at 30 GHz, suggesting possible straylight contamination due to an imperfect knowledge of the beam far sidelobes.
Context. The interaction of solar oscillations with near surface convection
is poorly understood. These interactions are likely the cause of several
problems in helio- and astero-seismology, including the so-called surface
effect and apparently unphysical travel time shifts as a function of center to
limb distance. There is thus a clear need for further theoretical understanding
and observational tests.
Aims. The aim is to determine how the observed modes are affected by the
convection.
Methods. I use HMI velocity and intensity images to construct k-$\omega$
diagrams showing how the oscillation amplitude and phase depend on the local
granulation intensity.
Results. There is a clear and significant dependence of the observed
properties of the oscillations on the local convection state.
We report the novel detection of complex high-column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (Av) of ~6-14. It has a slope of alpha=1.3-2 for the rho~r^-alpha profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above Av~40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with alpha>2. It correlates with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 10^4 cm^-3. This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: 1. rotation, which introduces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.
SDO/EVE provide rich information of the thermodynamic processes of solar activities, particularly of solar flares. Here, we develop a method to construct thermodynamic spectrum (TDS) charts based on the EVE spectral lines. Reading from the charts, we are able to easily recognize if there is a late phase following a main phase of a flare, and able to learn the begin, peak and end times of the flare as well as the drift of the temperature, i.e., the cooling rate, of the heated plasma during the flare. Through four M-class flares of different types, we illustrate which thermodynamic information can be revealed from the TDS charts. Further, we investigate the TDS charts of all the flares greater than M5.0, and some interesting results are achieved. First, there are two distinct drift patterns, called Type I and Type II. For Type I flares, the enhanced emission drifts from high to low temperture, whereas for Type II flares, the drift is somewhat reversed, suggesting a more violent and durable heating during Type II flares than Type I flares. Second, for late-phase flares, the peak intensity ratio of the late phase to the main phase read from the TDS chart is roughly correlated with the flare class identified by GOES SXR, and the flares with a stronger late phase are all confined. We believe that the re-deposition of the energy carried by a flux rope, that unsuccessfully erupts out, into thermal emissions is responsible for the strong late phase found in a confined flare. These results provide us new clues to advance our understanding of the thermodynamic processes of solar flares and associated solar eruptions, e.g., coronal mass ejections.
We report a new interpretation of the millimeter-wave polarization of the protoplanetary disk around HL Tau with self-polarization. We successfully reproduce the observed polarization signature with self-scattered light of dust grains. The detected polarization can be explained only if dust grains have a maximum size of around 150 ${\rm \mu m}$. This is a strong constraint on grain size in the early stage of a circumstellar disk. The obtained grain size contradicts to previously expected grain size, which is millimeter. The inferred grain size is too small to be trapped at gas pressure bumps, and therefore it requires planet formation theory to explain the mechanism to stop the grain growth or it might suggest the dust grains are no longer spherical but highly porous.
We present Keck NIRC2 high angular resolution adaptive optics observations of the microlensing event OGLE-2005-BLG-169, taken 8.21 years after the discovery of this planetary system. For the first time for a microlensing planetary event, the source and the lens are completely resolved, providing a precise measurement of their heliocentric relative proper motion, $\mu_{\rm{rel},\rm{helio}}=7.44 \pm 0.17$ mas yr$^{-1}$. This confirms and refines the initial model presented in the discovery paper and rules out a range of solutions that were allowed by the microlensing light curve. This is also the first time that parameters derived from a microlensing planetary signal are confirmed, both with Keck measurements, presented in this paper, and independent measurements obtained with the Hubble Space Telescope in I, V and B bands, presented in a companion paper. Hence, this new measurement of $\mu_{\rm{rel},\rm{helio}}$, as well as the measured brightness of the lens in H band, enabled the mass and distance of the system to be updated: a Uranus-mass planet ($m_\rm{p}=13.2\pm 1.3 M_\oplus$) orbiting a K5-type main sequence star ($M_*=0.65\pm 0.05 M_\odot$) separated by $a_\perp=3.4\pm 0.3$ AU, at the distance $D_\rm{L}=4.0\pm 0.4$ kpc from us.
Air showers induced by cosmic rays create nanosecond pulses detectable at radio frequencies. These pulses have been measured successfully in the past few years at the LOw- Frequency ARray (LOFAR) and are used to study the properties of cosmic rays. For a complete understanding of this phenomenon and the underlying physical processes, an absolute calibration of the detecting antenna system is needed. We present three approaches that were used to check and improve the antenna model of LOFAR and to provide an absolute calibration for air shower measurements. Two methods are based on calibrated reference sources and one on a calibration approach using the diffuse radio emission of the Galaxy, optimized for short data-sets. An accuracy of 35% in amplitude is reached. The absolute calibration is also compared to predictions from air shower simulations. These results are used to set an absolute energy scale for air shower measurements and can be used as a basis for an absolute scale for the measurement of astronomical transients with LOFAR.
Using data taken as part of the Bluedisk project we study the connection between neutral hydrogen (HI) in the environment of spiral galaxies and that in the galaxies themselves. We measure the total HI mass present in the environment in a statistical way by studying the distribution of noise peaks in the HI data cubes obtained for 40 galaxies observed with WSRT. We find that galaxies whose HI mass fraction is high relative to standard scaling relations have an excess HI mass in the surrounding environment as well. Gas in the environment consists of gas clumps which are individually below the detection limit of our HI data. These clumps may be hosted by small satellite galaxies and\or be the high-density peaks of a more diffuse gas distribution in the inter-galactic medium. We interpret this result as an indication for a picture in which the HI-rich central galaxies accrete gas from an extended gas reservoir present in their environment.
We reconstruct the projected mass distribution of a massive merging Hubble Frontier Fields cluster MACSJ0416 using the genetic algorithm based free-form technique called Grale. The reconstructions are constrained by 149 lensed images identified by Jauzac et al. using HFF data. No information about cluster galaxies or light is used, which makes our reconstruction unique in this regard. Using visual inspection of the maps, as well as galaxy-mass correlation functions we conclude that overall light does follow mass. Furthermore, the fact that brighter galaxies are more strongly clustered with mass is an important confirmation of the standard biasing scenario in galaxy clusters. On the smallest scales, approximately less than a few arcseconds the resolution afforded by 149 images is still not sufficient to confirm or rule out galaxy-mass offsets of the kind observed in ACO 3827. We also compare the mass maps of MACSJ0416 obtained by three different groups: Grale, and two parametric Lenstool reconstructions from the CATS and Sharon/Johnson teams. Overall, the three agree well; one interesting discrepancy between Grale and Lenstool galaxy-mass correlation functions occurs on scales of tens of kpc and may suggest that cluster galaxies are more biased tracers of mass than parametric methods generally assume.
The 2-point angular correlation function $w(\theta)$ (2PACF), where $\theta$ is the angular separation between pairs of galaxies, provides the transversal Baryon Acoustic Oscillation (BAO) signal almost model-independently. In this paper we use 409,337 luminous red galaxies in the redshift range $z = [0.440,0.555]$ obtained from the tenth data release of the Sloan Digital Sky Survey (SDSS DR10) to estimate $\theta_{\rm{BAO}}(z)$ from the 2PACF at six redshift {shells}. Since noise and systematics can hide the BAO signature in the $w - \theta$ plane, we also discuss some criteria to localize the acoustic bump. We identify two sources of model-dependence in the analysis, namely, the value of the acoustic scale from Cosmic Microwave Background (CMB) measurements and the correction in the $\theta_{\rm{BAO}}(z)$ position due to projection effects. Constraints on the dark energy equation-of-state parameter w$(z)$ from the $\theta_{\rm{BAO}}(z)$ diagram are derived, as well as from a joint analysis with current CMB measurements. We find that the standard $\Lambda$CDM model as well as some of its extensions are in good agreement with these $\theta_{\rm{BAO}}(z)$ measurements.
We report on Nuclear Spectroscopic Telescope Array (NuSTAR) hard X-ray observations of the young rotation-powered radio pulsar PSR B1509$-$58 in the supernova remnant MSH 15$-$52. We confirm the previously reported curvature in the hard X-ray spectrum, showing that a log parabolic model provides a statistically superior fit to the spectrum compared with the standard power law. The log parabolic model describes the NuSTAR data, as well as previously published gamma-ray data obtained with COMPTEL and AGILE, all together spanning 3 keV through 500 MeV. Our spectral modelling allows us to constrain the peak of the broadband high energy spectrum to be at 2.6$\pm$0.8 MeV, an improvement of nearly an order of magnitude over previous measurements. In addition, we calculate NuSTAR spectra in 26 pulse phase bins and confirm previously reported variations of photon indices with phase. Finally, we measure the pulsed fraction of PSR B1509$-$58 in the hard X-ray energy band for the first time. Using the energy resolved pulsed fraction results, we estimate that the pulsar's DC component has a photon index value between 1.26 and 1.96. Our results support a model in which the pulsar's lack of GeV emission is due to viewing geometry, with the X-rays originating from synchrotron emission from secondary pairs in the magnetosphere.
Ultra-high energy neutrinos are interesting messenger particles since, if detected, they can transmit exclusive information about ultra-high energy processes in the Universe. These particles, with energies above $10^{16}\mathrm{eV}$, interact very rarely. Therefore, detectors that instrument several gigatons of matter are needed to discover them. The ARA detector is currently being constructed at South Pole. It is designed to use the Askaryan effect, the emission of radio waves from neutrino-induced cascades in the South Pole ice, to detect neutrino interactions at very high energies. With antennas distributed among 37 widely-separated stations in the ice, such interactions can be observed in a volume of several hundred cubic kilometers. Currently 3 deep ARA stations are deployed in the ice of which two have been taking data since the beginning of the year 2013. In this publication, the ARA detector "as-built" and calibrations are described. Furthermore, the data reduction methods used to distinguish the rare radio signals from overwhelming backgrounds of thermal and anthropogenic origin are presented. Using data from only two stations over a short exposure time of 10 months, a neutrino flux limit of $3 \cdot 10^{-6} \mathrm{GeV} / (\mathrm{cm^2 \ s \ sr})$ is calculated for a particle energy of 10^{18}eV, which offers promise for the full ARA detector.
We present results from gamma-ray observations of the Coma cluster incorporating 6 years of Fermi-LAT data and the newly released "Pass 8" event-level analysis. Our analysis of the region reveals low-significance residual structures within the virial radius of the cluster that are too faint for a detailed investigation with the current data. Using a likelihood approach that is free of assumptions on the spectral shape we derive upper limits on the gamma-ray flux that is expected from energetic particle interactions in the cluster. We also consider a benchmark spatial and spectral template motivated by models in which the observed radio halo is mostly emission by secondary electrons. In this case, the median expected and observed upper limits for the flux above 100 MeV are 1.7x10-9 ph/cm2/s and 4.2x10-9 ph/cm2/s respectively. These bounds are comparable to several recent predictions for the hadronic gamma-ray emission in such secondary models, although direct comparisons are sensitive to assumptions regarding the magnetic field strength and other factors. The minimal expected gamma-ray flux from radio and star-forming galaxies within the Coma cluster is roughly an order of magnitude below the median sensitivity of our analysis.
We describe a general scenario, dubbed "Inflatable Dark Matter", in which the density of dark matter particles can be reduced through a short period of late-time inflation in the early universe. The overproduction of dark matter that is predicted within many otherwise well-motivated models of new physics can be elegantly remedied within this context, without the need to tune underlying parameters or to appeal to anthropic considerations. Thermal relics that would otherwise be disfavored can easily be accommodated within this class of scenarios, including dark matter candidates that are very heavy or very light. Furthermore, the non-thermal abundance of GUT or Planck scale axions can be brought to acceptable levels, without invoking anthropic tuning of initial conditions. A period of late-time inflation could have occurred over a wide range of scales from ~ MeV to the weak scale or above, and could have been triggered by physics within a hidden sector, with small but not necessarily negligible couplings to the Standard Model.
A possible solution for the problem of memory-size and computer-time, is the extrapolation of basis-set$^1$. This extrapolation has two exponents $\alpha$ and $\beta$, corresponding to the HF (reference energy) and the energy of correlations (EC), respectively. For a given system, the exponents are taken as constant$^2$, and potential energy surfaces (PES) are generated. We have found that the values of $\alpha$ and $\beta$ are not constant, but vary from position to position in the system. How to deal with such situation and get very accurate PES, is discussed.
Society needs to prepare for more severe space weather than it has experienced in the modern technological era. To enable that, we must both quantify extreme-event characteristics and analyze impacts of lesser events that are frequent yet severe enough to be informative. Exploratory studies suggest that economic impacts of a century-level space hurricane and of a century of lesser space-weather "gales" may turn out to be of the same order of magnitude. The economic benefits of effective mitigation of the impacts of space gales may substantially exceed the required investments, even as these investments provide valuable information to prepare for the worst possible storms.
We study spectral features in the gamma-ray emission from dark matter (DM) annihilation in the Next-to-Minimal Supersymmetric Standard Model (NMSSM), with either neutralino or right-handed (RH) sneutrino DM. We perform a series of scans over the NMSSM parameter space, compute the DM annihilation cross section into two photons and the contribution of box-shaped features, and compare them with the limits derived from the Fermi-LAT search for gamma-ray lines using the latest Pass 8 data. We implement the LHC bounds on the Higgs sector and on the masses of supersymmetric particles as well as the constraints on low-energy observables. We also consider the recent upper limits from the Fermi-LAT satellite on the continuum gamma-ray emission from dwarf spheroidal galaxies (dSphs). We show that in the case of the RH sneutrino the constraint on gamma-ray spectral features can be more stringent than the dSphs bounds. This is due to the Breit-Wigner enhancement near the ubiquitous resonances with a CP even Higgs and the contribution of scalar and pseudoscalar Higgs final states to box-shaped features. By contrast, for neutralino DM, the di-photon final state is only enhanced in the resonance with a $Z$ boson and box-shaped features are even more suppressed. Therefore, the observation of spectral features could constitute a discriminating factor between both models. In addition, we compare our results with direct DM searches, including the SuperCDMS and LUX limits on the elastic DM-nucleus scattering cross section and show that some of these scenarios would be accessible to next generation experiments. Thus, our findings strengthen the idea of complementarity among distinct DM search strategies.
Neutrino flavor oscillations in the presence of ambient neutrinos is nonlinear in nature which leads to interesting phenomenology that has not been well understood. It was recently shown that, in the two-dimensional, two-beam neutrino Line model, the inhomogeneous neutrino oscillation modes on small distance scales can become unstable at larger neutrino densities than the homogeneous mode does. We develop a numerical code to solve neutrino oscillations in the multi-angle/beam Line model with a continuous neutrino angular distribution. We show that the inhomogeneous oscillation modes can occur at even higher neutrino densities in the multi-angle model than in the two-beam model. We also find that the inhomogeneous modes on sufficiently small scales can be unstable at smaller neutrino densities with ambient matter than without, although a larger matter density does shift the instability region of the homogeneous mode to higher neutrino densities in the Line model as it does in the one-dimensional supernova Bulb model. Our results suggest that the inhomogeneous neutrino oscillation modes can be difficult to treat numerically because the problem of spurious oscillations becomes more severe for oscillations on smaller scales.
Links to: arXiv, form interface, find, astro-ph, recent, 1508, contact, help (Access key information)
We present spectral and timing analyses of the X-ray emission from the pulsar wind nebula DA 495 and its central object, J1952.2+2925, suggested to be the pulsar, using archival Chandra and XMM-Newton data. J1952.2+2925 has a pure thermal spectrum which is equally well fitted either by the blackbody model with a temperature of $\approx 215$ eV and an emitting area radius of $\approx 0.6$ km or by magnetized neutron star atmosphere models with temperatures of 80-90 eV. In the latter case the thermal emission can come from the entire neutron star surface which temperature is consistent with standard neutron star cooling scenarios. We place also an upper limit on the J1952.2+2925 nonthermal flux. The derived spectral parameters are generally compatible with published ones based only on the Chandra data, but they are much more accurate due to the inclusion of XMM-Newton data. No pulsations were found and we placed an upper limit for the J1952.2+2925 pulsed emission fraction of 40 per cent. Utilizing the interstellar absorption-distance relation, we estimated the distance to DA 495, which can be as large as 5 kpc if J1952.2+2925 emission is described by the atmosphere models. We compiled possible multi-wavelength spectra of the nebula including radio data; they depend on the spectral model of the central object. Comparing the results with other pulsar plus wind nebula systems we set reasonable constraints on the J1952.2+2925 spin-down luminosity and age. We suggest that the Fermi source 3FGL J1951.6+2926 is the likely $\gamma$-ray counterpart of J1952.2+2925.
We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar coordinates, and fully 3D cylindrical coordinates. Viscosity and diffusivity are implemented in the code to control the accretion rate in the disk and the rate of penetration of the disk matter through the magnetic field lines. The code has been utilized for the numerical investigations of a number of different astrophysical problems, several examples of which are shown.
We follow the premise that most intermediate luminosity optical transients (ILOTs) are powered by rapid mass accretion onto a main sequence star, and study the effects of jets launched by an accretion disk. The disk is formed due to large specific angular momentum of the accreted mass. The two opposite jets might expel some of the mass from the reservoir of gas that feeds the disk, hence reducing and shortening the mass accretion process. We argue that by this process ILOTs limit their luminosity and might even shut themselves off in this negative jet feedback mechanism (JFM). The group of ILOTs is a new member of a large family of astrophysical objects whose activity is regulated by the operation of the JFM.
We examine the scalings of galactic outflows with halo mass across a suite of twenty high-resolution cosmological zoom galaxy simulations covering halo masses from 10^9.5 - 10^12 M_sun. These simulations self-consistently generate outflows from the available supernova energy in a manner that successfully reproduces key galaxy observables including the stellar mass-halo mass, Tully-Fisher, and mass-metallicity relations. We quantify the importance of ejective feedback to setting the stellar mass relative to the efficiency of gas accretion and star formation. Ejective feedback is increasingly important as galaxy mass decreases; we find an effective mass loading factor that scales as v_circ^(-2.2), with an amplitude and shape that is invariant with redshift. These scalings are consistent with analytic models for energy-driven wind, based solely on the halo potential. Recycling is common: about half the outflow mass across all galaxy masses is later re-accreted. The recycling timescale is typically about 1 Gyr, virtually independent of halo mass. Recycled material is re-accreted farther out in the disk and with typically about 2-3 times more angular momentum. These results elucidate and quantify how the baryon cycle plausibly regulates star formation and alters the angular momentum distribution of disk material across the halo mass range where most of cosmic star formation occurs.
X-ray images of galaxy clusters and gas-rich elliptical galaxies show a wealth of small-scale features which reflect fluctuations in density and/or temperature of the intra-cluster medium. In this paper we study these fluctuations in M87/Virgo, to establish whether sound waves/shocks, bubbles or uplifted cold gas dominate the structure. We exploit the strong dependence of the emissivity on density and temperature in different energy bands to distinguish between these processes. Using simulations we demonstrate that our analysis recovers the leading type of fluctuation even in the presence of projection effects and temperature gradients. We confirm the isobaric nature of cool filaments of gas entrained by buoyantly rising bubbles, extending to 7' to the east and south-west, and the adiabatic nature of the weak shocks at 40" and 3' from the center. For features of 5--10 kpc, we show that the central 4'x 4' region is dominated by cool structures in pressure equilibrium with the ambient hotter gas while up to 30 percent of the variance in this region can be ascribed to adiabatic fluctuations. The remaining part of the central 14'x14' region, excluding the arms and shocks described above, is dominated by apparently isothermal fluctuations (bubbles) with a possible admixture (at the level of about 30 percent) of adiabatic (sound waves) and by isobaric structures. Larger features, of about 30 kpc, show a stronger contribution from isobaric fluctuations. The results broadly agree with an AGN feedback model mediated by bubbles of relativistic plasma.
We study the effect of surface brightness on the mass-metallicity relation using nearby galaxies whose gas content and metallicity profiles are available. Previous studies using fiber spectra indicated that lower surface brightness galaxies have systematically lower metallicity for their stellar mass, but the results were uncertain because of aperture effect. With stellar masses and surface brightnesses measured at WISE W1 and W2 bands, we re-investigate the surface brightness dependence with spatially-resolved metallicity profiles and find the similar result. We further demonstrate that the systematical difference cannot be explained by the gas content of galaxies. For two galaxies with similar stellar and gas masses, the one with lower surface brightness tends to have lower metallicity. Using chemical evolution models, we investigate the inflow and outflow properties of galaxies of different masses and surface brightnesses. We find that, on average, high mass galaxies have lower inflow and outflow rates relative to star formation rate. On the other hand, lower surface brightness galaxies experience stronger inflow than higher surface brightness galaxies of similar mass. The surface brightness effect is more significant for low mass galaxies. We discuss implications on the different inflow properties between low and high surface brightness galaxies, including star formation efficiency, environment and mass assembly history.
We consider the effectiveness of foreground cleaning in the recovery of Cosmic Microwave Background (CMB) polarization sourced by gravitational waves for tensor-to-scalar ratios in the range $0<r<0.1$. Using the planned survey area, frequency bands, and sensitivity of the Cosmology Large Angular Scale Surveyor (CLASS), we simulate maps of Stokes $Q$ and $U$ parameters at 40, 90, 150, and 220 GHz, including realistic models of the CMB, diffuse Galactic thermal dust and synchrotron foregrounds, and Gaussian white noise. We use linear combinations (LCs) of the simulated multifrequency data to obtain maximum likelihood estimates of $r$, the relative scalar amplitude $s$, and LC coefficients. We find that for 10,000 simulations of a CLASS-like experiment using only measurements of the reionization peak ($\ell\leq23$), there is a 95% C.L. upper limit of $r<0.017$ in the case of no primordial gravitational waves. For simulations with $r=0.01$, we recover at 68% C.L. $r=0.012^{+0.011}_{-0.006}$. The reionization peak corresponds to a fraction of the multipole moments probed by CLASS, and simulations including $30\leq\ell\leq100$ further improve our upper limits to $r<0.008$ at 95% C.L. ($r=0.01^{+0.004}_{-0.004}$ for primordial gravitational waves with $r=0.01$). In addition to decreasing the current upper bound on $r$ by an order of magnitude, these foreground-cleaned low multipole data will achieve a cosmic variance limited measurement of the E-mode polarization's reionization peak.
The main challenge for understanding the fuelling of supermassive black holes in active galactic nuclei is not to account for the source of fuel, but rather to explain its delivery from the boundaries of the black hole sphere of influence (10-100 pc) down to sub-parsec scales. In this work, we report on a series of numerical experiments aimed at exploring in further depth our model of "overlapping inflow events" as catalysts for rapid accretion, seeding a turbulent field in the infalling gas. We initially set a gaseous shell in non-equilibrium rotation around a supermassive black hole. After infall, the shell stalls in a disc-like structure. A second shell is then set in either co-rotation or counter-rotation with respect to the first and is let to impinge on the previously-formed disc. We find that combined turbulence and overlap significantly enhance accretion in counter-rotating inflows, while turbulence dominates for co-rotating inflows. The leftovers of overlapping inflows are warped nuclear discs, whose morphology depends on the relative orientation and angular momentum of the disc and the shell. Overlapping inflows leave observational signatures in the gas rotation curves.
We report the discovery of a relic Giant Radio Galaxy (GRG) J021659-044920 at redshift $z \sim 1.3$ that exhibits large-scale extended, nearly co-spatial, radio and X-ray emission from radio lobes, but no detection of AGN core, jets and hot-spots. The total angular extent of the GRG at the observed frame 0.325 GHz, using GMRT observations is found to be ${\sim}$ 2$^{\prime}$.4, that corresponds to a total projected linear size of $\sim$ 1.2 Mpc. The integrated radio spectrum between 0.240 GHz to 1.4 GHz shows high spectral curvature (${\alpha}_{\rm 0.610~GHz}^{\rm 1.4~GHz}$ - ${\alpha}_{\rm 0.240~GHz}^{\rm 0.325~GHz}$ $>$ 1.19) with sharp steepening above 0.325 GHz, consistent with relic radio emission that is $\sim$ 8 $\times$ 10$^{6}$ years old. The radio spectral index map between observed frame 0.325 and 1.4~GHz for the two lobes varies from 1.4 to 2.5 with the steepening trend from outer-end to inner-end, indicating backflow of plasma in the lobes. The extended X-ray emission characterized by an absorbed power law with photon index $\sim$ 1.86 favours inverse-Compton scattering of the Cosmic Microwave Background (ICCMB) photons as the plausible origin. Using both X-ray and radio fluxes under the assumption of ICCMB we estimate the magnetic field in the lobes to be 3.3 $\mu$G. The magnetic field estimate based on energy equipartition is $\sim$ 3.5 $\mu$G. Our work presents a case study of a rare example of a GRG caught in dying phase in the distant universe.
In the first paper of this series, we studied the effect of baryon acoustic oscillations (BAO), redshift space distortions (RSD) and weak lensing (WL) on measurements of angular cross-correlations in narrow redshift bins. Paper-II presented a multitracer forecast as Figures of Merit (FoM), combining a photometric and spectroscopic stage-IV survey. The uncertainties from galaxy bias, the way light traces mass, is an important ingredient in the forecast. Fixing the bias would increase our FoM equivalent to 3.3 times larger area for the combined constraints. This paper focus on how the modelling of bias affect these results. In the combined forecast, lensing both help and benefit from the improved bias measurements in overlapping surveys after marginalizing over the cosmological parameters. Adding a second lens population in counts-shear does not have a large impact on bias error, but removing all counts-shear information increases the bias error in a significant way. We also discuss the relative impact of WL, magnification, RSD and BAO, and how results change as a function of bias amplitude, photo-z error and sample density. By default we use one bias parameter per bin (with 72 narrow bins), but we show that the results do not change much when we use other parameterizations, with at least 3 parameters in total. Bias stochasticity, even when added as one new free parameter per bin, only produce moderate decrease in the FoM. In general, we find that the degradation in the figure of merit caused by the uncertainties in the knowledge of bias is significantly smaller for overlapping surveys.
We aimed to verify the nature and derive the basic parameters of the symbiotic star candidate [JD2002] 11. For this purpose, we obtained and analysed an X-Shooter spectrum of [JD2002] 11. We also used optical and infrared photometry available for the object. Emission-line diagnostic ratios are characteristic of a dusty type symbiotic star and reveal a two-component nebula (low- and high-density). The spectral energy distribution is well fitted with a two-component blackbody spectrum with the respective temperatures of 1150 K and 600 K. The total luminosity of $\rm 3000\,L_{\odot}$ is consistent with the expected luminosity of a typical Mira star, embedded in an optically thick dust shell. We conclude that [JD2002] 11 is the ninth symbiotic star in total and only the second dusty type symbiotic star discovered in the Small Magellanic Cloud.
Simulations of galaxy formation follow the gravitational and hydrodynamical interactions between gas, stars and dark matter through cosmic time. The huge dynamic range of such calculations severely limits strong scaling behaviour of the community codes in use, with load-imbalance, cache inefficiencies and poor vectorisation limiting performance. The new swift code exploits task-based parallelism designed for many-core compute nodes interacting via MPI using asynchronous communication to improve speed and scaling. A graph-based domain decomposition schedules interdependent tasks over available resources. Strong scaling tests on realistic particle distributions yield excellent parallel efficiency, and efficient cache usage provides a large speed-up compared to current codes even on a single core. SWIFT is designed to be easy to use by shielding the astronomer from computational details such as the construction of the tasks or MPI communication. The techniques and algorithms used in SWIFT may benefit other computational physics areas as well, for example that of compressible hydrodynamics. For details of this open-source project, see www.swiftsim.com
In a microlensing event, a large magnification occurs at caustic crossing and
provides an opportunity to obtain a stronger signal associated with the object.
In this paper we study the possibility of magnetic field detection in a
microlensing event through the Zeeman effect. We follow the prescription
introduced by \citep{Robinson:1980} which analyses the spectrum of a star in a
Fourier space to deconvolve other broadening mechanism from the actual Zeeman
effect. First we study magnification contrast between source and spot in terms
of spot size and then consider two distinct strategies using modern
spectrographs and perform a Monte Carlo simulation to find the detection
efficiency in magnification-magnetic field plan. The spectral resolution and
signal to noise ratio in each strategies, specify suitable places in this plane
to detect a magnetic field.
Apart from complexity of magnetic field on a star, we consider a simple model
and propose a fantastic method to detect magnetic field using spectrum of
source stars at caustic crossing.
The last 20 years have seen the development of new techniques in
Astroparticle Physics providing access to the highest end of the
electromagnetic spectrum. It has been shown that some sources emit photons up
to energies close to 100 TeV. Yet the fluxes of these photons are incredibly
low and new detection techniques are needed to go higher in energy.
A new technique that would use the new generation of Cherenkov Telescopes,
i.e., the Cherenkov Telescope Array (CTA), is proposed to push further the
energy frontier. It is based on the detection of the fluorescence radiation
emitted in extensive air showers, a successful method used in ultra-high-energy
cosmic ray experiments, like the Pierre Auger Observatory. It would complement
the standard imaging atmospheric Cherenkov technique with only minor
modifications of the hardware currently being developed for the CTA and would
not imply significant extra costs during its planned operation.
The measured fluxes of secondary particles produced by the interactions of cosmic rays with the astronomical environment represent a powerful tool to infer some properties of primary cosmic rays. In this work we investigate the production of secondary particles in inelastic hadronic interactions between several cosmic rays species of projectiles and different target nuclei of the interstellar medium. The yields of secondary particles have been calculated with the FLUKA simulation package, that provides with very good accuracy the energy distributions of secondary products in a large energy range. An application to the propagation and production of secondaries in the Galaxy is presented.
We present meter wave solar radio spectra of the highest spectrotemporal resolution achieved to date. The observations, obtained with the first station of the Long Wavelength Array (LWA1), show unprecedented detail of solar emissions across a wide bandwidth during a Type III/IIIb storm. Our flux calibration demonstrates that the LWA1 can detect Type III bursts much weaker than 1 SFU, much lower than previous observations, and that the distribution of fluxes in these bursts varies with frequency. The high sensitivity and low noise in the data provide strong constraints to models of this type of plasma emission. The continuous generation of electron beams in the corona revealed by the high density Type III storm is evidence for ubiquitous magnetic reconnection in the lower corona. Such an abundance of reconnection events not only contributes to the total coronal energy budget, but also provides an engine by which to form the populations of seed particles responsible for proton-rich solar energetic particle events. An active region with such levels of reconnection and the accompanying type III/IIIb storms is here proposed to be associated with an increase of SEP production if a CME erupts. The data's constraints on existing theories of type IIIIb production are used to make an association of the observed type IIIb storm to specific electron beam paths with increased inhomogeneities in density, temperature, and or turbulence. This scenario ties in the observed timing of III and IIIb storms, constrained theories of type III and IIIb emission, and the ability of the emitting AR to produce a strong SEP event. The result requires but a single observable to cement these ideas, the statistical correlation of type III/IIIb activity with SEP-productive AR.
We show that the standing accretion shock instability (SASI) that has been used to ease the shock revival in core collapse supernovae (CCSNe) neutrino-driven explosion models, might play a much more decisive role in supplying the stochastic angular momentum required to trigger an explosion with jittering jets. To play a minor role in neutrino-based explosion models, the kinetic energy of the gas inside the stalled shock associated with the transverse (non-radial) motion should be about more than ten percent of the energy of the accreted gas. We find that this implies a stochastic angular momentum that can reach about five percent of the Keplerian specific angular momentum around the newly born neutron star. Such an accretion flow leaves an open conical region along the poles with an average opening angle of about 5 degrees. The outflow from the open polar region powers an explosion according to the jittering-jets model.
Previous observations of the Ap star HD 32633 indicated that its magnetic field was unusually complex in nature and could not be characterised by a simple dipolar structure. Here we derive magnetic field maps and chemical abundance distributions for this star using full Stokes vector (Stokes $IQUV$) high-resolution observations obtained with the ESPaDOnS and Narval spectropolarimeters. Our maps, produced using the Invers10 magnetic Doppler imaging (MDI) code, show that HD 32633 has a strong magnetic field which features two large regions of opposite polarity but deviates significantly from a pure dipole field. We use a spherical harmonic expansion to characterise the magnetic field and find that the harmonic energy is predominately in the $\ell=1$ and $\ell=2$ poloidal modes with a small toroidal component. At the same time, we demonstrate that the observed Stokes parameter profiles of HD 32633 cannot be fully described by either a dipolar or dipolar plus quadrupolar field geometry. We compare the magnetic field topology of HD 32633 with other early-type stars for which MDI analyses have been performed, supporting a trend of increasing field complexity with stellar mass. We then compare the magnetic field topology of HD 32633 with derived chemical abundance maps for the elements Mg, Si, Ti, Cr, Fe, Ni and Nd. We find that the iron-peak elements show similar distributions, but we are unable to find a clear correlation between the location of local chemical enhancements or depletions and the magnetic field structure.
On the basis of the 30-year ago ground-based photometry and the recent Kepler space experiment there have been considered frequencies of occurrence and energetics of the solar-type stellar flares. It was concluded that frequencies of occurrence of such flares are proportional to sizes of stellar surfaces, and estimates of maximum flare radiation from the results of the ground-based photometry and space observations practically coincide.
A numerical study of a pseudoscalar inflation having an axion-photon-like coupling is performed by solving numerically the coupled differential equations of motion for inflaton and photon mode functions from the onset of inflation to the end of reheating. The backreaction due to particle production is also included self-consistently. We find that this particular inflation model realizes the idea of a warm inflation in which a steady thermal bath is established by the particle production. In most cases this thermal bath exceeds the amount of radiation released in the reheating process. In the strong coupling regime, the transition from the inflationary to the radiation-dominated phase does not involve either a preheating or reheating process. In addition, energy density peaks produced near the end of inflation may lead to the formation of primordial black holes.
We analyze new observations of superflares on G-stars discovered in the optical and near IR ranges with the Kepler mission. An evolution of solar-type activity is discussed. We give an estimate of the maximal total energy, $E_{tot} = 10^{34}\;\mbox{erg}$ of a flare that can occur on the young Sun at its age of 1 Gyr when the cycle was formed. We believe that the main source of the flare optical continuum is a low-temperature condensation forming in the course of the response of the chromosphere to an impulsive heating. For a superflare on the young Sun, we adopt the accelerated electron flux, $F_e (E>\mbox{20 keV}) = 3 \times 10^{11} \: \mbox{erg} \; \mbox{cm}^{-2} \; \mbox{s}^{-1}$, that is limited by the return current, and obtain the area of the optical continuum source on a G star, $S \approx 10^{19} \:\mbox{cm}^2$. This value is close to the area of the $H_\alpha$-ribbons in the largest solar flares, while the area of bright patches of a white-light flare on the contemporary Sun is smaller by about two orders of magnitude. At the same electron flux and the hard electron spectrum, the stellar flare of the similar energy should be accompanied by the microwave source of about 2 mJy at frequencies 10--100 GHz at distance 100 pc. We discuss the possible detection of the flare-produced lithium in the course of spallation reactions. The detection of the flare microwave source and the emission in the Li resonant line could demonstrate how effective can be particle acceleration on stars in the lower part of the main sequence.
Activity of the chromospheres and coronas of the stars of late (F, G, K) spectral classes and also their cyclic activity is analyzed on the observational data of solar-type stars from Mount Wilson "HK-project". The chromospheric and coronal activity of "HK-project stars" and the Sun was compared to the stars with active atmospheres observed in Planets Search Programs. It is shown that the cyclic activity begins to occur in stars, close by their spectral type to the Sun and becomes more pronounced at K-stars. The results of a comparative analysis of the characteristics of chromospheric and coronal activity of the Sun and F, G and K-stars are presented.
The life cycle of the population of interstellar polycyclic aromatic hydrocarbon (PAH) molecules depends partly on the photostability of the individual species. We have studied the dissociative photoionization of two ethynyl-substituted PAH species, namely, 9-ethynylphenanthrene and 1-ethynylpyrene. Their adiabatic ionization energy and the appearance energy of fragment ions have been measured with the photoelectron photoion coincidence (PEPICO) spectroscopy technique. The adiabatic ionization energy has been found at 7.84 +/- 0.02 eV for 9-ethynylphenanthrene and at 7.41 +/- 0.02 eV for 1-ethynylpyrene. These values are similar to those determined for the corresponding non-substituted PAH molecules phenanthrene and pyrene. The appearance energy of the fragment ion indicative of the loss of a H atom following photoionization is also similar for either ethynyl-substituted PAH molecule and its non-substituted counterpart. The measurements are used to estimate the critical energy for the loss of a H atom by the PAH cations and the stability of ethynyl-substituted PAH molecules upon photoionization. We conclude that these PAH derivatives are as photostable as the non-substituted species in HI regions. If present in the interstellar medium, they may play an important role in the growth of interstellar PAH molecules.
The appearance of debris disks around distant stars depends upon the scattering/phase function (SPF) of the material in the disk. However, characterizing the SPFs of these extrasolar debris disks is challenging because only a limited range of scattering angles are visible to Earth-based observers. By contrast, Saturn's tenuous rings can be observed over a much broader range of geometries, so their SPFs can be much better constrained. Since these rings are composed of small particles released from the surfaces of larger bodies, they are reasonable analogs to debris disks and so their SPFs can provide insights into the plausible scattering properties of debris disks. This work examines two of Saturn's dusty rings: the G ring (at 167,500 km from Saturn's center) and the D68 ringlet (at 67,600 km). Using data from the cameras onboard the Cassini spacecraft, we are able to estimate the rings' brightnesses at scattering angles ranging from 170 to 0.5 degrees. We find that both of the rings exhibit extremely strong forward-scattering peaks, but for scattering angles above 60 degrees their brightnesses are nearly constant. These SPFs can be well approximated by a linear combination of three Heyney-Greenstein functions, and are roughly consistent with the SPFs of irregular particles from laboratory measurements. Comparing these data to Fraunhofer and Mie models highlights several challenges involved in extracting information about particle compositions and size distributions from SPFs alone. The SPFs of these rings also indicate that the degree of forward scattering in debris disks may be greatly underestimated.
Comets and chondrites show non-monotonic behaviour of their Deuterium to
Hydrogen (D/H) ratio as a function of their formation location from the Sun.
This is difficult to explain with a classical protoplanetary disk model that
has a decreasing temperature structure with radius from the Sun.
We want to understand if a protoplanetary disc with a dead zone, a region of
zero or low turbulence, can explain the measured D/H values in comets and
chondrites.
We use time snapshots of a vertically layered disk model with turbulent
surface layers and a dead zone at the midplane. The disc has a non-monotonic
temperature structure due to increased heating from self-gravity in the outer
parts of the dead zone. We couple this to a D/H ratio evolution model in order
to quantify the effect of such thermal profiles on D/H enrichment in the
nebula.
We find that the local temperature peak in the disk can explain the diversity
in the D/H ratios of different chondritic families. This disk temperature
profile leads to a non-monotonic D/H enrichment evolution, allowing these
families to acquire their different D/H values while forming in close
proximity. The formation order we infer for these families is compatible with
that inferred from their water abundances. However, we find that even for very
young disks, the thermal profile reversal is too close to the Sun to be
relevant for comets.
In this paper my prime objective is to analyze the constraints on a sub-Planckian excursion of a single inflaton field within Effective Field Theory framework in a model independent fashion. For a generic single field inflationary potential, using the various parameterization of the primordial power spectrum I have derived the most general expression for the field excursion in terms of various inflationary observables, applying the observational constraints obtained from recent Planck 2015 and Planck 2015 +BICEP2/Keck Array data. By explicit computation I have reconstructed the structural form of the inflationary potential by constraining the Taylor expansion coefficients appearing in the generic expansion of the potential within the Effective Field Theory. Next I have explicitly derived, a set of higher order inflationary consistency relationships, which would help us to break the degeneracy between various class of inflationary models by differentiating them. I also provided two simple examples of Effective Theory of inflation- inflection-point model and saddle-point model to check the compatibility of the prescribed methodology in the light of Planck 2015 and Planck 2015 +BICEP2/Keck Array data. Finally, I have also checked the validity of the prescription by estimating the cosmological parameters and fitting the theoretical CMB TT, TE and EE angular power spectra with the observed data within the multipole range $2<l<2500$.
We present pre-perihelion infrared 8 to 31 micron spectrophotometric and imaging observations of comet C/2012 K1 (Pan-STARRS), a dynamically new Oort Cloud comet, conducted with NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) facility (+FORCAST) in 2014 June. As a "new" comet (first inner solar system passage), the coma grain population may be extremely pristine, unencumbered by a rime and insufficiently irradiated by the Sun to carbonize its surface organics. The comet exhibited a weak 10 micron silicate feature ~1.18 +/- 0.03 above the underlying best-fit 215.32 +/- 0.95 K continuum blackbody. Thermal modeling of the observed spectral energy distribution indicates that the coma grains are fractally solid with a porosity factor D = 3 and the peak in the grain size distribution, a_peak = 0.6 micron, large. The sub-micron coma grains are dominated by amorphous carbon, with a silicate-to-carbon ratio of 0.80 (+0.25) (- 0.20). The silicate crystalline mass fraction is 0.20 (+0.30) (-0.10), similar to with other dynamically new comets exhibiting weak 10 micron silicate features. The bolometric dust albedo of the coma dust is 0.14 +/- 0.01 at a phase angle of 34.76 degrees, and the average dust production rate, corrected to zero phase, at the epoch of our observations was Afrho ~ 5340~cm.
Many attempts have already been made for detecting exomoons around transiting exoplanets but the first confirmed discovery is still pending. The experience that have been gathered so far allow us to better optimize future space telescopes for this challenge, already during the development phase. In this paper we focus on the forthcoming CHaraterising ExOPlanet Satellite (CHEOPS),describing an optimized decision algorithm with step-by-step evaluation, and calculating the number of required transits for an exomoon detection for various planet-moon configurations that can be observable by CHEOPS. We explore the most efficient way for such an observation which minimizes the cost in observing time. Our study is based on PTV observations (photocentric transit timing variation, Szab\'o et al. 2006) in simulated CHEOPS data, but the recipe does not depend on the actual detection method, and it can be substituted with e.g. the photodynamical method for later applications. Using the current state-of-the-art level simulation of CHEOPS data we analyzed transit observation sets for different star-planet-moon configurations and performed a bootstrap analysis to determine their detection statistics. We have found that the detection limit is around an Earth-sized moon. In the case of favorable spatial configurations, systems with at least such a large moon and with at least Neptune-sized planet, 80\% detection chance requires at least 5-6 transit observations on average. There is also non-zero chance in the case of smaller moons, but the detection statistics deteriorates rapidly, while the necessary transit measurements increase fast. (abridged)
This article investigates the problem of estimating stellar atmospheric parameters from spectra. Feature extraction is a key procedure in estimating stellar parameters automatically. We propose a scheme for spectral feature extraction and atmospheric parameter estimation using the following three procedures: firstly, learn a set of basic structure elements (BSE) from stellar spectra using an autoencoder; secondly, extract representative features from stellar spectra based on the learned BSEs through some procedures of convolution and pooling; thirdly, estimate stellar parameters ($T_{eff}$, log$~g$, [Fe/H]) using a back-propagation (BP) network. The proposed scheme has been evaluated on both real spectra from Sloan Digital Sky Survey (SDSS)/Sloan Extension for Galactic Understanding and Exploration (SEGUE) and synthetic spectra calculated from Kurucz's new opacity distribution function (NEWODF) models. The best mean absolute errors (MAEs) are 0.0060 dex for log$~T_{eff}$, 0.1978 dex for log$~g$ and 0.1770 dex for [Fe/H] for the real spectra and 0.0004 dex for log$~T_{eff}$, 0.0145 dex for log$~g$ and 0.0070 dex for [Fe/H] for the synthetic spectra.
A new set of laws of motion for turbulent jets propagating in an intergalactic medium characterized by a decreasing density is derived by applying conservation of momentum flux both in the classical and relativistic framework. Two characteristic features of radio-galaxies, such as oscillations and curvature, are modeled by a classical helicoidal jet. A third feature of a radio-galaxy, the appearance of knots, is explained as an effect due to the theory of images.
HD 67044 currently assigned a B8 spectral type is one of the slowly rotating B stars situated in the northern hemisphere which we are currently observing. The selection criteria for this sample of stars are a declination higher than -15^o , spectral class B8 or B9, luminosity class V or IV, and a magnitude V brighter than 7.85. Most of the stars of this B8-9 sample have just recently been observed in December 2014. We are currently performing a careful abundance analysis study of high resolution high S/N ratio spectra of these objects and sort them out into chemically normal stars (ie. whose abundances do not depart more than +- 0.15 dex from solar), new spectroscopic binaries and new chemically peculiar B stars (CPs) which had remained unoticed so far. We present here new abundance determinations for HD 67044 which allow us to propose that this star is a new CP late B star. Monier et al. (2015) have recently published the discovery of 4 new HgMn stars (3 from this late-B stars sample and one from a sample of 47 early A types stars verifying the same criteria). Royer et al. (2014) have published the analysis of the sample of 47 early A stars having low apparent projected velocities in the northern hemisphere up to V=6.65 mag. A careful abundance analysis of high resolution high S/N ratio spectra of these objects has sorted out the sample into 17 chemically normal stars, 12 spectroscopic binaries and 13 Chemically Peculiar stars (CPs) among which 5 are new CP stars.
A scheme for estimating atmospheric parameters T$_{eff}$, log$~g$, and [Fe/H] is proposed on the basis of Least Absolute Shrinkage and Selection Operator (LASSO) algorithm and Haar wavelet. The proposed scheme consists of three processes. A spectrum is decomposed using the Haar wavelet transform and low-frequency components at the fourth level are considered as candidate features. Then, spectral features from the candidate features are detected using the LASSO algorithm to estimate the atmospheric parameters. Finally, atmospheric parameters are estimated from the extracted spectral features using the support-vector regression (SVR) method. The proposed scheme was evaluated using three sets of stellar spectra respectively from Sloan Digital Sky Survey (SDSS), Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST), and Kurucz's model, respectively. The mean absolute errors are as follows: for 40~000 SDSS spectra, 0.0062 dex for log~T$_{eff}$ (85.83 K for T$_{eff}$), 0.2035 dex for log$~g$ and 0.1512 dex for [Fe/H]; for 23963 LAMOST spectra, 0.0074 dex for log~T$_{eff}$ (95.37 K for T$_{eff}$), 0.1528 dex for log~$g$, and 0.1146 dex for [Fe/H]; and for 10469 synthetic spectra, 0.0010 dex for log T$_{eff}$(14.42K for T$_{eff}$), 0.0123 dex for log~$g$, and 0.0125 dex for [Fe/H].
We revisit a classical perturbative approach to the Hamiltonian related to the motions of Trojan bodies, in the framework of the Planar Circular Restricted Three-Body Problem (PCRTBP), by introducing a number of key new ideas in the formulation. In some sense, we adapt the approach of Garfinkel (1977) to the context of the normal form theory and its modern techniques. First, we make use of Delaunay variables for a physically accurate representation of the system. Therefore, we introduce a novel manipulation of the variables so as to respect the natural behavior of the model. We develop a normalization procedure over the fast angle which exploits the fact that singularities in this model are essentially related to the slow angle. Thus, we produce a new normal form, i.e. an integrable approximation to the Hamiltonian. We emphasize some practical examples of the applicability of our normalizing scheme, e.g. the estimation of the stable libration region. Finally, we compare the level curves produced by our normal form with surfaces of section provided by the integration of the non--normalized Hamiltonian, with very good agreement. Further precision tests are also provided. In addition, we give a step-by-step description of the algorithm, allowing for extensions to more complicated models.
We carry out a series of simulations of G2-type clouds interacting with the black hole Sgr A* at the galactic center. We determine that the accretion rate from the gas cloud onto Sgr A* for a range of simulation parameters, such as cloud structure, background structure, background density, grid resolution, and accretion radius. Regardless of the numerical considerations, the amount of cloud material accreted is small, both compared to the total cloud mass and the normal background accretion rate. The accretion rate will remain small for at least 30 years after periapsis. We also model Br-gamma, bolometric, and X-ray emission from the cloud with a variety of density profiles, and compare to observational data. In simulations, the bolometric and X-ray luminosity have a peak in luminosity lasting from about 1 year before to 1 year after periapsis, a feature not detected in observations. Simulated Br-gamma emission remains nearly flat with a small peak 1 month to 1 year before periapsis, depending on how centrally concentrated the cloud is. Br-gamma emission decreases rapidly as the cloud passes periapsis due to shock heating of the gas. Observations show an increase of the FWHM of the Br-gamma line velocity dispersion leading up to periapsis, consistent with our simulations. Reproducing observations of G2 likely requires two components for the object: an extended, cold gas cloud responsible for the Br-gamma emission, and a compact core or dusty stellar object dominating the bolometric luminosity. Any emission from a gaseous component of G2 should be undetectable by 1 year after periapsis, due to shock heating and expansion of the cloud. Any remaining emission should be from the compact component of G2.
The detection of 1.37$\, $s pulsations from NuSTAR J095551+6940.8 implies the existence of an accreting pulsar, which challenges the conventional understanding of ultraluminous X-ray source. This kind of sources are proposed to be massive pulsars in this paper. Considering the general relativistic effect, stronger gravity of massive pulsars could lead to a larger maximum luminosity, scaled as the Eddington luminosity. The pseudo-Newtonian potential is employed to simulate the gravitational field in general relativity, and the Eddington luminosity is calculated for self-bound stars (quark star and quark-cluster star) and for the Tolman IV solution. It is found that, for a massive pulsar with radius close to the Schwarszchild radius, the Eddington luminosity could be as high as $2\times10^{39} \, {\rm erg\, s}^{-1}$. It is able to account for the X-ray luminosity of NuSTAR J095551+6940.8 with reasonable beaming factor. It is also suggested that massive pulsar-like compact stars could form via this super-Eddington phase of ultraluminous X-ray source.
Was the Milky Way galaxy a low-luminosity active galactic nucleus (AGN) in the past? Can we find traces of remnant structures supporting this idea? What is the three-dimensional arrangement of matter around our central supermassive black hole? A number of fundamental questions concerning our own Galactic center remain controversial. To reveal the structure of the high-energy sky around our galactic core, a technique more sensitive to the morphology of the emitters than spectroscopy is needed. In this lecture note, I describe how X-ray polarimetry can open a new observational window by precisely measuring the three-dimensional position of the scattering material in the Galactic center. The observed polarization degree and polarization position angle would also determine unambiguously the primary source of emission and trace the centennial history of our supermassive black hole by detecting echoes of its past activity thanks to astrophysical mirrors. Finally, the synergy between X-ray polarimetry and infrared and radio observations can be used to better constrain the geometry in the first hundred parsecs around the Galactic center.
String Theory and Supergravity allow, in principle, to follow the transition of the inflaton from pre-inflationary fast roll to slow roll. This introduces an infrared depression in the primordial power spectrum that might have left an imprint in the CMB anisotropy, if it occurred at accessible wavelengths. We model the effect extending $\Lambda$CDM with a scale $\Delta$ related to the infrared depression and explore the constraints allowed by {\sc Planck} data, employing also more conservative, wider Galactic masks in the low resolution CMB likelihood. In an extended mask with $f_{sky}=39\%$, we thus find $\Delta = (0.351 \pm 0.114) \times 10^{-3} \, \mbox{Mpc}^{-1}$, at $99.4\%$ confidence level, to be compared with a nearby value at $88.5\%$ with the standard $f_{sky}=94\%$ mask. With about 65 $e$--folds of inflation, these values for $\Delta$ would translate into primordial energy scales ${\cal O}(10^{14})$ GeV.
[Abridged] We investigate and apply a method to incorporate signal models that allow an additional frequency harmonic in searches for gravitational waves from spinning neutron stars, and particularly known pulsars. We assume emission is given by the general model of Jones (2010), which describes a triaxial star with non-aligned core and crust angular momenta, whose waveform under certain conditions reduces to that of a biaxial precessing star, or a simple triaxial star with aligned core and crust. The first two models can produce emission at both the star's rotation frequency (f) and twice the rotation frequency (2f), whilst the latter only produces emission at 2f, and has been the standard model for previous searches. We performed Monte-Carlo simulations generating both noise and signals at f and 2f to assess detection efficiency and investigate how well we can distinguish between the three different models (assuming a source with a known position and frequency evolution) by evaluating the Bayesian evidence for each model given the data. We have found that for signal-to-noise ratios $\gtrsim 6$ there is no significant loss in efficiency if performing a search for a signal at both f and 2f when the source is only producing emission at 2f. However, for sources with emission at both f and 2f signals could be missed by a search only at 2f depending on the distribution of power between harmonics. We also find that for a triaxial aligned source, the triaxial aligned model is always favoured, but for a triaxial non-aligned source it can be hard to distinguish between the triaxial non-aligned model and the biaxial model, even at high SNR. Finally, we apply the method to a selection of known pulsars using data from the LIGO fifth science run. From this we produce the first upper limits on gravitational wave amplitude at both f and 2f and apply the model selection criteria on real data.
We investigate 3D atmosphere dynamics for tidally locked terrestrial planets
with an Earth-like atmosphere and irradiation for different rotation periods
($P_{rot}=1-100$ days) and planet sizes ($R_P=1-2 R_{Earth}$) with
unprecedented fine detail. We could precisely identify three climate state
transition regions that are associated with phase transitions in standing
tropical and extra tropical Rossby waves.
We confirm that the climate on fast rotating planets may assume multiple
states ($P_{rot}\leq 12$ days for $R_P=2 R_{Earth}$). Our study is, however,
the first to identify the type of planetary wave associated with different
climate states: The first state is dominated by standing tropical Rossby waves
with fast equatorial superrotation. The second state is dominated by standing
extra tropical Rossby waves with high latitude westerly jets with slower wind
speeds. For very fast rotations ($P_{rot}\leq 5$~days for $R_P=2 R_{Earth}$),
we find another climate state transition, where the standing tropical and extra
tropical Rossby wave can both fit on the planet. Thus, a third state with a
mixture of the two planetary waves becomes possible that exhibits three jets.
Different climate states may be observable, because the upper atmosphere's
hot spot is eastward shifted with respect to the substellar point in the first
state, westward shifted in the second state and the third state shows a
longitudinal 'smearing' of the spot across the substellar point.
We show, furthermore, that the largest fast rotating planet in our study
exhibits atmosphere features known from hot Jupiters like fast equatorial
superrotation and a temperature chevron in the upper atmosphere.
The planetary exospheres are poorly known in their outer parts, since the neutral densities are low compared with the instruments detection capabilities. The exospheric models are thus often the main source of information at such high altitudes. We present a new way to take into account analytically the additional effect of the stellar radiation pressure on planetary exospheres. In a series of papers, we present with an Hamiltonian approach the effect of the radiation pressure on dynamical trajectories, density profiles and escaping thermal flux. Our work is a generalization of the study by Bishop and Chamberlain (1989). In this third paper, we investigate the effect of the stellar radiation pressure on the Circular Restricted Three Body Problem (CR3BP), called also the photogravitational CR3BP, and its implication on the escape and the stability of planetary exospheres, especially for Hot Jupiters. In particular, we describe the transformation of the equipotentials and the location of the Lagrange points, and we provide a modified equation for the Hill sphere radius that includes the influence of the radiation pressure. Finally, an application to the hot Jupiter HD 209458b reveals the existence of a blow-off escape regime induced by the stellar radiation pressure.
Trojans are small bodies in planetary Lagrangian points. In our solar system, Jupiter has the largest number of such companions. Their existence is assumed for exoplanetary systems as well, but none has been found so far. We present an analysis by super-stacking $\sim4\times10^4$ Kepler planets with a total of $\sim9\times10^5$ transits, searching for an average trojan transit dip. Our result gives an upper limit to the average Trojan transiting area (per planet) corresponding to one body of radius $<460$km at $2\sigma$ confidence. We find a significant Trojan-like signal in a sub-sample for planets with more (or larger) Trojans for periods $>$60 days. Our tentative results can and should be checked with improved data from future missions like PLATO2.0, and can guide planetary formation theories.
In the framework of the ERTBP, we study an example of the influence of secondary resonances over the long term stability of Trojan motions. By the integration of ensembles of orbits, we find various types of chaotic diffusion, slow and fast. We show that the distribution of escape times is bi-modular, corresponding to two populations of short and long escape times. The objects with long escape times produce a power-law tail in the distribution.
In this work, we study the motions in the region around the equilateral Lagrangian equilibrium points L4 and L5, in the framework of the Circular Planar Restricted Three-Body Problem (hereafter, CPRTBP). We design a semi-analytic approach based on some ideas by Garfinkel in [4]: the Hamiltonian is expanded in Poincar\'e-Delaunay coordinates and a suitable average is performed. This allows us to construct (quasi) invariant tori that are moderately far from the Lagrangian points L4-L5 and approximate wide tadpole orbits. This construction provides the tools for studying optimal transfers in the neighborhood of the equilateral points, when instantaneous impulses are considered. We show some applications of the new averaged Hamiltonian for the Earth-Moon system, applied to the setting-up of some transfers which allow to enter in the stability region filled by tadpole orbits.
We measure the horizontal subsurface flow in a fast emerging active region (NOAA 11158) using the ring-diagram technique and the HMI high-spatial resolution Dopplergrams. This active region had a complex magnetic structure and displayed significant changes in the morphology during its disk passage. Over the period of six days from 2011 February 11 to 16, the temporal variation in the magnitude of total velocity is found to follow the trend of magnetic field strength. We further analyze regions of individual magnetic polarity within AR 11158 and find that the horizontal velocity components in these sub-regions have significant variation with time and depth. The leading and trailing polarity regions move faster than the mixed-polarity region. Further, both zonal and meridional components have opposite signs for trailing and leading polarity regions at all depths showing divergent flows within the active region. We also find a sharp decrease in the magnitude of total horizontal velocity in deeper layer around major flares. It is suggested that the re-organization of magnetic fields during flares combined with the sunspot rotation decreases the magnitude of horizontal flows or that the flow kinetic energy has been converted into the energy released by flares. After the decline in flare activity and the sunspot rotation, the flows tend to follow the pattern of the magnetic activity.We also observe less variation in the velocity components near the surface but these tend to increase with depth, further demonstrating that the deeper layers are more affected by the topology of active regions.
(abridged) Luminous blue compact galaxies are among the most active galaxies in the local universe in terms of their star formation rate per unit mass. They may be seen as the local analogs of higher redshift Lyman Break Galaxies. Studies of their kinematics is key to understanding what triggers their unusually active star formation In this work we investigate the kinematics of stars and ionised gas in Haro11, one of the most luminous blue compact galaxies in the local universe. Previous works have indicated that many such galaxies may be triggered by galaxy mergers. We have employed Fabry-Perot interferometry, long-slit spectroscopy and Integral Field Unit (IFU) spectroscopy to explore the kinematics of Haro11. We target the near infrared Calcium triplet to derive the stellar velocity field and velocity dispersion. Ionised gas is analysed through emission lines from hydrogen, [OIII] , and [SIII]. When spectral resolution and signal to noise allows we investigate the the line profile in detail and identify multiple velocity components when present. We find that to first order, the velocity field and velocity dispersions derived from stars and ionised gas agree. Hence the complexities reveal real dynamical disturbances providing further evidence for a merger in Haro11. Through decomposition of emission lines we find evidence for kinematically distinct components, for instance a tidal arm behind the galaxy. The ionised gas velocity field can be traced to large galactocentric radii, and shows significant velocity dispersion even far out in the halo. We discuss the origin of the line width, and interpreted as virial motions it indicates a mass of ~1E11 M_sun. Morphologically and kinematically Haro11 shows many resemblances with the famous Antennae galaxies, but is much denser which is the likely explanation for the higher star formation efficiency in Haro11.
In this note we propose a model independent framework for inflationary (p)reheating. Our approach is analogous to the Effective Field Theory of Inflation, however here the inflaton oscillations provide an additional source of (discrete) symmetry breaking. Using the Goldstone field that non-linearly realizes time diffeormorphism invariance we construct a model independent action for both the inflaton and reheating sectors. Utilizing the hierarchy of scales present during the reheating process we are able to recover known results in the literature in a simpler fashion, including the presence of oscillations in the primordial power spectrum. We also construct a class of models where the shift symmetry of the inflaton is preserved during reheating, which helps alleviate past criticisms of (p)reheating in models of Natural Inflation. Extensions of our framework suggest the possibility of analytically investigating non-linear effects (such as rescattering and back-reaction) during thermalization without resorting to lattice methods. By construction, the EFT relates the strength of many of these interactions to other operators in the theory, including those responsible for the efficiency of (p)reheating. We conclude with a discussion of the limitations and challenges for our approach.
We propose a scenario of brane cosmology in which the Peccei-Quinn field plays the role of the inflaton and solves simultaneously many cosmological and phenomenological issues such as the generation of a heavy Majorana mass for the right-handed neutrinos needed for seesaw mechanism, MSSM $\mu$-parameter, the right amount of baryon number asymmetry and dark matter relic density at the present universe, together with an axion solution to the strong CP problem without the domain wall obstacle. Interestingly, the scales of the soft SUSY-breaking mass parameter and that of the breaking of $U(1)_{\rm PQ}$ symmetry are lower bounded at $\mathcal{O}(10) {\mathrm TeV}$ and $\mathcal{O}(10^{11}) {\mathrm GeV}$, respectively.
We consider small perturbations about homogeneous backgrounds in dilatationally-invariant Galileon models. The issues we address are stability (absence of ghosts and gradient instabilities) and superluminality. We show that in Minkowski background, it is possible to construct the Lagrangian in such a way that any homogeneous Galileon background solution is stable and small perturbations about it are subluminal. On the other hand, in the case of FLRW backgrounds, for any Lagrangian functions there exist homogeneous background solutions to the Galileon equation of motion and time-dependence of the scale factor, such that the stability conditions are satisfied, but the Galileon perturbations propagate with superluminal speed. Thus, a popular class of the generalized Galileon models is plagued by superluminality.
In terms of Sturm's theorem, we reexamine a marginal stable circular orbit (MSCO) such as the innermost stable circular orbit (ISCO) of a timelike geodesic in any spherically symmetric and static spacetime. MSCOs for some of exact solutions to the Einstein's equation are discussed. Strum's theorem is explicitly applied to the Kottler (often called Schwarzschild-de Sitter) spacetime. Moreover, we analyze MSCOs for a spherically symmetric, static and vacuum solution in Weyl conformal gravity.
We investigate spherically symmetric cosmological models in Einstein-aether theory with a tilted (non-comoving) perfect fluid source. We use a 1+3 frame formalism and adopt the comoving aether gauge to derive the evolution equations, which form a well-posed system of first order partial differential equations in two variables. We then introduce normalized variables. The formalism is particularly well-suited for numerical computations and the study of the qualitative properties of the models, which are also solutions of Horava gravity. We study the local stability of the equilibrium points of the resulting dynamical system corresponding to physically realistic inhomogeneous cosmological models and astrophysical objects with values for the parameters which are consistent with current constraints. In particular, we consider dust models in ($\beta-$) normalized variables and derive a reduced (closed) evolution system and we obtain the general evolution equations for the spatially homogeneous Kantowski-Sachs models using appropriate bounded normalized variables. We then analyse these models, with special emphasis on the future asymptotic behaviour for different values of the parameters. Finally, we investigate static models for a mixture of a (necessarily non-tilted) perfect fluid with a barotropic equations of state and a scalar field.
We study physical processes around a rotating black hole in pure Gauss-Bonnet (GB) gravity. In pure GB gravity, gravitational potential has slower fall off as compared to the corresponding Einstein potential in the same dimension. It is therefore expected that the energetics of pure GB black hole would be weaker, and our analysis bears out that the efficiency of energy extraction by Penrose process is increased to $25.8\%$ and particle acceleration is increased to $55.28\%$, and optical shadow of the black hole is decreased. These are the distinguishing in principle observable features of pure GB black hole.
We present results from three-dimensional non-linear hydrodynamic simulations of a precession driven flow in cylindrical geometry. The simulations are motivated by a dynamo experiment currently under development at Helmholtz-Zentrum Dresden-Rossendorf (HZDR) in which the possibility of generating a magnetohydrodynamic dynamo will be investigated in a cylinder filled with liquid sodium and simultaneously rotating around two axes. In this study, we focus on the emergence of non-axisymmetric time-dependent flow structures in terms of inertial waves which - in cylindrical geometry - form so-called Kelvin modes. For a precession ratio ${\rm{Po}}=\Omega_p/\Omega_c=0.014$ the amplitude of the forced Kelvin mode reaches up to one fourth of the rotation velocity of the cylindrical container confirming that precession provides a rather efficient flow driving mechanism even at moderate values of ${\rm{Po}}$. More relevant for dynamo action might be free Kelvin modes with higher azimuthal wave number. These free Kelvin modes are triggered by non-linear interactions and may constitute a triadic resonance with the fundamental forced mode when the height of the container matches their axial wave lengths. Our simulations reveal triadic resonances at aspect ratios close to those predicted by the linear theory except around the primary resonance of the forced mode. In that regime we still identify various free Kelvin modes, however, all of them exhibit a retrograde drift around the symmetry axis of the cylinder and none of them can be assigned to a triadic resonance. The amplitudes of the free Kelvin modes always remain below the forced mode but may reach up to 6% of the of the container's angular velocity. The properties of the free Kelvin modes will be used in future simulations of the magnetic induction equation to investigate their ability to provide for dynamo action.
We study the excitation of axial quasi-normal modes of deformed non-rotating black holes by test-particles and we compare the associated gravitational wave signal with that expected in general relativity from a Schwarzschild black hole. Deviations from standard predictions are quantified by an effective deformation parameter, which takes into account deviations from both the Schwarzschild metric and the Einstein equations. We show that, at least in the case of non-rotating black holes, it is possible to test the metric around the compact object, in the sense that the measurement of the gravitational wave spectrum can constrain possible deviations from the Schwarzschild solution.
In an accelerating universe in General Relativity there is a maximum radius above which a shell of test particles cannot collapse, but is dispersed by the cosmic expansion. This radius could be used in conjunction with observations of large structures to constrain the equation of state of the universe. We extend the concept of turnaround radius to modified theories of gravity for which the gravitational slip is non-vanishing.
Links to: arXiv, form interface, find, astro-ph, recent, 1508, contact, help (Access key information)
We describe our implementation of a global-parameter optimizer and Square Root Information Filter (SRIF) into the asteroid-modelling software SHAPE. We compare the performance of our new optimizer with that of the existing sequential optimizer when operating on various forms of simulated data and actual asteroid radar data. In all cases, the new implementation performs substantially better than its predecessor: it converges faster, produces shape models that are more accurate, and solves for spin axis orientations more reliably. We discuss potential future changes to improve SHAPE's fitting speed and accuracy.
We report the results of a direct imaging survey of A- and F-type main sequence stars searching for giant planets. A/F stars are often the targets of surveys, as they are thought to have more massive giant planets relative to solar-type stars. However, most imaging is only sensitive to orbital separations $>$30 AU, where it has been demonstrated that giant planets are rare. In this survey, we take advantage of the high-contrast capabilities of the Apodizing Phase Plate coronagraph on NACO at the Very Large Telescope. Combined with optimized principal component analysis post-processing, we are sensitive to planetary-mass companions (2 to 12 $M_{\rm Jup}$) at Solar System scales ($\leq$30 AU). We obtained data on 13 stars in L'-band and detected one new companion as part of this survey: an M$6.0\pm0.5$ dwarf companion around HD 984. We re-detect low-mass companions around HD 12894 and HD 20385, both reported shortly after the completion of this survey. We use Monte Carlo simulations to determine new constraints on the low-mass ($<$80 $M_{\rm Jup}$) companion frequency, as a function of mass and separation. Assuming solar-type planet mass and separation distributions, normalized to the planet frequency appropriate for A-stars, and the observed companion mass-ratio distribution for stellar companions extrapolated to planetary masses, we derive a truncation radius for the planetary mass companion surface density of $<$135 AU at 95% confidence.
The final step of most large-scale structure analyses involves the comparison of power spectra or correlation functions to theoretical models. It is clear that the theoretical models have parameter dependence, but frequently the measurements and the covariance matrix depend upon some of the parameters as well. We show that a very simple interpolation scheme from an unstructured mesh allows for an efficient way to include this parameter dependence self-consistently in the analysis at modest computational expense. We describe two schemes for covariance matrices. The scheme which uses the geometric structure of such matrices performs roughly twice as well as the simplest scheme, though both perform very well.
Accreting black holes produce powerful relativistic plasma jets which emit radiation across all observable wavelengths but the details of the initial acceleration and confinement of the jet are uncertain. We apply an innovative new model that allows us to determine key properties of the acceleration zone via multi-frequency observations. The central component of the model is a relativistic steady-state fluid flow, and the emission from physically distinct regions can be seen to contribute to different energy bands in the overall spectrum. By fitting with unprecedented accuracy to 42 simultaneous multiwavelength blazar spectra we are able to constrain the location of the brightest synchrotron emitting region, and show that there must be a linear relation between the jet power and the radius of the brightest region of the jet. We also find a correlation between the length of the accelerating region and the maximum bulk Lorentz factor of the jet and find evidence for a bimodal distribution of accretion rates in the blazar population. This allows us to put constraints on the basic dynamical and structural properties of blazar jets and to understand the underlying physical differences which result in the blazar sequence.
The feedback from massive stars is important to super star cluster (SSC) evolution and the timescales on which it occurs. SSCs form embedded in thick material, and eventually, the cluster is cleared out and revealed at optical wavelengths -- however, this transition is not well understood. We are investigating this critical SSC evolutionary transition with a multi-wavelength observational campaign. Although previously thought to appear after the cluster has fully removed embedding natal material, we have found that SSCs may host large populations of Wolf-Rayet stars. These evolved stars provide ionization and mechanical feedback that we hypothesize is the tipping point in the combined feedback processes that drive a SSC to emerge. Utilizing optical spectra obtained with the 4m Mayall Telescope at Kitt Peak National Observatory and the 6.5m MMT, we have compiled a sample of embedded SSCs that are likely undergoing this short-lived evolutionary phase and in which we confirm the presence of Wolf-Rayet stars. Early results suggest that WRs may accelerate the cluster emergence.
We present optical and near-infrared light curves and optical spectra of SN 2013dx, associated with the nearby (redshift 0.145) gamma-ray burst GRB 130702A. The prompt isotropic gamma-ray energy released from GRB 130702A is measured to be $E_{\gamma,\mathrm{iso}} = 6.4_{-1.0}^{+1.3} \times 10^{50}$erg (1keV to 10MeV in the rest frame), placing it intermediate between low-luminosity GRBs like GRB 980425/SN 1998bw and the broader cosmological population. We compare the observed $g^{\prime}r^{\prime}i^{\prime}z^{\prime}$ light curves of SN 2013dx to a SN 1998bw template, finding that SN 2013dx evolves $\sim20$% faster (steeper rise time), with a comparable peak luminosity. Spectroscopically, SN 2013dx resembles other broad-lined Type Ic supernovae, both associated with (SN 2006aj and SN 1998bw) and lacking (SN 1997ef, SN 2007I, and SN 2010ah) gamma-ray emission, with photospheric velocities around peak of $\sim$21,000 km s$^{-1}$. We construct a quasi-bolometric ($g^{\prime}r^{\prime}i^{\prime}z^{\prime}yJH$) light curve for SN 2013dx, and, together with the photospheric velocity, we derive basic explosion parameters using simple analytic models. We infer a $^{56}$Ni mass of $M_{\mathrm{Ni}} = 0.38\pm 0.01$M$_{\odot}$, an ejecta mass of $M_{\mathrm{ej}} = 3.0 \pm 0.1$ M$_{\odot}$, and a kinetic energy of $E_{\mathrm{K}} = (8.2 \pm 0.40) \times 10^{51}$erg (statistical uncertainties only), consistent with previous GRB-associated SNe. When considering the ensemble population of GRB-associated SNe, we find no correlation between the mass of synthesized $^{56}$Ni and high-energy properties, despite clear predictions from numerical simulations that $M_{\mathrm{Ni}}$ should correlate with the degree of asymmetry. On the other hand, $M_{\mathrm{Ni}}$ clearly correlates with the kinetic energy of the supernova ejecta across a wide range of core-collapse events.
We use N-body simulations to model the 12 Gyr evolution of a suite of star clusters with identical initial stellar mass functions over a range of initial cluster masses, sizes, and orbits. Our models reproduce the distribution of present-day global stellar mass functions that is observed in the Milky Way globular cluster population. We find that the slope of a star cluster's stellar mass function is strongly correlated with the fraction of mass that the cluster has lost, independent of the cluster's initial mass, and nearly independent of its orbit and initial size. Thus, the mass function - initial mass relation can be used to determine a Galactic cluster's initial total stellar mass, if the initial stellar mass function is known. We apply the mass function - initial mass relation presented here to determine the initial stellar masses of 33 Galactic globular clusters, assuming an universal Kroupa initial mass function. Our study suggests that globular clusters had initial masses that were on average a factor of 10 times larger than their present day mass, with seven clusters showing evidence for initial total stellar masses > 10^7 Msun.
We show that diffusion due to chaotic mixing in the Neighbourhood of the Sun may not be as relevant as previously suggested in erasing phase space signatures of past Galactic accretion events. For this purpose, we analyse Solar Neighbourhood-like volumes extracted from cosmological simulations that naturally account for chaotic orbital behaviour induced by the strongly triaxial and cuspy shape of the resulting dark matter haloes, among other factors. In the approximation of an analytical static triaxial model, our results show that a large fraction of stellar halo particles in such local volumes have chaos onset times (i.e., the timescale at which stars commonly associated with chaotic orbits will exhibit their chaotic behaviour) significantly larger than a Hubble time. Furthermore, particles that do present a chaotic behaviour within a Hubble time do not exhibit significant diffusion in phase space.
We study Modified Gravity (MG) theories by modelling the redshifted matter power spectrum in a spherical Fourier-Bessel (sFB) basis. We use a fully non-linear description of the real-space matter power-spectrum and include the lowest-order redshift-space correction (Kaiser effect), taking into account some additional non-linear contributions. Ignoring relativistic corrections, which are not expected to play an important role for a shallow survey, we analyse two different modified gravity scenarios, namely the generalised Dilaton scalar-tensor theories and the $f({R})$ models in the large curvature regime. We compute the 3D power spectrum ${\cal C}^s_{\ell}(k_1,k_2)$ for various such MG theories with and without redshift space distortions, assuming precise knowledge of background cosmological parameters. Using an all-sky spectroscopic survey with Gaussian selection function $\varphi(r)\propto \exp(-{r^2 / r^2_0})$, $r_0 = 150 \, h^{-1} \, {\textrm{Mpc}}$, and number density of galaxies $\bar {\textrm{N}} =10^{-4}\;{\textrm{Mpc}}^{-3}$, we use a $\chi^2$ analysis, and find that the lower-order $(\ell \leq 25)$ multipoles of ${\cal C}^s_\ell(k,k')$ (with radial modes restricted to $k < 0.2 \, h \,{\textrm{Mpc}}^{-1}$) can constraint the parameter $f_{R_0}$ at a level of $2\times 10^{-5} (3\times 10^{-5})$ with $3 \sigma$ confidence for $n=1(2)$. Combining constraints from higher $\ell > 25$ modes can further reduce the error bars and thus in principle make cosmological gravity constraints competitive with solar system tests. However this will require an accurate modelling of non-linear redshift space distortions. Using a tomographic $\beta(a)$-$m(a)$ parameterization we also derive constraints on specific parameters describing the Dilaton models of modified gravity.
The recent ALMA observations of the disc surrounding HL Tau reveal a very complex dust spatial distribution. We present a radiative transfer model accounting for the observed gaps and bright rings as well as radial changes of the emissivity index. We find that the dust density is depleted by at least a factor 10 in the main gaps compared to the surrounding rings. Ring masses range from 10-100 M$_{\oplus}$ in dust, and, we find that each of the deepest gaps is consistent with the removal of up to 40 M$_{\oplus}$ of dust. If this material has accumulated into rocky bodies, these would be close to the point of runaway gas accretion. Our model indicates that the outermost ring is depleted in millimetre grains compared to the central rings. This suggests faster grain growth in the central regions and/or radial migration of the larger grains. The morphology of the gaps observed by ALMA - well separated and showing a high degree of contrast with the bright rings over all azimuths - indicates that the millimetre dust disc is geometrically thin (scale height $\approx$ 1 au at 100 au) and that a large amount of settling of large grains has already occurred. Assuming a standard dust settling model, we find that the observations are consistent with a turbulent viscosity coefficient of a few $10^{-4}$. We estimate the gas/dust ratio in this thin layer to be of the order of 5 if the initial ratio is 100. The HCO$^+$ and CO emission is consistent with gas in Keplerian motion around a 1.7 $M_\odot$ star at radii from $\leq 10 - 120\,$au.
We present a complete census of all 263 Herschel-detected sources within the HST Frontier Fields (HFF), a deep multi-filter HST programme covering six massive lensing clusters. We provide a robust legacy catalogue of Herschel fluxes, primarily based on imaging from the Herschel Lensing Survey (HLS) and PEP/HerMES Key Programmes. Photometry is derived via a simultaneous PSF-fit using priors from archival Spitzer imaging. We optimally combine Herschel, Spitzer and WISE infrared (IR) photometry with data from HST, VLA and ground-based observatories, identifying optical counterparts to gain source redshifts. Hence for each Herschel-detected source we also present magnification factor (mu), intrinsic IR luminosity and characteristic dust temperature, providing a comprehensive view of dust-obscured star formation within the HFF. We demonstrate the utility of our catalogues through an exploratory overview of HST morphologies for the IR-bright population. In particular we briefly describe the highest redshift (z>2.5) and most magnified (mu>4) sources in the gravitationally lensed background.
Coronal seismology is extensively used to estimate properties of the corona,
e.g. the coronal magnetic field strength are derived from oscillations observed
in coronal loops. We present a three-dimensional coronal simulation including a
realistic energy balance in which we observe oscillations of a loop in
synthesised coronal emission. We use these results to test the inversions based
on coronal seismology.
From the simulation of the corona above an active region we synthesise
extreme ultraviolet (EUV) emission from the model corona. From this we derive
maps of line intensity and Doppler shift providing synthetic data in the same
format as obtained from observations. We fit the (Doppler) oscillation of the
loop in the same fashion as done for observations to derive the oscillation
period and damping time.
The loop oscillation seen in our model is similar to imaging and
spectroscopic observations of the Sun. The velocity disturbance of the kink
oscillation shows an oscillation period of 52.5s and a damping time of 125s,
both being consistent with the ranges of periods and damping times found in
observation. Using standard coronal seismology techniques, we find an average
magnetic field strength of $B_{\rm kink}=79$G for our loop in the simulation,
while in the loop the field strength drops from some 300G at the coronal base
to 50G at the apex. Using the data from our simulation we can infer what the
average magnetic field derived from coronal seismology actually means. It is
close to the magnetic field strength in a constant cross-section flux tube that
would give the same wave travel time through the loop.
Our model produced not only a realistic looking loop-dominated corona, but
also provides realistic information on the oscillation properties that can be
used to calibrate and better understand the result from coronal seismology.
The HARPS spectrograph is showing an extreme stability close to the m s-1 level over more than ten years of data. However the radial velocities of some stars are contaminated by a spurious one-year signal with an amplitude that can be as high as a few m s-1 . This signal is in opposition of phase with the revolution of Earth around the Sun, and can be explained by the deformation of spectral lines crossing block stitchings of the CCD when the spectrum of an observed star is alternatively blue- and red-shifted due to the motion of Earth around the Sun. This annual perturbation can be supress by either removing those affected spectral lines from the correlation mask used by the cross correlation technique to derive precise radial velocities, or by simply fitting a yearly sinusoid to the RV data. This is mandatory if we want to detect long-period low-amplitude signals in the HARPS radial velocities of quiet solar-type stars.
We present results from the first solar full-disk F$_{10.7}$ (the radio flux at $10.7$ cm, $2.8$ GHz) image taken with the S-band receivers on the recently upgraded Karl G. Jansky Very Large Array (VLA) in order to assess the relationship between the F$_{10.7}$ index and solar extreme ultra-violet (EUV) emission. To identify the sources of the observed $2.8$ GHz emission, we calculate differential emission measures (DEMs) from EUV images collected by the Atmospheric Imaging Assembly (AIA) and use them to predict the bremsstrahlung component of the radio emission. By comparing the bremsstrahlung prediction and radio observation we find that $8.1\pm 0.5\%$ of the variable component of the F$_{10.7}$ flux is associated with the gyroresonance emission mechanism. Additionally, we identify optical depth effects on the radio limb which may complicate the use of F$_{10.7}$ time series as an EUV proxy. Our analysis is consistent with a coronal iron abundance that is $4$ times the photospheric level.
We present a detailed study of the early phases of the peculiar supernova 2011ay based on BVRI photometry obtained at Konkoly Observatory, Hungary, and optical spectra taken with the Hobby-Eberly Telescope at McDonald Observatory, Texas. The spectral analysis carried out with SYN++ and SYNAPPS confirms that SN 2011ay belongs to the recently defined class of SNe Iax, which is also supported by the properties of its light and color curves. The estimated photospheric temperature around maximum light, T_{phot} ~8,000 K, is lower than in most Type Ia SNe, which results in the appearance of strong Fe II features in the spectra of SN 2011ay, even during the early phases. We also show that strong blending with metal features (those of Ti II, Fe II, Co II) makes the direct analysis of the broad spectral features very difficult, and this may be true for all SNe Iax. We find two alternative spectrum models that both describe the observed spectra adequately, but their photospheric velocities differ by at least 3,000 km/s. The quasi-bolometric light curve of SN~2011ay has been assembled by integrating the UV-optical spectral energy distributions. Fitting a modified Arnett-model to L_{bol}(t), the moment of explosion and other physical parameters, i.e. the rise time to maximum, the ^{56}Ni mass and the total ejecta mass are estimated as t_{rise} ~14 +/-1 days, M_{Ni} ~0.22 +/- 0.01 M_{sol} and M_{ej} ~0.8 M_{sol}, respectively.
The richness of dynamical behavior exhibited by the rotational states of various solar system objects has driven significant advances in the theoretical understanding of their evolutionary histories. An important factor that determines whether a given object is prone to exhibiting non-trivial rotational evolution is the extent to which such an object can maintain a permanent aspheroidal shape, meaning that exotic behavior is far more common among the small body populations of the solar system. Gravitationally bound binary objects constitute a substantial fraction of asteroidal and TNO populations, comprising systems of triaxial satellites that orbit permanently deformed central bodies. In this work, we explore the rotational evolution of such systems with specific emphasis on quadrupole-quadrupole interactions, and show that for closely orbiting, highly deformed objects, both prograde and retrograde spin-spin resonances naturally arise. Subsequently, we derive capture probabilities for leading order commensurabilities and apply our results to the illustrative examples of (87) Sylvia and (216) Kleopatra asteroid systems. Cumulatively, our results suggest that spin-spin coupling may be consequential for highly elongated, tightly orbiting binary objects.
We study the fate of our Universe assuming that the present accelerated stage is due to a scalar field in a linear potential. Such a Universe would bounce and collapse in the future. We solve numerically and analytically the equations of motion for the scalar field and the scale factor. In particular, we relate the duration of the accelerated stage, the bounce and the collapse with the mass of the field and, thus, with the current value of equation of state $w$. We obtain an expression which predicts the age of the Universe for a given $w+1$. The present constraints on $w$ imply that the Universe will not collapse in the next $56$ billion years. Moreover, a cosmological solution to the coincidence problem favors a significant deviation of $w$ from $-1$ such that the Universe collapses in the not too distant future.
We report optical-infrared (IR) properties of faint 1.3 mm sources (S_1.3mm = 0.2-1.0 mJy) detected with the Atacama Large Millimeter/submillimeter Array (ALMA) in the Subaru/XMM-Newton Deep Survey (SXDS) field. We searched for optical/IR counterparts of 8 ALMA-detected sources (>=4.0 sigma, the sum of the probability of spurious source contamination is ~1) in a K-band source catalog. Four ALMA sources have K-band counterpart candidates within a 0.4" radius. Comparison between ALMA-detected and undetected K-band sources in the same observing fields shows that ALMA-detected sources tend to be brighter, more massive, and more actively forming stars. While many of the ALMA-identified submillimeter-bright galaxies (SMGs) in previous studies lie above the sequence of star-forming galaxies in stellar mass--star-formation rate plane, our ALMA sources are located in the sequence, suggesting that the ALMA-detected faint sources are more like `normal' star-forming galaxies rather than `classical' SMGs. We found a region where multiple ALMA sources and K-band sources reside in a narrow photometric redshift range (z ~ 1.3-1.6) within a radius of 5" (42 kpc if we assume z = 1.45). This is possibly a pre-merging system and we may be witnessing the early phase of formation of a massive elliptical galaxy.
We investigate the claim that the ratio {\beta} of radiation pressure force to gravitational force on a dust grain in our solar system can substantially exceed unity for some grain sizes, provided that grain porosity is high enough. For model grains consisting of random aggregates of silicate spherules, we find that the maximum value of {\beta} is almost independent of grain porosity, but for small (<0.3 {\mu}m) grains, {\beta} actually decreases with increasing porosity. These results affect the interpretation of the grain trajectories estimated from the Stardust mission, which were modeled assuming {\beta} values exceeding one. We find that radiation pressure effects are not large enough for particles Orion and Hylabrook captured by Stardust to be of interstellar origin given their reported impact velocities. We also investigate the effect of metallic iron inclusions in the dust grains, and find that metallic iron will increase {\beta}, but at least half the grain (by mass) must be iron in order to raise the maximum value of {\beta} as a function of grain size above one. We also investigate the effects of solar radiation on transverse velocities and grain spin, and show that radiation pressure introduces both transverse velocities and equatorial spin velocities of several hundred meters per second for incoming interstellar grains at 2 AU. Such spin rates may result in centrifugal disruption of aggregates.
We introduce a Bayesian solution to the problem of inferring the density profile of strong gravitational lenses when the lens galaxy may contain multiple dark or faint substructures. The source and lens models are based on a superposition of an unknown number of non-negative basis functions (or "blobs") whose form was chosen with speed as a primary criterion. The prior distribution for the blobs' properties is specified hierarchically, so the mass function of substructures is a natural output of the method. We use reversible jump Markov Chain Monte Carlo (MCMC) within Diffusive Nested Sampling (DNS) to sample the posterior distribution and evaluate the marginal likelihood of the model, including the summation over the unknown number of blobs in the source and the lens. We demonstrate the method on a simulated data set with a single substructure, which is recovered well with moderate uncertainties. We also apply the method to the g-band image of the "Cosmic Horseshoe" system, and find some hints of potential substructures. However, we caution that such results could also be caused by misspecifications in the model (such as the shape of the smooth lens component or the point spread function), which are difficult to guard against in full generality.
We have carried out a search for optically visible post-Asymptotic Giant Branch (post-AGB) stars in the Large Magellanic Cloud (LMC). First, we selected candidates with a mid-IR excess and then obtained their optical spectra. We disentangled contaminants with unique spectra such as M-stars, C-stars, planetary nebulae, quasi-stellar objects and background galaxies. Subsequently, we performed a detailed spectroscopic analysis of the remaining candidates to estimate their stellar parameters such as effective temperature, surface gravity (log g), metallicity ([Fe/H]), reddening and their luminosities. This resulted in a sample of 35 likely post-AGB candidates with late-G to late-A spectral types, low log g, and [Fe/H] < -0.5. Furthermore, our study confirmed the existence of the dusty post-Red Giant Branch (post-RGB) stars, discovered previously in our SMC survey, by revealing 119 such objects in the LMC. These objects have mid-IR excesses and stellar parameters (Teff, log g, [Fe/H]) similar to those of post-AGB stars except that their luminosities (< 2500 Lsun), and hence masses and radii, are lower. These post-RGB stars are likely to be products of binary interaction on the RGB. The post-AGB and post-RGB objects show SED properties similar to the Galactic post-AGB stars, where some have a surrounding circumstellar shell, while some others have a surrounding stable disc similar to the Galactic post-AGB binaries. This study also resulted in a new sample of 162 young stellar objects, identified based on a robust log g criterion. Other interesting outcomes include objects with an UV continuum and an emission line spectrum; luminous supergiants; hot main-sequence stars; and 15 B[e] star candidates, 12 of which are newly discovered in this study.
The life of a solar active prominence, one of the most remarkable objects on the Sun, is full of dynamics; after first appearing on the Sun the prominence continuously evolves with various internal motions and eventually produces a global eruption toward the interplane- tary space. Here we report that the whole life of an active prominence is successfully re- produced by performing as long-term a magnetohydrodynamic simulation of a magnetized prominence plasma as was ever done. The simulation reveals underlying dynamic processes that give rise to observed properties of an active prominence: invisible subsurface flows self- consistently produce the cancellation of magnetic flux observed at the photosphere, while observed and somewhat counterintuitive strong upflows are driven against gravity by en- hanced gas pressure gradient force along a magnetic field line locally standing vertical. The most highlighted dynamic event, transition into an eruptive phase, occurs as a natural con- sequence of the self-consistent evolution of a prominence plasma interacting with a magnetic field, which is obtained by seamlessly reproducing dynamic processes involved in the forma- tion and eruption of an active prominence.
X-ray fluorescence remote sensing technique plays a significant role in the chemical compositions research of the Moon. Here we describe the data analysis method for China's Chang'E-2 X-ray spectrometer (CE2XRS) in detail and present the preliminary results: the first global Mg/Si and Al/Si maps on the lunar surface. Our results show that the distributions of Mg/Si and Al/Si correlate well with the terrains of the Moon. The higher Mg/Si ratio corresponding to the mare regions while the lower value corresponding to the highland terrains. The map of Al/Si ratio shows a reverse relationship with the map of Mg/Si ratio.
SAGE1C\,J053634.78$-$722658.5 is a galaxy at redshift $z=0.14$, discovered behind the Large Magellanic Cloud in the {\it Spitzer} Space Telescope "Surveying the Agents of Galaxy Evolution" Spectroscopy survey (SAGE-Spec). It has very strong silicate emission at 10 $\mu$m but negligible far-IR and UV emission. This makes it a candidate for a bare AGN source in the IR, perhaps seen pole-on, without significant IR emission from the host galaxy. In this paper we present optical spectra taken with the Southern African Large Telescope (SALT) to investigate the nature of the underlying host galaxy and its AGN. We find broad H$\alpha$ emission characteristic of an AGN, plus absorption lines associated with a mature stellar population ($>9$ Gyr), and refine its redshift determination to $z=0.1428\pm0.0001$. There is no evidence for any emission lines associated with star formation. This remarkable object exemplifies the need for separating the emission from any AGN from that of the host galaxy when employing infrared diagnostic diagrams. We estimate the black hole mass, $M_{\rm BH}=3.5\pm0.8\times10^8$ M$_\odot$, host galaxy mass, $M_{\rm stars}=2.5^{2.5}_{1.2}\times10^{10}$ M$_\odot$, and accretion luminosity, $L_{\rm bol}({\rm AGN})=5.3\pm0.4\times10^{45}$ erg s$^{-1}$ ($\approx12$ per cent of the Eddington luminosity) and find the AGN to be more prominent than expected for a host galaxy of this modest size. The old age is in tension with the downsizing paradigm in which this galaxy would recently have transformed from a star-forming disc galaxy into an early-type, passively evolving galaxy.
I review the roles of jet-inflated bubbles in determining the evolution of different astrophysical objects. I discuss astrophysical systems where jets are known to inflate bubbles (cooling flow [CF] clusters; young galaxies; intermediate luminosity optical transients [ILOTs]; bipolar planetary nebulae [PNe]), and systems that are speculated to have jet-inflated bubbles (core collapse supernovae [CCSNe]; common envelope evolution [CEE]; grazing envelope evolution [GEE]). The jets in many of these cases act through a negative jet feedback mechanism (JFM). I discuss the outcomes when the JFM fizzle, or does not work at all. According to this perspective, some very interesting and energetic events owe their existence to the failure of the JFM, including stellar black holes, gamma ray bursts, and type Ia supernovae.
The available photometry from the online databases were used for the first light curve analysis of eight eclipsing binary systems EI Aur, XY Dra, BP Dra, DD Her, VX Lac, WX Lib, RZ Lyn, and TY Tri. All these stars are of Algol-type, having the detached components and the orbital periods from 0.92 to 6.8 days. For the systems EI Aur and BP Dra the large amount of the third light was detected during the light curve solution. Moreover, 468 new times of minima for these binaries were derived, trying to identify the period variations. For the systems XY Dra and VX Lac the third bodies were detected with the periods 17.7, and 49.3 years, respectively.
Observations of the WC9+OB system WR 65 in the infrared show variations of its dust emission consistent with a period near 4.8~yr, suggesting formation in a colliding-wind binary (CWB) having an elliptical orbit. If we adopt the IR maximum as zero phase, the times of X-ray maximum count and minimum extinction to the hard component measured by Oskinova & Hamann fall at phases 0.4--0.5, when the separation of the WC9 and OB stars is greatest. We consider WR 65 in the context of other WC8-9+OB stars showing dust emission.
The large-scale structure of the universe is comprised of virialized
blob-like clusters, linear filaments, sheet-like walls and huge near empty
three-dimensional voids. Characterizing the large scale universe is essential
to our understanding of the formation and evolution of galaxies. The density
range of clusters, walls and voids are relatively well separated, when compared
to filaments, which span a relatively larger range. The large scale filamentary
network thus forms an intricate part of the cosmic web.
In this paper, we describe Felix, a topology based framework for visual
exploration of filaments in the cosmic web. The filamentary structure is
represented by the ascending manifold geometry of the 2-saddles in the
Morse-Smale complex of the density field. We generate a hierarchy of
Morse-Smale complexes and query for filaments based on the density ranges at
the end points of the filaments. The query is processed efficiently over the
entire hierarchical Morse-Smale complex, allowing for interactive
visualization.
We apply Felix to computer simulations based on the heuristic Voronoi
kinematic model and the standard $\Lambda$CDM cosmology, and demonstrate its
usefulness through two case studies. First, we extract cosmic filaments within
and across cluster like regions in Voronoi kinematic simulation datasets. We
demonstrate that we produce similar results to existing structure finders.
Filaments that form the spine of the cosmic web, which exist in high density
regions in the current epoch, are isolated using Felix. Also, filaments present
in void-like regions are isolated and visualized. These filamentary structures
are often over shadowed by higher density range filaments and are not easily
characterizable and extractable using other filament extraction methodologies.
The Arcminute Microkelvin Imager Galactic Plane Survey (AMIGPS) provides mJy-sensitivity, arcminute-resolution interferometric images of the northern Galactic plane at $\approx$ 16 GHz. The first data release covered $76^{\circ} \lessapprox \ell \lessapprox 170^{\circ}$ between latitudes of $|b| \lessapprox 5^{\circ}$; here we present a second data release, extending the coverage to $53^{\circ} \lessapprox \ell \lessapprox 193^{\circ}$ and including high-latitude extensions to cover the Taurus and California giant molecular cloud regions, and the recently discovered large supernova remnant G159.6+7.3. The total coverage is now 1777 deg$^2$ and the catalogue contains 6509 sources. We also describe the improvements to the data processing pipeline which improves the positional and flux density accuracies of the survey.
Preliminary results are presented from spectroscopic data in the optical range of the Galactic ring nebulae NGC 6888, G2:4+1:4, RCW 58 and Sh2-308. Deep observations with long exposure times were carried out at the 6.5m Clay Telescope and at the 10.4m Gran Telescopio Canarias. In NGC 6888, recombination lines of C II, O II and N II are detected with signal-to-noise ratios higher than 8. The chemical content of NGC 6888 is discussed within the chemical enrichment predicted by evolution models of massive stars. For all nebulae, a forthcoming work will content in-depth details about observations, analysis and final results (Esteban et al. 2015, in prep.).
NIKA is a dual-band camera operating with 315 frequency multiplexed LEKIDs cooled at 100 mK. NIKA is designed to observe the sky in intensity and polarisation at 150 and 260 GHz from the IRAM 30-m telescope. It is a test-bench for the final NIKA2 camera. The incoming linear polarisation is modulated at four times the mechanical rotation frequency by a warm rotating multi-layer Half Wave Plate. Then, the signal is analysed by a wire grid and finally absorbed by the LEKIDs. The small time constant (< 1ms ) of the LEKID detectors combined with the modulation of the HWP enables the quasi-simultaneous measurement of the three Stokes parameters I, Q, U, representing linear polarisation. In this pa- per we present results of recent observational campaigns demonstrating the good performance of NIKA in detecting polarisation at mm wavelength.
Within the context of the collaboration "B fields in OB stars (BOB)", we used the FORS2 low-resolution spectropolarimeter to search for a magnetic field in 50 massive stars, including two reference magnetic massive stars. Because of the many controversies of magnetic field detections obtained with the FORS instruments, we derived the magnetic field values with two completely independent reduction and analysis pipelines. We compare and discuss the results obtained from the two pipelines. We obtained a general good agreement, indicating that most of the discrepancies on magnetic field detections reported in the literature are caused by the interpretation of the significance of the results (i.e., 3-4 sigma detections considered as genuine, or not), instead of by significant differences in the derived magnetic field values. By combining our results with past FORS1 measurements of HD46328, we improve the estimate of the stellar rotation period, obtaining P = 2.17950+/-0.00009 days. For HD125823, our FORS2 measurements do not fit the available magnetic field model, based on magnetic field values obtained 30 years ago. We repeatedly detect a magnetic field for the O9.7V star HD54879, the HD164492C massive binary, and the He-rich star CPD -57 3509. We obtain a magnetic field detection rate of 6+/-4%, while by considering only the apparently slow rotators we derive a detection rate of 8+/-5%, both comparable with what was previously reported by other similar surveys. We are left with the intriguing result that, although the large majority of magnetic massive stars is rotating slowly, our detection rate is not a strong function of the stellar rotational velocity.
We use a combination of simulated cosmological probes and astrophysical tests of the stability of the fine-structure constant $\alpha$, as expected from the forthcoming European Extremely Large Telescope (E-ELT), to constrain the class of string-inspired runaway dilaton models of Damour, Piazza and Veneziano. We consider three different scenarios for the dark sector couplings in the model and discuss the observational differences between them. We improve previously existing analyses investigating in detail the degeneracies between the parameters ruling the coupling of the dilaton field to the other components of the universe, and studying how the constraints on these parameters change for different fiducial cosmologies. We find that if the couplings are small (e.g., $\alpha_b=\alpha_V\sim0$) these degeneracies strongly affect the constraining power of future data, while if they are sufficiently large (e.g., $\alpha_b\gtrsim10^{-5}-\alpha_V\gtrsim0.05$, as in agreement with current constraints) the degeneracies can be partially broken. We show that E-ELT will be able to probe some of this additional parameter space.
To date, most simulations of the chemistry in protoplanetary discs have used 1+1D or 2D axisymmetric $\alpha$-disc models to determine chemical compositions within young systems. This assumption is inappropriate for non-axisymmetric, gravitationally unstable discs, which may be a significant stage in early protoplanetary disc evolution. Using 3D radiative hydrodynamics, we have modelled the physical and chemical evolution of a 0.17 M$_{\odot}$ self-gravitating disc over a period of 2000 yr. The 0.8 M$_{\odot}$ central protostar is likely to evolve into a solar-like star, and hence this Class 0 or early Class I young stellar object may be analogous to our early Solar System. Shocks driven by gravitational instabilities enhance the desorption rates, which dominate the changes in gas-phase fractional abundances for most species. We find that at the end of the simulation, a number of species distinctly trace the spiral structure of our relatively low-mass disc, particularly CN. We compare our simulation to that of a more massive disc, and conclude that mass differences between gravitationally unstable discs may not have a strong impact on the chemical composition. We find that over the duration of our simulation, successive shock heating has a permanent effect on the abundances of HNO, CN and NH$_3$, which may have significant implications for both simulations and observations. We also find that HCO$^+$ may be a useful tracer of disc mass. We conclude that gravitational instabilities induced in lower mass discs can significantly, and permanently, affect the chemical evolution, and that observations with high-resolution instruments such as ALMA offer a promising means of characterising gravitational instabilities in protosolar discs.
We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT source catalog (2FGL) and was among the top ten most significant unassociated sources in the recent third catalog (3FGL). PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of $8.9$ Hz, a characteristic age of $49$ kyr, and spin-down power $1.0 \times 10^{36}$ erg s$^{-1}$. In 2009 August it suffered one of the largest glitches detected from a gamma-ray pulsar ($\Delta f / f \approx 4.5\times10^{-6}$). Remaining undetected in dedicated radio follow-up observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this.
Neutrinos play a crucial role in the collapse and explosion of massive stars, governing the infall dynamics of the stellar core, triggering and fueling the explosion and driving the cooling and deleptonization of the newly formed neutron star. Due to their role neutrinos carry information from the heart of the explosion and, due to their weakly interacting nature, offer the only direct probe of the dynamics and thermodynamics at the center of a supernova. In this paper, we review the present status of modelling the neutrino physics and signal formation in collapsing and exploding stars. We assess the capability of current and planned large underground neutrino detectors to yield faithful information of the time and flavor dependent neutrino signal from a future Galactic supernova. We show how the observable neutrino burst would provide a benchmark for fundamental supernova physics with unprecedented richness of detail. Exploiting the treasure of the measured neutrino events requires a careful discrimination of source-generated properties from signal features that originate on the way to the detector. As for the latter, we discuss self-induced flavor conversions associated with neutrino-neutrino interactions that occur in the deepest stellar regions; matter effects that modify the pattern of flavor conversions in the dynamical stellar envelope; neutrino-oscillation signatures that result from structural features associated with the shock-wave propagation as well as turbulent mass motions in post-shock layers. Finally, we highlight our current understanding of the formation of the diffuse supernova neutrino background and we analyse the perspectives for a detection of this relic signal that integrates the contributions from all past core-collapse supernovae in the Universe.
We analyse the main physical parameters and the circumstellar environment of the young Herbig Be star HD 98922. We present AMBER/VLTI high spectral resolution (R =12000) interferometric observations across the Br$\gamma$ line, accompanied by UVES high-resolution spectroscopy and SINFONI-AO assisted near-infrared integral field spectroscopic data. To interpret our observations, we develop a magneto-centrifugally driven disc-wind model. Our analysis of the UVES spectrum shows that HD 98922 is a young (~5x10^5 yr) Herbig Be star (SpT=B9V), located at a distance of 440(+60-50) pc, with a mass accretion rate of ~9+/-3x10^(-7) M_sun yr^(-1). SINFONI K-band AO-assisted imaging shows a spatially resolved circumstellar disc-like region (~140 AU in diameter) with asymmetric brightness distribution. Our AMBER/VLTI UT observations indicate that the Br$\gamma$ emitting region (radius ~0.31+/-0.04 AU) is smaller than the continuum emitting region (inner dust radius ~0.7+/-0.2 AU), showing significant non-zero V-shaped differential phases (i.e. non S-shaped, as expected for a rotating disc). The value of the continuum-corrected pure Br$\gamma$ line visibility at the longest baseline (89 m) is ~0.8+/-0.1, i.e. the Br$\gamma$ emitting region is partially resolved. Our modelling suggests that the observed Br$\gamma$ line-emitting region mainly originates from a disc wind with a half opening angle of 30deg, and with a mass-loss rate of ~2x10(-7) M_sun yr^(-1). The observed V-shaped differential phases are reliably reproduced by combining a simple asymmetric continuum disc model with our Br$\gamma$ disc-wind model. The Br$\gamma$ emission of HD 98922 can be modelled with a disc wind that is able to approximately reproduce all interferometric observations if we assume that the intensity distribution of the dust continuum disc is asymmetric.
The Cherenkov Telescope Array (CTA) is a large collaborative effort aimed at the design and operation of an observatory dedicated to very high-energy gamma-ray astrophysics in the energy range from a few tens of GeV to above 100 TeV, which will yield about an order of magnitude improvement in sensitivity with respect to the current major arrays (H.E.S.S., MAGIC, and VERITAS). Within this framework, the Italian National Institute for Astrophysics is leading the ASTRI project, whose main goals are the design and installation on Mt. Etna (Sicily) of an end-to-end dual-mirror prototype of the CTA small size telescope (SST) and the installation at the CTA Southern site of a dual-mirror SST mini-array composed of nine units with a relative distance of about 300 m. The innovative dual-mirror Schwarzschild-Couder optical solution adopted for the ASTRI Project allows us to substantially reduce the telescope plate-scale and, therefore, to adopt silicon photo-multipliers as light detectors. The ASTRI mini-array is a wider international effort. The mini-array, sensitive in the energy range 1-100 TeV and beyond with an angular resolution of a few arcmin and an energy resolution of about 10-15%, is well suited to study relatively bright sources (a few $\times 10^{-12}$erg cm$^{-2}$s$^{-1}$ at 10 TeV) at very high energy. Prominent sources such as extreme blazars, nearby well-known BL Lac objects, Galactic pulsar wind nebulae, supernovae remnants, micro-quasars, and the Galactic Center can be observed in a previously unexplored energy range. The ASTRI mini-array will extend the current IACTs sensitivity well above a few tens of TeV and, at the same time, will allow us to compare our results on a few selected targets with those of current (HAWC) and future high-altitude extensive air-shower detectors.
A 1-parameter class of quadratic equations of state is confronted with the Hubble diagram of supernovae. The fit is found to be as good as the one using the standard LambdaCDM model. However this quadratic equation of state precludes objects with redshifts higher than z_max = 1.7. Adding a fair amount of cold baryons to the model increases z_max without spoiling the fit.
Deep ultraviolet (DUV) light sources are used to neutralise isolated test masses in highly sensitive space-based gravitational experiments. An example is the LISA Pathfinder charge management system, which uses low-pressure mercury lamps. A future gravitational-wave observatory such as eLISA will use UV light-emitting diodes (UV LEDs), which offer numerous advantages over traditional discharge lamps. Such devices have limited space heritage but are are now available from a number of commercial suppliers. Here we report on a test campaign that was carried out to quantify the general properties of three types of commercially available UV LEDs and demonstrate their suitability for use in space. Testing included general electrical and UV output measurements, spectral stability, pulsed performance, temperature dependence as well as thermal vacuum, radiation and vibration survivability.
Slowly pulsating B stars (SPB) and $\gamma$ Dor stars pulsate in high-order gravity (g-) modes. The frequencies of g-modes are sensitive to the detailed structure and evolution history of stars having convective cores. Receding convective cores in OB-type stars leave behind a chemically inhomogenous $\nabla_\mu>0$ radiative zone. Once a g-mode has radial nodes near the boundaries of these layers, the mode gets trapped and its period deviates from asymptotic period spacing. Careful study of such trapped modes allows constraining the extent of such layers by fitting individual pulsation frequencies. We employ 19 consecuitve dipole g-modes of a very rich Kepler SPB pulsator, KIC 10526294, to demonstrate the power of mode trapping in B-stars in studying the thermal and chemical stratification in the overshooting layer.
Black holes that accrete far below the Eddington limit are believed to accrete through a geometrically thick, optically thin, rotationally supported plasma that we will refer to as a radiatively inefficient accretion flow (RIAF). RIAFs are typically collisionless in the sense that the Coulomb mean free path is large compared to $GM/c^2$, and relativistically hot near the event horizon. In this paper we develop a phenomenological model for the plasma in RIAFs, motivated by the application to sources such as Sgr A* and M87. The model is derived using Israel-Stewart theory, which considers deviations up to second order from thermal equilibrium, but modified for a magnetized plasma. This leads to thermal conduction along magnetic field lines and a difference in pressure, parallel and perpendicular to the field lines (which is equivalent to anisotrotropic viscosity). In the non-relativistic limit, our model reduces to the widely used Braginskii theory of magnetized, weakly collisional plasmas. We compare our model to the existing literature on dissipative relativistic fluids, describe the linear theory of the plasma, and elucidate the physical meaning of the free parameters in the model. We also describe limits of the model when the conduction is saturated and when the viscosity implies a large pressure anisotropy. In future work, the formalism developed in this paper will be used in numerical models of RIAFs to assess the importance of non-ideal processes for the dynamics and radiative properties of slowly accreting black holes.
We present initial results of a new campaign of simulations focusing on the interaction of planetary winds with stellar environments using Adaptive Mesh Refinement methods. We have confirmed the results of Stone & Proga 2009 that an azimuthal flow structure is created in the planetary wind due to day/night temperatures differences. We show that a backflow towards the planet will occur with a strength that depends on the escape parameter. When a stellar outflow is included, we see unstable bow waves forming through the outflow's interaction with the planetary wind.
In a systematic study of data from a large number of nuclei with 50 < A < 248 (most taken with spin-0 targets) photon strength, level density and finally radiative capture processes are investigated with respect to their sensitivity on nuclear structure and especially nuclear deformation. A parameterization is proposed using a fit to GDR shapes by the sum of three Lorentzians; this TLO ansatz is based on the breaking of axiality and is global as it used one set of parameters for a range in nuclear mass number A from 50 to 250. This is reached by using calculated ground state deformations, obtained on the basis of a mapped collective Hamiltonian and the Gogny D1S interaction. To test the calculations they will also be compared to multiple Coulomb excitation results explicitly examining axial deformation and triaxiality in heavy nuclei. As level densities are also influenced strongly by deviations from sphericity, causing a significant collective enhancement, a coherent look at both, level density and photon strength, is presented. This finally allows a comparison with experimental Maxwellian averaged cross sections for neutron capture by even target nuclei from Se to Cm. The focus on spin-0 targets has the advantage of no additional spin-factors to be involved. Only one free global parameter enters, when ground state masses are used to quantify shell effects. The proposed parameterization of the electric dipole strength can also be considered global as again only two free fit parameters suffice to obtain good descriptions for all the nuclei, including odd isotopes. Thus a small set of global parameters tests the combination of photon strength and level density information and a reliable prediction for compound nuclear reactions also outside the valley of stability is expected.
We propose that the observed matter-antimatter asymmetry can be naturally produced as a byproduct of axion-driven slow-roll inflation by coupling the axion to standard-model neutrinos. We assume that GUT scale right-handed neutrinos are responsible for the masses of the standard model neutrinos and that the Higgs is a light field during inflation and develops a Hubble scale vacuum expectation value (VEV). In this set up, the rolling axion generates a helicity asymmetry in standard-model neutrinos. Following inflation, this helicity asymmetry becomes equal to a net lepton number as the Higgs VEV decays and is partially re-processed by the $SU(2)_{L}$ sphaleron into a net baryon number.
We study derivatively coupled fermions in axion-driven inflation, specifically $m_\phi^2\phi^2$ and monodromy inflation, and calculate particle production during the inflationary epoch and the post-inflationary axion oscillations. During inflation, the rolling axion acts as an effective chemical potential for helicity which biases the gravitational production of one fermion helicity over the other. This mechanism allows for efficient gravitational production of heavy fermion states that would otherwise be highly suppressed. Following inflation, the axion oscillates and fermions with both helicities are produced as the effective frequency of the fermion field changes non-adiabatically. For certain values of the fermion mass and axion-fermion coupling strength, the two helicity states are produced asymmetrically, resulting in unequal number-densities of left- and right-helicity fermions.
Links to: arXiv, form interface, find, astro-ph, recent, 1508, contact, help (Access key information)
Ground-based spectrometers have been developed to measure the concentration, velocity, and temperature of ozone in the mesosphere and lower thermosphere (MLT) using low-cost satellite television electronics to observe the 11.072 GHz spectral line of ozone. A two-channel spectrometer has been engineered to yield various performance improvements, including a doubling of the signal-to-noise ratio, improved data processing efficiency, and lower power consumption at 15 W. Following 2009 and 2012 observations of the seasonal and diurnal variations in ozone concentration near the mesopause, the ozone line was observed at an altitude near 95 km and latitude of 38 degrees north using three single-channel spectrometers located at the MIT Haystack Observatory (Westford, MA), Chelmsford High School (Chelmsford, MA), and Union College (Schenectady, NY) pointed south at 8 degrees. Observations from 2009 through 2014 are used to derive the nightly-averaged seasonal variation in meridional velocity, as well as the seasonally-averaged variation with local solar time. The results indicate a seasonal trend in which the winds at 95 km come from the north at about $10\,\text{m}\text{s}^{-1}$ in the summer of the northern hemisphere, and from the south at about $10\,\text{m}\text{s}^{-1}$ in the winter. Nighttime data from -5 to +5 hours local solar time show a gradual transition of the meridional wind velocity from about -$20\,\text{m}\text{s}^{-1}$ to +$20\,\text{m}\text{s}^{-1}$. These two trends correlate with nighttime wind measurements from the Millstone Hill High-Resolution F\'abry-Perot Interferometer (FPI) in Westford, MA, which uses the 557.7 nm green line nightglow from atomic oxygen centered at 95 km. The results have also been compared with average meridional winds measured with meteor radar.
Data re-sampling methods such as the delete-one jackknife are a common tool
for estimating the covariance of large scale structure probes. In this paper we
investigate the concepts of internal covariance estimation in the context of
cosmic shear two-point statistics. We demonstrate how to use log-normal
simulations of the convergence field and the corresponding shear field to carry
out realistic tests of internal covariance estimators and find that most
estimators such as jackknife or sub-sample covariance can reach a satisfactory
compromise between bias and variance of the estimated covariance.
In a forecast for the complete, 5-year DES survey we show that internally
estimated covariance matrices can provide a large fraction of the true
uncertainties on cosmological parameters in a 2D cosmic shear analysis. The
volume inside contours of constant likelihood in the $\Omega_m$-$\sigma_8$
plane as measured with internally estimated covariance matrices is on average
$\gtrsim 85\%$ of the volume derived from the true covariance matrix. The
uncertainty on the parameter combination $\Sigma_8 \sim \sigma_8
\Omega_m^{0.5}$ derived from internally estimated covariances is $\sim 90\%$ of
the true uncertainty.
The non-Gaussian nature of the Epoch of Reionization (EoR) 21-cm signal has a significant impact on the error variance of its power spectrum ${P_{\rm b}}({\bf k})$ (Mondal et al., 2015). Building on the previous work, we have used a large ensemble of semi-numerical simulations and an analytical model to estimate the effect of this non-Gaussianity on the entire error covariance matrix ${\mathcal{C}}_{ij}$. Our analytical model shows that ${\mathcal{C}}_{ij}$ has contributions from two sources. One is the usual variance for a Gaussian random field which scales inversely of the number of modes that goes into the estimation of ${P_{\rm b}}({\bf k})$. The other is the trispectrum of the signal. Using the simulated 21-cm signal ensemble, an ensemble of the randomized signal and ensembles of Gaussian random ensembles we have quantified the effect of the trispectrum on the error variance ${\mathcal{C}}_{ij}$. We find that its relative contribution is comparable to or larger than that of the Gaussian term for the $k$ range $0.3 \leq k \leq 1.0 \,{\rm Mpc}^{-1}$, and can be even $\sim 200$ times larger at $k \sim 5\, {\rm Mpc}^{-1}$. We also establish that the off-diagonal terms of ${\mathcal{C}}_{ij}$ have statistically significant non-zero values which arise purely from the trispectrum. This further signifies that the error in different $k$ modes are not independent. We find a strong correlation between the errors at large $k$ values $(\ge 0.5 \,{\rm Mpc}^{-1})$, and a weak correlation between the smallest and largest $k$ values. There is also a small anti-correlation between the errors in the smallest and intermediate $k$ values. These results are relevant for the $k$ range that will be probed by the current and upcoming EoR 21-cm experiments.
We present a self-consistent, absolute isochronal age scale for young (< 200 Myr), nearby (< 100 pc) moving groups, which is consistent with recent lithium depletion boundary ages for both the Beta Pic and Tucana-Horologium moving groups. This age scale was derived using a set of semi-empirical pre-main-sequence model isochrones that incorporate an empirical colour-Teff relation and bolometric corrections based on the observed colours of Pleiades members, with theoretical corrections for the dependence on log g. Absolute ages for young, nearby groups are vital as these regions play a crucial role in our understanding of the early evolution of low- and intermediate-mass stars, as well as providing ideal targets for direct imaging and other measurements of dusty debris discs, substellar objects and, of course, extrasolar planets.
The Kepler mission was launched in 2009 and has discovered thousands of
planet candidates. In a recent paper, Esteves et al. (2013) found a periodic
signal in the light curves of KOI-13 and HAT-P-7, with a frequency triple the
orbital frequency of a transiting planet. We found similar harmonics in many
systems with a high occurrence rate. At this time, the origins of the signal
are not entirely certain.
We look carefully at the possibility of errors being introduced through our
data processing routines but conclude that the signal is real. The harmonics on
multiples of the orbital frequency are a result of non-sinusoidal periodic
signals. We speculate on their origin and generally caution that these
harmonics could lead to wrong estimates of planet albedos, beaming mass
estimates, and ellipsoidal variations.
Radiation hydrodynamics (RHD) simulations are used to study many astrophysical phenomena, however they require the use of simplified radiation transport and thermal prescriptions to reduce computational cost. In this paper we present a systematic study of the importance of microphysical processes in RHD simulations using the example of D-type HII region expansion. We compare the simplest hydrogen-only models with those that include: ionisation of H, He, C, N, O, S and Ne, different gas metallicity, non-LTE metal line blanketed stellar spectral models of varying metallicity, radiation pressure, dust and treatment of photodissociation regions. Each of these processes are explicitly treated using modern numerical methods rather than parameterisation. In line with expectations, changes due to microphysics in either the effective number of ionising photons or the thermal structure of the gas lead to differences in D-type expansion. In general we find that more realistic calculations lead to the onset of D-type expansion at smaller radii and a slower subsequent expansion. Simulations of star forming regions using simplified microphysics are therefore likely overestimating the strength of radiative feedback. We find that both variations in gas metallicity and the inclusion of dust can affect the ionisation front evolution at the 10-20 per cent level over 500kyr, which could substantially modify the results of simplified 3D models including feedback. Stellar metallicity, radiation pressure and the inclusion of photodissociation regions are all less significant effects at the 1 per cent level or less, rendering them of minor importance in the modelling the dynamical evolution of HII regions.
In this work we present Spitzer 3.6 and 4.5 micron secondary eclipse observations of five new cool (<1200 K) transiting gas giant planets: HAT-P-19b, WASP-6b, WASP-10b, WASP-39b, and WASP-67b. We compare our measured eclipse depths to the predictions of a suite of atmosphere models and to eclipse depths for planets with previously published observations in order to constrain the temperature- and mass-dependent properties of gas giant planet atmospheres. We find that the dayside emission spectra of planets less massive than Jupiter require models with efficient circulation of energy to the night side and/or increased albedos, while those with masses greater than that of Jupiter are consistently best-matched by models with inefficient circulation and low albedos. At these relatively low temperatures we expect the atmospheric methane to CO ratio to vary as a function of metallicity, and we therefore use our observations of these planets to constrain their atmospheric metallicities. We find that the most massive planets have dayside emission spectra that are best-matched by solar metallicity atmosphere models, but we are not able to place strong constraints on metallicities of the smaller planets in our sample. Interestingly, we find that the ratio of the 3.6 and 4.5 micron brightness temperatures for these cool transiting planets is independent of planet temperature, and instead exhibits a tentative correlation with planet mass. If this trend can be confirmed, it would suggest that the shape of these planets' emission spectra depends primarily on their masses, consistent with the hypothesis that lower-mass planets are more likely to have metal-rich atmospheres.
We study the evolution of planet-induced vortices in radially stratified disks, with initial conditions allowing for radial buoyancy. For this purpose we run global two dimensional hydrodynamical simulations, using the PLUTO code. Planet-induced vortices are a product of the Rossby wave instability (RWI) triggered in the edges of a planetary gap. In this work we assess the influence of radial buoyancy for the development of the vortices. We found that radial buoyancy leads to smoother planetary gaps, which generates weaker vortices. This effect is less pronounced for locally isothermal and quasi-isothermal (very small cooling rate) disks. We observed the formation of two generations of vortices. The first generation of vortices is formed in the outer wall of the planetary gap. The merged primary vortex induces accretion, depleting the mass on its orbit. This process creates a surface density enhancement beyond the primary vortex position. The second generation of vortices arise in this surface density enhancement, indicating that the bump in this region is sufficient to trigger the RWI. The merged secondary vortex is a promising explanation for the location of the vortex in the Oph IRS48 system. Finally, we observed a nonmonotonic behavior for the vortex lifetimes as a function of the thermal relaxation timescale, agreeing with previous studies. The birth times of the secondary vortices also display a nonmonotonic behavior, which is correlated with the growth time of the primary vortex and its induced accretion.
We characterize magnetically driven accretion at radii between 1 au and 100 au in protoplanetary discs, using a series of local non-ideal magnetohydrodynamic (MHD) simulations. The simulations assume a Minimum Mass Solar Nebula (MMSN) disc that is threaded by a net vertical magnetic field of specified strength. Confirming previous results, we find that the Hall effect has only a modest impact on accretion at 30 au, and essentially none at 100 au. At 1-10 au the Hall effect introduces a pronounced bi-modality in the accretion process, with vertical magnetic fields aligned to the disc rotation supporting a strong laminar Maxwell stress that is absent if the field is anti-aligned. In the anti-aligned case, we instead find evidence for bursts of turbulent stress at 5-10 au, which we tentatively identify with the non-axisymmetric Hall-shear instability. The presence or absence of these bursts depends upon the details of the adopted chemical model, which suggests that appreciable regions of actual protoplanetary discs might lie close to the borderline between laminar and turbulent behaviour. Given the number of important control parameters that have already been identified in MHD models, quantitative predictions for disc structure in terms of only radius and accretion rate appear to be difficult. Instead, we identify robust qualitative tests of magnetically driven accretion. These include the presence of turbulence in the outer disc, independent of the orientation of the vertical magnetic fields, and a Hall-mediated bi-modality in turbulent properties extending from the region of thermal ionization to 10 au.
We present a collection of 250 anomalous Cepheids (ACs) discovered in the
OGLE-IV fields toward the Large (LMC) and Small Magellanic Cloud (SMC). The LMC
sample is an extension of the OGLE-III Catalog of ACs published in 2008, while
the SMC sample contains the first known bona fide ACs in this galaxy. The total
sample is composed of 141 ACs in the LMC and 109 ACs in the SMC. All these
stars pulsate in single modes: fundamental (174 objects) or first overtone (76
objects). Additionally, we report the discovery of four ACs located in the
foreground of the Magellanic Clouds. These are the first fundamental-mode ACs
known in the Galactic field.
We demonstrate that the coefficients phi_21 and phi_31 determined by the
Fourier light curve decomposition are useful discriminators between classical
Cepheids and ACs, at least in the LMC and in the field of the Milky Way. In the
SMC, the light curve shapes and mean magnitudes of short-period classical
Cepheids make them similar to ACs, which is a source of difficulties in the
discrimination of both classes of pulsators. The presence of unidentified ACs
in the catalogs of classical Cepheids may be partly responsible for the
observed non-linearity of the period-luminosity relation observed for
short-period Cepheids in the SMC. We compare spatial distributions of ACs,
classical Cepheids and RR Lyr stars. We show that the distribution of ACs
resembles that of old stars (RR Lyr variables), although in the LMC there are
visible structures typical for young population (classical Cepheids): the bar
and spiral arms. This may suggest that ACs are a mixture of relatively young
stars and mergers of very old stars.
Terrestrial planets are thought to be the result of a vast number of gravitational interactions and collisions between smaller bodies. We use numerical simulations to show that practically identical initial conditions result in a wide array of final planetary configurations. This highly chaotic behaviour questions the predictability of different scenarios for the formation and evolution of our solar system and planetary systems in general. However, multiple realisations of the same initial conditions can be used to predict certain global statistics. We present two sets of numerical experiments that quantify this behaviour. Firstly, we demonstrate that simulations with slightly displaced particles are completely divergent after ~500 years, irrespective of initial displacement, particle number, and code accuracy. If a single planetesimal is moved by less than one millimetre, then a different set of planets results -- this timescale for chaotic divergence decreases with increasing particle number. Secondly, we show final planetary configurations of initially similar simulations with and without giant planets after evolving them for ~148 Myr. We find that the same simulations including giant planets tend to generate higher mass planets at lower semi-major axes than simulations without gas giants. This prediction can be tested with forthcoming observational programs. By extracting outliers in the observations, we cautiously predict that Kepler-10, Kepler-9, 61 Vir, HD 134060, and HD 51608 may host as yet undetected giant planets.
With the availability of considerably more data, we revisit the question of how special our Solar System is, compared to observed exoplanetary systems. To this goal, we employ a mathematical transformation that allows for a meaningful, statistical comparison. We find that the masses and densities of the giant planets in our Solar System are very typical, as is the age of the Solar System. While the orbital location of Jupiter is somewhat of an outlier, this is most likely due to strong selection effects towards short-period planets. The eccentricities of the planets in our Solar System are relatively small compared to those in observed exosolar systems, but still consistent with the expectations for an 8-planet system (and could, in addition, reflect a selection bias towards high-eccentricity planets). The two characteristics of the Solar System that we find to be most special are the lack of super-Earths with orbital periods of days to months and the general lack of planets inside of the orbital radius of Mercury. Overall, we conclude that in terms of its broad characteristics our Solar System is not expected to be extremely rare, allowing for a level of optimism in the search for extrasolar life.
Asteroseismology of bright stars with well-determined properties from parallax measurements and interferometry can yield precise stellar ages and meaningful constraints on the composition. We substantiate this claim with an updated asteroseismic analysis of the solar-analog binary system 16 Cyg A & B using the complete 30-month data sets from the Kepler space telescope. An analysis with the Asteroseismic Modeling Portal (AMP), using all of the available constraints to model each star independently, yields the same age ($t=7.0 \pm 0.3$ Gyr) and composition ($Z=0.021 \pm 0.002$, $Y_i=0.25 \pm 0.01$) for both stars, as expected for a binary system. We quantify the accuracy of the derived stellar properties by conducting a similar analysis of a Kepler-like data set for the Sun, and we investigate how the reliability of asteroseismic inference changes when fewer observational constraints are available or when different fitting methods are employed. We find that our estimates of the initial helium mass fraction are probably biased low by 0.02-0.03 from neglecting diffusion and settling of heavy elements, and we identify changes to our fitting method as the likely source of small shifts from our initial results in 2012. We conclude that in the best cases reliable stellar properties can be determined from asteroseismic analysis even without independent constraints on the radius and luminosity.
We study the relationships between galaxy environments and galaxy properties related to disk (re)growth, considering two highly complete samples that are approximately baryonic mass limited into the high-mass dwarf galaxy regime, the Environmental COntext (ECO) catalog (data release herein) and the B-semester region of the REsolved Spectroscopy Of a Local VolumE (RESOLVE) survey. We quantify galaxy environments using both group identification and smoothed galaxy density field methods. We use by-eye and quantitative morphological classifications plus atomic gas content measurements and estimates. We find that blue early-type (E/S0) galaxies, gas-dominated galaxies, and UV-bright disk host galaxies all become distinctly more common below group halo mass ~10^11.5 Msun, implying that this low group halo mass regime may be a preferred regime for significant disk growth activity. We also find that blue early-type and blue late-type galaxies inhabit environments of similar group halo mass at fixed baryonic mass, consistent with a scenario in which blue early types can regrow late-type disks. In fact, we find that the only significant difference in the typical group halo mass inhabited by different galaxy classes is for satellite galaxies with different colors, where at fixed baryonic mass red early and late types have higher typical group halo masses than blue early and late types. More generally, we argue that the traditional morphology-environment relation (i.e., that denser environments tend to have more early types) can be largely attributed to the morphology-galaxy mass relation for centrals and the color-environment relation for satellites.
We report the results of a long-term spectroscopic monitoring of the FS\,CMa type object MWC\,728. We found that it is a binary system with a B5 Ve (T$_{\rm eff}$ = 14000$\pm$1000 K) primary and a G8 III type (T$_{\rm eff} \sim$ 5000 K) secondary. Absorption line positions of the secondary vary with a semi-amplitude of $\sim$20 km/s and a period of 27.5 days. The system's mass function is 2.3$\times10^{-2}$ M$_\odot$, and its orbital plane is $13^{\circ}-15^{\circ}$ tilted from the plane of the sky. The primary's $v \sin i \sim$110 km/s combined with this tilt implies that it rotates at a nearly breakup velocity. We detected strong variations of the Balmer and He I emission-line profiles on timescales from days to years. This points to a variable stellar wind of the primary in addition to the presence of a circum-primary gaseous disk. The strength of the absorption-line spectrum along with the optical and near-IR continuum suggest that the primary contributes $\sim$60% of the $V$--band flux, the disk contributes $\sim$30%, and the secondary $\sim$10%. The system parameters, along with the interstellar extinction, suggest a distance of $\sim$1 kpc, that the secondary does not fill its Roche lobe, and that the companions' mass ratio is $q \sim$0.5. Overall, the observed spectral variability and the presence of a strong IR-excess are in agreement with a model of a close binary system that has undergone a non-conservative mass-transfer.
We present the results of extensive multi-band intra-night optical monitoring
of BL Lacertae during 2010--2012. BL Lacertae was very active in this period
and showed intense variability in almost all wavelengths. We extensively
observed it for a total for 38 nights; on 26 of them observations were done
quasi-simultaneously in B, V, R and I bands (totaling 113 light curves), with
an average sampling interval of around 8 minutes. BL Lacertae showed
significant variations on hour-like timescales in a total of 19 nights in
different optical bands. We did not find any evidence for periodicities or
characteristic variability time-scales in the light curves.
The intranight variability amplitude is generally greater at higher
frequencies and decreases as the source flux increases.
We found spectral variations in BL Lacertae in the sense that the optical
spectrum becomes flatter as the flux increases but in several flaring states
deviates from the linear trend suggesting different jet components contributing
to the emission at different times.
The frequency ratios $r_{01}$ and $r_{10}$ of KIC 11081729 decrease firstly and then increase with frequency. These characteristics of the ratios can be directly reproduced by the models with the overshooting parameter $\delta_{\rm ov}$ between 1.7 and 1.8 and derive from the effects of both overshooting and the acoustic glitch between radiative and overshooting regions. The effects of the glitch on the ratios can be described by a function of sine, which depends on the frequency $\nu_{0}$ of the mode whose inner turning point is just located on the boundary between the radiative and overshooting regions. The value of the $\nu_{0}$ decreases with the increase in the radius of the overshooting region. The radius of the overshooting region increases with the increase in $\delta_{\rm ov}$. Thus the distributions of the ratios are sensitive to the value of $\delta_{\rm ov}$. The characteristics of the ratios provide a strict constraint on stellar models and aid in determining the size of the overshooting region. Observational constraints favor a star with $M=1.23\pm0.02$ $M_{\odot}$, $R=1.41\pm0.01$ $R_{\odot}$, $L=3.07\pm0.12$ $L_{\odot}$, $t=4.4\pm0.8$ Gyr, and $\delta_{\rm ov}$ between 1.7 and 1.8.
We observed the prototype blazar, BL Lacertae, extensively in optical and
radio bands during an active phase in the period 2010--2013 when the source
showed several prominent outbursts. We searched for possible correlations and
time lags between the optical and radio band flux variations using
multifrequency data to learn about the mechanisms producing variability. During
an active phase of BL Lacertae, we searched for possible correlations and time
lags between multifrequency light curves of several optical and radio bands. We
tried to estimate any possible variability timescales and inter-band lags in
these bands. We performed optical observations in B, V, R and I bands from
seven telescopes in Bulgaria, Georgia, Greece and India and obtained radio data
at 36.8, 22.2, 14.5, 8 and 4.8 GHz frequencies from three telescopes in
Ukraine, Finland and USA. Significant cross-correlations between optical and
radio bands are found in our observations with a delay of cm-fluxes with
respect to optical ones of ~250 days. The optical and radio light curves do not
show any significant timescales of variability. BL Lacertae showed many optical
'mini-flares' on short time-scales. Variations on longer term timescales are
mildly chromatic with superposition of many strong optical outbursts. In radio
bands, the amplitude of variability is frequency dependent. Flux variations at
higher radio frequencies lead the lower frequencies by days or weeks.
The optical variations are consistent with being dominated by a geometric
scenario where a region of emitting plasma moves along a helical path in a
relativistic jet. The frequency dependence of the variability amplitude
supports an origin of the observed variations intrinsic to the source.
Refraction and diffraction of incoming radio waves by the ionosphere induce time variability in the angular positions, peak amplitudes and shapes of radio sources, potentially complicating the automated cross-matching and identification of transient and variable radio sources. In this work, we empirically assess the effects of the ionosphere on data taken by the Murchison Widefield Array (MWA) radio telescope. We directly examine 51 hours of data observed over 10 nights under quiet geomagnetic conditions (global storm index Kp < 2), analysing the behaviour of short-timescale angular position and peak flux density variations of around ten thousand unresolved sources. We find that while much of the variation in angular position can be attributed to ionospheric refraction, the characteristic displacements (10-20 arcsec) at 154 MHz are small enough that search radii of 1-2 arcmin should be sufficient for cross-matching under typical conditions. By examining bulk trends in amplitude variability, we place upper limits on the modulation index associated with ionospheric scintillation of 1-3% for the various nights. For sources fainter than ~1 Jy, this variation is below the image noise at typical MWA sensitivities. Our results demonstrate that the ionosphere is not a significant impediment to the goals of time-domain science with the MWA at 154 MHz.
The ratios $r_{01}$ and $r_{10}$ of small to large separations of KIC 2837475 primarily exhibit an increase behavior in the observed frequency range. The calculations indicate that only the models with overshooting parameter $\delta_{\rm ov}$ between approximately 1.2 and 1.6 can reproduce the observed ratios $r_{01}$ and $r_{10}$ of KIC 2837475. The ratios $r_{01}$ and $r_{10}$ of the frequency separations of p-modes with inner turning points that are located in the overshooting region of convective core can exhibit an increase behavior. The frequencies of the modes that can reach the overshooting region decrease with the increase in $\delta_{\rm ov}$. Thus the ratio distributions are more sensitive to $\delta_{\rm ov}$ than to other parameters. The increase behavior of the KIC 2837475 ratios results from a direct effect of the overshooting of convective core. The characteristic of the ratios provides a strict constraint on stellar models. Observational constraints point to a star with $M=1.490\pm0.018$ $M_{\odot}$, $R=1.67\pm0.01$ $R_{\odot}$, age $=2.8\pm0.4$ Gyr, and $1.2\lesssim$ $\delta_{\rm ov}$ $\lesssim1.6$ for KIC 2837475.
The origin of the seeds which develop into the observed super-massive black holes at high redshifts may be hard to interpret in the context of the standard $\Lambda CDM$ of early universe cosmology based on Gaussian primordial perturbations. Here we consider the modification of the halo mass function obtained by introducing skewness and kurtosis of the primordial fluctuations. We show that such primordial non-Gaussianities constrained by the current observational bounds on the nonlinearity parameters of $f_{NL}$ and $g_{NL}$ are not effective at greatly increasing the number density of seeds which could develop into super-massive black holes at high redshifts. This is to be contrasted with the role which cosmic string loops could play in seeding super-massive black holes.
In our solar system, Mars-sized protoplanets frequently collided with each other during the last stage of terrestrial planet formation called the giant impact stage. Giant impacts eject a large amount of material from the colliding protoplanets into the terrestrial planet region, which may form debris disks with observable infrared excesses. Indeed, tens of warm debris disks around young solar-type stars have been observed. Here, we quantitatively estimate the total mass of ejected materials during the giant impact stages. We found that $\sim$0.4 times the Earth's mass is ejected in total throughout the giant impact stage. Ejected materials are ground down by collisional cascade until micron-sized grains are blown out by radiation pressure. The depletion timescale of these ejected materials is determined primarily by the mass of the largest body among them. We conducted high-resolution simulations of giant impacts to accurately obtain the mass of the largest ejected body. We then calculated the evolution of the debris disks produced by a series of giant impacts and depleted by collisional cascades to obtain the infrared excess evolution of the debris disks. We found that the infrared excess is almost always higher than the stellar infrared flux throughout the giant impact stage ($\sim$100 Myr) and is sometimes $\sim$10 times higher immediately after a giant impact. Therefore, giant impact stages would explain the infrared excess from most observed warm debris disks. The observed fraction of stars with warm debris disks indicates that the formation probability of our solar system-like terrestrial planets is approximately 10%.
I review experimental and observational constraints on a possible non-minimal coupling of a scalar field to electromagnetism (dilatonic coupling). Such a coupling is motivated from recent quasar spectrum observations that indicate a possible spatial and/or temporal variation of the fine-structure constant. I consider a dilatonic coupling of the form $B_F(\phi)=1+g\phi$. The strongest bounds on $g$ come from weak equivalence principle tests which impose the constraint $g<1.6 \times 10^{-17} GeV^{-1}$. This constraint is strong enough to rule out this class of models as a cause for an observable cosmological variation of the fine structure constant unless chameleon type mechanism is present.
For MeV gamma-ray astronomy, we have developed an electron-tracking Compton camera (ETCC) as a MeV gamma-ray telescope capable of rejecting the radiation background and attaining the high sensitivity of near 1 mCrab in space. Our ETCC comprises a gaseous time-projection chamber (TPC) with a micro pattern gas detector for tracking recoil electrons and a position-sensitive scintillation camera for detecting scattered gamma rays. After the success of a first balloon experiment in 2006 with a small ETCC (using a 10$\times$10$\times$15 cm$^3$ TPC) for measuring diffuse cosmic and atmospheric sub-MeV gamma rays (Sub-MeV gamma-ray Imaging Loaded-on-balloon Experiment I; SMILE-I), a (30 cm)$^{3}$ medium-sized ETCC was developed to measure MeV gamma-ray spectra from celestial sources, such as the Crab Nebula, with single-day balloon flights (SMILE-II). To achieve this goal, a 100-times-larger detection area compared with that of SMILE-I is required without changing the weight or power consumption of the detector system. In addition, the handling events are also expected to dramatically increase during observation. Here, we describe both the concept and the performance of the new data-acquisition system with this (30 cm)$^{3}$ ETCC to manage 100 times more data while satisfying the severe restrictions regarding the weight and power consumption imposed by a balloon-borne observation. In particular, to improve the detection efficiency of the fine tracks in the TPC from $\sim$10\% to $\sim$100\%, we introduce a new data-handling algorithm in the TPC. Therefore, for efficient management of such large amounts of data, we developed a data-acquisition system with parallel data flow.
The aim of this work is to solve the dispersion relations near the first excitation threshold of photon propagating along the magnetic field in the strong field limit. We have calculated the time damping of the photon in two particular cases: the degenerate gas as well as the diluted gas limit being both important from the Astrophysical point of view. In particular the diluted gas limit could describe the magnethosphere of neutron stars. The solutions have been used to obtain a finite Quantum Faraday angle in both limits. A resonant behavior for the Faraday angle is also obtained. To reproduce the semi-classical result for the Faraday rotation angle the weak field limit is considered.
Azimuthally asymmetric dust distributions observed with ALMA in transition disks have been interpreted as dust traps. We present VLA Ka band (34 GHz or 0.9 cm) and ALMA Cycle 2 Band 9 (680 GHz or 0.45 mm) observations at 0.2" resolution of the Oph IRS 48 disk, which suggest that larger particles could be more azimuthally concentrated than smaller dust grains, assuming an axisymmetric temperature field or optically thin 680 GHz emission. Fitting an intensity model to both data demonstrates that the azimuthal extent of the millimeter emission is 2.3 $\pm0.9$ times as wide as the centimeter emission, marginally consistent with the particle trapping mechanism under the above assumptions. The 34 GHz continuum image also reveals evidence for ionized gas emission from the star. Both the morphology and the spectral index variations are consistent with an increase of large particles in the center of the trap, but uncertainties remain due to the continuum optical depth at 680 GHz. Particle trapping has been proposed in planet formation models to allow dust particles to grow beyond millimeter sizes in the outer regions of protoplanetary disks. The new observations in the Oph IRS 48 disk provide support for the dust trapping mechanism for centimeter-sized grains, although additional data is required for definitive confirmation.
The recent discoveries of ultra-faint dwarf (UFD) galaxies in the vicinity of the Magellanic system supports the expectation from cosmological models that such faint objects exist and are numerous. By developing a mass model of the Local Group and backwards integrating the Magellanic Clouds' present kinematics, we find that the locations of these UFDs are consistent with those predicted if previously associated with the Large MC as part of a loose association. We further demonstrate how these satellites are likely to have been processed by the Galactic hot halo upon accretion, with the implication that ongoing detections of extremely gas-rich objects on the periphery of the Galaxy and without clear stellar counterparts are analogous to the progenitors of the gas-deficient UFDs. Our model allows us predict the locations of other putative Magellanic satellites, and propose how their distribution/kinematics provide a novel constraint on the dynamical properties of the Galaxy. We also predict that the stripped metal-poor HI, previously associated with these UFDs, lies coincident with but distinguishable from the extensive Magellanic Stream.
The interactions between radio-loud AGN and their environments play an
important r\^{o}le in galaxy and cluster evolution. Recent work has
demonstrated fundamental differences between High and Low Excitation Radio
Galaxies (HERGs and LERGs), and shown that they may have different
relationships with their environments. In the Chandra Large Project ERA
(Environments of Radio-loud AGN), we made the first systematic X-ray
environmental study of the cluster environments of radio galaxies at a single
epoch (z~0.5), and found tentative evidence for a correlation between radio
luminosity and cluster X-ray luminosity. We also found that this relationship
appeared to be driven by the LERG sub-population (Ineson et al. 2013).
We have now repeated the analysis with a low redshift sample (z~0.1), and
found strong correlations between radio luminosity and environment richness and
between radio luminosity and central density for the LERGs but not for the
HERGs. These results are consistent with models in which the HERGs are fuelled
from accretion discs maintained from local reservoirs of gas, while LERGs are
fuelled more directly by gas ingested from the intra-cluster medium.
Comparing the samples, we found that although the maximum environment
richness of the HERG environments is similar in both samples, there are poorer
HERG environments in the z~0.1 sample than in the z~0.5 sample. We have
therefore tentative evidence of evolution of the HERG environments. We found no
differences between the LERG sub-samples for the two epochs, as would be
expected if radio and cluster luminosity are related.
Galactic Gamma-Ray Bursts (GRBs) are copious sources of gamma-rays that can pose a threat to complex life. Using recent determinations of their rate and the probability of GRBs causing massive extinction, we explore what type of universes are most likely to harbour advanced forms of life. For this purpose we use cosmological N-body simulations to determine at what time and for what value of the cosmological constant ($\Lambda$) the chances of life being unaffected by cosmic explosions are maximised. We find that $\Lambda-$dominated universes favour the survival of life against GRBs. Within a $\Lambda$CDM cosmology, the parameters that govern the likelihood of life survival to GRBs are dictated by the value of $\Lambda$ and the age of the Universe. We find that we seem to live in a favorable point in this parameter phase space which minimises the exposure to cosmic explosions, yet maximises the number of main sequence (hydrogen-burning) stars around which advanced life forms can exist.
We examine the fraction of massive asymptotic giant branch (AGB) stars remaining bound in their parent star clusters and the effect of irradiation of these stars by intracluster ultraviolet (UV) field. We employ a set of N-body models of dynamical evolution of star clusters rotating in a galactic potential at the solar galactocentric radius. The cluster models are combined with stellar evolution formulae, a library of stellar spectra, and simple models for SiO photodissociation in circumstellar environment (CSE). The initial stellar masses of clusters are varied from $50\rm M_\odot$ to $10^{5}\rm M_\odot$. Results derived for individual clusters are combined using a mass distribution function for young star clusters. We find that about 30% of massive AGB stars initially born in clusters become members of the field population, while the rest evolves in star clusters. They are irradiated by strong intracluster UV radiation resulting in the decrease of the photodissociation radius of SiO molecules, in many stars down to the dust formation zone. In absence of dust shielding, the UV photons penetrate in the CSE deeper than $10R_*$ in 64% and deeper than $2 R_*$ in 42% of all massive AGB stars. If this suppresses following dust formation, the current injection rate of silicate dust from AGB stars in the local Galaxy decreases from $2.2 \times 10^{-4}\rm M_\odot\,kpc^{-2}\,Gyr^{-1}$ to $1.8 \times 10^{-4}\rm M_\odot\,kpc^{-2}\,Gyr^{-1}$ at most. A lower revised value of 40% for the expected fraction of presolar silicate grains from massive AGB stars is still high to explain the non-detection of these grains in meteorites.
We perform a detailed X-ray study of the filaments surrounding the brightest cluster galaxies in a sample of nearby galaxy clusters using deep Chandra observations, namely the Perseus, Centaurus and Virgo clusters, and Abell 1795. We compare the X-ray properties and spectra of the filaments in all of these systems, and find that their Chandra X-ray spectra are all broadly consistent with an absorbed two temperature thermal model, with temperature components at 0.75 and 1.7 keV. We find that it is also possible to model the Chandra ACIS filament spectra with a charge exchange model provided a thermal component is also present, and the abundance of oxygen is suppressed relative to the abundance of Fe. In this model, charge exchange provides the dominant contribution to the spectrum in the 0.5-1.0 keV band. However, when we study the high spectral resolution RGS spectrum of the filamentary plume seen in X-rays in Centaurus, the opposite appears to be the case. The properties of the filaments in our sample of clusters are also compared to the X-ray tails of galaxies in the Coma cluster and Abell 3627. In the Perseus cluster, we search for signs of absorption by a prominent region of molecular gas in the filamentary structure around NGC 1275. We do find a decrement in the X-ray spectrum below 2 keV, indicative of absorption. However the spectral shape is inconsistent with this decrement being caused by simply adding an additional absorbing component. We find that the spectrum can be well fit (with physically sensible parameters) with a model that includes both absorption by molecular gas and X-ray emission from the filament, which partially counteracts the absorption.
In this work we present laboratory measurements on the reduction of the threshold friction velocity necessary for lifting dust if the dust bed is illuminated. Insolation of a porous soil establishes a temperature gradient. At low ambient pressure this gradient leads to thermal creep gas flow within the soil. This flow leads to a sub-surface overpressure which supports lift imposed by wind. The wind tunnel was run with Mojave Mars Simulant and air at 3, 6 and 9 mbar, to cover most of the pressure range at martian surface levels. Our first measurements imply that the insolation of the martian surface can reduce the entrainment threshold velocity between 4 % and 19 % for the conditions sampled with our experiments. An insolation activated soil might therefore provide additional support for aeolian particle transport at low wind speeds.
The radio signature of a shock travelling through the solar corona is known as a type II solar radio burst. In rare cases these bursts can exhibit a fine structure known as `herringbones', which are a direct indicator of particle acceleration occurring at the shock front. However, few studies have been performed on herringbones and the details of the underlying particle acceleration processes are unknown. Here, we use an image processing technique known as the Hough transform to statistically analyse the herringbone fine structure in a radio burst at $\sim$20-90 MHz observed from the Rosse Solar-Terrestrial Observatory on 2011 September 22. We identify 188 individual bursts which are signatures of bi-directional electron beams continuously accelerated to speeds of 0.16$_{-0.10}^{+0.11} c$. This occurs at a shock acceleration site initially at a constant altitude of $\sim$0.6 R$_{\odot}$ in the corona, followed by a shift to $\sim$0.5 R$_{\odot}$. The anti-sunward beams travel a distance of 170$_{-97}^{+174}$ Mm (and possibly further) away from the acceleration site, while those travelling toward the sun come to a stop sooner, reaching a smaller distance of 112$_{-76}^{+84}$ Mm. We show that the stopping distance for the sunward beams may depend on the total number density and the velocity of the beam. Our study concludes that a detailed statistical analysis of herringbone fine structure can provide information on the physical properties of the corona which lead to these relatively rare radio bursts.
Measurements of the magnetic field in the interplanetary medium, of the sunspots area, and of the heliospheric current sheet position, reveal a possible North-South asymmetry in the magnetic field of the Sun. We study the North-South asymmetry as inferred from measurements of the deflection of polar coronal hole jets when they propagate throughout the corona. Since the corona is an environment where the magnetic pressure is greater than the kinetic pressure, we can assume that magnetic field controls the dynamics of plasma. On average, jets during their propagation follow the magnetic field lines, highlighting its local direction. The average jet deflection is studied both in the plane perpendicular to the line of sight, and, for a reduced number of jets, in three dimensional space. The observed jet deflection is studied in terms of an axisymmetric magnetic field model comprising dipole We measured the position angles at 1 rs and at 2 rs of the 79 jets from the catalogue of Nistico et al 2009., based on the STEREO ultraviolet and white-light coronagraph observations during the solar minimum period March 2007-April 2008. We found that the propagation is not radial, in agreement with the deflection due to magnetic field lines. Moreover, the amount of the deflection is different between jets over the north and those from the south pole. Comparison of jet deflections and field line tracing shows that a ratio g2/g1 ~ -0.5 for the quadrupole and a ratio g3/g1 ~ 1.6-2.0 for the esapole can describe the field. The presence of a non-negligible quadrupole moment. We find that the magnetic deflection of jets is larger in the North than in the South of the order of 25-40%, with an asymmetry which is consistent with a southward deflection of the heliospheric current sheet of the order of 10 deg, consistent with that inferred from other, independent, datasets and instruments.
We test here the first stage of a route of modifications to be applied to the public GADGET2 code for dynamically identifying accretion centers during the collision process of two adjacent and identical gas cores. Each colliding core has a uniform density profile and rigid body rotation; its mass and size have been chosen to represent the observed core $L1544$; for the thermal and rotational energy ratios with respect to the potential energy, we assume the values $\alpha=0.3$ and $\beta=0.1$, respectively. These values favor the gravitational collapse of the core. We here study cases of both -head-on and off-center collisions, in which the pre-collision velocity increases the initial sound speed of the barotropic gas by up to several times. In a simulation the accretion centers are formed by the highest density particles, so we here report their location and properties in order to realize the collision effects on the collapsing and colliding cores. In one of the models we observe a roughly spherical distribution of accretion centers located at the front wave of the collision. In a forthcoming publication we will apply the full modified GADGET code to study the collision of turbulent cores.
Context: This paper is part of a larger project in which we study the chemical abundances of extra-galactic post-AGB stars with the ultimate goal of improving our knowledge of the poorly understood AGB third dredge-up mixing processes and s-process nucleosynthesis. Aims: In this paper, we study two carefully selected post-AGB stars in the LMC. The combination of favourable atmospheric parameters for detailed abundance studies and their known distances make these objects ideal probes of the internal AGB third dredge-up and s-process nucleosynthesis in that they provide observational constraints for theoretical AGB models. Methods: We use high-resolution optical UVES spectra to determine accurate stellar parameters and perform detailed elemental abundance studies. Additionally, we use available photometric data to construct SEDs for reddening and luminosity determinations. We then estimate initial masses from theoretical post-AGB tracks. Results: Both stars show extreme s-process enrichment associated with relatively low C/O ratios of about 1.3. We could only derive upper limits of the lead (Pb) abundance which indicate no strong Pb overabundances with respect to other s-elements. Comparison with theoretical post-AGB evolutionary tracks in the HR-diagram reveals that both stars have low initial masses between 1.0 and 1.5 Msun. Conclusion: This study adds to the results obtained so far on a very limited number of s-process enriched post-AGB stars in the Magellanic Clouds. We find an increasing discrepancy between observed and predicted Pb abundances towards lower metallicities for all studied Magellanic Cloud post-AGB stars found so far, as well as moderate C/O ratios. We find that all s-process rich post-AGB stars in the LMC and SMC studied so far, cluster in the same region of the HR-diagram and are associated with low-mass stars with a low metallicity on average.
We study the synchrotron component of the spectral energy distribution (SED) on the sample of 877 blazars using the ASDC SED Builder Tool with available broadband data from the literature. Our sample includes 423 flat-spectrum radio sources (FSRQs), 361 BL Lac objects and candidates, and 93 blazars of uncertain type. We have made an estimation of the synchrotron peak frequency ($\nu_{peak}^{s}$) for the 875 objects and further classified them as high, intermediate and low synchrotron peaked sources (HSPs/ISPs/LSPs). There are 42 HSPs with $\nu_{peak}^{s} > 10^{16.5}$ Hz, 222 ISPs with $10^{14.5} < \nu_{peak}^{s} < 10^{16.5}$ Hz, and 611 LSPs with $\nu_{peak}^{s} < 10^{14.5}$ Hz in our sample. We have calculated an average value of $\nu_{peak}^{s}$ to be $10^{13.4 \pm 1.0}$ Hz for FSRQs and $10^{14.6 \pm 1.4}$ Hz for BL Lacs. We found out that $\nu_{peak}^{s}$ and the flux density at 4.8 GHz have a different distribution (as indicated by Kolmogorov-Smirnov test at significance level 0.05) for the FSRQ and BL Lac blazars, and for the RBL and XBL types of BL Lacs. Distribution of $\nu_{peak}^{s}$ values is broader for BL Lacs, than for FSRQs. There are no ultra-high energy peaked objects (with $\nu_{peak}^{s} > 10^{19}$ Hz) in our BL Lac sample according to our estimations. The significant part of FSRQs (41%) and small part of BL Lacs (9%) in our sample could be considered as candidates to the very-low synchrotron peaked blazars (with $\nu_{peak}^{s} < 10^{13}$ Hz). Our foundations confirm results of the previous studies made on samples with significantly smaller number of objects.
We propose a new way of looking at the Baryon Acoustic Oscillations in the Large Scale Structure clustering correlation function. We identify a scale s_LP that has two fundamental features: its position is insensitive to non-linear gravity, redshift space distortions, and scale-dependent bias at the 0.5% level; it is geometrical, i.e. independent of the power spectrum of the primordial density fluctuation parameters. These two properties together make s_LP, called the "linear point", an excellent cosmological standard ruler. The linear point is also appealing because it is easily identified irrespectively of how non-linearities distort the correlation function. Finally, the correlation function amplitude at s_LP is similarly insensitive to non-linear corrections to within a few percent. Hence, exploiting the particular Baryon features in the correlation function, we propose three new estimators for growth measurements. A preliminary analysis of s_LP in current data is encouraging.
High-energy neutrinos could be produced in the interaction of charged cosmic rays with matter or radiation surrounding astrophysical sources. Even with the recent detection of extraterrestrial high-energy neutrinos by the IceCube experiment, no astrophysical neutrino source has yet been discovered. Transient sources, such as gamma-ray bursts, core-collapse supernovae, or active galactic nuclei are promising candidates. Multi-messenger programs offer a unique opportunity to detect these transient sources. By combining the information provided by the ANTARES neutrino telescope with information coming from other observatories, the probability of detecting a source is enhanced, allowing the possibility of identifying a neutrino progenitor from a single detected event. A method based on optical and X-ray follow-ups of high-energy neutrino alerts has been developed within the ANTARES collaboration. This program, denoted as TAToO, triggers a network of robotic optical telescopes (TAROT and ROTSE) and the Swift-XRT with a delay of only a few seconds after a neutrino detection, and is therefore well-suited to search for fast transient sources. To identify an optical or X-ray counterpart to a neutrino signal, the images provided by the follow-up observations are analysed with dedicated pipelines. A total of 42 alerts with optical and 7 alerts with X-ray images taken with a maximum delay of 24 hours after the neutrino trigger have been analysed. No optical or X-ray counterparts associated to the neutrino triggers have been found, and upper limits on transient source magnitudes have been derived. The probability to reject the gamma-ray burst origin hypothesis has been computed for each alert.
As one of the best-characterized stellar-mass black holes, with good measurements of its mass, distance and inclination, V404 Cyg is the ideal candidate to study Eddington-limited accretion episodes. After a long quiescent period, V404 Cyg underwent a new outburst in June 2015. We obtained two Chandra HETG exposures of 20 ksec and 25 ksec. Many strong emission lines are observed; the ratio of Si He-like triplet lines gives an estimate for the formation region distance of $4\times10^{11}$ cm, while the higher ionization Fe XXV He-like triplet gives an estimate of $7\times10^9$ cm. A narrow Fe K$\alpha$ line is detected with an equivalent width greater than 1 keV in many epochs, signaling that we do not directly observe the central engine. Obscuration of the central engine and strong narrow emission lines signal that the outer disk may be illuminated, and its structure may help to drive the strong variability observed in V404 Cyg. In the highest flux phases, strong P-Cygni profiles consistent with a strong disk wind are observed. The kinetic power of this wind may be extremely high.
We show that a standard model gauge singlet fermion field, with mass of order keV or larger, and involved in the inverse seesaw mechanism of light neutrino mass generation, can be a good warm dark matter candidate. Our framework is based on B-L extension of the Standard Model. The construction ensures the absence of any mixing between active neutrinos and the aforementioned dark matter field. This circumvents the usual constraints on the mass of warm dark matter imposed by X-ray observations. We show that over-abundance of thermally produced warm dark matter (which nevertheless do not reach chemical equilibrium) can be reduced to an acceptable range in the presence of a moduli field decaying into radiation --- though only when the reheat temperature is low enough. Our warm dark matter candidate can also be produced directly from the decay of the moduli field during reheating. In this case, obtaining the right amount of relic abundance, while keeping the reheat temperature high enough as to be consistent with Big Bang nucleosynthesis bounds, places constraints on the branching ratio for the decay of the moduli field into dark matter.
In the context of neutron stars mergers, we study the gravitational wave spectrum of the merger remnant using numerical relativity simulations. Postmerger spectra are characterized by a main peak frequency $f_2$ related to the particular structure and dynamics of the remnant hot hypermassive neutron star. We show that $f_2$ is correlated with the tidal coupling constant $\kappa^T_2$ that characterizes the binary tidal interactions during the late-inspiral--merger. The relation $f_2(\kappa^T_2)$ depends very weakly on the binary total mass, mass-ratio, equation of state, and thermal effects. This observation opens up the possibility of developing a model of the gravitational spectrum of every merger unifying the late-inspiral and postmerger descriptions.
Real scalar fields are known to fragment into spatially localized and long-lived solitons called oscillons or $I$-balls. We prove the adiabatic invariance of the oscillons/$I$-balls for a potential that allows periodic motion even in the presence of non-negligible spatial gradient energy. We show that such potential is uniquely determined to be the quadratic one with a logarithmic correction, for which the oscillons/$I$-balls are absolutely stable. For slightly different forms of the scalar potential dominated by the quadratic one, the oscillons/$I$-balls are only quasi-stable, because the adiabatic charge is only approximately conserved. We check the conservation of the adiabatic charge of the $I$-balls in numerical simulation by slowly varying the coefficient of logarithmic corrections. This unambiguously shows that the longevity of oscillons/$I$-balls is due to the adiabatic invariance.
The performance of a radiatively cooled instrument is investigated in the context of optomechanical quantum experiments, where the environment of a macroscopic particle in a quantum-superposition has to be cooled to less than 20\,K in deep space. A heat-transfer analysis between the components of the instrument as well as a transfer-function analysis on thermal oscillations induced by the spacecraft interior and by dissipative sources is performed. The thermal behaviour of the instrument in an orbit around a Lagrangian point and in a highly elliptical Earth orbit is discussed. Finally, we investigate further possible design improvements aiming at lower temperatures of the environment of the macroscopic particle. These include a mirror-based design of the imaging system on the optical bench and the extension of the heat shields.
It is shown how quantum fluctuations of the radiation during the contraction era of a CBE (Comes Back Empty) cyclic cosmology can provide density fluctuations which re-enter the horizon during the subsequent expansion era and at lowest order are scale invariant, in a Harrison-Zel'dovich-Peebles sense, as necessary to be consistent with observations of large scale structure.
We discuss the possibility of devising cosmological observables which violate Bell's inequalities. Such observables could be used to argue that cosmic scale features were produced by quantum mechanical effects in the very early universe. As a proof of principle, we propose a somewhat elaborate inflationary model where a Bell inequality violating observable can be constructed.
Different mechanisms operate in various regions of the MSSM parameter space to bring the relic density of the lightest neutralino, neutralino_1, assumed here to be the LSP and thus the Dark Matter (DM) particle, into the range allowed by astrophysics and cosmology. These mechanisms include coannihilation with some nearly-degenerate next-to-lightest supersymmetric particle (NLSP) such as the lighter stau (stau_1), stop (stop_1) or chargino (chargino_1), resonant annihilation via direct-channel heavy Higgs bosons H/A, the light Higgs boson h or the Z boson, and enhanced annihilation via a larger Higgsino component of the LSP in the focus-point region. These mechanisms typically select lower-dimensional subspaces in MSSM scenarios such as the CMSSM, NUHM1, NUHM2 and pMSSM10. We analyze how future LHC and direct DM searches can complement each other in the exploration of the different DM mechanisms within these scenarios. We find that the stau_1 coannihilation regions of the CMSSM, NUHM1, NUHM2 can largely be explored at the LHC via searches for missing E_T events and long-lived charged particles, whereas their H/A funnel, focus-point and chargino_1 coannihilation regions can largely be explored by the LZ and Darwin DM direct detection experiments. We find that the dominant DM mechanism in our pMSSM10 analysis is chargino_1 coannihilation: {parts of its parameter space can be explored by the LHC, and a larger portion by future direct DM searches.
High-amplitude, chaotic/turbulent electromagnetic fluctuations are ubiquitous in high-energy-density laboratory and astrophysical plasmas, where they can be excited by various kinetic-streaming and/or anisotropy-driven instabilities, such as the Weibel instability. These fields typically exist on "sub-Larmor scales" -- scales smaller than the electron Larmor radius. Electrons moving through such magnetic fields undergo small-angle stochastic deflections of their pitch-angles, thus establishing diffusive transport on long time-scales. We show that this behavior, under certain conditions, is equivalent to Coulomb collisions in collisional plasmas. The magnetic pitch-angle diffusion coefficient, which acts as an effective "collision" frequency, may be substantial in these, otherwise, collisionless environments. We show that this effect, colloquially referred to as the plasma "quasicollisionality", may radically alter the expected radiative transport properties of candidate plasmas. We argue that the modified magneto-optic effects in these plasmas provide an attractive, novel radiative diagnostic tool for the exploration and characterization of small-scale magnetic turbulence, as well as affect inertial confinement fusion and other laser-plasma experiments.
Links to: arXiv, form interface, find, astro-ph, recent, 1508, contact, help (Access key information)
Planets form in the discs of gas and dust that surround young stars. It is not known whether gas giant planets on wide orbits form the same way as Jupiter or by fragmentation of gravitationally unstable discs. Here we show that a giant planet, which has formed in the outer regions of a protostellar disc, initially migrates fast towards the central star (migration timescale ~10,000 yr) while accreting gas from the disc. However, in contrast with previous studies, we find that the planet eventually opens up a gap in the disc and the migration is essentially halted. At the same time, accretion-powered radiative feedback from the planet, significantly limits its mass growth, keeping it within the planetary mass regime (i.e. below the deuterium burning limit) at least for the initial stages of disc evolution. Giant planets may therefore be able to survive on wide orbits despite their initial fast inward migration, shaping the environment in which terrestrial planets that may harbour life form.
We present the analysis of a Chandra High-Energy Transmission Grating (HETG) observation of the local Seyfert galaxy NGC 1365. The source, well known for its dramatic X-ray spectral variability, was caught in a reflection-dominated, Compton-thick state. The high spatial resolution afforded by Chandra allowed us to isolate the soft X-ray emission from the active nucleus, neglecting most of the contribution from the kpc-scale starburst ring. The HETG spectra thus revealed a wealth of He- and H-like lines from photoionized gas, whereas in larger aperture observations these are almost exclusively produced through collisional ionization in the circumnuclear environment. Once the residual thermal component is accounted for, the emission-line properties of the photoionized region close to the hard X-ray continuum source indicate that NGC 1365 has some similarities to the local population of obscured active galaxies. In spite of the limited overall data quality, several soft X-ray lines seem to have fairly broad profiles (~800-1300 km/s full-width at half maximum), and a range of outflow velocities (up to ~1600 km/s, but possibly reaching a few thousands km/s) appears to be involved. At higher energies, the K$\alpha$ fluorescence line from neutral iron is resolved with > 99 per cent confidence, and its width of ~3000 km/s points to an origin from the same broad-line region clouds responsible for eclipsing the X-ray source and likely shielding the narrow-line region.
We combine constraints on galaxy formation histories with planet formation models, yielding the Earth-like and giant planet formation histories of the Milky Way and the Universe as a whole. In the Hubble Volume (10^13 Mpc^3), we expect there to be ~10^20 Earth-like and ~10^20 giant planets; our own galaxy is expected to host ~10^9 and ~10^10 Earth-like and giant planets, respectively. Proposed metallicity thresholds for planet formation do not significantly affect these numbers. However, the metallicity dependence for giant planets results in later typical formation times and larger host galaxies than for Earth-like planets. The Solar System formed at the median age for existing giant planets in the Milky Way, and consistent with past estimates, formed after 80% of Earth-like planets. However, if existing gas within virialised dark matter haloes continues to collapse and form stars and planets, the Universe will form over 10 times more planets than currently exist. We show that this would imply at least a 92% chance that we are not the only civilisation the Universe will ever have, independent of arguments involving the Drake Equation.
We present spectroscopic confirmation of two new lensed quasars via data obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been selected from the Dark Energy Survey (DES) and WISE based on their multi-band photometry and extended morphology in DES images. Images of DES J0115-5244 show two blue point sources at either side of a red galaxy. Our long-slit data confirm that both point sources are images of the same quasar at $z_{s}=1.64.$ The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two blue components are visible in the DES and SDSS images. The SDSS fiber spectrum shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground lens galaxy. The long-slit Magellan spectra show that the blue components are resolved images of the same quasar. The Einstein Radius is $0.68$" corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other candidates were observed and rejected, two being low-redshift pairs of starburst galaxies, and one being a quasar behind a blue star. These first confirmation results provide an important empirical validation of the data-mining and model-based selection that is being applied to the entire DES dataset.
We present a model for the evolution of the galaxy ultraviolet (UV) luminosity function (LF) across cosmic time where star formation is linked to the assembly of dark matter halos under the assumption of a halo mass dependent, but redshift independent, star formation efficiency. This model improves on previous work by introducing a new self-consistent treatment of the halo star formation history, which allows us to make predictions at redshift $z>10$ (lookback time $\lesssim500$ Myr), when growth is rapid. With a calibration at a single redshift to set the stellar to halo mass ratio, and no further degrees of freedom, our model captures the evolution of the UV LF over all the available observations ($0\lesssim z\lesssim10$). The significant drop in the luminosity density of currently detectable galaxies beyond $z\sim8$ is explained by a shift of star formation toward less massive, fainter galaxies. Assuming that star formation proceeds down to atomic cooling halos, we derive a reionization optical depth $\tau = 0.056^{+0.007}_{-0.010}$ fully consistent with the latest Planck measurement, and implying that the universe is fully reionized at $z=7.84^{+0.65}_{-0.98}$. In addition, our model naturally produces smoothly rising star formation histories for galaxies with $L\lesssim L_*$ in agreement with observations and detailed hydrodynamical simulations. Before the epoch of reionization at $z>10$ we predict the LF to remain well-described by a Schechter function, but with an increasingly steep faint-end slope ($\alpha\sim-3.5$ at $z\sim16$). Finally, we construct detailed forecasts for surveys with JWST and WFIRST, including the boost from gravitational lensing magnification bias in blank fields, and predict that galaxies out to $z\sim14$ will be observed. However, galaxies at $z>15$ will likely be accessible to JWST and WFIRST only through the assistance of strong lensing magnification.
We summarize three recent efforts to constrain the first few moments of cosmic creation before and during the epoch of inflation. We consider two means to explain a slight dip in the power spectrum of the cosmic microwave background for multipoles in the range of $\ell= 10-30$ from both the {\it Planck} and {\it WMAP} data. We show that such a dip could be the result of resonant creation of a massive particle that couples to the inflaton field. For best-fit models, the epoch of resonant particle creation reenters the horizon at wave numbers of $k_* \sim 0.00011 \pm 0.0004 $ ($h$ Mpc$^{-1}$). The amplitude and location of these features correspond to the creation of a number of degenerate fermion species of mass $\sim 15/\lambda^{3/2} $ $m_{pl}$ during inflation where $\lambda$ is the coupling constant between the inflaton field and the created fermion species. Alternatively, one can explain the existence of such a dip as due to a jump in the inflation generating potential. We show that such a jump can also resolve the excessively large dark flow predicted from the M-theory landscape. Finally, we summarize our efforts to quantify constraints on the cosmic dark flow from a new analysis of the Type Ia supernova distance-redshift relation.
In this work we introduce a new way of binning sunspot group data with the
purpose of better understanding the impact of the solar cycle on sunspot
properties and how this defined the characteristics of the extended minimum of
cycle 23. Our approach assumes that the statistical properties of sunspots are
completely determined by the strength of the underlying large-scale field and
have no additional time dependencies. We use the amplitude of the cycle at any
given moment (something we refer to as activity level) as a proxy for the
strength of this deep-seated magnetic field.
We find that the sunspot size distribution is composed of two populations:
one population of groups and active regions and a second population of pores
and ephemeral regions. When fits are performed at periods of different activity
level, only the statistical properties of the former population, the active
regions, is found to vary.
Finally, we study the relative contribution of each component (small-scale
versus large-scale) to solar magnetism. We find that when hemispheres are
treated separately, almost every one of the past 12 solar minima reaches a
point where the main contribution to magnetism comes from the small-scale
component. However, due to asymmetries in cycle phase, this state is very
rarely reached by both hemispheres at the same time. From this we infer that
even though each hemisphere did reach the magnetic baseline, from a
heliospheric point of view the minimum of cycle 23 was not as deep as it could
have been.
We present the barium surface abundance of 12 blue stragglers (BSs) and 18 main-sequence (MS) stars in the intermediate-age open cluster NGC 6819 (2.5 Gyr) based on spectra obtained from the Hydra Multi-object Spectrograph on the WIYN 3.5 m telescope. For the MS stars we find [Fe/H] = $+$0.05 $\pm$ 0.04 and [Ba/Fe] = $-$0.01 $\pm$ 0.10. The majority of the BS stars are consistent with these values. We identify five BSs with significant barium enhancement. These stars most likely formed through mass transfer from an asymptotic giant branch star that polluted the surface of the BS with the nucleosynthesis products generated during thermal pulsations. This conclusion aligns with the results from the substantial work done on the BSs in old open cluster NGC 188 that identifies mass transfer as the dominant mechanism for BS formation in that open cluster. However, four of the BSs with enhanced barium show no radial-velocity evidence for a companion. The one star that is in a binary is a double-lined system, meaning the companion is not a white dwarf and not the remnant of a prior AGB star. In this paper we attempt to develop a consistent scenario to explain the origin of these five BSs.
We show that the Ratra model, where the inflaton is kinetically coupled to the photon, is a successful scenario of cosmic inflationary magnetogenesis, which is free from strong-coupling and backreaction problems.
Characterizing the physical properties of exoplanets, and understanding their formation and orbital evolution requires precise and accurate knowledge of the physical properties of their host stars. Accurately measuring stellar mass is particularly important because the masses of host stars likely influence planet occurrence and the architectures of planetary systems observed today. Single main-sequence stars typically have masses estimated from evolutionary tracks, which generally provide accurate results due to their extensive empirical calibration. However, the validity of this method for subgiants and giants has been called into question, with suggestions that the evolutionary models could contain systematic errors that would cause mass estimates of these evolved stars to be overestimated. We investigate these concerns using a sample of 59 benchmark evolved stars with model-independent masses (from binary systems or asteroseismology) obtained from the extant literature. We find very good agreement between these benchmark masses and the ones estimated using evolutionary tracks. The average fractional difference in the mass interval $\sim$0.7 - 4.5 $M_{\odot}$ is consistent with zero (-1.30 $\pm$ 2.42%), with no significant trends in the residuals relative to the input parameters. Six stars in our benchmark sample can be classified as Retired A Stars similar to the targets of various Doppler surveys. Our analysis of these small subset of stars reveals a systematic offset of -8.97 $\pm$ 8.93%, suggesting a possible underestimate rather than an overestimate of stellar masses in this region of the H-R diagram. Taken together, our results indicate that determination of masses of evolved stars using grids of evolutionary tracks is not significantly affected by systematic errors, and is thus valid for estimating the masses of isolated stars beyond the main sequence.
We report the discovery of an occultation event in the low-luminosity narrow-line Seyfert 1 galaxy WPVS 007 in 2015 February and March. In concert with longer timescale variability, these observations place strong constraints on the nature and location of the absorbing material. Swift monitoring has revealed a secular decrease since ~2010 accompanied by flattening of the optical and UV photometry that suggests variable reddening. Analysis of four Hubble Space Telescope COS observations since 2010, including a Director's Discretionary time observation during the occultation, shows that the broad-absorption-line velocity offset and the CIV emission-line width both decrease as the reddening increases. The occultation dynamical timescale, the BAL variability dynamical timescale, and the density of the BAL gas show that both the reddening material and the broad-absorption-line gas are consistent with an origin in the torus. These observations can be explained by a scenario in which the torus is clumpy with variable scale height, and the BAL gas is blown from the torus material like spray from the crest of a wave. As the obscuring material passes into our line of sight, we alternately see high-velocity broad absorption lines and a clear view to the central engine, or low-velocity broad absorption lines and strong reddening. WPVS 007 has a small black hole mass, and correspondingly short timescales, and so we may be observing behavior that is common in BALQSOs, but is not typically observable.
SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS u band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discovered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along ~25 degrees, consists of 20-30 BHBs with a total significance of around 6sigma over the background. Assuming we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs suggests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system.
Variability classes in the enigmatic black hole candidate GRS 1915+105 are known to be correlated with the variation of the Comptonizing Efficiency (CE) which is defined to be the ratio between the number of power-law (hard) photons and seed (soft) photons injected into the Compton cloud. Similarities of light curves of several variability classes of GRS 1915+105 and IGR 17091-3624, some of which are already reported in the literature, motivated us to compute CE for IGR 17091-3624 as well. We find that they are similar to what were reported earlier for GRS 1915+105, even though masses of these objects could be different. The reason is that the both the sizes of the sources of the seed photons and of the Comptonizing corona scale in the same way as the mass of the black hole. This indicates that characterization of variability classes based on CE is likely to be black hole mass independent, in general.
We present the SALT spectroscopy of a globular cluster in the center of the nearby isolated dSph galaxy KKs3 situated at a distance of 2.12 Mpc. Its heliocentric radial velocity is 316+-7 km/s that corresponds to V_{LG} = 112 km/s in the Local Group (LG) reference frame. We use its distance and velocity along with the data on other 35 field galaxies in the proximity of the LG to trace the local Hubble flow. Some basic properties of the local field galaxies: their morphology, absolute magnitudes, average surface brightnesses, specific star formation rates, and hydrogen mass-to-stellar mass ratios are briefly discussed. Surprisingly, the sample of the neighboring isolated galaxies displays no signs of compression under the influence of the expanding Local Void.
Using a stellar mass limited sample of $\sim 46,600$ galaxies ($M_* > 10^{9.1}\,M_{\odot}$) at $0.5 < z < 2$, we show that the tellar mass, rather than the environment, is the main parameter controlling quenching of star formation in galaxies with $M_* > 10^{10}\,M_{\odot}$ out to $z=2$. On the other hand, the environmental quenching becomes efficient at $z < 1$ regardless of galaxy mass, and it serves as a main star formation quenching mechanism for lower mass galaxies. Our result is based on deep optical and near-infrared imaging data over 2800 arcmin$^2$, enabling us to negate cosmic variance and identify 46 galaxy cluster candidates with $M \sim 10^{14}\,M_{\odot}$. From $M_* \sim 10^{9.5}$ to $10^{10.5}\,M_{\odot}$, the fraction of quiescent galaxies increases by a factor of $\sim 10$ over the entire redshift range, but the difference between cluster and field environment is negligible. Rapid evolution in the quiescent fraction is seen from $z=2$ to $z=1.3$ for massive galaxies suggesting a build-up of massive quiescent galaxies at $z > 1.3$. For galaxies with $M_* < 10^{10}\,M_{\odot}$ at $z < 1.0$, the quiescent fraction is found to be as much as a factor of 2 larger in clusters than in field, showing the importance of environmental quenching in low mass galaxies at low redshift. Most high mass galaxies are already quenched at $z > 1$, therefore environmental quenching does not play a significant role for them, although the environmental quenching efficiency is nearly identical between high and low mass galaxies.
Open clusters are known as excellent tracers of the structure and chemical evolution of the Galactic disk, however, the accuracy and reliability of open cluster parameters is poorly known. In recent years, several studies aimed to present homogeneous open cluster parameter compilations, which are based on some different approaches and photometric data. These catalogues are excellent sources to facilitate testing of the actual accuracy of open cluster parameters. We compare seven cluster parameter compilations statistically and with an external sample, which comprises the mean results of individual studies. Furthermore, we selected the objects IC 4651, NGC 2158, NGC 2383, NGC 2489, NGC 2627, NGC 6603, and Trumpler 14, with the main aim to highlight differences in the fitting solutions. We derived correction terms for each cluster parameter, using the external calibration sample. Most results by the compilations are reasonable scaled, but there are trends or constant offsets of different degree. We also identified one data set, which appears too erroneous to allow adjustments. After the correction, the mean intrinsic errors amount to about 0.2 dex for the age, 0.08 mag for the reddening, and 0.35 mag for the distance modulus. However, there is no study that characterises the cluster morphologies of all test cases in a correct and consistent manner. Furthermore, we found that the largest compilations probably include at least 20 percent of problematic objects, for which the parameters differ significantly. These could be among others doubtful or unlikely open clusters that do not facilitate an unambiguous fitting solution.
We present the results of archival coincidence analyses using public neutrino data from the 40-string configuration of IceCube (IC40) and contemporaneous public gamma-ray data from Fermi LAT. Our analyses have the potential to discover statistically significant coincidences between high-energy neutrino and gamma-ray signals, and hence, possible jointly-emitting neutrino/gamma-ray transients. This work is an example of more general multimessenger studies that the Astrophysical Multimessenger Observatory Network (AMON) aims to perform. AMON is currently under development and will link multiple running and future high-energy neutrino, cosmic ray and follow-up observatories as well as gravitational wave facilities. This single network will enable near real-time coincidence searches for multimessenger astrophysical transients and their electromagnetic counterparts. We will present the component high-energy neutrino and gamma-ray datasets, the statistical approaches that we used, and the results of analyses of the IC40+LAT datasets.
We report polarimetry results of a merging cluster of galaxies Abell 2256 with Karl G. Jansky Very Large Array (JVLA). We performed new observations with JVLA at S-band (2051-3947 MHz) and X-band (8051-9947 MHz) in the C array configuration, and detected significant polarized emissions from the radio relic, Source A, and Source B in this cluster. We calculated the total magnetic field strengths toward the radio relic using revised equipartition formula, which is 1.8-5.0 microG. With dispersions of Faraday rotation measure, magnetic-field strengths toward Sources A and B are estimated to be 0.63-1.26 microG and 0.11-0.21 microG, respectively. An extremely high degree of linear polarization, as high as ~ 35 %, about a half of the maximum polarization, was detected toward the radio relic, which indicates highly ordered magnetic lines of force over the beam sizes (~ 52 kpc).The fractional polarization of the radio relic decreases from ~ 35 % to ~ 20 % around 3 GHz as the frequency decreases and is nearly constant between 1.37 and 3 GHz. Both analyses with depolarization models and Faraday tomography suggest multiple depolarization components toward the radio relic and imply the existence of turbulent magnetic fields.
We present our analysis of the X-ray variability of two ultraluminous X-ray sources (ULXs) based on multiple XMM--Newton observations. We show the linear rms-flux relation is present in eight observations of NGC5408 X-1 and also in three observations of NGC6946 X-1, but data from other ULXs are generally not sufficient to constrain any rms-flux relation. The presence of this relation was previously reported in only two observations of NGC 5408X-1; our results show this is a persistent property of the variability of NGC5408 X-1 and extends to at least one other variable ULX. We speculate this is a ubiquitous property of ULX variability, as it is for X-ray variability in other luminous accreting sources. We also recover the time delay between hard and soft bands in NGC5408 X-1, with the soft band (<1 keV) delayed with respect to the hard band (>1 keV) by up to ~10 s (~0.2 rad) at frequencies above ~few mHz. For the first time, we extend the lag analysis to lower frequencies and find some evidence for a reversal of the lag, a hard lag of ~1 ks at frequencies of ~0.1 mHz. Our energy-resolved analysis shows the time delays are energy dependent. We argue that the lag is unlikely to be a result of reflection from an accretion disc (`reverberation') based on the lack of reflection features in the spectra, and the large size of the reflector inferred from the magnitude of the lag. We also argue that associating the soft lag with a quasi-periodic oscillation (QPO) in these ULXs - and drawing an analogy between soft lags in ULXs and soft lags seen in some low-frequency QPOs of Galactic X-ray binaries - is premature.
Polytropic transonic solutions of spherically symmetric and steady galactic winds in the gravitational potential of a dark matter halo (DMH) with a supermassive black hole (SMBH) are studied. The solutions are classified in terms of their topological features, and the gravitational potential of the SMBH adds a new branch to the transonic solutions generated by the gravity of the DMH. The topological types of the transonic solutions depend on the mass distribution, the amount of supplied energy, the polytropic index $\gamma$, and the slope $\alpha$ of the DMH mass distribution. When $\alpha$ becomes larger than a critical value $\alpha_\mathrm{c}$, the transonic solution types change dramatically. Further, our model predicts that it is possible for a slowly accelerating outflow to exist, even in quiescent galaxies with small $\gamma$. This slowly accelerating outflow differs from those considered in many of the previous studies focusing on supersonic outflows in active star-forming galaxies. In addition, our model indicates that outflows in active star-forming galaxies have only one transonic point in the inner region ($\sim$ 0.01 kpc). The locus of this transonic point does not strongly depend on $\gamma$. We apply the polytropic model incorporating mass flux supplied by stellar components to the Sombrero Galaxy, and conclude that it can reproduce the observed gas density and the temperature distribution well. This result differs significantly from the isothermal model, which requires an unrealistically large mass flux (Igarashi et al. 2014). Thus, we conclude that the polytropic model is more realistic than the isothermal model, and that the Sombrero Galaxy can have a slowly accelerating outflow.
HST images of proplyds in the Orion Nebula, as well as submillimeter/radio measurements, show that the dominant O7 star Theta1 Ori C photoevaporates nearby disks around pre-main sequence stars. Theory predicts that massive stars photoevaporate disks within distances of order 0.1 pc. These findings suggest that young, OB-dominated massive H II regions are inhospitable to the survival of protoplanetary disks, and subsequently to the formation and evolution of planets. In the current work, we test this hypothesis using large samples of pre-main sequence stars in 20 massive star-forming regions selected with X-ray and infrared photometry in the MYStIX survey. Complete disk destruction would lead to a deficit of cluster members with excess in JHKs and Spitzer/IRAC bands in the vicinity of O stars. In four MYStIX regions containing O stars and a sufficient surface density of disk-bearing sources to reliably test for spatial avoidance, we find no evidence for the depletion of inner disks around pre-main sequence stars in the vicinity of O-type stars, even very luminous O2-O5 stars. These results suggest that massive star-forming regions are not very hostile to the survival of protoplanetary disks and, presumably, to the formation of planets.
We report on the results of the extensive multi-wavelength campaign from optical to GeV gamma-rays of the 2014 periastron passage of PSR B1259-63, which is a unique high-mass gamma-ray emitting binary system with a young pulsar companion. Observations demonstrate the stable nature of the post-periastron GeV flare and prove the coincidence of the flare with the start of rapid decay of the H$\alpha$ equivalent width, usually interpreted as a disruption of the Be stellar disk. Intensive X-ray observations reveal changes in the X-ray spectral behaviour happening at the moment of the GeV flare. We demonstrate that these changes can be naturally explained as a result of synchrotron cooling of monoenergetic relativistic electrons injected into the system during the GeV flare.
The evolution of a pulsar wind nebula (PWN) inside a supernova remnant (SNR) is sensitive to properties of the central neutron star, pulsar wind, progenitor supernova, and interstellar medium. These properties are both difficult to measure directly and critical for understanding the formation of neutron stars and their interaction with the surrounding medium. In this paper, we determine these properties for PWN G54.1+0.3 by fitting its observed properties with a model for the dynamical and radiative evolution of a PWN inside an SNR. Our modeling suggests that the progenitor of G54.1+0.3 was an isolated ~15-20 Solar Mass star which exploded inside a massive star cluster, creating a neutron star initially spinning with period ~30-80ms. We also find that >99.9% of the pulsar's rotational energy is injected into the PWN as relativistic electrons and positrons whose energy spectrum is well characterized by a broken power-law. Lastly, we propose future observations which can both test the validity of this model and better determine the properties of this source -- in particular, its distance and the initial spin period of the central pulsar.
We report the discovery of a substantial stellar stream in the periphery of the Large Magellanic Cloud (LMC), found using public imaging from the first year of the Dark Energy Survey. The stream appears to emanate from the edge of the outer LMC disk at a radius $\approx 13.5$ degrees due north of its centre, and stretches more than $10$ kpc towards the east. It is roughly $1.5$ kpc wide and has an integrated $V$-band luminosity of at least $M_V = -7.4$. The stellar populations in the stream are indistinguishable from those in the outer LMC disk. We attempt to quantify the geometry of the outer disk using simple planar models, and find that only a disk with mild intrinsic ellipticity can simultaneously explain the observed stellar density on the sky and the azimuthal line-of-sight distance profile. We also see possible non-planar behaviour in the outer disk that may reflect a warp and/or flare. Based on all these observations, we conclude that the stream is most likely comprised of material that has been stripped from the outskirts of the LMC disk. We conduct a simple $N$-body simulation to show that this is plausibly due to the tidal force of the Milky Way, but we cannot rule out a recent close interaction between the LMC and the SMC as the source of the stripping. Finally, we observe tentative evidence for extremely diffuse LMC populations beyond the outer edge of the stream, at radii of up to $\sim 18.5$ kpc in the disk plane. These stars could serve as useful tracers of the total mass and orbital history of the LMC.
Reproducing the large Earth/Mars mass ratio requires a strong mass depletion in solids within the protoplanetary disk between 1 and 3 AU. The Grand Tack model invokes a specific migration history of the giant planets to remove most of the mass initially beyond 1 AU and to dynamically excite the asteroid belt. However, one could also invoke a steep density gradient created by inward drift and pile-up of small particles induced by gas-drag, as has been proposed to explain the formation of close-in super Earths. Here we show that the asteroid belt's orbital excitation provides a crucial constraint against this scenario for the Solar System. We performed a series of simulations of terrestrial planet formation and asteroid belt evolution starting from disks of planetesimals and planetary embryos with various radial density gradients and including Jupiter and Saturn on nearly circular and coplanar orbits. Disks with shallow density gradients reproduce the dynamical excitation of the asteroid belt by gravitational self-stirring but form Mars analogs significantly more massive than the real planet. In contrast, a disk with a surface density gradient proportional to $r^{-5.5}$ reproduces the Earth/Mars mass ratio but leaves the asteroid belt in a dynamical state that is far colder than the real belt. We conclude that no disk profile can simultaneously explain the structure of the terrestrial planets and asteroid belt. The asteroid belt must have been depleted and dynamically excited by a different mechanism such as, for instance, in the Grand Tack scenario.
It is well known that annihilations in the homogeneous fluid of dark matter (DM) can leave substantial imprints in the cosmic microwave background (CMB) anisotropy power spectrum. However, the relevance of DM annihilations in halos is still subject to debate, with previous works reaching different conclusions on this point. Furthermore, models of DM annihilations in halos have been invoked to solve the tension between WMAP measurement of the reionization optical depth and astrophysical Gunn-Peterson bound, requiring a significantly smaller value of the optical depth to reionization. This tension, although smaller, still exists in the new Planck data. In this work, we revisit these problems and aim at clarifying the situation, thanks to the most accurate treatment of DM annihilations in halos to this day. We find that the ionization fraction does exhibit a very particular (and potentially constraining) pattern, but the currently measurable reionization optical depth is left almost unchanged: For plausible halo models the modification of the signal with respect to the one coming from annihilation in the smooth background is tiny, below cosmic variance within currently allowed parameter space. We thus conclude that the impact of the virialised DM structures cannot be uncovered by CMB power spectra measurements, unless very peculiar models are invoked for the redshift evolution of the DM annihilation signal (e.g. via unconventional velocity dependence of the annihilation cross section). On the other hand, a precise measurement of the ionization fraction or of the temperature history of the universe (notably via the 21 cm signal) seems to be the most promising way for using halo formation as a tool in DM searches, improving over the current sensitivity of cosmological probes.
The term `solar tornadoes' has been used to describe apparently rotating magnetic structures above the solar limb, as seen in high resolution images and movies from the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO). These often form part of the larger magnetic structure of a prominence, however the links between them remain unclear. Here we present plasma diagnostics on a tornado-like structure and its surroundings, seen above the limb by the Extreme-ultraviolet Imaging Spectrometer (EIS) aboard the Hinode satellite. We aim to extend our view of the velocity patterns seen in tornado-like structures with EIS to a wider range of temperatures and to provide insight into the physical characteristics of the plasma. Using Gaussian fitting to fit and de-blend the spectral lines seen by EIS, we calculated line-of-sight velocities and non-thermal line widths. Along with information from the CHIANTI database, we used line intensity ratios to calculate electron densities at each pixel. Using a regularised inversion code we also calculated the differential emission measure (DEM) at different locations in the prominence. The split Doppler-shift pattern is found to be visible down to a temperature of around log(T) = 6.0. At temperatures lower than this, the pattern is unclear in this data set. We obtain an electron density of log(n_e) = 8.5 when looking towards the centre of the tornado structure at a plasma temperature of log(T) = 6.2, as compared to the surroundings of the tornado structure where we find log(n_e) to be nearer 9. Non-thermal line widths show broader profiles at the tornado location when compared to the surrounding corona. We discuss the differential emission measure in both the tornado and the prominence body, which suggests that there is more contribution in the tornado at temperatures below log(T) = 6.0 than in the prominence.
The FERMI observation of a $\gamma$-ray excess from the galactic-centre, as well as the PAMELA, AMS, and AMS-2 anti-proton excesses, and the recent claim of a FERMI excess in the Reticulum-2 dwarf galaxy have been put forward as possible detections of neutralino dark matter. These are of particular interest as the neutralino annihilation models which fit these observations might have observable consequences from radio to $\gamma$-ray emission. Since dark matter is expected to be a major matter constituents of cosmic structure, these multi-frequency consequences should also be common to structures across the mass spectrum. Thus, in this work we make predictions for the multi-frequency spectra of three well-known sources dominated by dark matter, e.g. the Coma cluster, the galaxy M81, and the Draco dwarf galaxy using models favoured by dark matter interpretations of the aforementioned observations. We pay special attention to the consequences for these models when their cross-sections are renormalised to reproduce the recent $\gamma$-ray excess observed in the Reticulum-2 dwarf galaxy, which throw a dark matter interpretation of this excess into doubt. We find that the multi-frequency data of Coma, M81 and Draco disfavour the dark matter interpretation of the AMS, Pamela and Fermi anti-proton excess. However, models derived from FERMI galactic centre observations present no such conflicts. We determine the detection prospects of the Square Kilometre Array, the Cherenkov Telescope Array, as well as the ASTROGAM and ASTRO-H satellites for the studied models. This demonstrates that ASTRO-H is well positioned to probe the X-ray emissions from neutralino annihilation. Thus, multi-frequency observation with the next generation experiments will allow for unprecedented sensitivity to the neutralino parameter space.
Observing a high-statistics neutrino signal from the supernova explosions in the Galaxy is a major goal of low-energy neutrino astronomy. The prospects for detecting all flavors of neutrinos and antineutrinos from the core-collapse supernova (ccSN) in operating and forthcoming large liquid scintillation detectors (LLSD) are widely discussed now. One of proposed LLSD is Baksan Large Volume Scintillation Detector (BLVSD). This detector will be installed at the Baksan Neutrino Observatory (BNO) of the Institute for Nuclear Research, Russian Academy of Sciences, at a depth of 4800 m.w.e. Low-energy neutrino astronomy is one of the main lines of research of the BLVSD.
A massive black hole (MBH) consumes stars whose orbits evolve into the small phase-space volume of unstable orbits, the "loss-cone", which take them directly into the MBH, or close enough to interact strongly with it. The resulting phenomena: tidal heating and tidal disruption, binary capture and hyper-velocity star ejection, gravitational wave (GW) emission by inspiraling compact remnants, or hydrodynamical interactions with an accretion disk, are of interest as they can produce observable signatures and thereby reveal the existence of the MBH, affect its mass and spin evolution, probe strong gravity, and provide information on stars and gas near the MBH. The continuous loss of stars and the processes that resupply them shape the central stellar distribution. We investigate relativistic stellar dynamics near the loss-cone of a non-spinning MBH in steady-state analytically and by Monte Carlo simulations of the diffusion of the orbital parameters. These take into account Newtonian mass precession due to enclosed stellar mass, in-plane precession due to general relativity, dissipation by GW, uncorrelated two-body relaxation, correlated resonant relaxation (RR) and adiabatic invariance due to secular precession, using a rigorously derived description of correlated post-Newtonian dynamics in the diffusion limit. We argue that general maximal entropy considerations strongly constrain orbital diffusion in steady-state, irrespective of the relaxation mechanism. We identify the exact phase-space separatrix between plunges and inspirals, predict their steady-state rates, and verify they are robust under a wide range of assumptions. We derive the dependence of the rates on the mass of the MBH, show that the contribution of RR is small, and discuss special cases where unquenched RR in restricted volumes of phase-space may affect the steady-state substantially.
This paper summarizes a presentation given on the occasion of the inauguration of the High Altitude Water Cherenkov (HAWC) Gamma-ray Observatory in Puebla, Mexico in March 2015. The inauguration of a new facility for the study of astrophysical gamma-rays provides an excellent opportunity to review the technical evolution and the scientific achievements of VERITAS (the Very Energetic Radiation Imaging Telescope Array System) since its own inauguration in 2007. HAWC and VERITAS are separated by only 14 degrees in longitude, and so can view much of the same sky at the same time. In combination with other ground-based facilities, and with the instruments onboard the Fermi Gamma-ray Space Telescope, VERITAS and HAWC will give an unprecedented view of the gamma-ray sky. We provide an overview of VERITAS, and discuss the complementarity of the two observatories for future gamma-ray observations.
Spectropolarimetric observations have been used to map stellar magnetic fields, many of which display strong bands of azimuthal fields that are toroidal. A number of explanations have been proposed to explain how such fields might be generated though none are definitive. In this paper, we examine the toroidal fields of a sample of 55 stars with magnetic maps, with masses in the range 0.1-1.5$\,{\rm M}_\odot$. We find that the energy contained in toroidal fields has a power law dependence on the energy contained in poloidal fields. However the power index is not constant across our sample, with stars less and more massive than 0.5$\,{\rm M}_\odot$ having power indices of 0.72$\pm$0.08 and 1.25$\pm$0.06 respectively. There is some evidence that these two power laws correspond to stars in the saturated and unsaturated regimes of the rotation-activity relation. Additionally, our sample shows that strong toroidal fields must be generated axisymmetrically. The latitudes at which these bands appear depend on the stellar rotation period with fast rotators displaying higher latitude bands than slow rotators. The results in this paper present new constraints for future dynamo studies.
A new formula for the gravitational potential of flattened systems is proposed. It is a modification of the Miyamoto-Nagai potential and should be applied to very flattened systems, exponential discs as a typical example. The resulting rotation curve agrees sufficiently well with that obtained by using special functions and the total masses remain the same. The functions contained in the new term can improve the agreement for the rotation curve and also reduce the effect of negative density values which appear off the midplane.
The first one started in 2012 July 16 (MJD 56124) and the second one in 2013 April 9 (MJD 56391). The multiwavelength data analysis shows that the $\gamma$-ray flare observed in 2012 was not detected in the hard-X ray bands. This result is usually interpreted as an "orphan" flare. In 2013, the analysis of the multiwavelength light curves shows that there are two very bright states detected in the optical R-band. The first one in 2013 April 9 (R =11.74 $\pm$ 0.04) and the second one in May 12 (R =11.62 $\pm$ 0.04). Also, high activity states were detected in the soft and hard X-rays. A discrete correlation function analysis of this last flare shows a strong correlation between the GeV $\gamma$-rays and the optical/hard-X ray emission. These results are discussed in terms of the more adequate standard scenarios that could explain the multiwavelength variations displayed by this blazar.
We investigate, in terms of production from pulsars and their nebulae, the cosmic ray positron and electron fluxes above $\sim10$ GeV, observed by the AMS-02 experiment up to 1 TeV. We concentrate on the Vela-X case. Starting from the gamma-ray photon spectrum of the source, generated via synchrotron and inverse Compton processes, we estimated the electron and positron injection spectra. Several features are fixed from observations of Vela-X and unknown parameters are borrowed from the Crab nebula. The particle spectra produced in the pulsar wind nebula are then propagated up to the Solar System, using a diffusion model. Differently from previous works, the omnidirectional intensity excess for electrons and positrons is obtained as a difference between the AMS-02 data and the corresponding local interstellar spectrum. An equal amount of electron and positron excess is observed and we interpreted this excess (above $\sim$100 GeV in the AMS-02 data) as a supply coming from Vela-X. The particle contribution is consistent with models predicting the gamma-ray emission at the source. The input of a few more young pulsars is also allowed, while below $\sim$100 GeV more aged pulsars could be the main contributors.
We address the thermal history of the Earth after the Moon-forming impact, taking tidal heating and thermal blanketing by the atmosphere into account. The atmosphere sets an upper bound of ~100 W/m^2 on how quickly the Earth can cool. The liquid magma ocean cools over 2-10 Myrs, with longer times corresponding to high angular-momentum events. Tidal heating is focused mostly in mantle materials that are just beginning to freeze. The atmosphere's control over cooling sets up a negative feedback between viscosity-dependent tidal heating and temperature-dependent viscosity of the magma ocean. While the feedback holds, evolution of the Moon's orbit is limited by the modest radiative cooling rate of Earth's atmosphere. Orbital evolution is orders of magnitude slower than in conventional constant Q models, which promotes capture by resonances. The evection resonance is encountered early, when the Earth is molten. Capture by the evection resonance appears certain but unlikely to generate much eccentricity because it is encountered early when the Earth is molten and Q_Earth >> Q_Moon. Tidal dissipation in the Earth becomes more efficient (Q_Earth << Q_Moon) later when the Moon is between ~20 R_Earth and ~40 R_Earth. If lunar eccentricity grew great, this was when it did so, perhaps setting the table for some other process to leave its mark on the inclination of the Moon.
The detection of high-energy astrophysical neutrinos of extraterrestrial origin by the IceCube neutrino observatory in Antarctica has opened a unique window to the cosmos that may help to probe both the distant Universe and our cosmic backyard. The arrival directions of these high-energy events have been interpreted as uniformly distributed on the celestial sphere. Here, we revisit the topic of the putative isotropic angular distribution of these events applying Monte Carlo techniques to investigate a possible anisotropy. A modest evidence for anisotropy is found. An excess of events appears projected towards a section of the Local Void, where the density of galaxies with radial velocities below 3000 km/s is rather low, suggesting that this particular group of somewhat clustered sources are located either very close to the Milky Way or perhaps beyond 40 Mpc. The results of further analyses of the subsample of southern hemisphere events favour an origin at cosmological distances with the arrival directions of the events organized in a fractal-like structure. Although a small fraction of closer sources is possible, remote hierarchical structures appear to be the main source of these very energetic neutrinos. Some of the events may have their origin at the IBEX ribbon.
Dissipation of magnetohydrodynamic (MHD) wave energy has been proposed as a viable heating mechanism in the solar chromospheric plasma. Here, we use a simplified one-dimensional model of the chromosphere to theoretically investigate the physical processes and the spatial scales that are required for the efficient dissipation of Alfv\'en waves and slow magnetoacoustic waves. We consider the governing equations for a partially ionized hydrogen-helium plasma in the single-fluid MHD approximation and include realistic wave damping mechanisms that may operate in the chromosphere, namely Ohmic and ambipolar magnetic diffusion, viscosity, thermal conduction, and radiative losses. We perform an analytic local study in the limit of small amplitudes to approximately derive the lengthscales for critical damping and efficient dissipation of MHD wave energy. We find that the critical dissipation lengthscale for Alfv\'en waves depends strongly on the magnetic field strength and ranges from 10~m to 1~km for realistic field strengths. The damping of Alfv\'en waves is dominated by Ohmic diffusion for weak magnetic field and low heights in the chromosphere, and by ambipolar diffusion for strong magnetic field and medium/large heights in the chromosphere. Conversely, the damping of slow magnetoacoustic waves is less efficient, and spatial scales shorter than 10~m are required for critical damping. Thermal conduction and viscosity govern the damping of slow magnetoacoustic waves and play an equally important role at all heights. These results indicate that the spatial scales at which strong wave heating may work in the chromosphere are currently unresolved by observations.
We present a dynamical analysis of the merging galaxy cluster system Abell 2146 using spectroscopy obtained with the Gemini Multi-Object Spectrograph on the Gemini North telescope. As revealed by the Chandra X-ray Observatory, the system is undergoing a major merger and has a gas structure indicative of a recent first core passage. The system presents two large shock fronts, making it unique amongst these rare systems. The hot gas structure indicates that the merger axis must be close to the plane of the sky and that the two merging clusters are relatively close in mass, from the observation of two shock fronts. Using 63 spectroscopically determined cluster members, we apply various statistical tests to establish the presence of two distinct massive structures. With the caveat that the system has recently undergone a major merger, the virial mass estimate is M_vir = 8.5 +4.3 -4.7 x 10 ^14 M_sol for the whole system, consistent with the mass determination in a previous study using the Sunyaev-Zeldovich signal. The newly calculated redshift for the system is z = 0.2323. A two-body dynamical model gives an angle of 13-19 degrees between the merger axis and the plane of the sky, and a timescale after first core passage of 0.24-0.28 Gyr.
Over the past decades, General Relativity and the concordance $\Lambda$CDM model have been successfully tested using several different astrophysical and cosmological probes based on large datasets ({\it precision cosmology}). Despite their successes, some shortcomings emerge due to the fact that General Relativity should be revised at infrared and ultraviolet limits and to the fact that the fundamental nature of Dark Matter and Dark Energy is still a puzzle to be solved. In this perspective, $f(R)$ gravity have been extensively investigated being the most straightforward way to modify General Relativity and to overcame some of the above shortcomings. In this paper, we review various aspects of $f(R)$ gravity at extragalactic and cosmological levels. In particular, we consider cluster of galaxies, cosmological perturbations, and N-Body simulations, focusing on those models that satisfy both cosmological and local gravity constraints. The perspective is that some classes of $f(R)$ models can be consistently constrained by Large Scale Structure.
The assumption of collisionless cold dark matter on its own cannot reconcile several astrophysical discrepancies (cusp-vs-core problem, missing satellite problem, too-big-to-fail problem). Self-interacting dark matter provides a promising framework for solving all these problems, and self-interaction cross sections are duly constrained in the literature. Following the work of Tulin, Yu, and Zurek [1], we can constrain the dark matter mass and the mass of a light mediator assuming a generic scalar Yukawa-type interaction. In particular, we constrain the strongly coupled inflationary dark matter of the luminogenesis model, a unification model with the gauge group $SU(3)_C \times SU(6) \times U(1)_Y$, which breaks to the Standard Model with an extra gauge group for dark matter when the inflaton rolls into its true vacuum. The luminogenesis model is additionally subject to constraints on inflation, and we find an upper bound on the scale of symmetry breaking of the inflaton and the decoupling scale $M_1$ of certain representations of the gauge group. We emphasize that the luminogenesis model enables a unique connection between astrophysical constraints, the nature of dark matter, and inflation.
We describe two aspects of the physics of hybrid stars that have a sharp
interface between a core of quark matter and a mantle of nuclear matter.
Firstly, we analyze the mass-radius relation. We describe a generic "Constant
Speed of Sound" (CSS) parameterization of the quark matter equation of state
(EoS), in which the speed of sound is independent of density. In terms of the
three parameters of the CSS EoS we obtain the phase diagram of possible forms
of the hybrid star mass-radius relation, and we show how observational
constraints on the maximum mass and typical radius of neutron stars can be
expressed as constraints on the CSS parameters.
Secondly, we propose a mechanism for the damping of density oscillations,
including r-modes, in hybrid stars with a sharp interface. The dissipation
arises from the periodic conversion between quark matter and nuclear matter
induced by the pressure oscillations in the star. We find the damping grows
nonlinearly with the amplitude of the oscillation and is powerful enough to
saturate an r-mode at very low saturation amplitude, of order $10^{-10}$, which
is compatible with currently-available observations of neutron star spin
frequencies and temperatures.
The largest geomagnetic storms of solar cycle 24 so far occurred on 2015 March 17 and June 22 with $D_{\rm st}$ minima of $-223$ and $-195$ nT, respectively. Both of the geomagnetic storms show a multi-step development. We examine the plasma and magnetic field characteristics of the driving coronal mass ejections (CMEs) in connection with the development of the geomagnetic storms. A particular effort is to reconstruct the in situ structure using a Grad-Shafranov technique and compare the reconstruction results with solar observations, which gives a larger spatial perspective of the source conditions than one-dimensional in situ measurements. Key results are obtained concerning how the plasma and magnetic field characteristics of CMEs control the geomagnetic storm intensity and variability: (1) a sheath-ejecta-ejecta mechanism and a sheath-sheath-ejecta scenario are proposed for the multi-step development of the 2015 March 17 and June 22 geomagnetic storms, respectively; (2) two contrasting cases of how the CME flux-rope characteristics generate intense geomagnetic storms are found, which indicates that a southward flux-rope orientation is not a necessity for a strong geomagnetic storm; and (3) the unexpected 2015 March 17 intense geomagnetic storm resulted from the interaction between two successive CMEs plus the compression by a high-speed stream from behind, which is essentially the "perfect storm" scenario proposed by \citet[][i.e., a combination of circumstances results in an event of unusual magnitude]{liu14a}, so the "perfect storm" scenario may not be as rare as the phrase implies.
An electron-tracking Compton camera (ETCC) is a detector that can determine the arrival direction and energy of incident sub-MeV/MeV gamma-ray events on an event-by-event basis. It is a hybrid detector consisting of a gaseous time projection chamber (TPC), that is the Compton-scattering target and the tracker of recoil electrons, and a position-sensitive scintillation camera that absorbs of the scattered gamma rays, to measure gamma rays in the environment from contaminated soil. To measure of environmental gamma rays from soil contaminated with radioactive cesium (Cs), we developed a portable battery-powered ETCC system with a compact readout circuit and data-acquisition system for the SMILE-II experiment. We checked the gamma-ray imaging ability and ETCC performance in the laboratory by using several gamma-ray point sources. The performance test indicates that the field of view (FoV) of the detector is about 1$\;$sr and that the detection efficiency and angular resolution for 662$\;$keV gamma rays from the center of the FoV is $(9.31 \pm 0.95) \times 10^{^-5}$ and $5.9^{\circ} \pm 0.6^{\circ}$, respectively. Furthermore, the ETCC can detect 0.15$\;\mu\rm{Sv/h}$ from a $^{137}$Cs gamma-ray source with a significance of 5$\sigma$ in 13 min in the laboratory. In this paper, we report the specifications of the ETCC and the results of the performance tests. Furthermore, we discuss its potential use for environmental gamma-ray measurements.
We propose a new scalar-tensor model which induces significant deviation from general relativity inside dense objects like neutron stars, while passing solar-system and terrestrial experiments, extending a model proposed by Damour and Esposito-Farese. Unlike their model, we employ a massive scalar field dubbed asymmetron so that it not only realizes proper cosmic evolution but also can account for the cold dark matter. In our model, asymmetron undergoes spontaneous scalarization inside dense objects, which results in reduction of the gravitational constant by a factor of order unity. This suggests that observational tests of constancy of the gravitational constant in high density phase are the effective ways to look into the asymmetron model.
We show how the correction to the calculation of the mass in the original relativistic model of a rotating star by Hartle [6], found recently [10], appears in the Newtonian limit, and that the correcting term is indeed present, albeit hidden, in the original Newtonian approach by Chandrasekhar [2].
We use crossing statistics and its generalization to determine the anisotropic direction imposed on a stochastic fields in $(2+1)$Dimension. This approach enables us to examine not only the rotational invariance of morphology but also we can determine the Gaussianity of underlying stochastic field in various dimensions. Theoretical prediction of up-crossing statistics (crossing with positive slope at a given threshold $\alpha$ of height fluctuation), $\nu^+_{\diamond}(\alpha)$, and generalized roughness function, $N^{\diamond}_{tot}(q)$, for correlation length ($\xi_{\diamond}$) and/with scaling exponent ($\gamma_{\diamond}$) anisotropies are calculated. The strategy to examine the anisotropy nature and to determine its direction is as follows: we consider a set of normal axes, and sign them $||$ (parallel) and $ \bot$ (normal) with respect to unknown anisotropic direction. Then we determine $\nu_{\diamond}^+ (\alpha)$ and $N^{\diamond}_{tot}(q)$ in both directions. The directional dependency of difference between computed results in mentioned directions are clarify. Finally we systematically recognize the anisotropy direction at $3\sigma$ confidence interval using P-value approach. In order to distinguish between nature of anisotropies, after applying a typical method in determining the scaling exponents in both mentioned directions with respect to the recognized anisotropy direction using up-crossing statistics, the kind and the ratio of correlation length anisotropy are specified. Our algorithm can be mounted with a simple software on various instruments for surface analysis, such as AFM, STM and etc.
We discuss simple models which predict the existence of significant gamma ray fluxes from dark matter annihilation. In this context the dark matter candidate is a Majorana fermion with velocity-suppressed tree-level annihilation into Standard Model fermions but unsuppressed annihilation into photons. These gamma lines can easily be distinguished from the continuum and allow us to test these models in the near future.
We suggest a structure for the vacuum comprised of a network of tightly knotted/linked flux tubes formed in a QCD-like cosmological phase transition and show that such a network can drive cosmological inflation. As the network can be topologically stable only in three space dimensions, this scenario provides a dynamical explanation for the existence of exactly three large spatial dimensions in our Universe.
We show that an inflationary slow-roll potential can be derived as an IR limit of the non-perturbative exact renormalisation group equation for a scalar field within the mean-field approximation. The result follows without having to specify a Lagrangian for the UV theory at the Planck scale. All we assume is that the theory contains a scalar mode with suppressed coupling to other UV fields. The resulting effective potential gives rise to slow-roll inflation, which is fully consistent with the recent observations.
The implications of the recent classical nonlocal generalization of Einstein's theory of gravitation for gravitational physics in the Solar System are investigated. In this theory, the nonlocal character of gravity simulates dark matter. Nonlocal gravity in the Newtonian regime involves a reciprocal kernel with three spatial parameters, of which two have already been determined from the rotation curves of spiral galaxies and the internal dynamics of clusters of galaxies. However, the short-range parameter a_0 remains to be determined. In this connection, the nonlocal contribution to the perihelion precession of a planetary orbit is estimated and a preliminary lower limit on a_0 is determined.
Links to: arXiv, form interface, find, astro-ph, recent, 1508, contact, help (Access key information)