Sco X-1 has been the subject of many multi-wavelength studies in the past, being the brightest persistent extra-solar X-ray source ever observed. Here we revisit Sco X-1 with simultaneous short cadence Kepler optical photometry and MAXI X-ray photometry over a 78 day period, as well as optical spectroscopy obtained with HERMES. We find Sco X-1 to be highly variable in all our datasets. The optical fluxes are clearly bimodal, implying the system can be found in two distinct optical states. These states are generally associated with the known flaring/normal branch X-ray states, although the flux distributions associated with these states overlap. Furthermore, we find that the optical power spectrum of Sco X-1 differs substantially between optical luminosity states. Additionally we find rms-flux relations in both optical states, but only find a linear relation during periods of low optical luminosity. The full optical/X-ray discrete correlation function displays a broad ~12.5 hour optical lag. However during the normal branch phase the X-ray and optical fluxes are anti-correlated, whilst being correlated during the flaring branch. We also performed a Cepstrum analysis on the full Kepler light curve to determine the presence of any echoes within the optical light curve alone. We find significant echo signals, consistent with the optical lags found using the discrete cross-correlation. We speculate that whilst some of the driving X-ray emission is reflected by the disk, some is absorbed and re-processed on the thermal timescale, giving rise to both the observed optical lags and optical echoes.
Astrophysical discs are often warped, that is, their orbital planes change with radius. This occurs whenever there is a non-axisymmetric force acting on the disc, for example the Lense-Thirring precession induced by a misaligned spinning black hole, or the gravitational pull of a misaligned companion. Such misalignments appear to be generic in astrophysics. The wide range of systems that can harbour warped discs - protostars, X-ray binaries, tidal disruption events, quasars and others - allows for a rich variety in the disc's response. Here we review the basic physics of warped discs and its implications.
We investigate the potential of angular auto- and cross-correlation power spectra in identifying sterile neutrino dark matter in the cosmic X-ray background. We take as reference the performance of the soon-to-be-launched eROSITA satellite. The main astrophysical background sources against sterile neutrino decays are active galactic nuclei, galaxies powered by X-ray binaries, and clusters of galaxies. While sterile neutrino decays are always subdominant in the auto-correlation power spectra, they can be efficiently enhanced when cross-correlating with tracers of the dark matter distribution such as galaxies in the 2MASS catalogues. We show that the planned four-years eROSITA all-sky survey will provide a large enough photon statistics to potentially yield very stringent constraints on the decay lifetime, enabling to firmly test the recently claimed 3.56-keV X-ray line found towards several clusters and galaxies and its decaying dark matter interpretation. However, we also show that in order to fully exploit the potential of eROSITA for dark matter searches, it is vital to overcome the shot-noise limitations inherent to galaxy catalogues as tracers for the dark matter distribution.
(abridged) Motivated by the possibility that an optical coronagraph will be put on WFIRST/AFTA, we present an exploration of the general character of the direct detectability of extrasolar giant planets (EGPs) in the optical. We quantify a planet's direct detectability by the fraction of its orbit for which it is in an observable configuration--defined to be its observability fraction ($f_{obs}$). Using a suite of Monte Carlo experiments, we study the dependence of $f_{obs}$ upon various technological and astrophysical parameters, including the inner working angle (IWA) and minimum achievable contrast ($C_{min}$) of the direct-imaging observatory; the planet's scattering phase function, geometric albedo, single-scattering albedo, radius, and distance from Earth; and the semi-major axis distribution of EGPs. We assume cloud-free, homogeneous atmospheres and calculate phase functions for a given geometric or single-scattering albedo, assuming various scattering mechanisms. We find that the often-assumed Lambertian phase function can predict significantly larger $f_{obs}$ values with respect to the more physically motivated Rayleigh phase function. For observations made with WFIRST/AFTA's baseline coronagraphic capabilities ($C_{min}\sim10^{-9}$, IWA$\sim0.2''$), Jupiter-like planets orbiting stars within 10, 30, and 50 pc from Earth have volume-averaged (assuming a uniform distribution of stars) observability fractions of ${\sim}$12%, 3%, and 0.5%, respectively. Using a plausible estimate for the occurrence rate of EGPs, we find that, in all but the most optimistic configurations, the probability a blind search will lead to a detection is low ($< 5\%$). However, with orbital parameter constraints from long-term radial-velocity campaigns and Gaia astrometry, the tools we develop in this work can be used to determine both the most promising systems to target and when to observe them.
We present an optimised variant of the halo model, designed to produce accurate matter power spectra well into the non-linear regime for a wide range of cosmological models. To do this, we introduce physically-motivated free parameters into the halo-model formalism and fit these to data from high-resolution N-body simulations. For a variety of $\Lambda$CDM and $w$CDM models the halo-model power is accurate to $\simeq 5$ per cent for $k\leq 10h\,\mathrm{Mpc}^{-1}$ and $z\leq 2$. We compare our results with recent revisions of the popular HALOFIT model and show that our predictions are more accurate. An advantage of our new halo model is that it can be adapted to account for the effects of baryonic feedback on the power spectrum. We demonstrate this by fitting the halo model to power spectra from the OWLS hydrodynamical simulation suite via parameters that govern halo internal structure. We are able to fit all feedback models investigated at the 5 per cent level using only two free parameters, and we place limits on the range of these halo parameters for feedback models investigated by the OWLS simulations. Accurate predictions to high-$k$ are vital for weak lensing surveys, and these halo parameters could be considered nuisance parameters to marginalise over in future analyses to mitigate uncertainty regarding the details of feedback, the limits we find on these parameters provide a prior. Finally we investigate how lensing observables predicted by our model compare to those from simulations and from HALOFIT for a range of $k$-cuts and feedback models and quantify the angular scales at which these effects become important. Code to calculate power spectra from the model presented in this paper can be found at https://github.com/alexander-mead/hmcode.
We present 3D-Barolo, a new code that derives rotation curves of galaxies from emission-line observations. This software fits 3D tilted-ring models to spectroscopic data-cubes and can be used with a variety of observations: from HI and molecular lines to optical/IR recombination lines. We describe the structure of the main algorithm and show that it performs much better than the standard 2D approach on velocity fields. A number of successful applications, from high to very low spatial resolution data are presented and discussed. 3D-Barolo can recover the true rotation curve and estimate the intrinsic velocity dispersion even in barely resolved galaxies (about 2 resolution elements) provided that the signal to noise of the data is larger that 2-3. It can also be run automatically thanks to its source-detection and first-estimate modules, which make it suitable for the analysis of large 3D datasets. These features make 3D-Barolo a uniquely useful tool to derive reliable kinematics for both local and high-redshift galaxies from a variety of different instruments including the new-generation IFUs, ALMA and the SKA pathfinders.
A growing body of evidence indicates that the star formation rate per unit stellar mass (sSFR) decreases with increasing mass in normal "main-sequence" star forming galaxies. Many processes have been advocated as responsible for such a trend (also known as mass quenching), e.g., feedback from active galactic nuclei (AGNs), and the formation of classical bulges. We determine a refined star formation versus stellar mass relation in the local Universe. To this aim we use the Halpha narrow-band imaging follow-up survey (Halpha3) of field galaxies selected from the HI Arecibo Legacy Fast ALFA Survey (ALFALFA) in the Coma and Local superclusters. By complementing this local determination with high-redshift measurements from the literature, we reconstruct the star formation history of main-sequence galaxies as a function of stellar mass from the present epoch up to z=3. In agreement with previous studies, our analysis shows that quenching mechanisms occur above a threshold stellar mass M_knee that evolves with redshift as propto (1+z)^{2}. Moreover, visual morphological classification of individual objects in our local sample reveals a sharp increase in the fraction of visually-classified strong bars with mass, hinting that strong bars may contribute to the observed downturn in the sSFR above M_knee. We test this hypothesis using a simple but physically-motivated numerical model for bar formation, finding that strong bars can rapidly quench star formation in the central few kpc of field galaxies. We conclude that strong bars contribute significantly to the red colors observed in the inner parts of massive galaxies, although additional mechanisms are likely required to quench the star formation in the outer regions of massive spiral galaxies. Intriguingly, when we extrapolate our model to higher redshifts, we successfully recover the observed redshift evolution for M_knee.
Context. Young loose nearby associations are unique samples of close (<150
pc), young (approx 5-100 Myr) pre-main sequence (PMS) stars. A significant
number of members of these associations have been identified in the SACY
collaboration. We can use the proximity and youth of these members to
investigate key ingredients in star formation processes, such as multiplicity.
Aims. We present the statistics of identified multiple systems from 113
confirmed SACY members. We derive multiplicity frequencies, mass-ratio, and
physical separation distributions in a consistent parameter space, and compare
our results to other PMS populations and the field.
Methods. We have obtained adaptive-optics assisted near-infrared observations
with NACO (ESO/VLT) and IRCAL (Lick Observatory) for at least one epoch of all
113 SACY members. We have identified multiple systems using co-moving
proper-motion analysis and using contamination estimates. We have explored
ranges in projected separation and mass-ratio of a [3-1000 au], and q [0.1-1],
respectively.
Results. We have identified 31 multiple systems (28 binaries and 3 triples).
We derive a multiplicity frequency (MF) of MF_(3-1000au)=28.4 +4.7, -3.9% and a
triple frequency (TF) of TF_(3-1000au)=2.8 +2.5, -0.8% in the separation range
of 3-1000 au. We do not find any evidence for an increase in the MF with
primary mass. The estimated mass-ratio of our statistical sample (with
power-law index gamma=-0.04 +/- 0.14) is consistent with a flat distribution
(gamma = 0).
Conclusions. We show further similarities (but also hints of discrepancies)
between SACY and the Taurus region: flat mass-ratio distributions and
statistically similar MF and TF values. We also compared the SACY sample to the
field (in the separation range of 19-100 au), finding that the two
distributions are indistinguishable, suggesting a similar formation mechanism.
Hydrodynamic simulations of galaxies that include stellar feedback are now able to comfortably reproduce observations of the cool halo gas around star-forming galaxies. Without implementing quasar feedback, however, these simulations under-predict the cool gas content of quasar host halos by more than a factor of two. To better understand the source of this tension, we have compiled an exhaustive sample of 195 quasars at z=1 with constraints on chemically enriched, cool gas traced by MgII absorption in background quasar spectra from the Sloan Digital Sky Survey. We find a strong correlation between quasar luminosity and cool gas covering fraction. In particular, low-luminosity quasars exhibit a mean gas covering fraction comparable to inactive galaxies of similar masses, but more luminous quasars exhibit excess cool gas approaching what is reported previously at z=2.2. Moreover, 40% of the MgII absorption occurs at radial velocities of |v|>300 km/s from the quasar systemic redshift, inconsistent with gas gravitationally bound to the quasar host halo. We discuss possible scenarios to explain the large velocity offsets and observed luminosity dependence of the cool gas near quasars including gas arising in: (1) neighboring halos correlated through large-scale structure at Mpc scales, (2) feedback from luminous quasars, and (3) debris from the mergers thought to trigger luminous quasars. The first of these scenarios can be ruled out by the lack of correlation between quasar luminosity and clustering while the latter two each make distinct predictions that can be tested with additional observations.
We develop a Bayesian hierarchical modelling approach for cosmic shear power spectrum inference, jointly sampling from the posterior distribution of the cosmic shear field and its (tomographic) power spectra. Inference of the shear power spectrum is a powerful intermediate product for a cosmic shear analysis, since it requires very few model assumptions and can be used to perform inference on a wide range of cosmological models \emph{a posteriori} without loss of information. We show that joint posterior for the shear map and power spectrum can be sampled effectively by Gibbs sampling, iteratively drawing samples from the map and power spectrum, each conditional on the other. This approach neatly circumvents difficulties associated with complicated survey geometry and masks that plague frequentist power spectrum estimators, since the power spectrum inference provides prior information about the field in masked regions at every sampling step. We demonstrate this approach for inference of tomographic shear $E$-mode, $B$-mode and $EB$-cross power spectra from a simulated galaxy shear catalogue with a number of important features; galaxies distributed on the sky and in redshift with photometric redshift uncertainties, realistic random ellipticity noise for every galaxy and a complicated survey mask. The obtained posterior distributions for the tomographic power spectrum coefficients recover the underlying simulated power spectra for both $E$- and $B$-modes, where the latter are recovered at a level of $1$-$2$ orders of magnitude below the ellipticity noise level.
Upcoming space-based surveys such as Euclid and WFIRST-AFTA plan to measure Baryonic Acoustic Oscillations (BAOs) in order to study dark energy. These surveys will use IR slitless grism spectroscopy to measure redshifts of a large number of galaxies over a significant redshift range. In this paper, we use the WFC3 Infrared Spectroscopic Parallel Survey (WISP) to estimate the expected number of Halpha (Ha) emitters observable by these future surveys. WISP is an ongoing HST slitless spectroscopic survey, covering the 0.8-1.65micron wavelength range and allowing the detection of Ha emitters up to z~1.5 and [OIII] emitters to z~2.3. We derive the Ha-[OIII] bivariate line luminosity function for WISP galaxies at z~1 using a maximum likelihood estimator that properly accounts for uncertainties in line luminosity measurement, and demonstrate how it can be used to derive the Ha luminosity function from exclusively fitting [OIII] data. Using the z~2 [OIII] line luminosity function, and assuming that the relation between Ha and [OIII] luminosity does not change significantly over the redshift range, we predict the Ha number counts at z~2 - the upper end of the redshift range of interest for the future surveys. For the redshift range 0.7<z<2, we expect ~3000 galaxies/deg^2 for a flux limit of 3x10^{-16} ergs/s/cm^2 (the proposed depth of Euclid galaxy redshift survey) and ~20,000 galaxies/deg^2 for a flux limit of ~10^{-16} ergs/s/cm^2 (the baseline depth of WFIRST galaxy redshift survey).
We study the characteristic size and shape of idealized blazar-induced cascade halos in the $1-100 \, {\rm GeV}$ energy range assuming various non-helical and helical configurations for the intergalactic magnetic field (IGMF). While the magnetic field creates an extended halo, the helicity provides the halo with a twist. Under simplifying assumptions, we assess the parameter regimes for which it is possible to measure the size and shape of the halo from a single source and then to deduce properties of the IGMF. We find that blazar halo measurements with an experiment similar to Fermi-LAT are best suited to probe a helical magnetic field with strength and coherence length today in the ranges $10^{-17} \lesssim B_{0} / {\rm Gauss} \lesssim 10^{-13}$ and $10 \, {\rm Mpc} \lesssim \lambda \lesssim 10 \, {\rm Gpc}$ where $\mathcal{H} \sim B_0^2 / \lambda$ is the magnetic helicity density. Stronger magnetic fields or smaller coherence scales can still potentially be investigated, but the connection between the halo morphology and the magnetic field properties is more involved. Weaker magnetic fields or longer coherence scales require high photon statistics or superior angular resolution.
Astro-H will be able for the first time to map gas velocities and detect turbulence in galaxy clusters. One of the best targets for turbulence studies is the Coma cluster, due to its proximity, absence of a cool core, and lack of a central active galactic nucleus. To determine what constraints Astro-H will be able to place on the Coma velocity field, we construct simulated maps of the projected gas velocity and compute the second-order structure function, an analog of the velocity power spectrum. We vary the injection scale, dissipation scale, slope, and normalization of the turbulent power spectrum, and apply measurement errors and finite sampling to the velocity field. We find that even with sparse coverage of the cluster, Astro-H will be able to measure the Mach number and the injection scale of the turbulent power spectrum--the quantities determining the energy flux down the turbulent cascade and the diffusion rate for everything that is advected by the gas (metals, cosmic rays, etc). Astro-H will not be sensitive to the dissipation scale or the slope of the power spectrum in its inertial range, unless they are outside physically-motivated intervals. We give the expected confidence intervals for the injection scale and the normalization of the power spectrum for a number of possible pointing configurations, combining the structure function and velocity dispersion data. Importantly, we also determine that measurement errors on the line shift will bias the velocity structure function upward, and show how to correct this bias.
This thesis assesses the influence of astronomical phenomena on the Earth's biosphere and climate. I examine in particular the relevance of both the path of the Sun through the Galaxy and the evolution of the Earth's orbital parameters in modulating non-terrestrial mechanisms. I build models to predict the extinction rate of species, the temporal variation of the impact cratering rate and ice sheet deglaciations, and then compare these models with other models within a Bayesian framework. I find that the temporal distribution of mass extinction events over the past 550 Myr can be explained just as well by a uniform random distribution as by other models, such as variations in the stellar density local to the Sun arising from the Sun's orbit. Given the uncertainties in the Galaxy model and the Sun's current phase space coordinates, as well as the errors in the geological data, it is not possible to draw a clear connection between terrestrial extinction and the solar motion. In a separate study, I find that the solar motion, which modulates the Galactic tidal forces imposed on Oort cloud comets, does not significantly influence this cratering rate. My dynamical models, together with the solar apex motion, can explain the anisotropic perihelia of long period comets without needing to invoke the existence of a Jupiter-mass solar companion. Finally, I find that variations in the Earth's obliquity play a dominant role in triggering terrestrial deglaciations over the past 2 Myr. The precession of the equinoxes, in contrast, only becomes important in pacing large deglaciations after the transition from the 100-kyr dominant periodicity in the ice coverage to a 41-kyr dominant periodicity, which occurred 0.7 Myr ago.
New mid-infrared period-luminosity (PL) relations are presented for \rrl{} variables in the globular cluster M4 (NGC 6121). Accurate photometry was obtained for 37 \rrl{} variables using observations from the Infrared Array Camera onboard the Spitzer Space Telescope. The dispersion of M4's PL relations is 0.056, and the uncertainty in the slope is 0.11 mag. Additionally, weestablished calibrated PL relations at 3.6 and 4.5~\micron{} using published Hubble Space Telescopegeometric parallaxes of five Galactic \rrl{} stars. The resulting band-averaged distance modulus for M4 is $ \mu = 11.399 \pm 0.007 \textrm{(stat)} \pm 0.080 \textrm{(syst)} \pm 0.015 \textrm{(cal)} \pm 0.020 \textrm{(ext)}$. The systematic uncertainty will be greatly reduced when parallaxes of more stars become available from the GAIA mission. Optical and infrared period-color (PC) relations are also presented, and the lack of a MIR PC relation suggests that \rrl{} stars are not affected by CO absorptionin the 4.5~\micron{} band.
Absorption of infrared and visible radiation from stellar emission spectra by the negatively charged hydrogen ions H$^{-}$ is considered. The explicit formula for the photodetachment cross-section of the negatively charged hydrogen ion(s) is derived. Photodetachemnt cross-sections of the ${}^{\infty}$H$^{-}$, ${}^{3}$H$^{-}$ (or T$^{-}$), ${}^{2}$H$^{-}$ (or D$^{-}$) and ${}^{1}$H$^{-}$ ions are determined to high accuracy and for a large number of photo-electron momenta/energies. We introduce criteria which can be used to evaluate the overall quality of highly accurate wave functions of the hydrogen ion(s). One of these criteria is based on highly accurate calculations of the lowest order QED corrections in the negatively charged hydrogen ions, including ${}^{1}$H$^{-}$ (protium), ${}^{2}$H$^{-}$ (deuterium), ${}^{3}$H$^{-}$ (tritium) and model ion with the infinitely heavy nucleus ${}^{\infty}$H$^{-}$. An effective approach has been developed to calculate three-body integrals with the Bessel functions of different orders. Some preliminary evaluations of the phototdetachment cross-sections of the negatively charged hydrogen ions are performed. Inverse bremsstrahlung in the field of the neutral hydrogen atom is briefly discussed.
We present the measurements and modelling of the small-to-intermediate scale (0.1--25 Mpc/h) projected and three-dimensional (3D) redshift-space two-point correlation functions (2PCFs) of local galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7. We find a clear dependence of galaxy clustering on luminosity in both projected and redshift spaces, generally being stronger for more luminous samples. The measurements are successfully interpreted within the halo occupation distribution (HOD) framework with central and satellite velocity bias parameters to describe galaxy kinematics inside haloes and to model redshift-space distortion (RSD) effects. In agreement with previous studies, we find that more luminous galaxies reside in more massive haloes. Including the redshift-space 2PCFs helps tighten the HOD constraints. Moreover, we find that luminous central galaxies are not at rest at the halo centres, with the velocity dispersion about 30% that of the dark matter. Such a relative motion may reflect the consequence of galaxy and halo mergers, and we find that central galaxies in lower mass haloes tend to be more relaxed with respect to their host haloes. The motion of satellite galaxies in luminous samples is consistent with their following that of the dark matter. For faint samples, satellites tends to have slower motion, with velocity dispersion inside haloes about 85% that of the dark matter. We discuss possible applications of the velocity bias constraints on studying galaxy evolution and cosmology. In the appendix, we characterize the distribution of galaxy redshift measurement errors, which is well described by a Gaussian-convolved double exponential distribution.
We present near-infrared observations of the environments around three radio-loud sources (MG1 J0442+0202, 3C 068.2, and MS 1426.9+1052) at redshifts z=1.10,1.57, and 1.83 (respectively), that are surrounded by near-infrared galaxy overdensities. Overdensities with respect to field counts were found to be significant up to 19-sigma, with twelve times the expected number of galaxies within the inner regions of the densest proto-cluster. Color-magnitude relations are constructed in K_s, J-K_s, with each candidate cluster exhibiting a feature consistent with the beginnings of a red sequence. Galaxy models based on the redshift of the radio source are used to compare expected color-magnitude relations for a given formation epoch with the observed red sequence of each candidate, and are found to be consistent with an old (z_f > 5) formation epoch for a few bright, red galaxies on the red sequence.
We show that interferometry can be applied to study irregular, rapidly rotating structures, as are expected in the turbulent accretion flow near a black hole. Specifically, we analyze the lagged covariance between interferometric baselines of similar lengths but slightly different orientations. We demonstrate that the peak in the lagged covariance indicates the direction and angular velocity of the flow. Importantly, measuring the direction of the flow as clockwise or counterclockwise on the sky breaks a degeneracy in accretion disk inclinations when analyzing time-averaged images alone. We explore the potential efficacy using three-dimensional, general relativistic magnetohydrodynamic (GRMHD) simulations, and we highlight several baseline pairs for the Event Horizon Telescope (EHT) that are well-suited to this application. These results indicate that the EHT is capable of determining the direction and angular velocity of the emitting material near Sgr A*, even for highly-inclined flows, and they suggest that a rotating flow may even be utilized to improve imaging capabilities.
Cosmological simulations of the low-density intergalactic medium exhibit a strikingly tight power-law relation between temperature and density that holds over two decades in density. It is found that this relation should roughly apply Delta z ~ 1-2 after a reionization event, and this limiting behavior has motivated the power-law parameterizations used in most analyses of the Ly-alpha forest. This relation has been explained by using equations linearized in the baryonic overdensity (which does not address why a tight power-law relation holds over two decades in density) or by equating the photoheating rate with the cooling rate from cosmological expansion (which we show is incorrect). Previous explanations also did not address why recombination cooling and Compton cooling off of the cosmic microwave background, which are never negligible, do not alter the character of this relation. We provide an understanding for why a tight power-law relation arises for unshocked gas at all densities for which collisional cooling is unimportant. We also use our results to comment on (1) how quickly fluctuations in temperature redshift away after reionization processes, (2) how much shock heating occurs in the low-density intergalactic medium, and (3) how the temperatures of collapsing gas parcels evolve.
We present results from ab initio simulations of liquid water-hydrogen mixtures in the range from 2 to 70 GPa and from 1000 to 6000 K, covering conditions in the interiors of ice giant planets and parts of the outer envelope of gas giant planets. In addition to computing the pressure and the internal energy, we derive the Gibbs free energy by performing a thermodynamic integration. For all conditions under consideration, our simulations predict hydrogen and water to mix in all proportions. The thermodynamic behavior of the mixture can be well described with an ideal mixing approximation. We suggest a substantial fraction of water and hydrogen in giant planets may occur in homogeneously mixed form rather than in separate layers. The extend of mixing depends on the planet's interior dynamics and its conditions of formation, in particular on how much hydrogen was present when icy planetesimals were delivered. Based on our results, we do not predict water-hydrogen mixtures to phase separate during any stage of the evolution of giant planets. We also show that the hydrogen content of an exoplanet is much higher if the mixed interior is assumed.
Halos, filaments, sheets and voids in the cosmic web can be defined in terms of the eigenvalues of the smoothed shear tensor and a threshold $\lambda_{\rm th}$. Using analytic methods, we construct mean maps centered on each of these four types of structures for Gaussian random fields corresponding to cosmological initial conditions. Each map also requires a choice of shear at the origin; we consider two choices motivated by average shear calculations and find characteristic sizes, shapes and other properties of the basic objects in the mean maps. We explore how these mean map properties change with varying the threshold and smoothing scale, i.e. with varying the separation of the cosmic web into different kinds of components. The mean maps become increasingly complex as the threshold $\lambda_{\rm th}$ decreases to zero. We also describe scatter around these mean maps, subtleties which can arise in their construction, and some comparisons between them and the full ensemble of corresponding structures at final times.
In the second paper of this series, we present results from cosmological simulations on the demographics of flyby interactions to gauge their potential impact on galaxy evolution. In a previous paper, we demonstrated that flybys -- an interaction where two independent halos inter-penetrate but detach at a later time and do not merge -- occur much more frequently than previously believed. In particular, we found that the frequency of flybys increases at low redshift and is comparable to or even greater than the frequency of mergers for halos $\gtrsim 10^{11} M_\odot/h$. In this paper, we classify flybys according to their orbits and the level of perturbation exacted on both the halos involved. We find that the majority of flybys penetrate deeper than $\sim R_{half}$ of the primary and have an initial relative speed $\sim 1.6\times V_{vir}$ of the primary. The typical flyby mass-ratio is $\sim 0.1$ at high $z$ for all halos, while at low $z$, massive primary halos undergo flybys with small secondary halos. We estimate the perturbation from the flyby on both the primary and the secondary and find that a typical flyby is mostly non-perturbative for the primary halo. However, since a massive primary experiences so many flybys at any given time, they are nearly continually a victim of a perturbative event. In particular, we find flybys that cause $\sim 1\%$ change in the binding energy of a primary halo occurs $\gtrsim 1 $ Gyr$^{-1}$ for halos $> 10^{10} M_\odot/h$ for $z \lesssim 4$. Secondary halos, on the other hand, are highly perturbed by the typical encounter, experiencing a change in binding energy of nearly order unity. Our results imply that flybys can drive a significant part of galaxy transformation at moderate to lower redshifts ($z \lesssim 4$). We touch on implications for observational surveys, mass-to-light ratios, and galaxy assembly bias.
Recent discoveries of recombining plasmas (RPs) in supernova remnants (SNRs) have dramatically changed our understanding of SNR evolution. To date, a dozen of RP SNRs have been identified in the Galaxy. Here we present Suzaku deep observations of four SNRs in the Large Magellanic Cloud (LMC), N49, N49B, N23, and DEM L71, for accurate determination of their plasma state. Our uniform analysis reveals that only N49 is in the recombining state among them, which is the first robust discovery of a RP from an extra-galactic SNR. Given that RPs have been identified only in core-collapse SNRs, our result strongly suggests a massive star origin of this SNR. On the other hand, no clear evidence for a RP is confirmed in N23, from which detection of recombination lines and continua was previously claimed. Comparing the physical properties of the RP SNRs identified so far, we find that all of them are categorized into the "mixed-morphology" class and interacting with surrounding molecular clouds. This might be a key to solve formation mechanisms of the RPs.
Wrapping around the Milky Way, the Sagittarius stream is the dominant substructure in the halo. Our statistical selection method has allowed us to identify 106 highly likely members of the Sagittarius stream. Spectroscopic analysis of metallicity and kinematics of all members provides us with a new mapping of the Sagittarius stream. We find correspondence between the velocity distribution of stream stars and those computed for a triaxial model of the Milky Way dark matter halo. The Sagittarius trailing arm exhibits a metallicity gradient, ranging from $-0.59$ dex to $-0.97$ dex over 142$^{\circ}$. This is consistent with the scenario of tidal disruption from a progenitor dwarf galaxy that possessed an internal metallicity gradient. We note high metallicity dispersion in the leading arm, causing a lack of detectable gradient and possibly indicating orbital phase mixing. We additionally report on a potential detection of the Sextans dwarf spheroidal in our data.
To study the usefulness of \CFP\ as a tracer of the regions where C\p\ and \HH\ coexist in the interstellar medium. We used the Plateau de Bure Interferometer to synthesize \CFP\ J=1-0 absorption at 102.6 GHz toward the core of the distant HII region W49N at l = 43.2\degr, b=0.0\degr, and we modeled the fluorine chemistry in diffuse/translucent molecular gas. We detected \CFP\ absorption over a broad range of velocity showing that \CFP\ is widespread in the \HH-bearing Galactic disk gas. Originally detected in dense gas in the Orion Bar and Horsehead PDR, \CFP\ was subsequently detected in absorption from diffuse and translucent clouds seen toward \bll\ and 3C111. Here we showed that \CFP\ is distributed throughout the diffuse and translucent molecular disk gas with N(\CFP)/N(\HH) $= 1.5-2.0\times10^{-10}$, increasing to N(\CFP)/N(\HH) $= 3.5\times10^{-10}$ in one cloud at 39 \kms\ having higher N(\HH) $\approx 3\times10^{21}\pcc$. Models of the fluorine chemistry reproduce the observed column densities and relative abundance of HF, from which \CFP\ forms, but generally overpredict the the column density of \CFP\ by factors of 1.4-4. We show that a free space photodissociation rate $\Gamma \ga 10^{-9}\ps$, comparable to that of CH, might account for much of the discrepancy but a recent calculation finds a value about ten times smaller. In the heavily blended and kinematically complex spectra seen toward W49, \CFP\ absorption primarily traces the peaks of the \HH\ distribution.
We perform a discrimination procedure with the spectral index diagram of TiO5 and CaH2+CaH3 to separate M giants from M dwarfs. Using the M giant spectra identified from the LAMOST DR1 with high signal-to-noise ratio (SNR), we have successfully assembled a set of M giant templates, which show more reliable spectral features. Combining with the M dwarf/subdwarf templates in Zhong et al. (2015), we present an extended M-type templates library which includes not only M dwarfs with well-defined temperature and metallicity grid but also M giants with subtype from M0 to M6. Then, the template-fit algorithm were used to automatically identify and classify M giant stars from the LAMOST DR1. The result of M giant stars catalog is cross-matched with 2MASS JHKs and WISE W1/W2 infrared photometry. In addition, we calculated the heliocentric radial velocity of all M giant stars by using the cross-correlation method with the template spectrum in a zero-velocity restframe. Using the relationship between the absolute infrared magnitude MJ and our classified spectroscopic subtype, we derived the spectroscopic distance of M giants with uncertainties of about 40%. A catalog of 8639 M giants is provided. As an additional search result, we also present 101690 M dwarfs/subdwarfs catalog which were classified by our classification pipeline.
We describe the data processing pipeline developed to reduce the pointing observation data of Lunar-based Ultraviolet Telescope (LUT), which belongs to the Chang'e-3 mission of the Chinese Lunar Exploration Program. The pointing observation program of LUT is dedicated to variable objects monitoring in a near-ultraviolet (245-345 nm) band. LUT works in lunar daytime for sufficient power supply, so some special data processing strategies have been developed in the pipeline. The procedures of the pipeline mainly include stray light removing, astrometry, flat fielding employing superflat technique, source extraction and cosmic rays rejection, aperture and PSF photometry, aperture correction, and catalogues archiving. It has been intensively tested and works smoothly with observation data. The photometric accuracy is typically ~0.02 mag for LUT 10 mag stars (30s exposure), with errors from background noises, residuals of stray light removing, and flat fielding. The accuracy degrades to be ~0.2 mag for stars of 13.5 mag which is the 5{\sigma} detection limit of LUT.
Both the Cosmic Ray Flux (CRF) and Solar Energetic Particles (SEPs) have left an imprint on SOHO technical systems. While the solar array efficiency degraded irreversibly down to ~77% of its original level over roughly 1 1/2 solar cycles, Single Event Upsets (SEUs) in the solid state recorder (SSR) have been reversed by the memory protection mechanism. We compare the daily CRF observed by the Oulu station with the daily SOHO SEU rate and with the Degradation curve of the solar arrays. The Oulu CRF and the SOHO SSR SEU rate are both modulated by the solar cycle and are highly correlated, except for sharp spikes in the SEU rate, caused by isolated SEP events, which also show up as discontinuities in the otherwise slowly decreasing solar ray efficiency. This allows to discriminate between effects with solar and non-solar origin and to compare the relative strength of both. We find that during solar cycle 23 (1996 Apr 1 -- 2008 Aug 31) only 6% of the total number of SSR SEUs were caused by SEPs; the remaining 94% were due to galactic cosmic rays. During the maximum period of cycle 23 (2000 Jan 1 -- 2003 Dec 31), the SEP contribution increased to 22%, and during 2001, the year with the highest SEP rate, to 30%. About 40% of the total solar array degradation during the 17 years from Jan 1996 through Feb 2013 can be attributed to proton events, i.e. the effect of a series of short-lived, violent SEP events is comparable to the cycle-integrated damage by cosmic rays.
A very important issue concerning protostellar jets is the mechanism behind their formation. Obtaining information on the region at the base of a jet can shed light into the subject and some years ago this has been done through a search for a rotational signature at the jet line spectrum. The existence of such signatures, however, remains controversial. In order to contribute to the clarification of this issue, in this paper we show that the Principal Component Analysis (PCA) can potentially help to distinguish between rotation and precession effects in protostellar jet images. We apply the PCA to synthetic spectro-imaging datacubes generated as an output of numerical simulations of protostellar jets. In this way we generate a benchmark to which a PCA diagnostics of real observations can be confronted. Using the computed emission line profiles for [O I]6300A and [S II]6716A, we recover and analyze the effects of rotation and precession in tomograms generated by PCA. We show that different combinations of the eigenvectors can be used to enhance and to identify the rotation features present in the data. Our results indicate that the PCA can be useful for disentangling rotation from precession in jets with an inclination of the jet with respect to the plane of the sky as high as 45 degrees. We have been able to recover the initially imposed rotation jet profile for models at moderate inclination angle (< 15 degrees) and without precession (abridged).
Astrometric monitoring of the nearby early-L dwarf DE0823$-$49 has revealed a low-mass companion in a 248-day orbit that was announced in an earlier work. Here, we present new astrometric and spectroscopic observations that allow us to characterise the system in detail. The optical spectrum shows LiI-absorption indicative of a young age and/or substellar mass for the primary component. The near-infrared spectrum is best reproduced by a binary system of brown dwarfs with spectral types of L1.5 $+$ L5.5 and effective temperatures of $2150\pm100$ K and $1670\pm140$ K. To conform with the photocentric orbit size measured with astrometry and the current understanding of substellar evolution, the system must have an age in the 80--500 Myr range. Evolutionary models predict component masses in the ranges of $M_1\simeq0.028-0.063\,M_\odot$ and $M_2\simeq0.018-0.045\,M_\odot$ with a mass ratio of $q\simeq0.64-0.74$. Multi-epoch radial velocity measurements unambiguously establish the three-dimensional orbit of the system and allow us to investigate its kinematic properties. DE0823$-$49 emerges as a rare example of a nearby brown dwarf binary with orbit, component properties, and age that are characterised well. It is a juvenile resident of the solar neighbourhood, but does not appear to belong to a known young association or moving group.
In-situ spacecraft observations recently suggested that the transport of energetic particles accelerated at heliospheric shocks can be anomalous, i.e. the mean square displacement can grow non-linearly in time. In particular, a new analysis technique has permitted the study of particle transport properties from energetic particle time profiles upstream of interplanetary shocks. Indeed, the time/spatial power laws of the differential intensity upstream of several shocks are indicative of superdiffusion. A complete determination of the key parameters of superdiffusive transport comprises the power-law index, the superdiffusion coefficient, the related transition scale at which the energetic particle profiles turn to decay as power laws, and the energy spectral index of the shock accelerated particles. Assuming large-scale spatial homogeneity of the background plasma, the power-law behaviour can been derived from both a (microscopic) propagator formalism and a (macroscopic) fractional transport equation. We compare the two approaches and find a relation between the diffusion coefficients used in the two formalisms. Based on the assumption of superdiffusive transport, we quantitatively derive these parameters by studying energetic particle profiles observed by the Ulysses and Voyager 2 spacecraft upstream of shocks in the heliosphere, for which a superdiffusive particle transport has previously been observed. Further, we have jointly studied the electron energy spectra, comparing the values of the spectral indices observed with those predicted by the standard diffusive shock acceleration theory and by a model based on superdiffusive transport. For a number of interplanetary shocks and for the solar wind termination shock, for the first time we obtain the anomalous diffusion constants and the scale at which the probability of particle free paths changes to a power-law...
EX Lupi is a young star, prototype of EXor variables. Its spectrum is very rich in emission lines, including many metallic lines. It has been also proposed to have a close companion. We use the metallic emission lines to study the accretion structures and to test the companion hypothesis. We analyse 54 spectra taken in 5 years of quiescence time. We study the line profile variability and the radial velocity of the metallic emission lines. We use the velocity signatures of different species with various excitation conditions and their time dependency to track the dynamics associated to accretion. We observe periodic velocity variations in the line components consistent with rotational modulation. The modulation is stronger for lines with higher excitation potentials. We propose that the narrow line components are produced in the post-shock region, while the broad components originate in the more extended, pre-shock material. All the emission lines suffer velocity modulation due to the rotation of the star. The broad components are responsible for the line-dependent veiling observed in EX Lupi. Rotationally-modulated line-dependent veiling can explain the radial velocity signatures, making the close-in companion hypothesis unnecessary. The accretion structure is locked to the star and very stable during the 5 years of observations. Not all stars with similar spectral types and accretion rates show the same metallic emission lines, maybe related to differences in temperature and density in their accretion structure(s). The contamination of photospheric lines by accretion processes can be turned into a very useful tool to determine the innermost details of the accretion channels in the proximities of the star. Emission lines from very stable accretion columns will nevertheless be a very strong limitation for the detection of companions by radial velocity in young stars. (Abridged)
There has been growing observational evidence for the presence of condensates in the atmospheres and/or comet-like tails of extrasolar planets. As a result, systematic and homogeneous tables of dust properties are useful in order to facilitate further observational and theoretical studies. In this paper we present calculations and analysis of non-isotropic phase functions, asymmetry parameter (mean cosine of the scattering angle), absorption and scattering opacities, single scattering albedos, equilibrium temperatures, and radiative accelerations of dust grains relevant for extrasolar planets. Our assumptions include spherical grain shape, Deirmendjian particle size distribution, and Mie theory. We consider several species: corundum/alumina, iron, olivines with 0% and 50% iron content, pyroxenes with 0%, 20% and 60% iron content, carbon at two different temperatures, water ice, liquid water, and ammonia. The presented tables cover the wavelength range of 0.2 to 500 micron and modal particle radii from 0.01 micron to 100 micron. Equilibrium temperatures and radiative accelerations assume irradiation by a non-black-body source of light with temperatures from 7000K to 700K seen at solid angles from 2$\pi$ to $10^{-6}$ sr. The tables are provided to the community together with a simple code which allows for an optional, finite, angular dimension of the source of light (star) in the phase function.
We present a description of the pipeline used to calibrate the Planck Low Frequency Instrument (LFI) timelines into thermodynamic temperatures for the Planck 2015 data release, covering 4 years of uninterrupted operations. As in the 2013 data release, our calibrator is provided by the spin-synchronous modulation of the CMB dipole, exploiting both the orbital and solar components. Our 2015 LFI analysis provides an independent Solar dipole estimate in excellent agreement with that of HFI and within $1\sigma$ (0.3 % in amplitude) of the WMAP value. This 0.3 % shift in the peak-to-peak dipole temperature from WMAP and a global overhaul of the iterative calibration code increases the overall level of the LFI maps by 0.45 % (30 GHz), 0.64 % (44 GHz), and 0.82 % (70 GHz) in temperature with respect to the 2013 Planck data release, thus reducing the discrepancy with the power spectrum measured by WMAP. We estimate that the LFI calibration uncertainty is at the level of 0.20 % for the 70 GHz map, 0.26 % for the 44 GHz map, and 0.35 % for the 30 GHz map. We provide a detailed description of the impact of all the changes implemented in the calibration since the previous data release.
The Gaia-ESO Survey (GES) is a large public spectroscopic survey that aims at observing with FLAMES@VLT the main stellar components of our Galaxy. The study of the population of open clusters is one of the main objectives of GES. We present some results from the first 18 months of observations, among them, a preliminary view of the radial metallicity gradient as traced by open cluster data and a comparison of the chemical patterns of clusters located in different parts of the Galactic disk.
We compute trajectories of dust grains starting from a homogeneous surface activity-profile. Despite the homogeneous initial distribution a collimation in jet-like structures becomes visible. The fine structure is caused by topographical features with similar bundles of normal vectors. Gravitational forces are accurately determined from the triangular surface mesh. For the comet 67P/Churyumov-Gerasimenko, we find several areas of good agreement between the homogeneous dust emission model and Rosetta observation of dust jets.
In 2014 the Royal Astronomical Society carried out a survey of its membership, finding that we are both more and less diverse than UK society as a whole. Robert Massey summarizes the findings and what they mean for the Society in future.
4C +21.35 is a flat-spectrum-radio-quasar-type blazar, in which the rapid variability of very high energy (VHE, $E_{\gamma}\gtrsim 100$\,GeV) emission as short as $\sim$ 10 minutes was observed by MAGIC Cherenkov telescopes, and the VHE spectrum extends up to at least 400\,GeV. In this paper, by using a flat broad-line region (BLR) structure, we study the location and properties of $\gamma$-ray emitting region of 4C +21.35 under the constraints of multiwavelength data. We fit three quasi-simultaneous spectral energy distributions (SEDs) using homogeneous one-zone leptonic model, in which the flat BLR with the aperture angle of $\alpha=25^{\circ}$ and a spherically symmetric hot dusty torus with the temperature of $T_{\rm sub}=1200$\,K, are assumed. The results show that the jet structure of 4C +21.35 is nearly conical with a half-opening angle of $\theta_{\rm j}\simeq 0.29^{\circ}-0.6^{\circ}$. Furthermore, the emitting region is located within the BLR clouds and approaches to outer radius of the BLR during the flaring states, while it is well beyond the dusty torus in quiescent state. The quiescent high-energy emission is dominated by synchrotron self-Compton (SSC) process, the high-energy emission during the flaring periods is dominated by Compton scattering of BLR and dusty torus photons with the value of Compton-dominance parameter is about 30. Moreover, the fit to optical/ultraviolet data provides a further support that the central black hole (BH) mass of 4C +21.35 is $6\times 10^{8}$\,$M_{\odot}$.
We report a detailed investigation of the bulk motions of the nearby Galactic stellar disk, based on three samples selected from the LSS-GAC DR2: a global sample containing 0.57 million FGK dwarfs out to $\sim$ 2 kpc, a local subset of the global sample consisting $\sim$ 5,400 stars within 150 pc, and an anti-center sample containing $\sim$ 4,400 AFGK dwarfs and red clump stars within windows of a few degree wide centered on the Galactic anti-center. The global sample is used to construct a three-dimensional map of bulk motions of the Galactic disk from the solar vicinity out to $\sim$ 2 kpc with a spatial resolution of $\sim$ 250 pc. Typical values of the radial and vertical components of bulk motion range from $-$15 km s$^{-1}$ to 15 km s$^{-1}$, while the lag behind the circular speed dominates the azimuthal component by up to $\sim$ 15 km s$^{-1}$. The map reveals spatially coherent, kpc-scale stellar flows in the disk, with typical velocities of a few tens km s$^{-1}$. Bending- and breathing-mode perturbations are clearly visible, and vary smoothly across the disk plane. Our data also reveal higher-order perturbations, such as breaks and ripples, in the profiles of vertical motion versus height. From the local sample, we find that stars of different populations exhibit very different patterns of bulk motion. Finally, the anti-center sample reveals a number of peaks in stellar number density in the line-of-sight velocity versus distance distribution, with the nearer ones apparently related to the known moving groups. The "velocity bifurcation" reported by Liu et al. (2012) at Galactocentric radii 10--11 kpc is confirmed. However, just beyond this distance, our data also reveal a new triple-peaked structure.
The search for weakly interacting massive particles (WIMPs) by direct detection faces an encroaching background due to coherent neutrino-nucleus scattering. As the sensitivity of these experiments improves, the question of how to best distinguish a dark matter signal from neutrinos will become increasingly important. A proposed method of overcoming this so-called 'neutrino floor' is to utilize the directional signature that both neutrino and dark matter induced recoils possess. We show that directional experiments can indeed probe WIMP-nucleon cross-sections below the neutrino floor with little loss in sensitivity due to the neutrino background. In particular we find at low WIMP masses (around 6 GeV) the discovery limits for directional detectors penetrate below the non-directional limit by several orders of magnitude. For high WIMP masses (around 100 GeV), the non-directional limit is overcome by a factor of a few. Furthermore we show that even for directional detectors which can only measure 1- or 2-dimensional projections of the 3-dimensional recoil track, the discovery potential is only reduced by a factor of 3 at most. We also demonstrate that while the experimental limitations of directional detectors, such as sense recognition and finite angular resolution, have a detrimental effect on the discovery limits, it is still possible to overcome the ultimate neutrino background faced by non-directional detectors.
We present a catalog of 908 objects observed with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) in the vicinity fields of M31 and M33, targeted as globular clusters (GCs) and candidates. The targets include known GCs and candidates selected from the literature, as well as new candidates selected from the Sloan Digital Sky Survey (SDSS). Analysis shows that 356 of them are likely GCs of various degree of confidence, while the remaining ones turn out to be background galaxies and quasars, stars and HII regions in M31 or foreground Galactic stars. The 356 likely GCs include 298 bona fide GCs and 26 candidates known in the literature. Three candidates selected from the Revised Bologna Catalog of M31 GCs and candidates (RBC) and one possible cluster from Johnson et al. are confirmed to be bona fide clusters. We search for new GCs in the halo of the M31 amongst the new candidates selected from the SDSS photometry. Based on radial velocities yielded by LAMSOT spectra and visual examination of the SDSS images, we find 28 objects, 5 bona fide and 23 likely GCs. Amongst the five bona fide GCs, three have been recently discovered independently by others, the remaining 25 are our new identifications, including two bona fide ones. The new identified objects fall at projected distances ranging from 13 to 265 kpc from M31. Of the two newly discovered bona fide GCs, one is located near M33, probably a GC belonging to M33. The other bona fide GC falls on the Giant Stream with a projected distance of 78 kpc from M31. Of the 23 newly identified likely GCs, one has a projected distance of about 265 kpc from M31 and could be an intergalactic cluster.
We use 297 042 main sequence turn-off stars selected from the LSS-GAC to determine the radial and vertical gradients of stellar metallicity of the Galactic disk in the anti-center direction. We determine ages of those turn-off stars by isochrone fitting and measure the temporal variations of metallicity gradients. Our results show that the gradients, both in the radial and vertical directions, exhibit significant spatial and temporal variations. The radial gradients yielded by stars of oldest ages (>11 Gyr) are essentially zero at all heights from the disk midplane, while those given by younger stars are always negative. The vertical gradients deduced from stars of oldest ages (>11Gyr) are negative and show only very weak variations with the Galactocentric distance in the disk plane, $R$, while those yielded by younger stars show strong variations with $R$. After being essentially flat at the earliest epochs of disk formation, the radial gradients steepen as age decreases, reaching a maxima (steepest) at age 7-8 Gyr, and then they flatten again. Similar temporal trends are also found for the vertical gradients. We infer that the assemblage of the Milky Way disk may have experienced at least two distinct phases. The earlier phase is probably related to a slow, pressure-supported collapse of gas, when the gas settles down to the disk mainly in the vertical direction. In the later phase, there are significant radial flows of gas in the disk, and the rate of gas inflow near the solar neighborhood reaches a maximum around a lookback time of 7-8 Gyr. The transition of the two phases occurs around a lookback time between 8 and 11 Gyr. The two phases may be responsible for the formation of the Milky Way thick and thin disks, respectively. And, as a consequence, we recommend that stellar age is a natural, physical criterion to distinguish thin and thick disk stars. ... (abridged)
In this work, we report new quasars discovered in the vicinity fields of the Andromeda (M31) and Triangulum (M33) galaxies with the LAMOST (Large Sky Area Multi-Object Fiber Spectroscopic Telescope, also named Guoshoujing Telescope) during the 2013 observational season, the second year of Regular Survey. In total, 1330 new quasars are discovered in an area of $\sim$133 deg$^2$ around M31 and M33. With $i$ magnitudes ranging from 14.79 to 20.0, redshifts from 0.08 to 4.85, the 1330 new quasars represent a significant increase of the number of identified quasars in the vicinity fields of M31 and M33. Up to the moment, there are in total 1870 quasars discovered by LAMOST in this area (see also Huo et al. 2010, 2013). The much enlarged sample of known quasars in this area can potentially be utilized to construct a precise astrometric reference frame for the measurement of the minute proper motions of M31, M33 and the associated substructures, vital for the understanding of the formation and evolution of M31, M33 and the Local Group of galaxies.Meanwhile, amongst the sample, there are in total 45, 98 and 225 quasars with $i$ magnitudes brighter than 17.0, 17.5 and 18.0 respectively. In the aforementioned brightness bins, 15, 35 and 84 quasars are reported here for the first time, 6, 21 and 81 are reported in Huo et al. (2010, 2013), while 0, 1 and 6 are from the Sloan Digital Sky Survey, and 24, 41 and 54 are from the NED database. These bright quasars provide an invaluable sample for the kinematics and chemistry study of the interstellar/intergalactic medium of the Local Group.
Using a sample of over 70, 000 red clump (RC) stars with $5$-$10$% distance accuracy selected from the LAMOST Spectroscopic Survey of the Galactic Anti-center (LSS-GAC), we study the radial and vertical gradients of the Galactic disk(s) mainly in the anti-center direction, covering a significant disk volume of projected Galactocentric radius $7 \leq R_{\rm GC} \leq 14$ kpc and height from the Galactic midplane $0 \leq |Z| \leq 3$ kpc. Our analysis shows that both the radial and vertical metallicity gradients are negative across much of the disk volume probed, and exhibit significant spatial variations. Near the solar circle ($7 \leq R_{\rm GC} \leq 11.5$ kpc), the radial gradient has a moderately steep, negative slope of $-0.08$ dex kpc$^{-1}$ near the midplane ($|Z| < 0.1$ kpc), and the slope flattens with increasing $|Z|$. In the outer disk ($11.5 < R_{\rm GC} \leq 14$ kpc), the radial gradients have an essentially constant, much less steep slope of $-0.01$ dex kpc $^{-1}$ at all heights above the plane, suggesting that the outer disk may have experienced an evolution path different from that of the inner disk. The vertical gradients are found to flatten largely with increasing $R_{\rm GC}$. However, the vertical gradient of the lower disk ($0 \leq |Z| \leq 1$ kpc) is found to flatten with $R_{\rm GC}$ quicker than that of the upper disk ($1 < |Z| \leq 3$ kpc). Our results should provide strong constraints on the theory of disk formation and evolution, as well as the underlying physical processes that shape the disk (e.g. gas flows, radial migration, internal and external perturbations).
The multi-parameter character of supersymmetric dark-matter models implies the combination of their experimental studies with astrophysical and cosmological probes. The physics of the early Universe provides nontrivial effects of non-equilibrium particles and primordial cosmological structures. Primordial black holes (PBHs) are a profound signature of such structures that may arise as a cosmological consequence of supersymmetric (SUSY) models. SUSY-based mechanisms of baryosynthesis can lead to the possibility of antimatter domains in a baryon asymmetric Universe. In the context of cosmoparticle physics, which studies the fundamental relationship of the micro- and macro-worlds, the development of SUSY illustrates the main principles of this approach, as the physical basis of the modern cosmology provides cross-disciplinary tests in physical and astronomical studies.
Polarimetric studies of astrophysical sources can make important contributions to resolve the geometry of the emitting region and determine the photon emission mechanism. PoGOLite is a balloon-borne polarimeter operating in the hard X-ray band (25-240 keV), with a Pathfinder mission focussing on Crab observations. Within the polarimeter, the distribution of Compton scattering angles is used to determine the polarisation fraction and angle of incident photons. To assure an unbiased measurement of the polarisation during a balloon flight it is crucial to characterise the performance of the instrument before the launch. This paper presents the results of the PoGOLite calibration tests and simulations performed before the 2013 balloon flight. The tests performed confirm that the polarimeter does not have any intrinsic asymmetries and therefore does not induce bias into the measurements. Generally, good agreement is found between results from test data and simulations which allows the polarimeter performance to be estimated for Crab observations.
Titan's polar surface is dotted with hundreds of lacustrine depressions. Based on the hypothesis that they are karstic in origin, we aim at determining the efficiency of surface dissolution as a landshaping process on Titan, in a comparative planetology perspective with the Earth as reference. Our approach is based on the calculation of solutional denudation rates and allow inference of formation timescales for topographic depressions developed by chemical erosion on both planetary bodies. The model depends on the solubility of solids in liquids, the density of solids and liquids, and the average annual net rainfall rates. We compute and compare the denudation rates of pure solid organics in liquid hydrocarbons and of minerals in liquid water over Titan and Earth timescales. We then investigate the denudation rates of a superficial organic layer in liquid methane over one Titan year. At this timescale, such a layer on Titan would behave like salts or carbonates on Earth depending on its composition, which means that dissolution processes would likely occur but would be 30 times slower on Titan compared to the Earth due to the seasonality of precipitation. Assuming an average depth of 100 m for Titan's lacustrine depressions, these could have developed in a few tens of millions of years at polar latitudes higher than 70{\deg} N and S, and a few hundreds of million years at lower polar latitudes. The ages determined are consistent with the youth of the surface (<1 Gyr) and the repartition of dissolution-related landforms on Titan.
The Spitzer Space Telescope mapped the Perseus molecular cloud complex with IRAC and MIPS as part of the c2d Spitzer Legacy project. This paper combines the observations from both instruments giving an overview of low-mass star formation across Perseus from 3.6 to 70 micron. We provide an updated list of young stellar objects with new classifications and source fluxes from previous works, identifying 369 YSOs in Perseus with the Spitzer dataset. By synthesizing the IRAC and MIPS maps of Perseus and building on the work of previous papers in this series (Jorgensen et al. 2006, Rebull et al. 2007), we present a current census of star formation across the cloud and within smaller regions. 67% of the YSOs are associated with the young clusters NGC 1333 and IC 348. The majority of the star formation activity in Perseus occurs in the regions around the clusters, to the eastern and western ends of the cloud complex. The middle of the cloud is nearly empty of YSOs despite containing regions of high visual extinction. The western half of Perseus contains three-quarters of the total number of embedded YSOs (Class 0+I and Flat SED sources) in the cloud and nearly as many embedded YSOs as Class II and III sources. Class II and III greatly outnumber Class 0+I objects in eastern Perseus and IC 348. These results are consistent with previous age estimates for the clusters. Across the cloud, 56% of YSOs and 91% of the Class 0+I and Flat sources are in areas where Av > 5 mag, indicating a possible extinction threshold for star formation.
In preparing the way for the Square Kilometre Array and its pathfinders, there is a pressing need to begin probing the transient sky in a fully robotic fashion using the current generation of radio telescopes. Effective exploitation of such surveys requires a largely automated data-reduction process. This paper introduces an end-to-end automated reduction pipeline, AMIsurvey, used for calibrating and imaging data from the Arcminute Microkelvin Imager Large Array. AMIsurvey makes use of several component libraries which have been packaged separately for open-source release. The most scientifically significant of these is chimenea, which implements a telescope agnostic algorithm for automated imaging of pre-calibrated multi-epoch radio-synthesis data, making use of CASA subroutines for the underlying image-synthesis operations. At a lower level, AMIsurvey relies upon two libraries, drive-ami and drive-casa, built to allow use of mature radio-astronomy software packages from within Python scripts. These packages may be of wider techical interest, since they demonstrate use of the Python library pexpect to emulate terminal interaction with an external process. This approach allows for rapid development of a Python interface to any legacy or externally maintained pipeline which accepts command-line input, without requiring alterations to the original code.
We develop and apply the position angle-only shear estimator of Whittaker et al. (2014) to realistic galaxy images. This is done by demonstrating the method on the simulations of the GREAT3 challenge (Mandelbaum et al. 2014b), which include contributions from anisotropic PSFs. We measure the position angles of the galaxies using three distinct methods - the integrated light method, quadrupole moments of surface brightness, and using model-based ellipticity measurements provided by IM3SHAPE. A weighting scheme is adopted to address biases in the position angle measurements which arise in the presence of an anisotropic PSF. Biases on the shear estimates, due to measurement errors on the position angles and correlations between the measurement errors and the true position angles, are corrected for using simulated galaxy images and an iterative procedure. The properties of the simulations are estimated using the deep field images provided as part of the challenge. A method is developed to match the distributions of galaxy fluxes and half-light radii from the deep fields to the corresponding distributions in the field of interest. We recover angle-only shear estimates with a performance close to current well-established model and moments-based methods for all three angle measurement techniques. The Q-values for all three methods are found to be Q~400. The code is freely available online at this http URL
Fractal grains have large surface area, which leads to more efficient condensation. The special limit case where the volume-area ratio is constant (corresponding to, e.g., a very rough grain surface or non-compacts aggregates) is particularly interesting, as well as convenient, from a mathematical point of view. If dust grains from AGB stars have `rough surfaces', it may have important implications for our understanding of dust and wind formation in AGB stars.
It is shown that high-energy astrophysical neutrinos observed in the IceCube experiment can be produced by protons accelerated in extragalactic Type IIn supernova remnants by shocks propagating in the dense circumstellar medium. The nonlinear diffusive shock acceleration model is used for description of particle acceleration.
FRBs are puzzling, millisecond, energetic radio transients with no discernible source; observations show no counterparts in other frequency bands. The birth of a quark star from a parent neutron star experiencing a quark nova - previously thought undetectable when born in isolation - provides a natural explanation for the emission characteristics of FRBs. The generation of unstable r-process elements in the quark nova ejecta provides millisecond exponential injection of electrons into the surrounding strong magnetic field at the parent neutron star's light cylinder via $\beta$-decay. This radio synchrotron emission has a total duration of hundreds of milliseconds and matches the observed spectrum while reducing the inferred dispersion measure by approximately 200 cm$^{-3}$ pc. The model allows indirect measurement of neutron star magnetic fields and periods in addition to providing astronomical measurements of $\beta$-decay chains of unstable neutron rich nuclei. Using this model, we can calculate expected FRB average energies ($\sim$ 10$^{41}$ ergs), spectra shapes and provide a theoretical framework for determining distances.
We present results for the estimation of gravity, effective temperature, and radial velocity of poorly studied chemically peculiar stars recently observed with the spectropolarimeter Echelle SpectroPolarimetric Device for Observations of Stars at the Canada-France-Hawaii Telescope in the frame of the Vertical Stratification of Element Abundances project. The effective temperature and surface gravity values are determined for the very first time for four of the stars from our sample (HD23878, HD83373, HD95608, and HD164584). Grids of stellar atmosphere models with the corresponding fluxes have been calculated using version 15 of the PHOENIX code for effective temperatures in the range of 5000-15,000 K, for the logarithm of surface gravities in the range of 3.0-4.5 and for the metallicities from -1.0 to +1.5. We used these fluxes to fit the Balmer line profiles employing the code FITSB2 that produces estimates of the effective temperature, gravity, and radial velocity for each star. When possible, our results are compared to those previously published. The physical characteristics of 16 program stars are discussed with the future aim to study the abundance anomalies of chemical species and the possible vertical abundance stratification in their stellar atmosphere.
A modified gravitational theory is developed in which the gravitational coupling constants $G$ and $Q$ and the effective mass $m_\phi$ of a repulsive vector field run with momentum scale $k$ or length scale $\ell =1/k$, according to a renormalization group flow. The theory can explain cosmological early universe data with a dark hidden photon and late time galaxy and cluster dynamics without dark matter. The theory agrees with solar system and binary pulsar observations.
With a view to developing a formalism that will be applicable at second perturbative order, we devise a new practical scheme for computing the gravitational self-force experienced by a point mass moving in a curved background spacetime. Our method works in the frequency domain and employs the effective-source approach, in which a distributional source for the retarded metric perturbation is replaced with an effective source for a certain regularized self-field. A key ingredient of the calculation is the analytic determination of an appropriate puncture field from which the effective source and regularized residual field can be calculated. In addition to its application in our effective-source method, we also show how this puncture field can be used to derive tensor-harmonic mode-sum regularization parameters that improve the efficiency of the traditional mode-sum procedure. To demonstrate the method, we calculate the first-order-in-the-mass-ratio self-force and redshift invariant for a point mass on a circular orbit in Schwarzschild spacetime.
We propose new large field inflation scenarios built on the framework of F-term axion monodromy. Our setup is based on string compactifications where D-branes create potentials for closed string axions via F-terms. Because the source of the axion potential is different from the standard sources of moduli stabilisation, it is possible to lower the inflaton mass as compared to other massive scalars. We discuss a particular class of models based on type IIA flux compactifications with D6-branes. In the small field regime they describe supergravity models of quadratic chaotic inflation with a stabiliser field. In the large field regime the inflaton potential displays a flattening effect due to Planck suppressed corrections, allowing to easily fit the cosmological parameters of the model within current experimental bounds.
An epoch of Higgs relaxation may occur in the early universe during or immediately following postinflationary reheating. It has recently been pointed out that leptogenesis may occur in minimal extensions of the Standard Model during this epoch. We analyse Higgs relaxation taking into account the effects of perturbative and non-perturbative decays of the Higgs condensate, and we present a detailed derivation of the relevant kinetic equations and of the relevant particle interaction cross sections. We identify the parameter space in which a sufficiently large asymmetry is generated.
We discuss the detection of gravitational-wave backgrounds in the context of Bayesian inference and suggest a practical definition of what it means for a signal to be considered stochastic---namely, that the Bayesian evidence favors a stochastic signal model over a deterministic signal model. A signal can further be classified as Gaussian-stochastic if a Gaussian signal model is favored. In our analysis we use Bayesian model selection to choose between several signal and noise models for simulated data consisting of uncorrelated Gaussian detector noise plus a superposition of sinusoidal signals from an astrophysical population of gravitational-wave sources. For simplicity, we consider co-located and co-aligned detectors with white detector noise, but the method can be extended to more realistic detector configurations and power spectra. The general trend we observe is that a deterministic model is favored for small source numbers, a non-Gaussian stochastic model is preferred for intermediate source numbers, and a Gaussian stochastic model is preferred for large source numbers. However, there is very large variation between individual signal realizations, leading to fuzzy boundaries between the three regimes. We find that a hybrid, trans-dimensional model comprised of a deterministic signal model for individual bright sources and a Gaussian-stochastic signal model for the remaining confusion background outperforms all other models in most instances.
Major discrepancies have been noted for some time between fossil ages and molecular divergence dates for a variety of taxa. Recently, systematic trends within avian clades have been uncovered. The trends show that the disparity is much larger for mitochondrial DNA than for nuclear DNA; also that it is larger for crown fossil dates than stem fossil dates. It was argued that this pattern is largely inconsistent with incompleteness of the fossil record as the principal driver of the disparity. A case is presented that given the expected mutations from a fluctuating background of astrophysical radiation from such sources as supernovae, the rate of molecular clocks is variable and should increase into the past. This is a possible explanation for the disparity. One test of this hypothesis is to look for an acceleration of molecular clocks 2 to 2.5 Ma due to a probable moderately nearby supernova at that time. Another is to look for reduced disparity in benthic organisms of the deep ocean.
Interactions of different types of topological defects can play an important role in the aftermath of a phase transition. We study interactions of fundamental magnetic monopoles and stable domain walls in a Grand Unified theory in which $SU(5) \times Z_2$ symmetry is spontaneously broken to $SU(3)\times SU(2)\times U(1)/Z_6$. We find that there are only two distinct outcomes depending on the relative orientation of the monopole and the wall in internal space. In one case, the monopole passes through the wall, while in the other it unwinds on hitting the wall.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
Radio emission from non-magnetic cataclysmic variables (CVs, accreting white
dwarfs) could allow detailed studies of outflows and possibly accretion flows
in these nearby, numerous and non-relativistic compact accretors. Up to now,
however, very few CVs have been detected in the radio.
We have conducted a VLA pilot survey of four close and optically-bright
novalike CVs at 6 GHz, detecting three, and thereby doubling the number of
radio detections of these systems. RW Sex, V603 Aql and the old nova TT Ari
were detected in both of the epochs, while V1084 Her was not detected (to a
$3\sigma$ upper-limit of 7.8 $\mu\rm{Jy}\,\rm{beam}^{-1}$). These observations
clearly show that the sensitivity of previous surveys was typically too low to
detect these objects and that non-magnetic CVs can indeed be significant radio
emitters.
The three detected sources show a range of properties, including flaring and
variability on both short ($\sim$200 s) and longer-term (days) time-scales, as
well as circular polarization levels of up to 100\%. The spectral indices range
from steep to inverted; TT Ari shows a spectral turnover at $\sim$6.5 GHz,
while the spectral index of V603 Aql flattened from $\alpha=0.54\pm0.05$ to
$0.16\pm0.08$ ($F_{\nu}\propto \nu^{\alpha}$) in the week between observations.
This range of properties suggests that more than one emission process can be
responsible for the radio emission in non-magnetic CVs. In this sample we find
that individual systems are consistent with optically thick synchrotron
emission, gyrosynchrotron emission or cyclotron maser emission.
There has been significant controversy over the mechanisms responsible for forming compact stellar systems like ultra compact dwarfs (UCDs), with suggestions that UCDs are simply the high mass extension of the globular cluster (GC) population, or alternatively, the liberated nuclei of galaxies tidally stripped by larger companions. Definitive examples of UCDs formed by either route have been difficult to find, with only a handful of persuasive examples of stripped-nucleus type UCDs being known. In this paper we present very deep Gemini/GMOS spectroscopic observations of the suspected stripped nucleus UCD NGC 4546-UCD1 taken in good seeing conditions (< 0.7"). With these data we examine the spatially resolved kinematics and star formation history of this unusual object. We find no evidence of a rise in the central velocity dispersion of the UCD, suggesting that this UCD lacks a massive central black hole like those found in some other compact stellar systems, a conclusion confirmed by detailed dynamical modelling. Finally we are able to use our extremely high signal to noise spectrum to detect a temporally extended star formation history for this UCD. We find that the UCD was forming stars since the earliest epochs until at least 1-2 Gyr ago. Taken together these observations confirm that NGC 4546-UCD1 is the remnant nucleus of a nucleated dwarf galaxy that was tidally destroyed by NGC 4546 within the last 1-2 Gyr.
We perform uniformly sampled large-scale cosmological simulations including magnetic fields with the moving mesh code AREPO. We run two sets of MHD simulations: one including adiabatic gas physics only; the other featuring the fiducial feedback model of the Illustris simulation. In the adiabatic case, the magnetic field amplification follows the $B \propto \rho^{2/3}$ scaling derived from `flux-freezing' arguments, with the seed field strength providing an overall normalisation factor. At high baryon overdensities the amplification is enhanced by shear flows and turbulence. Feedback physics and the inclusion of radiative cooling change this picture dramatically. Gas collapses to much larger densities and the magnetic field is amplified strongly, reaching saturation and losing memory of the initial seed field. At lower densities a dependence on the seed field strength and orientation, which in principle can be used to constrain models of cosmological magnetogenesis, is still present. Inside the most massive haloes magnetic fields reach values of $\sim 10-100\,\mu G$, in agreement with galaxy cluster observations. The topology of the field is tangled and gives rise to rotation measure signals in reasonable agreement with the observations. However, the rotation measure signal declines too rapidly towards larger radii compared to observational data.
Knots or blobs observed in astrophysical jets are commonly interpreted as shock waves moving along them. Long time observations of the HST-1 knot inside the jet of the galaxy M87 have produced detailed multi-wavelength light curves. In this article, we model these light curves using the semi-analytical approach developed by Mendoza et al. (2009). This model was developed to account for the light curves of working surfaces moving along relativistic jets. These working surfaces are generated by periodic oscillations of the injected flow velocity and mass ejection rates at the base of the jet. Using genetic algorithms to fit the parameters of the model, we are able to explain the outbursts observed in the light curves of the HST-1 knot with an accuracy greater than a 2 sigma statistical confidence level.
We present simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR ) and Suzaku observations of the X-ray binary Cygnus X-1 in the hard state. This is the first time this state has been observed in Cyg X-1 with NuSTAR, which enables us to study the reflection and broad-band spectra in unprecedented detail. We confirm that the iron line cannot be fit with a combination of narrow lines and absorption features, and instead requires a relativistically blurred profile in combination with a narrow line and absorption from the companion wind. We use the reflection models of Garcia et al. (2014) to simultaneously measure the black hole spin, disk inner radius, and coronal height in a self-consistent manner. Detailed fits to the iron line profile indicate a high level of relativistic blurring, indicative of reflection from the inner accretion disk. We find a high spin, a small inner disk radius, and a low source height, and rule out truncation to greater than three gravitational radii at the 3{\sigma} confidence level. In addition, we find that the line profile has not changed greatly in the switch from soft to hard states, and that the differences are consistent with changes in the underlying reflection spectrum rather than the relativistic blurring. We find that the blurring parameters are consistent when fitting either just the iron line or the entire broad-band spectrum, which is well modelled with a Comptonized continuum plus reflection model.
We used recently produced Solar System ephemeris, which incorporate two years of ranging observations to the MESSENGER spacecraft, to extract the secular orbital elements for Mercury and associated uncertainties. As Mercury is in a stable 3:2 spin-orbit resonance these values constitute an important reference for the planet's measured rotational parameters, which in turn strongly bear on physical interpretation of Mercury's interior structure. In particular, we derive an mean orbital period of 87.96934962 $\pm$ 0.00000037 days and (assuming the perfect resonance) a spin rate of 6.138506839 $\pm$ 0.000000028 degree/day. The difference between this rotation rate and the currently adopted rotation rate (Archinal et al, 2011) corresponds to a longitudinal displacement of approx. 67 m per year at the equator. Moreover, we present a basic approach for the calculation of the orientation of the instantaneous Laplace and Cassini planes of Mercury. The analysis allows us to assess the uncertainties in physical parameters of the planet when derived from observations of Mercury's rotation.
In this paper, we present the FATS (Feature Analysis for Time Series) library. FATS is a Python library which facilitates and standardizes feature extraction for time series data. In particular, we focus on one application: feature extraction for astronomical light curve data, although the library is generalizable for other uses. We detail the methods and features implemented for light curve analysis, and present examples for its usage.
Dark matter in the Milky Way may annihilate directly into gamma rays, producing a monoenergetic spectral line. Therefore, detecting such a signature would be strong evidence for dark matter annihilation or decay. We search for spectral lines in the Fermi Large Area Telescope observations of the Milky Way halo in the energy range 200 MeV to 500 GeV using analysis methods from our most recent line searches. The main improvements relative to previous works are our use of 5.8 years of data reprocessed with the Pass 8 event-level analysis and the additional data resulting from the modified observing strategy designed to increase exposure of the Galactic center region. We searched in five sky regions selected to optimize sensitivity to different theoretically-motivated dark matter scenarios and find no significant detections. In addition to presenting the results from our search for lines, we also investigate the previously reported tentative detection of a line at 133 GeV using the new Pass 8 data.
We report on the determination of the lithium abundance [A(Li)] of 52 solar-like stars. For 41 objects the A(Li) here presented corresponds to the first measurement. We have measured the equivalent widths of the 6708\AA\ lithium feature in high-resolution spectroscopic images ($R \sim 80\,000$), obtained at the Observatorio Astrof\'isico Guillermo Haro (Sonora, Mexico), as part of the first scientific observations of the revitalized Lunar and Planetary Laboratory (LPL) Echelle Spectrograph, now known as the Cananea High-resolution Spectrograph (CanHiS). Lithium abundances were derived with the Fortran code MOOG, using as fundamental input a set of atmospheric parameters recently obtained by our group. With the help of an additional small sample with previous A(Li) determinations, we demonstrate that our lithium abundances are in agreement, to within uncertainties, with other works. Two target objects stand out from the rest of the sample. The star BD+47 3218 ($T_{\rm eff}$ = 6050$\pm$52 K, A(Li) = 1.86$\pm$ 0.07 dex) lies inside the so-called lithium desert in the the A(Li)--$T_{\rm eff}$ plane. The other object, BD+28 4515, has an A(Li) = 3.05$\pm$0.07 dex, which is the highest of our sample and compatible with the expected abundances of relatively young stars.
We present ultraviolet (UV) follow-up of a sample of potential dwarf galaxy candidates selected for their neutral hydrogen (HI) properties, taking advantage of the low UV background seen by the GALEX satellite and its large and publicly available imaging footprint. The HI clouds, which are drawn from published GALFA-HI and ALFALFA HI survey compact cloud catalogs, are selected to be galaxy candidates based on their spatial compactness and non-association with known high-velocity cloud complexes or Galactic HI emission. Based on a comparison of their UV characteristics to those of known dwarf galaxies, half (48%) of the compact HI clouds have at least one potential stellar counterpart with UV properties similar to those of nearby dwarf galaxies. If galaxies, the star formation rates, HI masses, and star formation efficiencies of these systems follow the trends seen for much larger galaxies. The presence of UV emission is an efficient method to identify the best targets for spectroscopic follow-up, which is necessary to prove that the stars are associated with the compact HI. Further, searches of this nature help to refine the salient HI properties of likely dwarfs (even beyond the Local Group). In particular, HI compact clouds considered to be velocity outliers relative to their neighbor HI clouds have the most significant detection rate of single, appropriate UV counterparts. Correcting for the sky coverage of the two all-Arecibo sky surveys yielding the compact HI clouds, these results may imply the presence of potentially hundreds of new tiny galaxies across the entire sky.
The irregular satellites of the giant planets are believed to have been captured during the evolution of the solar system. Knowing their physical parameters, such as size, density, and albedo is important for constraining where they came from and how they were captured. The best way to obtain these parameters are observations in situ by spacecrafts or from stellar occultations by the objects. Both techniques demand that the orbits are well known. We aimed to obtain good astrometric positions of irregular satellites to improve their orbits and ephemeris. We identified and reduced observations of several irregular satellites from three databases containing more than 8000 images obtained between 1992 and 2014 at three sites (Observat\'orio do Pico dos Dias, Observatoire de Haute-Provence, and European Southern Observatory - La Silla). We used the software PRAIA (Platform for Reduction of Astronomical Images Automatically) to make the astrometric reduction of the CCD frames. The UCAC4 catalog represented the International Celestial Reference System in the reductions. Identification of the satellites in the frames was done through their ephemerides as determined from the SPICE/NAIF kernels. Some procedures were followed to overcome missing or incomplete information (coordinates, date), mostly for the older images. We managed to obtain more than 6000 positions for 18 irregular satellites: 12 of Jupiter, 4 of Saturn, 1 of Uranus (Sycorax), and 1 of Neptune (Nereid). For some satellites the number of obtained positions is more than 50\% of what was used in earlier orbital numerical integrations. Comparison of our positions with recent JPL ephemeris suggests there are systematic errors in the orbits for some of the irregular satellites. The most evident case was an error in the inclination of Carme.
Cosmic-ray (CR) electrons and nuclei interact with the Galactic interstellar gas and produce high-energy $\gamma$ rays. The $\gamma$-ray emission rate per hydrogen atom, called emissivity, provides a unique indirect probe of the CR flux. We present the measurement and the interpretation of the emissivity in the solar neighborhood for $\gamma$-ray energy from 50~MeV to 50~GeV. We analyzed a subset of 4 years of observations from the Large Area Telescope (LAT) aboard the {\it Fermi Gamma-ray Space Telescope} ({\it Fermi}) restricted to absolute latitudes $10^o<|b| <70^o$. From a fit to the LAT data including atomic, molecular and ionized hydrogen column density templates as well as a dust optical depth map we derived the emissivities, the molecular hydrogen to CO conversion factor $X_{CO}=(0.902\pm0.007) \times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ and the dust-to-gas ratio $X_{DUST}=(41.4\pm0.3) \times 10^{20}$ cm$^{-2}$ mag$^{-1}$. Moreover we detected for the first time $\gamma$-ray emission from ionized hydrogen. We compared the extracted emissivities to those calculated from $\gamma$-ray production cross-sections and to CR spectra measured in the heliosphere. We observed that the experimental emissivities are reproduced only if the solar modulation is accounted for. This provides a direct detection of solar modulation observed previously through the anticorrelation between CR fluxes and solar activity. Finally we fitted a parametrized spectral form to the heliospheric CR observations as well as to the {\it Fermi}-LAT emissivity and obtained compatible local interstellar spectra for proton and Helium kinetic energy per nucleon between between 1 and 100~GeV and for electron-positrons between 0.1 and 100~GeV.
We present the particle-in-cell (PIC) simulation results of the interaction of a high-energy lepton plasma flow with background electron-proton plasma and focus on the acceleration processes of the protons. It is found that the acceleration follows a two-stage processes. In the first stage, protons are accelerated transversely (perpendicular to the lepton flow) by the turbulent magnetic field "islands" generated via the strong Weibel-type instabilities. The accelerated protons shows a perfect inverse-power energy spectrum. As the interaction continues, a shockwave structure forms and the protons in front of the shockwave are reflected at twice of the shock speed, resulting in a quasi-monoenergetic peak located near 200MeV under the simulation parameters.
We report a novel radio autocorrelation (AC) search for extraterrestrial intelligence (SETI). For selected frequencies across the terrestrial microwave window (1-10 GHz) observations were conducted at the Allen Telescope Array to identify artificial non-sinusoidal periodic signals with radio bandwidths greater than 1 kHz, which are capable of carrying substantial messages with symbol-rates from 10-10e6 Hz. Out of 243 observations, about half (101) were directed toward sources with known continuum flux greater than 1 Jy (quasars, pulsars, supernova remnants and masers), based on the hypothesis that they might harbor heretofore undiscovered natural or artificial, repetitive, phase or frequency modulation. The rest of the targets were mostly toward exoplanet stars and similarly interesting targets from the standpoint of SETI. This campaign rules out several previously untested hypotheses relating to the number of artificially modulated "natural" sources. Since we are using a phase sensitive detector, these observations break new ground on this topic. We conclude that the maximum probability that future observations like the ones described here will reveal repetitively modulated emissions from a wide variety of sources, including quasars, supernova remnants and bright stars, is no more than 15-30 percent, depending on source type. The paper concludes by describing an approach to expanding this survey to many more targets and much greater sensitivity using archived and real-time data from interferometers all over the world.
The primary science driver for 3D galaxy surveys is their potential to constrain cosmological parameters. Forecasts of these surveys' effectiveness typically assume Gaussian statistics for the underlying matter density, despite the fact that the actual distribution is decidedly non-Gaussian. To quantify the effect of this assumption, we employ an analytic expression for the power spectrum covariance matrix to calculate the Fisher information for BAO-type model surveys. We find that for typical number densities, at $k_\mathrm{max} = 0.5 h$ Mpc$^{-1}$, Gaussian assumptions significantly overestimate the information on all parameters considered, in some cases by up to an order of magnitude. However, after marginalizing over a six-parameter set, the form of the covariance matrix (dictated by $N$-body simulations) causes the majority of the effect to shift to the "amplitude-like" parameters, leaving the others virtually unaffected. We find that Gaussian assumptions at such wavenumbers can underestimate the dark energy parameter errors by well over 50 per cent, producing dark energy figures of merit almost 3 times too large. Thus, for 3D galaxy surveys probing the non-linear regime, proper consideration of non-Gaussian effects is essential.
In this work we present a procedure to infer the mass of progenitors and remnants of Gamma Ray Bursts (GRB), starting from the observed energy $E_{iso}^{GRB}$ emitted isotropically and considering the associated emission of Gravitational Waves (GW) $ E_{iso}^{GW}$ in the different phases. We assume that the GW energy of the progenitor $E_{PROG}^{GW}$ is emitted partially during a star collapse, and the residual energy is related to the GW energy emitted by the remnant. We take a sample of $237$ Long GRB, and use an hybrid Montecarlo procedure to explore, for each of them, a region of possible solutions of $ E_{iso}^{GW}$ as a function of the masses, radii, oblateness, rotation frequencies of progenitor and remnant and the fraction of energy $k$ emitted as GW by the GRB. We discriminate between a Neutron Star (NS) or Black Hole (BH) for the remnant and obtain interesting values for the GW emitted by the remnant NS or BH, for the conversion factor $k$ of and for the masses and radii of GRB progenitor stars. We also obtain remnant populations with mean masses, mean GW frequencies and GRB frequency of GW emission in agreement with the most accepted models.
We investigate the multiple stellar populations in one of peculiar globular clusters (GCs) M22 using new ground-based wide-field Ca by and HST/WFC3 photometry with equivalent passbands, confirming our previous result that M22 has a distinctive red-giant branch (RGB) split mainly due to difference in metal abundances. We also make use of radial velocity measurements of the large number of cluster membership stars by other. Our main results are followings. (i) The RGB and the subgiant branch (SGB) number ratios show that the calcium-weak (Ca-w) group is the dominant population of the cluster. However, an irreconcilable difference can be seen in the rather simple two horizontal branch (HB) classification by other. (ii) Each group has its own CN-CH anti-correlation. However, the alleged CN-CH positive correlation is likely illusive. (iii) The location of the RGB bump of the calcium-strong (Ca-s) group is significantly fainter, which may pose a challenge to the helium enhancement scenario in the Ca-s group. (iv) The positions of the center are similar. (v) The Ca-w group is slightly more centrally concentrated, while the the Ca-s is more elongated at larger radii. (vi) The mean radial velocities for both groups are similar, but the Ca-s group has a larger velocity dispersion. (vii) The Ca-s group rotates faster. The plausible scenario for the formation of M22 is that it formed via a merger of two GCs in a dwarf galaxy environment and accreted later to our Galaxy.
The cosmic distance duality relation is a milestone of cosmology involving the luminosity and angular diameter distances. Any departure of the relation points to new physics or systematic errors in the observations, therefore tests of the relation are extremely important to build a consistent cosmological framework. Here, two new tests are proposed based on galaxy clusters observations (angular diameter distance and gas mass fraction) and $H(z)$ measurements. By applying Gaussian Processes, a non-parametric method, we are able to derive constraints on departures of the relation where no evidence of deviation is found in both methods, reinforcing the cosmological and astrophysical hypotheses adopted so far.
Large-area surveys operating at mid-infrared wavelengths have proven to be a valuable means of discovering and characterizing minor planets. Through the use of radiometric models, it is possible to derive physical properties such as diameters, albedos, and thermal inertia for large numbers of objects. Modern detector array technology has resulted in a significant improvement in spatial resolution and sensitivity compared with previous generations of space-based infrared telescopes, giving rise to a commensurate increase in the number of objects that have been observed at these wavelengths. Space-based infrared surveys of asteroids therefore offer an effective means of rapidly gathering information about small body populations' orbital and physical properties. The AKARI, WISE/NEOWISE, Spitzer, and Herschel missions have significantly increased the number of minor planets with well-determined diameters and albedos.
The X-ray and mid-IR emission from active galactic nuclei (AGN) are strongly correlated. However, while various published parameterizations of this correlation are consistent with the low-redshift, local Seyfert galaxy population, extrapolations of these relations to high luminosity differ by an order of magnitude at nuL(nu)(6um) = 1e47 erg/s. Using data from the Wide-field Infrared Survey Explorer, we determine the mid-IR luminosities of the most luminous quasars from the Sloan Digital Sky Survey and present a revised formulation of the X-ray to mid-IR relation of AGN which is appropriate from the Seyfert regime to the powerful quasar regime.
The Tully$-$Fisher relation (TFR) expresses the connection between rotating galaxies and the dark matter haloes they inhabit, and therefore contains a wealth of information about galaxy formation. We construct a general framework to investigate whether models based on halo abundance matching are able to reproduce the observed stellar mass TFR and mass$-$size relation (MSR), and use the data to constrain galaxy formation parameters. Our model tests a range of plausible scenarios, differing in the response of haloes to disc formation, the relative angular momentum of baryons and dark matter, the impact of selection effects, and the abundance matching parameters. We show that agreement with the observed TFR puts an upper limit on the scatter between galaxy and halo properties, requires weak or reversed halo contraction, and favours selection effects that preferentially eliminate fast-rotating galaxies. The MSR constrains the ratio of the disc to halo specific angular momentum to be approximately in the range 0.6$-$1.2. We identify and quantify two problems that models of this nature face: 1) They predict too large an intrinsic scatter for the MSR, and 2) they predict too strong an anticorrelation between the TFR and MSR residuals. We argue that resolving these problems requires introducing a correlation between stellar surface density and enclosed dark matter mass, which could be achieved by correlating abundance matching scatter with halo spin, or by performing a second abundance match between galaxy size and halo mass or concentration. Finally, we explore the expected difference between the TFRs of central and satellite galaxies, and find that in the favoured models this difference should be detectable in a sample of $\sim$700 galaxies.
PolyChord is a novel nested sampling algorithm tailored for high-dimensional parameter spaces. This paper coincides with the release of PolyChord v1.3, and provides an extensive account of the algorithm. PolyChord utilises slice sampling at each iteration to sample within the hard likelihood constraint of nested sampling. It can identify and evolve separate modes of a posterior semi-independently, and is parallelised using openMPI. It is capable of exploiting a hierarchy of parameter speeds such as those present in CosmoMC and CAMB, and is now in use in the CosmoChord and ModeChord codes. PolyChord is available for download at: this http URL
A computed line list for formaldehyde, H$_2{}^{12}$C$^{16}$O, applicable to temperatures up to $T=1500$~K is presented. An empirical potential energy and {\it ab initio} dipole moment surfaces are used as the input to nuclear motion program TROVE. The resulting line list, referred to as \textit{AYTY}, contains 10.3 million rotational-vibrational states and around 10 billion transition frequencies. Each transition includes associated Einstein-$A$ coefficients and absolute transition intensities, for wavenumbers below 10~000 cm\(^{-1}\) and rotational excitations up to \(J=70\). Room-temperature spectra are compared with laboratory measurements and data currently available in the HITRAN database. These spectra show excellent agreement with experimental spectra and highlight the gaps and limitations of the HITRAN data. The full line list is available from the CDS database as well as at \url{www.exomol.com}.
Accurate and complete rotational, rotational-vibrational and rotational-vibrational-electronic line lists are calculated for sodium hydride: both the NaH and NaD isotopologues are considered. These line lists cover all ro-vibrational states of the ground ($X$~$^1\Sigma^+$) and first excited ($A$~$^1\Sigma^+$) electronic states. The calculations use available spectroscopically-determined potential energy curves and new high-quality, \textit{ab initio} dipole moment curves. Partition functions for both isotopologues are calculated and the effect of quasibound states is considered. The resulting line lists are suitable for temperatures up to about 7000~K and are designed for studies of exoplanet atmospheres, brown dwarfs and cool stars. In particular, the NaH $A-X$ band is found to show a broad absorption feature at about 385 nm which should provide a signature for the molecule. All partition functions, lines and transitions are available as Supplementary Information to this article and at \url{www.exomol.com}.
There are very few confirmed black holes with a mass that could be $\sim\! 4\, M_\odot$ and no neutron stars with masses greater than $\sim\! 2\, M_\odot$, creating a gap in the observed distribution of compact star masses. Some black holes with masses between 2 and $4\, M_\odot$ might be hiding among other X-ray sources, whose masses are difficult to measure. We present new high-speed optical photometry of the low-mass X-ray binary V1408 Aql (= 4U 1957+115), which is a persistent X-ray source thought to contain a black hole. The optical light curve of V1408~Aql shows a nearly sinusoidal modulation at the orbital period of the system superimposed on large night-to-night variations in mean intensity. We combined the new photometry with previously-published photometry to derive a more precise orbital period, $P = 0.388893(3)$\ d, and to better define the orbital light curve and night-to-night variations. The orbital light curve agrees well with a model in which the modulation is caused entirely by the changing aspect of the heated face of the secondary star. The lack of eclipses rules out orbital inclinations greater than $65^{\circ}$. Our best models for the orbital light curve favor inclinations near $13^{\circ}$ and black hole masses near $3\, M_\odot$ with a 90\% upper bound of $6.2\, M_\odot$, and a lower bound of $2.0\, M_\odot$ imposed solely by the maximum mass of neutron stars. We favor a black hole primary over a neutron star primary based on evidence from the X-ray spectra, the high spin of the compact object, and the fact that a type I X-ray burst has not been observed for this system. Although uncertainties in the data and the models allow higher masses, possibly much higher masses, the compact star in V1408~Aql is a viable candidate for a black hole lying in the mass gap.
We present timing solutions and spin properties of the young pulsar PSR B0540-69 from analysis of 15.8 yr of data from the Rossi X-Ray Timing Explorer. We perform a partially phase-coherent timing analysis in order to mitigate the pronounced effects of timing noise in this pulsar. We also perform fully coherent timing over large subsets of the data set in order to arrive at a more precise solution. In addition to the previously reported first glitch undergone by this pulsar, we find a second glitch, which occurred at MJD 52927 $\pm$ 4, with fractional changes in spin frequency $\Delta\nu/\nu = (1.64 \pm 0.05) \times 10^{-9}$ and spin-down rate $\Delta\dot{\nu}/\dot{\nu} = (0.930 \pm 0.011) \times 10^{-4}$ (taken from our fully coherent analysis). We measure a braking index that is consistent over the entire data span, with a mean value $n = 2.129 \pm 0.012$, from our partially coherent timing analysis. We also investigated the emission behavior of this pulsar, and have found no evidence for significant flux changes, flares, burst-type activity, or pulse profile shape variations. While there is strong evidence for the much-touted similarity of PSR B0540-69 to the Crab pulsar, they nevertheless differ in several aspects, including glitch activity, where PSR B0540-69 can be said to resemble other very young pulsars. It seems clear that the specific processes governing the formation, evolution, and interiors of this population of recently born neutron stars can vary significantly, as reflected in their observed properties.
Rotation curves of galaxies show a wide range of shapes, which can be paramaterized as scatter in Vrot(1kpc)/Vmax i.e.the ratio of the rotation velocity measured at 1kpc and the maximum measured rotation velocity. We examine whether the observed scatter can be accounted for by combining scatters in disc scale-lengths, the concentration-halo mass relation, and the Mstar-Mhalo relation. We use these scatters to create model galaxy populations; when housed within dark matter halos that have universal, NFW density profiles, the model does not match the lowest observed values of Vrot(1kpc)/Vmax and has too little scatter in Vrot(1kpc)/Vmax compared to observations. By contrast, a model using a mass dependant dark matter profile, where the inner slope is determined by the ratio of Mstar/Mhalo, produces galaxies with low values of Vrot(1kpc)/Vmax and a much larger scatter, both in agreement with observation. We conclude that the large observed scatter in Vrot(1kpc)/Vmax favours density profiles that are significantly affected by baryonic processes. Alternative dark matter core formation models are also disfavoured by the large scatter in Vrot(1kpc)/Vmax, implying that baryonic physics would be required to affect density profiles in such models, which defeats a large part of their motivation without ruling them out.
We utilise ground-based, balloon-born and satellite climatology data to
reconstruct site and season-dependent vertical profiles of precipitable water
vapour (PWV). We use these profiles to numerically solve radiative transfer
through the atmosphere, and derive atmospheric brightness temperature ($T_{\rm
atm}$) and optical depth ($\tau$) at the centimetre wavelengths.
We validate the reconstruction by comparing the model column PWV, with
photometric measurements of PWV, performed in the clear sky conditions towards
the Sun. Based on the measurements, we devise a selection criteria to filter
the climatology data to match the PWV levels to the expectations of the clear
sky conditions.
We apply the reconstruction to the location of the Polish 32-metre radio
telescope, and characterise $T_{\rm atm}$ and $\tau$ year-round, at selected
frequencies. We also derive the zenith distance dependence for these
parameters, and discuss shortcomings of using planar, single-layer, and
optically thin atmospheric model approximations, in continuum radio-source
flux-density measurement calibrations.
We obtain PWV-$T_{\rm atm}$ and PWV-$\tau$ scaling relations in the clear sky
conditions, and constrain limits to which the actual $T_{\rm atm}$ and $\tau$
can deviate from those derived solely from the climatological data.
Finally, we suggest a statistical method to detect clear sky that involves
ground-level measurements of relative humidity. The accuracy is tested using a
local climatological data. The method may be useful to constrain cloud cover in
cases when no other (and more robust) climatological data are available.
We study orbit of a pressure-confined cloud in the broad-line region (BLR) of active galactic nuclei (AGNs) when the combined effects of the central gravity and anisotropic radiation pressure and the drag force are considered. Physical properties of the intercloud gas such as its pressure and dynamic viscosity are defined as power-law functions of the radial distance. For a drag force proportional to the relative velocity of a cloud and the background gas, a detailed analysis of the orbits is performed for different values of the input parameters. We also present analytical solutions for a situation where the intercloud pressure is uniform and the viscosity is proportional to the inverse square of the radial distance. Our analytical and numerical solutions demonstrate decay of the orbits because of considering the drag force so that a cloud will eventually fall onto the central region after so called time-of-flight. We found that time-of-flight of a BLR cloud is proportional to the inverse of the dimensionless drag coefficient. We discuss if time-of-flight becomes shorter than the life time of the whole system, then existence of mechanisms for continually forming BLR clouds is needed.
The present investigation relates the orbital radii of regular satellites of Uranus, Jupiter, Neptune, and Saturn to photon energies in the spectra of atomic and molecular hydrogen. To explain these observations a model is developed involving stimulated radiative molecular association (SRMA) reactions among the photons and atoms in the protosatellite disks of the planets. In this model thermal energy is extracted from each disk due to a resonance at radii where there is a match between the temperature in the disk and a photon energy. Matter accumulates at these radii, and satellites and rings are ultimately formed. Orbital radii of satellites of Uranus, Jupiter, and Neptune are related to photon energies ($E_{PM}$ values) in the spectrum of molecular hydrogen. Orbital radii of satellites of Saturn are related to photon energies ($E_{PA}$ values) in the spectrum of atomic hydrogen. The first hint that such relationships exist is found in the linearity of the graphs of orbital radii of uranian satellites vs. orbital radii of jovian satellites, as well as in the graphs of orbital radii of uranian satellites vs. orbital radii of neptunian satellites. An expression is determined which gives the temperature in protosatellite disks where the evolution of each satellite begins. This expression is used to find temperature distributions in the disks, which are found to be similar to distributions calculated by other investigators.
We present the photometric properties of a sample of infrared (IR) bright dust obscured galaxies (DOGs). Combining wide and deep optical images obtained with the Hyper Suprime-Cam (HSC) on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer (WISE), we discovered 48 DOGs with $i - K_\mathrm{s} > 1.2$ and $i - [22] > 7.0$, where $i$, $K_\mathrm{s}$, and [22] represent AB magnitude in the $i$-band, $K_\mathrm{s}$-band, and 22 $\mu$m, respectively, in the GAMA 14hr field ($\sim$ 9 deg$^2$). Among these objects, 31 ($\sim$ 65 %) show power-law spectral energy distributions (SEDs) in the near-IR (NIR) and MIR regime, while the remainder show a NIR bump in their SEDs. Assuming that the redshift distribution for our DOGs sample is Gaussian, with mean and sigma $z$ = 1.99 $\pm$ 0.45, we calculated their total IR luminosity using an empirical relation between 22 $\mu$m luminosity and total IR luminosity. The average value of the total IR luminosity is (3.5 $\pm$ 1.1) $\times$ $10^{13}$ L$_{\odot}$, which classifies them as hyper-luminous infrared galaxies (HyLIRGs). We also derived the total IR luminosity function (LF) and IR luminosity density (LD) for a flux-limited subsample of 18 DOGs with 22 $\mu$m flux greater than 3.0 mJy and with $i$-band magnitude brighter than 24 AB magnitude. The derived space density for this subsample is log $\phi$ = -6.59 $\pm$ 0.11 [Mpc$^{-3}$]. The IR LF for DOGs including data obtained from the literature is well fitted by a double-power law. The derived lower limit for the IR LD for our sample is $\rho_{\mathrm{IR}}$ $\sim$ 3.8 $\times$ 10$^7$ [L$_{\odot}$ Mpc$^{-3}$] and its contributions to the total IR LD, IR LD of all ultra-luminous infrared galaxies (ULIRGs), and that of all DOGs are $>$ 3 %, $>$ 9 %, and $>$ 15 %, respectively.
We numerically study the detailed evolutionary features of the wave-like disturbance and its propagation in the eruption. This work is a follow-up to Wang et al., using significantly upgraded new simulations. We focus on the contribution of the velocity vortices and the fast shock reflection and refraction in the solar corona to the formation of the EUV waves. Following the loss of equilibrium in the coronal magnetic structure, the flux rope exhibits rapid motions and invokes the fast-mode shock forward of the rope, which then produces the type II radio burst. The expansion of the fast shock, which is associated with outward motion, takes place in various directions, and the downward expansion shows the reflection and the refraction as a result of the non-uniform background plasma. The reflected component of the fast shock propagates upward and the refracted component propagates downward. As the refracted component reaches the boundary surface, a weak echo is excited. The Moreton wave is invoked as the fast shock touches the bottom boundary, so the Moreton wave lags the type II burst. A secondary echo occurs in the area where reflection of the fast shock encounters the slow-mode shock, and the nearby magnetic field lines are further distorted because of the interaction between the secondary echo and the velocity vortices. Our results indicate that the EUV wave may arise from various processes that are revealed in the new simulations.
The knee phenomenon of the cosmic ray spectrum, which plays an important role in studying the acceleration mechanism of cosmic rays, is still an unsolved mystery. We try to reconcile the knee spectra measured by ARGO-YBJ and Tibet-III. A simple broken power-law model fails to explain the experimental data. Therefore a modified broken power-law model with non-linear acceleration effects is adopted, which can describe the sharp knee structure. This model predicts that heavy elements dominate at the knee.
High-precision asteroseismic data provides an excellent opportunity to test theories of stellar evolution and new physics, such as the properties of the dark matter (DM) of the Universe. Here we will show that some models of DM lead to changes in the classical scenario of stellar evolution. The accumulation of DM in the core of low-mass stars reduces their central temperatures and inhibits the formation of small convective cores in 1.1-1.3 Ms stars. We review the asteroseismic constraints that have been set to the characteristics of the DM particles, obtained comparing the oscillations of the star Alpha Cen B with modified stellar models. To conclude, we discuss the prospects to use CoRoT and Kepler data on main-sequence stars and red giants to further constrain the nature of DM.
The analysis of the vertical velocity dispersion of disc stars is the most direct astronomical means of estimating the local dark matter density, $\rho_{DM}$. Current estimates based on the mid-plane dynamic density use a local baryonic correction that ignores the non-local effects of spiral structure and significantly underestimates the amount of dynamically relevant gas; the additional gas plus the remaining uncertainties make it practically impossible to measure $\rho_{DM}$ from mid-plane kinematics alone. The sampling of inhomogeneous tracer populations with different scale-heights and scale-lengths results in a systematic increase in the observed dispersion gradients and changes in the nominal density distributions that, if not properly considered, can be misinterpreted as a sign of more dark matter. If the disc gravity is modelled using an infinite disc, the local variation in the vertical gravity due to the globally exponential disc components results in an underestimation of the baryonic contribution by as much as ~40% Given only the assumptions of stationarity, an axially and vertically symmetric disc, doubly exponential tracer and mass-component density profiles, a phenomenologically justified model for the cross-dispersion component $\sigma_{Rz}$, and a realistic model for $g_z$, it is possible to solve the full vertical Jeans equation analytically for the vertical dispersion $\sigma_{z}(z)$ and hence test the robustness of previous attempts at measuring $\rho_{DM}$. When the model parameters are estimated from SEGUE G dwarf star data, it is still not possible to explain the difference in behaviour seen in the simple thick- and thin-disc datasets reported by Buedenbender et al.. Rather than being a fundamental problem with the kinematical model, this effect appears to be a further sign of the difficulty of defining and handling kinematically homogeneous tracer populations.
The semiregular variable Z UMa experienced a particularly deep minimum at the end of March 2014 followed by a three month rise to maximum. The rise and subsequent decline were followed with low resolution (R~1000) spectroscopy, possibly for the first time. The spectral type of Z UMa varied between M7III at minimum and M4III at maximum while the bright hydrogen emission lines reported by observers during the last century now appear relatively weak and transitory.
We describe observations of the apparently empty ring galaxy ESO 474-G040 obtained with the Southern African Large Telescope (SALT). The observations, consisting of imaging, long-slit spectroscopy and Fabry-Perot mapping of the H-alpha line, allow determining the ring kinematics as well as estimating the metallicity of the ring and the stellar population composition in its various parts. We propose that the object could best be understood as being the result of a past merger of disk galaxies, which formed a gas ring that subsequently disrupted via the bead instability and is presently forming stars.
Due to the renewed interest in dark matter after the upgrade of the large hadron collider and its dedication to dark matter research it is timely to reassess the whole problem. Dark matter is one way to reconcile the discrepancy between the velocity of matter in the outer regions of galaxies and the observed galactic mass. So far, no credible candidate for dark matter has been identified. Here, we develop a model accounting for observations by rotations and interactions between rotating objects analogous to magnetic fields and interactions with moving charges. The magnitude of these fields is described by a fundamental constant of the order 10-41kg-1. The same interactions can be observed in the solar system where they lead to small changes of planetary orbits.
We model the dust evolution in protoplanetary disks with full 3D, Smoothed Particle Hydrodynamics (SPH), two-phase (gas+dust) hydrodynamical simulations. The gas+dust dynamics, where aerodynamic drag leads to the vertical settling and radial migration of grains, is consistently treated. In a previous work, we characterized the spatial distribution of non-growing dust grains of different sizes in a disk containing a gap-opening planet and investigated the gap's detectability with the Atacama Large Millimeter/submillimeter Array (ALMA). Here we take into account the effects of grain growth and fragmentation and study their impact on the distribution of solids in the disk. We show that rapid grain growth in the two accumulation zones around planet gaps is strongly affected by fragmentation. We discuss the consequences for ALMA observations.
We have carried out multiwavelength observations of the near-by ($z=0.046$) rich, merging galaxy cluster Abell 3376 with the Murchison Widefield Array (MWA). As a part of the GaLactic and Extragalactic All-sky MWA survey (GLEAM), this cluster was observed at 88, 118, 154, 188 and 215 MHz. The known radio relics, towards the eastern and western peripheries of the cluster, were detected at all the frequencies. The relics, with a linear extent of $\sim$ 1 Mpc each, are separated by $\sim$ 2 Mpc. Combining the current observations with those in the literature, we have obtained the spectra of these relics over the frequency range 80 -- 1400 MHz. The spectra follow power laws, with $\alpha$ = $-1.17\pm0.06$ and $-1.37\pm0.08$ for the west and east relics, respectively ($S \propto \nu^{\alpha}$). Assuming the break frequency to be near the lower end of the spectrum we estimate the age of the relics to be $\sim$ 0.4 Gyr. No diffuse radio emission from the central regions of the cluster (halo) was detected. The upper limit on the radio power of any possible halo that might be present in the cluster is a factor of 35 lower than that expected from the radio power and X-ray luminosity correlation for cluster halos. From this we conclude that the cluster halo is very extended ($>$ 500 kpc) and/or most of the radio emission from the halo has decayed. The current limit on the halo radio power is a factor of ten lower than the existing upper limits with possible implications for models of halo formation.
A two dimensional hydrochemical hybrid code, KM2, is constructed to deal with astrophysical problems that would require coupled hydrodynamical and chemical evolution. The code assumes axisymmetry in cylindrical coordinate system, and consists of two modules: a hydrodynamics module and a chemistry module. The hydrodynamics module solves hydrodynamics using a Godunov-type finite volume scheme and treats included chemical species as passively advected scalars. The chemistry module implicitly solves non-equilibrium chemistry and change of the energy due to thermal processes with transfer of external ultraviolet radiation. Self-shielding effects on photodissociation of CO and H$_2$ are included. In this introductory paper, the adopted numerical method is presented, along with code verifications using the hydrodynamics modules, and a benchmark on the chemistry module with reactions specific to a photon-dominated region (PDR). Finally, as an example of the expected capability, the hydrochemical evolution of a PDR is presented based on the PDR benchmark.
Using the anelastic approximation of linearised hydrodynamic equations, we investigate the development of axially symmetric small perturbations in thin Keplerian discs. The sixth-order dispersion equation is derived and numerically solved for different values of relevant physical parameters (viscosity, heat conductivity, disc semi-thickness and vertical structure). The analysis reveals the appearance of two overstable modes which split out from the classical Rayleigh inertial modes in a wide range of the parameters in both ionized and neutral gases. These modes have a viscous-convective nature and can serve as a seed for turbulence in astrophysical discs even in the absence of magnetic fields.
A study of the star KIC5892969 observed by the Kepler satellite is presented. Its three highest amplitude modes present a strong amplitude modulation. The aim of this work is to investigate amplitude variations in this star and their possible cause. Using the 4 years-long observations available, we obtained the frequency content of the full light curve. Then, we studied the amplitude and phase variations with time using shorter time stamps. The results obtained are compared with the predicted ones for resonant mode coupling of an unstable mode with lower frequency stable modes. Our conclusion is that resonant mode coupling is consistent as an amplitude limitation mechanism in several modes of KIC5892969 and we discuss to which extent it might play an important role for other delta Scuti stars.
Statistical isotropy (SI) is one of the fundamental assumptions made in cosmological model building. This assumption is now being rigorously tested using the almost full sky measurements of the CMB anisotropies. A major hurdle in any such analysis is to handle the large biases induced due to the process of masking. We have developed a new method of analysis, using the bipolar spherical harmonic basis functions, in which we semi-analytically evaluate the modifications to SI violation induced by the mask. The method developed here is generic and can be potentially used to search for any arbitrary form of SI violation. We specifically demonstrate the working of this method by recovering the Doppler boost signal from a set of simulated, masked CMB skies.
From a large sample of $\approx 170,000$ local SDSS galaxies we find that the Fundamental Metallicity Relation (FMR) has an overabundance of outliers, compared to what would be expected from a Gaussian distribution of residuals, with significantly lower metallicities than predicted from their stellar mass and star formation rate (SFR). This low-metallicity population has lower stellar masses, bimodial specific SFRs with enhanced star formation within the aperture and smaller half-light radii than the general sample, and is hence a physically distinct population. We show that they are consistent with being galaxies that are merging or have recently merged with a satellite galaxy. In this scenario, low-metallicity gas flows in from large radii, diluting the metallicity of star-forming regions and enhancing the specific SFR until the inflowing gas is processed and the metallicity has recovered. We introduce a simple model in which mergers with a mass ratio larger than a minimum dilute the central galaxy's metallicity by an amount that is proportional to the stellar mass ratio for a constant time, and show that it provides an excellent fit to the distribution of FMR residuals. We find the dilution time-scale to be $\tau=1.568_{-0.027}^{+0.029}$ Gyr, the average metallicity depression caused by a 1:1 merger to be $\alpha=0.2480_{-0.0020}^{+0.0017}$ dex and the minimum mass ratio merger that can be discerned from the intrinsic Gaussian scatter in the FMR to be $\xi_\text{min}=0.2030_{-0.0095}^{+0.0127}$ (these are statistical errors only). From this we derive that the average metallicity depression caused by a merger with mass ratio between 1:5 and 1:1 is 0.114 dex.
We present a new extensive analysis of the old problem of finding a satisfactory calibration of the relation between the geometric albedo and some measurable polarization properties of the asteroids. To achieve our goals, we use all polarimetric data at our disposal. For the purposes of calibration, we use a limited sample of objects for which we can be confident to know the albedo with good accuracy, according to previous investigations of other authors. We find a new set of updated calibration coefficients for the classical slope - albedo relation, but we generalize our analysis and we consider also alternative possibilities, including the use of other polarimetric parameters, one being proposed here for the first time, and the possibility to exclude from best-fit analyzes the asteroids having low albedos. We also consider a possible parabolic fit of the whole set of data.
There is a forbidden region in the parameter space of quasinormal modes of black holes in general relativity. Using both inspiral and ringdown phases of gravitational waves from binary black holes, we propose two methods to test general relativity. We also evaluate how our methods will work when we apply them to Pop III black-hole binaries with typical masses. Adopting simple mean of the estimated range of the event rate, we have the expected rate of 500 ${\rm yr^{-1}}$. Then, the rates of events with SNR $>20$ and SNR $>50$ are 32 ${\rm yr^{-1}}$ and 2 ${\rm yr^{-1}}$, respectively. Therefore, there is a good chance to confirm (or refute) the Einstein theory in the strong gravity region by observing the expected quasinormal modes.
We monitored BL Lacertae for 13 nights in optical B, V, R, and I bands during October and November 2014 including quasi-simultaneous observations in V and R bands using two optical telescopes in India. We have studied multi-band optical flux variations, colour variation and spectral changes in this blazar. Source was found to be active during the whole monitoring period and showed significant intraday variability on 3 nights in V and R filters while displayed hints of variability on 6 other dates in R passband and on 2 nights in V filter. From the colour-magnitude analysis of the source we found that the spectra of the target gets flatter as it becomes brighter on intra-night timescale. Using discrete correlation technique, we found that intraday light curves in both V and R filters are almost consistent and well correlated with each other. We also generated spectral energy distribution (SED) of the target using the B, V, R, and I data sets for all 13 nights which could help us investigate the physical process responsible for the observed variations in BL Lacertae objects. We also discuss possible physical causes of the observed spectral variability.
We present a high-resolution elemental-abundance analysis for a sample of 23 very metal-poor (VMP; [Fe/H] < -2.0) stars, 12 of which are extremely metal-poor (EMP; [Fe/H] < -3.0), and 4 of which are ultra metal-poor (UMP; [Fe/H] < -4.0). These stars were targeted to explore differences in the abundance ratios for elements that constrain the possible astrophysical sites of element production, including Li, C, N, O, the alpha-elements, the iron-peak elements, and a number of neutron-capture elements. This sample substantially increases the number of known carbon-enhanced metal-poor (CEMP) and nitrogen-enhanced metal-poor (NEMP) stars -- our program stars include eight that are considered "normal" metal-poor stars, six CEMP-no stars, five CEMP-s stars, two CEMP-r stars, and two CEMP-r/s stars. One of the CEMP-$r$ stars and one of the CEMP-r/s stars are possible NEMP stars. We detect lithium for three of the six CEMP-no stars, all of which are Li-depleted with respect to the Spite plateau. The majority of the CEMP stars have [C/N] > 0. The stars with [C/N] < 0 suggest a larger degree of mixing; the few CEMP-no stars that exhibit this signature are only found at [Fe/H] < -3.4, a metallicity below which we also find the CEMP-no stars with large enhancements in Na, Mg, and Al. We confirm the existence of two plateaus in the absolute carbon abundances of CEMP stars, as suggested by Spite et al. We also present evidence for a "floor" in the absolute Ba abundances of CEMP-no stars at A(Ba)~ -2.0.
We have undertaken a nearly simultaneous optical/UV and X-ray variability study of the flat spectrum radio quasar, 3C 273 using data available from the XMM$-$Newton satellite mission from June 2000 to July 2012. Here we focus on the multi-wavelength flux variability on both intra-day and long time scales of this very well known radio-loud source. We found high flux variability over long time scales in all bands for which observations were made. The optical/UV variability amplitude was more than twice than that in the X-ray bands. There is some frequency dependence of the variability in optical/UV bands in the sense that the variability amplitude increases with increasing frequency; however, the X-ray emissions disagree with this trend as the variability amplitude decreases from soft to hard X-ray bands. On intra-day time scales 3C 273 showed small amplitude variability in X-ray bands. A hardness ratio analysis in the X-ray regime indicates that the particle acceleration mechanism dominates the cooling mechanism during most of the ~12 year span of these observations.
We develop a method to search for pair halos around active galactic nuclei (AGN) through a temporal analysis of gamma-ray data. The basis of our method is an analysis of the spatial distributions of photons coming from AGN flares and from AGN quiescent states and a further comparison of these two spatial distributions. This method can also be used for a reconstruction of a point spread function (PSF). We found no evidence for a pair halo component through this method by applying it to the Fermi-LAT data in the energy bands of 4.5-6, 6-10, and >10 GeV and set upper limits on the fraction of photons attributable to a pair halo component. An illustration of how to reconstruct the PSF of Fermi-LAT is given.
To search for optical variability on a wide range of timescales, we have carried out photometric monitoring of two flat spectrum radio quasars, 3C 454.3 and 3C 279, plus one BL Lac, S5 0716+714, all of which have been exhibiting remarkably high activity and pronounced variability at all wavelengths. CCD magnitudes in B, V, R and I pass-bands were determined for $\sim$ 7000 new optical observations from 114 nights made during 2011 - 2014, with an average length of $\sim$ 4 h each, at seven optical telescopes: four in Bulgaria, one in Greece, and two in India. We measured multiband optical flux and colour variations on diverse timescales. Discrete correlation functions were computed among B, V, R, and I observations, to search for any time delays. We found weak correlations in some cases with no significant time lags. The structure function method was used to estimate any characteristic time-scales of variability. We also investigated the spectral energy distribution of the three blazars using B, V, R, I, J and K pass-band data. We found that the sources almost always follows a bluer-when-brighter trend. We discuss possible physical causes of the observed spectral variability.
Ongoing studies with XMM-Newton have shown that powerful accretion disc winds, as revealed through highly-ionised Fe\,K-shell absorption at E>=6.7 keV, are present in a significant fraction of Active Galactic Nuclei (AGN) in the local Universe (Tombesi et al. 2010). In Gofford et al. (2013) we analysed a sample of 51 Suzaku-observed AGN and independently detected Fe K absorption in ~40% of the sample, and we measured the properties of the absorbing gas. In this work we build upon these results to consider the properties of the associated wind. On average, the fast winds (v_out>0.01c) are located <r>~10^{15-18} cm (typically ~10^{2-4} r_s) from their black hole, their mass outflow rates are of the order <M_out>~0.01-1 Msun/yr or ~(0.01-1) M_edd and kinetic power is constrained to <L_k> ~10^{43-45} erg/s, equivalent to ~(0.1-10%) L_edd. We find a fundamental correlation between the source bolometric luminosity and the wind velocity, with v_out \propto L_bol^{\alpha} and \alpha=0.4^{+0.3}_{-0.2}$ (90% confidence), which indicates that more luminous AGN tend to harbour faster Fe K winds. The mass outflow rate M_out, kinetic power L_k and momentum flux P_out of the winds are also consequently correlated with L_bol, such that more massive and more energetic winds are present in more luminous AGN. We investigate these properties in the framework of a continuum-driven wind, showing that the observed relationships are broadly consistent with a wind being accelerated by continuum-scattering. We find that, globally, a significant fraction (~85%) of the sample can plausibly exceed the L_k/L_bol~0.5% threshold thought necessary for feedback, while 45% may also exceed the less conservative ~5% of L_bol threshold as well. This suggests that the winds may be energetically significant for AGN--host-galaxy feedback processes.
We consider perfect fluid bodies ("stars") in general relativity that are axisymmetric, asymptotically flat, and that admit a maximal hypersurface. We show that configurations that extremize the total entropy at fixed ADM mass, ADM angular momentum, and total particle number are stationary with circular flow. For such stars, this establishes that thermodynamic equilibrium implies dynamic equilibrium.
We consider generic rotating axially symmetric "dirty" (surrounded by matter) black holes. Near-horizon circular equatorial orbits are examined in two different cases of near-extremal (small surface gravity $\kappa $) and exactly extremal black holes. This has a number of qualitative distinctions. In the first case, it is shown that such orbits can lie as close to the horizon as one wishes on suitably chosen slices of space-time when $\kappa \rightarrow 0$. This generalizes observation of T.\ Jacobson Class. Quantum Grav. 28 187001 (2011) made for the Kerr metric. If a black hole is extremal ($\kappa =0$), circular on-horizon orbits are impossible for massive particles but, in general, are possible in its vicinity. The corresponding black hole parameters determine also the rate with which a fine-tuned particle on the noncircular near-horizon orbit asymptotically approaches the horizon. Properties of orbits under discussion are also related to the Ba% \~{n}ados-Silk-West effect of high energy collisions near black holes. Impossibility of the on-horizon orbits in question is manifestation of kinematic censorship that forbids infinite energies in collisions.
In this paper we investigate the evolution of a Jordan-Brans-Dicke scalar field, $\Phi$, with a power-law potential in the presence of a second scalar field, $\phi$, with an exponential potential, in both the Jordan and the Einstein frames. We prove that in Jordan-Brans-Dicke theory, the de Sitter solution is not a natural attractor. Instead, we show that the attractor in the Jordan frame corresponds to an "intermediate accelerated" solution of the form $a(t)\simeq e^{\alpha_1 t^{p_1}}$ as $t\rightarrow \infty$ where $\alpha_1>0$ and $0<p_1<1$ for a wide range of parameters. Furthermore, when we work in the Einstein frame we get that the attractor is also an "intermediate accelerated" solution of the form $\mathfrak{a}(\mathfrak{t})\simeq e^{\alpha_2 \mathfrak{t}^{p_2}}$ as $\mathfrak{t}\rightarrow \infty$ where $\alpha_2>0$ and $0<p_2<1$ for the same conditions on the parameter space as in the Jordan frame. In the special case of a quadratic potential in the Jordan frame, or for a constant potential in Einstein's frame, the above intermediate solutions are of saddle type. These results were proved using the center manifold theorem, which is not based on linear approximation.
In the context of massive gravity, bi-gravity and multi-gravity non-minimal matter couplings via a specific composite effective metric were investigated recently. Even if these couplings generically reintroduce the Boulware-Deser ghost, this composite metric is unique in the sense that the ghost reemerges only beyond the decoupling limit and the matter quantum loop corrections do not detune the potential interactions. We consider non-minimal {\it derivative} couplings of the composite metric to matter fields for a specific subclass of Horndeski scalar-tensor interactions. We first explore these couplings in the mini-superspace and investigate in which scenario the ghost remains absent. We further study these non-minimal derivative couplings in the decoupling-limit of the theory and show that the equation of motion for the helicity-0 mode remains second order in derivatives. Finally, we discuss preliminary implications for cosmology.
Links to: arXiv, form interface, find, astro-ph, recent, 1506, contact, help (Access key information)
Detailed studies of stellar populations in the halos of the Milky Way and the Andromeda (M 31) galaxies have shown increasing numbers of tidal streams and dwarf galaxies, attesting to a complicated and on-going process of hierarchical structure formation. The most prominent feature in the halo of M 31 is the Giant Stellar Stream, a structure ~4.5 degrees in extent along the sky, which is close to, but not coincident with the galaxy's minor axis. The stars that make up this stream are kinematically and chemically distinct from the other stars in the halo. Here, we present HST/COS high-resolution ultraviolet absorption spectra of three Active Galactic Nuclei sight lines which probe the M 31 halo, including one that samples gas in the main southwestern portion of the Giant Stream. We see two clear absorption components in many metal species at velocities typical of the M 31 halo and a third, blue-shifted component which arises in the stream. Photoionization modeling of the column density ratios in the different components shows gas in an ionization state typical of that seen in other galaxy halo environments and suggests solar to slightly super-solar metallicity, consistent with previous findings from stellar spectroscopy.
GRB 130925A is one of the recent additions to the growing family of ultra-long GRBs (T90$ \gtrsim 1000$ s). While the X-ray emission of ultra-long GRBs have been studied extensively in the past, no comprehensive radio dataset has been obtained so far. We report here the early discovery of an unusual radio afterglow associated with the ultra-long GRB 130925A. The radio emission peaks at low-frequencies ($\sim 7$ GHz) at early times, only $2.2$ days after the burst occurred. More notably, the radio spectrum at frequencies above $10$ GHz exhibits a rather steep cut-off, compared to other long GRB radio afterglows. This cut-off can be explained if the emitting electrons are either mono-energetic or originate from a rather steep, $dN/dE \propto E^{-4}$, power-law energy distribution. An alternative electron acceleration mechanism may be required to produce such an electron energy distribution. Furthermore, the radio spectrum exhibits a secondary underlying and slowly varying component. This may hint that the radio emission we observed is comprised of emission from both a reverse and a forward shock. We discuss our results in comparison with previous works that studied the unusual X-ray spectrum of this event and discuss the implications of our findings on progenitor scenarios.
Multiple stellar populations are a widespread phenomenon among Galactic globular clusters. Even though the origin of the enriched material from which new generations of stars are produced remains unclear, it is likely that self-enrichment will be feasible only in clusters massive enough to retain this enriched material. We searched for multiple populations in the low mass (M~1.4 x 10^4 M_sun) globular cluster E 3, analyzing SOAR/Goodman multi-object spectroscopy centered on the blue CN absorption features of 23 red giant branch stars. We find that the CN abundance does not present the typical bimodal behavior seen in clusters hosting multi stellar populations, but rather a unimodal distribution that indicates the presence of a genuine single stellar population, or a level of enrichment much lower than in clusters that show evidence for two populations from high-resolution spectroscopy. E 3 would be the first bona fide Galactic old globular cluster where no sign of self-enrichment is found.
Variations in the stellar initial mass function (IMF) have been invoked to explain the spectroscopic and dynamical properties of early-type galaxies. However, no observations have yet been able to disentangle the physical driver. We analyse here a sample of 24 early-type galaxies drawn from the CALIFA survey, deriving in a homogeneous way their stellar population and kinematic properties. We find that the local IMF is tightly related to the local metallicity, becoming more bottom-heavy towards metal-rich populations. Our result, combined with the galaxy mass-metallicity relation, naturally explains previous claims of a galaxy mass-IMF relation, derived from non-IFU spectra. If we assume that - within the star formation environment of early-type galaxies - metallicity is the main driver of IMF variations, a significant revision of the interpretation of galaxy evolution observables is necessary.
We compare Halpha emission profiles of 12 dynamically confirmed black holes (BHs) and 2 neutron star X-ray transients (SXTs) in quiescence with those of a sample of 43 Cataclysmic Variables (CVs), also quiescent. The full-width-half maximum (FWHM) of the Halpha line in SXTs is tightly correlated with the velocity semi-amplitude of the donor star K2=0.233(13) FWHM. This new correlation, when combined with orbital periods (i.e. through photometric light curves), opens the possibility to estimate compact object mass functions from single integration, low-resolution spectroscopy. On the other hand, CVs above the period gap are found to follow a flatter correlation, a likely consequence of their larger mass ratios. We also find that the FWHM traces the disc velocity at ~42% R_L1, independently of binary mass ratio. In addition, for a given FWHM, BHs tend to have lower EWs than CVs. This might be explained by the fact that CVs must be seen at higher inclinations to mimic the same projected disc velocities as BH SXTs. For the same reason CVs with FWHM>~1500 km/s are mostly eclipsing while none of our sample BHs are. Further, we show that there is a vacant/unoccupied region for CVs in the FWHM-EW plane defined by FWHM > 2568 ( 1 - (9/EW)**2 )**0.5. Both the FWHM-K2 correlation and the FWHM-EW plane can be exploited, together with photometric light curves, to efficiently discover quiescent BHs in deep Halpha surveys of the Galactic Plane.
The CALIFA team has recently found that the stellar angular momentum and concentration of late-type spiral galaxies are incompatible with those of lenticular galaxies (S0s), concluding that fading alone cannot satisfactorily explain the evolution from spirals into S0s. Here we explore whether major mergers can provide an alternative way to transform spirals into S0s by analysing the spiral-spiral major mergers from the GalMer database that lead to realistic, relaxed S0-like galaxies. We find that the change in stellar angular momentum and concentration can explain the differences in the $\lambda_\mathrm{Re}$--$R_{90}/R_{50}$ plane found by the CALIFA team. Major mergers thus offer a feasible explanation for the transformation of spirals into S0s.
Extremely large surveys with future experiments like Euclid and the SKA will soon allow us to access perturbation modes close to the Hubble scale, with wavenumbers $k \sim {\cal H}$. If a modified gravity theory is responsible for cosmic acceleration, the Hubble scale is a natural regime for deviations from General Relativity (GR) to become manifest. The majority of studies to date have concentrated on the consequences of alternative gravity theories for the subhorizon, quasi-static regime, however. We investigate how modifications to the gravitational field equations affect perturbations around the Hubble scale, and how this translates into deviations of ultra large-scale relativistic observables from their GR behaviour. Adopting a model-independent ethos that relies only on the broad physical properties of gravity theories, we find that the deviations of the observables are small unless modifications to GR are drastic. The angular dependence and redshift evolution of the deviations is highly parameterisation- and survey-dependent, however, and so they are possibly a rich source of modified gravity phenomenology if they can be measured.
This paper uses dynamical invariants to describe the evolution of collisionless systems subject to time-dependent gravitational forces without resorting to maximum-entropy probabilities. We show that collisionless relaxation can be viewed as a special type of diffusion process in the integral-of-motion space. In time-varying potentials with a fixed spatial symmetry the diffusion coefficients are closely related to virial quantities, such as the specific moment of inertia, the virial factor and the mean kinetic and potential energy of microcanonical particle ensembles. The non-equilibrium distribution function (DF) is found by convolving the initial DF with the Green function that solves Einstein's equation for freely diffusing particles. Such a convolution also yields a natural solution to the Fokker-Planck equations in the energy space. Our mathematical formalism can be generalized to potentials with a time-varying symmetry, where diffusion extends over multiple dimensions of the integral-of-motion space. The new probability theory is in many ways analogous to stochastic calculus, with two significant differences: (i) the equations of motion that govern the trajectories of particles are fully deterministic, and (ii) the diffusion coefficients can be derived self-consistently from microcanonical phase-space averages without relying on ergodicity assumptions. For illustration we follow the cold collapse of $N$-body models in a time-dependent logarithmic potential. Comparison between the analytical and numerical results shows excellent agreement in regions where the potential evolution does not depart too strongly from the adiabatic regime.
The modified geodetic brane cosmology (MGBC) is tested with observational data. The MGBC is derived from the geodetic brane gravity action corrected by the extrinsic curvature of the braneworld. The density parameter coming from this additional term produces an accelerated expansion of geometrical origin. Subject to the Supernovae Ia, Observable Hubble parameter, Baryon Acoustic Oscillations and Cosmic Microwave Background probes, the obtained fit provides enough evidence in the sense that the extrinsic curvature effect is able to reproduce the accelerated expansion of the universe without need of invoking dark energy, exotic matter or cosmological constant. Moreover the MGBC is free of the problems present in other braneworld models.
We compare the spatial distribution of stars which form in hydrodynamical simulations to the spatial distribution of the gas, using the $\mathcal{Q}$-parameter. The $\mathcal{Q}$-parameter enables a self-consistent comparison between the stars and gas because it uses a pixelated image of the gas as a distribution of points, in the same way that the stars (sink particles in the simulations) are a distribution of points. We find that, whereas the stars have a substructured, or hierarchical spatial distribution ($\mathcal{Q} \sim 0.4 - 0.7$), the gas is dominated by a smooth, concentrated component and typically has $\mathcal{Q} \sim 0.9$. We also find no statistical difference between the structure of the gas in simulations that form with feedback, and those that form without, despite these two processes producing visually different distributions. These results suggest that the link between the spatial distributions of gas, and the stars which form from them, is non-trivial.
We combine SDSS and WISE photometry for the full SDSS spectroscopic galaxy sample, creating SEDs that cover lambda=0.4-22 micron for an unprecedented large and comprehensive sample of 858,365 present-epoch galaxies. Using MAGPHYS we then model simultaneously and consistently both the attenuated stellar SED and the dust emission at 12 micron and 22 micron, producing robust new calibrations for monochromatic mid-IR star formation rate proxies. These modeling results provide the first mid-IR-based view of the bi-modality in star formation activity among galaxies, exhibiting the sequence of star-forming galaxies (main sequence) with a slope of dlogSFR/dlogM*=0.80 and a scatter of 0.39 dex. We find that these new star-formation rates along the SF main sequence are systematically lower by a factor of 1.4 than those derived from optical spectroscopy. We show that for most present-day galaxies the 0.4-22 micron SED fits can exquisitely predict the fluxes measured by Herschel at much longer wavelengths. Our analysis also illustrates that the majorities of stars in the present-day universe is formed in luminous galaxies (~L*) in and around the green valley of the color-luminosity plane. We make the matched photometry catalog and SED modeling results publicly available.
We present the study of a set of N-body+SPH simulations of a Milky Way-like system produced by the radiative cooling of hot gas embedded in a dark matter halo. The galaxy and its gaseous halo evolve for 10 Gyr in isolation, which allows us to study how internal processes affect the evolution of the system. We show how the morphology, the kinematics and the evolution of the galaxy are affected by the input supernova feedback energy E$_{\rm SN}$, and we compare its properties with those of the Milky Way. Different values of E$_{\rm SN}$ do not significantly affect the star formation history of the system, but the disc of cold gas gets thicker and more turbulent as feedback increases. Our main result is that, for the highest value of E$_{\rm SN}$ considered, the galaxy shows a prominent layer of extra-planar cold (log(T)<4.3) gas extended up to a few kpc above the disc at column densities of $10^{19}$ cm$^{-2}$. The kinematics of this material is in agreement with that inferred for the HI halos of our Galaxy and NGC 891, although its mass is lower. Also, the location, the kinematics and the typical column densities of the hot (5.3<log(T)<5.7) gas are in good agreement with those determined from the O$_{\rm VI}$ absorption systems in the halo of the Milky Way and external galaxies. In contrast with the observations, however, gas at log(T)<5.3 is lacking in the circumgalactic region of our systems.
The velocity function derived from large scale surveys can be compared with the predictions of LCDM cosmology, by matching the measured rotation velocities Vrot of galaxies to the maximum circular velocity of dark matter (DM) halos Vmax. For Vrot<50km/s, a major discrepancy arises between the observed and LCDM velocity functions. However, the manner in which different observational measures of Vrot are associated with Vmax is not straight forward in dwarf galaxies. We instead relate galaxies to DM halos using the empirical baryon- mass to halo-mass relation, and show that different observational measures of Vrot result in very different velocity functions. We show how the W50 velocity function, i.e. using the HI profile line width at 50% of peak HI flux to measure Vrot, can be reconciled with a LCDM cosmology. Our semi-empirical methodology allows us to determine the region of rotation curves that are probed by HI measurements (RHI), and shows that the Vrot of dwarfs are generally measured at a fraction of Rmax, explaining their tendency to have rising rotation curves. We provide fitting formulae for relating RHI and Reff (the effective radius) to the virial radius of DM halos. To continue to use velocity functions as a probe of LCDM cosmology, it is necessary to be precise about how the different measures of rotation velocity are probing the mass of the DM halos, dropping the assumption that any measure of rotational velocity can be equally used as a proxy for Vmax.
We complement a recent ApJ Letter by Luo et al. by comparing the fraction of starburst galaxies which are interacting with the overall fraction of interacting galaxies in the nearby galaxy population (within 40 Mpc). We confirm that in starburst galaxies the fraction of interacting galaxies is enhanced, by a factor of around 2, but crucially we do so by studying a sample of almost 1500 of the nearest galaxies, including many dwarfs and irregulars. We discuss how adjusting the starburst definition influences the final result and conclude that our result is stable. We find significantly lower fractions of interacting galaxies than Luo et al. did from their larger but more distant sample of galaxies, and argue that the difference is most likely due to various biases in the sample selection, with a representative sample of the nearest galaxies, such as the one used here, being the best possible representation of a general picture. Our overall conclusion is that interactions can and do increase the number of starburst galaxies, and that this result is extremely robust. By far most starburst galaxies, however, show no evidence of a present interaction.
The first observations of Saturn's visible-wavelength aurora were made by the Cassini camera. The aurora was observed between 2006 and 2013 in the northern and southern hemispheres. The color of the aurora changes from pink at a few hundred km above the horizon to purple at 1000-1500 km above the horizon. The spectrum observed in 9 filters spanning wavelengths from 250 nm to 1000 nm has a prominent H-alpha line and roughly agrees with laboratory simulated auroras. Auroras in both hemispheres vary dramatically with longitude. Auroras form bright arcs between 70 and 80 degree latitude north and between 65 and 80 degree latitude south, which sometimes spiral around the pole, and sometimes form double arcs. A large 10,000-km-scale longitudinal brightness structure persists for more than 100 hours. This structure rotates approximately together with Saturn. On top of the large steady structure, the auroras brighten suddenly on the timescales of a few minutes. These brightenings repeat with a period of about 1 hour. Smaller, 1000-km-scale structures may move faster or lag behind Saturn's rotation on timescales of tens of minutes. The persistence of nearly-corotating large bright longitudinal structure in the auroral oval seen in two movies spanning 8 and 11 rotations gives an estimate on the period of 10.65 $\pm$0.15 h for 2009 in the northern oval and 10.8$\pm$ 0.1 h for 2012 in the southern oval. The 2009 north aurora period is close to the north branch of Saturn Kilometric Radiation (SKR) detected at that time.
We present an analysis of 123 Gamma-ray bursts (GRBs) with known redshifts possessing an afterglow plateau phase. We reveal that $L_a-T^{*}_a$ correlation between the X-ray luminosity $L_a$ at the end of the plateau phase and the plateau duration, $T^*_a$, in the GRB rest frame has a power law slope different, within more than 2 $\sigma$, from the slope of the prompt $L_{f}-T^{*}_{f}$ correlation between the isotropic pulse peak luminosity, $L_{f}$, and the pulse duration, $T^{*}_{f}$, from the time since the GRB ejection. Analogously, we show differences between the prompt and plateau phases in the energy-duration distributions with the afterglow emitted energy being on average $10\%$ of the prompt emission. Moreover, the distribution of prompt pulse versus afterglow spectral indexes do not show any correlation. In the further analysis we demonstrate that the $L_{peak}-L_a$ distribution, where $L_{peak}$ is the peak luminosity from the start of the burst, is characterized with a considerably higher Spearman correlation coefficient, $\rho=0.79$, than the one involving the averaged prompt luminosity, $L_{prompt}-L_a$, for the same GRB sample, yielding $\rho=0.60$. Since some of this correlation could result from the redshift dependences of the luminosities, namely from their cosmological evolution we use the Efron-Petrosian method to reveal the intrinsic nature of this correlation. We find that a substantial part of the correlation is intrinsic. We apply a partial correlation coefficient to the new de-evolved luminosities showing that the intrinsic correlation exists.
We present a new survey of HCN(1-0) emission, a tracer of dense molecular gas, focused on the little-explored regime of normal star-forming galaxy disks. Combining HCN, CO, and infrared (IR) emission, we investigate the role of dense gas in Star Formation (SF), finding systematic variations in both the apparent dense gas fraction and the apparent SF efficiency (SFE) of dense gas. The latter may be unexpected, given the popularity of gas density threshold models to explain SF scaling relations. We used the IRAM 30-m telescope to observe HCN(1-0) across 29 nearby disk galaxies whose CO(2-1) emission has previously been mapped by the HERACLES survey. Because our observations span a range of galactocentric radii, we are able to investigate the properties of the dense gas as a function of local conditions. We focus on how the IR/CO, HCN/CO, and IR/HCN ratios (observational cognates of the SFE, dense gas fraction, and dense gas SFE) depend on the stellar surface density and the molecular/atomic ratio. The HCN/CO ratio correlates tightly with these two parameters across a range of 2.1 dex and increases in the high surface density parts of galaxies. Simultaneously, the IR/HCN ratio decreases systematically with these same parameters and is ~6-8 times lower near galaxy centers than in the outer regions. For fixed line-mass conversion factors, these results are incompatible with a simple model in which SF depends only on the gas mass above some density threshold. Only a specific set of environment-dependent conversion factors can render our observations compatible with such a model. Whole cloud models, such as the theory of turbulence regulated SF, do a better job of matching our data. We explore one such model in which variations in the Mach number and in the mean density would respectively drive the trends within galaxy disks and the differences between disk and merging galaxies (abridged).
When searching for deviations of statistical isotropy in CMB, a popular strategy is to write the two-point correlation function (2pcf) as the most general function of four spherical angles (i.e., two unit vectors) in the celestial sphere. Then, using a basis of bipolar spherical harmonics, statistical anisotropy will show up if and only if any coefficient of the expansion with non-trivial bipolar momentum is detected -- although this detection will not in general elucidate the origin of the anisotropy. In this work we show that two new sets of four angles and basis functions exist which completely specifies the 2pcf, while, at the same time, offering a clearer geometrical interpretation of the mechanisms generating the signal. Since the coefficients of these expansions are zero if and only if isotropy holds, they act as a simple and geometrically motivated null test of statistical isotropy, with the advantage of allowing cosmic variance to be controlled in a systematic way. We report the results of the application of these null tests to the latest temperature data released by the Planck collaboration.
Estimates of the total particulate mass of the plumes of Enceladus are important to constrain theories of particle formation and transport at the surface and interior of the satellite. We revisit the calculations of Ingersoll and Ewald (2011), who estimated the particulate mass of the Enceladus plumes from strongly forward scattered light in Cassini ISS images. We model the plume as a combination of spherical particles and irregular aggregates resulting from the coagulation of spherical monomers, the latter of which allows for plumes of lower particulate mass. Though a continuum of solutions are permitted by the model, the best fits to the ISS data consist either of low mass plumes composed entirely of small aggregates or high mass plumes composed of large aggregates and spheres. The high mass plumes can be divided into a population of large aggregates with total particulate mass of 116 +/- 12 X 10^3 kg, and a mixed population of spheres and aggregates consisting of a few large monomers that has a total plume particulate mass of 166 +/- 42 X 10^3 kg, consistent with the results of Ingersoll and Ewald (2011). Meanwhile, the low particulate mass aggregate plumes have masses of 25 +/- 4 X 10^3 kg, leading to a solid to vapor mass ratio of 0.07 +/- 0.01 for the plume. If indeed the plumes are made of such aggregates, then a vapor-based origin for the plume particles is possible. The process of aggregate formation by the coagulation of monomers, which depends on the bulk monomer number density inside the plume vents, requires a total plume vent cross sectional area of at most 1800 m^2 to allow for the aggregates to form before the monomers are ejected into space. Differentiation between the high mass and low mass solutions may be possible if forward scattering observations are taken at scattering angles <2 degrees, or else an independent plume particulate mass measurement becomes available.
Massive gravity provides a natural solution for the dark energy problem of cosmology and is also a candidate for resolving the dark matter problem. I demonstrate that, assuming reasonable scaling relations, massive gravity can provide for Milgrom's law of gravity (or "modified Newtonian dynamics") which is known to remove the need for particle dark matter from galactic dynamics. Milgrom's law comes with a characteristic acceleration, Milgrom's constant, which is observationally constrained to $a_0\approx1.1\times10^{-10}$ m/s$^2$. In the derivation presented here, this constant arises naturally from the cosmologically required mass of gravitons like $a_0\propto c\sqrt{\Lambda}\propto c H_0\sqrt{3\Omega_{\Lambda}}$, with $\Lambda$, $H_0$, and $\Omega_{\Lambda}$ being the cosmological constant, the Hubble constant, and the third cosmological parameter, respectively. My derivation suggests that massive gravity could be the mechanism behind both, dark matter and dark energy.
HESS J1841$-$055 is a diffuse unidentified gamma-ray source with the size of $\sim\,1.3^{\circ}\times\,1^{\circ}$. No conclusive counterpart in other wavelengths has so far detected. To search for X-rays responsible for the TeV emission, the Suzaku observations were conducted, which covered a half region of the HESS source. In the soft band (0.5-2.0 keV), we discovered a diffuse emission, Suzaku J1840.2$-$0552, with the size of $\sim10'$. Since its spectrum was fitted by an optically thin thermal plasma model, Suzaku J1840.2$-$0552 is likely to be a supernova remnant. We also discovered an extended source, Suzaku J1840.2$-$0544, in the hard band (2.0-8.0 keV) with an emission line at 6.1 keV. From the spectral feature and large interstellar absorption, this source is likely to be a cluster of galaxies behind the Galactic plane at the red-shift of $\sim$0.09. The other diffuse source spatially overlaps with the SNR candidate G26.6$-$0.2, which shows a non-thermal dominant spectrum. Since no other candidate is found in the hard X-ray band, we infer that these largely extended sources could be possible counterparts of HESS J1841$-$055.
We have reconstructed the sunspot group count, not by comparisons with other reconstructions and correcting those where they were deemed to be deficient, but by a re-assessment of original sources. The resulting series is a pure solar index and does not rely on input from other proxies, e.g. radionuclides, auroral sightings, or geomagnetic records. 'Backboning' the data sets, our chosen method, provides substance and rigidity by using long-time observers as a stiffness character. Solar activity, as defined by the Group Number, appears to reach and sustain for extended intervals of time the same level in each of the last three centuries since 1700 and the past several decades do not seem to have been exceptionally active, contrary to what is often claimed.
A common idea for the origin of the Galactic diffuse X-ray emission, particularly that of the iron lines from neutral and highly ionized atoms, is a superposition of many cataclysmic variables and coronally active binaries. In this scenario, the flux should symmetrically distribute between the east and west on the plane with respect to Sagittarius A$^*$ because the stellar mass distribution determined by infrared observations is nearly symmetric. This symmetry is confirmed for the highly ionized iron line as well as the continuum emission. However, a clear excess of the neutral iron line in the near east of the Galactic center compared to the near-west side is found. The flux distribution of the excess emission well correlates with molecular column density. The X-ray spectrum of the excess emission is described by a power-law continuum plus a 6.4 keV line with the large equivalent width of $\sim1.3$ keV, which is hardly explained by the low-energy electron bombardment scenario. The longitudinal and latitudinal distribution of the excess emission disfavors the X-ray irradiation, neither by Sagittarius A$^*$ nor by nearby X-ray binaries. Then the low-energy proton bombardment is the most probable origin although the high energy density $\sim 80$ eV cm$^{-3}$ in 0.1-1000 MeV is required and there is no conventional proton source in the vicinity.
We present the results of a Monte Carlo technique to calculate the absolute magnitudes (H) and slope parameters (G) of about 240,000 asteroids observed by the Pan-STARRS1 telescope during the first 15 months of its 3-year all-sky survey mission. The system's exquisite photometry with photometric errors < 0.04 mags, and well-defined filter and photometric system, allowed us to derive accurate H and G even with a limited number of observations and restricted range in phase angles. Our Monte Carlo method simulates each asteroid's rotation period, amplitude and color to derive the most-likely H and G, but its major advantage is in estimating realistic statistical+systematic uncertainties and errors on each parameter. The method was confirmed by comparison with the well-established and accurate results for about 500 asteroids provided by Pravec et al. (2012) and then applied to determining H and G for the Pan-STARRS1 asteroids using both the Muinonen et al. (2010) and Bowell et al. (1989) phase functions.
We have compared the stellar parameters, temperature, gravity and
metallicity, between the LAMOST-DR2 and SDSS-DR12/APOGEE database for stars in
common. It is found that the LAMOST database provides a better red-clump
feature than the APOGEE database in the Teff versus logg diagram. With this
advantage, we have separated red clump stars from red giant stars, and attempt
to establish the calibrations between the two datasets for the two groups of
stars respectively.
It shows that there is a good consistency in temperature with a calibration
close to the one-to-one line, and we can establish a satisfied metallicity
calibration of [Fe/H]_APO=1.18[Fe/H]_APO+0.11 with a scatter of 0.08 dex for
both red clump and red giant branch samples. For gravity, there is no any
correlation for red clump stars between the two databases, and scatters around
the calibrations of red giant stars are substantial. We found two main sources
of the scatters of logg for red giant stars. One is a group of stars with
$0.00253*Teff-8.67<logg<2.6 locating at the forbidden region, and the other is
the contaminated red clump stars, which could be picked out from the unmatched
region where stellar metallicity is not consistent with the position in the
Teff versus logg diagram. After excluding stars at the two regions, we have
established two calibrations for red giant stars,
logg_APO=0.000615*Teff_LAMO+0.697*logg_LAMO-2.208$ (sigma=0.150) for [Fe/H]>-1
and logg_APO=0.000874*Teff_LAMO+0.588*logg_LAMO-3.117 (sigma=0.167) for
[Fe/H]<-1. The calibrations are valid for stars with Teff=3800-5400 K and
logg=0-3.8 dex, and are useful in a work aiming to combine the LAMOST and
APOGEE databases in the future study. In addition, we find that an SVM method
based on seismic logg is a good way to greatly improve the accuracy of gravity
for these two regions at least in the LAMOST database.
NGC 247 is a nearby late-type bulgeless spiral galaxy that contains an inactive nucleus. We report a serendipitous discovery of an X-ray flare from the galaxy center with a luminosity up to 2*10^39 erg/s in the 0.3-10 keV band with XMM-Newton. A Chandra observation confirms that the new X-ray source is spatially coincident with the galaxy nucleus. The XMM-Newton data revealed a hard power-law spectrum with a spectral break near 3-4 keV, no pulsations on timescales longer than 150 ms, and a flat power spectrum consistent with Poisson noise from 1 mHz to nearly 10 Hz. Follow-up observations with Swift detected a second flux peak followed by a luminosity drop by factor of almost 20. The spectral and temporal behaviors of the nuclear source are well consistent with the scenario that the flare was due to an outburst of a low-mass X-ray binary that contains a stellar-mass black hole emitting near its Eddington limit at the peak. However, it cannot be ruled out that the sudden brightening in the nucleus was due to accretion onto a possible low-mass nuclear black hole, fed by a tidally disrupted star or a gas cloud; the MAXI observations limit the peak luminosity of the flare to less than ~10^43 erg/s, suggesting that it is either a low mass black hole or an inefficient tidal disruption event (TDE).
A sample of 70 E+A galaxies are selected from 37, 206 galaxies in the second data release (DR2) of Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) according to the criteria for E+A galaxies defined by Goto, and each of these objects is further visually identified. In this sample, most objects are low redshift E+A galaxies with z < 0.25, and locate in the high latitude sky area with magnitude among 14 to 18 mag in g, r and i bands. A stellar population analysis for the whole sample indicates that the E+A galaxies are characterized by both young and old stellar populations (SPs), and the metal-rich SPs have relatively higher contributions than the metal-poor ones. Additionally, a morphological classification for these objects is performed based on the images taken from the Sloan Digital Sky Survey (SDSS).
We know that magnetic fields are pervasive across all scales in the Universe and over all of cosmic time and yet our understanding of many of the properties of magnetic fields is still limited. We do not yet know when, where or how the first magnetic fields in the Universe were formed, nor do we fully understand their role in fundamental processes such as galaxy formation or cosmic ray acceleration or how they influence the evolution of astrophysical objects. The greatest challenge to addressing these issues has been a lack of deep, broad bandwidth polarimetric data over large areas of the sky. The Square Kilometre Array will radically improve this situation via an all-sky polarisation survey that delivers both high quality polarisation imaging in combination with observations of 7-14 million extragalactic rotation measures. Here we summarise how this survey will improve our understanding of a range of astrophysical phenomena on scales from individual Galactic objects to the cosmic web.
Titanium dioxide, TiO$_2$, is a refractory species that could play a crucial role in the dust-condensation sequence around oxygen-rich evolved stars. To date, gas phase TiO$_2$ has been detected only in the complex environment of the red supergiant VY CMa. We aim to constrain the distribution and excitation of TiO$_2$ around VY CMa in order to clarify its role in dust formation. We analyse spectra and channel maps for TiO$_2$ extracted from ALMA science verification data. We detect 15 transitions of TiO$_2$, and spatially resolve the emission for the first time. The maps demonstrate a highly clumpy, anisotropic outflow in which the TiO$_2$ emission likely traces gas exposed to the stellar radiation field. A roughly east-west oriented, accelerating bipolar-like structure is found, of which the blue component runs into and breaks up around a solid continuum component. A distinct tail to the south-west is seen for some transitions, consistent with features seen in the optical and near-infrared. We find that a significant fraction of TiO$_2$ remains in the gas phase outside the dust-formation zone and suggest that this species might play only a minor role in the dust-condensation process around extreme oxygen-rich evolved stars like VY CMa.
Solar energetic particles (SEPs) are one of the extreme space weather phenomena. A huge SEP event increases the radiation dose received by aircrews, who should be warned of such events as early as possible. We developed a warning system for aviation exposure to SEPs. This article describes one component of the system, which calculates the temporal evolution of the SEP intensity and the spectrum immediately outside the terrestrial magnetosphere. To achieve this, we performed numerical simulations of SEP transport in interplanetary space, in which interplanetary SEP transport is described by the focused transport equation. We developed a new simulation code to solve the equation using a set of stochastic differential equations. In the code, the focused transport equation is expressed in a magnetic field line coordinate system, which is a non-orthogonal curvilinear coordinate system. An inverse Gaussian distribution is employed as the injection profile of SEPs at an inner boundary located near the Sun. We applied the simulation to observed SEP events as a validation test. The results show that our simulation can closely reproduce observational data for the temporal evolution of particle intensity. By employing the code, we developed the WArning System for AVIation Exposure to Solar energetic particles (WASAVIES).
Analyzing the available photometry from the Kepler satellite and other databases, we performed detailed light curve modeling of 10 eclipsing binary systems that were found to exhibit a periodic modulation of their orbital periods. All of the selected systems are detached Algol type, with orbital periods from 0.9 to 2.9 days. In total, 9448 times of minimum for these binaries were analyzed in an attempt to identify the period variations caused by the third bodies in these systems. The well-known method of the light-travel time effect was used for the analysis. The orbital periods of the outer bodies were found to be between 1 and 14 years. This hypothesis makes such systems interesting for future prospective detections of these components, despite their low predicted masses. Considering the dynamical interaction between the orbits, the system KIC 3440230 seems to be the most interesting, in which one would expect the detection of some effects (i.e., changing the inclination) even after a few years or decades of observations.
Astrophysical dynamo theories provide various mechanisms for magnetic field amplification inside galaxies, where weak initial fields grow exponentially on various timescales. We investigate the particular role played by stellar feedback mechanisms in creating strong fluid turbulence, allowing for a magnetic dynamo to emerge. Performing magnetohydrodynamic simulations of isolated cooling halos, for both dwarf and Milky Way sized objects, we compare the magnetic field evolution for various initial field topologies and various stellar feedback mechanisms. We find that feedback can indeed drive strong gas turbulence which gives rise to a fast exponential magnetic field growth. Our simulations feature typical properties of Kolmogorov turbulence with a $k ^{-5/3}$ kinetic energy spectrum, as well as the characteristic properties of a small-scale dynamo, with a $k^{3/2}$ magnetic energy spectrum as predicted by Kazantsev dynamo theory. In these feedback-dominated galaxies, stellar feedback provides forcing on large scales close to the halo scale radius, providing thus exponential field growth on all scales within the galaxy. We also investigate simulations with a final quiescent phase by manually turning off the feedback. As turbulence decreases, the galactic fountain settles into a thin, rotationally supported disk. The magnetic field develops a large-scale, well-ordered structure with quadrupole symmetry, irrespective of the initial field topology, which is in good agreement with magnetic field observations of nearby spirals. Our findings suggest that weak initial seed fields were first amplified by a small-scale dynamo during a violent, feedback-dominated early phase in the galaxy formation history, followed by a more quiescent evolution, where the fields have slowly decayed or were maintained via large-scale dynamo action.
We identify 4 unusually bright (H<~25.5) galaxies from HST and Spitzer CANDELS data with probable redshifts z~7-9. These identifications constitute the brightest-known galaxies to date at z>~7.5. As Y-band observations are not available over the full CANDELS program to perform a standard Lyman-break selection of z>7 galaxies, here we employ an alternate strategy using the deep Spitzer/IRAC data. We identify z~7.1-9.1 galaxies by selecting z>~6 galaxies from the HST CANDELS data that show quite red IRAC [3.6]-[4.5] colors, indicating a strong [OIII] line in the 4.5mu band. This selection strategy was validated using a modest sample for which we have deep Y-band coverage. Here we focus on using this criterion to select the brightest z>~7 sources. Applying the IRAC criteria to all HST-selected optical-dropout galaxies over the full ~900 arcmin**2 of the 5 CANDELS fields revealed four unusually bright z~7.1, 7.6, 7.9 and 8.6 candidates. The median [3.6]-[4.5] color of our selected z~7.1-9.1 sample is consistent with rest-frame [OIII]+Hbeta EWs of ~1600A in the [4.5] band. Keck/MOSFIRE spectroscopy has already been reported for one of our selected sources EGS-zs8-1, showing Lyalpha at a redshift of 7.7302$\pm$0.0006. We present similar Keck/MOSFIRE spectroscopy for a second selected galaxy with a probable 4.7sigma Lyalpha-line at a redshift of 7.4770+/-0.0008. Both have H-band magnitudes of ~25 mag and are ~0.5 mag more luminous (M(UV)~-22.0) than any previously discovered z~8 galaxy, with important implications for the UV LF. Our 3 brightest, highest redshift z>~7 galaxies all lie within the CANDELS EGS field, providing a dramatic illustration of the potential impact of field-to-field variance.
We present a theory of large-scale dynamo action in a turbulent flow that has stochastic, zero-mean fluctuations of the $\alpha$ parameter. We extend the Kraichnan-Moffatt model to explore effects of finite memory of $\alpha$ fluctuations, in a spirit similar to that of Sridhar & Singh (2014), hereafter SS14. Using the first-order smoothing approximation, we derive a linear integro-differential equation governing the dynamics of the large-scale magnetic field, which is non-perturbative in the $\alpha$-correlation time $\tau_{\alpha}$. We recover earlier results in the exactly solvable white-noise (WN) limit where the Moffatt drift does not contribute to the dynamo growth/decay. To study finite memory effects, we reduce the integro-differential equation to a partial differential equation by assuming that the $\tau_{\alpha}$ be small but nonzero and the large-scale magnetic field is slowly varying. We derive the dispersion relation and provide explicit expression for the growth rate as a function of four independent parameters. When $\tau_{\alpha}\neq 0$, we find that: (i) in the absence of the Moffatt drift, but with finite Kraichnan diffusivity, only strong $\alpha$-fluctuations can enable a mean-field dynamo (this is qualitatively similar to the WN case); (ii) in the general case when also the Moffatt drift is nonzero, both, weak or strong $\alpha$ fluctuations, can lead to a large-scale dynamo; and (iii) there always exists a wavenumber ($k$) cutoff at some large $k$ beyond which the growth rate turns negative, irrespective of weak or strong $\alpha$ fluctuations. Thus we show that a finite Moffatt drift can always facilitate large-scale dynamo action if sufficiently strong, even in case of weak $\alpha$ fluctuations, and the maximum growth occurs at intermediate wavenumbers ($\sim$ few times $k_{\alpha}$).
Standard inflationary models yield a characteristic signature of a primordial power spectrum with a red tensor and scalar tilt. Nevertheless, Cannone et al recently suggested that, by breaking the assumption of spatial diffeomorphism invariance in the context of the effective field theory of inflation, a blue tensor spectrum can be achieved without violating the Null Energy Condition. In this context, we explore in which cases a blue tensor tilt can be obtained along with a red tilt in the scalar spectrum. Ultimately, we analyze under which conditions this model can reproduce the specific consistency relation of String Gas Cosmology.
We propose a novel method to constrain the spatial extent of dust around galaxies through the measurement of dust temperature. Our method combines the dust emission of galaxies from far-infrared (FIR) image stacking analysis and the quasar reddening due to the dust absorption around galaxies. As a specific application of our method, we use the stacked FIR emission profiles of SDSS photometric galaxies over the IRAS 100$\mu$m map, and the recent measurement of the SDSS galaxy-quasar cross-correlation. If we adopt a single-temperature dust model, the resulting temperature is around 18K, which is consistent with a typical dust temperature for a central part of galaxies. If we assume an additional dust component with much lower temperature, the current data imply the temperature of the galactic dust needs to be higher, 20K to 30K. Since the model of the density and temperature distribution of dust adopted in the current paper is very simple, we cannot draw any strong conclusion at this point. Nevertheless our novel method with the elaborated theoretical model and multi-band measurement of dust will offer an interesting constraint on the statistical nature of galactic dust.
N2 is abundant in Pluto's atmosphere and on its surface, but the supply is depleted by prodigious atmospheric escape. We demonstrate that cometary impacts could not have delivered enough N2 mass to resupply Pluto's atmospheric escape over time; thus Pluto's N2 is likely endogenous, and therefore was either acquired early in its history or created by chemistry inside/on Pluto. We find that cratering could excavate a considerable amount of N2 to resupply the atmosphere against escape if the near-surface N2 reservoir is deep. However, we find that this process likely falls short of that necessary to resupply the atmosphere against escape by at least an order of magnitude. We conclude that either the escape of N2 from Pluto's atmosphere was on average much lower than the predictions for the current epoch, or that internal activity could be necessary to bring N2 to the surface and resupply escape losses. Observations made by the New Horizons spacecraft in mid-2015 will yield further constraints on the provenance and evolution of Pluto's surface and atmospheric N2, and could reveal evidence for past or present internal activity.
In General Relativity, the constraint equation relating metric and density perturbations is inherently nonlinear, leading to an effective non-Gaussianity in the density field on large scales -- even if the primordial metric perturbation is Gaussian. This imprints a relativistic signature in large-scale structure which is potentially observable, for example via a scale-dependent galaxy bias. The effect has been derived and then confirmed by independent calculations, using second-order perturbation theory. Recently, the physical reality of this relativistic effect has been disputed. The counter-argument is based on the claim that a very long wavelength curvature perturbation can be removed by a coordinate transformation. We argue that while this is true locally, the large-scale curvature cannot be removed by local coordinate transformations. The transformation itself contains the long-wavelength modes and thus includes the correlation. We show how the separate universe approach can be used to understand this correlation, confirming the results of perturbation theory.
We investigate the internal and relative motions of the Taurus and Ophiuchus star-forming regions using a sample of young stars with accurately measured radial velocities and proper motions. We find no evidence for expansion or contraction of the Taurus complex, but a clear indication for a global rotation, resulting in velocity gradients, this suggests a common origin, possibly related to that of Gould's Belt.
The solar photospheric abundance of oxygen is still a matter of debate. For about ten years some determinations have favoured a low oxygen abundance which is at variance with the value inferred by helioseismology. Among the oxygen abundance indicators, the forbidden line at 630nm has often been considered the most reliable even though it is blended with a NiI line. In Papers I and Paper II of this series we reported a discrepancy in the oxygen abundance derived from the 630nm and the subordinate [OI] line at 636nm in dwarf stars, including the Sun. Here we analyse several, in part new, solar observations of the the centre-to-limb variation of the spectral region including the blend at 630nm in order to separate the individual contributions of oxygen and nickel. We analyse intensity spectra observed at different limb angles in comparison with line formation computations performed on a CO5BOLD 3D hydrodynamical simulation of the solar atmosphere. The oxygen abundances obtained from the forbidden line at different limb angles are inconsistent if the commonly adopted nickel abundance of 6.25 is assumed in our local thermodynamic equilibrium computations. With a slightly lower nickel abundance, A(Ni)~6.1, we obtain consistent fits indicating an oxygen abundance of A(O)=8.73+/-0.05. At this value the discrepancy with the subordinate oxygen line remains. The derived value of the oxygen abundance supports the notion of a rather low oxygen abundance in the solar hotosphere. However, it is disconcerting that the forbidden oxygen lines at 630 and 636nm give noticeably different results, and that the nickel abundance derived here from the 630nm blend is lower than expected from other nickel lines.
The paper presents an interactive module created through the Wolfram Demonstrations Project that visualizes the Zeeman effect for the small magnetic field strengths present in the interstellar medium. The paper provides an overview of spectral lines and a few examples of strong and weak Zeeman splitting before discussing the module in depth. Student discovery is aided with example situations to investigate using the interactive module, which is targeted at the upper undergraduate or early graduate level. This module (this http URL), which uses free software, can be used in classroom activities or as a means of introducing students to the Wolfram Demonstrations Project as a learning resource.
We present a time series of high resolution spectra of the Type Ia supernova 2014J, which exploded in the nearby galaxy M82. The spectra were obtained with the HEROS echelle spectrograph installed at the 1.2 m TIGRE telescope. We present a series of 33 spectra with a resolution of R = 20, 000, which covers the important bright phases in the evolution of SN 2014J during the period from January 24 to April 1 of 2014. The spectral evolution of SN 2014J is derived empirically. The expansion velocities of the Si II P-Cygni features were measured and show the expected decreasing behaviour, beginning with a high velocity of 14,000 km/s on January 24. The Ca II infrared triplet feature shows a high velocity component with expansion velocities of > 20, 000 km/s during the early evolution apart from the normal component showing similar velocities as Si II. Further broad P-Cygni profiles are exhibited by the principal lines of Ca II, Mg II and Fe II. The TIGRE SN 2014J spectra also resolve several very sharp Na I D doublet absorption components. Our analysis suggests interesting substructures in the interstellar medium of the host galaxy M82, as well as in our Milky Way, confirming other work on this SN. We were able to identify the interstellar absorption of M82 in the lines of Ca II H & K at 3933 and 3968 A as well as K I at 7664 and 7698 A. Furthermore, we confirm several Diffuse Interstellar Bands, at wavelengths of 6196, 6283, 6376, 6379 and 6613 A and give their measured equivalent widths.
In the framework of metric $f(R)$ gravity, we find the dispersion relation for the propagation of tightly wound spiral density waves in the surface of rotating, self-gravitating disks. Also, new Toomre-like stability criteria for differentially rotating disks has been derived for both fluid and stellar disks.
We present a comprehensive analysis of the globular cluster (GC) system of the Local Group dwarf irregular galaxy NGC 6822. Our study is based on homogeneous optical and near-IR photometry, as well as long-slit spectroscopic observations which are used to determine new radial velocities for 6 GCs, two of which had no previous spectroscopic information. We construct optical-near IR colour-colour diagrams and through comparison to simple stellar population models infer that the GCs have old ages consistent with being 9 Gyr or older, while their metallicities are in the range between -1.6 < [Fe/H] < -0.4. We conduct a kinematic analysis of the GC population and find tentative evidence for weak net rotation of the GC system, in the same sense as that exhibited by the underlying spheroid. The most likely amplitude of rotation is ~10 km/s, approximately half the magnitude of the observed velocity dispersion. Finally, we use the GCs to estimate the dynamical mass of NGC 6822 within 11 kpc and we formally find it to be in the range between (3-4)10^9 Msun. This implies an overall mass-to-light ratio in the range of ~ 30-40 and indicates that NGC 6822 is highly dark matter dominated. The mass and the corresponding mass-to-light ratio estimates are affected by various additional systematic effects due to limitations of the data and the model that are not necessary reflected in the formal uncertainties.
The velocity, position, and action variable evolution of a tidal stream drawn out of a star cluster in a triaxial isochrone potential containing a sub-halo population reproduces many of the orbital effects of more general cosmological halos but allows easy calculation of orbital actions. We employ a spherical shell code which we show accurately reproduces the results of a tree gravity code for a collisionless star cluster. Streams from clusters on high eccentricity orbits, $e\gtrsim 0.6$, can spread out so much that the amount of material at high enough surface density to stand out on the sky may be only a few percent of the stream's total mass. Low eccentricity streams remain more spatially coherent, but sub-halos both broaden the stream and displace the centerline with details depending on the orbits allowed within the potential. Overall, the majority of stream particles have changes in their total actions of only 1-2\%, leaving the mean stream relatively undisturbed. A halo with 1\% of the mass in sub-halos typically spreads the velocity distribution about a factor of two wider than would be expected for a smooth halo. Strong density variations, "gaps", along with mean velocity offsets, are clearly detected in low eccentricity streams for even a 0.2\% sub-halo mass fraction. Around one hundred velocity measurements per kiloparsec of stream will enable tests for the presence of a local sub-halo density as small as 0.2-0.5\% of the local mass density, with about 1\% predicted for 30 kiloparsec orbital radii streams.
Inflationary models can correlate small-scale density perturbations with the long-wavelength gravitational waves (GW) in the form of the Tensor-Scalar-Scalar (TSS) bispectrum. This correlation affects the mass-distribution in the Universe and leads to the off-diagonal correlations of the density field modes in the form of the quadrupole anisotropy. Interestingly, this effect survives even after the tensor mode decays when it re-enters the horizon, known as the fossil effect. As a result, the off-diagonal correlation function between different Fourier modes of the density fluctuations can be thought as a way to probe the large-scale GW and the mechanism of inflation behind the fossil effect. Models of single field slow roll inflation generically predict a very small quadrupole anisotropy in TSS while in models of multiple fields inflation this effect can be observable. Therefore this large scale quadrupole anisotropy can be thought as a spectroscopy for different inflationary models. In addition, in models of anisotropic inflation there exists quadrupole anisotropy in curvature perturbation power spectrum. Here we consider TSS in models of anisotropic inflation and show that the shape of quadrupole anisotropy is different than in single field models. In addition in these models the quadrupole anisotropy is projected into the preferred direction and its amplitude is proportional to $g_* N_e$ where $N_e$ is the number of e-folds and $g_*$ is the amplitude of quadrupole anisotropy in curvature perturbation power spectrum. We use this correlation function to estimate the large scale GW as well as the preferred direction and discuss the detectability of the signal in the galaxy surveys like Euclid and 21 cm surveys.
Galactic HI (neutral hydrogen) shells are central to our understanding of the interstellar medium (ISM), which plays a key role in the development and evolution of galaxies, including our own. Several models involving supernovae and stellar winds have contributed to our broad understanding, but a complete, detailed picture remains elusive. To extend existing Galactic shell catalogs, we visually examined the SETHI (Search for Extraterrestrial HI) database to identify shell-like structures. This high-sensitivity 21-cm radio survey covering the Arecibo sky uniquely provides high-resolution data on shells at a wide range of Galactic latitudes. We present basic information (location, radial velocity, angular size, shape) for 74 previously unidentified HI shells. Due to limitations of coverage and data quality, and the biases inherent in search techniques, our catalog is not a complete sample of Galactic shells. We discuss the catalog completeness, and comment on the new shells' relationship with known interstellar structure as warranted. Unlike many previous catalogs, this sample is not biased towards expanding shells. Where possible we also estimate the kinematic distances, physical sizes, expansion velocities, and energies of these shells. Overall, they are relatively large and old, each the result of multiple supernovae. Unlike previous surveys, we do not find that the shells in our sample are preferentially aligned relative to the Galactic plane.
Following the discovery of high-energy (HE; $E>10\,{\rm MeV}$) and very-high-energy (VHE; $E>100\,{\rm GeV}$) $\gamma$-ray emission from the low-frequency-peaked BL~Lac (LBL) object AP Librae, its electromagnetic spectrum is studied over 60 octaves in energy. Contemporaneous data in radio, optical and UV together with the $\gamma$-ray data are used to construct the most precise spectral energy distribution of this source. The data have been found to be modeled with difficulties with single zone homogeneous leptonic synchrotron self-Compton (SSC) radiative scenarios due to the unprecedented width of the high-energy component when compared to the lower-energy component. The two other LBL objects also detected at VHE appear to have similar modeling difficulties. Nevertheless, VHE $\gamma$ rays produced in the extended jet could account for the VHE flux observed by H.E.S.S.
The Fermi GBM Catalog has been recently published. Previous classification analyses of the BATSE, RHESSI, BeppoSAX, and Swift databases found three types of gamma-ray bursts. Now we analyzed the GBM catalog to classify the GRBs. PCA and Multiclustering analysis revealed three groups. Validation of these groups, in terms of the observed variables, shows that one of the groups coincides with the short GRBs. The other two groups split the long class into a bright and dim part, as defined by the peak flux. Additional analysis is needed to determine whether this splitting is only a mathematical byproduct of the analysis or has some real physical meaning.
So far, multiple stellar systems harbor more than 130 extra solar planets. Dynamical simulations show that the outcome of planetary formation process can lead to various planetary architecture (i.e. location, size, mass and water content) when the star system is single or double. In the late phase of planetary formation, when embryo-sized objects dominate the inner region of the system, asteroids are also present and can provide additional material for objects inside the habitable zone (hereafter HZ). In this study, we make a comparison of several binary star systems and their efficiency to move icy asteroids from beyond the snow-line into orbits crossing the HZ. We modeled a belt of 10000 asteroids (remnants from the late phase of planetary formation process) beyond the snow-line. The planetesimals are placed randomly around the primary star and move under the gravitational influence of the two stars and a gas giant. As the planetesimals do not interact with each other, we divided the belt into 100 subrings which were separately integrated. In this statistical study, several double star configurations with a G-type star as primary are investigated. Our results show that small bodies also participate in bearing a non-negligible amount of water to the HZ. The proximity of a companion moving on an eccentric orbit increases the flux of asteroids to the HZ, which could result into a more efficient water transport on a short timescale, causing a heavy bombardment. A comparison with similar simulations in a single star system indicates also a more efficient flux of asteroids in binary stars.
We reanalyze the behavior of Friedmann-Lema\^itre-Robertson-Walker cosmologies in the recently proposed quasidilaton massive-gravity model, and discover that the background dynamics present hitherto unreported features that require unexpected fine-tuning of the additional fundamental parameters of the theory for an observationally consistent background cosmology. We also identify new allowed regions in the parameters space and exclude some of the previously considered ones. The evolution of the mass of gravitational waves reveals non-trivial behavior, exhibiting a mass squared that may be negative in the past, and that presently, while positive, is larger than the square of the Hubble parameter. These properties of the gravity-wave mass have the potential to lead to observational tests of the theory. While quasidilaton massive gravity is known to have issues with stability at short distances, the current analysis is a first step toward the investigation of the more stable extended quasidilaton massive gravity theory, with some expectation that both the fine-tuning of parameters and the interesting behavior of the gravity-wave mass will persist.
HH 30 is a T Tauri star. In this thesis photometric and polarimetric observations through the source are reported. The observations were carried out using the 84cm telescope of the National Observatory of Mexico (OAN-SPM). This thesis also present a model of the brightness of the source using a Monte Carlo code of Watson and Henney (2001) that solves the radiative transfer equation. I modified this code in order to include the polarization. I used the parameters of Wood y Whitney (1998) to calculate the polarimetric variability observed in HH 30.
Resolved submillimeter imaging of transitional disks is increasingly revealing the complexity of disk structure. Here we present the first high-resolution submillimeter image of a recently identified transitional disk around IRAS 04125+2902 in the Taurus star-forming region. We measure an inner disk hole of ~20 AU around IRAS 04125+2902 by simultaneously modeling new 880 micron Submillimeter Array (SMA) data along with an existing spectral energy distribution supplemented by new Discovery Channel Telescope (DCT) photometry. We also constrain the outer radius of the dust disk in IRAS~04125+2902 to ~50-60 AU. Such a small dust disk could be attributed to initial formation conditions, outward truncation by an unseen companion, or dust evolution in the disk. Notably, the dust distribution of IRAS 04125+2902 resembles a narrow ring (delta R ~ 35 AU) composed of large dust grains at the location of the disk wall. Such narrow dust rings are also seen in other transitional disks and may be evidence of dust trapping in pressure bumps, possibly produced by planetary companions. More sensitive submillimeter observations of the gas are necessary to further probe the physical mechanisms at work in shaping the spatial distribution of large dust in this disk. Interestingly, the IRAS 04125+2902 disk is significantly fainter than other transitional disks that have been resolved at submillimeter wavelengths, hinting that more objects with large disk holes may exist at the faint end of the submillimeter luminosity distribution that await detection with more sensitive imaging telescopes.
Ekpyrotic bouncing cosmologies have been proposed as alternatives to inflation. In these scenarios, the universe is smoothed and flattened during a period of slow contraction preceding the bounce while quantum fluctuations generate nearly scale-invariant super-horizon perturbations that seed structure in the post-bounce universe. An analysis by Tolley and Wesley (2007) showed that, for a wide range of ekpyrotic models, generating a scale-invariant spectrum of adiabatic or entropic fluctuations is only possible if the cosmological background is unstable, in which case the scenario is highly sensitive to initial conditions. In this paper, we analyze an important counterexample: a simple action that generates a Gaussian, scale-invariant spectrum of entropic perturbations during ekpyrotic contraction without requiring fine-tuned initial conditions. Based on this example, we discuss some generalizations.
This paper examines the properties of `lattice universes' wherein point masses are arranged in a regular lattice on space-like hypersurfaces; open, flat, and closed universes are considered. The universes are modelled using the Lindquist-Wheeler (LW) approximation scheme, which approximates the space-time in each lattice cell by Schwarz\-schild geometry. Extending Lindquist and Wheeler's work, we derive cosmological scale factors describing the evolution of all three types of universes, and we use these scale factors to show that the universes' dynamics strongly resemble those of Friedmann-Lema\^itre-Robertson-Walker (FLRW) universes. In particular, we use the scale factors to make more salient the resemblance between Clifton and Ferreira's Friedmann-like equations for the LW models and the actual Friedmann equations of FLRW space-times. Cosmological redshifts for such universes are then determined numerically, using a modification of Clifton and Ferreira's approach; the redshifts are found to closely resemble their FLRW counterparts, though with certain differences attributable to the `lumpiness' in the underlying matter content. Most notably, the LW redshifts can differ from their FLRW counterparts by as much as 30\%, even though they increase linearly with FLRW redshifts, and they exhibit a non-zero integrated Sachs-Wolfe effect, something which would not be possible in matter-dominated FLRW universes without a cosmological constant.
We derive a new equation of state (EoS) for neutron stars (NS) from the outer crust to the core based on modern microscopic Brueckner-Hartree-Fock (BHF) calculations using the Argonne $v_{18}$ potential plus three-body forces computed with the Urbana model. To deal with the inhomogeneous structures of matter in the NS crust, we use the recent Barcelona-Catania-Paris-Madrid (BCPM) nuclear energy density functional that is directly based on the same microscopic BHF calculations, and which is able to reproduce the ground-state properties of nuclei along the periodic table. The EoS of the outer crust requires the masses of neutron-rich nuclei, which are obtained through Hartree-Fock-Bogoliubov calculations with the BCPM functional when they are unknown experimentally. To compute the inner crust, Thomas-Fermi calculations in Wigner-Seitz cells are performed with the same functional. Existence of nuclear pasta is predicted in a range of average baryon densities between $\simeq$0.067 fm$^{-3}$ and $\simeq$0.0825 fm$^{-3}$, where the transition to the core takes place. The NS core is computed from the nuclear EoS of the BHF calculation assuming non-exotic constituents (core of $npe\mu$ matter). In each region of the star, we discuss the comparison of the new EoS with previous EoSes for the complete NS structure, in particular, with the Lattimer-Swesty EoS and with the Shen et al. EoS widely used in astrophysical calculations. The new microscopically derived EoS fulfills at the same time a NS maximum mass of 2~$M_\odot$ with a radius of 10 km, and a 1.5~$M_\odot$ NS with a radius of 11.7 km.
According to Rayleigh's criterion, rotating flows are linearly stable when their specific angular momentum increases radially outward. The celebrated magnetorotational instability opens a way to destabilize those flows, as long as the angular velocity is decreasing outward. Using a short-wavelength approximation we demonstrate that even flows with very steep positive shear can be destabilized by azimuthal magnetic fields which are current-free within the fluid. We illustrate the transition of this instability to a rotationally enhanced kink-type instability in case of a homogeneous current in the fluid, and discuss the prospects for observing it in a magnetized Taylor-Couette flow.
We calculate static and spherically symmetric solutions for the Rastall modification of gravity to describe Neutron Stars (NS). The key feature of the Rastall gravity is the non-conservation of the energy-momentum tensor proportionally to the space-time curvature. Using realistic equations of state for the NS interior we place a bound on the non-GR behaviour of the Rastall theory which should be $\lesssim 0.1\%$ level. This work presents the more stringent contraints on the deviations of GR caused by the Rastall proposal.
In a metastable de Sitter space any object has a finite life expectancy beyond which it undergoes vacuum decay. However, by spreading into different parts of the universe which will fall out of causal contact of each other in future, a civilization can increase its collective life expectancy, defined as the average time after which the last settlement disappears due to vacuum decay. We study in detail the collective life expectancy of two comoving objects in de Sitter space as a function of the initial separation, the horizon radius and the vacuum decay rate. We find that even with a modest initial separation, the collective life expectancy can reach a value close to the maximum possible value of 1.5 times that of the individual object if the decay rate is less than 1% of the expansion rate. Our analysis can be generalized to any number of objects, general trajectories not necessarily at rest in the comoving coordinates and general FRW space-time. As part of our analysis we find that in the current state of the universe dominated by matter and cosmological constant, the vacuum decay rate is increasing as a function of time due to accelerated expansion of the volume of the past light cone. Present decay rate is about 3.7 times larger than the average decay rate in the past and the final decay rate in the cosmological constant dominated epoch will be about 56 times larger than the average decay rate in the past. This considerably weakens the lower bound on the half-life of our universe based on its current age.
Inflationary $\alpha$-attractor models in supergravity, which provide excellent fits to the latest observational data, are based on the Poincare disk hyperbolic geometry. We refine these models by constructing Kahler potentials with built-in inflaton shift symmetry and by making a canonical choice of the goldstino Kahler potential. The refined models are stable with respect to all scalar fields at all $\alpha$, no additional stabilization terms are required. The scalar potential V has a nearly Minkowski minimum at small values of the inflaton field $\varphi$, and an infinitely long dS valley of constant depth and width at large $\varphi$. Because of the infinite length of this shift-symmetric valley, the initial value of the inflaton field at the Planck density is expected to be extremely large. We show that the inflaton field $\varphi$ does not change much until all fields lose their energy and fall to the bottom of the dS valley at large $\varphi$. This provides natural initial conditions for inflation driven by the inflaton field slowly rolling along the dS valley towards the minimum of the potential at small $\varphi$. A detailed description of this process is given for $\alpha$-attractors in supergravity, but we believe that our general conclusions concerning naturalness of initial conditions for inflation are valid for a broad class of inflationary models with sufficiently flat potentials.
We construct a new isentropic equation of state (EOS) at finite temperature "CRover" on the basis of the hadron-quark crossover at high density. By using the new EOS, we study the structure of hot neutron stars at birth with the typical lepton fraction ($Y_l=0.3-0.4$) and the typical entropy per baryon ($S=1-2$). Due to the gradual appearance of quark degrees of freedom at high density, the temperature T and the baryon density at the center of the hot neutron stars with the hadron-quark crossover are found to be smaller than those without the crossover by a factor of 2 or more. Typical energy release due to the contraction of a hot neutron star to a cold neutron star with 1.4 solarmass is shown to be about 0.04 solarmass with the spin-up rate about 14%.
We present new traversable wormhole and non-singular black hole solutions in pure, scale-free $R^2$ gravity. These exotic solutions require no null energy condition violating or "exotic" matter and are supported only by the vacuum of the theory. It is well known that $f(R)$ theories of gravity may be recast as dual theories in the Einstein frame. The solutions we present are found when the conformal transformation required to move to the dual frame is singular. In the specific case of purely quadratic $R^2$ gravity, the required conformal factor is identically zero for spacetimes with $R=0$. Solutions in this case are argued to arise in the strong coupling limit of General Relativity.
Links to: arXiv, form interface, find, astro-ph, recent, 1506, contact, help (Access key information)
Dark matter cores within galaxy haloes can be formed by energy feedback from star forming regions: an energy balance suggests that the maximum core formation efficiency arises in galaxies with M$_{\star}\sim10^{8.5}$M$_{\odot}$. We show that a model population of galaxies, in which the density profile has been modified by such baryonic feedback, is able to explain the observed galaxy velocity function and Tully-Fisher relations significantly better than a model in which a universal cuspy density profile is assumed. Alternative models, namely warm or self-interacting dark matter, also provide a better match to these observed relations than a universal profile model does, but make different predictions for how halo density profiles vary with mass compared to the baryonic feedback case. We propose that different core formation mechanisms may be distinguished based on the imprint they leave on galaxy populations over a wide range of mass. Within the current observational data we find evidence of the expected signatures of the mass dependence of core formation generated by baryonic feedback.
We present Keck/DEIMOS spectroscopy of stars in the recently discovered Milky Way satellites Hydra II, Pisces II, and Laevens 1. We measured a velocity dispersion of 5.4 (+3.6 -2.4) km/s for Pisces II, but we did not resolve the velocity dispersions of Hydra II or Laevens 1. We marginally resolved the metallicity dispersions of Hydra II and Pisces II but not Laevens 1. Furthermore, Hydra II and Pisces II obey the luminosity-metallicity relation for Milky Way dwarf galaxies (<[Fe/H]> = -2.02 +/- 0.08 and -2.45 +/- 0.07, respectively), whereas Laevens 1 does not (<[Fe/H]> = -1.68 +/- 0.05). The kinematic and chemical properties suggest that Hydra II and Pisces II are dwarf galaxies, and Laevens 1 is a globular cluster. We determined that two of the previously observed blue stars near the center of Laevens 1 are not members of the cluster. A third blue star has ambiguous membership. If it is a member, we suggest that it could be a Type II Cepheid variable. Hydra II has a radial velocity <v_helio> = 303.1 +/- 1.4 km/s, similar to the leading arm of the Magellanic stream. The mass-to-light ratio for Pisces II is 370 (+310 -240) M_sun/L_sun. It is not among the most dark matter-dominated dwarf galaxies, but it is still worthy of inclusion in the search for gamma rays from dark matter self-annihilation.
In the "Grand Tack" (GT) scenario for the young solar system, Jupiter formed beyond 3.5 AU from the Sun and migrated as close as 1.5 AU until it encountered an orbital resonance with Saturn. Both planets then supposedly migrated outward for several $10^5$ yr, with Jupiter ending up at ~5 AU. The initial conditions of the GT and the timing between Jupiter's migration and the formation of the Galilean satellites remain unexplored. We study the formation of Ganymede and Callisto, both of which consist of ~50% water and rock, respectively, in the GT scenario. We examine why they lack dense atmospheres, while Titan is surrounded by a thick nitrogen envelope. We model an axially symmetric circumplanetary disk (CPD) in hydrostatic equilibrium around Jupiter. The CPD is warmed by viscous heating, Jupiter's luminosity, accretional heating, and the Sun. The position of the water ice line in the CPD, which is crucial for the formation of massive moons, is computed at various solar distances. We assess the loss of Galilean atmospheres due to high-energy radiation from the young Sun. Ganymede and Callisto cannot have accreted their water during Jupiter's supposed GT, because its CPD (if still active) was too warm to host ices and much smaller than Ganymede's contemporary orbit. From a thermal perspective, the Galilean moons might have had significant atmospheres, but these would probably have been eroded during the GT in < $10^5$ yr by solar XUV radiation. Jupiter and the Galilean moons formed beyond 4.5 (+/-0.5) AU and prior to the proposed GT. Thereafter, Jupiter's CPD would have been dry, and delayed accretion of planetesimals should have created water-rich Io and Europa. While Galilean atmospheres would have been lost during the GT, Titan would have formed after Saturn's own tack, because Saturn still accreted substantially for ~$10^6$ yr after its closest solar approach, ending up at about 7 AU.
Using high-quality spectra of the twin stars in the XO-2 binary system, we have detected significant differences in the chemical composition of their photospheres. The differences correlate strongly with the elements' dust condensation temperature. In XO-2N, volatiles are enhanced by about 0.015 dex and refractories are overabundant by up to 0.090 dex. On average, our error bar in relative abundance is 0.012 dex. We present an early metal-depletion scenario in which the formation of the gas giant planets known to exist around these stars is responsible for a 0.015 dex offset in the abundances of all elements while 20 M_Earth of non-detected rocky objects that formed around XO-2S explain the additional refractory-element difference. An alternative explanation involves the late accretion of at least 20 M_Earth of planet-like material by XO-2N, allegedly as a result of the migration of the hot Jupiter detected around that star. Dust cleansing by a nearby hot star as well as age or Galactic birthplace effects can be ruled out as valid explanations for this phenomenon.
The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits.
The Cosmic X-ray Background (CXB) is the total emission from past AGN activity and peaks in the hard X-ray band (30~keV). In this paper, we consider how the energetics of the CXB are dependent on black hole spin. We show that the higher radiative efficiency of rapidly-spinning black holes will over-represent high spin in the X-ray background, meaning that such sources could be a dominant contributor to the CXB. Using a simple bimodal spin distribution, we demonstrate that a population consisting of 15 per cent maximally-spinning black holes can produce 50 per cent of the CXB. The effect is even more pronounced for flux-limited surveys: 7 per cent of sources with maximally-spinning black holes can produce half of the source counts in such surveys. A host of observations are consistent with an over-representation of high-spin black holes in the samples accessible to us, such as obtained from X-ray reverberation studies, detailed modelling of the iron K$\alpha$ line at 6.4~keV, stacked AGN X-ray spectra and calculation of accretion efficiences from optical--to--UV AGN SED modelling. Future NuSTAR and ASTRO-H hard X-ray surveys will provide the best constraints on the role of spin within the AGN population and in producing the CXB spectrum.
We present a detailed analysis of a large-scale galactic outflow in the CGM of a massive (M_h ~ 10^12.5 Msun), star-forming (6.9 Msun/yr), sub-L* (0.5 L_B*) galaxy at z=0.39853 that exhibits a wealth of metal-line absorption in the spectra of the background quasar Q 0122-003 at an impact parameter of 163 kpc. The galaxy inclination angle (i=63 degree) and the azimuthal angle (Phi=73 degree) imply that the QSO sightline is passing through the projected minor-axis of the galaxy. The absorption system shows a multiphase, multicomponent structure with ultra-strong, wide velocity spread OVI (logN = 15.16\pm0.04, V_{90} = 419 km/s) and NV (logN = 14.69\pm0.07, V_{90} = 285 km/s) lines that are extremely rare in the literature. The highly ionized absorption components are well explained as arising in a low density (10^{-4.2} cm^{-3}), diffuse (10 kpc), cool (10^4 K) photoionized gas with a super-solar metallicity ([X/H] > 0.3). From the observed narrowness of the Lyb profile, the non-detection of SIV absorption, and the presence of strong CIV absorption we rule out equilibrium/non-equilibrium collisional ionization models. The low-ionization photoionized gas with a density of 10^{-2.5} cm^{-3} and a metallicity of [X/H] > -1.4 is possibly tracing recycled halo gas. We estimate an outflow mass of ~2x10^{10} Msun, a mass-flow rate of ~54 Msun/yr, a kinetic luminosity of ~9x10^{41} erg/s, and a mass loading factor of ~8 for the outflowing high-ionization gas. These are consistent with the properties of "down-the-barrel" outflows from infrared-luminous starbursts as studied by Rupke et al. Such powerful, large-scale, metal-rich outflows are the primary means of sufficient mechanical and chemical feedback as invoked in theoretical models of galaxy formation and evolution.
We compare the measured angular cross-correlation between the Fermi-LAT gamma-ray sky and catalogues of extra-galactic objects with the expected signal induced by weakly interacting massive particle (WIMP) dark matter (DM). We include a detailed description of the contribution of astrophysical gamma-ray emitters such as blazars, misaligned AGN and star forming galaxies, and perform a global fit to the measured cross-correlation. Five catalogues are considered: SDSS-DR6 quasars, 2MASS galaxies, NVSS radio galaxies, SDSS-DR8 Luminous Red Galaxies and SDSS-DR8 main galaxy sample. To model the cross-correlation signal we use the halo occupation distribution formalism to estimate the number of galaxies of a given catalogue in DM halos and their spatial correlation properties. We discuss uncertainties in the predicted cross-correlation signal arising from the DM clustering and WIMP microscopic properties, which set the DM gamma-ray emission. The use of different catalogues probing objects at different redshifts reduces significantly, though not completely, the degeneracy among the different gamma-ray components. We find that the presence of a significant WIMP DM signal is allowed by the data but not significantly preferred by the fit, although this is mainly due to a degeneracy with the misaligned AGN component. With modest substructure boost, the sensitivity of this method excludes thermal annihilation cross sections at 95% C.L. for WIMP masses up to few tens of GeV. Constraining the low-redshift properties of astrophysical populations with future data will further improve the sensitivity to DM.
[Abridged] Modern analyses of structure formation predict a universe tangled in a cosmic web of dark matter and diffuse baryons. These theories further predict that by the present day, a significant fraction of the baryons will be shock-heated to $T \sim 10^{5}-10^{7}$K yielding a warm-hot intergalactic medium (WHIM), but whose actual existence has eluded a firm observational confirmation. We have designed a novel experiment to search for signatures of the WHIM, by targeting the putative filaments connecting galaxy clusters. Here, we detail the experimental design and report on our first study of a remarkable QSO sightline, that passes within $\Delta d < 3$ Mpc from $7$ independent cluster-pair axes at redshifts $0.1<z<0.5$. We observed this unique QSO using HST/COS, and conducted a survey of diffuse gas. We find excesses of total HI, narrow HI (NLA; Doppler parameters $b<50$ km/s), broad HI (BLA; $b>50$ km/s) and OVI absorption lines within $\Delta v < 1000$ km/s from the cluster-pairs redshifts, corresponding to $\sim 2$, $\sim 2$, $\sim 6$ and $\sim 4$ times their field expectations, respectively. We also report on covering fractions, $f_c$, of gas close to cluster-pairs, and find that the $f_c$ of BLAs are $\sim 4-7$ times higher than the random expectation (at the $\sim 2 \sigma$ c.l.), whereas the $f_c$ of NLAs and OVI are not significantly enhanced. We argue that a larger relative excess of BLAs compared to those of NLAs is a clear signature of the WHIM, even when accounting for turbulence. We also conclude that most of the reported excesses of NLAs and BLAs cannot be attributed to individual galaxy halos but rather to truly intergalactic material. These results suggest an ubiquitous presence of the WHIM in inter-cluster filaments. The HST archive will enable to extend the present analysis to tens of sightlines, eventually leading towards the firm detection of the WHIM.
Semiconvection - mixing that occurs in regions that are stable when considering compositional gradients, but unstable when ignoring them - is shown to have the greatest potential impact on main sequence stars with masses in the range 1.2 - 1.7 solar masses. We present the first stellar evolution calculations using a prescription for semiconvection derived from extrapolation of direct numerical simulations of double-diffusive mixing down to stellar parameters. The dominant mode of semiconvection in stars is layered semiconvection, where the layer height is an adjustable parameter analogous to the mixing length in convection. The rate of mixing across the semiconvective region is sensitively dependent on the layer height. We find that there is a critical layer height that separates weak semiconvective mixing (where evolution is well-approximated by using the Ledoux criterion) from strong semiconvective mixing (where evolution is well-approximated by using the Schwarzschild criterion). This critical layer height is much smaller than the minimum layer height expected from simulations so we predict that for realistic layer heights, the evolution is nearly the same as a model ran with the Schwarzschild criterion. We also investigate the effects of compositional gradient smoothing, finding that it causes convective cores to artificially shrink in the absence of additional mixing beyond the convective boundary. Layered semiconvection with realistic layer heights provides enough such mixing that stars will still evolve as if the Schwarzschild criterion is employed. Finally, we discuss the potential of detecting such semiconvection and its implication on convective core sizes in solar-like oscillators.
The deep, wide-area (900 arcmin**2) near-IR/WFC3/IR + Spitzer/IRAC observations over the CANDELS program represents a significant resource for constraining the bright end of the z~9 and z~10 luminosity functions (LFs) in the UV. We recently reported the discovery of 6 luminous z~9-10 candidates over the GOODS-North+South fields, but extending this search to the full CANDELS program was limited due to the lack of HST-depth 1.05-micron observations in the other 3 CANDELS fields. Here we attempt to significantly realize the potential of CANDELS for z=9-10 science by combining a search over all 5 fields with results from a new HST program (B9-CANDELS) designed to follow up the highest-probability z~9-10 galaxy candidates with observations at 1.05 microns. The targeted z~9-10 candidates are preselected by taking advantage of the full HST, Spitzer/IRAC S-CANDELS observations, and the deepest-available ground-based optical+near-IR observations. With our follow-up program now 91% complete, we identify 4 new high-probability z~9-10 galaxies. This brings our total sample of bright z~9-10 galaxies to 14, including several other new sources from the CANDELS GOODS + ERS fields.Through extensive simulations, we replicate the selection process for our sample (both the preselection and follow-up) and obtain an accurate estimate of the volume density of bright galaxies (M_{UV,AB}<-20) at both z~9 and z~10. The volume density of bright z~9 and z~10 galaxies that we find is 4.5(-1.3)(+2.5)x and 8(-3)(+7)x lower than found at z~8. When compared with the best-fit evolution in the UV luminosities densities from z~8 to z~4 (>0.4L*), the luminosity densities we find at z~9 and z~10 are ~2x lower than the extrapolated trends. We would expect significant additional gains in these results from the on-going 500-orbit BoRG[z910] program and by obtaining additional follow-up observations over the CANDELS-WIDE fields.
We present an integral field study of the internal structure, kinematics and stellar population of the almost edge-on, intermediate luminosity ($L_ {*}$) elliptical galaxy NGC~4697. We build extended 2-dimensional (2D) maps of the stellar kinematics and line-strengths of the galaxy up to $\sim 0.7 $ effective radii (R$_{eff}$) using a mosaic of 8 VIMOS (VIsible Multi-Objects Spectrograph on the VLT) integral-field unit pointings. We find clear evidence for a rotation-supported structure along the major axis from the 2D kinematical maps, confirming the previous classification of this system as a `fast-rotator'. We study the correlations between the third and fourth Gauss-Hermite moments of the line-of-sight velocity distribution (LOSVD) $h_3$ and $h_4$ with the rotation parameter ($V/\sigma$), and compare our findings to hydrodynamical simulations. We find remarkable similarities to predictions from gas-rich mergers. Based on photometry, we perform a bulge/disk decomposition and study the stellar population properties of the two components. The bulge and the disk show different stellar populations, with the stars in the bulge being older (age$_{\rm bulge}=13.5^{+1.4}_{-1.4}$ Gyr, age$_{\rm disk}=10.5^{+1.6}_{-2.0}$Gyr) and more metal-poor ($\mathrm{[M/H]_{bulge}} = -0.17^{+0.12}_{-0.1}$, $\mathrm{[M/H]_{disk}}=-0.03^{+0.02}_{-0.1}$). The evidence of a later-formed, more metal-rich disk embedded in an older, more metal-poor bulge, together with the LOSVD structure, supports a mass assembly scenario dominated by gas-rich minor mergers and possibly with a late gas-rich major merger that left a previously rapidly rotating system unchanged. The bulge and the disk do not show signs of different stellar Initial Mass Function slopes, and both match well with a Milky Way-like IMF.
The density wave theory for the grand-design two-armed spiral pattern in galaxies is successful in explaining several observed features. However, the long-term persistence of this spiral structure is a serious problem since the group transport would destroy it within about a billion years as shown in a classic paper by Toomre. In this paper we include the low velocity dispersion component, namely gas, on an equal footing with stars in the formulation of the density wave theory, and obtain the dispersion relation for this coupled system. We show that the inclusion of gas makes the group transport slower by a factor of few, thus allowing the pattern to persist longer - for several billion years. Though still less than the Hubble time, this helps in making the spiral structure more long-lived. Further we show that addition of gas is essential to get a stable wave for the observed pattern speed for the Galaxy, which otherwise is not possible for a one-component stellar disc.
The results of the searching on the bases of GPS UKIDSS survey's data of dense compact stellar clusters in the vicinity of 19 YSOs with high and middle masses are presented. Totally we have revealed clusters in 12 areas. Around 5 objects (IRAS 18151-1208, IRAS 18316-0602, IRAS 19110+1045, IRAS 19213+1723, IRAS 20056+3350) they are newly detected. The clusters associated with IRAS 05168+3634, IRAS 20188+3928, IRAS 19374+2352 and IRAS 19388+2357 sources have been already revealed on the less depth data than GPS UKIDSS survey images. The compact groups of stars located in the vicinity of IRAS 05358+3543, IRAS 18507+0121 and IRAS 20198+3716 sources belong to the more extensive clusters. The radii and stellar density have significant gradient: from 0.2 to 2.7 pc and from 3 to 1000 stars/arcmin^2 respectively. In the vicinity of 7 IRAS sources (IRAS 18174-1612, IRAS 18360-0537, IRAS 18385-0512, IRAS 18517+0437, IRAS 19092+0841, IRAS 19410+2336, IRAS 20126+4104) the stellar clusters were not revealed.
The detection of astrophysical very high energy (VHE) neutrinos in the range of TeV-PeV energies by the IceCube observatory has opened a new season in high energy astrophysics. Energies ~PeV imply that the neutrinos are originated from sources where cosmic rays (CRs) can be accelerated up to ~ 10^{17}eV. Recently, we have shown that the observed TeV gamma-rays from radio-galaxies may have a hadronic origin in their nuclear region and in such a case this could lead to neutrino production. In this paper we show that relativistic protons accelerated by magnetic reconnection in the core region of these sources may produce VHE neutrinos via the decay of charged pions produced by photo-meson process. We have also calculated the diffuse flux of VHE neutrinos and found that it can be associated to the IceCube data.
Peak counts have been shown to be an excellent tool to extract the non-Gaussian part of the weak lensing signal. Recently, we developped a fast stochastic forward model to predict weak-lensing peak counts. Our model is able to reconstruct the underlying distribution of observables for analyses. In this work, we explore and compare various strategies for constraining parameter using our model, focusing on the matter density $\Omega_\mathrm{m}$ and the density fluctuation amplitude $\sigma_8$. First, we examine the impact from the cosmological dependency of covariances (CDC). Second, we perform the analysis with the copula likelihood, a technique which makes a weaker assumption compared to the Gaussian likelihood. Third, direct, non-analytic parameter estimations are applied using the full information of the distribution. Fourth, we obtain constraints with approximate Bayesian computation (ABC), an efficient, robust, and likelihood-free algorithm based on accept-reject sampling. We find that neglecting the CDC effect enlarges parameter contours by 22%, and that the covariance-varying copula likelihood is a very good approximation to the true likelihood. The direct techniques work well in spite of noisier contours. Concerning ABC, the iterative process converges quickly to a posterior distribution that is in an excellent agreement with results from our other analyses. The time cost for ABC is reduced by two orders of magnitude. The stochastic nature of our weak-lensing peak count model allows us to use various techniques that approach the true underlying probability distribution of observables, without making simplifying assumptions. Our work can be generalized to other observables where forward simulations provide samples of the underlying distribution.
Cosmological analyses can be accelerated by approximating slow calculations using a training set, which is either precomputed or generated dynamically. However, this approach is only safe if the approximations are well understood and controlled. This paper surveys issues associated with the use of machine-learning based emulation strategies for accelerating cosmological parameter estimation. We describe a learn-as-you-go algorithm that is implemented in the Cosmo++ code and (1) trains the emulator while simultaneously estimating posterior probabilities; (2) identifies unreliable estimates, computing the exact numerical likelihoods if necessary; and (3) progressively learns and updates the error model as the calculation progresses. We explicitly describe and model the emulation error and show how this can be propagated into the posterior probabilities. We apply these techniques to the Planck likelihood and the calculation of $\Lambda$CDM posterior probabilities. The computation is significantly accelerated without a pre-defined training set and uncertainties in the posterior probabilities are subdominant to statistical fluctuations. We have obtained a speedup factor of $6.5$ for Metropolis-Hastings and $3.5$ for nested sampling. Finally, we discuss the general requirements for a credible error model and show how to update them on-the-fly.
We present WHFast, a fast and accurate implementation of a Wisdom-Holman symplectic integrator for long-term orbit integrations of planetary systems. WHFast is significantly faster and conserves energy better than all other Wisdom-Holman integrators tested. We achieve this by significantly improving the Kepler-solver and ensuring numerical stability of coordinate transformations to and from Jacobi coordinates. These refinements allow us to remove the linear secular trend in the energy error that is present in other implementations. For small enough timesteps we achieve Brouwer's law, i.e. the energy error is dominated by an unbiased random walk due to floating-point round-off errors. We implement symplectic correctors up to order eleven that significantly reduce the energy error. We also implement a symplectic tangent map for the variational equations. This allows us to efficiently calculate two widely used chaos indicators the Lyapunov characteristic number (LCN) and the Mean Exponential Growth factor of Nearby Orbits (MEGNO). WHFast is freely available as a flexible C package, as a shared library, and as an easy-to-use python module.
We explore the detailed and broad properties of carbon burning in Super Asymptotic Giant Branch (SAGB) stars with 2755 MESA stellar evolution models. The location of first carbon ignition, quenching location of the carbon burning flames and flashes, angular frequency of the carbon core, and carbon core mass are studied as a function of the ZAMS mass, initial rotation rate, and mixing parameters such as convective overshoot, semiconvection, thermohaline and angular momentum transport. In general terms, we find these properties of carbon burning in SAGB models are not a strong function of the initial rotation profile, but are a sensitive function of the overshoot parameter. We quasi-analytically derive an approximate ignition density, $\rho_{ign} \approx 2.1 \times 10^6$ g cm$^{-3}$, to predict the location of first carbon ignition in models that ignite carbon off-center. We also find that overshoot moves the ZAMS mass boundaries where off-center carbon ignition occurs at a nearly uniform rate of $\Delta M_{\rm ZAMS}$/$\Delta f_{\rm{ov}}\approx$ 1.6 $M_{\odot}$. For zero overshoot, $f_{\rm{ov}}$=0.0, our models in the ZAMS mass range $\approx$ 8.9 to 11 $M_{\odot}$ show off-center carbon ignition. For canonical amounts of overshooting, $f_{\rm{ov}}$=0.016, the off-center carbon ignition range shifts to $\approx$ 7.2 to 8.8 $M_{\odot}$. Only systems with $f_{\rm{ov}}$ $\geq 0.01$ and ZAMS mass $\approx$ 7.2-8.0 $M_{\odot}$ show carbon burning is quenched a significant distance from the center. These results suggest a careful assessment of overshoot modeling approximations on claims that carbon burning quenches an appreciable distance from the center of the carbon core.
We explore the assembly history of the M31 bulge within a projected major-axis radius of 180" (~680 pc) by studying its stellar populations in Hubble Space Telescope WFC3 and ACS observations. Colors formed by comparing near-ultraviolet vs. optical bands are found to become bluer with increasing major-axis radius, which is opposite to that predicted if the sole sources of near-ultraviolet light were old extreme horizontal branch stars with a negative radial gradient in metallicity. Spectral energy distribution fits require a metal-rich intermediate-age stellar population (300 Myr to 1 Gyr old, ~solar metallicity) in addition to the dominant old population. The radial gradients in age and metallicity of the old stellar population are consistent with those in previous works. For the intermediate-age population, we find an increase in age with radius and a mass fraction that increases up to 2% at 680 pc away from the center. We exclude contamination from the M31 disk and/or halo as the main origin for this population. Our results thus suggest that intermediate-age stars exist beyond the central 5" (19 pc) of M31 and contribute ~1% of the total stellar mass in the bulge. These stars could be related to the secular growth of the M31 bulge.
Understanding the nonlinear property of redshift-space distortion, i.e., Fingers-of-God (FoG) effect, is important for the redshift-space distortion studies to test gravity models. FoG effect has been usually attributed to the random motion of galaxies inside the clusters. In this paper, we demonstrate the importance of the coherent infalling motion of satellite galaxies toward the cluster center. We analytically derive the satellite velocity distribution due to the infall motion together with the random motion and show that the velocity distribution becomes far from Maxwellian when the infalling motion is dominant. We use simulated subhalo catalogs to find that the contribution of infall motion is important for massive subhalos and that their velocity distribution has top-hat like shape as expected from our analytic model.We also study the FoG effect due to infall motion on the redshift-space power spectrum. Using the simulated subhalo catalogs based on the halo occupation distribution of luminous red galaxies, we find that the redshift-space power spectra significantly differs from the expectations when the infall motion is ignored. We discuss the impacts on the estimation of the satellite velocity bias.
We use frequency-resolved spectroscopy to examine the energy spectra of the prominent low frequency Quasi-Periodic Oscillation (QPO) and its harmonic in GX 339-4. We track the evolution of these spectra as the source makes a transition from a bright low/hard to hard intermediate state. The time averaged spectrum from each observation can be fit by a disc, and two comptonisation components (one quite soft and one much harder) together with their reflection. The spectrum of the QPO is consistent with being dominated by an even harder comptonisation spectrum than the hard Compton component seen in the time averaged spectrum, together with its reflection. Conversely, the harmonic is consistent with being dominated by the soft Compton component. Simultaneously fitting the time averaged spectrum together with the harmonic leads to significantly better constraints on this component than can be derived from the time averaged spectrum alone. We point to additional evidence from the power spectrum and spectral lags which also argue for the presence of this soft comptonisation component. This adds to the growing evidence for inhomogeneous comptonisation in black hole binaries, and we caution that this more complex continuum model may change the derived broad line profile. The dramatic difference in this object between the spectrum of the QPO and its harmonic disfavours models where the harmonic is solely produced from the angular dependence of Compton scattering in a single region as this would produce a similar spectrum for both harmonic and fundamental. Instead, we suggest that the soft Compton region is located predominantly above the disc while the hard Compton is from the hotter inner flow.
We present global radiative transfer models for heated relativistic jets. The simulations include all relevant radiative processes, starting deep in the opaque zone and following the evolution of radiation to and beyond the photosphere of the jet. The transfer models are compared with three gamma-ray bursts GRB 990123, GRB 090902B, and GRB 130427A, which have well-measured and different spectra. The models provide good fits to the observed spectra in all three cases. The fits give estimates for the jet magnetization parameter $\varepsilon_{\rm B}$ and the Lorentz factor $\Gamma$. In the small sample of three bursts, $\varepsilon_{\rm B}$ varies between 0.01 and 0.1, and $\Gamma$ varies between 340 and 1200.
Asteroseismology provides a powerful way to constrain stellar parameters. Solar-like oscillations have been observed on subgiant stars with the \emph{Kepler\/} mission. The continuous and high-precision time series enables us to carry out a detailed asteroseismic study for these stars. We carry out data processing of two subgiants of spectral type G: KIC 6442183 and KIC 11137075 observed with the \emph{Kepler} mission, and perform seismic analysis for the two evolved stars. We estimate the values of global asteroseismic parameters: $\Delta\nu=64.9\pm 0.2 $ $\mu$Hz and $\nu_{\rm max}=1225 \pm 17$ $\mu$Hz for KIC 6442183, $\Delta\nu=65.5\pm 0.2 $ $\mu$Hz and $\nu_{\rm max}=1171 \pm 8$ $\mu$Hz for KIC 11137075, respectively. In addition, we extract the individual mode frequencies of the two stars. We compare stellar models and observations, including mode frequencies and mode inertias. The mode inertias of mixed modes, which are sensitive to the stellar interior, are used to constrain stellar models. We define a quantity $d\nu_{\rm m-p}$ that measures the difference between the mixed modes and the expected pure pressure modes, which is related to the inertia ratio of mixed modes to radial modes. Asteroseismic together with spectroscopic constraints provide the estimations of the stellar parameters: $M = 1.04_{-0.04}^{+0.01} M_{\odot}$, $R = 1.66_{-0.02}^{+0.01} R_{\odot}$ and $t=8.65_{-0.06}^{+1.12}$ Gyr for KIC 6442183, and $M = 1.00_{-0.01}^{+0.01} M_{\odot}$, $R = 1.63_{-0.01}^{+0.01} R_{\odot}$ and $t=10.36_{-0.20}^{+0.01}$ Gyr for KIC 11137075. Either mode inertias or $d\nu_{\rm m-p}$ could be used to constrain stellar models.
An enigmatic and rare type of young stellar object is the FU Orionis class. The members are interpreted as "outbursting," that is, currently in a state of enhanced accretion by several orders of magnitude relative to the more modest disk-to-star accretion rates measured in typical T Tauri stars. They are key to our understanding of the history of stellar mass assembly and pre-main sequence evolution, as well as critical to consider in the chemical and physical evolution of the circumstellar environment -- where planets form. A common supposition is that *all* T Tauri stars undergo repeated such outbursts, more frequently in their earlier evolutionary stages when the disks are more massive, so as to build up the requisite amount of stellar mass on the required time scale. However, the actual data supporting this traditional picture of episodically enhanced disk accretion are limited, and the observational properties of the known sample of FU Ori objects quite diverse. To improve our understanding of these rare objects, we outline the logic for meaningfully constraining the rate of FU Ori outbursts and present numbers to guide parameter choices in the analysis of time domain surveys.
XTE J1908+094 is an X-ray transient black hole candidate in the Galactic plane that was observed in outburst in 2002 and 2013. Here we present multi-frequency radio and X-ray data, including radio polarimetry, spanning the entire period of the 2013 outburst. We find that the X-ray behaviour of XTE J1908+094 traces the standard black hole hardness-intensity path, evolving from a hard state, through a soft state, before returning to a hard state and quiescence. Its radio behaviour is typical of a compact jet that becomes quenched before discrete ejecta are launched during the late stages of X-ray softening. The radio and X-ray fluxes, as well as the light curve morphologies, are consistent with those observed during the 2002 outburst of this source. The polarisation angle during the rise of the outburst infers a jet orientation in agreement with resolved observations but also displays a gradual drift, which we associate with observed changes in the structure of the discrete ejecta. We also observe an unexpected 90deg rotation of the polarisation angle associated with a second component.
There is much evidence that planet formation is occurring in the disk around the Herbig Be star HD100546. To learn more about the processes occurring in this disk we conducted high resolution imaging at 43/45 GHz with the Australia Telescope Compact Array (ATCA). Multiple array configurations were used, providing a best spatial resolution of $\sim$ 0.15 arcsec, or 15 AU at HD100546's distance of $\sim$ 100 pc. Significant structure is revealed, but its precise form is dependent on the $u-v$ plane sampling used for the image reconstruction. At a resolution of $\leq$ 30 AU we detected an inner gap in the disk with a radius of $\sim$ 25 AU and a position angle approximately along the known disk major axis. With different weighting, and an achieved resolution of $\sim$ 15 AU, emission appears at the centre and the disk takes on the shape of an incomplete ring, much like a horseshoe, again with a gap radius of $\sim$ 25 AU. The position angle of the disk major axis and its inclination from face-on are determined to be $140^{\circ}\pm5^{\circ}$ and $40^{\circ}\pm5^{\circ}$ respectively. The $\sim$ 25 AU gap radius is confirmed by a null in the real part of the binned visibilities at 320$\pm$10 k$\lambda$, whilst the non-axisymmetric nature is also confirmed through significant structure in the imaginary component. The emission mechanism at the central peak is most likely to be free-free emission from a stellar or disk wind. Overall our data support the picture of at least one, but probably several, giant planets orbiting HD100546 within 25 AU.
We present an unprecedented measurement of the disc stability and local instability scales in the luminous infrared Seyfert 1 host, NGC7469, based on ALMA observations of dense gas tracers and with a synthesized beam of 165 x 132 pc. While we confirm that non-circular motions are not significant in redistributing the dense interstellar gas in this galaxy, we find compelling evidence that the dense gas is a suitable tracer for studying the origin of its intensely high-mass star forming ring-like structure. Our derived disc stability parameter accounts for a thick disc structure and its value falls below unity at the radii in which intense star formation is found. Furthermore, we derive the characteristic instability scale and find a striking agreement between our measured scale of ~ 180 pc, and the typical sizes of individual complexes of young and massive star clusters seen in high-resolution images.
In the late stages of nuclear burning for massive stars ($M$> 10 $M_{sun}$), the production of neutrino-antineutrino pairs through various processes becomes the dominant mechanism of stellar cooling. As the star evolves, the energy of these neutrinos increases and in the days preceding the supernova a significant fraction of emitted electron anti-neutrinos exceeds the threshold for inverse beta decay on free hydrogen. This is the golden channel for liquid scintillator detectors because the coincidence signature allows for significant reductions in background signals. We find that the kiloton-scale liquid scintillator detector KamLAND can detect these pre-supernova neutrinos from a star with a mass of 25 $M_{sun}$ at a distance less than 660 pc with 3{\sigma} significance before the supernova. This limit is dependent on the neutrino mass hierarchy and background levels. KamLAND takes data continuously and can provide an alarm for supernovae to the community.
Context. The diffusive shock acceleration mechanism has been widely accepted as the acceleration mechanism for galactic cosmic rays. While self-consistent hybrid simulations have shown how power-law spectra are produced, detailed information on the interplay of diffusive particle motion and the turbulent electromagnetic fields responsible for repeated shock crossings are still elusive. Aims. The framework of test-particle theory is applied to investigate the effect of diffusive shock acceleration by inspecting the obtained cosmic-ray energy spectra. The resulting energy spectra can be obtained this way from the particle motion and, depending on the prescribed turbulence model, the influence of stochastic acceleration through plasma waves can be studied. Methods. A numerical Monte-Carlo simulation code is extended to include collisionless shock waves. This allows one to trace the trajectories of test particle while they are being accelerated. In addition, the diffusion coefficients can be obtained directly from the particle motion, which allows for a detailed understanding of the acceleration process. Results. The classic result of an energy spectrum with $E^{-2}$ is only reproduced for parallel shocks, while, for all other cases, the energy spectral index is reduced depending on the shock obliqueness. Qualitatively, this can be explained in terms of the diffusion coefficients in the directions that are parallel and perpendicular to the shock front.
We investigated the absolute age of the Galactic globular cluster M71 (NGC 6838) by using optical ground-based images (u',g',r',i',z') collected with the MegaCam camera at the Canada-France-Hawaii-Telescope (CFHT). We performed a robust selection of field and cluster stars by applying a new method based on the 3D (r',u'-g',g'-r') Color-Color-Magnitude-Diagram. The comparison between the Color-Magnitude-Diagram of the candidate cluster stars and a new set of isochrones, at the locus of the Main Sequence Turn Off (MSTO), suggests an absolute age of 12+/-2 Gyr. The absolute age was also estimated using the difference in magnitude between the MSTO and the so-called main sequence knee, a well defined bending occurring in the lower main sequence. This feature was originally detected in the near-infrared (NIR) bands and explained as a consequence of an opacity mechanism (collisionally induced absorption of molecular hydrogen) in the atmosphere of cool low-mass stars (Bono et al. 2010). The same feature was also detected in the r',u'-g' and in the r',g'-r' CMD, thus supporting previous theoretical predictions by Borysow et al. (1997). The key advantage in using the Delta^(Knee)_(TO) as an age diagnostic is that it is independent of uncertainties affecting the distance, the reddening and the photometric zero-point. We found an absolute age of 12+/-1 Gyr that agrees, within the errors, with similar age estimates, but the uncertainty is on average a factor of two smaller. We also found that the Delta^(Knee)_(TO) is more sensitive to the metallicity than the MSTO, but the dependence becomes vanishing using the difference in color between the MSK and the MSTO.
Massive black holes in active galaxies are immersed in huge concentrations of late type stars in the galactic bulges and also early type massive stars in the nuclear stellar clusters which are additionally surrounded by quasi-spherical several kpc scale halos containing from a few hundred up to several thousand globular clusters (GCs). It is expected that significant numbers of red giant stars, massive stars and also GCs can move through the jet expelled from the central engine of active galaxy. We consider collisions of stars from the galactic bulge, nuclear cluster and globular clusters with the jet plasma. As a result of such collisions, multiple shocks are expected to appear in the jet around these compact objects. Therefore, the plasma in the kpc scale jet can be significantly disturbed. We show that particles can be accelerated on these shocks up to the multi-TeV energies. TeV leptons emit synchrotron radiation, extending up to the X-ray energies, and also comptonize radiation produced in a stellar cluster and also the Microwave Background Radiation to TeV $\gamma$-ray energies. We show that such non-thermal radiation is likely to be detectable from the intermediate scale jets of the nearby active galaxies for reasonable number of stars and GCs immersed within the jet. As an example, we calculate the expected non-thermal emission in the X-ray and gamma-ray energies from the nearby radio galaxy Cen A from which a steady gamma-ray emission with the complex spectrum has been recently reported by the Fermi and the HESS Observatories.
Recent precise observations of cosmic rays (CRs) by AMS-02 experiment clearly show (1) harder spectra of helium and carbon compared to protons by $\propto R^{0.08}$, and (2) concave breaks in proton and helium spectra at a rigidity $R \sim 300$ GV. In particular the helium and carbon spectra are exactly similar, pointing to the same acceleration site. We examine possible interpretations of these features and identify a chemically enriched region, that is, superbubbles as the most probable origin of Galactic CRs. We further show that CRs originate primarily from the supernova ejecta in the superbubble core, mixed with negligible interstellar medium, and predict similar spectra for heavy nuclei.
Using V band photometry of the WINGS survey, we derive galaxy luminosity functions (LF) in nearby clusters. This sample is complete down to Mv=-15.15, and it is homogeneous, thus allowing the study of an unbiased sample of clusters with different characteristics. We constructed the photometric LF for 72 out of the original 76 WINGS clusters, excluding only those without a velocity dispersion estimate. For each cluster we obtained the LF for galaxies in a region of radius=0.5 x r200, and fitted them with single and double Schechter's functions. We also derive the composite LF for the entire sample, and those pertaining to different morphological classes. Finally we derive the spectroscopic cumulative LF for 2009 galaxies that are cluster members. The double Schechter fit parameters are neither correlated with the cluster velocity dispersion, nor with the X-ray luminosity. Our median values of the Schechter's fit slope are, on average, in agreement with measurements of nearby clusters, but are less steep that those derived from large surveys, such as the SDSS. Early--type galaxies outnumber late-types at all magnitudes, but both early and late types contribute equally to the faint end of the LF. Finally, the spectroscopic LF is in excellent agreement with the ones derived for A2199, A85 and Virgo, and with the photometric one at the bright magnitudes (where both are available). There is a large spread in the LF of different clusters. However, this spread is not caused by correlation of the LF shape with cluster characteristics such as X--ray luminosity or velocity dispersions. The faint end is flatter than what previously derived (alpha_f=-1.7) at odds with what predicted from numerical simulations.
Aims. The Alice far-ultraviolet spectrograph onboard Rosetta is designed to observe emissions from various atomic and molecular species from within the coma of comet 67P/ Churyumov-Gerasimenko and to determine their spatial distribution and evolution with time and heliocentric distance. Methods. Following orbit insertion in August 2014, Alice made observations of the inner coma above the limbs of the nucleus of the comet from cometocentric distances varying between 10 and 80 km. Depending on the position and orientation of the slit relative to the nucleus, emissions of atomic hydrogen and oxygen were initially detected. These emissions are spatially localized close to the nucleus and spatially variable with a strong enhancement above the comet's neck at northern latitudes. Weaker emission from atomic carbon and CO were subsequently detected. Results. Analysis of the relative line intensities suggests photoelectron impact dissociation of H2O vapor as the source of the observed H I and O I emissions. The electrons are produced by photoionization of H2O. The observed C I emissions are also attributed to electron impact dissociation, of CO2, and their relative brightness to H I reflects the variation of CO2 to H2O column abundance in the coma.
IGR J17454-2919 is a hard X-ray transient discovered by INTEGRAL on 2014 September 27. We report on our 20ks Chandra observation of the source, performed about five weeks after the discovery, as well as on INTEGRAL and Swift monitoring long-term observations. X-ray broad-band spectra of the source are compatible with an absorbed power-law, $\Gamma\sim$1.6-1.8, ${\rm N_H}\sim$(10-12)$\times 10^{22}\,{\rm cm}^{-2}$, with no trace of a cut-off in the data up to about 100keV, and with an average absorbed 0.5-100keV flux of about (7.1-9.7)${\times 10^{-10}~erg~cm^{-2}~s^{-1}}$. With Chandra, we determine the most accurate X-ray position of IGR J17454-2919, $\alpha_{J2000}$=17$^{h}$ 45$^{m}$ 27$^{s}$.69, $\delta_{J2000}$= $-$29$^{\circ}$ 19$^{\prime}$ 53$^{\prime \prime}$.8 (90% uncertainty of 0$^{\prime\prime}$.6), consistent with the NIR source 2MASS J17452768-2919534. We also include NIR investigations from our observations of the source field on 2014 October 6 with GROND. With the multi-wavelength information at hand, we discuss the possible nature of IGR J17454-2919.
The detection of mixed oscillation modes offers a unique insight into the internal structure of core helium burning (CHeB) stars. The stellar structure during CHeB is very uncertain because the growth of the convective core, and/or the development of a semiconvection zone, is critically dependent on the treatment of convective boundaries. In this study we calculate a suite of stellar structure models and their non-radial pulsations to investigate why the predicted asymptotic g-mode $\ell = 1$ period spacing $\Delta\Pi_1$ is systematically lower than is inferred from Kepler field stars. We find that only models with large convective cores, such as those calculated with our newly proposed "maximal-overshoot" scheme, can match the average $\Delta\Pi_1$ reported. However, we also find another possible solution that is related to the method used to determine $\Delta\Pi_1$: mode trapping can raise the observationally inferred $\Delta\Pi_1$ well above its true value. Even after accounting for these two proposed resolutions to the discrepancy in average $\Delta\Pi_1$, models still predict more CHeB stars with low $\Delta\Pi_1$ ($ < 270$ s) than are observed. We establish two possible remedies for this: i) there may be a difficulty in determining $\Delta\Pi_1$ for early CHeB stars (when $\Delta\Pi_1$ is lowest) because of the effect that the sharp composition profile at the hydrogen burning shell has on the pulsations, or ii) the mass of the helium core at the flash is higher than predicted. Our conclusions highlight the need for the reporting of selection effects in asteroseismic population studies in order to safely use this information to constrain stellar evolution theory.
Observations by the Fermi Large Area Telescope of gamma-ray millisecond pulsar light curves imply copious pair production in their magnetospheres, and not exclusively in those of younger pulsars. Such pair cascades may be a primary source of Galactic electrons and positrons, contributing to the observed enhancement in positron flux above ~10 GeV. Fermi has also uncovered many new millisecond pulsars, impacting Galactic stellar population models. We investigate the contribution of Galactic millisecond pulsars to the flux of terrestrial cosmic-ray electrons and positrons. Our population synthesis code predicts the source properties of present-day millisecond pulsars. We simulate their pair spectra invoking an offset-dipole magnetic field. We also consider positrons and electrons that have been further accelerated to energies of several TeV by strong intrabinary shocks in black widow and redback systems. Since millisecond pulsars are not surrounded by pulsar wind nebulae or supernova shells, we assume that the pairs freely escape and undergo losses only in the intergalactic medium. We compute the transported pair spectra at Earth, following their diffusion and energy loss through the Galaxy. The predicted particle flux increases for non-zero offsets of the magnetic polar caps. Pair cascades from the magnetospheres of millisecond pulsars are only modest contributors around a few tens of GeV to the lepton fluxes measured by AMS-02, PAMELA, and Fermi, after which this component cuts off. The contribution by black widows and redbacks may, however, reach levels of a few tens of percent at tens of TeV, depending on model parameters.
The primary cosmic ray sources are searched by means of CHRONOTRON - KLARA separate array of the P.N. Lebedev Physical Institute Tien Shan station. It was done on the base of 35 millions observed PCR extensive air showers from 5x10**13 to 5x10**14 eV energies. The data analysis was carried on the method of the direct selection of local areas in equatorial coordinates where the deviation of event numbers exceeded the definite value from normal Gaussian standard. These directions are compared with other arrays observed results and with coordinates of astrophysical sources.
To reconcile the observed unusual high luminosity of periodic source M82X-2 of the first NuSTAR ultraluminous X-ray pulsations with the most extreme violation of the Eddington limit, and in view that a persistent X-ray radiation from M82X-2 ultimately precludes the possibility of typical pulsars, we tackle the problem by the implications of "microscopic theory of black hole", the preceding developments of which are of vital interest for the physics of ultra-high energy (UHE) cosmic-rays. Replacing a central singularity by the infrastructures inside event horizon, subject to certain rules, MTBH explains the origin of ZeV-neutrinos which are of vital interest for the source of UHE- particles. Withal, M82X-2 is assumed to be a spinning intermediate mass black hole resided in final stage of growth. Then the thermal blackbody X-ray emission, arisen due to the rotational kinetic energy of black hole, escapes from event horizon through the vista to outside world that detected as ultraluminous X-ray pulsations. The M82X-2 indeed releases $99.59\%$ of its pulsed radiative energy predominantly in the X-ray bandpass $0.3-30$ keV, while the pulsed radiation over the band $0.3-3$ keV was not detected yet because it is seemingly suppressed by more powerful persistent radiation over this band. We derive a pulse profile and give a quantitative account of energetics and orbital parameters of the semi-detached X-ray binary containing a primary accretor M82X-2 of inferred mass $M\simeq 138.5-226\,M_{\odot}$ and secondary massive, $M_{2}> 48.3- 64.9\,M_{\odot}$, O/B-type donor star with radius of $R> 22.1- 25.7\,R_{\odot}$, respectively.
Twisted magnetic fields should be ubiquitous in the solar corona. The magnetic energy contained in such twisted fields can be released during solar flares and other explosive phenomena. Reconnection in helical magnetic coronal loops results in plasma heating and particle acceleration distributed within a large volume, including the lower coronal and chromospheric sections of the loops, and can be a viable alternative to the standard flare model, where particles are accelerated only in a small volume located in the upper corona. The goal of this study is to investigate the observational signatures of plasma heating and particle acceleration in kink-unstable twisted coronal loops using combination of MHD simulations and test-particle methods. The simulations describe the development of kink instability and magnetic reconnection in twisted coronal loops using resistive compressible MHD, and incorporate atmospheric stratification and large-scale loop curvature. The resulting distributions of hot plasma let us estimate thermal X-ray emission intensities. Test-particle trajectories combined with the density distributions let us deduce synthetic hard X-ray bremsstrahlung intensities. Our simulations emphasise that the geometry of the emission patterns produced in flaring twisted coronal loops can differ from the actual geometry of the underlying magnetic fields. The twist angles revealed by the soft X-ray thermal emission (SXR) are lower than the field-line twist present at the onset of the kink-instability. Hard X-ray (HXR) emission due to the collisions of energetic electrons with the stratified background are concentrated at the loop foot-points, even though the electrons are accelerated everywhere within the coronal volume of the loop. The HXR light-curve are approximately proportional to the temporal derivative of the SXR light-curve.
Algol ($\beta$ Persei) is the prototypical semi-detached eclipsing binary and a hierarchical triple system. From 2006 to 2010 we obtained 121 high-resolution and high-S/N \'{e}chelle spectra of this object. Spectral disentangling yields the individual spectra of all three stars, and greatly improved elements both the inner and outer orbits. We find masses of $M_{\rm A} = 3.39\pm0.06$ M$_\odot$, $M_{\rm B} = 0.770\pm0.009$ M$_\odot$ and $M_{\rm C} = 1.58\pm0.09$ M$_\odot$. The disentangled spectra also give the light ratios between the components in the $B$ and $V$ bands. Atmospheric parameters for the three stars are determined, including detailed elemental abundances for Algol A and Algol C. We find the following effective temperatures: $T_{\rm A} = 12\,550\pm120$ K, $T_{\rm B} = 4900\pm300$ K and $T_{\rm C} = 7550\pm250$ K. The projected rotational velocities are $v_{\rm A} \sin i_{\rm A} = 50.8\pm0.8$ km/s, $v_{\rm B} \sin i_{\rm B} = 62\pm2$ km/s and $v_{\rm C} \sin i_{\rm C} = 12.4\pm0.6$ km/s. This is the first measurement of the rotational velocity for Algol B, and confirms that it is synchronous with the orbital motion. The abundance patterns of components A and C are identical to within the measurement errors, and are basically solar. They can be summarised as mean metal abundances: [M/H]$_{\rm A} = -0.03\pm0.08$ and [M/H]$_{\rm C} = 0.04\pm0.09$. A carbon deficiency is confirmed for Algol A, with tentative indications for a slight overabundance of nitrogen. The ratio of their abundances is (C/N)$_{\rm A} = 2.0\pm0.4$, half of the solar value of (C/N)$_{\odot} = 4.0\pm0.7$. The new results derived in this study, including detailed abundances and metallicities, will enable tight constraints on theoretical evolutionary models for this complex system.
The new spectropolarimetric observations of HD 22920 with ESPaDOnS at CFHT reveal a strong variability of its spectral line profiles with the phase of stellar rotation. We have obtained Teff = 13640 K, logg=3.72 for this star from the best fit of its nine Balmer line profiles. The respective model of stellar atmosphere was calculated to perform abundance analysis of HD 22920 using the spectra obtained for three different phases of stellar rotation. We have found that silicon and chromium abundances appear to be vertically stratified in the atmosphere of HD 22920. Meanwhile, silicon shows hints for a possible variability of vertical abundance stratification with rotational phase.
The BESS-Polar Collaboration measured the energy spectra of cosmic-ray protons and helium during two long-duration balloon flights over Antarctica in December 2004 and December 2007, at substantially different levels of solar modulation. Proton and helium spectra probe the origin and propagation history of cosmic rays in the galaxy, and are essential to calculations of the expected spectra of cosmic-ray antiprotons, positrons, and electrons from interactions of primary cosmic-ray nuclei with the interstellar gas, and to calculations of atmospheric muons and neutrinos. We report absolute spectra at the top of the atmosphere for cosmic-ray protons in the kinetic energy range 0.2-160 GeV and helium nuclei 0.2-80 GeV/nucleon. The corresponding magnetic rigidity ranges are 0.6-160 GV for protons and 1.3-160 GV for helium. These spectra are compared to measurements from previous BESS flights and from AMS-01, ATIC-2, PAMELA, and AMS-02. We also report the ratio of the proton and helium fluxes from 1.3 GV to 160 GV and compare to ratios from PAMELA and AMS-02.
A new research project on spectral analysis that aims to characterize the vertical stratification of element abundances in stellar atmospheres of chemically peculiar (CP) stars is discussed in detail. Some results on detection of vertical abundance stratification in several slowly rotating main sequence CP stars are presented and considered as an indicator of the effectiveness of the atomic diffusion mechanism responsible for the observed peculiarities of chemical abundances. This study is carried out in the frame of Project VeSElkA (Vertical Stratification of Elements Abundance) for which 34 slowly rotating CP stars have been observed with the ESPaDOnS spectropolarimeter at CFHT.
Photosynthesis is a highly efficient mechanism developed by terrestrial life to utilize the energy from photons of solar origin for biological use. Subsurface regions are isolated from the photosphere, and consequently are incapable of utilizing this energy. This opens up the opportunity for life to cultivate alternative mechanisms in order to take advantage of other available energy sources. Studies have shown that in subsurface environments, life can use energy generated from geochemical and geothermal processes to sustain a minimal metabolism. Another mechanism is radiolysis, in which particles emitted by radioactive substances are indirectly utilized for metabolism. One such example is the bacterium fueled by radiation, found 2 miles deep in a South African mine, which consumes hydrogen formed from particles emitted by radioactive U, Th and K present in rock. An additional source of radiation in the subsurface environments is secondary particles, such as muons generated by Galactic Cosmic Rays (GCRs). It is a steady source of a small amount of energy, and the possibility of a slow metabolizing life flourishing on it cannot be ruled out. Muon-induced radiolysis can produce H2 which is used by methanogens for abiotic hydrocarbon synthesis. We propose three mechanisms through which GCR-induced secondary particles, which are able to penetrate in deep subsurface environments, can be utilized for biological use. (1) GCRs injecting energy in the environment through muon-induced radiolysis, (2) organic synthesis from GCR secondaries interacting with the medium and (3) direct capture of radiation with the help of pigments such as melanin. We discuss the implications of these mechanisms on finding life in the Solar System and elsewhere in the Universe.
Recently some pessimism has been expressed about our lack of progress in understanding quasars over the 50+ year since their discovery. It is worthwhile to look back at some of the progress that has been made - but still lies under the radar - perhaps because few people are working on optical/UV spectroscopy in this field. Great advances in understanding quasar phenomenology have emerged using eigenvector techniques. The 4D eigenvector 1 context provides a surrogate H-R Diagram for quasars with a source main sequence driven by Eddington ratio convolved with line-of-sight orientation. Appreciating the striking differences between quasars at opposite ends of the main sequence (so-called population A and B sources) opens the door towards a unified model of quasar physics, geometry and kinematics. We present a review of some of the progress that has been made over the past 15 years, and point out unsolved issues.
In the present study we investigate magnetic reconnection in twisted magnetic fluxtubes with different initial configurations. In all considered cases, energy release is triggered by the ideal kink instability, which is itself the result of applying footpoint rotation to an initially potential field. The main goal of this work is to establish the influence of the field topology and various thermodynamic effects on the energy release process. Specifically, we investigate convergence of the magnetic field at the loop footpoints, atmospheric stratification, as well as thermal conduction. In all cases, the application of vortical driving at the footpoints of an initally potential field leads to an internal kink instability. With the exception of the curved loop with high footpoint convergence, the global geometry of the loop change little during the simulation. Footpoint convergence, curvature and atmospheric structure clearly influences the rapidity with which a loop achieves instability as well as the size of the subsequent energy release. Footpoint convergence has a stabilising influence and thus the loop requires more energy for instability, which means that the subsequent relaxation has a larger heating effect. Large-scale curvature has the opposite result: less energy is needed for instability and so the amount of energy released from the field is reduced. Introducing a stratified atmosphere gives rise to decaying wave phenomena during the driving phase, and also results in a loop that is less stable.
The Spitzer-Cosmic Assembly Deep Near-Infrared Extragalactic Legacy Survey (S-CANDELS; PI G. Fazio) is a Cycle 8 Exploration Program designed to detect galaxies at very high redshifts (z > 5). To mitigate the effects of cosmic variance and also to take advantage of deep coextensive coverage in multiple bands by the Hubble Space Telescope Multi-Cycle Treasury Program CANDELS, S-CANDELS was carried out within five widely separated extragalactic fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the HST Deep Field North, and the Extended Groth Strip. S-CANDELS builds upon the existing coverage of these fields from the Spitzer Extended Deep Survey (SEDS) by increasing the integration time from 12 hours to a total of 50 hours but within a smaller area, 0.16 square degrees. The additional depth significantly increases the survey completeness at faint magnitudes. This paper describes the S-CANDELS survey design, processing, and publicly-available data products. We present IRAC dual-band 3.6+4.5 micron catalogs reaching to a depth of 26.5 AB mag. Deep IRAC counts for the roughly 135,000 galaxies detected by S-CANDELS are consistent with models based on known galaxy populations. The increase in depth beyond earlier Spitzer/IRAC surveys does not reveal a significant additional contribution from discrete sources to the diffuse Cosmic Infrared Background (CIB). Thus it remains true that only roughly half of the estimated CIB flux from COBE/DIRBE is resolved.
HATS-8b is a low density transiting super-Neptune discovered as part of the HATSouth project. The planet orbits its solar-like G dwarf host (V=14.03 $\pm$ 0.10 and T$_{eff}$ =5679 $\pm$ 50 K) with a period of 3.5839 d. HATS-8b is the third lowest mass transiting exoplanet to be discovered from a wide-field ground based search, and with a mass of 0.138 $\pm$ 0.019 M$_J$ it is approximately half-way between the masses of Neptune and Saturn. However HATS-8b has a radius of 0.873 (+0.123,-0.075) R$_J$, resulting in a bulk density of just 0.259 $\pm$ 0.091 g.cm$^{-3}$. The metallicity of the host star is super-Solar ([Fe/H]=0.210 $\pm$ 0.080), arguing against the idea that low density exoplanets form from metal-poor environments. The low density and large radius of HATS-8b results in an atmospheric scale height of almost 1000 km, and in addition to this there is an excellent reference star of near equal magnitude at just 19 arcsecond separation on the sky. These factors make HATS-8b an exciting target for future atmospheric characterization studies, particularly for long-slit transmission spectroscopy.
The structure and properties of the diffuse interstellar medium (ISM) on small scales, sub-au to 1 pc, are poorly understood. We compare interstellar absorption-lines, observed towards a selection of O- and B-type stars at two or more epochs, to search for variations over time caused by the transverse motion of each star combined with changes in the structure in the foreground ISM. Two sets of data were used: 83 VLT- UVES spectra with approximately 6 yr between epochs and 21 McDonald observatory 2.7m telescope echelle spectra with 6 - 20 yr between epochs, over a range of scales from 0 - 360 au. The interstellar absorption-lines observed at the two epochs were subtracted and searched for any residuals due to changes in the foreground ISM. Of the 104 sightlines investigated with typically five or more components in Na I D, possible temporal variation was identified in five UVES spectra (six components), in Ca II, Ca I and/or Na I absorption-lines. The variations detected range from 7\% to a factor of 3.6 in column density. No variation was found in any other interstellar species. Most sightlines show no variation, with 3{\sigma} upper limits to changes of the order 0.1 - 0.3 dex in Ca II and Na I. These variations observed imply that fine-scale structure is present in the ISM, but at the resolution available in this study, is not very common at visible wavelengths. A determination of the electron densities and lower limits to the total number density of a sample of the sightlines implies that there is no striking difference between these parameters in sightlines with, and sightlines without, varying components.
We present a new, fully generative model of optical telescope image sets, along with a variational procedure for inference. Each pixel intensity is treated as a Poisson random variable, with a rate parameter dependent on latent properties of stars and galaxies. Key latent properties are themselves random, with scientific prior distributions constructed from large ancillary data sets. We check our approach on synthetic images. We also run it on images from a major sky survey, where it exceeds the performance of the current state-of-the-art method for locating celestial bodies and measuring their colors.
The `standard' model of cosmology is founded on the basis that the expansion rate of the universe is accelerating at present --- as was inferred originally from the Hubble diagram of Type Ia supernovae. There exists now a much bigger database of supernovae so we can perform rigorous statistical tests to check whether these `standardisable candles' indeed indicate cosmic acceleration. Taking account of the empirical procedure by which corrections are made to their absolute magnitudes to allow for the varying shape of the light curve and extinction by dust, we find, rather surprisingly, that the data are still quite consistent with a constant rate of expansion.
We investigate how the lensing potential can be constrained with future galaxy surveys using their number counts. Such a measurement is an independent test of the standard LCDM framework and can be used to discern modified gravity models. We perform a Fisher matrix forecast based on galaxy angular power spectra, assuming specifications consistent with future photometric Euclid-like surveys, for which we provide a computation and fitting formula of magnification bias, and spectroscopic SKA-like surveys. The analysis suggests that the amplitude of the lensing potential can be constrained at the same level as other standard LCDM cosmological parameters.
Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfv\'{e}n waves, which interact and rapidly generate strong MHD turbulence. The total energy is conserved and the rate of turbulent energy decay is very similar in both codes, although the fluid code has numerical dissipation whereas the kinetic code has kinetic dissipation. The inertial range power spectrum index is similar in both the codes. The fluid code shows a perpendicular wavenumber spectral slope of $k_{\perp}^{-1.3}$. The kinetic code shows a spectral slope of $k_{\perp}^{-1.5}$ for smaller simulation domain, and $k_{\perp}^{-1.3}$ for larger domain. We estimate that collisionless damping mechanisms in the kinetic code can account for the dissipation of the observed nonlinear energy cascade. Current sheets are geometrically characterized. Their lengths and widths are in good agreement between the two codes. The length scales linearly with the driving scale of the turbulence. In the fluid code, their thickness is determined by the grid resolution as there is no explicit diffusivity. In the kinetic code, their thickness is very close to the skin-depth, irrespective of the grid resolution. This work shows that kinetic codes can reproduce the MHD inertial range dynamics at large scales, while at the same time capturing important kinetic physics at small scales.
We have calculated the third order moments of scalar curvature perturbations in configuration space for different inflationary models. We developed a robust numerical technique to compute the bispectrum for different models that have some features in the inflationary potential. From the bispectrum we evaluated moments analytically in the slow roll regime while we devised a numerical mechanism to calculated these moments for non slow roll single field inflationary models with standard kinetic term that are minimally coupled to gravity. With help of these third order moments one can directly predict many non-Gaussian and geometrical measures of CBM distributions in the configuration space. Thus, we have devised a framework to calculate different third order moments and geometrical measures, e.g. Minkowski functionals or skeleton statistic, generated by different single field models of inflation.
We compute corrections to the inflationary potential due to conformally coupled non-relativistic matter. We find that under certain conditions of the matter coupling, inflation may be interrupted abruptly. We display this in the superconformal Starobinsky model, where matter is conformally coupled to the Einstein frame metric. These corrections may easily stop inflation provided that there is an initial density of non-relativistic matter. Since these additional heavy degrees of freedom generically occur in higher dimension theories, for example as Kaluza-Klein modes, this effect can arise in multiple scenarios.
It is generally believed that a cosmological arrow of time must be associated with entropy production. Indeed, in his seminal work on cyclic cosmology, Tolman introduced a viscous fluid in order to make successive expansion/contraction cycles larger than previous ones, thereby generating an arrow of time. However, as we demonstrate in this letter, the production of entropy is not the only means by which a cosmological arrow of time may emerge. Remarkably, systems which are dissipationless may nevertheless demonstrate a preferred direction of time provided they possess attractors. An example is provided by a homogeneous scalar-field driven cyclic cosmology where the presence of cosmological hysteresis causes an arrow of time to emerge in a system which is formally dissipationless.
We emphasize the importance of effects from heavy fields on supergravity models of inflation. We study, in particular, the backreaction of stabilizer fields and geometric moduli in the presence of supersymmetry breaking. Many effects do not decouple even if those fields are much heavier than the inflaton field. We apply our results to successful models of Starobinsky-like inflation and natural inflation. In most scenarios producing a Starobinsky potential it proves difficult to retain the flatness of the potential after backreactions are taken into account. Some of them are incompatible with non-perturbative moduli stabilization. In natural inflation there exist a number of models which are not constrained by backreactions at all. In those cases the correction terms from heavy fields have the same inflaton-dependence as the uncorrected potential, so that inflation may be possible even for very large gravitino masses.
We provide a unified description of cosmological $\alpha$-attractors and late-time acceleration, in excellent agreement with the latest Planck data. Our construction involves two superfields playing distinctive roles: one is the dynamical field and its evolution determines inflation and dark energy, the other is nilpotent and responsible for a landscape of vacua and supersymmetry breaking. We prove that the attractor nature of the theory is enhanced when combining the two sectors: cosmological attractors are very stable with respect to any possible value of the cosmological constant and, interestingly, to any generic coupling of the inflationary sector with the field responsible for uplifting. Finally, as related result, we show how specific couplings generate an arbitrary inflaton potential in a supergravity framework with varying Kahler curvature.
Links to: arXiv, form interface, find, astro-ph, recent, 1506, contact, help (Access key information)
The current hierarchical merging paradigm and $\Lambda$CDM predict that the $z \sim 4-8$ universe should be a time in which the most massive galaxies are transitioning from their initial halo assembly to the later baryonic evolution seen in star-forming galaxies and quasars. However, no evidence of this transition has been found in many high redshift galaxy surveys including CFHTLS, CANDELS and SPLASH, the first studies to probe the high-mass end at these redshifts. Indeed, if halo mass to stellar mass ratios estimated at lower-redshift continue to $z \sim 6-8$, CANDELS and SPLASH report several orders of magnitude more $M \sim 10^{12-13} M_\odot$ halos than are possible to have formed by those redshifts, implying these massive galaxies formed impossibly early. We consider various systematics in the stellar synthesis models used to estimate physical parameters and possible galaxy formation scenarios in an effort to reconcile observation with theory. Although known uncertainties can greatly reduce the disparity between recent observations and cold dark matter merger simulations, even taking the most conservative view of the observations, there remains considerable tension with current theory.
We carry out a direct search for bar-like non-circular flows in intermediate-inclination, gas-rich disk galaxies with a range of morphological types and photometric bar classifications from the first data release (DR1) of the CALIFA survey. We use the DiskFit algorithm to apply rotation only and bisymmetric flow models to H$\alpha$ velocity fields for 49/100 CALIFA DR1 systems that meet our selection criteria. We find satisfactory fits for a final sample of 37 systems. DiskFit is sensitive to the radial or tangential components of a bar-like flow with amplitudes greater than $15\,$km$\,$s$^{-1}$ across at least two independent radial bins in the fit, or ~2.25 kpc at the characteristic final sample distance of ~75 Mpc. The velocity fields of 25/37 $(67.6^{+6.6}_{-8.5}\%)$ galaxies are best characterized by pure rotation, although only 17/25 $(68.0^{+7.7}_{-10.4}\%)$ of them have sufficient H$\alpha$ emission near the galaxy centre to afford a search for non-circular flows. We detect non-circular flows in the remaining 12/37 $(32.4^{+8.5}_{-6.6}\%)$ galaxies. We conclude that the non-circular flows detected in 11/12 $(91.7^{+2.8}_{-14.9}\%)$ systems stem from bars. Galaxies with intermediate (AB) bars are largely undetected, and our detection thresholds therefore represent upper limits to the amplitude of the non-circular flows therein. We find 2/23 $(8.7^{+9.6}_{-2.9}\%)$ galaxies that show non-circular motions consistent with a bar-like flow, yet no photometric bar is evident. This suggests that in ~10% of galaxies either the existence of a bar may be missed completely in photometry or other processes may drive bar-like flows and thus secular galaxy evolution.
We study the rest-frame ultra-violet sizes of massive (~0.8 x 10^11 M_Sun) galaxies at 3.4<z<4.2, selected from the FourStar Galaxy Evolution Survey (ZFOURGE), by fitting single Sersic profiles to HST/WFC3/F160W images from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). Massive quiescent galaxies are very compact, with a median circularized half-light radius r_e = 0.63 +/- 0.18 kpc. Removing 5/16 (31%) sources with signs of AGN activity does not change the result. Star-forming galaxies have r_e = 2.0 +/- 0.60 kpc, 3.2 +/- 1.3 x larger than quiescent galaxies. Quiescent galaxies at z~4 are on average 6.0 +\- 0.17 x smaller than at z~0 and 1.9 +/- 0.7 x smaller than at z~2. Star-forming galaxies of the same stellar mass are 2.4 +/- 0.7 x smaller than at z~0. Overall, the size evolution at 0<z<4 is well described by a powerlaw, with r_e = 5.08 +/- 0.28 (1+z)^(-1.44+/-0.08) kpc for quiescent and r_e = 6.02 +/- 0.28 (1+z)^(-0.72+/-0.05) kpc for star-forming galaxies. Compact star-forming galaxies are rare in our sample: we find only 1/14 (7%) with r_e / (M / 10^11 M_Sun)^0.75 < 1.5, whereas 13/16 (81%) of the quiescent galaxies is compact. The number density of compact quiescent galaxies at z~4 is 1.8 +/- 0.8 x 10^-5 Mpc^-3 and increases rapidly, by >5 x, between 2<z<4. The paucity of compact star-forming galaxies at z~4 and their large rest-frame ultra-violet median sizes suggest that the formation phase of compact cores is very short and/or highly dust obscured.
We present spectral analysis of five NuSTAR and Swift observations of GX 339-4 taken during a failed outburst in summer 2013. These observations cover Eddington luminosity fractions in the range ~0.9-6%. Throughout this outburst, GX 339-4 stayed in the hard state, and all five observations show similar X-ray spectra with a hard power-law with a photon index near 1.6 and significant contribution from reflection. Using simple reflection models we find unrealistically high iron abundances. Allowing for different photon indices for the continuum incident on the reflector relative to the underlying observed continuum results in a statistically better fit and reduced iron abundances. With a photon index around 1.3, the input power-law on the reflector is significantly harder than that which is directly observed. We study the influence of different emissivity profiles and geometries and consistently find an improvement when using separate photon indices. The inferred inner accretion disk radius is strongly model dependent, but we do not find evidence for a truncation radius larger than 100 r_g in any model. The data do not allow independent spin constraints but the results are consistent with the literature (i.e., a>0). Our best-fit models indicate an inclination angle in the range 40-60 degrees, consistent with limits on the orbital inclination but higher than reported in the literature using standard reflection models. The iron line around 6.4 keV is clearly broadened, and we detect a superimposed narrow core as well. This core originates from a fluorescence region outside the influence of the strong gravity of the black hole and we discuss possible geometries.
We consider the dynamical evolution of two planets orbiting in the vicinity of a first order mean motion reso- nance while simultaneously undergoing eccentricity damping and convergent migration. Following Goldreich & Schlichting (2014), we include a coupling between the dissipative semimajor axis evolution and the damping of the eccentricities. In agreement with past studies, we find that this coupling can lead to overstability of the resonance and that for a certain range of parameters capture into resonance is only temporary. Using a more general model, we show that whether overstable motion can occur depends in a characteristic way on the mass ratio between the two planets as well as their relative eccentricity damping timescales. Moreover, we show that even when escape from resonance does occur, the timescale for escape is long enough such at any given time a pair of planets is more likely to be found in a resonance rather than migrating between them. Thus, we argue that overstability of resonances cannot singlehandedly reconcile convergent migration with the observed lack of Kepler planet pairs found near resonances. However, it is possible that overstable motion in combination with other effects such as large scale orbital instability could produce the observed period ratio distribution.
Clusters of galaxies are embedded in halos of optically thin, gravitationally stratified, weakly magnetized plasma at the system's virial temperature. Due to radiative cooling and anisotropic heat conduction, such intracluster medium (ICM) is subject to local instabilities, which are combinations of the thermal, magnetothermal and heat-flux-driven buoyancy instabilities. If the ICM rotates significantly, its stability properties are substantially modified and, in particular, also the magnetorotational instability (MRI) can play an important role. We study simple models of rotating cool-core clusters and we demonstrate that the MRI can be the dominant instability over significant portions of the clusters, with possible implications for the dynamics and evolution of the cool cores. Our results give further motivation for measuring the rotation of the ICM with future X-ray missions such as ASTRO-H and ATHENA.
We assess the impact of starspots on the evolution of late-type stars during the pre-main sequence (pre-MS) using a modified stellar evolution code. We find that heavily spotted models of mass 0.1-1.2\msun\ are inflated by up to $10$% during the pre-MS, and up to 4% and 9% for fully- and partially-convective stars at the zero-age MS, consistent with measurements from active eclipsing binary systems. Spots similarly decrease stellar luminosity and $T_{\rm eff}$, causing isochrone-derived masses to be under-estimated by up to a factor of $2 \times$, and ages to be under-estimated by a factor of 2-10$\times$, at 3 Myr. Consequently, pre-MS clusters and their active stars are systematically older and more massive than often reported. Cluster ages derived with the lithium depletion boundary technique are erroneously young by $\sim 15$% and $10$% at $30$ and $100$ Myr respectively, if 50% spotted stars are interpreted with un-spotted models. Finally, lithium depletion is suppressed in spotted stars with radiative cores, leading to a fixed-temperature lithium dispersion on the MS if a range of spot properties are present on the pre-MS. Such dispersions are large enough to explain Li abundance spreads seen in young open clusters, and imply a range of radii at fixed mass and age during the pre-MS Li burning epoch. By extension, this implies that mass, composition, and age do not uniquely specify the HR diagram location of pre-MS stars.
We present an analysis of deep HST/WFC3 near-IR (NIR) imaging data of the globular cluster M4. The best-photometry NIR colour-magnitude diagram (CMD) clearly shows the main sequence extending towards the expected end of the Hydrogen burning limit and going beyond this point towards fainter sources. The white dwarf sequence can be identified. As such, this is the deepect NIR CMD of a globular cluster to date. Archival HST optical data were used for proper-motion cleaning of the CMD and for distinguishing the WDs from BD candidates. Detection limits in the NIR are around F110W approx. 26.5 mag and F160W approx. 27 mag, and in the optical around F775W approx. 28 mag. Comparing our observed CMDs with theoretical models, we conclude that we have reached beyond the H-burning limit in our NIR CMD and are probably just above or around this limit in our optical-NIR CMDs. Thus, faint NIR sources that have no optical counterpart are potential Brown Dwarf candidates. We visually inspected the positions of NIR sources fainter than H-burning limit in F110W and without optical counterparts in the optical images. We found one such source which is therefore a very good BD candidate. For a further four faint NIR sources upper optical magnitude limits of a possible optical counterpart are estimated, resulting in three additional good BD candidates and one WD candidate.
The goal of this investigation is to reinterpret and upgrade the astrometric and other data on comet C/1945 X1, the least prominent among the Kreutz system sungrazers discovered from the ground in the 20th century. The central issue is to appraise the pros and cons of a possibility that this object is --- despite its brightness reported at discovery --- a dwarf Kreutz sungrazer. We confirm Marsden's (1989) conclusion that C/1945 X1 has a common parent with C/1882 R1 and C/1965 S1, in line with the Sekanina-Chodas (2004) scenario of their origin in the framework of the Kreutz system's evolution. We integrate the orbit of C/1882 R1 back to the early 12th century and then forward to around 1945 to determine the nominal direction of the line of apsides and perform a Fourier analysis to get insight into effects of the indirect planetary perturbations. To better understand the nature of C/1945 X1, its orbital motion, fate, and role in the hierarchy of the Kreutz system as well as to attempt detecting the comet's possible terminal outburst shortly after perihelion and answer the question in the title of this investigation, we closely examined the relevant Boyden Observatory logbooks and identified both the photographs with the comet's known images and nearly 20 additional patrol plates, taken both before and after perihelion, on which the comet or traces of its debris will be searched for, once the process of their digitization, currently conducted as part of the Harvard College Observatory's DASCH Project, has been completed and the scanned copies made available to the scientific community.
Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo finders. We compare five halo finding algorithms' (AHF, HBT, Rockstar, SubFind, and VELOCIraptor) recovery of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z>1). This raises significant uncertainties for theoretical models that require major mergers for, e.g., galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented.
We model particle growth in a turbulent, viscously evolving protoplanetary nebula, incorporating sticking, bouncing, fragmentation, and mass transfer at high speeds. We treat small particles using a moments method and large particles using a traditional histogram binning, including a probability distribution function of collisional velocities. The fragmentation strength of the particles depends on their composition (icy aggregates are stronger than silicate aggregates). The particle opacity, which controls the nebula thermal structure, evolves as particles grow and mass redistributes. While growing, particles drift radially due to nebula headwind drag. Particles of different compositions evaporate at "evaporation fronts" (EFs) where the midplane temperature exceeds their respective evaporation temperatures. We track the vapor and solid phases of each component, accounting for advection and radial and vertical diffusion. We present characteristic results in evolutions lasting $2 \times 10^5$ years. In general, (a) mass is transferred from the outer to inner nebula in significant amounts, creating strong radial concentrations of solids at EFs; (b) particle sizes are limited by a combination of fragmentation, bouncing, and drift; (c) "lucky" large particles never represent a significant amount of mass; and (d) restricted radial zones just outside each EF become compositionally enriched in the associated volatiles. We point out implications for mm-submm SEDs and inference of nebula mass, radial banding, the role of opacity on new mechanisms for generating turbulence, enrichment of meteorites in heavy oxygen isotopes, variable and nonsolar redox conditions, primary accretion of silicate and icy planetesimals, and the makeup of Jupiter's core.
Previous studies have shown that the energy release mechanism of some solar flares follow the Standard magnetic-reconnection model, but the detailed properties of high-energy electrons produced in the flare are still not well understood. We conducted a unique, multi-wavelength study that discloses the spatial, temporal and energy distributions of the accelerated electrons in the X2.2 solar flare on 2011, Feb. 15. We studied the source locations of seven distinct temporal peaks observed in hard X-ray (HXR) and microwave (MW) lightcurves using the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) in 50 to 75 keV channels and Nobeyama Radioheliograph (NoRH) in 34 GHz, respectively. We found that the seven emission peaks did not come from seven spatially distinct sites in HXR and MW, but rather in HXR we observed a sudden change in location only between the second and the third peak, with the same pattern occurring, but evolving more slowly in MW. Comparison between the HXR lightcurve and the temporal variations in intensity in the two MW source kernels also confirmed that the seven peaks came predominantly from two sources, each with multiple temporal peaks. In addition, we studied the polarization properties of MW sources, and time delay between HXR and MW. We discuss our results in the context of the tether-cutting model.
We present a sub-100 pc-scale analysis of the CO molecular gas emission and kinematics of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large Millimetre/submillimetre Array (ALMA) science verification data and a visibility-plane lens reconstruction technique. We find clear evidence for an excitation dependent structure in the unlensed molecular gas distribution, with emission in CO (5-4) being significantly more diffuse and structured than in CO (8-7). The intrinsic line luminosities are 1.04+/-0.13x10^10 and 0.32+/-0.07x10^10 K km/s pc^2 for CO (5-4) and CO (8-7), respectively, after correcting for the differential magnification. A kinemetry analysis of the velocity-fields shows evidence for a perturbed disk with multiple velocity components. Source reconstructions from ALMA and the Hubble Space Telescope are combined to investigate the morphological structure of the stellar, molecular gas and dust components of SDP.81. Together with Karl G. Jansky Very Large Array CO (1-0) data, they provide corroborative evidence for a complex ~2 kpc-scale starburst that is embedded within a larger structure that is ~15 kpc in size.
Light curves of microlensing events involving stellar binaries and planetary systems can provide information about the orbital elements of the system due to orbital modulations of the caustic structure. Accurately measuring the orbit in either the stellar or planetary case requires detailed modeling of subtle deviations in the light curve. At the same time, the natural, Cartesian parameterization of a microlensing binary is partially degenerate with the microlens parallax. Hence, it is desirable to perform independent tests of the predictions of microlens orbit models using radial velocity time series of the lens binary system. To this end, we present 3.5 years of RV monitoring of the binary lens system OGLE-2009-BLG-020L, for which Skowron et al. (2011) constrained all internal parameters of the 200--700 day orbit. Our RV measurements reveal an orbit that is consistent with the predictions of the microlens light curve analysis, thereby providing the first confirmation of orbital elements inferred from microlensing events.
We demonstrate that observations of glitches in the Vela pulsar can be used to investigate the strength of the crust-core coupling in a neutron star, and suggest that recovery from the glitch is dominated by torque exerted by the re-coupling of superfluid components of the core that were decoupled from the crust during the glitch. Assuming that the recoupling is mediated by mutual friction between the superfluid neutrons and the charged components of the core, we use the observed magnitudes and timescales of the shortest timescale components of the recoveries from two recent glitches in the Vela pulsar to infer the fraction of the core that is coupled to the crust during the glitch, and hence spun up by the glitch event. Within the framework of a two-fluid hydrodynamic model of glitches, we analyze whether crustal neutrons alone are sufficient to drive the glitch activity observed in the Vela pulsar. We use two sets of neutron star equations of state (EOSs), both of which span crust and core consistently and cover a range of the slope of the symmetry energy at saturation density $30 < L <120$ MeV. One set produces maximum masses $\approx$2.0$M_{\odot}$, the second $\approx$2.6$M_{\odot}$. We also include the effects of entrainment of crustal neutrons by the superfluid lattice. We find that for medium to stiff EOSs, observations imply $>70\%$ of the moment of inertia of the core is coupled to the crust during the glitch, though for softer EOSs $L\approx 30$MeV as little as $5\%$ could be coupled. No EOS is able to reproduce the observed glitch activity with crust neutrons alone, but extending the region where superfluid vortices are strongly pinned into the core by densities as little as 0.016fm$^{-3}$ above the crust-core transition density restores agreement with the observed glitch activity.
There is ample evidence from rotation curves that dark matter halo's around disk galaxies have nontrivial dynamics. Of particular significance are: a) the cored dark matter profile of disk galaxies, b) correlations of the shape of rotation curves with baryonic properties, and c) the Tully-Fisher relation. Dark matter halo's around disk galaxies may have nontrivial dynamics if dark matter is strongly self interacting and dissipative. Multicomponent hidden sector dark matter featuring a massless `dark photon' (from an unbroken dark $U(1)$ gauge interaction) which kinetically mixes with the ordinary photon provides a concrete example of such dark matter. The kinetic mixing interaction facilitates halo heating by enabling ordinary supernovae to be a source of these `dark photons'. Dark matter halo's can expand and contract in response to the heating and cooling processes, but for a sufficiently isolated halo should have evolved to a steady state or `equilibrium' configuration where heating and cooling rates locally balance. This dynamics allows the dark matter density profile to be related to the distribution of ordinary supernovae in the disk of a given galaxy. In a previous paper a simple and predictive formula was derived encoding this relation. Here we improve on previous work by modelling the supernovae distribution via the measured UV and $H\alpha$ fluxes. The resulting dark matter halo profile is then tested against the rotation curve data of all 26 dwarf galaxies in the LITTLE THINGS sample. The dissipative dark matter concept is further developed and some conclusions drawn.
Modulated (p)reheating is thought to be an alternative mechanism for producing super-horizon curvature perturbations in CMB. But large non-gaussianity and iso-curvature perturbations produced by this mechanism rule out its acceptability as the sole process behind creating CMB perturbations. We explore the situation where CMB perturbations are mostly generated by usual quantum fluctuations of inflaton during inflation, but a modulated coupling constant between inflaton and a secondary scalar affect the preheating process and produces some extra curvature perturbations. We show that to avoid fine tuning of the coupling constant the contribution of modulated preheating in producing curvature perturbations should be less than 15\% of the total observed CMB perturbations.
The long-term (over more than one decade) X-ray emission from two massive stellar systems known to be particle accelerators is investigated using XMM-Newton. Their X-ray properties are interpreted taking into account recent information about their multiplicity and orbital parameters. The two targets, HD168112 and HD167971 appear to be overluminous in X-rays, lending additional support to the idea that a significant contribution of the X-ray emission comes from colliding-wind regions. The variability of the X-ray flux from HD168112 is interpreted in terms of varying separation expected to follow the 1/D rule for adiabatic shocked winds. For HD167971, marginal decrease of the X-ray flux in September 2002 could tentatively be explained by a partial wind eclipse in the close pair. No long-term variability could be demonstrated despite the significant difference of separation between 2002 and 2014. This suggests the colliding-wind region in the wide orbit does not contribute a lot to the total X-ray emission, with a main contribution coming from the radiative shocked winds in the eclipsing pair. The later result provides evidence that shocks in a colliding-wind region may be efficient particle accelerators even in the absence of bright X-ray emission, suggesting particle acceleration may operate in a wide range of conditions. Finally, in hierarchical triple O-type systems, thermal X-rays do not necessarily constitute an efficient tracer to detect the wind-wind interaction in the long period orbit.
The regular or chaotic character of orbits of stars moving in the meridional plane (R,z) of an axially symmetric elliptical galaxy with a dense, massive spherical nucleus and a dark matter halo component is under investigation. In particular, we explore how the flattening of an elliptical galaxy influences the overall orbital structure of the system, by computing in each case the percentage of chaotic orbits, as well as the percentages of orbits composing the main regular families. In an attempt to discriminate safely and with certainty between regular and chaotic motion, we use the Smaller ALingment Index (SALI) method to extensive samples of orbits obtained by integrating numerically the basic equations of motion as well as the variational equations. In addition, a technique which is based mainly on the field of spectral dynamics that utilizes the Fourier transform of the time series of each coordinate is used for classifying the regular orbits into different families and also to recognize the secondary resonances that usually bifurcate from them. Three cases are considered in our work: (i) the case where the elliptical galaxy is prolate (ii) the case where a spherically symmetric elliptical galaxy is present and (iii) the case where the elliptical galaxy has an oblate shape. Comparison between the current results and early related work is also made.
This paper presents the results of the analysis of the very first dedicated X-ray observation with XMM-Newton of WR106. This carbon-rich WC9d Wolf-Rayet star belongs to the category of persistent dust makers (WCd stars). The issue of the multiplicity of these dust makers is pivotal to understand the dust formation process, and in this context X-ray observations may allow to reveal an X-ray emission attributable to colliding-winds in a binary system. The main result of this analysis is the lack of detection of X-rays coming from WR106. Upper limits on the X-ray flux are estimated, but the derived numbers are not sufficient to provide compelling constraints on the existence or not of a colliding-wind region. Detailed inspection of archive data bases reveals that persistent dust makers have been poorly investigated by the most sensitive X-ray observatories. Certainly, the combination of several approaches to indirectly constrain their multiplicity should be applied to lift a part of the veil on the nature of these persistent dust makers.
We present a direct comparison between the observed star formation rate functions (SFRF) and the state-of-the-art predictions of semi-analytic models (SAM) of galaxy formation and evolution. We use the PACS Evolutionary Probe Survey (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) data-sets in the COSMOS and GOODS-South fields, combined with broad-band photometry from UV to sub-mm, to obtain total (IR+UV) instantaneous star formation rates (SFRs) for individual Herschel galaxies up to z~4, subtracted of possible active galactic nucleus (AGN) contamination. The comparison with model predictions shows that SAMs broadly reproduce the observed SFRFs up to z~2, when the observational errors on the SFR are taken into account. However, all the models seem to under-predict the bright-end of the SFRF at z>2. The cause of this underprediction could lie in an improper modelling of several model ingredients, like too strong (AGN or stellar) feedback in the brighter objects or too low fall-back of gas, caused by weak feedback and outflows at earlier epochs.
Modified gravity models with screening in local environments appear in three different guises: chameleon, K-mouflage and Vainshtein mechanisms. We propose to look for differences between these classes of models by considering cosmological observations at low redshift. In particular, we analyse the redshift dependence of the fine structure constant and the proton to electron mass ratio in each of these scenarios. When the absorption lines belong to unscreened regions of space such as dwarf galaxies, a time variation would be present for chameleons. For both K-mouflage and Vainshtein mechanisms, the cosmological time variation of the scalar field is not suppressed in both unscreened and screened environments, therefore enhancing the variation of constants and their detection prospect. We also consider the time variation of the redshift of distant objects using their spectrocopic velocities. We find that models of the K-mouflage and Vainshtein types have very different spectroscopic velocities as a function of redshift and that their differences with the $\Lambda$-CDM template should be within reach of the future ELT- HIRES observations.
Supermassive black hole binaries are one of the primary targets for gravitational wave searches using pulsar timing arrays. Gravitational wave signals from such systems are well represented by parametrized models, allowing the standard Generalized Likelihood Ratio Test (GLRT) to be used for their detection and estimation. However, there is a dichotomy in how the GLRT can be implemented for pulsar timing arrays: there are two possible ways in which one can split the set of signal parameters for semi-analytical and numerical extremization. The straightforward extension of the method used for continuous signals in ground-based gravitational wave searches, where the so-called pulsar phase parameters are maximized numerically, was addressed in an earlier paper (Wang et al. 2014). In this paper, we report the first study of the performance of the second approach where the pulsar phases are maximized semi-analytically. This approach is scalable since the number of parameters left over for numerical optimization does not depend on the size of the pulsar timing array. Our results show that, for the same array size (9 pulsars), the new method performs somewhat worse in parameter estimation, but not in detection, than the previous method where the pulsar phases were maximized numerically. The origin of the performance discrepancy is likely to be in the ill-posedness that is intrinsic to any network analysis method. However, scalability of the new method allows the ill-posedness to be mitigated by simply adding more pulsars to the array. This is shown explicitly by taking a larger array of pulsars.
Variations in the observed dust emission and extinction indicate a systematic evolution of grain properties in the transition from the diffuse interstellar medium (ISM) to denser molecular clouds. The differences in the dust spectral energy distribution (SED) observed from the diffuse ISM to denser regions, namely an increase in the spectral index at long wavelengths, an increase in the FIR opacity, and a decrease in temperature, are usually assumed to be the result of changes in dust properties. We investigate if evolutionary processes, such as coagulation and accretion, are able to change the dust properties of grains in a way that is consistent with observations. We use a core-mantle grain model to describe diffuse ISM-type grains, and using DDA we calculate how the accretion of mantles and coagulation into aggregates vary the grain optical properties. We calculate the dust SED and extinction using DustEM and the radiative transfer code CRT. We show that the accretion of an aliphatic carbon mantle on diffuse ISM-type dust leads to an increase in the FIR opacity by a factor of about 2 and in the FIR/submm spectral index from 1.5 to 1.8, and to a decrease in the temperature by about 2 K. We also show that the coagulation of these grains into aggregates further decreases the temperature by 3 K and increases the spectral index up to a value of $\sim$2. The FIR opacity is increased by a factor of 3 (7) for these aggregates (with an additional ice-mantle) compared to the diffuse ISM-dust. Dust evolution in the ISM resulting from coagulation and accretion, leads to significant changes in the optical properties of the grains that can explain the observed variations in the dust SED in the transition from the diffuse ISM to denser regions.
The High-Altitude Ballooning programme began at Indian Institute of Astrophysics, Bangalore, in the year 2011 with the primary purpose of developing and flying low-cost scientific payloads on a balloon-borne platform. Some of the science goals are studies of the phenomena occurring in the upper atmosphere, of airglow and zodiacal light, and observations of extended astronomical objects such as, for example, comets, from near space (20 to 30 km). A brief summary and results of the tethered flights carried out at CREST campus are given in Ref.~1. Here we present a complete overview of the 9 free-flying balloon experiments conducted from March 2013 to November 2014. We describe the launch procedures, payloads, methods of tracking and recovery of the payloads. Since we fall in the light/medium balloon category, the weight of the payload is limited to less than 5 kg --- we use a 3-D printer to fabricate lightweight boxes and structures for our experiments. We are also developing in-house lightweight sensors and controllers to use in the balloon flights. The flight and scientific data obtained from the different launches, and future plans for the development of a fully-fledged 3-axis pointing and stabilization system and a low-cost star camera-cum-sensor for forthcoming balloon flights are briefly discussed.
We present the detection of a giant radio halo (GRH) in the Sunyaev-Zel'dovich (SZ)-selected merging galaxy cluster ACT-CL J0256.5+0006 ($z = 0.363$), observed with the Giant Metrewave Radio Telescope at 325 MHz and 610 MHz. We find this cluster to host a faint ($S_{610} = 5.6 \pm 1.4$ mJy) radio halo with an angular extent of 2.6 arcmin, corresponding to 0.8 Mpc at the cluster redshift, qualifying it as a GRH. J0256 is one of the lowest-mass systems, $M_{\rm 500,SZ} = (5.0 \pm 1.2) \times 10^{14} M_\odot$, found to host a GRH. We measure the GRH at lower significance at 325 MHz ($S_{325} = 10.3 \pm 5.3$ mJy), obtaining a spectral index measurement of $\alpha^{610}_{325} = 1.0^{+0.7}_{-0.9}$. This result is consistent with the mean spectral index of the population of typical radio halos, $\alpha = 1.2 \pm 0.2$. Adopting the latter value, we determine a 1.4 GHz radio power of $P_{1.4\text{GHz}} = (1.0 \pm 0.3) \times 10^{24}$ W Hz$^{-1}$, placing this cluster within the scatter of known scaling relations. Various lines of evidence, including the ICM morphology, suggest that ACT-CL J0256.5+0006 is composed of two subclusters. We determine a merger mass ratio of 7:4, and a line-of-sight velocity difference of $v_\perp = 1880 \pm 280$ km s$^{-1}$. We construct a simple merger model to infer relevant time-scales in the merger. From its location on the $P_{\rm 1.4GHz}{-}L_{\rm X}$ scaling relation, we infer that we observe ACT-CL J0256.5+0006 approximately 500 Myr before first core crossing.
An Equation of State (\textit{EoS}) closes the set of fluid equations. Although an ideal EoS with a constant \textit{adiabatic index} $\Gamma$ is the preferred choice due to its simplistic implementation, many astrophysical fluid simulations may benefit from a more sophisticated treatment that can account for diverse chemical processes. Here, we first review the basic thermodynamic principles of a gas mixture in terms of its thermal and caloric EoS by including effects like ionization, dissociation as well as temperature dependent degrees of freedom such as molecular vibrations and rotations. The formulation is revisited in the context of plasmas that are either in equilibrium conditions (local thermodynamic- or collisional excitation- equilibria) or described by non-equilibrium chemistry coupled to optically thin radiative cooling. We then present a numerical implementation of thermally ideal gases obeying a more general caloric EoS with non-constant adiabatic index in Godunov-type numerical schemes.We discuss the necessary modifications to the Riemann solver and to the conversion between total energy and pressure (or vice-versa) routinely invoked in Godunov-type schemes. We then present two different approaches for computing the EoS.The first one employs root-finder methods and it is best suited for EoS in analytical form. The second one leans on lookup table and interpolation and results in a more computationally efficient approach although care must be taken to ensure thermodynamic consistency. A number of selected benchmarks demonstrate that the employment of a non-ideal EoS can lead to important differences in the solution when the temperature range is $500-10^4$ K where dissociation and ionization occur. The implementation of selected EoS introduces additional computational costs although using lookup table methods can significantly reduce the overhead by a factor $3\sim 4$.
We present a study of Hen 2-155 and Hen 2-161, two planetary nebulae which
bear striking morphological similarities to other planetary nebulae known to
host close-binary central stars. Both central stars are revealed to be
photometric variables while spectroscopic observations confirm that Hen 2-155
is host to a double-eclipsing, post-common-envelope system with an orbital
period of 3h33m making it one of the shortest period binary central stars
known. The observations of Hen 2-161 are found to be consistent with a
post-common-envelope binary of period ~1 day.
A detailed model of central star of Hen 2-155, is produced, showing the
nebular progenitor to be a hot, post-AGB remnant of approximately 0.62 Msol,
consistent with the age of the nebula, and the secondary star to be an M dwarf
whose radius is almost twice the expected ZAMS radius for its mass. In spite of
the small numbers, all main-sequence companions, of planetary nebulae central
stars, to have had their masses and radii constrained by both photometric and
spectroscopic observations have also been found to display this "inflation".
The cause of the "inflation" is uncertain but is probably related to rapid
accretion, immediately before the recent common-envelope phase, to which the
star has not yet thermally adjusted.
The chemical composition of both nebulae is also analysed, showing both to
display elevated abundance discrepancy factors. This strengthens the link
between elevated \abundance discrepancy factors and close binarity in the
nebular progenitor.
We report results of our search for fast oscillations in lightcurve of one of the brightest accretion powered pulsars on the sky V0332+53 with the help of data of the PCA spectrometer of the RXTE observatory. In course of this search we have carefully explored complications appearing if one uses only sub-bands of the total bandpass of the PCA spectrometer. We show that lightcurves collected in the soft sub-band of the PCA spectrometer contains an additional instrumental noise, lightcurves of harder sub-bands lack some fraction of the anticipated Poisson noise. We show that this noise is caused by a cross-talk of energy bands, which lasts up to ~200usec. One hypothesis is that these effects are caused by temporarily drop of the PCA detector gain after any occurred event due to slowly moving ions in the detector volume. In order to avoid this effect we searched for fast oscillations in flux of V0332+53 only in the total bandpass of the PCA spectrometer 2-60 keV. We have not detected any quasi-periodic oscillations in lightcurve of the source with an upper limit at the level of 0.5% in the Fourier frequency range 200-1500 Hz.
Context: The solution of the nonlocal thermodynamical equilibrium (non-LTE)
radiative transfer equation usually relies on stationary iterative methods,
which may falsely converge in some cases. Furthermore, these methods are often
unable to handle large-scale systems, such as molecular spectra emerging from,
for example, cool stellar atmospheres.
Aims: Our objective is to develop a new method, which aims to circumvent
these problems, using nonstationary numerical techniques and taking advantage
of parallel computers.
Methods: The technique we develop may be seen as a generalization of the
coupled escape probability method. It solves the statistical equilibrium
equations in all layers of a discretized model simultaneously. The numerical
scheme adopted is based on the generalized minimum residual method.
Result:. The code has already been applied to the special case of the water
spectrum in a red supergiant stellar atmosphere. This demonstrates the fast
convergence of this method, and opens the way to a wide variety of
astrophysical problems.
Here we report on the significant role of a so far overlooked dynamical
aspect, namely a secular resonance between the dwarf planet Ceres and other
asteroids. We demonstrate that this type of secular resonance can be the
dominant dynamical factor in certain regions of the main asteroid belt.
Specifically, we performed a dynamical analysis of the asteroids belonging to
the (1726) Hoffmeister family. To identify which dynamical mechanisms are
actually at work in this part of the main asteroid belt, i.e. to isolate the
main perturber(s), we study the evolution of this family in time. The study is
accomplished using numerical integrations of test particles performed within
different dynamical models. The obtained results reveal that the post-impact
evolution of the Hoffmeister asteroid family is a direct consequence of the
nodal secular resonance with Ceres.
This leads us to the conclusion that similar effects must exist in other
parts of the asteroid belt. In this respect, the obtained results shed light on
an important and entirely new aspect of the long-term dynamics of small bodies.
Ceres' fingerprint in asteroid dynamics, expressed through the discovered
secular resonance effect, completely changes our understanding of the way in
which perturbations by Ceres-like objects affect the orbits of nearby bodies.
We present the results of a synthetic spectral analysis of the far ultraviolet archival IUE, HST and FUSE observations of the fast old nova V603 Aql, obtained some 90 years after its 1918 nova outburst. Our analysis utilizes the new Hubble FGS parallax distance for this nearly face-on old nova, a high white dwarf mass and a low reddening. Our analysis includes non-truncated optically thick accretion disks since V603 Aql is neither a polar nor an intermediate polar. Our synthetic spectral modeling of the FUSE and HST spectra analyzed separately indicate a mass transfer rate of 1.5-2.2xe-9 solar mass per year for the FUSE and HST spectra respectively, assuming a WD mass of 1.2Msun. The mass accretion rate also depends on the assumed WD mass, and increases by a factor of two for a WD mass of 0.8Msun. Combining the FUSE and HST spectra together lead to the same results. Potential implications are discussed.
We present an algorithm that can reconstruct the full distributions of metric components within the class of spherically symmetric dust universes that may include a cosmological constant. The algorithm is capable of confronting this class of solutions with arbitrary data. In this work we use luminosity and age data to constrain the geometry of the universe up to a redshift of $z = 1.75$. We show that the current data are perfectly compatible with homogeneous models of the universe and that these models seem to be favoured at low redshift.
The black hole candidate Swift J1753.5-0127 went into outburst in 2005 June. Rather than fade into quiescence as most black-hole-candidate transients do, it has remained in a low hard spectral state for most of the $\sim$ 9 years after outburst. The persistent emission while in a hard state is reminiscent of the black hole Cyg X-1 and the black hole candidates 1E 1740.7-2942 and GRS 1758-258. Thus far hard X-ray/soft gamma-ray results have focused mainly on the 2005 flare, with a few additional observations during 2007. Here, we present results from INTEGRAL/SPI observations from 2005 - 2010 spanning the 22 - 650 keV energy range. Spectral analysis shows a weak high-energy excess ( $\sim$ 2.9 $\sigma$) above a cutoff powerlaw model that is well fit by a powerlaw suggesting an additional spectral component. Observations of Cyg X-1, 1E 1740.7-2942, and GRS 1758-258 have shown similar spectra requiring an additional high-energy component The SPI results are compared to previously reported results for Swift J1753.5-0127 as well as observations of other sources.
Jupiter's tropospheric composition and cloud structure are studied using Cassini VIMS 4.5-5.1 {\mu}m thermal emission spectra from the 2000-2001 flyby. We make use of both nadir and limb darkening observations on the planet's nightside, and compare these with dayside observations. Although there is significant spatial variability in the 5-{\mu}m brightness temperatures, the shape of the spectra remain very similar across the planet, suggesting the presence of a spectrally-flat, spatially inhomogeneous cloud deck. We find that a simple cloud model consisting of a single, compact cloud is able to reproduce both nightside and dayside spectra, subject to the following constraints: (i) the cloud base is located at pressures of 1.2 bar or lower; (ii) the cloud particles are highly scattering; (iii) the cloud is sufficiently spectrally flat. Using this cloud model, we search for global variability in the cloud opacity and the phosphine deep volume mixing ratio. We find that the vast majority of the 5-{\mu}m inhomogeneity can be accounted for by variations in the thickness of the cloud decks, with huge differences between the cloudy zones and the relatively cloud-free belts. The relatively low spectral resolution of VIMS limits reliable retrievals of gaseous species, but some evidence is found for an enhancement in the abundance of phosphine at high latitudes.
The next generation of multi-object spectrographs (MOS) will deliver comprehensive surveys of the Galaxy, Magellanic Clouds and nearby dwarfs. These will provide us with the vast samples, spanning the full extent of the Hertzsprung-Russell diagram, that are needed to explore the chemistry, history and dynamics of their host systems. Further ahead, the Extremely Large Telescopes (ELTs) will have sufficient sensitivity and angular resolution to extend stellar spectroscopy well beyond the Local Group, opening-up studies of the chemical evolution of galaxies across a broad range of galaxy types and environments. In this contribution I briefly reflect on current and future studies of stellar populations, and introduce plans for the MOSAIC instrument for the European ELT.
Zinc in stars is an important reference element because it is a proxy to Fe in studies of damped Lyman-alpha systems, permitting a comparison of chemical evolution histories of bulge stellar populations and DLAs. In terms of nucleosynthesis, it behaves as an alpha element because it is enhanced in metal-poor stars. The aim of this work is to derive the iron-peak element Zn abundances in 56 bulge giants from high resolution spectra. These results are compared with data from other bulge samples, as well as from disk and halo stars, and damped Lyman-alpha systems, in order to better understand the chemical evolution in these environments. High-resolution spectra were obtained using FLAMES+UVES on the Very Large Telescope. We find [Zn/Fe]=+0.24+-0.02 in the range -1.3 < [Fe/H] < -0.5 and [Zn/Fe]=+0.06+-0.02 in the range -0.5 < [Fe/H] < -0.1, whereas for [Fe/H] > -0.1, it shows a spread of -0.60 < [Zn/Fe] < +0.15, with most of these stars having low [Zn/Fe]<0.0. These low zinc abundances at the high metallicity end of the bulge define a decreasing trend in [Zn/Fe] with increasing metallicities. A comparison with Zn abundances in DLA systems is presented, where a dust-depletion correction was applied for both Zn and Fe. Finally, we present a chemical evolution model of Zn enrichment in massive spheroids, representing a typical classical bulge.
Binary stars hosting exoplanets are a unique laboratory where chemical tagging can be performed to measure with high accuracy the elemental abundances of both stellar components, with the aim to investigate the formation of planets and their subsequent evolution. Here, we present a high-precision differential abundance analysis of the XO-2 wide stellar binary based on high resolution HARPS-N@TNG spectra. Both components are very similar K-dwarfs and host planets. Since they formed presumably within the same molecular cloud, we expect they should possess the same initial elemental abundances. We investigate if the presence of planets can cause some chemical imprints in the stellar atmospheric abundances. We measure abundances of 25 elements for both stars with a range of condensation temperature $T_{\rm C}=40-1741$ K, achieving typical precisions of $\sim 0.07$ dex. The North component shows abundances in all elements higher by $+0.067 \pm 0.032$ dex on average, with a mean difference of +0.078 dex for elements with $T_{\rm C} > 800$ K. The significance of the XO-2N abundance difference relative to XO-2S is at the $2\sigma$ level for almost all elements. We discuss the possibility that this result could be interpreted as the signature of the ingestion of material by XO-2N or depletion in XO-2S due to locking of heavy elements by the planetary companions. We estimate a mass of several tens of $M_{\oplus}$ in heavy elements. The difference in abundances between XO-2N and XO-2S shows a positive correlation with the condensation temperatures of the elements, with a slope of $(4.7 \pm 0.9) \times 10^{-5}$ dex K$^{-1}$, which could mean that both components have not formed terrestrial planets, but that first experienced the accretion of rocky core interior to the subsequent giant planets.
We present atmospheric parameters for 51 nearby FG dwarfs uniformly distributed over the -2.60 < [Fe/H] < +0.20 metallicity range that is suitable for the Galactic chemical evolution research. Lines of iron, Fe I and Fe II, were used to derive a homogeneous set of effective temperatures, surface gravities, iron abundances, and microturbulence velocities. We used high-resolution (R>60000) Shane/Hamilton and CFHT/ESPaDOnS observed spectra and non-local thermodynamic equilibrium (NLTE) line formation for Fe I and Fe II in the classical 1D model atmospheres. The spectroscopic method was tested with the 20 benchmark stars, for which there are multiple measurements of the infrared flux method (IRFM) Teff and their Hipparcos parallax error is < 10%. We found NLTE abundances from lines of Fe I and Fe II to be consistent within 0.06 dex for every benchmark star, when applying a scaling factor of S_H = 0.5 to the Drawinian rates of inelastic Fe+H collisions. The obtained atmospheric parameters were checked for each program star by comparing its position in the log g-Teff plane with the theoretical evolutionary track in the Yi et al. (2004) grid. Our final effective temperatures lie in between the T_IRFM scales of Alonso et al. (1996) and Casagrande et al. (2011), with a mean difference of +46 K and -51 K, respectively. NLTE leads to higher surface gravity compared with that for LTE. The shift in log g is smaller than 0.1 dex for stars with either [Fe/H] > -0.75, or Teff < 5750 K, or log g > 4.20. NLTE analysis is crucial for the VMP turn-off and subgiant stars, for which the shift in log g between NLTE and LTE can be up to 0.5 dex. The obtained atmospheric parameters will be used in the forthcoming papers to determine NLTE abundances of important astrophysical elements from lithium to europium and to improve observational constraints on the chemo-dynamical models of the Galaxy evolution.
One of the most exciting near-term prospects in physics is the potential discovery of gravitational waves by the advanced LIGO and Virgo detectors. To maximise both the confidence of the detection and the science return, it is essential to identify an electromagnetic counterpart. This is not trivial, as the events are expected to be poorly localised, particularly in the near-term, with error regions covering hundreds or even thousands of square degrees. In this paper we discuss the prospects for finding an X-ray counterpart to a gravitational wave trigger with the Swift X-ray Telescope, using the assumption that the trigger is caused by a binary neutron star merger which also produces a short gamma-ray burst. We show that it is beneficial to target galaxies within the GW error region, highlighting the need for substantially complete galaxy catalogues out to distances of 300 Mpc. We also show that nearby, on-axis short GRBs are either extremely rare, or are systematically less luminous than those detected to date. We consider the prospects for detecting afterglow emission from an an off-axis GRB which triggered the GW facilities, finding that the detectability, and the best time to look, are strongly dependent on the characteristics of the burst such as circumburst density and our viewing angle.
This work by Feng and papers which published its conclusions do not cite nor do they deal with objections by the author published in 2013-2014. There are many fundamental problems. We will summarize here the principal problems, as published by Melott and Bambach which render irrelevant most of the work presented by Feng.
Context. Asteroids with semi-major axes very close to that of a host planet
can avoid node crossings when their nodal points are at perihelion and at
aphelion. This layout protects the asteroids from close encounters, and
eventual collisions, with the host planet.
Aims. Here, we study the short-term dynamical evolution of four recently
discovered near-Earth asteroids (NEAs) --2012 FC71, 2014 EK24, 2014 QD364, and
2014 UR-- that follow very Earth-like orbits.
Methods. Our analysis is based on results of direct N-body calculations that
use the most updated ephemerides and include perturbations by the eight major
planets, the Moon, the barycentre of the Pluto-Charon system, and the three
largest asteroids.
Results. These four NEAs exhibit an orbital evolution unlike any other known
near-Earth object (NEO). Beyond horseshoe, tadpole, or quasi-satellite
trajectories, they follow co-orbital passing orbits relative to the Earth
within the Kozai domain. Our calculations show that secular interactions induce
librations of their relative argument of perihelion with respect to our planet
but also to Venus, Mars, and Jupiter. Secular chaos is also present. The size
of this transient population is probably large.
Conclusions. Although some of these NEAs can remain orbitally stable for many
thousands of years, their secular dynamics is substantially more complicated
than commonly thought and cannot be properly described within the framework of
the three-body problem alone due to the overlapping of multiple secular
resonances. Objects in this group are amongst the most atypical NEOs regarding
favourable visibility windows because these are separated in time by many
decades or even several centuries.
Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. We perform the first systematic search for gamma-ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV from the ensemble of 147 SNe Type IIn exploding in dense CSM. We search for a gamma-ray excess at each SNe location in a one year time window. In order to enhance a possible weak signal, we simultaneously study the closest and optically brightest sources of our sample in a joint-likelihood analysis in three different time windows (1 year, 6 months and 3 months). For the most promising source of the sample, SN 2010jl (PTF10aaxf), we repeat the analysis with an extended time window lasting 4.5 years. We do not find a significant excess in gamma rays for any individual source nor for the combined sources and provide model-independent flux upper limits for both cases. In addition, we derive limits on the gamma-ray luminosity and the ratio of gamma-ray-to-optical luminosity ratio as a function of the index of the proton injection spectrum assuming a generic gamma-ray production model. Furthermore, we present detailed flux predictions based on multi-wavelength observations and the corresponding flux upper limit at 95% confidence level (CL) for the source SN 2010jl (PTF10aaxf).
We fit model spectral energy distributions to each pixel in 67 nearby (<z>=0.0057) galaxies using broadband photometry from the Sloan Digital Sky Survey and GALEX. For each galaxy, we compare the stellar mass derived by summing the mass of each pixel to that found from fitting the entire galaxy treated as an unresolved point source. We find that, while the pixel-by-pixel and unresolved masses of galaxies with low specific star formation rates (such as ellipticals and lenticulars) are in rough agreement, the unresolved mass estimate for star-forming galaxies is systematically lower then the measurement from spatially-resolved photometry. The discrepancy is strongly correlated with sSFR, with the highest sSFRs in our sample having masses underestimated by 25% (0.12 dex) when treated as point sources. We found a simple relation to statistically correct mass estimates derived from unresolved broad-band SED fitting to the resolved mass estimates: m_{resolved} = m_{unresolved}/(-0.057log(sSFR) + 0.34) where sSFR is in units of yr^{-1}. We study the effect of varying spatial resolution by degrading the image resolution of the largest images and find a sharp decrease in the pixel-by-pixel mass estimate at a physical scale of approximately 3 kpc, which is comparable to spiral arm widths. The effects we observe are consistent with the "outshining" idea which posits that the youngest stellar populations mask more massive, older -- and thus fainter -- stellar populations. Although the presence of strong dust lanes can also lead to a drastic difference between resolved and unresolved mass estimates (up to 45% or 0.3 dex) for any individual galaxy, we found that resolving dust does not affect mass estimates on average. The strong correlation between mass discrepancy and sSFR is thus most likely due to the outshining systematic bias.
The Trojan asteroids provide a unique perspective on the history of Solar System. As a large population of small bodies, they record important gravitational interactions and dynamical evolution of the Solar System. In the past decade, significant advances have been made in understanding physical properties, and there has been a revolution in thinking about the origin of Trojans. The ice and organics generally presumed to be a significant part of Trojan compositions have yet to be detected directly, though low density of the binary system Patroclus (and possibly low density of the binary/moonlet system Hektor) is consistent with an interior ice component. By contrast, fine-grained silicates that appear to be similar to cometary silicates in composition have been detected, and a color bimodality may indicate distinct compositional groups among the Trojans. Whereas Trojans had traditionally been thought to have formed near 5 AU, a new paradigm has developed in which the Trojans formed in the proto-Kuiper Belt, and they were scattered inward and captured in the Trojan swarms as a result of resonant interactions of the giant planets. Whereas the orbital and population distributions of current Trojans are consistent with this origin scenario, there are significant differences between current physical properties of Trojans and those of Kuiper Belt objects. These differences may be indicative of surface modification due to the inward migration of objects that became the Trojans, but understanding of appropriate modification mechanisms is poor and would benefit from additional laboratory studies. Many open questions remain, and the future promises significant strides in our understanding of Trojans. The time is ripe for a spacecraft mission to the Trojans, to turn these objects into geologic worlds that can be studied in detail to unravel their complex history.
The basic structure of the solar system is set by the presence of low-mass terrestrial planets in its inner part and giant planets in its outer part. This is the result of the formation of a system of multiple embryos with approximately the mass of Mars in the inner disk and of a few multi-Earth-mass cores in the outer disk, within the lifetime of the gaseous component of the protoplanetary disk. What was the origin of this dichotomy in the mass distribution of embryos/cores? We show in this paper that the classic processes of runaway and oligarchic growth from a disk of planetesimals cannot explain this dichotomy, even if the original surface density of solids increased at the snowline. Instead, the accretion of drifting pebbles by embryos and cores can explain the dichotomy, provided that some assumptions hold true. We propose that the mass-flow of pebbles is two-times lower and the characteristic size of the pebbles is approximately ten times smaller within the snowline than beyond the snowline (respectively at heliocentric distance $r<r_{ice}$ and $r>r_{ice}$, where $r_{ice}$ is the snowline heliocentric distance), due to ice sublimation and the splitting of icy pebbles into a collection of chondrule-size silicate grains. In this case, objects of original sub-lunar mass would grow at drastically different rates in the two regions of the disk. Within the snowline these bodies would reach approximately the mass of Mars while beyond the snowline they would grow to $\sim 20$ Earth masses. The results may change quantitatively with changes to the assumed parameters, but the establishment of a clear dichotomy in the mass distribution of protoplanets appears robust, provided that there is enough turbulence in the disk to prevent the sedimentation of the silicate grains into a very thin layer.
Context: Several hundred candidate hybrid pulsators of type A-F have been
identified from space-based observations. Their large number allows both
statistical analyses and detailed investigations of individual stars. This
offers the opportunity to study the full interior of the genuine hybrids, in
which both low-radial-order p- and high-order g-modes are self-excited at the
same time. However, a few other physical processes can also be responsible for
the observed hybrid nature, related to binarity or to surface inhomogeneities.
The finding that most delta Scuti stars also show long-period light variations
represents a real challenge for theory.
Methods: Fourier analysis of all the available Kepler light curves.
Investigation of the frequency and period spacings. Determination of the
stellar physical parameters from spectroscopic observations. Modelling of the
transit events.
Results: The Fourier analysis of the Kepler light curves revealed 55
significant frequencies clustered into two groups, which are separated by a gap
between 15 and 27 c/d. The light variations are dominated by the beating of two
dominant frequencies located at around 4 c/d. The amplitudes of these two
frequencies show a monotonic long-term trend. The frequency spacing analysis
revealed two possibilities: the pulsator is either a highly inclined moderate
rotator (v~70 km/s, i > 70 deg) or a fast rotator (v~200 km/s) with i~20 deg.
The transit analysis disclosed that the transit events which occur with a ~197
c/d period may be caused by a 1.6 R_Jup body orbiting a fainter star, which
would be spatially coincident with KIC 9533489.
We complement the PARSEC database of stellar evolutionary tracks with new models of massive stars, from the pre-main sequence phase to the central carbon ignition. We consider a broad range of metallicities, 0.0001$\leq Z \leq$0.04 and initial masses up to $M_{\rm ini}=350\,M_\odot$. The main difference with respect to our previous models of massive stars is the adoption of a recent formalism accounting for the mass-loss enhancement when the ratio of the stellar to the Eddington luminosity, $\Gamma_e$, approaches unity. With this new formalism, the models are able to reproduce the Humphreys-Davidson limit observed in the Galactic and LMC colour-magnitude diagrams, without an ad hoc mass-loss enhancement. We also follow the predictions of recent wind models indicating that the metallicity dependence of the mass-loss rates becomes shallower when $\Gamma_e$ approaches unity. We thus find that the more massive stars may suffer from substantial mass-loss even at low metallicity. We also predict that the Humphreys-Davidson limit should become brighter at decreasing metallicity. We supplement the evolutionary tracks with new tables of theoretical bolometric corrections, useful to compare tracks and isochrones with the observations. For this purpose, we homogenize existing stellar atmosphere libraries of hot and cool stars (PoWR, ATLAS9 and Phoenix) and we add, where needed, new atmosphere models computed with WM-basic. The mass, age and metallicity grids are fully adequate to perform detailed investigations of the properties of very young stellar systems, both in local and distant galaxies. The new tracks supersede the previous old Padova models of massive stars.
In this work, we investigate the thermophysical properties of OSIRIS-REx target asteroid (101955) Bennu (hereafter, Bennu), where thermal inertia plays an important role in understanding the nature of the asteroid's surface, and will definitely provide substantial information for the sampling return mission. Using a thermophysical model incorporating the recently updated 3D radar-derived shape model \citep{Nolan2013} and mid-infrared observations of Spitzer-PUI, Spitzer-IRAC, Herschel/PACS and ESO VLT/VISIR \citep{Muller2012,Emery2014}, we derive the surface thermophysical properties of Bennu. The asteroid has an effective diameter of $510^{+6}_{-40}$ m, a geometry albedo of $0.047^{+0.0083}_{-0.0011}$, a roughness fraction of $0.04^{+0.26}_{-0.04}$, and thermal inertia of $240^{+440}_{-60}\rm~Jm^{-2}s^{-0.5}K^{-1}$ for a best-fit solution at 1$\sigma$ level. The best-estimate thermal inertia indicates that fine-grained regolith may cover a large area of Bennu's surface, with a grain size that may range from $1.3$ to $31$~mm, and our outcome further supports that Bennu would be a suitable target for the OSIRIS-REx mission to return samples from the asteroid to Earth.
The spectral profiles of the coronal Ne viii line at 77 nm have different shapes in quiet-Sun regions and coronal holes (CHs). A single Gaussian fit of the line profile provides an adequate approximation in quiet-Sun areas, whereas a strong shoulder on the long-wavelength side is a systematic feature in CHs. Although this has been noticed since 1999, no physical reason for the peculiar shape could be given. In an attempt to identify the cause of this peculiarity, we address three problems that could not be conclusively resolved in a review article by a study team of the International Space Science Institute (ISSI; Wilhelm et al. 2011) : (1) The physical processes operating at the base and inside of plumes as well as their interaction with the solar wind (SW). (2) The possible contribution of plume plasma to the fast SW streams. (3) The signature of the first-ionization potential (FIP) effect between plumes and inter-plume regions (IPRs). Before the spectroscopic peculiarities in IPRs and plumes in polar coronal holes (PCHs) can be further investigated with the instrument Solar Ultraviolet Measurements of Emitted Radiation (SUMER) aboard the Solar and Heliospheric Observatory (SOHO), it is mandatory to summarize the results of the review to place the spectroscopic observations into context. Finally, a plume model is proposed that satisfactorily explains the plasma flows up and down the plume field lines and leads to the shape of the neon line in PCHs.
We detect emission from [Fe XXI] $\lambda$1354.1, which is a tracer of $10^7$ K gas, in archival HST-COS spectra from the centers of the well-known elliptical galaxies M87 and NGC 4696. The detections are at moderate significance, with S/N of 4.9 and 4.1 respectively. Using this line, we measure the kinematics of the hot gaseous halos in these galaxies, which are stirred by turbulence and bulk flows. The hot gas has a mean velocity which is consistent with zero relative to each galaxy, although in the case of M87 spatial broadening by the off-axis nucleus may be introducing a slight artificial blueshift. In both systems we measure velocity dispersions for this line, which are likely contaminated by spatial broadening. We estimate the effect of spatial broadening and infer turbulent line-of sight velocities of $105^{+28}_{-22}$ km/s and $85^{+22}_{-18}$ km/s, corresponding to turbulent pressures of $7^{+4}_{-3}$% and $5\pm2$% of the total thermal pressure in these respective galaxies. These uncertainties include measurement errors only; the spatial broadening is somewhat uncertain as well so we also provide 90\% upper limits on the turbulent pressure, for which we do not subtract the spatial broadening. These very conservative upper limits are 21% and 17% of the total thermal pressure. We perform emission measure analysis for both sightlines, constraining the characteristic path length and column density of the $\sim10^7$ K gas.
We propose that 'standard pings', brief broadband radio impulses, can be used to study the three-dimensional clustering of matter in the Universe even in the absence of redshift information. The dispersion of radio waves as they travel through the intervening plasma can, like redshift, be used as a cosmological distance measure. Due to inhomogeneities in the electron density along the line-of-sight, dispersion is an imperfect proxy for radial distance and we show that this leads to calculable dispersion-space distortions in the apparent clustering of sources. Fast radio bursts (FRBs) are a new class of radio transients that are the prototypical standard ping, and due to their high observed dispersion, have been interpreted as originating at cosmological distances. The rate of fast radio bursts has been estimated to be several thousand over the whole sky per day, and if cosmological, the sources of these events should trace the large-scale structure of the Universe. We calculate the dispersion-space power spectra for a simple model where electrons and FRBs are biased tracers of the large-scale structure of the Universe and show that the clustering signal could be measured using as few as 10000 events. Such a survey is in-line with what may be achieved with upcoming wide-field radio telescopes.
With a basic varying space-time cutoff $\tilde\ell$, we study a regularized and quantized Einstein-Cartan gravitational field theory and its domains of ultraviolet-unstable fixed point $g_{\rm ir}\gtrsim 0$ and ultraviolet-stable fixed point $g_{\rm uv}\approx 4/3$ of the gravitational gauge coupling $g=(4/3)G/G_{\rm Newton}$. Because the fundamental operators of quantum gravitational field theory are dimension-2 area operators, the cosmological constant is inversely proportional to the squared correlation length $\Lambda\propto \xi^{-2}$. The correlation length $\xi$ characterizes an infrared size of a causally correlate patch of the universe. The cosmological constant $\Lambda$ and the gravitational constant $G$ are related by a generalized Bianchi identity. As the basic space-time cutoff $\tilde\ell$ decreases and approaches to the Planck length $\ell_{\rm pl}$, the universe undergoes inflation in the domain of the ultraviolet-unstable fixed point $g_{\rm ir}$, then evolves to the low-redshift universe in the domain of ultraviolet-stable fixed point $g_{\rm uv}$. We give the quantitative description of the low-redshift universe in the scaling-invariant domain of the ultraviolet-stable fixed point $g_{\rm uv}$, and its deviation from the $\Lambda$CDM can be examined by low-redshift $(z\lesssim 1)$ cosmological observations, such as supernova Type Ia.
We extend the solid-state nudged elastic band method to handle a non-conserved order parameter - in particular, magnetization, that couples to volume and leads to many observed effects in magnetic systems. We apply this formalism to the well-studied magneto-volume collapse during the pressure-induced transformation in iron - from ferromagnetic body-centered cubic (bcc) austenite to hexagonal close-packed (hcp) martensite. We find a bcc-hcp equilibrium coexistence pressure of 8.4 GPa, with the transition-state enthalpy of 156 meV/Fe at this pressure. A discontinuity in magnetization and coherent stress occurs at the transition state, which has a form of a cusp on the potential-energy surface (yet all the atomic and cell degrees of freedom are continuous); the calculated pressure jump of 25 GPa is related to the observed 25 GPa spread in measured coexistence pressures arising from martensitic and coherency stresses in samples. Our results agree with experiments, but necessarily differ from those arising from drag and restricted parametrization methods having improperly constrained or uncontrolled degrees of freedom.
We explore the implications of a new minimal supersymmetric hybrid inflation model in which the MSSM $\mu$ term arises from a coupling of the Higgs doublets to the inflaton sector and plays an essential role both during and after inflation. Successful inflation with the scalar spectral index $n_s = 0.96-0.97$, is followed by a relatively high reheat temperature, $T_{RH} \gtrsim 10^{12}$ GeV, in the presence of this new coupling. Consistency with big bang nucleosynthesis favors a gravitino mass $m_G \gtrsim 5 \times 10^7$ GeV, so that the gravitino decays before the LSP neutralino freezes out. With $\mu \gtrsim m_G \sim 5 \times 10^7$ GeV, and soft scalar masses of the same order, the correct value for the SM-like Higgs boson mass is realized for $\tan\beta \simeq 1.7$. An LSP wino with mass $\simeq$ 2 TeV turns out to be the simplest dark matter candidate. The tensor-to-scalar ratio $r$, a canonical measure of gravity waves, can be as high as $0.001$.
The basic observational properties of the 'coronal partings'--the special type of the coronal magnetic structures, identified by a comparison of the coronal X-ray images and solar magnetograms--are considered. They represent channels inside the unipolar large-scale magnetic fields, formed by the rows of magnetic arcs directed to the neighboring fields of opposite polarity. The most important characteristics of the partings are revealed. It is found that--from the evolutionary and spatial point of view--the partings can transform to the coronal holes and visa versa. The classes of global, intersecting, and complex partings are identified.
We entertain a new paradigm according to which the observed matter-antimatter asymmetry is generated as a large-scale quantum fluctuation over the baryon-symmetric state that occurred during the cosmic inflation.
We formulate $R^2$ pure supergravity as a scale invariant theory built only in terms of superfields describing the geometry of curved superspace. The standard supergravity duals are obtained in both "old" and "new" minimal formulations of auxiliary fields. These theories have massless fields in de Sitter space as they do in their non supersymmetric counterpart. Remarkably, the dual theory of $R^2$ supergravity in the new minimal formulation is an extension of the Freedman model, describing a massless gauge field and a massless chiral multiplet in de Sitter space, with inverse radius proportional to the Fayet-Iliopoulos term. This model can be interpreted as the "de-Higgsed" phase of the dual companion theory of $R+R^2$ supergravity.
We propose a new theory of massive gravity with only two propagating degrees of freedom. After defining the theory in the unitary gauge in the vielbein language, we shall perform a Hamiltonian analysis to count the number of physical degrees of freedom, and then study some phenomenologies. While the homogeneous and isotropic background cosmology and the tensor linear perturbations around it are described by exactly the same equations as those in the de Rham - Gabadadze - Tolley (dRGT) massive gravity, the scalar and vector gravitational degrees of freedom are absent in the new theory at the fully nonlinear level. Hence the new theory provides a stable nonlinear completion of the self-accelerating cosmological solution that was originally found in the dRGT theory.
We recently reported the discovery of an unpublished manuscript by Albert Einstein in which he attempted a 'steady-state' model of the universe, i.e., a cosmic model in which the expanding universe remains essentially unchanged due to a continuous formation of matter from empty space. The manuscript was apparently written in early 1931, many years before the steady-state models of Fred Hoyle, Hermann Bondi and Thomas Gold. We compare Einstein's steady-state cosmology with that of Hoyle, Bondi and Gold and consider the reasons Einstein abandoned his model. The relevance of steady-state models for today's cosmology is briefly reviewed.
This work reviews basic features of both internally heated (IH) convection and Rayleigh-B\'enard (RB) convection, along with findings on IH convection from laboratory experiments and numerical simulations. In the first chapter, six canonical models of convection are described: three configurations of IH convection driven by constant and uniform volumetric heating, and three configurations of RB convection driven by the boundary conditions. The IH models are distinguished by differing pairs of thermal boundary conditions: top and bottom boundaries of equal temperature, an insulating bottom with heat flux fixed at the top, and an insulating bottom with temperature fixed at the top. The RB models also are distinguished by whether temperatures or heat fluxes are fixed at the top and bottom boundaries. Integral quantities important to heat transport are discussed, including the mean fluid temperature, the mean temperature difference between the boundaries, and the mean convective heat transport. Integral relations and bounds are presented, and further bounds are conjectured for the IH cases. The second chapter presents results that can be derived mathematically from the governing equations: linear and nonlinear stability thresholds of static states, and parameter-dependent bounds. Known bounds are on mean temperatures in IH convection and on convective transport in RB convection. The third chapter reviews numerical simulations and laboratory experiments on IH convection, emphasizing quantitative results.
We study general gauge-dependent dynamical equations describing homogeneous
isotropic cosmologies coupled to a scalar field $\psi$ (scalaron). For flat
cosmologies ($k=0$), we analyze in detail the previously proposed
gauge-independent equation describing the differential, $\chi(\alpha)$, of the
map of the metric $\alpha$ to the scalaron $\psi$, which is the main
mathematical characteristic (`portrait') of cosmologies in $\alpha$-version. In
a more habitual $\psi$-version, the similar equation for the differential of
the inverse map, $\bar{\chi}(\psi)$, can be solved asymptotically or for
special scalaron potentials $v(\psi)$.
In the $\alpha$-version the whole dynamical system is explicitly integrable
for $k\neq 0$ and any `potential' $\bar{v}(\alpha)$ replacing $v(\psi)$. There
is no \textit{a priori} relation between the two potentials before deriving
$\chi$, which depends on the potential itself, though relations between the two
pictures can be found in asymptotic regions. An alternative proposal is to
specify a cosmology by assuming a characteristic solution or its phase portrait
and then finding the potentials from the solutions of the dynamical equations.
Our main subject is the mathematical structure of cosmologies, but possible
applications of the results are briefly discussed.
We develop four-parameter supergravity models of inflation and dark energy, constrained so that ${\delta\rho\over \rho}$, $n_s$ and the cosmological constant $\Lambda $ take their known observable values, but where the mass of gravitino $m_{3/2}$ and the tensor-to-scalar ratio $r$ are free parameters. We focus on generalized cosmological $\alpha$-attractor models, with logarithmic Kahler potentials, a nilpotent goldstino and spontaneously broken supersymmetry at the de Sitter minimum. The future data on B-modes will specify the parameter $\alpha$, measuring the geometry of the Kahler, manifold. The string landscape idea for dark energy is supported in these models via an incomplete cancellation of the universal positive goldstino and negative gravitino contribution. The scale of SUSY breaking M related to the mass of gravitino in our models is a controllable parameter, independent on the scale of inflation, it will be constrained by LHC data and future collider Energy-frontier experiments.
Links to: arXiv, form interface, find, astro-ph, recent, 1506, contact, help (Access key information)