Recent observational studies have demonstrated that most satellite galaxies tend to orbit their hosts on highly flattened, vast, possibly co-rotating planes. Two nearly parallel planes of satellites have been confirmed around the M31 galaxy and around the Centaurus A galaxy, while the Milky Way also sports a plane of satellites. It has been argued that such an alignment of satellites on vast planes is unexpected in the standard ($\Lambda$CDM) model of cosmology if not even in contradiction to its generic predictions. Guided by $\Lambda$CDM numerical simulations, which suggest that satellites are channeled towards hosts along the axis of the slowest collapse as dictated by the ambient velocity shear tensor, we re-examine the planes of local satellites systems within the framework of the local shear tensor derived from the Cosmicflows-2 dataset. The analysis reveals that the Local Group and Centaurus A reside in a filament stretched by the Virgo cluster and compressed by the expansion of the Local Void. Four out of five thin planes of satellite galaxies are indeed closely aligned with the axis of compression induced by the Local Void. Being the less massive system, the moderate misalignment of the Milky Way's satellite plane can likely be ascribed to its greater susceptibility to tidal torques, as suggested by numerical simulations. The alignment of satellite systems in the local universe with the ambient shear field is thus in agreement with predictions of the $\Lambda$CDM model.
We have conducted the first parallax and proper motion measurements of 6.7 GHz methanol maser emission using the Australian Long Baseline Array (LBA). The parallax of G339.884$-$1.259 measured from five epochs of observations is 0.48$\pm $0.08 mas, corresponding to a distance of $2.1^{+0.4}_{-0.3}$ kpc, placing it in the Scutum spiral arm. This is consistent (within the combined uncertainty) with the kinematic distance estimate for this source at 2.5$\pm $0.5 kpc using the latest Solar and Galactic rotation parameters. We find from the Lyman continuum photon flux that the embedded core of the young star is of spectral type B1, demonstrating that luminous 6.7 GHz methanol masers can be associated with high-mass stars towards the lower end of the mass range.
Building on our previous cross-correlation analysis (Xia et al. 2011) between the isotropic gamma-ray background (IGRB) and different tracers of the large-scale structure of the universe, we update our results using 60-months of data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. We perform a cross-correlation analysis between the IGRB and objects that may trace the astrophysical sources of the IGRB: SDSS-DR6 QSOs, the SDSS-DR8 Main Galaxy Sample, Luminous Red Galaxies (LRGs) in the SDSS catalog, 2MASS galaxies, and radio NVSS galaxies. The benefit of correlating the Fermi-LAT signal with catalogs of objects at various redshifts is to provide tomographic information on the IGRB which is crucial to separate the various contributions and to clarify its origin. We observe a significant (>3.5 sigma) cross-correlation signal on angular scales smaller than 1 deg in the NVSS, 2MASS and QSO cases and, at lower statistical significance (~3.0 sigma), with SDSS galaxies. These results are robust against the choice of the statistical estimator, estimate of errors, map cleaning procedure and instrumental effects. Finally, we test the hypothesis that the IGRB observed by Fermi-LAT originates from the summed contributions of three types of unresolved extragalactic sources: BL Lacs, FSRQs and Star-Forming Galaxies (SFGs). We find that a model in which the IGRB is mainly produced by SFGs ($72^{+23}_{-37}$% with 2 sigma errors), with BL Lacs and FSRQs giving a minor contribution, provides a good fit to the data. We also consider a possible contribution from Misaligned Active Galactic Nuclei, and we find that, depending on the details of the model and its uncertainty, they can also provide a substantial contribution, partly degenerate with the SFG one. (abridged)
Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the accuracy of present-day $N$-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used $N$-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at $k\leq1$ $h\,\rm Mpc^{-1}$ and to within three percent at $k\leq10$ $h\,\rm Mpc^{-1}$. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code {\tt Pkdgrav3}. We demonstrate that both a minimum box size of $L=500$ $h^{-1}\rm Mpc$ and a maximum particle mass of $M_{\rm p}=10^{9}$ $h^{-1}\rm M_{\odot}$ are required to obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.
If dark matter (DM) is composed by particles which are non-gravitationally coupled to ordinary matter, their annihilations or decays in cosmic structures can result in detectable radiation. We show that the most powerful technique to detect a particle DM signal outside the Local Group is to study the angular cross-correlation of non-gravitational signals with low-redshift gravitational probes. This method allows to enhance signal-to-noise from the regions of the Universe where the DM-induced emission is preferentially generated. We demonstrate the power of this approach by focusing on GeV-TeV DM and on the recent cross-correlation analysis between the 2MASS galaxy catalogue and the Fermi-LAT gamma-ray maps. We show that this technique is more sensitive than other extragalactic gamma-ray probes, such as the energy spectrum and angular autocorrelation of the extragalactic background, and emission from clusters of galaxies. Intriguingly, we find that the measured cross-correlation can be well fitted by a DM component, with thermal annihilation cross section and mass between 10 and 100 GeV, depending on the small-scale DM properties and gamma-ray production mechanism. This solicits further data collection and dedicated analyses.
We present a systematic analysis of the rotation curves of 187 galaxies with masses greater than 10^10 M_sol, with atomic gas masses from the GALEX Arecibo Sloan Survey (GASS), and with follow-up long-slit spectroscopy from the MMT. Our analysis focuses on stellar rotation curves derived by fitting stellar template spectra to the galaxy spectra binned along the slit. In this way, we are able to obtain accurate rotation velocity measurements for a factor of 2 more galaxies than possible with the Halpha line. Galaxies with high atomic gas mass fractions are the most dark-matter dominated galaxies in our sample and have dark matter halo density profiles that are well fit by Navarro, Frenk & White profiles with an average concentration parameter of 10. The inner slopes and of the rotation curves correlate more strongly with stellar population age than with galaxy mass or structural parameters. At fixed stellar mass, the rotation curves of more actively star-forming galaxies have steeper inner slopes than less actively star-forming galaxies. The ratio between the galaxy specific angular momentum and the total specific angular momentum of its dark matter halo, R_j, correlates strongly with galaxy mass, structure and gas content. Low mass, disk-dominated galaxies with atomic gas mass fractions greater than 20% have median values of R_j of around 1, but massive, bulge-dominated galaxies have R_j=0.2-0.3. We argue that these trends can be understood in a picture where gas inflows triggered by disk instabilities lead to the formation of passive, bulge-dominated galaxies with low specific angular momentum.
The presence of multiple populations in globular clusters has been well established thanks to high-resolution spectroscopy. It is widely accepted that distinct populations are a consequence of different stellar generations: intra-cluster pollution episodes are required to produce the peculiar chemistry observed in almost all clusters. Unfortunately, the progenitors responsible have left an ambiguous signature and their nature remains unresolved. To constrain the candidate polluters, we have measured lithium and aluminium abundances in more than 180 giants across three systems: NGC~1904, NGC~2808, and NGC~362. The present investigation along with our previous analysis of M12 and M5 affords us the largest database of simultaneous determinations of Li and Al abundances. Our results indicate that Li production has occurred in each of the three clusters. In NGC~362 we detected an M12-like behaviour, with first and second-generation stars sharing very similar Li abundances favouring a progenitor that is able to produce Li, such as AGB stars. Multiple progenitor types are possible in NGC~1904 and NGC~2808, as they possess both an intermediate population comparable in lithium to the first generation stars and also an extreme population, that is enriched in Al but depleted in Li. A simple dilution model fails in reproducing this complex pattern. Finally, the internal Li variation seems to suggest that the production efficiency of this element is a function of the cluster's mass and metallicity - low-mass or relatively metal-rich clusters are more adept at producing Li.
Much evidence points towards that the photosphere in the relativistic outflow in GRBs plays an important role in shaping the observed MeV spectrum. However, it is unclear whether the spectrum is fully produced by the photosphere or whether a substantial part of the spectrum is added by processes far above the photosphere. Here we make a detailed study of the $\gamma-$ray emission from single pulse GRB110920A which has a spectrum that becomes extremely narrow towards the end of the burst. We show that the emission can be interpreted as Comptonisation of thermal photons by cold electrons in an unmagnetised outflow at an optical depth of $\tau \sim 20$. The electrons receive their energy by a local dissipation occurring close to the saturation radius. The main spectral component of GRB110920A and its evolution is thus, in this interpretation, fully explained by the emission from the photosphere including localised dissipation at high optical depths.
Two-point correlation functions of cosmic microwave background polarization provide a physically independent probe of the surprising suppression of correlations in the cosmic microwave background temperature anisotropies at large angular scales. We investigate correlation functions constructed from both the Q and U Stokes parameters and from the E and B polarization components. The dominant contribution to these correlation functions comes from local physical effects at the last scattering surface or from the epoch of reionization at high redshift, so all should be suppressed if the temperature suppression is due to an underlying lack of correlations in the cosmological metric perturbations larger than a given scale. We evaluate the correlation functions for the standard $\Lambda$CDM cosmology constrained by the observed temperature correlation function, and compute statistics characterizing their suppression on large angular scales. Future full-sky polarization maps with minimal systematic errors on large angular scales will provide strong tests of whether the observed temperature correlation function is a statistical fluke or reflects a fundamental shortcoming of the standard cosmological model.
The collimation and evolution of relativistic outflows in $\gamma$-ray bursts (GRBs) are determined by their interaction with the stellar envelope through which they travel before reaching the much larger distance where the energy is dissipated and $\gamma$-rays are produced. We consider the case of a Poynting flux dominated relativistic outflow and show that it suffers strong inverse-Compton (IC) scattering drag near the stellar surface and the jet is slowed down to sub-relativistic speed if its initial magnetization parameter ($\sigma_0$) is larger than about 10$^5$. If the temperature of the cocoon surrounding the jet were to be larger than about 10 keV, then an optically thick layer of electrons and positrons forms at the interface of the cocoon and the jet, and one might expect this pair screen to protect the interior of the jet from IC drag. However, the pair screen turns out to be ephemeral, and instead of shielding the jet it speeds up the IC drag on it. Although a high $\sigma_0$ jet might not survive its passage through the star, a fraction of its energy is converted to 1-100 MeV radiation that escapes the star and appears as a bright flash lasting for about 10 s.
We present an analysis of archival Spitzer images and new ground-based and Hubble Space Telescope (HST) near-infrared (IR) and optical images of the field of M83 with the goal of identifying rare, dusty, evolved massive stars. We present point source catalogs consisting of 3778 objects from $Spitzer$ Infrared Array Camera (IRAC) Band 1 (3.6 $\mu$m) and Band 2 (4.5 $\mu$m), and 975 objects identified in Magellan 6.5m FourStar near-IR $J$ and $K_{\rm s}$ images. A combined catalog of coordinate matched near- and mid-IR point sources yields 221 objects in the field of M83. Using this photometry we identify 185 massive evolved stellar candidates based on their location in color-magnitude and color-color diagrams. We estimate the background contamination to our stellar candidate lists and further classify candidates based on their appearance in $HST$ Wide Field Camera 3 (WFC3) observations of M83. We find 49 strong candidates for massive stars which are very promising objects for spectroscopic follow-up. Based on their location in a $B-V$ versus $V-I$ diagram, we expect at least 24, or roughly 50%, to be confirmed as red supergiants.
Spatial variations in the distribution of galaxy luminosities, estimated from redshifts as distance proxies, are correlated with the peculiar velocity field. Comparing these variations with the peculiar velocities inferred from galaxy redshift surveys is a powerful test of gravity and dark energy theories on cosmological scales. Using ~ 2 $\times$ 10$^{5}$ galaxies from the SDSS Data Release 7, we perform this test in the framework of gravitational instability to estimate the normalized growth rate of density perturbations f$\sigma_{8}$ = 0.37 +/- 0.13 at z ~ 0.1, which is in agreement with the $\Lambda$CDM scenario. This unique measurement is complementary to those obtained with more traditional methods, including clustering analysis. The estimated accuracy at z ~ 0.1 is competitive with other methods when applied to similar datasets.
We perform a systematic study to explore the accuracy with which the parameters of intermediate-mass black-hole binary systems can be measured from their gravitational wave (GW) signatures using second-generation GW detectors. We make use of the most recent reduced-order models containing inspiral, merger and ringdown signals of aligned-spin effective-one-body waveforms (SEOBNR) to significantly speed up the calculations. We explore the phenomenology of the measurement accuracies for binaries with total masses between 50 and 500 $M_\odot$ and mass ratios between 0.1 and 1. We find that (i) at total masses below ~200 $M_\odot$, where the signal-to-noise-ratio is dominated by the inspiral portion of the signal, the chirp mass parameter can be accurately measured; (ii) at higher masses, the information content is dominated by the ringdown, and total mass is measured more accurately; (iii) the mass of the lower-mass companion is poorly estimated, especially at high total mass and more extreme mass ratios; (iv) spin cannot be accurately measured for our injection set with non-spinning components. Most importantly, we find that for binaries with non-spinning components at all values of the mass ratio in the considered range and at network signal-to-noise ratio of 15, analyzed with spin-aligned templates, the presence of an intermediate-mass black hole with mass >100 $M_\odot$ can be confirmed with 95% confidence in any binary that includes a component with a mass of 130 $M_\odot$ or greater.
Particle methods are a ubiquitous tool for solving the Vlasov-Poisson equation in comoving coordinates, which is used to model the gravitational evolution of dark matter in an expanding universe. However, these methods are known to produce poor results on idealized test problems, particularly at late times, after the particle trajectories have crossed. To investigate this, we have performed a series of one- and two-dimensional "Zel'dovich Pancake" calculations using the popular Particle-in-Cell (PIC) method. We find that PIC can indeed converge on these problems provided the following modifications are made. The first modification is to regularize the singular initial distribution function by introducing a small but finite artificial velocity dispersion. This process is analogous to artificial viscosity in compressible gas dynamics, and, as with artificial viscosity, the amount of regularization can be tailored so that its effect outside of a well-defined region - in this case, the high-density caustics - is small. The second modification is the introduction of a particle remapping procedure that periodically re-expresses the dark matter distribution function using a new set of particles. We describe a remapping algorithm that is third-order accurate and adaptive in phase space. This procedure prevents the accumulation of numerical errors in integrating the particle trajectories from growing large enough to significantly degrade the solution. Once both of these changes are made, PIC converges at second order on the Zel'dovich Pancake problem, even at late times, after many caustics have formed. Furthermore, the resulting scheme does not suffer from the unphysical, small-scale "clumping" phenomenon known to occur on the Pancake problem when the perturbation wave vector is not aligned with one of the Cartesian coordinate axes.
We present results from multiwavelength observations of the galaxy NGC5005. We use new neutral hydrogen (HI) observations from the Very Large Array to examine the neutral gas morphology and kinematics. We find an HI disk with a well-behaved flat rotation curve in the radial range 20\arcsec-140\arcsec. Ionized gas observations from the SparsePak integral field unit on the WIYN 3.5m telescope provide kinematics for the central 70\arcsec. We use both the SparsePak and HI velocity fields to derive a rotation curve for NGC5005. Deep 3.6{\mu}m observations from the Spitzer Space Telescope probe the faint extended stellar population of NGC5005. The images reveal a large stellar disk with a high surface brightness component that transitions to a low surface brightness component at a radius nearly 1.6 times farther than the extent of the gas disk detected in HI. The 3.6{\mu}m image is also decomposed into bulge and disk components to account for the stellar light distribution. Optical broadband B and R and narrowband H{\alpha} from the WIYN 0.9m telescope complement the 3.6{\mu}m data by providing information about the dominant stellar population and current star formation activity. The neutral and ionized gas rotation curve is used along with the stellar bulge and disk light profiles to decompose the mass distributions in NGC5005 and determine a dark matter halo model. The maximum stellar disk contribution to the total rotation curve is only about 70\%, suggesting that dark matter makes a significant contribution to the dynamics at all radii.
Transition edge sensor microcalorimeters can measure x-ray and gamma-ray energies with very high energy resolution and high photon-collection efficiency. For this technology to reach its full potential in future x-ray observatories, each sensor must be able to measure hundreds or even thousands of photon energies per second. Current "optimal filtering" approaches to achieve the best possible energy resolution work only for photons well isolated in time, a requirement in direct conflict with the need for high-rate measurements. We describe a new analysis procedure to allow fitting for the pulse height of all photons even in the presence of heavy pulse pile-up. In the limit of isolated pulses, the technique reduces to the standard optimal filtering with long records. We employ reasonable approximations to the noise covariance function in order to render multi-pulse fitting computationally viable even for very long data records. The technique is employed to analyze x-ray emission spectra at 600 eV and 6 keV at rates up to 250 counts per second in microcalorimeters having exponential signal decay times of approximately 1.2 ms.
We study the Cosmic Microwave Background using the three-scale framework of Hu et al. to derive the dependence of the CMB temperature anisotropy spectrum on the fundamental constants. We show that, as expected, the observed spectrum depends only on \emph{dimensionless} combinations of the constants, and we emphasize the points that make this generally true for cosmological observations. Our analysis suggests that the CMB spectrum shape is mostly determined by $\alpha^2m_e/m_p$ and the proton-CDM-particle mass ratio, $m_p/\mchi$, with a sub-dominant dependence on $(G\mchi m_e/\hbar c)\alpha^\beta$ with $\beta\sim -7$. The distance to the last-scattering surface depends on $Gm_p\mchi/\hbar c$, so published CMB observational limits on time variations of the constants, besides making assumptions about the form of the dark-energy, implicitly assume the time-independence of this quantity. On the other hand, low-redshift $H_0$, BAO and large-scale structure data can be combined with the \emph{shape} of the CMB spectrum to give information that is largely independent of the dark-energy model. In particular we show that the pre-recombination values of $G\mchi^2/\hbar c$ and $\alpha^2m_e/\mchi$ could not have differed from their present values by more than of order 25\%.
Ellipsoidal red giant binaries are close binary systems where an unseen, relatively close companion distorts the red giant, leading to light variations as the red giant moves around its orbit. These binaries are likely to be the immediate evolutionary precursors of close binary planetary nebula and post-asymptotic giant branch and post-red giant branch stars. Due to the MACHO and OGLE photometric monitoring projects, the light variability nature of these ellipsoidal variables has been well studied. However, due to the lack of radial velocity curves, the nature of their masses, separations, and other orbital details has so far remained largely unknown. In order to improve this situation, we have carried out spectral monitoring observations of a large sample of 80 ellipsoidal variables in the Large Magellanic Cloud and we have derived radial velocity curves. At least 12 radial velocity points with good quality were obtained for most of the ellipsoidal variables. The radial velocity data are provided with this paper. Combining the photometric and radial velocity data, we present some statistical results related to the binary properties of these ellipsoidal variables.
Hubble Space Telescope Advanced Camera for Surveys has been used to determine accurate distances for the spiral galaxy NGC2683 and 12 other galaxies in a zone of the "local velocity anomaly" from measurements of the luminosities of the brightest red giant branch stars. These galaxies lie in the Leo Spur, the nearest filament beyond our home Local Sheet. The new, accurate distance measurements confirm that galaxies along the Leo Spur are more distant than expected from uniform cosmic expansion, hence have large peculiar velocities toward us. The motions are generally explained by a previously published model that posits that the Local Sheet is descending at 259 km/s toward the south supergalactic pole due to expansion of the Local Void and being attracted at 185 km/s toward the Virgo Cluster. With the standard $\Lambda$CDM cosmology an empty void expands at 16 km/s/Mpc so a motion of 259 km/s requires the Local Void to be impressively large and empty. Small residuals from the published model can be attributed to an upward push toward the north supergalactic pole by expansion of the Gemini-Leo Void below the Leo Spur. The Leo Spur is sparsely populated but among its constituents there are two associations that contain only dwarf galaxies.
Acoustic modes of oscillation are affected by stellar activity, however it is unclear how starspots contribute to these changes. Here we investigate the non-magnetic effects of starspots on global modes with angular degree $\ell \leq 2$ in highly active stars, and characterize the spot seismic signature on synthetic light curves. We perform 3D time-domain simulations of linear acoustic waves to study their interaction with a model starspot. We model the spot as a 3D change in the sound speed stratification with respect to a convectively stable stellar background, built from solar Model S. We perform a parametric study by considering different depths and perturbation amplitudes. Exact numerical simulations allow investigation of the wavefield-spot interaction beyond first order perturbation theory. The interaction of the axisymmetric modes with the starspot is strongly nonlinear. As mode frequency increases, the frequency shifts for radial modes exceed the value predicted by linear theory, while the shifts for the $\ell=2, m=0$ modes are smaller than predicted by linear theory, with avoided-crossing-like patterns forming between the $m=0$ and $m=1$ mode frequencies. The nonlinear behavior increases with increasing spot amplitude and/or decreasing depth. Linear theory still reproduces the correct shifts for non-axisymmetric modes. In the nonlinear regime the mode eigenfunctions are not pure spherical harmonics, but rather a mixture of different spherical harmonics. This mode mixing, together with the frequency changes, may lead to misidentification of the modes in the observed acoustic power spectra.
We present the redshift-space generalization of the equal-time angular-averaged consistency relations between $(\ell+n)$- and $n$-point polyspectra of the cosmological matter density field. Focusing on the case of $\ell=1$ large-scale mode and $n$ small-scale modes, we use an approximate symmetry of the gravitational dynamics to derive explicit expressions that hold beyond the perturbative regime, including both the large-scale Kaiser effect and the small-scale fingers-of-god effects. We explicitly check these relations, both perturbatively, for the lowest-order version that applies to the bispectrum, and nonperturbatively, for all orders but for the one-dimensional dynamics. Using a large ensemble of $N$-body simulations, we find that our squeezed bispectrum relation is valid to better than $20\%$ up to $1h$Mpc$^{-1}$, for both the monopole and quadrupole at $z=0.35$, in a $\Lambda$CDM cosmology. Additional simulations done for the Einstein-de Sitter background suggest that these discrepancies mainly come from the breakdown of the approximate symmetry of the gravitational dynamics. For practical applications, we introduce a simple ansatz to estimate the new derivative terms in the relation using only observables. Although the relation holds worse after using this ansatz, we can still recover it within $20\%$ up to $1h$Mpc$^{-1}$, at $z=0.35$ for the monopole. On larger scales, $k = 0.2 h\mathrm{Mpc}^{-1}$, it still holds within the statistical accuracy of idealized simulations of volume $\sim8h^{-3}\mathrm{Gpc}^3$ without shot-noise error.
We investigate the improvements in Shack-Hartmann wavefront sensor image processing that can be realised using total variation minimisation techniques to remove noise from these images. We perform Monte-Carlo simulation to demonstrate that at certain signal-to-noise levels, sensitivity improvements of up to one astronomical magnitude can be realised. We also present on-sky measurements taken with the CANARY adaptive optics system that demonstrate an improvement in performance when this technique is employed, and show that this algorithm can be implemented in a real-time control system. We conclude that total variation minimisation can lead to improvements in sensitivity of up to one astronomical magnitude when used with adaptive optics systems.
We present measurements of relative elemental abundances in plumes and interplumes. Plumes are bright, narrow structures in coronal holes that extend along open magnetic field lines far out into the corona. Previous work has found that in some coronal structures the abundances of elements with a low first ionization potential (FIP) < 10 eV are enhanced relative to their photospheric abundances. This coronal-to-photospheric abundance ratio, commonly called the FIP bias, is typically 1 for element with a high-FIP (> 10 eV). We have used EIS spectroscopic observations made on 2007 March 13 and 14 over an ~24 hour period to characterize abundance variations in plumes and interplumes. To assess their elemental composition, we have used a differential emission measure (DEM) analysis, which accounts for the thermal structure of the observed plasma. We have used lines from ions of iron, silicon, and sulfur. From these we have estimated the ratio of the iron and silicon FIP bias relative to that for sulfur. From the results, we have created FIP-bias-ratio maps. We find that the FIP-bias ratio is sometimes higher in plumes than in interplumes and that this enhancement can be time dependent. These results may help to identify whether plumes or interplumes contribute to the fast solar wind observed in situ and may also provides constraints on the formation and heating mechanisms of plumes.
We have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.
In this work we present the results of a suite of dynamical simulations following the orbital evolution of 8,022 hypothetical Centaur objects. These Centaurs begin our integrations on orbits in the vicinity of the 2:1 mean motion resonance with Neptune, and we follow their dynamical evolution for a period of 3 Myr under the gravitational influence of a motionless Sun and the four Jovian planets. The great majority of the test particles studied rapidly escaped from the vicinity of the 2:1 mean motion resonance of Neptune and diffused throughout the Solar System. The average libration time of Centaurs in the vicinity of 2:1 mean motion resonance of Neptune was found to be just 27 kyr. Although two particles did remain near the resonance for more than 1 Myr. Upon leaving the vicinity of the 2:1 resonance, the majority of test particles evolved by a process of random walk in semi-major axis, due to repeated close encounters with the giant planets.
The Pseudo-Nambu-Goldstone Boson (PNGB) potential, defined through the amplitude $M^4$ and width $f$ of its characteristic potential $V(\phi) = M^4[1 + \cos(\phi~ /~ f)]$, is one of the best-suited models for the study of thawing quintessence. We analyse its present observational constraints by direct numerical solution of the scalar field equation of motion. Observational bounds are obtained using data from Union 2.1 for Supernovae, cosmic microwave background temperature anisotropies from Planck plus WMAP polarization data, and baryon acoustic oscillations data. We find the parameter ranges for which PNGB quintessence remains a viable theory for dark energy. We compare the direct potential analysis and use of an approximate equation-of-state parameterization for thawing theories; this comparison highlights a strong prior dependence to the outcome coming from the choice of modelling methodology, which current data are not sufficient to override.
We have conducted a detailed object-by-object study of a mass-complete (M*>10^11 M_sun) sample of 56 galaxies at 1.4 < z < 2 in the GOODS-South field, showing that an accurate de-blending in MIPS/24um images is essential to properly assign to each galaxy its own star formation rate (SFR), whereas an automatic procedure often fails. This applies especially to galaxies with SFRs below the Main Sequence (MS) value, which may be in their quenching phase. After that, the sample splits evenly between galaxies forming stars within a factor of 4 of the MS rate (~45%), and sub-MS galaxies with SFRs ~10-1000 times smaller (~55%). We did not find a well defined class of intermediate, transient objects below the MS, suggesting that the conversion of a massive MS galaxy into a quenched remnant may take a relatively short time (<1 Gyr), though a larger sample should be analyzed in the same way to set precise limits on the quenching timescale. X-ray detected AGNs represent a ~30% fraction of the sample, and are found among both star-forming and quenched galaxies. The morphological analysis revealed that ~50% of our massive objects are bulge-dominated, and almost all MS galaxies with a relevant bulge component host an AGN. We also found sub-MS SFRs in many bulge-dominated systems, providing support to the notion that bulge growth, AGN activity and quenching of star formation are closely related to each other.
The problem of cosmological production of (massless) gravitons is discussed in the framework of an expanding, spatially homogeneous and isotropic FRW type Universe with decaying vacuum energy density ($\Lambda \equiv \Lambda(H(t))$) described by general relativity theory. The gravitational wave equation is established and its time-dependent part has analytically been solved for different epochs in the case of a flat geometry. Unlike the standard $\Lambda$CDM cosmology (no interacting vacuum), we show that massless gravitons can be produced during the radiation era. However, high frequency modes are damped out even faster than in the standard cosmology both in the radiation and matter-vacuum dominated epoch. The formation of the stochastic background of gravitons and the remnant power spectrum generated at different cosmological eras are also explicitly evaluated.
New insights on stellar evolution and stellar interior physics are being made possible by asteroseismology, the study of stars by the observation of their natural, resonant oscillations. Asteroseismology is making significant contributions to our understanding of solar-type stars, in great part due to the exquisite data that have been made available by NASA's Kepler space telescope. Of particular interest is the synergy between asteroseismology and exoplanetary science. Herein I will review recent contributions from asteroseismology to the determination of fundamental properties of Kepler exoplanet-host stars and stress its potential in constraining the spin-orbit alignment of exoplanet systems.
We show Akari data, Herschel data and data from the SCUBA2 camera on JCMT, of molecular clouds. We focus on pre-stellar cores within the clouds. We present Akari data of the L1147-1157 ring in Cepheus and show how the data indicate that the cores are being externally heated. We present SCUBA2 and Herschel data of the Ophiuchus region and show how the environment is also affecting core evolution in this region. We discuss the effects of the magnetic field in the Lupus I region, and how this lends support to a model for the formation and evolution of cores in filamentary molecular clouds.
Interplanetary Coronal Mass Ejections are the manifestation of solar transient eruptions, which can significantly modify the plasma and magnetic conditions in the heliosphere. They are often preceded by a shock, and a magnetic flux rope is detected in situ in a third to half of them. The main aim of this study is to obtain the best quantitative shape for the flux rope axis and for the shock surface from in situ data obtained during spacecraft crossings of these structures. We first compare the orientation of the flux ropes axes and shock normals obtained from independent data analyses of the same events, observed in situ at 1AU from the Sun. Then, we carry out an original statistical analysis of axes/shock normals by deriving the statistical distributions of their orientations. We fit the observed distributions using the distributions derived from several synthetic models describing these shapes. We show that the distributions of axis/shock orientations are very sensitive to their respective shape. One classical model, used to analyze interplanetary imager data, is incompatible with the in situ data. Two other models are introduced, for which the results for axis and shock normals lead to very similar shapes; the fact that the data for MCs and shocks are independent strengthen this result. The model which best fit all the data sets has an ellipsoidal shape with similar aspect ratio values for all the data sets. These derived shapes for the flux rope axis and shock surface have several potential applications. First, these shapes can be used to construct a consistent ICME model. Second, these generic shapes can be used to develop a quantitative model to analyze imager data, as well as constraining the output of numerical simulations of ICMEs. Finally, they will have implications for space weather forecasting, in particular for forecasting the time arrival of ICMEs at the Earth.
Context. he study of prestellar cores is essential to understanding the initial stages of star formation. With $Herschel$ more cold clumps have been detected than ever before. For this study we have selected 21 cold clumps from 20 $Herschel$ fields observed as a follow-up on original $Planck$ detections. We have observed these clumps in $^{13}$CO (1-0), C$^{18}$O (1-0), and N$_2$H$^+$ (1-0) lines. Aims. Our aim is to find out if these cold clumps are prestellar. We have examined to what extent independent analysis of the dust and the molecular lines lead to similar conclusions about the masses of these objects. Methods. We calculate the clump masses and densities from the dust continuum and molecular line observations and compare these to each other and to the virial and Bonnor-Ebert masses calculated for each clump. Finally we examine two of the fields with radiative transfer models to estimate CO abundances. Results. When excitation temperatures could be estimated, the column densities derived from molecular line observations were comparable to those from dust continuum data. The median column density estimates are 4.2$\times 10^{21}$cm$^{-2}$ and 5.5$\times 10^{21}$cm$^{-2}$ for the line and dust emission data, respectively. The calculated abundances, column densities, volume densities, and masses all have large uncertainties and one must be careful when drawing conclusions. Abundance of $^{13}$CO was found in modeling the two clumps in the field G131.65$+$9.75 to be close to the usual value of 10$^{-6}$. The abundance ratio of $^{13}$CO and C$^{18}$O was $\sim$10. Molecular abundances could only be estimated with modeling, relying on dust column density data. Conclusions. The results indicate that most cold clumps, even those with dust color temperatures close to 11 K, are not necessarily prestellar.
The shock interaction and evolution of nova ejecta with a wind from a red giant star in a symbiotic binary system are investigated via three-dimensional hydrodynamics simulations. We specifically model the March 2010 outburst of the symbiotic recurrent nova V407~Cygni from the quiescent phase to its eruption phase. The circumstellar density enhancement due to wind-white dwarf interaction is studied in detail. It is found that the density-enhancement efficiency depends on the ratio of the orbital speed to the red giant wind speed. Unlike another recurrent nova, RS~Ophiuchi, we do not observe a strong disk-like density enhancement, but instead observe an aspherical density distribution with $\sim 20\%$ higher density in the equatorial plane than at the poles. To model the 2010 outburst, we consider several physical parameters, including the red giant mass loss rate, nova eruption energy, and ejecta mass. A detailed study of the shock interaction and evolution reveals that the interaction of shocks with the red giant wind generates strong Rayleigh-Taylor instabilities. In addition, the presence of the companion and circumstellar density enhancement greatly alter the shock evolution during the nova phase. The ejecta speed after sweeping out most of the circumstellar medium decreases to $\sim 100-300$ km-s$^{-1}$, depending on model, which is consistent with the observed extended redward emission in [N~II] lines in April 2011.
The CHIME Pathfinder is a new interferometric radio telescope that uses a hybrid FPGA/GPU FX correlator. The GPU-based X-engine of this correlator processes over 819 Gb/s of 4+4-bit complex astronomical data from N=256 inputs across a 400 MHz radio band. A software framework is presented to manage this real-time data flow, which allows each of 16 processing servers to handle 51.2 Gb/s of astronomical data, plus 8 Gb/s of ancillary data. Each server receives data in the form of UDP packets from an FPGA F-engine over the eight 10 GbE links, combines data from these packets into large (32MB-256MB) buffered frames, and transfers them to multiple GPU co-processors for correlation. The results from the GPUs are combined and normalized, then transmitted to a collection server, where they are merged into a single file. Aggressive optimizations enable each server to handle this high rate of data; allowing the efficient correlation of 25 MHz of radio bandwidth per server. The solution scales well to larger values of N by adding additional servers.
We report the presence of scattered light from dust grains located in the giant planet formation region of the circumbinary disk orbiting the ~20-Myr-old close (~0.045 AU separation) binary system V4046 Sgr AB based on observations with the new Gemini Planet Imager (GPI) instrument. These GPI images probe to within ~7 AU of the central binary with linear spatial resolution of ~3 AU, and are thereby capable of revealing dust disk structure within a region corresponding to the giant planets in our solar system. The GPI imaging reveals a relatively narrow (FWHM ~10 AU) ring of polarized near-infrared flux whose brightness peaks at ~14 AU. This ~14 AU radius ring is surrounded by a fainter outer halo of scattered light extending to ~45 AU, which coincides with previously detected mm-wave thermal dust emission. The presence of small grains that efficiently scatter starlight well inside the mm-wavelength disk cavity supports current models of planet formation that suggest planet-disk interactions can generate pressure traps that impose strong radial variations in the particle size distribution throughout the disk.
The aim of this work is to determine the total integrated flux of cosmic radiation which a commercial aircraft is exposed to along specific flight trajectories. To study the radiation background during a flight and its modulation by effects such as altitude, latitude, exposure time and transient magnetospheric events, we perform simulations based on Magnetocosmics and CORSIKA codes, the former designed to calculate the geomagnetic effects on cosmic rays propagation and the latter allows us to simulate the development of extended air showers in the atmosphere. In this first work, by considering the total flux of cosmic rays from 5 GeV to 1 PeV, we obtained the expected integrated flux of secondary particles on board of a commercial airplane during the Bogot\'a-Buenos Aires trip by point-to-point numerical integration.
We present the design and implementation of a custom GPU-based compute cluster that provides the correlation X-engine of the CHIME Pathfinder radio telescope. It is among the largest such systems in operation, correlating 32,896 baselines (256 inputs) over 400MHz of radio bandwidth. Making heavy use of consumer-grade parts and a custom software stack, the system was developed at a small fraction of the cost of comparable installations. Unlike existing GPU backends, this system is built around OpenCL kernels running on consumer-level AMD GPUs, taking advantage of low-cost hardware and leveraging packed integer operations to double algorithmic efficiency. The system achieves the required 105TOPS in a 10kW power envelope, making it among the most power-efficient X-engines in use today.
Interferometric radio telescopes often rely on computationally expensive O(N^2) correlation calculations; fortunately these computations map well to massively parallel accelerators such as low-cost GPUs. This paper describes the OpenCL kernels developed for the GPU based X-engine of a new hybrid FX correlator. Channelized data from the F-engine is supplied to the GPUs as 4-bit, offset-encoded real and imaginary integers. Because of the low bit width of the data, two values may be packed into a 32-bit register, allowing multiplication and addition of more than one value with a single fused multiply-add instruction. With this data and calculation packing scheme, as many as 5.6 effective tera-operations per second (TOPS) can be executed on a 4.3 TOPS GPU. The kernel design allows correlations to scale to large numbers of input elements, limited only by maximum buffer sizes on the GPU. This code is currently working on-sky with the CHIME Pathfinder Correlator in BC, Canada.
The recent discovery of ten new dwarf galaxy candidates by the Dark Energy Survey (DES) and the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) could increase the Fermi Gamma-Ray Space Telescope's sensitivity to annihilating dark matter particles, potentially enabling a definitive test of the dark matter interpretation of the long-standing Galactic Center gamma-ray excess. In this paper, we compare the previous analyses of Fermi data from the directions of the new dwarf candidates (including the relatively nearby Reticulum II) and perform our own analysis, with the goal of establishing the statistical significance of any gamma-ray signal from these sources. We confirm the presence of an excess from Reticulum II, with a spectral shape that is compatible with the Galactic Center signal. The significance of this emission is greater than that observed from 99.84% of randomly chosen high-latitude blank-sky locations, corresponding to a local detection significance of 3.2 sigma. We improve upon the standard blank-sky calibration approach through the use of multi-wavelength catalogs, which allow us to avoid regions that are likely to contain unresolved gamma-ray sources.
An excess in gamma-rays from the galactic center observed by the Fermi Large Area Telescope has been proposed as a possible signal of dark matter annihilation. Recently, the Fermi collaboration showed that systematic errors broaden the range of spectral shapes for this excess. We demonstrate fits to this range for (1) flavor-violating annihilations to top-charm pairs and (2) annihilations to on-shell bosonic mediators which decay to Standard Model quarks in a boosted frame. Annihilation of 40 - 100 GeV DM to pairs of spin-1 mediators provide a good fit to the Fermi-LAT spectrum with a normalization consistent with a thermal relic. Top-charm modes and annihilation to three pseudoscalar mediators can fit the spectral shape but typically require non-thermal annihilation cross sections.
By applying fractional calculus to the equation proposed by M. Planck in 1900, we obtain a new blackbody radiation law described by a Mittag-Leffler (ML) function. We have analyzed NASA COBE data by means of a non-extensive formula with a parameter $(q-1)$, a formula proposed by Ertik et al. with a fractional parameter $(\alpha-1)$, and our new formula including a parameter $(p-1)$, as well as the Bose-Einstein distribution with a dimensionless chemical potential $\mu$. It can be said that one role of the fractional parameter $(p-1)$ is almost the same as that of chemical potential $(\mu)$ as well as that of the parameter $(q-1)$ in the non-extensive approach.
We consider the Damour-Dyson analysis of the sensitivity of neutron resonance energies $E_i$ to changes in the fine structure constant $\alpha$. We point out that, with more appropriate choices of nuclear parameters, their result for ${}^{150}$Sm is increased by a factor of 2.5. We go on to identify and compute excitation, Coulomb and deformation corrections. To this end, we use deformed Fermi density distributions fitted to the output of HF+BCS calculations, the energetics of the surface diffuseness of nuclei, and thermal properties of their deformation; we also invoke the eigenstate thermalization hypothesis, performing the requisite microcanonical averages with phenomenological level densities which include the effect of increased surface diffuseness. We find that the corrections diminish the revised ${}^{150}$Sm sensitivity but not by more than 25\%. More precisely, we establish that $\alpha dE_i / d\alpha < (-2.2\pm 0.3)\,\text{MeV}$ (similar inequalities are also obtained for ${}^{156}$Gd and ${}^{158}$Gd). Subject to a weak and testable restriction on the change in $m_q/\Lambda$ (relative to the change in $\alpha$) since the time when the Oklo reactors were active ($m_q$ is the average of the $\text{u}$ and $\text{d}$ current quark masses, and $\Lambda$ is the mass scale of quantum chromodynamics), we deduce that $|\alpha_{\text{Oklo}} - \alpha_{\text{now}}| < 1.3\times 10^{-8} \alpha_{\text{now}}$. This bound is comparable to existing Oklo-based limits, but has a stronger theoretical basis which resolves uncertainties that plagued earlier treatments.
In 'hilltop inflation', inflation takes place when the inflaton field slowly rolls from close to a maximum of its potential (i.e. the 'hilltop') towards its minimum. When the inflaton potential is associated with a phase transition, possible topological defects produced during this phase transition, such as domain walls, are efficiently diluted during inflation. It is typically assumed that they also do not reform after inflation, i.e. that the inflaton field stays on its side of the 'hill', finally performing damped oscillations around the minimum of the potential. In this paper we study the linear and the non-linear phases of preheating after hilltop inflation. We find that the fluctuations of the inflaton field during the tachyonic oscillation phase grow strong enough to allow the inflaton field to form regions in position space where it crosses 'over the top of the hill' towards the 'wrong vacuum'. We investigate the formation and behaviour of these overshooting regions using lattice simulations: Rather than durable domain walls, these regions form oscillon-like structures (i.e. localized bubbles that oscillate between the two vacua) which should be included in a careful study of preheating in hilltop inflation.
We describe a class of modified gravity theories that deform general
relativity in a way that breaks time reversal invariance and, very mildly,
locality. The algebra of constraints, local physical degrees of freedom, and
their linearized equations of motion, are unchanged, yet observable effects may
be present on cosmological scales, which have implications for the early
history of the universe.
This is achieved in the Hamiltonian framework, in a way that requires the
constant mean curvature gauge conditions and is, hence, inspired by shape
dynamics.
We attempt to simultaneously explain the recently observed 3.55 keV X-ray line in the analysis of XMM-Newton telescope data and the galactic center gamma ray excess observed by the Fermi gamma ray space telescope within an abelian gauge extension of standard model. We consider a two component dark matter scenario with a mass difference 3.55 keV such that the heavier one can decay into the lighter one and a photon with energy 3.55 keV. The lighter dark matter candidate is protected from decaying into the standard model particles by a remnant $Z_2$ symmetry into which the abelian gauge symmetry gets spontaneously broken. If the mass of the dark matter particle is chosen to be within $31-40$ GeV, then this model can also explain the galactic center gamma ray excess if the dark matter annihilation into $b\bar{b}$ pairs has a cross section of $\langle \sigma v \rangle \simeq (1.4-2.0) \times 10^{-26} \; \text{cm}^3/\text{s}$. We constrain the model from the requirement of producing correct dark matter relic density, 3.55 keV X-ray line flux and galactic center gamma ray excess. We also impose the bounds coming from dark matter direct detection experiments as well as collider limits on additional gauge boson mass and coupling. We also briefly discuss how this model can give rise to sub-eV neutrino masses at tree level as well as one loop level while keeping the dark matter mass at few tens of GeV. We also show the natural origin of keV mass splitting between two electroweak scale dark matter particles at one loop level in this model.
There is a growing appreciation that the environmental conditions that we call space weather impact the technological infrastructure that powers the coupled economies around the world. With that comes the need to better shield society against space weather by improving forecasts, environmental specifications, and infrastructure design. [...] advanced understanding of space weather requires a coordinated international approach to effectively provide awareness of the processes within the Sun-Earth system through observation-driven models. This roadmap prioritizes the scientific focus areas and research infrastructure that are needed to significantly advance our understanding of space weather of all intensities and of its implications for society. Advancement of the existing system observatory through the addition of small to moderate state-of-the-art capabilities designed to fill observational gaps will enable significant advances. Such a strategy requires urgent action: key instrumentation needs to be sustained, and action needs to be taken before core capabilities are lost in the aging ensemble. We recommend advances through priority focus (1) on observation-based modeling throughout the Sun-Earth system, (2) on forecasts more than 12 hrs ahead of the magnetic structure of incoming coronal mass ejections, (3) on understanding the geospace response to variable solar-wind stresses that lead to intense geomagnetically-induced currents and ionospheric and radiation storms, and (4) on developing a comprehensive specification of space climate, including the characterization of extreme space storms to guide resilient and robust engineering of technological infrastructures. The roadmap clusters its implementation recommendations by formulating three action pathways, and outlines needed instrumentation and research programs and infrastructure for each of these. [...]
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
We have discovered that SDSSJ105213.51+442255.7 (T0.5$\pm$1.0) is a binary in Keck laser guide star adaptive optics imaging, displaying a large J-to-K-band flux reversal ($\Delta$J = -0.45$\pm$0.09 mag, $\Delta$K = 0.52$\pm$0.05 mag). We determine a total dynamical mass from Keck orbital monitoring (88$\pm$5 $M_{\rm Jup}$) and a mass ratio by measuring the photocenter orbit from CFHT/WIRCam absolute astrometry ($M_B/M_A$ = 0.78$\pm$0.07). Combining these provides the first individual dynamical masses for any field L or T dwarfs, 49$\pm$3 $M_{\rm Jup}$ for the L6.5$\pm$1.5 primary and 39$\pm$3 $M_{\rm Jup}$ for the T1.5$\pm$1.0 secondary. Such a low mass ratio for a nearly equal luminosity binary implies a shallow mass$-$luminosity relation over the L/T transition ($\Delta$log$L_{\rm bol}$/$\Delta$log$M = 0.6^{+0.6}_{-0.8}$). This provides the first observational support that cloud dispersal plays a significant role in the luminosity evolution of substellar objects. Fully cloudy models fail our coevality test for this binary, giving ages for the two components that disagree by 0.2 dex (2.0$\sigma$). In contrast, our observed masses and luminosities can be reproduced at a single age by "hybrid" evolutionary tracks where a smooth change from a cloudy to cloudless photosphere around 1300 K causes slowing of luminosity evolution. Remarkably, such models also match our observed JHK flux ratios and colors well. Overall, it seems that the distinguishing features SDSSJ1052+4422AB, like a J-band flux reversal and high-amplitude variability, are normal for a field L/T binary caught during the process of cloud dispersal, given that the age (1.11$^{+0.17}_{-0.20}$ Gyr) and surface gravity (log$g$ = 5.0$-$5.2) of SDSSJ1052+4422AB are typical for field ultracool dwarfs.
Weak gravitational lensing is a powerful cosmological probe, with non--Gaussian features potentially containing the majority of the information. We examine constraints on the parameter triplet $(\Omega_m,w,\sigma_8)$ from non-Gaussian features of the weak lensing convergence field, including a set of moments (up to $4^{\rm th}$ order) and Minkowski functionals, using publicly available data from the 154deg$^2$ CFHTLenS survey. We utilize a suite of ray--tracing N-body simulations spanning 91 points in $(\Omega_m,w,\sigma_8)$ parameter space, replicating the galaxy sky positions, redshifts and shape noise in the CFHTLenS catalogs. We then build an emulator that interpolates the simulated descriptors as a function of $(\Omega_m,w,\sigma_8)$, and use it to compute the likelihood function and parameter constraints. We employ a principal component analysis to reduce dimensionality and to help stabilize the constraints with respect to the number of bins used to construct each statistic. Using the full set of statistics, we find $\Sigma_8\equiv\sigma_8(\Omega_m/0.27)^{0.55}=0.75\pm0.04$ (68% C.L.), in agreement with previous values. We find that constraints on the $(\Omega_m,\sigma_8)$ doublet from the Minkowski functionals suffer a strong bias. However, high-order moments break the $(\Omega_m,\sigma_8)$ degeneracy and provide a tight constraint on these parameters with no apparent bias. The main contribution comes from quartic moments of derivatives.
We present the discovery of a new dwarf galaxy, Hydra II, found serendipitously within the data from the ongoing Survey of the MAgellanic Stellar History (SMASH) conducted with the Dark Energy Camera on the Blanco 4m Telescope. The new satellite is compact (r_h = 68 +/- 11 pc) and faint (M_V = -4.8 +/- 0.3), but well within the realm of dwarf galaxies. The stellar distribution of HydraII in the color-magnitude diagram is well-described by a metal-poor ([Fe/H] = -2.2) and old (13 Gyr) isochrone and shows a distinct blue horizontal branch, some possible red clump stars, and faint stars that are suggestive of blue stragglers. At a heliocentric distance of 134 +/- 10 kpc, Hydra II is located in a region of the Galactic halo that models have suggested may host material from the leading arm of the Magellanic Stream. A comparison with N-body simulations hints that the new dwarf galaxy could be or could have been a satellite of the Magellanic Clouds.
The next generation gravitational waves (GW) detectors are most sensitive to GW emitted by compact (neutron star/black hole) binary mergers. If one of those is a neutron star the merger will also emit electromagnetic radiation via three possible channels: Gamma-ray bursts and their (possibly orphan) afterglows (Eichler et al. 1989), Li-Paczynski Macronovae (Li & Paczynski 1998) and radio flares (Nakar & Piran 2011). This accompanying electromagnetic radiation is vitally important in confirming the GW detections (Kochanek & Piran 1993). It could also reveal a wealth of information regarding the merger and will open a window towards multi-messenger astronomy. Identifying and characterizing these counterparts is therefore of utmost importance. In this work we explore late time radio flares emitted by the dynamically ejected outflows. We build upon previous work and consider the effect of the outflow's non-trivial geometry. Using an approximate method we estimate the radio light-curves for several ejected matter distributions obtained in numerical simulations. Our method provides an upper limit to the effect of non-sphericity. Together with the spherical estimates the resulting light curves bound the actual signal. We find that while non-spherical geometries can in principle lead to an enhanced emission, in most cases they result in an increase in the timescale compared with a corresponding spherical configuration. This would weaken somewhat these signals and might decrease the detection prospects.
We present the first pointed X-ray observations of ten candidate fossil galaxy groups and clusters. With these Suzaku observations, we determine global temperatures and bolometric X-ray luminosities of the intracluster medium (ICM) out to $r_{500}$ for six systems in our sample. The remaining four systems show signs of significant contamination from non-ICM sources. For the six objects with successfully determined $r_{500}$ properties, we measure global temperatures between $2.8 \ \mathrm{keV} \leq T_{\mathrm{X}} \leq 5.3 \ \mathrm{keV}$, bolometric X-ray luminosities of $0.6 \times 10^{44} \ \mathrm{ergs} \ \mathrm{s}^{-1} \leq L_{\mathrm{X,bol}} \leq 7.2\times 10^{44} \ \mathrm{ergs} \ \mathrm{s}^{-1}$, and estimate masses, as derived from $T_{\mathrm{X}}$, of $M_{500} \gtrsim 10^{14} \ \mathrm{M}_{\odot}$. Scaling relations are constructed for an assembled sample of fossil and non-fossil systems using global X-ray luminosities, temperatures, optical luminosities, and velocity dispersions. The fit of the scaling relations for fossil systems is found to be consistent with the fit of the relations for normal groups and clusters. We find fossil systems have global ICM X-ray properties similar to those of comparable mass non-fossil systems.
At z>1, the distinction between merging and 'normal' star-forming galaxies based on single band morphology is often hampered by the presence of large clumps which result in a disturbed, merger-like appearance even in rotationally supported disks. In this paper we discuss how a classification based on canonical, non-parametric structural indices measured on resolved stellar mass maps, rather than on single-band images, reduces the misclassification of clumpy but not merging galaxies. We calibrate the mass-based selection of mergers using the MIRAGE hydrodynamical numerical simulations of isolated and merging galaxies which span a stellar mass range of $10^{9.8}$-$10^{10.6}M_{sun}$ and merger ratios between 1:1-1:6.3. These simulations are processed to reproduce the typical depth and spatial resolution of observed HUDF data. We test our approach on a sample of real z~2 galaxies with kinematic classification into disks or mergers and on ~100 galaxies in the HUDF field with photometric/spectroscopic redshift between 1.5$\leqslant z \leqslant$3 and $M>10^{9.4}M_{sun}$. We find that a combination of the asymmetry $A_{MASS}$ and $M_{20, MASS}$ indices measured on the stellar mass maps can efficiently identify real (major) mergers with $\lesssim$20% contamination from clumpy disks in the merger sample. This mass-based classification cannot be reproduced in star-forming galaxies by $H-$band measurements alone, which instead result in a contamination from clumpy galaxies that can be as high as 50%. Moreover, we find that the mass-based classification always results in a lower contamination from clumpy galaxies than an $H-$band classification, regardless of the depth of the imaging used (e.g., CANDELS vs. HUDF).
NGC 7129 is a bright reflection nebula located in the molecular cloud complex near l=105.4, b=+9.9, about 1.15 kpc distant. Embedded within the reflection nebula is a young cluster dominated by a compact grouping of four early-type stars: BD+65 1638 (B3V), BD+65 1637 (B3e), SVS 13 (B5e), and LkH-alpha 234 (B8e). About 80 H-alpha emission sources brighter than V~23 are identified in the region, many of which are presumably T Tauri star members of the cluster. We also present deep (V~23), optical (VRI) photometry of a field centered on the reflection nebula and spectral types for more than 130 sources determined from low dispersion, optical spectroscopy. The narrow pre-main sequence evident in the color-magnitude diagram suggests that star formation was rapid and coeval. A median age of about 1.8 Myr is inferred for the H-alpha and literature-identified X-ray emission sources having established spectral types, using pre-main sequence evolutionary models. Our interpretation of the structure of the molecular cloud and the distribution of young stellar objects is that BD+65 1638 is primarily responsible for evacuating the blister-like cavity within the molecular cloud. LkH-alpha 234 and several embedded sources evident in near infrared adaptive optics imaging have formed recently within the ridge of compressed molecular gas. The compact cluster of low-mass stars formed concurrently with the early-type members, concentrated within a central radius of ~0.7 pc. Star formation is simultaneously occurring in a semi-circular arc some ~3 pc in radius that outlines remaining dense regions of molecular gas. High dispersion, optical spectra are presented for BD+65 1638, BD+65 1637, SVS 13, LkH-alpha 234, and V350 Cep. These spectra are discussed in the context of the circumstellar environments inferred for these stars.
We investigated the copper abundances for $64$ late-type stars in the Galactic disk and halo with effective temperatures from $5400$ K to $6700$ K and [Fe/H] from $-1.88$ to $-0.17$. For the first time, the copper abundances are derived using both local thermodynamic equilibrium (LTE) and non-local thermodynamic equilibrium (non-LTE) calculations. High resolution ($R > 40,000$), high signal-to-noise ratio ($S/N > 100$) spectra from the FOCES spectrograph are used. The atmospheric models are calculated based on the MAFAGS opacity sampling code. All the abundances are derived using the spectrum synthesis methods. Our results indicate that the non-LTE effects of copper are important for metal-poor stars, showing a departure of $\sim 0.17$ dex at the metallicity $\sim -1.5$. We also find that the copper abundances derived from non-LTE calculations are enhanced compared with those from LTE. The enhancements show clear dependence on the metallicity, which gradually increase with decreasing [Fe/H] for our program stars, leading to a flatter distribution of [Cu/Fe] with [Fe/H] than previous work. There is a hint that the thick- and thin-disk stars have different behaviors in [Cu/Fe], and a bending for disk stars may exist.
In the summer of 2012, during a Pulsar Search Collaboratory workshop, two high-school students discovered J1930$-$1852, a pulsar in a double neutron star (DNS) system. Most DNS systems are characterized by short orbital periods, rapid spin periods and eccentric orbits. However, J1930$-$1852 has the longest spin period ($P_{\rm spin}\sim$185 ms) and orbital period ($P_{\rm b}\sim$45 days) yet measured among known, recycled pulsars in DNS systems, implying a shorter than average and/or inefficient recycling period before its companion went supernova. We measure the relativistic advance of periastron for J1930$-$1852, $\dot{\omega}=0.00078$(4) deg/yr, which implies a total mass (M$_{\rm{tot}}=2.59$(4) M$_{\odot}$) consistent with other DNS systems. The $2\sigma$ constraints on M$_{\rm{tot}}$ place limits on the pulsar and companion masses ($m_{\rm p}<1.32$ M$_{\odot}$ and $m_{\rm c}>1.30$ M$_{\odot}$ respectively). J1930$-$1852's spin and orbital parameters challenge current DNS population models and make J1930$-$1852 an important system for further investigation.
We report accurate measurements of the physical properties (mass, radius, temperature) of components of the G+M eclipsing binary V530 Ori. The M-type secondary shows a larger radius and a cooler temperature than predicted by standard stellar evolution models, as has been found for many other low-mass stars and ascribed to the effects of magnetic activity and/or spots. We show that models from the Dartmouth series that incorporate magnetic fields are able to match the observations with plausible field strengths of 1-2 kG, consistent with a rough estimate we derive for that star.
Magnetars are neutron stars in which a strong magnetic field is the main energy source. About two dozens of magnetars, plus several candidates, are currently known in our Galaxy and in the Magellanic Clouds. They appear as highly variable X-ray sources and, in some cases, also as radio and/or optical pulsars. Their spin periods (2-12 s) and spin-down rates (~10^{-13}-10^{-10} s/s) indicate external dipole fields of ~10^{13-15} G, and there is evidence that even stronger magnetic fields are present inside the star and in non-dipolar magnetospheric components. Here we review the observed properties of the persistent emission from magnetars, discuss the main models proposed to explain the origin of their magnetic field and present recent developments in the study of their evolution and connection with other classes of neutron stars.
The launch of the gamma-ray telescope Fermi Large Area Telescope (Fermi-LAT) started a pivotal period in indirect detection of dark matter. By outperforming expectations, for the first time a robust and stringent test of the paradigm of weakly interacting massive particles (WIMPs) is within reach. In this paper, we discuss astrophysical targets for WIMP detection and the challenges they present, review the analysis tools which have been employed to tackle these challenges, and summarize the status of constraints on and the claimed detections in the WIMP parameter space. Methods and results will be discussed in comparison to Imaging Air Cherenkov Telescopes. We also provide an outlook on short term and longer term developments.
AA Tau, a classical T Tauri star in the Taurus cloud, has been the subject of intensive photometric monitoring for more than two decades due to its quasi-cyclic variation in optical brightness. Beginning in 2011, AA Tau showed another peculiar variation -- its median optical though near-IR flux dimmed significantly, a drop consistent with a 4-mag increase in visual extinction. It has stayed in the faint state since.Here we present 4.7um CO rovibrational spectra of AA Tau over eight epochs, covering an eleven-year time span, that reveal enhanced 12CO and 13CO absorption features in the $J_{\rm low}\leqslant$13 transitions after the dimming. These newly appeared absorptions require molecular gas along the line of sight with T~500 K and a column density of log (N12CO)~18.5 cm^{-2}, with line centers that show a constant 6 km s$^{-1}$ redshift. The properties of the molecular gas confirm an origin in the circumstellar material. We suggest that the dimming and absorption are caused by gas and dust lifted to large heights by a magnetic buoyancy instability. This material is now propagating inward, and on reaching the star within a few years will be observed as an accretion outburst.
We report on the weak lensing detection of a filament between two galaxy clusters at $z=0.55$, CL0015.9+1609 and RX J0018.3+1618. We conduct weak lensing analysis of deep multi-band Subaru/Suprime-Cam images with $Lensfit$. The weak lensing signals from the filament are contaminated by signals from the adjacent massive clusters and we statistically subtract the cluster component using two different methods. Both methods yield consistent shear profiles on the filament with $\gtrsim2\sigma$ significance and the average surface mass density of the filament is $<\Sigma>=(3.20\pm0.10)\times10^{14}h$ M$_\odot$Mpc$^{-2}$, which is in broad agreement with previous studies. On-going surveys such as Hyper Suprime-Cam will identify more filaments, which will serve as a new probe of structure formation in the Universe.
The use of type Ic Super Luminous Supernovae (SLSN Ic) to examine the cosmological expansion introduces a new standard ruler with which to test theoretical models. The sample suitable for this kind of work now includes 11 SLSNe Ic, which have thus far been used solely in tests involving $\Lambda$CDM. In this paper, we broaden the base of support for this new, important cosmic probe by using these observations to carry out a one-on-one comparison between the $R_{\rm h}=ct$ and $\Lambda$CDM cosmologies. We individually optimize the parameters in each cosmological model by minimizing the $\chi^{2}$ statistic. We also carry out Monte Carlo simulations based on these current SLSN Ic measurements to estimate how large the sample would have to be in order to rule out either model at a $\sim 99.7\%$ confidence level. The currently available sample indicates a likelihood of $\sim$$70-80\%$ that the $R_{\rm h}=ct$ Universe is the correct cosmology versus $\sim$$20-30\%$ for the standard model. These results are suggestive, though not yet compelling, given the current limited number of SLSNe Ic. We find that if the real cosmology is $\Lambda$CDM, a sample of $\sim$$240$ SLSNe Ic would be sufficient to rule out $R_{\rm h}=ct$ at this level of confidence, while $\sim$$480$ SLSNe Ic would be required to rule out $\Lambda$CDM if the real Universe is instead $R_{\rm h}=ct$. This difference in required sample size reflects the greater number of free parameters available to fit the data with $\Lambda$CDM. If such SLSNe Ic are commonly detected in the future, they could be a powerful tool for constraining the dark-energy equation of state in $\Lambda$CDM, and differentiating between this model and the $R_{\rm h}=ct$ Universe.
A bright fireball with an absolute magnitude of -10.5 $\pm$ 0.5 was recorded over the South of Spain on August 15, 2012. The atmospheric trajectory, radiant and heliocentric orbit of this event are calculated. These data show that the parent meteoroid belonged to the \kappa-Cygnid meteoroid stream. The emission spectrum of this bolide, which was obtained in the wavelength range between 350 and 800 nm, suggests a chondritic nature for the progenitor meteoroid. Besides, the spectrum of the meteoric afterglow was also recorded for about 0.7 seconds. The evolution with time of the intensity of the main emission lines identified in this signal is discussed.
We present a far-infrared all-sky atlas from a sensitive all-sky survey using the Japanese $AKARI$ satellite. The survey covers $> 99$% of the sky in four photometric bands centred at 65 $\mu$m, 90 $\mu$m, 140 $\mu$m, and 160 $\mu$m with spatial resolutions ranging from 1 to 1.5 arcmin. These data provide crucial information for the investigation and characterisation of the properties of dusty material in the Interstellar Medium (ISM), since significant portion of its energy is emitted between $\sim$50 and 200 $\mu$m. The large-scale distribution of interstellar clouds, their thermal dust temperatures and column densities, can be investigated with the improved spatial resolution compared to earlier all-sky survey observations. In addition to the point source distribution, the large-scale distribution of ISM cirrus emission, and its filamentary structure, are well traced. We have made the first public release of the full-sky data to provide a legacy data set for use by the astronomical community.
We present an initial analysis of the properties of the all-sky image obtained by the Far-Infrared Surveyor (FIS) onboard the AKARI satellite, at 65~$\mu$m (N60), 90~$\mu$m (WIDE-S), 140~$\mu$m (WIDE-L),and 160~$\mu$m (N160). Absolute flux calibration was determined by comparing the data with the COBE/DIRBE data sets, and the intensity range was as wide as from a few MJy~sr$^{-1}$ to $>$1~GJy~sr$^{-1}$. The uncertainties are considered to be the standard deviations with respect to the DIRBE data, and they are less than 10\% for intensities above 10, 3, 25, and 26~MJy~sr$^{-1}$ at the N60, WIDE-S, WIDE-L, and N160 bands, respectively. The characteristics of point sources in the image were also determined by stacking maps centred on photometric standard stars. The full width at half maxima of the point spread functions (PSFs) were 63$"$, 78$"$, and 88$"$ at the N60, WIDE-S, and WIDE-L bands, respectively. The PSF at the N160 band was not obtained due to the sensitivity, but it is thought to be the same as that of the WIDE-L one.
G349.7$+$0.2 is a supernova remnant (SNR) expanding in a dense medium of molecular clouds and interacting with clumps of molecular material emitting gamma rays. We analyzed the gamma-ray data of Large Area Telescope on board Fermi Gamma Ray Space Telescope and detected G349.7$+$0.2 in the energy range of 0.2$-$300 GeV with a significance of $\sim$13$\sigma$ showing no extended morphology. Modeling of the gamma-ray spectrum revealed that the GeV gamma-ray emission dominantly originates from the decay of neutral pions, where the protons follow a broken power-law distribution with a spectral break at $\sim$12 GeV. To search for features of radiative recombination continua in the eastern and western regions of the remnant, we analyzed the Suzaku data of G349.7$+$0.2 and found no evidence for overionized plasma. In this paper we discuss possible scenarios to explain the hadronic gamma-ray emission in G349.7$+$0.2 and the mixed morphology nature of this SNR.
We present an analysis of the physical conditions in an extreme molecular cloud in the Antennae merging galaxies. This cloud has properties consistant with those required to form a globular cluster. We have obtained ALMA CO and 870$\mu$m observations of the Antennae galaxy system with $\sim 0".5$ resolution. This cloud stands out in the data with a radius of $\lesssim 24$~pc and mass of $>5\times 10^6$~M$_\odot$. The cloud appears capable of forming a globular cluster, but the lack of associated thermal radio emission indicates that star formation has not yet altered the environment. The lack of thermal radio emission places the cloud in an early stage of evolution, which we expect to be short-lived ($\lesssim 1$~Myr) and thus rare. Given its mass and kinetic energy, for the cloud to be confined (as its appearance strongly suggests) it must be subject to an external pressure of P/$k_B \gtrsim 10^8$~K~cm$^{-3}$ -- 10,000 times higher than typical interstellar pressure. This would support theories that high pressures are required to form globular clusters and may explain why extreme environments like the Antennae are preferred environments for generating such objects. Given the cloud temperature of $\sim 25$~K, the internal pressure must be dominated by non-thermal processes, most likely turbulence. We expect the molecular cloud to collapse and begin star formation in $\lesssim 1$~Myr.
We consider a model of the early universe which consists of two scalar fields: the inflaton, and a second field which drives the stabilisation of the Planck mass (or gravitational constant). We show that the non-minimal coupling of this second field to the Ricci scalar sources a non-adiabatic pressure perturbation. By performing a fully numerical calculation we find, in turn, that this boosts the amplitude of the primordial power spectrum after inflation.
The ratio of deuterium to hydrogen (D/H ratio) of Solar System bodies is an important clue to their formation histories. Here we fit a Neptunian atmospheric model to Gemini Near Infrared Spectrograph (GNIRS) high spectral resolution observations and determine the D/H ratio in methane absorption in the infrared H-band ($\sim$ 1.6 {\mu}m). The model was derived using our radiative transfer software VSTAR (Versatile Software for the Transfer of Atmospheric Radiation) and atmospheric fitting software ATMOF (ATMOspheric Fitting). The methane line list used for this work has only become available in the last few years, enabling a refinement of earlier estimates. We identify a bright region on the planetary disc and find it to correspond to an optically thick lower cloud. Our preliminary determination of CH$_{\rm 3}$D/CH$_{\rm 4}$ is 3.0$\times10^{-4}$, which is in line with the recent determination of Irwin et al. (2014) of 3.0$^{+1.0}_{-0.9}\sim\times10^{-4}$, made using the same model parameters and line list but different observational data. This supports evidence that the proto-solar ice D/H ratio of Neptune is much less than that of the comets, and suggests Neptune formed inside its present orbit.
We present the deepest optical photometry for any dwarf elliptical (dE) galaxy based on Hubble Space Telescope ACS observations of the Local Group dE galaxies NGC 147 and NGC 185. The resulting F606W and F814W color-magnitude diagrams are the first to reach below the main sequence turnoff in a dE galaxy, allowing us to determine full star formation histories in these systems. The ACS fields are located ~1.5 effective radii from the galaxy center to avoid photometric crowding. While our ACS pointings in both dEs show unambiguous evidence for old and intermediate age stars, the mean age in NGC 147 is ~ 4 Gyr younger as compared to NGC 185. In NGC 147, only 40% of stars were in place 12.5 Gyrs ago (z~5), with the bulk of the remaining stellar population forming between 5 to 7 Gyr. In contrast, 70% of stars were formed in NGC 185 field more than12.5 Gyr ago with the majority of the remaining population forming between 8 to 10 Gyr. Star formation ceased in both ACS fields at least 3 Gyr ago. Previous observations in the central regions of NGC 185 show evidence for star formation as recent as 100 Myr ago and a strong metallicity gradient with radius. We suggest that the orbit of NGC 185 has a larger pericenter as compared to NGC 147, allowing it to preserve radial gradients and maintain a small central reservoir of recycled gas. We interpret the inferred differences in star formation histories to imply an earlier infall time into the M31 environment for NGC 185 as compared to NGC 147.
We have designed, fabricated, and measured a 5-channel prototype spectrometer pixel operating in the WR10 band to demonstrate a novel moderate-resolution (R=f/{\Delta}f~100), multi-pixel, broadband, spectrometer concept for mm and submm-wave astronomy. Our design implements a transmission line filter bank using waveguide resonant cavities as a series of narrow-band filters, each coupled to an aluminum kinetic inductance detector (KID). This technology has the potential to perform the next generation of spectroscopic observations needed to drastically improve our understanding of the epoch of reionization (EoR), star formation, and large-scale structure of the universe. We present our design concept, results from measurements on our prototype device, and the latest progress on our efforts to develop a 4-pixel demonstrator instrument operating in the 130-250 GHz band.
A pulsar-like compact star is the rump left behind after a supernova where normal baryonic matter is intensely compressed by gravity, but the real state of such compressed baryonic matter is still not well understood because of the non-perturbative nature of the fundamental color interaction. We argue that pulsars could be of condensed matter of quark clusters, i.e., "quark-cluster stars" which distinguish from conventional neutron and quark stars. In comparison with 2-flavour symmetric micro-nuclei, a quark-cluster star could simply be considered as a macro-nucleus with 3-flavour symmetry. New research achievements both theoretical and observational are briefly presented.
The Spitzer Survey of Stellar Structure in Galaxies (S$^4$G, Sheth et. al. 2010) is a deep 3.6 and 4.5 $\mu$m imaging survey of 2352 nearby ($< 40$ Mpc) galaxies. We describe the S$^4$G data analysis pipeline 4, which is dedicated to 2-dimensional structural surface brightness decompositions of 3.6 $\mu$m images, using GALFIT3.0 \citep{peng2010}. Besides automatic 1-component S\'ersic fits, and 2-component S\'ersic bulge + exponential disk fits, we present human supervised multi-component decompositions, which include, when judged appropriate, a central point source, bulge, disk, and bar components. Comparison of the fitted parameters indicates that multi-component models are needed to obtain reliable estimates for the bulge S\'ersic index and bulge-to-total light ratio ($B/T$), confirming earlier results \citep{laurikainen2007, gadotti2008, weinzirl2009}. In this first paper, we describe the preparations of input data done for decompositions, give examples of our decomposition strategy, and describe the data products released via IRSA and via our web page ({\bf \tt www.oulu.fi/astronomy/S4G\_PIPELINE4/MAIN}). These products include all the input data and decomposition files in electronic form, making it easy to extend the decompositions to suit specific science purposes. We also provide our IDL-based visualization tools (GALFIDL) developed for displaying/running GALFIT-decompositions, as well as our mask editing procedure (MASK\_EDIT) used in data preparation. In the second paper we will present a detailed analysis of the bulge, disk, and bar parameter derived from multi-component decompositions.
We present a general relativistic (GR) model of jet variability in active galactic nuclei due to orbiting blobs in helical motion along a funnel or cone shaped magnetic surface anchored to the accretion disk near the black hole. Considering a radiation pressure driven flow in the inner region, we find that it stabilizes the flow yielding Lorentz factors ranging between $1.1 - 7$ at small radii for reasonable initial conditions. Assuming these as inputs, simulated light curves (LCs) for the funnel model include Doppler and gravitational shifts, aberration, light bending and time delay. These LCs are studied for quasi-periodic oscillations (QPOs) and the power spectral density (PSD) shape and yield an increased amplitude ($\sim$ 12 %); a beamed portion and a systematic phase shift with respect to that from a previous special relativistic model. The results strongly justify implementing a realistic magnetic surface geometry in a GR framework to describe effects on emission from orbital features in the jet close to the horizon radius. A power law shaped PSD with a typical slope of $-2$ and QPOs with timescales in the range $(1.37 - 130.7)$ days consistent with optical variability in Blazars, emerges from the simulations for black hole masses $M_{\bullet} = (0.5 - 5) \times 10^8 M_{\odot}$ and initial Lorentz factors $\gamma_{jet,i} = 2 - 10$. The models presented here can be applied to explain radio, optical and X-ray variability from a range of jetted sources including AGN, X-ray binaries and neutron stars.
CK Vulpeculae was observed in outburst in 1670-16721, but no counterpart was seen until 1982, when a bipolar nebula was found at its location. Historically, CK Vul has been considered to be a nova (Nova Vul 1670), but a similarity to 'red transients', which are more luminous than classical nova and thought to be the result of stellar collisions, has re-opened the question of CK Vul's status. Red transients cool to resemble late M-type stars, surrounded by circumstellar material rich in molecules and dust. No stellar source has been seen in CK Vul, though a radio continuum source was identified at the expansion centre of the nebula. Here we report CK Vul is surrounded by chemically rich molecular gas with peculiar isotopic ratios, as well as dust. The chemical composition cannot be reconciled with a nova or indeed any other known explosion. In addition, the mass of the surrounding gas is too high for a nova, though the conversion from observations of CO to a total mass is uncertain. We conclude that CK Vul is best explained as the remnant of a merger of two stars.
Although the inflationary paradigm is the most widely accepted explanation for the current cosmological observations, it does not necessarily correspond to what actually happened in the early stages of our Universe. To decide on this issue, two paths can be followed: first, all the possible predictions it makes must be derived thoroughly and compared with available data, and second, all imaginable alternatives must be ruled out. Leaving the first task to all other contributors of this volume, we concentrate here on the second option, focusing on the bouncing alternatives and their consequences.
The formation of massive stars play a crucial role in galaxies from numerous points of view. The protostar generates a strong ultraviolet radiation field that ionizes its surroundings, and it drives powerful shock waves in the neighbouring medium in the form of jets and bipolar outflows, whose structure can be partially organized by local, strong magnetic field. Such an ejection activity locally modifies the interstellar chemistry, contributing to the cycle of matter. It also significantly participates to the energetic balance of galaxies. In the latter stages of massive star formation, the protostar is surrounded by an ultra-compact HII region, and irradiates its bipolar outflows, where an intrinsically strong magnetic field structure is associated to the generally high densities. In the HII region, or in the bipolar outflows, the question of in situ cosmic rays acceleration can then be raised by the simultaneous presence of strong magnetic fields, significant ionization of the matter, and mechanical energy available in large quantities. In this contribution, we will only summarise the results of our study of potential in situ cosmic rays acceleration the irradiated shocks in the W28 A2 massive star forming region, based on observations from the APEX, IRAM 30m, and Herschel telescopes, and presented in Gusdorf et al., recently submitted to A&A.
Supernova remnants (SNRs) are considered as being the sources of galactic cosmic rays. In order to understand the origin, acceleration, and composition of these cosmic rays, detailed knowledge of the physical conditions in the local interstellar medium is needed. The shock interaction of SNRs with molecular clouds that gives rise to strong molecular emission in the far-IR and sub-mm wavelength regimes can be used as a highly valuable tracer of these conditions. The application of MHD shock models in the interpretation of the resulting line emission can yield information on the energetic and chemical impact of supernova remnants. We have mapped two regions in the supernova remnant W44 with the APEX telescope in ${}^{12}$CO (3-2), (4-3), (6-5), (7-6) and ${}^{13}$CO (3-2). The extraction of integrated intensities on five different positions, corresponding to local maxima of CO emission, allows to compare these intensities to the outputs of a grid of models, which combine an MHD shock code with a radiative transfer module based on the 'large velocity gradient' approximation. We find that the observed CO line emission is compatible with non-stationary shocks and a pre-shock density of $10^4$ cm${}^{-3}$. Our models furthermore allow to constrain shock ages, velocities, the pre-shock magnetic field strength components perpendicular to the line-of-sight, and the full ladder of CO transitions. Finally, our analysis can be used to estimate the contribution of such SNRs to, e.g. the galactic energy balance and the momentum-injection into the surrounding interstellar medium.
Few high-redshift, radio-loud quasars are known to date. The extremely luminous, radio-bright quasar, SDSS J013127.34-032100.1 was recently discovered at a redshift of $z=5.18$. We observed the source with high resolution very long baseline interferometry (VLBI) at 1.7 GHz with the European VLBI Network (EVN) and found a single compact radio component. We estimated a lower limit to the brightness temperature of the detected radio component, T_B~10^{11} K. Additionaly, when compared to archival radio data, the source showed significant flux density variation. These two findings are indicative of the blazar nature of the source.
The chemical evolution in high-mass star-forming regions is still poorly constrained. Studying the evolution of deuterated molecules allows to differentiate between subsequent stages of high-mass star formation regions due to the strong temperature dependence of deuterium isotopic fractionation. We observed a sample of 59 sources including 19 infrared dark clouds, 20 high-mass protostellar objects, 11 hot molecular cores and 9 ultra-compact HII regions in the (3-2) transitions of the four deuterated molecules, DCN, DNC, DCO+ and N2D+ as well as their non-deuterated counterpart. The overall detection fraction of DCN, DNC and DCO+ is high and exceeds 50% for most of the stages. N2D+ was only detected in a few infrared dark clouds and high-mass protostellar objects. It can be related to problems in the bandpass at the frequency of the transition and to low abundances in the more evolved, warmer stages. We find median D/H ratios of ~0.02 for DCN, ~0.005 for DNC, ~0.0025 for DCO+ and ~0.02 for N2D+. While the D/H ratios of DNC, DCO+ and N2D+ decrease with time, DCN/HCN peaks at the hot molecular core stage. We only found weak correlations of the D/H ratios for N2D+ with the luminosity of the central source and the FWHM of the line, and no correlation with the H2 column density. In combination with a previously observed set of 14 other molecules (Paper I) we fitted the calculated column densities with an elaborate 1D physico-chemical model with time-dependent D-chemistry including ortho- and para-H2 states. Good overall fits to the observed data have been obtained the model. It is one of the first times that observations and modeling have been combined to derive chemically based best-fit models for the evolution of high-mass star formation including deuteration.
We present the first Swift Ultraviolet/Optical Telescope Serendipitous Source Catalogue (UVOTSSC). The catalogue was compiled from 23,059 Swift datasets taken within the first five years of observations with the Swift UVOT. A purpose-built processing pipeline, based around the standard Swift processing tools, was employed. The catalogue contains positions, photometry in three UV and three optical bands, morphological information and data quality flags. In total, the catalogue contains 6,200,016 unique sources of which more than 2 million have multiple observations in the catalogue.
Distinguishing $\gamma$-ray and hadron initiated extensive air showers (EAS) based on lateral distribution of electrons has been studied by detailed Monte Carlo (MC) simulations. The possibility of using the local age parameter (LAP) of EAS for the gamma-hadron separation has been explored. It is found that separating $\gamma$-ray and hadron induced EAS on the basis of LAP can be useful for surface detector experiments those have no reliable muon measurement facilities.
Despite the fact that the Ly$\alpha$ line is an important diagnostic for star formation at high redshift, the interpretation of its flux and line profile is difficult due to its resonance nature. Trends between the escape of Ly$\alpha$ photons and dust and ISM properties have been found, but detailed comparisons between Ly$\alpha$ emission and the properties of the gas in local high redshift analogues are vital to understand the relation between Ly$\alpha$ emission and galaxy properties. For the first time we can directly infer the properties of the ionized gas at the same location and comparable spatial scales of the extended Ly$\alpha$ halo around ESO 338-IG04. VLT/MUSE integral field spectra have been obtained. We use ionization parameter mapping of the [SII/OIII] line ratio and the kinematics of H$\alpha$ to study the ionization state and kinematics of the interstellar medium of ESO 338. The velocity map reveals two outflows. The entire central area of the galaxy is highly ionized by photons leaking from the HII regions around the youngest star clusters. Three highly ionized cones have been identified, of which one is associated with an outflow detected in the H$\alpha$. We propose a scenario where the outflows are created by mechanical feedback of the older clusters, while the highly ionized gas is caused by the hard ionizing photons emitted by the youngest clusters. A comparison with the Ly$\alpha$ map shows that the (approximately bipolar) asymmetries observed in the Ly$\alpha$ emission are consistent with the base of the outflows detected in H$\alpha$. No clear correlation with the ionization cones is found. The mechanical and ionization feedback of star clusters significantly changes the state of the ISM by creating ionized cones and outflows. The comparison with Ly$\alpha$ suggests that especially the outflows could facilitate the escape of Ly$\alpha$ photons.
Aims: To predict the cruise radiation environment related to future human missions to Mars, the correlation between solar modulation potential and the dose rate measured by the Radiation Assessment Detector (RAD) has been analyzed and empirical models have been employed to quantify this correlation. Methods: The instrument RAD, onboard Mars Science Laboratory's (MSL) rover Curiosity, measures a broad spectrum of energetic particles along with the radiation dose rate during the 253-day cruise phase as well as on the surface of Mars. With these first ever measurements inside a spacecraft from Earth to Mars, RAD observed the impulsive enhancement of dose rate during solar particle events as well as a gradual evolution of the galactic cosmic ray (GCR) induced radiation dose rate due to the modulation of the primary GCR flux by the solar magnetic field, which correlates with long-term solar activities and heliospheric rotation. Results: We analyzed the dependence of the dose rate measured by RAD on solar modulation potentials and estimated the dose rate and dose equivalent under different solar modulation conditions. These estimations help us to have approximate predictions of the cruise radiation environment, such as the accumulated dose equivalent associated with future human missions to Mars. Conclusions: The predicted dose equivalent rate during solar maximum conditions could be as low as one-fourth of the current RAD cruise measurement. However, future measurements during solar maximum and minimum periods are essential to validate our estimations.
We explore the properties of high-redshift Lyman-alpha emitters (LAE), and their link with the Lyman-Break galaxy population (LBG), using a semi-analytic model of galaxy formation that takes into account resonant scattering of Lya photons in gas outflows. We can reasonably reproduce the abundances of LAEs and LBGs from redshift 3 to 7, as well as most UV LFs of LAEs. The stronger dust attenuation for (resonant) Lya photons compared to UV continuum photons in bright LBGs provides a natural interpretation to the increase of the LAE fraction in LBG samples, X_LAE, towards fainter magnitudes. The redshift evolution of X_LAE seems however very sensitive to UV magnitudes limits and EW cuts. In spite of the apparent good match between the statistical properties predicted by the model and the observations, we find that the tail of the Lya equivalent width distribution (EW > 100 A) cannot be explained by our model, and we need to invoke additional mechanisms. We find that LAEs and LBGs span a very similar dynamical range, but bright LAEs are about 4 times rarer than LBGs in massive halos. Moreover, massive halos mainly contain weak LAEs in our model, which might introduce a bias towards low-mass halos in surveys which select sources with high EW cuts. Overall, our results are consistent with the idea that LAEs and LBGs make a very similar galaxy population. Their apparent differences seem mainly due to EW selections, UV detection limits, and a decreasing Lya-to-UV escape fraction ratio in high SFR galaxies.
We present results based on YJKs photometry of star clusters located in the outermost, eastern region of the Small Magellanic Cloud (SMC). We analysed a total of 51 catalogued clusters whose colour--magnitude diagrams (CMDs), having been cleaned from field-star contamination, were used to assess the clusters' reality and estimate ages of the genuine systems. Based on CMD analysis, 15 catalogued clusters were found to be possible non-genuine aggregates. We investigated the properties of 80% of the catalogued clusters in this part of the SMC by enlarging our sample with previously obtained cluster ages, adopting a homogeneous scale for all. Their spatial distribution suggests that the oldest clusters, log(t yr-1) >= 9.6, are in general located at greater distances to the galaxy's centre than their younger counterparts -- 9.0 <= log(t yr-1) <= 9.4 -- while two excesses of clusters are seen at log(t yr-1) ~ 9.2 and log(t yr-1) ~ 9.7. We found a trail of younger clusters which follow the Wing/Bridge components. This long spatial sequence does not only harbour very young clusters, log(t yr-1) ~ 7.3, but it also hosts some of intermediate ages, log(t yr-1) ~ 9.1. The derived cluster and field-star formation frequencies as a function of age are different. The most surprising feature is an observed excess of clusters with ages of log(t yr-1) < 9.0, which could have been induced by interactions with the LMC.
The GAMMA-400 gamma-ray telescope is intended to measure the fluxes of gamma rays and cosmic-ray electrons and positrons in the energy range from 100 MeV to several TeV. Such measurements concern with the following scientific goals: search for signatures of dark matter, investigation of gamma-ray point and extended sources, studies of the energy spectra of Galactic and extragalactic diffuse emission, studies of gamma-ray bursts and gamma-ray emission from the active Sun, as well as high-precision measurements of spectra of high-energy electrons and positrons, protons, and nuclei up to the knee. The main components of cosmic rays are protons and helium nuclei, whereas the part of lepton component in the total flux is ~10E-3 for high energies. In present paper, the capability of the GAMMA-400 gamma-ray telescope to distinguish electrons and positrons from protons in cosmic rays is investigated. The individual contribution to the proton rejection is studied for each detector system of the GAMMA-400 gamma-ray telescope. Using combined information from all detector systems allow us to provide the proton rejection from electrons with a factor of ~4x10E5 for vertical incident particles and ~3x10E5 for particles with initial inclination of 30 degrees. The calculations were performed for the electron energy range from 50 GeV to 1 TeV.
We present new magnetic field measurements of the K2 main-sequence star $\epsilon$ Eridani based on principal components analysis (PCA) line-profile reconstructions. The aim of this paper is to quantify the surface-averaged magnetic field and search for possible variations. A total of 338 optical echelle spectra from our robotic telescope facility STELLA with a spectral resolution of 55,000 were available for analysis. This time-series was used to search for the small line-profile variations due to a surface magnetic field with the help of a PCA. Evidence for a spatial and temporal inhomogeneous magnetic field distribution is presented. The mean, surface averaged, magnetic field strength was found to be $\langle B \rangle = 186 \pm 47\,\mathrm{G}$ in good agreement with previous Zeeman-broadening measurements. Clear short-term variations of the surface averaged magnetic field of up to few tens Gauss were detected together with evidence for a three-year cycle in the surface-averaged magnetic field of $\epsilon$ Eridani.
Aims: We develop a method for estimating the properties of stellar winds for
low-mass main-sequence stars between masses of 0.4 and 1.1 solar masses at a
range of distances from the star.
Methods: We use 1D thermal pressure driven hydrodynamic wind models run using
the Versatile Advection Code. Using in situ measurements of the solar wind, we
produce models for the slow and fast components of the solar wind. We consider
two radically different methods for scaling the base temperature of the wind to
other stars: in Model A, we assume that wind temperatures are fundamentally
linked to coronal temperatures, and in Model B, we assume that the sound speed
at the base of the wind is a fixed fraction of the escape velocity. In Paper II
of this series, we use observationally constrained rotational evolution models
to derive wind mass loss rates.
Results: Our model for the solar wind provides an excellent description of
the real solar wind far from the solar surface, but is unrealistic within the
solar corona. We run a grid of 1200 wind models to derive relations for the
wind properties as a function of stellar mass, radius, and wind temperature.
Using these results, we explore how wind properties depend on stellar mass and
rotation.
Conclusions: Based on our two assumptions about the scaling of the wind
temperature, we argue that there is still significant uncertainty in how these
properties should be determined. Resolution of this uncertainty will probably
require both the application of solar wind physics to other stars and detailed
observational constraints on the properties of stellar winds. In the final
section of this paper, we give step by step instructions for how to apply our
results to calculate the stellar wind conditions far from the stellar surface.
Theories which modify general relativity to explain the accelerated expansion of the Universe often use screening mechanisms to satisfy constraints on Solar System scales. We investigate the effects of the cosmic web and the local environmental density of dark matter halos on the screening properties of the Vainshtein and chameleon screening mechanisms. We compare the cosmic web morphology of dark matter particles, mass functions of dark matter halos, mass and radial dependence of screening, velocity dispersions and peculiar velocities, and environmental dependence of screening mechanisms in $f(R)$ and nDGP models. Using the ORIGAMI cosmic web identification routine we find that the Vainshtein mechanism depends on the cosmic web morphology of dark matter particles, since these are defined according to the dimensionality of their collapse, while the chameleon mechanism shows no morphology dependence. The chameleon screening of halos and their velocity dispersions depend on halo mass, and small halos and subhalos can be environmentally screened in the chameleon mechanism. On the other hand, the screening of halos in the Vainshtein mechanism does not depend on mass nor environment, and their velocity dispersions are suppressed. The peculiar velocities of halos in the Vainshtein mechanism are enhanced because screened objects can still feel the fifth force generated by external fields, while peculiar velocities of chameleon halos are suppressed when the halo centers are screened.
I have created, for public use, a Python code allowing the simulation of
light curves with any given power spectral density and any probability density
function, following the algorithm described in Emmanoulopoulos et al. 2013. The
simulated products have exactly the same variability and statistical properties
as the observed light curves. The code and its documentation are available at:
https://github.com/samconnolly/DELightcurveSimulation
Note that a Mathematica code of the algorithm is given in Emmanoulopoulos et
al. [arXiv:1305.0304]
We provide a status report on the determination of stellar ages from asteroseismology for stars of various masses and evolutionary stages. The ability to deduce the ages of stars with a relative precision of typically 10 to 20% is a unique opportunity for stellar evolution and also of great value for both galactic and exoplanet studies. Further, a major uncalibrated ingredient that makes stellar evolution models uncertain, is the stellar interior rotation frequency $\Omega(r)$ and its evolution during stellar life. We summarize the recent achievements in the derivation of $\Omega(r)$ for different types stars, offering stringent observational constraints on theoretical models. Core-to-envelope rotation rates during the red giant stage are far lower than theoretical predictions, pointing towards the need to include new physical ingredients that allow strong and efficient coupling between the core and the envelope in the models of low-mass stars in the evolutionary phase prior to the core helium burning. Stars are subject to efficient mixing phenomena, even at low rotation rates. Young massive stars with seismically determined interior rotation frequency reveal low core-to-envelope rotation values.
The results of Mopra molecular spectral line observations towards the supernova remnant HESSJ1731-347 (G353.6-0.7) and the unidentified gamma-ray source HESSJ1729-345 are presented. Dense molecular gas in three different velocity-bands (corresponding to three Galactic arms) are investigated using the CS(1-0) line. The CS-traced component provides information about the dense target material in a hadronic scenario for gamma-ray production (cosmic rays interacting with gas) and an understanding of the dynamics. Furthermore, the effects of cosmic ray diffusion into dense gas may alter the gamma-ray spectrum to cause a flattening of spectra towards such regions. Dense molecular gas mass at a level of ~10^5 Mo was revealed in this survey, with mass of the order of ~10^3 Mo towards HESSJ1729-345 in each coincident Galactic arm, but no significant detection of dense molecular gas towards HESSJ1731-347 at the currently-preferred distance of ~5.2-6.2 kpc was discovered.
Recent extensive observations of Type Ia Supernovae (SNe Ia) have revealed the existence of a diversity of SNe Ia, including SN 2002cx-like objects (also called SN Iax). We introduce two possible channels in the single degenerate scenario: 1) double detonations in sub-Chandrasekhar (Ch) mass CO white dwarfs (WDs), where a thin He envelope is developed with relatively low accretion rates after He novae even at low metallicities, and 2) carbon deflagrations in Ch-mass possibly hybrid C+O+Ne WDs, where WD winds occur at [Fe/H] ~ -2.5 at high accretion rates. These subclasses of SNe Ia are rarer than `normal' SNe Ia and do not affect the chemical evolution in the solar neighborhood, but can be very important in metal-poor systems with stochastic star formation. In dwarf spheroidal galaxies in the Local Group, the decrease of [\alpha/Fe] ratios at [Fe/H] ~ -2 to -1.5 can be produced depending on the star formation history. SNe Iax give high [Mn/Fe], while sub-Ch-mass SNe Ia give low [Mn/Fe], and thus a model including a mix of the two is favoured by the available observations.
The proximity profile in the spectra of z~3 quasars, where fluxes extend
blueward of the He II Lya wavelength 304 (1+z) A, is one of the most important
spectral features in the study of the intergalactic medium. Based on the HST
spectra of 24 He II quasars, we find that the majority of them display a
proximity profile, corresponding to an ionization radius as large as 20 Mpc in
the source's rest frame. In comparison with those in the H i spectra of the
quasars at z~6, the He II proximity effect is more prominent and is observed
over a considerably longer period of reionization. The He II proximity zone
sizes decrease at higher redshifts, particularly at z > 3.3. This trend is
similar to that for H I, signaling an onset of He II reionization at z~4.
For quasar SDSS1253+6817 (z=3.48), the He II absorption trough displays a
gradual decline and serves a good case for modeling the He II reionization. To
model such a broad profile requires a quasar radiation field whose distribution
between 4 and 1 Rydberg is considerably harder than normally assumed. The UV
continuum of this quasar is indeed exceptionally steep, and the He II
ionization level in the quasar vicinity is higher than the average level in the
intergalactic medium. These results are evidence that a very hard EUV continuum
from this quasar produces a large ionized zone around it.
Distinct exceptions are the two brightest He II quasars at z~2.8, for which
no significant proximity profile is present, possibly implying that they are
young.
We discuss how recent advances in observations, theory and numerical simulations have allowed the stellar community to progress in its understanding of stellar convection, rotation and magnetism and to assess the degree to which the Sun and other stars share similar dynamical properties. Ensemble asteroseismology has become a reality with the advent of large time domain studies, especially from space missions. This new capability has provided improved constraints on stellar rotation and activity, over and above that obtained via traditional techniques such as spectropolarimetry or CaII H&K observations. New data and surveys covering large mass and age ranges have provided a wide parameter space to confront theories of stellar magnetism. These new empirical databases are complemented by theoretical advances and improved multi-D simulations of stellar dynamos. We trace these pathways through which a lucid and more detailed picture of magnetohydrodynamics of solar-like stars is beginning to emerge and discuss future prospects.
Well-calibrated spectropolarimetry studies at resolutions of $R>$10,000 with signal-to-noise ratios (SNRs) better than 0.01\% across individual line profiles, are becoming common with larger aperture telescopes. Spectropolarimetric studies require high SNR observations and are often limited by instrument systematic errors. As an example, fiber-fed spectropolarimeters combined with advanced line-combination algorithms can reach statistical error limits of 0.001\% in measurements of spectral line profiles referenced to the continuum. Calibration of such observations is often required both for cross-talk and for continuum polarization. This is not straightforward since telescope cross-talk errors are rarely less than $\sim$1\%. In solar instruments like the Daniel K. Inouye Solar Telescope (DKIST), much more stringent calibration is required and the telescope optical design contains substantial intrinsic polarization artifacts. This paper describes some generally useful techniques we have applied to the HiVIS spectropolarimeter at the 3.7m AEOS telescope on Haleakala. HiVIS now yields accurate polarized spectral line profiles that are shot-noise limited to 0.01\% SNR levels at our full spectral resolution of 10,000 at spectral sampling of $\sim$100,000. We show line profiles with absolute spectropolarimetric calibration for cross-talk and continuum polarization in a system with polarization cross-talk levels of essentially 100\%. In these data the continuum polarization can be recovered to one percent accuracy because of synchronized charge-shuffling model now working with our CCD detector. These techniques can be applied to other spectropolarimeters on other telescopes for both night and day-time applications such as DKIST, TMT and ELT which have folded non-axially symmetric foci.
We explore the kinematics (both the radial velocity and the proper motion) of the vertical X-shaped feature in the Milky Way with an N-body bar/bulge model. From the solar perspective, the distance distribution of particles is double-peaked in fields passing through the X-shape. The separation and amplitude ratio between the two peaks qualitatively match the observed trends towards the Galactic bulge. We confirm clear signatures of cylindrical rotation in the pattern of mean radial velocity across the bar/bulge region. We also find possible imprints of coherent orbital motion inside the bar structure in the radial velocity distribution along l=0 degree, where the near and far sides of the bar/bulge show excesses of approaching and receding particles. The coherent orbital motion is also reflected in the slight displacement of the zero-velocity-line in the mean radial velocity, and the displacement of the maximum/minimum in the mean longitudinal proper motion across the bulge region. We find some degree of anisotropy in the stellar velocity within the X-shape, but the underlying orbital family of the X-shape cannot be clearly distinguished. Two potential applications of the X-shape in previous literature are tested, i.e., bulge rotation and Galactic center measurements. We find that the proper motion difference between the two sides of the X-shape can be used to estimate the mean azimuthal streaming motion of the bulge, but not the pattern speed of the bar. We also demonstrate that the Galactic center can be located with the X-shape, but the accuracy depends on the fitting scheme, the number of fields, and their latitudinal coverage.
This work studies jet-like electromagnetic configurations surrounding a slowly-spinning black-hole immersed in a uniformly magnetized force-free plasma. In the first part of this work, we present a family of stationary solutions that are jet-capable. While these solutions all satisfy the force-free equations and the appropriate boundary conditions, our numerical experiments show a unique relaxed state starting from different initial data, and so one member of the family is likely preferred over the others. In the second part of this work, we analyze the perturbations of this family of jet-like solutions, and show that the perturbative modes exhibit a similar split into the trapped and traveling categories previously found for perturbed Blandford-Znajek solutions. In the eikonal limit, the trapped modes can be identified with the fast magnetosonic waves in the force-free plasma and the traveling waves are essentially the Alfven waves. Moreover, within the scope of our analysis, we have not seen signs of unstable modes at the light-crossing timescale of the system, within which the numerical relaxation process occurs. This observation disfavors mode instability as the selection mechanism for picking out a preferred solution. Consequently, our analytical study is unable to definitively select a particular solution out of the family to serve as the aforementioned preferred final state. This remains an interesting open problem.
We analyze the parametric space of the constrained minimal supersymmetric standard model with mu>0 supplemented by a generalized asymptotic Yukawa coupling quasi-unification condition which yields acceptable masses for the fermions of the third family. We impose constraints from the cold dark matter abundance in the universe and its direct detection experiments, the B-physics, as well as the masses of the sparticles and the lightest neutral CP-even Higgs boson. Fixing the mass of the latter to its central value from the LHC and taking 40<=tanbeta<=50, we find a relatively wide allowed parameter space with -11<=A_0/M_{1/2}<=15 and mass of the lightest sparticle in the range (0.09-1.1) TeV. This sparticle is possibly detectable by the present cold dark matter direct search experiments. The required fine-tuning for the electroweak symmetry breaking is much milder than the one needed in the neutralino-stau coannihilation region of the same model.
In this present work, we try to build up a cosmological model using a non-canonical scalar field within the framework of a spatially flat FRW space-time. In this context, we have considered four different parametrizations of the equation of state parameter of the non- canonical scalar field. Under this scenario, an analytical solution for the various cosmological parameters have been found out. It has been found that the deceleration parameter shows a smooth transition from a positive value to some negative value which indicates that the universe was undergoing an early deceleration followed by late time acceleration which is essential for the structure formation of the universe. With these four parametrizations, the future evolution of the models are also discussed. We have also shown that the two models mimic as the concordance $\Lambda$CDM in the near future, whereas the other two models diverge due to the future singularity. Finally, we have studied these theoretical models with the Union2.1 SN Ia dataset.
In the approach of the effective field theory of modified gravity, we derive the equations of motion for linear perturbations in the presence of a barotropic perfect fluid on the flat isotropic cosmological background. In a simple version of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, which is the minimum extension of Horndeski theories, we show that a slight deviation of the tensor propagation speed squared $c_{\rm t}^2$ from 1 generally leads to the large modification to the propagation speed squared $c_{\rm s}^2$ of a scalar degree of freedom $\phi$. This problem persists whenever the kinetic energy $\rho_X$ of the field $\phi$ is much smaller than the background energy density $\rho_m$, which is the case for most of dark energy models in the asymptotic past. Since the scaling solution characterized by the constant ratio $\rho_X/\rho_m$ is one way out for avoiding such a problem, we study the evolution of perturbations for a scaling dark energy model in the framework of GLPV theories in the Jordan frame. Provided the oscillating mode of scalar perturbations is fine-tuned so that it is initially suppressed, the anisotropic parameter $\eta=-\Phi/\Psi$ between the two gravitational potentials $\Psi$ and $\Phi$ significantly deviates from 1 for $c_{\rm t}^2$ away from 1. For other general initial conditions, the deviation of $c_{\rm t}^2$ from 1 gives rise to the large oscillation of $\Psi$ with the frequency related to $c_{\rm s}^2$. In both cases, the model can leave distinct imprints for the observations of CMB and weak lensing.
We report on numerical calculations of orbital and epicyclic frequencies in nearly circular orbits around rotating neutron stars and strange quark stars. The FPS equation of state was used to describe the neutron star structure while the MIT bag model was used to model the equation of state of strange quark stars. The uniformly rotating stellar configurations were computed in full general relativity. We find that the vertical epicyclic frequency is very sensitive to the oblateness of the rotating star. For models of rotating neutron stars of moderate mass, as well as for strange quark star models, the sense of the nodal precession of test particle orbits close to the star changes at a certain stellar rotation rate. These findings may have implications for models of kHz QPOs.
We investigate the evolution of the baryon asymmetry of the Universe (BAU) in its symmetric phase before the electroweak phase transition (EWPT) induced by leptogenesis in the hypermagnetic field of an arbitrary structure and with a maximum hypermagnetic helicity density. The novelty of this work is that the BAU has been calculated for a continuous hypermagnetic helicity spectrum. The observed BAU $B_{obs} = 10^{-10}$ that can be in large-scale hypermagnetic fields satisfying the wave number inequality $k \leq k_{max}$ grows with increasing $k_{max}$. We will also show that the initial right-handed electron asymmetry $\xi_{eR}(\eta_0)$ used in our leptogenesis model as a free parameter cannot take too large values, $\xi_{eR}(\eta_0) = 10^{-4}$, because this leads to a negative BAU by the EWPT time. In contrast, a sufficiently small initial right-handed electron asymmetry, $\xi_{eR}(\eta_0)$, provides its further growth and the corresponding BAU growth from zero to some positive value,including the observed $B_{obs} = 10^{-10}$.
We derive the primordial power spectra and spectral indexes of the density fluctuations and gravitational waves in the framework of loop quantum cosmology (LQC) with holonomy and inverse-volume corrections, by using the uniform asymptotic approximation method to its third-order, at which the upper error bounds are $\lesssim 0.15\%$, accurate enough for the current and forthcoming cosmological observations. Then, using the Planck, BAO and SN data we obtain new constraints on quantum gravitational effects from LQC corrections, and find that such effects could be well within the detection of the current and forthcoming experiments.
We search for self tuning solutions to the Einstein-scalar field equations for the simplest class of `Fab-Four' models with constant potentials. We first review the conditions under which self tuning occurs in a cosmological spacetime, and by introducing a small modification to the original theory - introducing the second and third Galileon terms - show how one can obtain de Sitter states where the expansion rate is independent of the vacuum energy. We then consider whether the same self tuning mechanism can persist in a spherically symmetric inhomogeneous spacetime. We show that there are no asymptotically flat solutions to the field equations in which the vacuum energy is screened, other than the trivial one (Minkowski space). We then consider the possibility of constructing Schwarzschild de Sitter spacetimes for the modified Fab Four plus Galileon theory. We argue that the only model that can successfully screen the vacuum energy in both an FLRW and Schwarzschild de Sitter spacetime is one containing `John' $\sim G^{\mu}{}_{\nu} \partial_{\mu}\phi\partial^{\nu}\phi$ and a canonical kinetic term $\sim \partial_{\alpha}\phi \partial^{\alpha}\phi$. This behaviour was first observed in (Babichev&Charmousis,2013). The screening mechanism, which requires redundancy of the scalar field equation in the `vacuum', fails for the `Paul' term in an inhomogeneous spacetime.
The cosmological models called $\alpha$-attractors provide an excellent fit to the latest observational data. Their predictions $n_{s} = 1-2/N$ and $r = 12\alpha/N^{2}$ are very robust with respect to the modifications of the inflaton potential. An intriguing interpretation of $\alpha$-attractors is based on a geometric moduli space with a boundary: a Poincare disk model of a hyperbolic geometry with the radius $\sqrt{3\alpha}$, beautifully represented by the Escher's picture Circle Limit IV. In such models, the amplitude of the gravitational waves is proportional to the square of the radius of the Poincare disk.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
Whether a short gamma-ray burst (GRB) is caused by a black hole (BH) or a neutron star (NS) after the merger of a NS binary is a crucial problem. We propose a BH model that explains short GRBs with long-lasting activities such as extended emission and plateau emission up to $\sim10000$ s. To extract the BH rotational energy, the topological evolution of the magnetic field should accompany the mass ejection, mass fallback, and magnetic field reconnection. The observations suggest the magnetic field decay from $\sim10^{14}$ G to $\sim10^{13} - 10^{11}$ G at the BH, bounded below by the pre-merger strength and kept constant while the luminosity is constant, and the fallback mass of $\sim10^{-4} - 10^{-2} M_{\odot}$, comparable to the ejecta mass implied by the macronova in GRB 130603B. The BH model has implications for gravitational waves and the equation of state of NS matter.
We investigate the variation of the ratio of the equivalent widths of the FeII$\lambda$2600 line to the MgII$\lambda\lambda$2796,2803 doublet as a function of redshift in a large sample of absorption lines drawn from the JHU-SDSS Absorption Line Catalog. We find that despite large scatter, the observed ratio shows a trend where the equivalent width ratio $\mathcal{R}\equiv W_{\rm FeII}/W_{\rm MgII}$ decreases monotonically with increasing redshift $z$ over the range $0.55 \le z \le 1.90$. Selecting the subset of absorbers where the signal-to-noise ratio of the MgII equivalent width $W_{\rm MgII}\ge$3 and modeling the equivalent width ratio distribution as a gaussian, we find that the mean of the gaussian distribution varies as $\mathcal{R}\propto (-0.045\pm0.005)z$. We discuss various possible reasons for the trend. A monotonic trend in the Fe/Mg abundance ratio is predicted by a simple model where the abundances of Mg and Fe in the absorbing clouds are assumed to be the result of supernova ejecta and where the cosmic evolution in the SNIa and core-collapse supernova rates is related to the cosmic star-formation rate. If the trend in $\mathcal{R}$ reflects the evolution in the abundances, then it is consistent with the predictions of the simple model.
Collisions between giant molecular clouds are a potential mechanism for triggering the formation of massive stars, or even super star clusters. The trouble is identifying this process observationally and distinguishing it from other mechanisms. We produce synthetic position-velocity diagrams from models of: cloud-cloud collisions, non-interacting clouds along the line of sight, clouds with internal radiative feedback and a more complex cloud evolving in a galactic disc, to try and identify unique signatures of collision. We find that a broad bridge feature connecting two intensity peaks, spatially correlated but separated in velocity, is a signature of a high velocity cloud-cloud collision. We show that the broad bridge feature is resilient to the effects of radiative feedback, at least to around 2.5Myr after the formation of the first massive (ionising) star. However for a head on 10km/s collision we find that this will only be observable from 20-30 per cent of viewing angles. Such broad-bridge features have been identified towards M20, a very young region of massive star formation that was concluded to be a site of cloud-cloud collision by Torii et al (2011), and also towards star formation in the outer Milky Way by Izumi et al (2014).
Collider, direct and indirect searches for dark matter have typically little or no sensitivity to weakly interacting massive particles (WIMPs) with masses above a few TeV. This rather unexplored regime can however be probed through the search for distinctive gamma-ray spectral features produced by the annihilation of WIMPs at very high energies. Here we present a dedicated search for gamma-ray boxes -- sharp spectral features that cannot be mimicked by astrophysical sources -- with the upcoming Cherenkov Telescope Array (CTA). Using realistic projections for the instrument performance and detailed background modelling, a profile likelihood analysis is implemented to derive the expected upper limits and sensitivity reach after 100 h of observations towards a $2^\circ\times2^\circ$ region around the Galactic centre. Our results show that CTA will be able to probe gamma-ray boxes down to annihilation cross sections of $10^{-27}-10^{-26}\,\text{cm}^3\text{/s}$ up to tens of TeV. We also identify a number of concrete particle physics models providing thermal dark matter candidates that can be used as target benchmarks in future search campaigns. This constitutes a golden opportunity for CTA to either discover or rule out multi-TeV thermal dark matter in a corner of parameter space where all other experimental efforts are basically insensitive.
In the Local Group, nearly all of the dwarf galaxies (M_star < 10^9 M_sun) that are satellites within 300 kpc (the virial radius) of the Milky Way (MW) and Andromeda (M31) have quiescent star formation and little-to-no cold gas. This contrasts strongly with comparatively isolated dwarf galaxies, which are almost all actively star-forming and gas-rich. This near dichotomy implies a rapid transformation after falling into the halos of the MW or M31. We combine the observed quiescent fractions for satellites of the MW and M31 with the infall times of satellites from the ELVIS suite of cosmological simulations to determine the typical timescales over which environmental processes within the MW/M31 halos remove gas and quench star formation in low-mass satellite galaxies. The quenching timescales for satellites with M_star < 10^8 M_sun are short, < 2 Gyr, and quenching is more rapid at lower M_star. These quenching timescales can be 1 - 2 Gyr longer if environmental preprocessing in lower-mass groups prior to MW/M31 infall is important. We compare with quenching timescales for more massive satellites from previous works and synthesize the nature of satellite quenching across the observable range of M_star = 10^{3-11} M_sun. The environmental quenching timescale increases rapidly with satellite M_star, peaking at ~9.5 Gyr for M_star ~ 10^9 M_sun, and rapidly decreases at higher M_star to < 5 Gyr at M_star > 5x10^9 M_sun. Overall, galaxies with M_star ~ 10^9 M_sun, similar to the Magellanic Clouds, exhibit the longest quenching timescales, regardless of environmental or internal mechanisms.
To better understand the nature of the multiphase material found in outflowing galaxies, we study the evolution of cold clouds embedded in flows of hot and fast material. Using a suite of adaptive-mesh refinement simulations that include radiative cooling, we investigate both cloud mass loss and cloud acceleration under the full range of conditions observed in galaxy outflows. The simulations are designed to track the cloud center of mass, enabling us to study the cloud evolution at long disruption times. For supersonic flows, a Mach cone forms around the cloud, which damps the Kelvin-Helmholtz instability but also establishes a streamwise pressure gradient that stretches the cloud apart. If time is expressed in units of the cloud crushing time, both the cloud lifetime and the cloud acceleration rate are independent of cloud radius, and we find simple scalings for these quantities as a function of the Mach number of the external medium. A resolution study suggests that our simulations have sufficient resolution to accurately describe the evolution of cold clouds in the absence of thermal conduction and magnetic fields, physical processes whose roles will be studied in forthcoming papers.
The vast majority of dwarf satellites orbiting the Milky Way and M31 are quenched, while comparable galaxies in the field are gas-rich and star-forming. Assuming that this dichotomy is driven by environmental quenching, we use the ELVIS suite of N-body simulations to constrain the characteristic timescale upon which satellites must quench following infall into the virial volumes of their hosts. The high satellite quenched fraction observed in the Local Group demands an extremely short quenching timescale (~ 2 Gyr) for dwarf satellites in the mass range Mstar ~ 10^6-10^8 Msun. This quenching timescale is significantly shorter than that required to explain the quenched fraction of more massive satellites (~ 8 Gyr), both in the Local Group and in more massive host halos, suggesting a dramatic change in the dominant satellite quenching mechanism at Mstar < 10^8 Msun. Combining our work with the results of complementary analyses in the literature, we conclude that the suppression of star formation in massive satellites (Mstar ~ 10^8 - 10^11 Msun) is broadly consistent with being driven by starvation, such that the satellite quenching timescale corresponds to the cold gas depletion time. Below a critical stellar mass scale of ~ 10^8 Msun, however, the required quenching times are much shorter than the expected cold gas depletion times. Instead, quenching must act on a timescale comparable to the dynamical time of the host halo. We show that ram-pressure stripping can naturally explain this behavior, with the critical mass (of Mstar ~ 10^8 Msun) corresponding to halos with gravitational restoring forces that are too weak to overcome the drag force encountered when moving through an extended, hot circumgalactic medium.
We consider the formation of satellites around the Pluto-Charon binary. An early collision between the two partners likely produced the binary and a narrow ring of debris, out of which arose the moons -- Styx, Nix, Kerberos and Hydra. Yet how the satellites emerged from the compact ring is uncertain. Here we show that a particle ring can spread from physical collisions and collective gravitational scattering, similar to migration. Around a binary, these processes take place in the reference frames of `most circular' orbits, akin to circular ones in a Keplerian potential. Ring particles can damp to these orbits, avoiding destructive collisions. Damping and diffusion can also help particles survive dynamical instabilities driven by resonances with the binary. In some situations, particles get trapped near resonances that sweep outward with the tidal evolution of the Pluto-Charon binary. With simple models and numerical experiments, we show how the Pluto-Charon impact ring may have expanded into a broad disk, out of which grew the circumbinary moons. In some scenarios, the ring can spread well beyond the orbit of Hydra, the most distant moon, to form a handful of smaller satellites. If these small moons exist, New Horizons will find them.
Interactions between a supermassive black-hole binary and the surrounding accretion disc can both assist the binary inspiral and align the black-hole spins to the disc angular momentum. While binary migration is due to angular-momentum transfer within the circumbinary disc, the spin-alignment process is driven by the mass accreting onto each black hole. Mass transfer between different disc components thus couples the inspiral and the alignment process together. Mass is expected to leak through the cavity cleared by the binary, and preferentially accretes onto the lighter (secondary) black-hole which orbits closer to the disc edge. Low accretion rate onto the heavier (primary) black hole slows the alignment process down. We revisit the problem and develop a semi-analytical model to describe the coupling between gas-driven inspiral and spin alignment, finding that binaries with mass ratio q<~0.2 approach the gravitational-wave driven inspiral in differential misalignment: light secondaries prevent primaries from aligning. Binary black holes with misaligned primaries are ideal candidates for precession effects in the strong-gravity regime and may suffer from moderately large (~1500 km/s) recoil velocities.
We present an HCO$^{+}$ J=3-2 survey of Class 0+I and Flat SED young stellar objects (YSOs) found in the Gould Belt clouds by surveys with Spitzer. Our goal is to provide a uniform Stage 0+I source indicator for these embedded protostar candidates. We made single point HCO$^{+}$ J = 3-2 measurements toward the source positions at the CSO and APEX of 546 YSOs (89% of the Class 0+I + Flat SED sample). Using the criteria from van Kempen et al. (2009), we classify sources as Stage 0+I or bona fide protostars and find that 84% of detected sources meet the criteria. We recommend a timescale for the evolution of Stage 0+I (embedded protostars) of 0.54 Myr. We find significant correlations of HCO$^{+}$ integrated intensity with ${\alpha}$ and $T_{bol}$ but not with $L_{bol}$. The detection fraction increases smoothly as a function of ${\alpha}$ and $L_{bol}$, while decreasing smoothly with $T_{bol}$. Using the Stage 0+I sources tightens the relation between protostars and high extinction regions of the cloud; 89% of Stage I sources lie in regions with $A_{V}$ >8 mag. Class 0+I and Flat SED YSOs that are not detected in HCO$^{+}$ have, on average, a factor of ~2 higher $T_{bol}$ and a factor of ~5 lower $L_{bol}$ than YSOs with HCO$^{+}$ detections. We find less YSO contamination, defined as the number of undetected YSOs divided by the total number surveyed, for sources with $T_{bol}\lesssim$ 600 K and $L_{bol} \gtrsim$ 1 $L_{\odot}$. The contamination percentage is >90% at $A_{V}$< 4 mag and decreases as $A_{V}$ increases.
As long as magnetic fields remain frozen into the gas, the magnetic braking prevents the formation of protostellar discs. This condition is subordinate to the ionisation fraction characterising the inmost parts of a collapsing cloud. The ionisation level is established by the number and the energy of the cosmic rays able to reach these regions. Adopting the method developed in our previous studies, we computed how cosmic rays are attenuated as a function of column density and magnetic field strength. We applied our formalism to low- and high-mass star formation models obtained by numerical simulations of gravitational collapse that include rotation and turbulence. In general, we found that the decoupling between gas and magnetic fields, condition allowing the collapse to go ahead, occurs only when the cosmic-ray attenuation is taken into account with respect to a calculation in which the cosmic-ray ionisation rate is kept constant. We also found that the extent of the decoupling zone also depends on the dust grain size distribution and is larger if large grains (of radius about 0.1 microns) are formed by compression and coagulation during cloud collapse. The decoupling region disappears for the high-mass case due to magnetic field diffusion that is caused by turbulence and that is not included in the low-mass models. We infer that a simultaneous study of the cosmic-ray propagation during the cloud's collapse may lead to values of the gas resistivity in the innermost few hundred AU around a forming protostar that is higher than generally assumed.
We use the most extensive integral field spectroscopic map of a local galaxy, NGC 628, combined with gas and stellar mass surface density maps, to study the distribution of metals in this galaxy out to 3 effective radii ($\rm R_e$). At each galactocentric distance, we compute the metal budget and thus constrain the mass of metals lost. We find that in the disc about half of the metals have been lost throughout the lifetime of the galaxy. The fraction of metals lost is higher in the bulge ($\sim$70%) and decreases towards the outer disc ($\rm \sim 3 \ R_e$). In contrast to studies based on the gas kinematics, which are only sensitive to ongoing outflow events, our metal budget analysis enables us to infer the average outflow rate during the galaxy lifetime. By using simple physically motivated models of chemical evolution we can fit the observed metal budget at most radii with an average outflow loading factor of order unity, thus clearly demonstrating the importance of outflows in the evolution of disc galaxies. The observed gas phase metallicity is higher than expected from the metal budget and suggests late-time accretion of enriched gas, likely raining onto the disc from the metal-enriched halo.
The galaxy power spectrum encodes a wealth of information about cosmology and the matter fluctuations. Its unbiased and optimal estimation is therefore of great importance. In this paper we generalise the framework of Feldman et al. (1994) to take into account the fact that galaxies are not simply a Poisson sampling of the underlying dark matter distribution. Besides finite survey-volume effects and flux-limits, our optimal estimation scheme incorporates several of the key tenets of galaxy formation: galaxies form and reside exclusively in dark matter haloes; a given dark matter halo may host several galaxies of various luminosities; galaxies inherit part of their large-scale bias from their host halo. Under these broad assumptions, we prove that the optimal weights "do not" explicitly depend on galaxy luminosity, other than through defining the maximum survey volume and effective galaxy density at a given position. Instead, they depend on the bias associated with the host halo; the first and second factorial moments of the halo occupation distribution; a selection function, which gives the fraction of galaxies that can be observed in a halo of mass M at position r in the survey; and an effective number density of galaxies. If one wishes to reconstruct the matter power spectrum, then, provided the model is correct, this scheme provides the only unbiased estimator. The practical challenges with implementing this approach are also discussed.
We investigate the limits of applicability of the quasi-static approximation in cosmologies featuring general models of dark energy or modified gravity. We show that the quasi-static approximation always breaks down outside of the sound horizon of the dark-energy, rather than the cosmological horizon as is frequently assumed. When the sound speed of dark energy is significantly below that of light, the quasi-static limit is only valid in a limited range of observable scales and this must be taken into account when computing effects on observations in such models. In particular, in the analysis of data from today's weak-lensing and peculiar-velocity surveys, dark energy can be modelled as quasi-static only if the sound speed is larger than order 1% of that of light. In upcoming surveys, such as Euclid, it should only be used when the sound speed exceeds around 10% of the speed of light.
The high-mass star-forming site G009.62-00.20E hosts the 6.7 GHz methanol maser source with the greatest flux density in the Galaxy which has been flaring periodically over the last ten years. We performed high-resolution astrometric measurements of the CH3OH, H2O, and OH maser emission and 7 mm continuum in the region. The radio continuum emission was resolved in two sources separated by 1300 AU. The CH3OH maser cloudlets are distributed along two north-south ridges of emission to the east and west of the strongest radio continuum component. This component likely pinpoints a massive young stellar object which heats up its dusty envelope, providing a constant IR pumping for the Class II CH3OH maser transitions. We suggest that the periodic maser activity may be accounted for by an independent, pulsating, IR radiation field provided by a bloated protostar in the vicinity of the brightest masers. We also report about the discovery of an elliptical distribution of CH3OH maser emission in the region of periodic variability.
We present a morphological study of nebulae around Wolf-Rayet (WR) stars using archival narrow-band optical and Wide-field Infrared Survey Explorer (WISE) infrared images. The comparison among WISE images in different bands and optical images proves to be a very efficient procedure to identify the nebular emission from WR nebulae, and to disentangle it from that of the ISM material along the line of sight. In particular, WR nebulae are clearly detected in the WISE W4 band at 22 $\mu$m. Analysis of available mid-IR Spitzer spectra shows that the emission in this band is dominated by thermal emission from dust spatially coincident with the thin nebular shell or most likely with the leading edge of the nebula. The WR nebulae in our sample present different morphologies that we classified into well defined WR bubbles (bubble ${\cal B}$-type nebulae), clumpy and/or disrupted shells (clumpy/disrupted ${\cal C}$-type nebulae), and material mixed with the diffuse medium (mixed ${\cal M}$-type nebulae). The variety of morphologies presented by WR nebulae shows a loose correlation with the central star spectral type, implying that the nebular and stellar evolutions are not simple and may proceed according to different sequences and time-lapses. We report the discovery of an obscured shell around WR35 only detected in the infrared.
The statistics of extrasolar planetary systems indicate that the default mode of planet formation generates planets with orbital periods shorter than 100 days, and masses substantially exceeding that of the Earth. When viewed in this context, the Solar System is unusual. Here, we present simulations which show that a popular formation scenario for Jupiter and Saturn, in which Jupiter migrates inward from a > 5 AU to a ~ 1.5 AU before reversing direction, can explain the low overall mass of the Solar System's terrestrial planets, as well as the absence of planets with a < 0.4 AU. Jupiter's inward migration entrained s ~ 10-100 km planetesimals into low-order mean-motion resonances, shepherding and exciting their orbits. The resulting collisional cascade generated a planetesimal disk that, evolving under gas drag, would have driven any pre-existing short-period planets into the Sun. In this scenario, the Solar System's terrestrial planets formed from gas-starved mass-depleted debris that remained after the primary period of dynamical evolution.
We explore the linear stability of astrophysical discs exhibiting vertical shear, which arises when there is a radial variation in the temperature or entropy. Such discs are subject to a "vertical-shear instability", which recent nonlinear simulations have shown to drive hydrodynamic activity in the MRI-stable regions of protoplanetary discs. We first revisit locally isothermal discs using the quasi-global reduced model derived by Nelson et al. (2013). This analysis is then extended to global axisymmetric perturbations in a cylindrical domain. We also derive and study a reduced model describing discs with power law radial entropy profiles ("locally polytropic discs"), which are somewhat more realistic in that they possess physical (as opposed to numerical) surfaces. In all cases the fastest growing modes have very short wavelengths and are localised at the disc surfaces (if present), where the vertical shear is maximal. An additional class of modestly growing vertically global body modes is excited, corresponding to destabilised classical inertial waves ("r-modes"). We discuss the properties of both types of modes, and stress that those that grow fastest occur on the shortest available length scales (determined either by the numerical grid or the physical viscous length). This ill-posedness makes simulations of the instability difficult to interpret. We end with some brief speculation on the nonlinear saturation and resulting angular momentum transport.
In this paper, we estimate the X-ray emission from close-in exoplanets. We show that the Solar/Stellar Wind Charge Exchange Mechanism (SWCX) which produces soft X-ray emission is very effective for hot Jupiters. In this mechanism, X-ray photons are emitted as a result of the charge exchange between heavy ions in the solar wind and the atmospheric neutral particles. In the Solar System, comets produce X-rays mostly through the SWCX mechanism, but it has also been shown to operate in the heliosphere, in the terrestrial magnetosheath, and on Mars, Venus and Moon. Since the number of emitted photons is proportional to the solar wind mass flux, this mechanism is not very effective for the Solar system giants. Here we present a simple estimate of the X-ray emission intensity that can be produced by close-in extrasolar giant planets due to charge exchange with the heavy ions of the stellar wind. Using the example of HD~209458b, we show that this mechanism alone can be responsible for an X-ray emission of $\approx 10^{22}$~erg~s$^{-1}$, which is $10^6$ times stronger than the emission from the Jovian aurora. We discuss also the possibility to observe the predicted soft X-ray flux of hot Jupiters and show that despite high emission intensities they are unobservable with current facilities.
The Telescope Array (TA) is an experiment to observe Ultra-High Energy Cosmic Rays (UHECRs). TA's recent results, the energy spectrum and anisotropy based on the 6-year surface array data, and the primary composition obtained from the shower maximum Xmax are reported. The spectrum demonstrates a clear dip and cutoff. The shape of the spectrum is well described by the energy loss of extra-galactic protons interacting with the cosmic microwave background (CMB). Above the cutoff, a medium-scale (20 degrees radius) flux enhancement was observed near the Ursa-Major. A chance probability of creating this hotspot from the isotropic flux is 4.0 sigma. The measured Xmax is consistent with the primary being proton or light nuclei for energies 10^18.2 eV - 10^19.2 eV.
We model the full-Stokes spectrum of the brown dwarf LSR J1835+3259 in the bands of the diatomic molecules CrH, TiO, and FeH in order to infer its magnetic properties. The models are then compared to the observational data obtained with the Low Resolution Imaging Spectrograph (LRISp) at the Keck observatory. Our preliminary analysis shows that the brown dwarf considered possesses a magnetic field of the order of 2-3 kG.
The VISTA survey of the Magellanic Clouds system (VMC) is a public survey project of the European Southern Observatory. It is collecting multi-band near-infrared data across large areas of the Large and Small Magellanic Clouds, the Magellanic Bridge and a few fields in the Magellanic Stream. The combination of great sensitivity to stars below the old main sequence turn off, and the multiplicity at Ks band (at least 12 epochs) make VMC highly suitable for the determination of the spatially resolved star formation history and three-dimensional geometry, using variable stars such as Cepheids and RR Lyrae stars. The VMC observations are progressing well and will be completed in 2018. The VMC survey has a high legacy value and many science results have already been published, e.g. the most detailed star formation history map of the Small Magellanic Cloud, and others are in preparation, e.g. a comprehensive investigation of classical Cepheids and a study of the proper motion in the foreground 47 Tuc cluster.
We identified three samples of ROSAT sources with Sloan Digital Sky Survey spectra, one at redshift z = 0.1, a second one at z = 0.2 and a third one at z = 0.3. 812 sources in total were examined. We determined the nature and spectral sub-types of the sources by visual inspection. The fraction of each sub-type at each of the three redshifts are then calculated. We consider selection biases caused by the luminosity cut-off threshold to determine whether any systematic trends in AGN type are evident with increasing redshift. We hence probe if an evolution effect is detected in our sample.
We report the discovery of a group of apparently young CoRoT red-giant stars exhibiting enhanced [alpha/Fe] abundance ratios (as determined from APOGEE spectra) with respect to Solar values. Their existence is not explained by standard chemical evolution models of the Milky Way, and shows that the chemical-enrichment history of the Galactic disc is more complex. We find similar stars in previously published samples for which isochrone-ages could be robustly obtained, although in smaller relative numbers, which could explain why these stars have not received prior attention. The young [alpha/Fe]-rich stars are much more numerous in the CoRoT-APOGEE (CoRoGEE) inner-field sample than in any other high-resolution sample available at present, as only CoRoGEE can explore the inner-disc regions and provide ages for its field stars. The kinematic properties of the young [$\alpha$/Fe]-rich stars are not clearly thick-disc like, despite their rather large distances from the Galactic mid-plane. Our tentative interpretation of these and previous intriguing observations in the Milky Way is that these stars were formed close to the end of the Galactic bar, near corotation -- a region where gas can be kept inert for longer times, compared to other regions shocked more frequently by the passage of spiral arms. Moreover, that is where the mass return from older inner-disc stellar generations should be maximal (according to an inside-out disc-formation scenario), further diluting the in-situ gas. Other possibilities to explain these observations (e.g., a recent gas-accretion event) are also discussed.
We present a Monte Carlo radiative transfer technique for calculating synthetic spectropolarimetry for multi-dimensional supernova explosion models. The approach utilises "virtual-packets" that are generated during the propagation of the Monte Carlo quanta and used to compute synthetic observables for specific observer orientations. Compared to extracting synthetic observables by direct binning of emergent Monte Carlo quanta, this virtual-packet approach leads to a substantial reduction in the Monte Carlo noise. This is vital for calculating synthetic spectropolarimetry (since the degree of polarisation is typically very small) but also useful for calculations of light curves and spectra. We first validate our approach via application of an idealised test code to simple geometries. We then describe its implementation in the Monte Carlo radiative transfer code ARTIS and present test calculations for simple models for Type Ia supernovae. Specifically, we use the well-known one-dimensional W7 model to verify that our scheme can accurately recover zero polarisation from a spherical model, and to demonstrate the reduction in Monte Carlo noise compared to a simple packet-binning approach. To investigate the impact of aspherical ejecta on the polarisation spectra, we then use ARTIS to calculate synthetic observables for prolate and oblate ellipsoidal models with Type Ia supernova compositions.
Since 2009 we are running a project to identify flashes produced by the impact of meteoroids on the surface of the Moon. For this purpose we are employing small telescopes and high-sensitivity CCD video cameras. To automatically identify these events a software package called MIDAS was developed and tested. This package can also perform the photometric analysis of these flashes and estimate the value of the luminous efficiency. Besides, we have implemented in MIDAS a new method to establish which is the likely source of the meteoroids (known meteoroid stream or sporadic background). The main features of this computer program are analyzed here, and some examples of lunar impact events are presented.
Limb-darkening is fundamental in determining transit lightcurve shapes, and is typically modeled by a variety of laws that parametrize the intensity profile of the star that is being transited. Confronted with a transit lightcurve, some authors fix the parameters of these laws, the so-called limb-darkening coefficients (LDCs), while others prefer to let them float in the lightcurve fitting procedure. Which of these is the best strategy, however, is still unclear, as well as how and by how much each of these can bias the retrieved transit parameters. In this work we attempt to clarify those points by first re-calculating these LDCs, comparing them to measured values from Kepler transit lightcurves using an algorithm that takes into account uncertainties in both the geometry of the transit and the parameters of the stellar host. We show there are significant departures from predicted model values, suggesting that our understanding of limb-darkening still needs to improve. Then, we show through simulations that if one uses the quadratic limb-darkening law to parametrize limb-darkening, fixing and fitting the LDCs can lead to significant biases -up to $\sim 3\%$ and $\sim 1\%$ in $R_p/R_*$, respectively-, which are important for several confirmed and candidate exoplanets. We conclude that, in this case, the best approach is to let the LDCs be free in the fitting procedure. Strategies to avoid biases in data from present and future missions involving high precision measurements of transit parameters are described.
The model is presented of H$_2$O maser in the $1_{10}-1_{01}$ line within the first excited vibrational state of the molecule around oxygen-rich asymptotic giant branch stars. It is suggested that the maser cloud is located in the inner layers of the circumstellar envelope where intense dust formation takes place. The calculations took into account rotational levels belonging to the five lowest vibrational states of the H$_2$O molecule. The model predicts the gain values of the 658-GHz maser about $10^{-14} -10^{-13}$ cm$^{-1}$ at H$_2$ molecule concentrations $10^9 - 10^{11}$ cm$^{-3}$ and at high ortho-H$_2$O concentrations $ \gtrsim 10^5$ cm$^{-3}$. The gas temperatures $1000 - 1500$ K are considered to be a necessary condition for the effective maser operation. Results are presented for other maser transitions of the excited vibrational states of the molecule.
(abridged) The heating mechanism at high densities during M dwarf flares is poorly understood. Spectra of M dwarf flares in the optical and near-ultraviolet wavelength regimes have revealed three continuum components during the impulsive phase: 1) an energetically dominant blackbody component with a color temperature of T $\sim$ 10,000 K in the blue-optical, 2) a smaller amount of Balmer continuum emission in the near-ultraviolet at lambda $<$ 3646 Angstroms and 3) an apparent pseudo-continuum of blended high-order Balmer lines. These properties are not reproduced by models that employ a typical "solar-type" flare heating level in nonthermal electrons, and therefore our understanding of these spectra is limited to a phenomenological interpretation. We present a new 1D radiative-hydrodynamic model of an M dwarf flare from precipitating nonthermal electrons with a large energy flux of $10^{13}$ erg cm$^{-2}$ s$^{-1}$. The simulation produces bright continuum emission from a dense, hot chromospheric condensation. For the first time, the observed color temperature and Balmer jump ratio are produced self-consistently in a radiative-hydrodynamic flare model. We find that a T $\sim$ 10,000 K blackbody-like continuum component and a small Balmer jump ratio result from optically thick Balmer and Paschen recombination radiation, and thus the properties of the flux spectrum are caused by blue light escaping over a larger physical depth range compared to red and near-ultraviolet light. To model the near-ultraviolet pseudo-continuum previously attributed to overlapping Balmer lines, we include the extra Balmer continuum opacity from Landau-Zener transitions that result from merged, high order energy levels of hydrogen in a dense, partially ionized atmosphere. This reveals a new diagnostic of ambient charge density in the densest regions of the atmosphere that are heated during dMe and solar flares.
Measuring the morphological parameters of galaxies is a key requirement for
studying their formation and evolution. Surveys such as the Sloan Digital Sky
Survey (SDSS) have resulted in the availability of very large collections of
images, which have permitted population-wide analyses of galaxy morphology.
Morphological analysis has traditionally been carried out mostly via visual
inspection by trained experts, which is time-consuming and does not scale to
large ($\gtrsim10^4$) numbers of images.
Although attempts have been made to build automated classification systems,
these have not been able to achieve the desired level of accuracy. The Galaxy
Zoo project successfully applied a crowdsourcing strategy, inviting online
users to classify images by answering a series of questions. Unfortunately,
even this approach does not scale well enough to keep up with the increasing
availability of galaxy images.
We present a deep neural network model for galaxy morphology classification
which exploits translational and rotational symmetry. It was developed in the
context of the Galaxy Challenge, an international competition to build the best
model for morphology classification based on annotated images from the Galaxy
Zoo project.
For images with high agreement among the Galaxy Zoo participants, our model
is able to reproduce their consensus with near-perfect accuracy ($> 99\%$) for
most questions. Confident model predictions are highly accurate, which makes
the model suitable for filtering large collections of images and forwarding
challenging images to experts for manual annotation. This approach greatly
reduces the experts' workload without affecting accuracy. The application of
these algorithms to larger sets of training data will be critical for analysing
results from future surveys such as the LSST.
The equilibria formed by the self-gravitating, collisionless collapse of simple initial conditions have been investigated for decades. We present the results of our attempts to describe the equilibria formed in $N$-body simulations using thermodynamically-motivated models. Previous work has suggested that it is possible to define distribution functions for such systems that describe maximum entropy states. These distribution functions are used to create radial density and velocity distributions for comparison to those from simulations. A wide variety of $N$-body code conditions are used to reduce the chance that results are biased by numerical issues. We find that a subset of initial conditions studied lead to equilibria that can be accurately described by these models, and that direct calculation of the entropy shows maximum values being achieved.
In this review I discuss the ultra-high energy neutrinos (UHEN) originated
from Cosmic-Rays propogation (GZK neutrinos) and from Gamma Ray Bursts (GRBs),
and discuss their detectability in kilometers scale detectors like ARA and
IceCube.
While GZK neutrinos are expected from cosmic ray interactions on the CMB, the
GRB neutrinos depend on the physics inside the sources. GRBs are predicted to
emit UHEN in the prompt and in the later 'after-glow' phase.
I discuss the constraints on the hadronic component of GRBs derived from the
search of four years of IceCube data for a prompt neutrino fux from gamma-ray
bursts (GRBs) and more in general I present the results of the search for
high-energy neutrinos interacting within the IceCube detector between 2010 and
2013.
We report the discovery of an extremely close white dwarf plus F dwarf main-sequence star in a 12 hour binary identified by combining data from the RAdial Velocity Experiment (RAVE) survey and the Galaxy Evolution Explorer (GALEX) survey. A combination of spectral energy distribution fitting and optical and Hubble Space Telescope ultraviolet spectroscopy allowed us to place fairly precise constraints on the physical parameters of the binary. The system, TYC 6760-497-1, consists of a hot Teff~21,500K, M~0.65Ms white dwarf and an F8 star (M~1.23Ms, R~1.35Rs) seen at a low inclination (i~35 deg). The system is likely the descendent of a binary that contained the F star and a ~2Ms A-type star that filled its Roche-lobe on the second asymptotic giant branch, initiating a common envelope phase. The F star is extremely close to Roche-lobe filling and there is likely to be a short phase of thermal timescale mass-transfer onto the white dwarf. During this phase it will grow in mass by up to 20 per cent, until the mass ratio reaches close to unity, at which point it will appear as a standard cataclysmic variable star. Therefore, TYC 6760-497-1 is the first known progenitor of a super-soft source system, but will not undergo a supernova Ia explosion. Once an accurate distance to the system is determined by Gaia, we will be able to place very tight constraints on the stellar and binary parameters.
We derive the particle asymmetry due to inflationary baryogenesis involving a complex inflaton, obtaining a different result to that in the literature. While asymmetries are found to be significantly smaller than previously calculated, in certain parameter regions baryogenesis can still be achieved.
We consider a minimally coupled scalar field with a monomial potential and a perfect fluid in flat FLRW cosmology. We apply local and global dynamical systems techniques to a new three-dimensional dynamical systems reformulation of the field equations on a compact state space. This leads to a visual global description of the solution space and asymptotic behavior. At late times we employ averaging techniques to prove statements about how the relationship between the equation of state of the fluid and the monomial exponent of the scalar field affects asymptotic source dominance and asymptotic manifest self-similarity breaking. We also situate the `attractor' solution in the three-dimensional state space and show that it corresponds to the one-dimensional unstable center manifold of a de Sitter fixed point, located on an unphysical boundary associated with the dynamics at early times. By deriving a center manifold expansion we obtain approximate expressions for the attractor solution. We subsequently improve the accuracy and range of the approximation by means of Pad\'e approximants and compare with the slow-roll approximation.
The next-to-minimal supersymmetric standard model predicts the formation of domain walls due to the spontaneous breaking of the discrete $Z_3$-symmetry at the electroweak phase transition, and they collapse before the epoch of big bang nucleosynthesis if there exists a small bias term in the potential which explicitly breaks the discrete symmetry. Signatures of gravitational waves produced from these unstable domain walls are estimated and their parameter dependence is investigated. It is shown that the amplitude of gravitational waves becomes generically large in the decoupling limit, and that their frequency is low enough to be probed in future pulsar timing observations.
We consider the motion of classical spinning test particles in Schwarzschild and Kerr metrics and investigate innermost stable circular orbits (ISCO). The main goal of this work is to find analytically the small-spin corrections for the parameters of ISCO (radius, total angular momentum, energy) of spinning test particles in the case of vectors of black hole spin, particle spin and orbital angular momentum being collinear to each other. We analytically derive the small-spin linear corrections for arbitrary Kerr parameter $a$. The cases of Schwarzschild, slowly rotating and extreme Kerr black hole are considered in details. For a slowly rotating black hole the ISCO parameters are obtained up to quadratic in $a$ and particle's spin $s$ terms. From the formulae obtained it is seen that the spin-orbital coupling has attractive character when spin and angular momentum are parallel and repulsive when they are antiparallel. For the case of the extreme Kerr black hole with co-rotating particle we succeed to find the exact (on spin) analytical solution for the limiting ISCO parameters. It has been shown that the limiting value of ISCO radius does not depend on the particle's spin. We have also considered circular orbits of arbitrary radius and have found small-spin linear corrections for the total angular momentum and energy at given radius. System of equations for numerical calculation of ISCO parameters for arbitrary $a$ and $s$ is also explicitly written.
The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing Cosmic Ray Observatory (Tibet, People's Republic of China, 4300 m a.s.l.). The high altitude, combined with the full-coverage technique, allows the detection of extensive air showers in a wide energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum down to the TeV region, where direct balloon/space-borne measurements are available. The detector has been in stable data taking in its full configuration from November 2007 to February 2013. In this paper the measurement of the cosmic ray proton plus helium energy spectrum is presented in the region 3-300 TeV by analyzing the full collected data sample. The resulting spectral index is $\gamma = -2.64 \pm 0.01$. These results demonstrate the possibility of performing an accurate measurement of the spectrum of light elements with a ground based air shower detector.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
We present Keck-I MOSFIRE spectroscopy in the Y and H bands of GDN-8231, a massive, compact, star-forming galaxy (SFG) at a redshift $z\sim1.7$. Its spectrum reveals both H$_{\alpha}$ and [NII] emission lines and strong Balmer absorption lines. The H$_{\alpha}$ and Spitzer MIPS 24 $\mu$m fluxes are both weak, thus indicating a low star formation rate of SFR $\lesssim5-10$ M$_{\odot}$ yr$^{-1}$. This, added to a relatively young age of $\sim700$ Myr measured from the absorption lines, provides the first direct evidence for a distant galaxy being caught in the act of rapidly shutting down its star formation. Such quenching allows GDN-8231 to become a compact, quiescent galaxy, similar to 3 other galaxies in our sample, by $z\sim1.5$. Moreover, the color profile of GDN-8231 shows a bluer center, consistent with the predictions of recent simulations for an early phase of inside-out quenching. Its line-of-sight velocity dispersion for the gas, $\sigma^{\rm{gas}}_{\!_{\rm LOS}}=127\pm32$ km s$^{-1}$, is nearly 40% smaller than that of its stars, $\sigma^{\star}_{\!_{\rm LOS}}=215\pm35$ km s$^{-1}$. High-resolution hydro-simulations of galaxies explain such apparently colder gas kinematics of up to a factor of $\sim1.5$ with rotating disks being viewed at different inclinations and/or centrally concentrated star-forming regions. A clear prediction is that their compact, quiescent descendants preserve some remnant rotation from their star-forming progenitors.
We present the detection of potassium in the atmosphere of HAT-P-1b using optical transit narrowband photometry. The results are obtained using the 10.4 m Gran Telescopio Canarias (GTC) together with the OSIRIS instrument in tunable filter imaging mode. We observed four transits, two at continuum wavelengths outside the potassium feature, at 6792 {\AA} and 8844 {\AA}, and two probing the potassium feature in the line wing at 7582.0 {\AA} and the line core at 7664.9 {\AA} using a 12 {\AA} filter width (R~650). The planet-to-star radius ratios in the continuum are found to be $R_{\rm{pl}}/R_{\star}$ = 0.1176 $\pm$ 0.0013 at 6792 {\AA} and $R_{\rm{pl}}/R_{\star}$ = 0.1168 $\pm$ 0.0022 at 8844 {\AA}, significantly lower than the two observations in the potassium line: $R_{\rm{pl}}/R_{\star}$ = 0.1248 $\pm$ 0.0014 in the line wing at 7582.0 {\AA} and $R_{\rm{pl}}/R_{\star}$ = 0.1268 $\pm$ 0.0012 in the line core at 7664.9 {\AA}. With a weighted mean of the observations outside the potassium feature $R_{\rm{pl}}/R_{\star}$ = 0.1174 $\pm$ 0.0010, the potassium is detected as an increase in the radius ratio of {\Delta}$R_{\rm{pl}}/R_{\star}$ = 0.0073 $\pm$ 0.0017 at 7582.0 {\AA} and {\Delta}$R_{\rm{pl}}/R_{\star}$ = 0.0094 $\pm$ 0.0016 at 7664.9 {\AA} (a significance of 4.3 and 6.1 $\sigma$ respectively). We hypothesise that the strong detection of potassium is caused by a large scale height, which can be explained by a high-temperature at the base of the upper atmosphere. A lower mean molecular mass caused by the dissociation of molecular hydrogen into atomic hydrogen by the EUV flux from the host star may also partly explain the amplitude of our detection.
We have carried out a joint photometric and structural analysis of red sequence galaxies in four clusters at a mean redshift of z ~ 1.25 using optical and near-IR HST imaging reaching to at least 3 magnitudes fainter than $M^*$. As expected, the photometry and overall galaxy sizes imply purely passive evolution of stellar populations in red sequence cluster galaxies. However, the morphologies of red sequence cluster galaxies at these redshifts show significant differences to those of local counterparts. Apart from the most massive galaxies, the high redshift red sequence galaxies are significantly diskier than their low redshift analogues. These galaxies also show significant colour gradients, again not present in their low redshift equivalents, most straightforwardly explained by radial age gradients. A clear implication of these findings is that red sequence cluster galaxies originally arrive on the sequence as disk-dominated galaxies whose disks subsequently fade or evolve secularly to end up as high S\'ersic index early-type galaxies (classical S0s or possibly ellipticals) at lower redshift. The apparent lack of growth seen in a comparison of high and low redshift red sequence galaxies implies that any evolution is internal and is unlikely to involve significant mergers. While significant star formation may have ended at high redshift, the cluster red sequence population continues to evolve (morphologically) for several Gyrs thereafter.
We investigate the recent claim of 'photon underproduction crisis' by Kollmeier et al. (2014) which suggests that the known sources of ultra-violet (UV) radiation may not be sufficient to generate the inferred hydrogen photoionization rate ($\Gamma_{\rm HI}$) in the low redshift inter-galactic medium. Using the updated QSO emissivities from the recent studies and our radiative transfer code developed to estimate the UV background, we show that the QSO contribution to $\Gamma_{\rm HI}$ is higher by a factor ~2 as compared to the previous estimates. Using self-consistently computed combinations of star formation rate density and dust attenuation, we show that a typical UV escape fraction of 4% from star forming galaxies should be sufficient to explain the inferred $\Gamma_{\rm HI}$ by Kollmeier et al. (2014). Interestingly, we find that the contribution from QSOs alone can explain the recently inferred $\Gamma_{\rm HI}$ by Shull et al. (2015) which used the same observational data but different simulation. Therefore, we conclude that the crisis is not as severe as it was perceived before and there seems no need to look for alternate explanations such as low luminosity hidden QSOs or decaying dark matter particles.
Telescope networks are gaining traction due to their promise of higher resource utilization than single telescopes and as enablers of novel astronomical observation modes. However, as telescope network sizes increase, the possibility of scheduling them completely or even semi-manually disappears. In an earlier paper, a step towards software telescope scheduling was made with the specification of the Reservation formalism, through the use of which astronomers can express their complex observation needs and preferences. In this paper we build on that work. We present a solution to the discretized version of the problem of scheduling a telescope network. We derive a solvable integer linear programming (ILP) model based on the Reservation formalism. We show computational results verifying its correctness, and confirm that our Gurobi-based implementation can address problems of realistic size. Finally, we extend the ILP model to also handle the novel observation requests that can be specified using the more advanced Compound Reservation formalism.
We study dynamical capture binary black hole-neutron star mergers focusing on the effects of the neutron star spin. These events may arise in dense stellar regions, such as globular clusters, where the majority of neutron stars are expected to be rapidly rotating. Using general-relativistic hydrodynamical simulations, we find that even moderate spins can significantly increase the amount of mass in unbound material. In some of the more extreme cases, there can be up to a third of a solar mass in unbound matter. Similarly, large amounts of tidally stripped material can remain bound and eventually accrete onto the black hole --- as much as a tenth of a solar mass in some cases. These simulations demonstrate that it is important to treat neutron star spin in order to make reliable predictions of the gravitational wave and electromagnetic transient signals accompanying these sources.
Dust formation in supernova ejecta is currently the leading candidate to explain the large quantities of dust observed in the distant, early Universe. However, it is unclear whether the ejecta-formed dust can survive the hot interior of the supernova remnant (SNR). We present infrared observations of ~0.02 $M_\odot$ of warm (~100 K) dust seen near the center of the ~10,000 yr-old Sgr A East SNR at the Galactic center. Our findings signify the detection of dust within an older SNR that is expanding into a relatively dense surrounding medium ($n_e$ ~ 100 $\mathrm{cm}^{-3}$) and has survived the passage of the reverse shock. The results suggest that supernovae may indeed be the dominant dust production mechanism in the dense environment of early Universe galaxies.
Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are recognized as the most promising magnetar candidates as indicated by their energetic bursts and rapid spin-downs. It is expected that the strong magnetic field leaves distinctive imprints on the emergent radiation both by affecting the radiative processes in atmospheres of magnetars and by scattering in the upper magnetospheres. We construct a self-consistent physical model that incorporates emission from the magnetar surface and its reprocessing in the three-dimensional (3D) twisted magnetosphere using a Monte Carlo technique. The synthetic spectra are characterized by four parameters: surface temperature $kT$, surface magnetic field strength $B$, magnetospheric twist angle $\Delta\phi$, and the normalized electron velocity $\beta$. We also create a tabular model (STEMS3D) and apply it to a large sample of XMM-Newton spectra of magnetars. The model successfully fits nearly all spectra, and the obtained magnetic field for the 7 out of 11 sources are consistent with the values inferred from the spin-down rates. We conclude that the continuum-fitting with our model is a robust method to measure the magnetic field strength and magnetospheric configuration of AXPs and SGRs. Investigating the multiple observations of variable sources, we also study the mechanism of their spectral evolution. Our results suggest that the magnetospheres in these sources are highly twisted ($\Delta\phi > 1$), and the behavior of magnetospheric twisting and untwisting is revealed in the 2002 outburst of 1E 2259+586.
A powerful method to measure the mass profile of a galaxy is through the velocities of tracer particles distributed through its halo. Transforming this kind of data accurately to a mass profile M(r), however, is not a trivial problem. In particular, limited or incomplete data may substantially affect the analysis. In this paper we develop a Bayesian method to deal with incomplete data effectively; we have a hybrid-Gibbs sampler that treats the unknown velocity components of tracers as parameters in the model. We explore the effectiveness of our model using simulated data, and then apply our method to the Milky Way using velocity and position data from globular clusters and dwarf galaxies. We find that in general, missing velocity components have little effect on the total mass estimate. However, the results are quite sensitive to the outer globular cluster Pal 3. Using a basic Hernquist model with an isotropic velocity dispersion, we obtain credible regions for the cumulative mass profile M(r) of the Milky Way, and provide estimates for the model parameters with 95 percent Bayesian credible intervals. The mass contained within 260 kpc is 1.37x10^12 solar masses, with a 95 percent credible interval of (1.27,1.51)x10^12 solar masses. The Hernquist parameters for the total mass and scale radius are 1.55 (+0.18/-0.13)x10^12 solar masses and 16.9 (+4.8/-4.1) kpc, where the uncertainties span the 95 percent credible intervals. The code we developed for this work, Galactic Mass Estimator (GME), will be available as an open source package in the R Project for Statistical Computing.
We demonstrate a new procedure to derive accurate and precise surface gravities from high resolution spectra without the use of external constraints. Our analysis utilizes Spectroscopy Made Easy (SME) with robust spectral line constraints and uses an iterative process to mitigate degeneracies in the fitting process. We adopt an updated radiative transfer code, a new treatment for neutral perturber broadening, a line list with multiple gravity constraints and separate fitting for global stellar properties and abundance determinations. To investigate the sources of temperature dependent trends in determining log g noted in previous studies, we obtained Keck HIRES spectra of 42 Kepler asteroseismic stars. In comparison to asteroseismically determined log g our spectroscopic analysis has a constant offset of 0.01 dex with a root mean square (RMS) scatter of 0.05 dex. We also analyzed 30 spectra which had published surface gravities determined using the $a/R_*$ technique from planetary transits and found a constant offset of 0.06 dex and RMS scatter of 0.07 dex. The two samples covered effective temperatures between 5000K and 6700K with log g between 3.7 and 4.6.
We present an initial study of the mass and evolutionary state of a massive and distant cluster, RCS2 J232727.6-020437. This cluster, at z=0.6986, is the richest cluster discovered in the RCS2 project. The mass measurements presented in this paper are derived from all possible mass proxies: X-ray measurements, weak-lensing shear, strong lensing, Sunyaev Zel'dovich effect decrement, the velocity distribution of cluster member galaxies, and galaxy richness. While each of these observables probe the mass of the cluster at a different radius, they all indicate that RCS2 J232727.6-020437 is among the most massive clusters at this redshift, with an estimated mass of M_200 ~3 x10^15 h^-1 Msun. In this paper, we demonstrate that the various observables are all reasonably consistent with each other to within their uncertainties. RCS2 J232727.6-020437 appears to be well relaxed -- with circular and concentric X-ray isophotes, with a cool core, and no indication of significant substructure in extensive galaxy velocity data.
The transiting hot Jupiter planet Qatar-1 b was presented to exhibit variations in transit times that could be of perturbative nature. A hot Jupiter with a planetary companion on a nearby orbit would constitute an unprecedented planetary configuration, important for theories of formation and evolution of planetary systems. We performed a photometric follow-up campaign to confirm or refute transit timing variations. We extend the baseline of transit observations by acquiring 18 new transit light curves acquired with 0.6-2.0 m telescopes. These photometric time series, together with data available in the literature, were analyzed in a homogenous way to derive reliable transit parameters and their uncertainties. We show that the dataset of transit times is consistent with a linear ephemeris leaving no hint for any periodic variations with a range of 1 min. We find no compelling evidence for the existence of a close-in planetary companion to Qatar-1 b. This finding is in line with a paradigm that hot Jupiters are not components of compact multi-planetary systems. Based on dynamical simulations, we place tighter constraints on a mass of any fictitious nearby planet in the system. Furthermore, new transit light curves allowed us to redetermine system parameters with the precision better than that reported in previous studies. Our values generally agree with previous determinations.
Recent observations have discovered the presence of a Box/Peanut or X-shape structure in the Galactic bulge. Such Box/Peanut structures are common in external disc galaxies, and are well-known in N-body simulations where they form following the buckling instability of a bar. From studies of analytical potentials and N-body models it has been claimed in the past that Box/Peanut bulges are supported by "bananas", or x1v1 orbits. We present here a set of N-body models where instead the peanut bulge is mainly supported by brezel-like orbits, allowing strong peanuts to form with short extent relative to the bar length. This shows that stars in the X-shape do not necessarily stream along banana orbits which follow the arms of the X-shape. The brezel orbits are also found to be the main orbital component supporting the peanut shape in our recent Made-to-Measure dynamical models of the Galactic bulge. We also show that in these models the fraction of stellar orbits that contribute to the X-structure account for 40-45% of the stellar mass.
Sulfur appears to be depleted by an order of magnitude or more from its elemental abundance in star-forming regions. In the last few years, numerous observations and experiments have been performed in order to to understand the reasons behind this depletion without providing a satisfactory explanation of the sulfur chemistry towards high-mass star-forming cores. Several sulfur-bearing molecules have been observed in these regions, and yet none are abundant enough to make up the gas-phase deficit. Where, then, does this hidden sulfur reside? This paper represents a step forward in our understanding of the interactions among the various S-bearing species. We have incorporated recent experimental and theoretical data into a chemical model of a hot molecular core in order to see whether they give any indication of the identity of the sulfur sink in these dense regions. Despite our model producing reasonable agreement with both solid-phase and gas-phase abundances of many sulfur-bearing species, we find that the sulfur residue detected in recent experiments takes up only ~6 per cent of the available sulfur in our simulations, rather than dominating the sulfur budget.
The first stars continue to elude modern telescopes, but much has been accomplished in observing the glow of the first galaxies. As detection capabilities improve we will eventually resolve these galaxies, but hopes of observing an individual star remains dim for the foreseeable future. However, our first view of an individual first star may be possible due to its explosion. In this work, we present evolution calculations for Population III (Pop III) stars and their subsequent supernovae explosions. Our evolution models include a mass range of 15 - 100 M$_\odot$, each with initial heavy element abundance Z = 10-14. Our models are evolved from pre-main sequence through formation of an iron core, and thus near to core collapse. We find that modelling the evolution of these stars is very sensitive to the choice of convection criterion; here we provide evolution results using both the Schwarzschild and Ledoux criteria. We also use the final structure from our models for numerical simulation of their supernovae light curves using a radiation hydrodynamics code. In doing so, we estimate a lower bound of initial model mass that may be possible to observe in the near future. We also find that our 40 M$_\odot$ Schwarzschild evolution model produces the brightest supernova peak and statistically should be the most frequently observed. At our highest redshift z= 15, only the 60 M$_\odot$ Schwarzschild model at peak magnitude starts to rival the 40 M$_\odot$ model in brightness.
Whereas it is generally accepted that calcium-aluminum-rich inclusions (CAIs) from chondritic meteorites formed in a hot environment in the solar protoplanetary disk, the conditions of their formation remain debated. Recent laboratory studies of CAIs have provided new kind of data: their size distributions. We show that size distributions of CAIs measured in laboratory from sections of carbonaceous chondrites have a power law size distribution with cumulative size exponent between -1.7 and -1.9, which translates into cumulative size exponent between -2.5 and -2.8 after correction for sectioning. To explain these observations, numerical simulations were run to explore the growth of CAIs from micrometer to centimeter sizes, in a hot and turbulent protoplanetary disk through the competition of coagulation and fragmentation. We show that the size distributions obtained in growth simulations are in agreement with CAIs size distributions in meteorites. We explain the CAI sharp cut-off of their size distribution at centimeter sizes as the direct result from the famous fragmentation barrier, provided that CAI fragment for impact velocities larger than 10 m/s. The growth/destruction timescales of millimeter- and centimeter-sized CAIs is inversely proportional to the local dust/gas ratio and is about 10 years at 1300 K and up to 104 years at 1670K. This implies that the most refractory CAIs are expected to be smaller in size owing to their long growth timescale compared to less refractory CAIs. Conversely, the least refractory CAIs could have been recycled many times during the CAI production era which may have profound consequences for their radiometric age.
Here we analyze radio, optical, and X-ray data for a peculiar cluster Abell 578. This cluster is not fully relaxed and consists of two merging sub-systems. The brightest cluster galaxy, CGPG 0719.8+6704, is a pair of interacting ellipticals with projected separation $\sim$10 kpc, the brighter of which hosts the radio source 4C +67.13. The Fanaroff-Riley type-II radio morphology of 4C +67.13 is unusual for central radio galaxies in local Abell clusters. Our new optical spectroscopy revealed that both nuclei of the CGPG 0719.8+6704 pair are active, albeit at low accretion rates corresponding to the Eddington ratio $\sim10^{-4}$ (for the estimated black hole masses of $\sim 3 \times 10^8\,M_\odot$ and $\sim 10^9 \, M_\odot$). The gathered X-ray ({\it Chandra}) data allowed us to confirm and to quantify robustly the previously noted elongation of the gaseous atmosphere in the dominant sub-cluster, as well as a large spatial offset ($\sim 60$\,kpc projected) between the position of the brightest cluster galaxy and the cluster center inferred from the modeling of the X-ray surface brightness distribution. Detailed analysis of the brightness profiles and temperature revealed also that the cluster gas in the vicinity of 4C\,+67.13 is compressed (by a factor of about $\sim 1.4$) and heated (from $\simeq 2.0$\,keV up to 2.7\,keV), consistent with the presence of a weak shock (Mach number $\sim 1.3$) driven by the expanding jet cocoon. This would then require the jet kinetic power of the order of $\sim 10^{45}$\,erg\,s$^{-1}$, implying either a very high efficiency of the jet production for the current accretion rate, or a highly modulated jet/accretion activity in the system.
The Milky Way bulge offers a unique opportunity to investigate in detail the role that different processes such as dynamical instabilities, hierarchical merging, and dissipational collapse may have played in the history of the Galaxy formation and evolution based on its resolved stellar population properties. Large observation programmes and surveys of the bulge are providing for the first time a look into the global view of the Milky Way bulge that can be compared with the bulges of other galaxies, and be used as a template for detailed comparison with models. The Milky Way has been shown to have a box/peanut (B/P) bulge and recent evidence seems to suggest the presence of an additional spheroidal component. In this review we summarise the global chemical abundances, kinematics and structural properties that allow us to disentangle these multiple components and provide constraints to understand their origin. The investigation of both detailed and global properties of the bulge now provide us with the opportunity to characterise the bulge as observed in models, and to place the mixed component bulge scenario in the general context of external galaxies. When writing this review, we considered the perspectives of researchers working with the Milky Way and researchers working with external galaxies. It is an attempt to approach both communities for a fruitful exchange of ideas.
To characterize the signals registered by the different types of water Cherenkov detectors (WCD) used by the Latin American Giant Observatory (LAGO) Project, it is necessary to develop detailed simulations of the detector response to the flux of secondary particles at the detector level. These particles are originated during the interaction of cosmic rays with the atmosphere. In this context, the LAGO project aims to study the high energy component of gamma rays bursts (GRBs) and space weather phenomena by looking for the solar modulation of galactic cosmic rays (GCRs). Focus in this, a complete and complex chain of simulations is being developed that account for geomagnetic effects, atmospheric reaction and detector response at each LAGO site. In this work we shown the first steps of a GEANT4 based simulation for the LAGO WCD, with emphasis on the induced effects of the detector internal diffusive coating.
We conducted astrometric VLBI observations of water-vapor maser emission in the massive star forming region IRAS 21379+5106 to measure the annual parallax and proper motion, using VERA. The annual parallax was measured to be $0.262 \pm 0.031$ mas corresponding to a trigonometric distance of $3.82^{+0.51}_{-0.41}$ kpc. The proper motion was $(\mu_\alpha\cos{\delta}, \mu_\delta)=(-2.74 \pm 0.08, -2.87 \pm 0.18)$ mas yr$^{-1}$. Using this result, the Galactic rotational velocity was estimated to be $V_\theta=218\pm 19$ km s$^{-1}$ at the Galactocentric distance $R=9.22\pm0.43$ kpc, when we adopted the Galactic constants $R_0=8.05\pm 0.45$ kpc and $V_0=238\pm 14$ km s$^{-1}$. With newly determined distance, {the bolometric luminosity of the central young stellar object was re-evaluated to $(2.15\pm 0.54)\times 10^3 L_\odot$, which corresponds to spectral type of} B2--B3. Maser features were found to be distributed along a straight line from south-west to north-east. In addition, a vector map of the internal motions constructed from the residual proper motions implies that maser features trace a bipolar flow and that it cannot be explained by simple ballistic motion.
The Kepler mission has provided high quality light curves for more than 2000 eclipsing binaries. Tertiary companions to these binaries can be detected if they transit one or both stars in the binary or if they perturb the binary enough to cause deviations in the observed times of the primary and secondary eclipses (in a few cases both effects are observed in the same eclipsing binary). From the study of eclipse timing variations, it is estimated that 15 to 20% of the Kepler eclipsing binaries have close-in tertiary companions. I will give an overview of recent results and discuss some specific systems of interest.
GRAVITY is the second generation Very Large Telescope Interferometer instrument for precision narrow-angle astrometry and interferometric imaging in the Near Infra-Red (NIR). It shall provide precision astrometry of order 10 microarcseconds, and imaging capability at a few milliarcsecond resolution, and hence will revolutionise dynamical measurements of celestial objects. GRAVITY is currently in the last stages of its integration and tests in Garching at MPE, and will be delivered to the VLT Interferometer (VLTI) in 2015. We present here the instrument, with a particular focus on the components making use of fibres: integrated optics beam combiners, polarisation rotators, fibre differential delay lines, and the metrology.
Chondrules are millimeter-sized spherules that dominate primitive meteorites (chondrites) originating from the asteroid belt. The incorporation of chondrules into asteroidal bodies must be an important step in planet formation, but the mechanism is not understood. We show that the main growth of asteroids can result from gas-drag-assisted accretion of chondrules. The largest planetesimals of a population with a characteristic radius of 100 km undergo run-away accretion of chondrules within ~3 Myr, forming planetary embryos up to Mars sizes along with smaller asteroids whose size distribution matches that of main belt asteroids. The aerodynamical accretion leads to size-sorting of chondrules consistent with chondrites. Accretion of mm-sized chondrules and ice particles drives the growth of planetesimals beyond the ice line as well, but the growth time increases above the disk life time outside of 25 AU. The contribution of direct planetesimal accretion to the growth of both asteroids and Kuiper belt objects is minor. In contrast, planetesimal accretion and chondrule accretion play more equal roles for the formation of Moon-sized embryos in the terrestrial planet formation region. These embryos are isolated from each other and accrete planetesimals only at a low rate. However, the continued accretion of chondrules destabilizes the oligarchic configuration and leads to the formation of Mars-sized embryos and terrestrial planets by a combination of direct chondrule accretion and giant impacts.
HH 212 is a nearby (400 pc) highly collimated protostellar jet powered by a Class 0 source in Orion. We have mapped the inner 80" (~ 0.16 pc) of the jet in SiO (J=8-7) and CO (J=3-2) simultaneously at ~ 0.5 resolution with the Atacama Millimeter/Submillimeter Array at unprecedented sensitivity. The jet consists of a chain of knots, bow shocks, and sinuous structures in between. As compared to that seen in our previous observations with the Submillimeter Array, it appears to be more continuous, especially in the northern part. Some of the knots are now seen associated with small bow shocks, with their bow wings curving back to the jet axis, as seen in pulsed jet simulations. Two of them are reasonably resolved, showing kinematics consistent with sideways ejection, possibly tracing the internal working surfaces formed by a temporal variation in the jet velocity. In addition, nested shells are seen in CO around the jet axis connecting to the knots and bow shocks, driven by them. The proper motion of the jet is estimated to be ~ 115+-50 km/s, comparing to our previous observations. The jet has a small semi-periodical wiggle, with a period of ~ 93 yrs. The amplitude of the wiggle first increases with the distance from the central source and then stays roughly constant. One possible origin of the wiggle could be the kink instability in a magnetized jet.
We recently reported the first identification of rotational transitions of methyl formate (HCOOCH$_{3}$) in the second torsionally excited state toward Orion Kleinmann-Low (KL) observed with the Nobeyama 45 m telescope. In combination with the identified transitions of methyl formate in the ground state and the first torsional excited state, it was found that there is a difference in rotational temperature and vibrational temperature, where the latter is higher. In this study, high spatial resolution analysis by using Atacama Large Millimeter/Submillimeter Array (ALMA) science verification data was carried out to verify and understand this difference. Toward the Compact Ridge, two different velocity components at 7.3 and 9.1 km s$^{-1}$ were confirmed, while a single component at 7.3 km s$^{-1}$ was identified towards the Hot Core. The intensity maps in the ground, first, and second torsional excited states have quite similar distributions. Using extensive ALMA data, we determined the rotational and vibrational temperatures for the Compact Ridge and Hot Core by the conventional rotation diagram method. The rotational temperature and vibrational temperatures agree for the Hot Core and for one component of the Compact Ridge. At the 7.3 km s$^{-1}$ velocity component for the Compact Ridge, the rotational temperature was found to be higher than the vibrational temperature. This is different from what we obtained from the results by using the single-dish observation. The difference might be explained by the beam dilution effect of the single-dish data and/or the smaller number of observed transitions within the limited range of energy levels ($\leq$30 K) of $E_u$ in the previous study.
The Planck-HFI all-sky survey from 353 to 857GHz combined with the 100 microns IRAS show that the dust properties vary in the diffuse ISM at high Galactic latitude (1e19<NH<2.5e20 H/cm2). Our aim is to explain these variations with changes in the ISM properties and grain evolution. Our starting point is the latest core-mantle dust model. It consists of small aromatic-rich carbon grains, larger amorphous carbon grains with aliphatic-rich cores and aromatic-rich mantles, and amorphous silicates with Fe/FeS nano-inclusions covered by aromatic-rich carbon mantles. We explore whether variations in the radiation field or in the gas density distribution in the diffuse ISM could explain the observations. The dust properties are also varied in terms of mantle thickness, Fe/FeS inclusions, carbon abundance, and size distribution. Variations in the radiation field intensity and gas density distribution cannot explain the observed variations but radiation fields harder than the standard ISRF may participate in creating part of them. We further show that variations in the grain mantle thickness coupled with changes in the grain size distribution can reproduce most of the observations. We put a limit on the mantle thickness of the silicates (~10-15nm), and find that aromatic-rich mantles are needed for the carbon grains (at least 5-7.5nm thick). We also find that changes in the carbon abundance in the grains could explain part of the observed variations. Finally, we show that varying the composition of Fe/FeS inclusions in the silicates cannot account for the variations. With small variations in the dust properties, we are able to explain most of the variations in the dust emission observed by Planck-HFI in the diffuse ISM. We also find that the small realistic changes in the dust properties that we consider almost perfectly match the anti-correlation and scatter in the observed beta-T relation.
We discuss the conditions to satisfy the Higuchi bound and to avoid gradient instabilities in the scalar sector for cosmological solutions in singly coupled bimetric gravity theories. We find that in expanding universes the ratio of the scale factors of the reference and observable metric has to increase at all times. This automatically implies a ghost-free helicity-2 sector and enforces a phantom Dark Energy. Furthermore, the condition for the absence of gradient instabilities in the scalar sector will be analyzed. Finally, we discuss whether cosmological solutions, including exotic evolutions like bouncing cosmologies, can exist, in which both the Higuchi ghost and scalar instabilities are absent at all times.
One of the most energetic gamma-ray burst GRB 110731A was observed from optical to GeV energy range. Previous analysis on the prompt phase revealed similarities with the Large Area Telescope (LAT) bursts observed by Fermi: i) a delayed onset of the high-energy emission ($> 100$ MeV), ii) a short-lasting bright peak at later times and iii) a temporally extended component from this phase and lasting hundreds of seconds. Additionally to the prompt phase, multiwavelength observations over different epochs showed that the spectral energy distribution was better fitted by a wind afterglow model. We present a leptonic model based on an early afterglow that evolves in a stellar wind of its progenitor. We apply this model to interpret the temporally extended LAT emission and the brightest LAT peak exhibited by the prompt phase of GRB 110731A. Additionally, using the same set of parameters, we describe the multiwavelength afterglow observations. The origin of the temporally extended LAT, X-ray and optical flux is explained through synchrotron radiation from the forward shock and the brightest LAT peak is described evoking the synchrotron self-Compton emission from the reverse shock. The bulk Lorentz factor required in this model ($\Gamma\simeq520$) lies in the range of values demanded for most LAT-detected gamma-ray bursts. We show that the strength of the magnetic field in the reverse-shock region is $\sim$ 50 times stronger than in the forward-shock region. This result suggests that for GRB 110731A, the central engine is likely entrained with strong magnetic fields.
<abridged> NASA's Asteroid Redirect Mission (ARM) has been proposed with the
aim to capture a small asteroid a few meters in size and redirect it into an
orbit around the Moon. There it can be investigated at leisure by astronauts
aboard an Orion or other spacecraft. The target for the mission has not yet
been selected, and there are very few potential targets currently known. Though
sufficiently small near-Earth asteroids (NEAs) are thought to be numerous, they
are also difficult to detect and characterize with current observational
facilities. Here we collect the most up-to-date information on the smallest
known near-Earth asteroids to outline the properties of these small NEAs as
currently understood, in order to examine what the eventual ARM target might be
like. Observational biases certainly mean that our sample is not an ideal
representation of the true population of small NEAs. However our sample is
representative of the eventual target list for the ARM mission, which will be
compiled under very similar observational constraints unless dramatic changes
are made to the way near-Earth asteroids are searched for and studied.
We find that the typical rotation period is 40 minutes. The mean and median
axis ratios were 1.43 and 1.29. Rotation rates much faster than the spin
barrier are seen, reaching below 30 seconds, and implying that most of these
bodies are monoliths. Non-principal axis rotation is uncommon. Axial ratios
often reach values as high as two, though no undisputed results reach above
three. We find little correlation of axis ratio with size. The most common
spectral type in the sample of small NEAs is S-type (> 90%), with only a
handful of C and X types known.
Discovering the mass of neutrinos is a principle goal in high energy physics and cosmology. In addition to cosmological measurements based on two-point statistics, the neutrino mass can also be estimated by observations of neutrino wakes resulting from the relative motion between dark matter and neutrinos. Such a detection relies on an accurate reconstruction of the dark matter-neutrino relative velocity which is affected by non-linear structure growth and galaxy bias. We investigate our ability to reconstruct this relative velocity using large N-body simulations where we evolve neutrinos as distinct particles alongside the dark matter. We find that the dark matter velocity power spectrum is overpredicted by linear theory whereas the neutrino velocity power spectrum is underpredicted. The magnitude of the relative velocity observed in the simulations is found to be lower than what is predicted in linear theory. Since neither the dark matter nor the neutrino velocity fields are directly observable from galaxy or 21 cm surveys, we test the accuracy of a reconstruction algorithm based on halo density fields and linear theory. Assuming prior knowledge of the halo bias, we find that the reconstructed relative velocities are highly correlated with the simulated ones with correlation coefficients of 0.94, 0.93, 0.91 and 0.88 for neutrinos of mass 0.05, 0.1, 0.2 and 0.4 eV. We confirm that the relative velocity field reconstructed from large scale structure observations such as galaxy or 21 cm surveys can be accurate in direction and, with appropriate scaling, magnitude.
We present evolution of galaxy effective radius r_e obtained from the HST samples of ~190,000 galaxies at z=0-10. Our HST samples consist of 176,152 photo-z galaxies at z=0-6 from the 3D-HST+CANDELS catalogue and 10,454 LBGs at z=4-10 identified in CANDELS, HUDF09/12, and HFF parallel fields, providing the largest data set to date for galaxy size evolution studies. We derive r_e with the same technique over the wide-redshift range of z=0-10, evaluating the optical-to-UV morphological K-correction and the selection bias of photo-z galaxies+LBGs as well as the cosmological surface brightness dimming effect. We find that r_e values at a given luminosity significantly decrease towards high-z, regardless of statistics choices. For star-forming galaxies, there is no evolution of the power-law slope of the size-luminosity relation and the median Sersic index (n~1.5). Moreover, the r_e-distribution is well represented by log-normal functions whose standard deviation \sigma_{\ln{r_e}} does not show significant evolution within the range of \sigma_{\ln{r_e}}~0.45-0.75. We calculate the stellar-to-halo size ratio from our r_e measurements and the dark-matter halo masses estimated from the abundance matching study, and obtain a nearly constant value of r_e/r_vir=1.0-3.5% at z=0-8. The combination of the r_e-distribution shape+standard deviation, the constant r_e/r_vir, and n~1.5 suggests a picture that typical high-z star-forming galaxies have disk-like stellar components in a sense of dynamics and morphology over cosmic time of z~0-6. If high-z star-forming galaxies are truly dominated by disks, the r_e/r_vir value and the disk formation model indicate that the specific angular momentum of the disk normalized by the host halo is j_d/m_d=0.5-1. These are statistical results for galaxies' major stellar components, and the detailed study of clumpy sub-components is presented in the paper II.
In the context of single field inflation, models with a quadratic potential and models with a natural potential with subplanckian decay constant are in tension with the Planck data. We show that, when embedded in a two-field model with an additional super massive field, they can become consistent with observations. Our results follow if the inflaton is the phase of a complex field (or an angular variable) protected by a mildly broken U(1) symmetry, and the radial component, whose mass is much greater than the Hubble scale, is stabilized at subplanckian values. The presence of the super massive field, besides modifying the effective single field potential, causes a reduction in the speed of sound of the inflaton fluctuations, which drives the prediction for the primordial spectrum towards the allowed experimental values. We discuss these effects also for the linear potential, and show that this model increases its agreement with data as well
Aims: We study the evolution of stellar rotation and wind properties for
low-mass main-sequence stars. Our aim is to use rotational evolution models to
constrain the mass loss rates in stellar winds and to predict how their
properties evolve with time on the main-sequence.
Methods: We construct a rotational evolution model that is driven by observed
rotational distributions of young stellar clusters. Fitting the free parameters
in our model allows us to predict how wind mass loss rate depends on stellar
mass, radius, and rotation. We couple the results to the wind model developed
in Paper I of this series to predict how wind properties evolve on the
main-sequence.
Results: We estimate that wind mass loss rate scales with stellar parameters
as $\dot{M}_\star \propto R_\star^2 \Omega_\star^{1.33} M_\star^{-3.36}$. We
estimate that at young ages, the solar wind likely had a mass loss rate that is
an order of magnitude higher than that of the current solar wind. This leads to
the wind having a higher density at younger ages; however, the magnitude of
this change depends strongly on how we scale wind temperature. Due to the
spread in rotation rates, young stars show a large range of wind properties at
a given age. This spread in wind properties disappears as the stars age.
Conclusions: There is a large uncertainty in our knowledge of the evolution
of stellar winds on the main-sequence, due both to our lack of knowledge of
stellar winds and the large spread in rotation rates at young ages. Given the
sensitivity of planetary atmospheres to stellar wind and radiation conditions,
these uncertainties can be significant for our understanding of the evolution
of planetary environments.
We describe a new model for the "stripes" of synchrotron radiation seen in the remnant of Tycho's supernova. In our picture, cosmic rays streaming ahead of the forward shock generate parallel (with respect to the local magnetic field direction) circularly polarized Alfven waves that are almost free of dissipation, and due to being circularly polarized exhibit no spatial variation of magnetic field strength. Following interaction with the SNR shock with nonzero obliquity, these parallel waves become obliquely propagating, due the the wave refraction (different in principle for the different plane wave components), and dissipation sets in. The magnetosonic polarization decays faster, due to transit time damping, leaving only the Alfven mode. This surviving mode now exhibits a spatial variation of the magnetic field, leading to local maxima and minima in the synchrotron emission, i.e. the stripes. We attribute the initial wave generation to the Bell instability, which in contrast to the resonant generation of upstream Alfven waves, gives rise to a preferred wavelength, and hence the single wave period at which the stripes are seen. Based on estimates for damping rates due to turbulent cascade and transit time damping, we estimate the dependence of the visibility of the stripes on the shock obliquity, and determine a maximum cosmic ray energy in Tycho's SNR in the range $6\times 10^{14} - 1\times 10^{15}$ eV.
The following is a comment on the recent letter by Iocco et al. (2015, arXiv:1502.03821) where the authors claim to have found "...convincing proof of the existence of dark matter...". The letter in question presents a compilation of recent rotation curve observations for the Milky Way, together with Newtonian rotation curve estimates based on recent baryonic matter distribution measurements. A mismatch between the former and the latter is then presented as "evidence for dark matter". Here we show that the reported discrepancy is the well known gravitational anomaly which consistently appears when dynamical accelerations approach the critical Milgrom acceleration a_0 = 1.2 \times 10^{-10} m / s^2. Further, using a simple modified gravity force law, the baryonic models presented in Iocco et al. (2015), yield dynamics consistent with the observed rotation values.
We describe the new spectroscopic data reduction pipeline for the multi-object MMT/Magellan Infrared Spectrograph. The pipeline is implemented in idl as a stand-alone package and is publicly available in both stable and development versions. We describe novel algorithms for sky subtraction and correction for telluric absorption. We demonstrate that our sky subtraction technique reaches the Poisson limit set by the photon statistics. Our telluric correction uses a hybrid approach by first computing a correction function from an observed stellar spectrum, and then differentially correcting it using a grid of atmosphere transmission models for the target airmass value. The pipeline provides a sufficient level of performance for real time reduction and thus enables data quality control during observations. We reduce an example dataset to demonstrate the high data reduction quality.
Using multi-epoch broad-band observations of the GRB 121024A afterglow, we measure the three characteristic break frequencies of the synchrotron spectrum. We use 6 epochs of combined XRT and GROND data to constrain the temporal slope, the dust extinction and the spectral slope with high accuracy. Two further epochs of combined data from XRT, GROND, APEX, CARMA and EVLA are used to set constraints on the break frequencies and therefore on the micro-physical and dynamical parameters. The XRT and GROND light curves show a simultaneous break at around 42 ks. No spectral evolution is observed between the afterglow SEDs before and after the break. As a result, the crossing of the synchrotron cooling break is not suitable as an explanation for the break in the light curve. The multi-wavelength data give us a unique opportunity to discern between two plausible scenarios explaining the break: the end of energy injection and a jet break. The observations are explained by two possible scenarios, a jet break and an energy injection model. The jet break model has been suggested by previous analysis of the observed linear and circular polarisation although it requires a flat electron spectrum, a very low cooling break and a non-spreading jet. The energy injection avoids an extremely flat spectrum for the shock-accelerated electrons, the very low cooling break frequency and the extreme prompt emission efficiency. However some atypical values for the micro-physics of particle acceleration arise in this model. It is consistent with the correlation between luminosity and end-time of the plateau in the light curve, reportedfrom large sample studies of XRT data.
We set conservative, robust constraints on the annihilation and decay of dark matter into various Standard Model final states under various assumptions about the distribution of the dark matter in the Milky Way halo. We use the inclusive photon spectrum observed by the Fermi Gamma-ray Space Telescope through its main instrument, the Large-Area Telescope (LAT). We use simulated data to first find the "optimal" regions of interest in the gamma-ray sky, where the expected dark matter signal is largest compared with the expected astrophysical foregrounds. We then require the predicted dark matter signal to be less than the observed photon counts in the a priori optimal regions. This yields a very conservative constraint as we do not attempt to model or subtract astrophysical foregrounds. The resulting limits are competitive with other existing limits, and, for some final states with cuspy dark-matter distributions in the Galactic Center region, disfavor the typical cross section required during freeze-out for a weakly interacting massive particle (WIMP) to obtain the observed relic abundance.
We study moduli stabilization in combination with inflation in heterotic orbifold compactifications in the light of a large Hubble scale and the favored tensor-to-scalar ratio $r \approx 0.05$. To account for a trans-Planckian field range we implement aligned natural inflation. Although there is only one universal axion in heterotic constructions, further axions from the geometric moduli can be used for alignment and inflation. We argue that such an alignment is rather generic on orbifolds, since all non-perturbative terms are determined by modular weights of the involved fields and the Dedekind $\eta$ function. We present two setups inspired by the mini-landscape models of the $\mathbb Z_{6-\text{II}}$ orbifold which realize aligned inflation and stabilization of the relevant moduli. One has a supersymmetric vacuum after inflation, while the other includes a gaugino condensate which breaks supersymmetry at a high scale.
We study the left-right symmetric extensions of the Standard Model (LRSM) and point out that the discovery of a right-handed charged gauge boson $W_R^\pm$ with mass in the TeV range will have profound consequences for leptogenesis. We consider the LRSM with both triplet and doublet Higgs scalars, and in both the cases we find that if the $W_R^\pm$ with mass of around a few TeV is found, for example through a signal of two leptons and two jets that has been reported by CMS to have a 2.8$\sigma$ local excess over the Standard Model background at the LHC, then it will rule out all possibilities of leptogenesis. In that case, the baryon asymmetry of the universe has to be generated after the electroweak phase transition and from the baryon number violation that can give rise to neutron-antineutron oscillation or ($B - L$) violating proton decay.
Many modified gravity theories are under consideration in cosmology as the source of the accelerated expansion of the universe and linear perturbation theory, valid on the largest scales, has been examined in many of these models. However, smaller non-linear scales offer a richer phenomenology with which to constrain modified gravity theories. Here, we consider the Hu-Sawicki form of $f(R)$ gravity and apply the post-Friedmann approach to derive the leading order equations for non-linear scales, i.e. the equations valid in the Newtonian-like regime. We reproduce the standard equations for the scalar field, gravitational slip and the modified Poisson equation in a coherent framework. In addition, we derive the equation for the leading order correction to the Newtonian regime, the vector potential. We measure this vector potential from $f(R)$ N-body simulations at redshift zero and one, for two values of the $f_{R_0}$ parameter. We find that the vector potential at redshift zero in $f(R)$ gravity can be close to 50\% larger than in GR on small scales for $|f_{R_0}|=1.289\times10^{-5}$, although this is less for larger scales, earlier times and smaller values of the $f_{R_0}$ parameter. Similarly to in GR, the small amplitude of this vector potential suggests that the Newtonian approximation is highly accurate for $f(R)$ gravity, and also that the non-linear cosmological behaviour of $f(R)$ gravity can be completely described by just the scalar potentials and the $f(R)$ field.
The inversion of temperature at the solar corona is hard to understand from classical physics, and the coronal heating mechanism remains unclear. The heating in the quiet region seems contradicting with the thermodynamics and is a keen problem for physicists. A new mechanism for the coronal heating based on the neutrino radiative transition unique in the corona region is studied. The probability is enormously amplified by an electroweak Chern-Simons form and overlapping waves, and the sufficient energy is transfered. Thus the coronal heating is understood from the quantum effects of the solar neutrino.
The deformability of a compact object induced by a perturbing tidal field is encoded in the tidal Love numbers, which depend sensibly on the object's internal structure. These numbers are known only for static, spherically-symmetric objects. As a first step to compute the tidal Love numbers of a spinning compact star, here we extend powerful perturbative techniques to compute the exterior geometry of a spinning object distorted by an axisymmetric tidal field to second order in the angular momentum. The spin of the object introduces couplings between electric and magnetic deformations and new classes of induced Love numbers emerge. For example, a spinning object immersed in a quadrupolar, electric tidal field can acquire some induced mass, spin, quadrupole, octupole and hexadecapole moments to second order in the spin. The deformations are encoded in a set of inhomogeneous differential equations which, remarkably, can be solved analytically in vacuum. We discuss certain subtleties in defining the multipole moments of the central object, which are due to the difficulty in separating the tidal field from the linear response of the object in the solution. By extending the standard procedure to identify the linear response in the static case, we prove analytically that the Love numbers of a Kerr black hole remain zero to second order in the spin. As a by-product, we provide the explicit form for a slowly-rotating, tidally-deformed Kerr black hole to quadratic order in the spin, and discuss its geodesic and geometrical properties.
f(R) gravity is one of the simplest generalizations of general relativity, which may explain accelerated cosmic expansion without introducing a cosmological constant. Transformed into the Einstein frame, a new scalar degree of freedom appears and it couples with matter fields. In order for f(R) theories to pass the local tests of general relativity, it has been known that the chameleon mechanism with a so-called thin-shell solution must operate. If the thin-shell constraint is applied to a cosmological situation, it has been claimed that the equation of state parameter of dark energy w must be extremely close to -1. We argue this is due to incorrect use of the Poisson's equation which is valid only in the static case. By solving the correct Klein-Gordon equation perturbatively, we show that a thin-shell solution exists even if w deviates from -1 appreciably.
A non-conformally invariant coupling between the inflaton and the photon in the minimal Lorentz-violating standard model extension is analyzed. For specific forms of the Lorentz-violating background tensor, the strong coupling and backreaction problems of magnetogenesis in de Sitter inflation with scale $\sim 10^{16}$GeV are evaded, the electromagnetic-induced primordial spectra of (Gaussian and non-Gaussian) scalar and tensor curvature perturbations are compatible with cosmic microwave background observations, and the inflation-produced magnetic field directly accounts for cosmic magnetic fields.
We study the constraints on dark energy equation of state $\omega^{X}$ and the deceleration parameter $q$ from the recent observational data including Hubble data and the cosmic microwave background (CMB) radiation by using a model-independent deceleration parameter $q(z)=1/2-a/(1+z)^b$ and dark energy equation of state $\omega^{X}=\omega_{0}+\omega_{1}z/(1+z)$ in the scope of anisotropic Bianchi type I space-time. For the cases of Hubble dataset, CMB data, and their combination, our results indicate that the constraints on transition redshift $z_{\ast}$ are $0.62^{+1.45}_{-0.56}$, $0.34^{+0.13}_{-0.06}$, and $0.60^{+0.20}_{-0.10}$ respectively.
For quark matter studies in astrophysics the thermodynamic bag model (tdBAG) has been widely used. Despite its success it fails to account for various phenomena expected from Quantum-Chromo-Dynamics (QCD). We suggest a straightforward extension of tdBAG in order to take the dynamical breaking of chiral symmetry and the influence of vector interactions explicitly into account. As for tdBAG the model mimics confinement in a phenomenological approach. It is based on an analysis of the Nambu--Jona-Lasinio (NJL) model at finite density. Furthermore, we demonstrate how NJL and bag models in this regime follow from the more general and QCD based framework of Dyson-Schwinger (DS) equations in medium by assuming a simple gluon contact interaction. Based on our simple and novel model, we construct quark hadron hybrid equations of state (EoS) and study systematically chiral and deconfinement phase transitions, the appearance of $s$-quarks and the role of vector interaction. We further study these aspects for matter in beta-equilibrium at zero temperature, with particular focus on the current ~2 solar masses maximum mass constraint for neutron stars. Our approach indicates that the currently only theoretical evidence for the hypothesis of stable strange matter is an artifact of tdBAG and results from neglecting the dynamical breaking of chiral symmetry.
We address the problem of noise regression in the output of gravitational-wave (GW) interferometers, using data from the physical environmental monitors (PEM). The objective of the regression analysis is to predict environmental noise in the gravitational-wave channel from the PEM measurements. One of the most promising regression method is based on the construction of Wiener-Kolmogorov filters. Using this method, the seismic noise cancellation from the LIGO GW channel has already been performed. In the presented approach the Wiener-Kolmogorov method has been extended, incorporating banks of Wiener filters in the time-frequency domain, multi-channel analysis and regulation schemes, which greatly enhance the versatility of the regression analysis. Also we presents the first results on regression of the bi-coherent noise in the LIGO data.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
We present a new reduction pipeline for the VIRCAM@VISTA detector and describe the method developed to obtain high-precision astrometry with the VISTA Variables in the V\'ia L\'actea (VVV) data set. We derive an accurate geometric-distortion correction using as calibration field the globular cluster NGC 5139, and showed that we are able to reach a relative astrometric precision of about 8 mas per coordinate per exposure for well-measured stars over a field of view of more than 1 square degree. This geometric-distortion correction is made available to the community. As a test bed, we chose a field centered around the globular cluster NGC 6656 from the VVV archive and computed proper motions for the stars within. With 45 epochs spread over four years, we show that we are able to achieve a precision of 1.4 mas/yr and to isolate each population observed in the field (cluster, Bulge and Disk) using proper motions. We used proper-motion-selected field stars to measure the motion difference between Galactic disk and bulge stars. Our proper-motion measurements are consistent with UCAC4 and PPMXL, though our errors are much smaller. Models have still difficulties in reproducing the observations in this highly-reddened Galactic regions.
We present new wide-field photometry and spectroscopy of the globular clusters (GCs) around NGC 4649 (M60), the third brightest galaxy in the Virgo cluster. Imaging of NGC 4649 was assembled from a recently-obtained HST/ACS mosaic, and new Subaru/Suprime-Cam and archival CFHT/MegaCam data. About 1200 sources were followed up spectroscopically using combined observations from three multi-object spectrographs: Keck/DEIMOS, Gemini/GMOS and MMT/Hectospec. We confirm 431 unique GCs belonging to NGC 4649, a factor of 3.5 larger than previous datasets and with a factor of 3 improvement in velocity precision. We confirm significant GC colour bimodality and find that the red GCs are more centrally concentrated, while the blue GCs are more spatially extended. We infer negative GC colour gradients in the innermost 20 kpc and flat gradients out to large radii. Rotation is detected along the galaxy major axis for all tracers: blue GCs, red GCs, galaxy stars and planetary nebulae. We compare the observed properties of NGC 4649 with galaxy formation models. We find that formation via a major merger between two gas-poor galaxies, followed by satellite accretion, can consistently reproduce the observations of NGC 4649 at different radii. We find no strong evidence to support an interaction between NGC 4649 and the neighbouring spiral galaxy NGC 4647. We identify interesting GC kinematic features in our data, such as counter-rotating subgroups and bumpy kinematic profiles, which encode more clues about the formation history of NGC 4649.
Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse structure of the progenitor, which in combination allow for a clear separation of exploding and non-exploding cases with only few exceptions (~1--2.5%) in our set of 621 investigated stellar models. One parameter is M4, defining the enclosed mass for a dimensionless entropy per nucleon of s = 4, and the other is mu4 = dm/dr|_{s=4}, being the mass-derivative at this location. The two parameters mu4 and M4*mu4 can be directly linked to the mass-infall rate, Mdot, of the collapsing star and the electron-type neutrino luminosity of the accreting proto-NS, L_nue ~ M_ns*Mdot, which play a crucial role in the "critical luminosity" concept for the theoretical description of neutrino-driven explosions as runaway phenomenon of the stalled accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino-driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a gray transport solver, while the uncertain neutrino emission of the 1.1 Msun proto-NS core is parametrized by an analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce the observables of Supernova 1987A for five different progenitor models.
Over the past 15 years, examples of exotic radio-quiet quasars with intrinsically weak or absent broad emission line regions (BELRs) have emerged from large-scale spectroscopic sky surveys. Here, we present spectroscopy of seven such weak emission line quasars (WLQs) at moderate redshifts (z=1.4-1.7) using the X-shooter spectrograph, which provides simultaneous optical and near-infrared spectroscopy covering the rest-frame ultraviolet through optical. These new observations effectively double the number of WLQs with spectroscopy in the optical rest-frame, and they allow us to compare the strengths of (weak) high-ionization emission lines (e.g., CIV) to low-ionization lines (e.g., MgII, Hb, Ha) in individual objects. We detect broad Hb and Ha emission in all objects, and these lines are generally toward the weaker end of the distribution expected for typical quasars (e.g., Hb has rest-frame equivalent widths ranging from 15-40 Ang.). However, these low-ionization lines are not exceptionally weak, as is the case for high-ionization lines in WLQs. The X-shooter spectra also display relatively strong optical FeII emission, Hb FWHM <4000 km/s, and significant CIV blueshifts (1000-5500 km/s) relative to the systemic redshift; two spectra also show elevated ultraviolet FeII emission, and an outflowing component to their (weak) MgII emission lines. These properties suggest that WLQs are exotic versions of "wind-dominated" quasars. Their BELRs either have unusual high-ionization components, or their BELRs are in an atypical photoionization state because of an unusually soft continuum.
In 2006 February, shortly after its launch, Swift began monitoring the center of the Milky Way with the onboard X-Ray Telescope using short 1-ks exposures performed every 1-4 days. Between 2006 and 2014, over 1200 observations have been obtained, amounting to ~1.2 Ms of exposure time. This has yielded a wealth of information about the long-term X-ray behavior of the supermassive black hole Sgr A*, and numerous transient X-ray binaries that are located within the 25'x25' region covered by the campaign. In this review we highlight the discoveries made during these first nine years, which includes 1) the detection of seven bright X-ray flares from Sgr A*, 2) the discovery of the magnetar SGR J1745-29, 3) the first systematic analysis of the outburst light curves and energetics of the peculiar class of very-faint X-ray binaries, 4) the discovery of three new transient X-ray sources, 5) exposing low-level accretion in otherwise bright X-ray binaries, and 6) the identification of a candidate X-ray binary/millisecond radio pulsar transitional object. We also reflect on future science to be done by continuing this Swift's legacy campaign of the Galactic center, which includes high-cadence monitoring of how the interaction between the gaseous object `G2' and Sgr A* plays out in the future.
Results from exoplanet surveys indicate that small planets (super-Earth size and below) are abundant in our Galaxy. However, little is known about their interiors and atmospheres. There is therefore a need to find small planets transiting bright stars, which would enable a detailed characterisation of this population of objects. We present the results of a search for the transit of the Earth-mass exoplanet Alpha Centauri Bb with the Hubble Space Telescope (HST). We observed Alpha Centauri B twice in 2013 and 2014 for a total of 40 hours. We achieve a precision of 115 ppm per 6-s exposure time in a highly-saturated regime, which is found to be consistent across HST orbits. We rule out the transiting nature of Alpha Centauri Bb with the orbital parameters published in the literature at 96.6% confidence. We find in our data a single transit-like event that could be associated to another Earth-size planet in the system, on a longer period orbit. Our program demonstrates the ability of HST to obtain consistent, high-precision photometry of saturated stars over 26 hours of continuous observations.
In this paper, we show that the supernova remnant (SNR) masses cumulated from core-collapse supernovae along the star formation history of two powerful z=3.8 radio galaxies 4C41.17 and TN J2007-1316 reach up to > 10^9 Msun, comparable with supermassive black hole (SMBH) masses measured from the SDSS sample at similar redshifts. The SNR mass is measured from the already exploded supernova mass after subtraction of ejecta at the galaxy age where the mass of still luminous stars fits at best the observed spectral energy distribution (SED), continuously extended to the optical-Spitzer-Herschel-submm domains, with the help of the galaxy evolution model P\'egase.3. For the recent and old stellar populations, SNR masses vary on 10^(9 to 10) Msun and the SNR-to-star mass ratio between 1 and 0.1 percent is comparable to the observed low-z SMBH-to-star mass ratio. For the template radio galaxy 4C41.17, SNR and stellar population masses estimated from large aperture (>4arcsec=30kpc) observations are compatible, within one mass order, with the total mass of multiple optical HST (~700pc) structures, associated with VLA radio emissions, both at 0.1 arcsec. Probing the SNR accretion by central black holes is a simple explanation for SMBH growth, requiring physics on star formation, stellar and galaxy dynamics with consequences on various processes (quenching, mergers, negative feedback) and a key to the relation bulge-SMBH.
Several scenarios have been proposed to explain the presence of multiple stellar populations in globular clusters. Many of them invoke multiple generations of stars to explain the observed chemical abundance anomalies, but it has also been suggested that self-enrichment could occur via accretion of ejecta from massive stars onto the circumstellar disc of low-mass pre-main sequence stars. These scenarios imply different initial conditions for the kinematics of the various stellar populations. Given some net angular momentum initially, models for which a second generation forms from gas that collects in a cooling flow into the core of the cluster predict an initially larger rotational amplitude for the polluted stars compared to the pristine stars. This is opposite to what is expected from the accretion model, where the polluted stars are the ones crossing the core and are on preferentially radial (low-angular momentum) orbits, such that their rotational amplitude is lower. Here we present the results of a suite of $N$-body simulations with initial conditions chosen to capture the distinct kinematic properties of these pollution scenarios. We show that initial differences in the kinematics of polluted and pristine stars can survive to the present epoch in the outer parts of a large fraction of Galactic globular clusters. The differential rotation of pristine and polluted stars is identified as a unique kinematic signature that could allow us to distinguish between various scenarios, while other kinematic imprints are generally very similar from one scenario to the other.
Halo Abundance Matching has been used to construct a one-parameter mapping between galaxies and dark matter haloes by assuming that halo mass and galaxy luminosity (or stellar mass) are monotonically related. While this approach has been reasonably successful, it is known that galaxies must be described by at least two parameters, as can be seen from the two-parameter Fundamental Plane on which massive early-type galaxies lie. In this paper, we derive a connection between initial dark matter density perturbations in the early universe and present-day virialized dark matter haloes by assuming simple spherical collapse combined with conservation of mass and energy. We find that $z = 0$ halo concentration, or alternatively the inner slope of the halo density profile $\alpha$, is monotonically and positively correlated with the collapse redshift of the halo. This is qualitatively similar to the findings of some previous works based on numerical simulations, with which we compare our results. We then describe how the halo mass and concentration (or inner slope $\alpha$) can be used as two halo parameters in combination with two parameters of early-type galaxies to create an improved abundance matching scheme.
Here we present a simple, parameter-free, non-perturbative algorithm that gives low-redshift cosmological particle realizations accurate to few-Megaparsec scales, called muscle (MUltiscale Spherical ColLapse Evolution). It has virtually the same cost as producing N-body-simulation initial conditions, since it works with the 'stretch' parameter {\psi}, the Lagrangian divergence of the displacement field. It promises to be useful in quickly producing mock catalogs, and to simplify computationally intensive reconstructions of galaxy surveys. muscle applies a spherical-collapse prescription on multiple Gaussian-smoothed scales. It achieves higher accuracy than perturbative schemes (Zel'dovich and 2LPT), and, by including the void-in-cloud process (voids in large-scale collapsing regions), solves problems with a single-scale spherical-collapse scheme. Additionally, we show the behavior of {\psi} for different morphologies (voids, walls, filaments, and haloes). A Python code to produce these realizations is available at this http URL
For the first time a proper comparison of the average depth of shower maximum ($X_{\rm max}$) published by the Pierre Auger and Telescope Array Observatories is presented. The $X_{\rm max}$ distributions measured by the Pierre Auger Observatory were fit using simulated events initiated by four primaries (proton, helium, nitrogen and iron). The primary abundances which best describe the Auger data were simulated through the Telescope Array (TA) Middle Drum (MD) fluorescence and surface detector array. The simulated events were analyzed by the TA Collaboration using the same procedure as applied to their data. The result is a simulated version of the Auger data as it would be observed by TA. This analysis allows a direct comparison of the evolution of $\langle X_{\rm max} \rangle$ with energy of both data sets. The $\langle X_{\rm max} \rangle$ measured by TA-MD is consistent with a preliminary simulation of the Auger data through the TA detector and the average difference between the two data sets was found to be $(2.9 \pm 2.7\;(\text{stat.}) \pm 18\;(\text{syst.}))~\text{g/cm}^2$.
The NEA 2001 SN263 is the target of the ASTER MISSION - First Brazilian Deep Space Mission. Araujo et al. (2012), characterized the stable regions around the components of the triple system for the planar and prograde cases. Knowing that the retrograde orbits are expected to be more stable, here we present a complementary study. We now considered particles orbiting the components of the system, in the internal and external regions, with relative inclinations between $90^{\circ}< I \leqslant180^{\circ}$, i.e., particles with retrograde orbits. Our goal is to characterize the stable regions of the system for retrograde orbits, and then detach a preferred region to place the space probe. For a space mission, the most interesting regions would be those that are unstable for the prograde cases, but stable for the retrograde cases. Such configuration provide a stable region to place the mission probe with a relative retrograde orbit, and, at the same time, guarantees a region free of debris since they are expected to have prograde orbits. We found that in fact the internal and external stable regions significantly increase when compared to the prograde case. For particles with $e=0$ and $I=180^{\circ}$, we found that nearly the whole region around Alpha and Beta remain stable. We then identified three internal regions and one external region that are very interesting to place the space probe. We present the stable regions found for the retrograde case and a discussion on those preferred regions. We also discuss the effects of resonances of the particles with Beta and Gamma, and the role of the Kozai mechanism in this scenario. These results help us understand and characterize the stability of the triple system 2001 SN263 when retrograde orbits are considered, and provide important parameters to the design of the ASTER mission.
We investigate the relationship between the rest-frame equivalent width (EW) of the C IV \lambda1549 broad-emission line, monochromatic luminosity at rest-frame 5100 A, and the Hbeta-based Eddington ratio in a sample of 99 ordinary quasars across the widest possible ranges of redshift (0 < z < 3.5) and bolometric luminosity (10^{44} <~ L <~ 10^{48} erg s^{-1}). We find that EW(C IV) is primarily anti-correlated with the Eddington ratio, a relation we refer to as a modified Baldwin effect (MBE), an extension of the result previously obtained for quasars at z < 0.5. Based on the MBE, weak emission line quasars (WLQs), typically showing EW(C IV) <~ 10 A, are expected to have extremely high Eddington ratios. By selecting all WLQs with archival Hbeta and C IV spectroscopic data, nine sources in total, we find that their Hbeta-based Eddington ratios are typical of ordinary quasars with similar redshifts and luminosities. Four of these WLQs can be accommodated by the MBE, but the other five deviate significantly from this relation, at the >~3 \sigma\ level, by exhibiting C IV lines much weaker than predicted from their Hbeta-based Eddington ratios. Assuming the supermassive black-hole masses in all quasars can be determined reliably using the single-epoch Hbeta-method, our results indicate that EW(C IV) cannot depend solely on the Eddington ratio. We briefly discuss a strategy for further investigation into the roles that basic physical properties play in controlling the relative strengths of broad-emission lines in quasars.
We present Giant Metrewave Radio Telescope (GMRT) 240 MHz observations of the nearby luminous FR I radio source 3C 270, in the group-central elliptical NGC 4261. Combining these data with reprocessed Very Large Array (VLA) 1.55 and 4.8 GHz observations, we produce spectral index maps that reveal a constant spectral index along the jets and a gradual steepening from the ends of the jets through the lobes towards the nucleus. A Jaffe & Perola (JP) model fitted to the integrated spectrum of the source gives an asymptotic low-frequency index of $\alpha_{inj}=0.53_{-0.02}^{+0.01}$, while JP models fitted to the observed spectral index trend along the lobes allow us to estimate radiative ages of $\sim29$ Myr and $\sim37$ Myr for the west and east lobes respectively. Our age estimates are a factor of two lower than the 75-Myr upper limit derived from X-ray data (O'Sullivan et al. 2011). We find unlikely the scenario of an early supersonic phase in which the lobe expanded into the ISM at approximately Mach 6 (3500 km s$^{-1}$), and suggest that either the source underwent multiple AGN outbursts with possible large changes in jet power, or possibly that the source age that we find is due to a backflow that transports young electrons from the jet tips through the lobes toward the nucleus relatively quickly. We calculate that in the lobes the energy ratio of non-radiating to radiating particles is $\sim4-24$ indicating significant gas entrainment. If the lobes are in pressure balance with their surroundings, the total energy required to heat the entrained material is $10^{58}$ erg, $\sim$40% of the total enthalpy of the lobes.
Power suppression of the cosmic microwave background on the largest observable scales could provide valuable clues about the particle physics underlying inflation. Here we consider the prospect of power suppression in the context of the multifield landscape. Based on the assumption that our observable universe emerges from a tunnelling event and that the relevant features originate purely from inflationary dynamics, we find that the power spectrum not only contains information on single-field dynamics, but also places strong con- straints on all scalar fields present in the theory. We find that the simplest single-field models giving rise to power suppression do not generalise to multifield models in a straightforward way, as the resulting superhorizon evolution of the curvature perturbation tends to erase any power suppression present at horizon crossing. On the other hand, multifield effects do present a means of generating power suppression which to our knowledge has so far not been considered. We propose a mechanism to illustrate this, which we dub flume inflation.
We present new 0.9-2.45 $\mu$m spectroscopy ($R \sim 1000$), and $Y$, $J$, $H$, $K_s$, $L^\prime$ photometry, obtained at Gemini North, of three low-mass brown dwarf companions on wide orbits around young stars of the Upper Scorpius OB association: HIP 78530 B, [PGZ2001] J161031.9-191305 B, and GSC 06214-00210 B. We use these data to assess the companions' spectral type, temperature, surface gravity and mass, as well as the ability of the BT-Settl and Drift-Phoenix atmosphere models to reproduce the spectral features of young substellar objects. For completeness, we also analyze the archival spectroscopy and photometry of the Upper Scorpius planetary mass companion 1RXS J160929.1-210524 b. Based on a comparison with model spectra we find that the companions, in the above order, have effective temperatures of 2700, 2500, 2300 and 1700 K. These temperatures are consistent with our inferred spectral types, respectively M7 $\beta$, M9 $\gamma$, M9 $\gamma$, and L4 $\gamma$. From bolometric luminosities estimated from atmosphere model spectra adjusted to our photometry, and using evolution models at 5-10 Myr, we estimate masses of 21-25, 28-70, 14-17 and 7-12 $M_{\rm Jup}$, respectively. J1610-1913 B appears significantly over-luminous for its inferred temperature, which explains its higher mass estimate. Synthetic spectra based on the BT-Settl and Drift-Phoenix atmosphere models generally offer a good fit to our observed spectra, although our analysis has highlighted a few problems. For example, the best fits in the individual near-infrared bands occur at different model temperatures. Also, temperature estimates based on a comparison of the broadband magnitudes and colors of the companions to synthetic magnitudes from the models are systematically lower than the temperature estimates based on a comparison with synthetic spectra.
A planned rapid submillimeter (submm) Gamma Ray Burst (GRBs) follow-up observations conducted using the Greenland Telescope (GLT) is presented. The GLT is a 12-m submm telescope to be located at the top of the Greenland ice sheet, where the high-altitude and dry weather porvides excellent conditions for observations at submm wavelengths. With its combination of wavelength window and rapid responding system, the GLT will explore new insights on GRBs. Summarizing the current achievements of submm GRB follow-ups, we identify the following three scientific goals regarding GRBs: (1) systematic detection of bright submm emissions originating from reverse shock (RS) in the early afterglow phase, (2) characterization of forward shock and RS emissions by capturing their peak flux and frequencies and performing continuous monitoring, and (3) detections of GRBs as a result of the explosion of first-generation stars result of GRBs at a high redshift through systematic rapid follow ups. The light curves and spectra calculated by available theoretical models clearly show that the GLT could play a crucial role in these studies.
Thanks to the recent very high-precision photometry of red giants from satellites such as Kepler, precise mass and radius values as well as accurate information of evolutionary stages are already established by asteroseismic approach for a large number of G-K giants. Based on the high-dispersion spectra of selected such 55 red giants in the Kepler field with precisely known seismic parameters (among which parallaxes are available for 9 stars), we checked the accuracy of the determination method of stellar parameters previously applied to many red giants by Takeda et al. (2008, PASJ, 60, 781), since it may be possible to discriminate their complex evolutionary status by using the surface gravity vs. mass diagram. We confirmed that our spectroscopic gravity and the seismic gravity satisfactorily agree with each other (to within ~0.1 dex) without any systematic difference. However, the mass values of He-burning red clump giants derived from stellar evolutionary tracks (~2-3 Msun) were found to be markedly larger by ~50% compared to the seismic values (~1-2 Msun) though such discrepancy is not seen for normal giants in the H-burning phase, which reflects the difficulty of mass determination from intricately overlapping tracks on the luminosity vs. effective temperature diagram. This consequence implies that the mass results of many red giants in the clump region determined by Takeda et al. (2008) are likely to be significantly overestimated. We also compare our spectroscopically established parameters with recent literature values, and further discuss the prospect of distinguishing the evolutionary status of red giants based on the conventional (i.e., non-seismic) approach.
Evolution in the measured rest frame ultraviolet spectral slope and ultraviolet to optical flux ratios indicate a rapid evolution in the dust obscuration of galaxies during the first 3 billion years of cosmic time (z>4). This evolution implies a change in the average interstellar medium properties, but the measurements are systematically uncertain due to untested assumptions, and the inability to measure heavily obscured regions of the galaxies. Previous attempts to directly measure the interstellar medium in normal galaxies at these redshifts have failed for a number of reasons with one notable exception. Here we report measurements of the [CII] gas and dust emission in 9 typical (~1-4L*) star-forming galaxies ~1 billon years after the big bang (z~5-6). We find these galaxies have >12x less thermal emission compared with similar systems ~2 billion years later, and enhanced [CII] emission relative to the far-infrared continuum, confirming a strong evolution in the interstellar medium properties in the early universe. The gas is distributed over scales of 1-8 kpc, and shows diverse dynamics within the sample. These results are consistent with early galaxies having significantly less dust than typical galaxies seen at z<3 and being comparable to local low-metallicity systems.
We report a detailed modeling of a mass profile of a $z = 0.2999$ massive elliptical galaxy using 30 milli-arcsecond resolution 1-mm Atacama Large Millimeter/submillimeter Array (ALMA) images of the galaxy-galaxy lensing system SDP.81. The detailed morphology of the lensed multiple images of the $z = 3.042$ infrared-luminous galaxy, which is found to consist of tens of $\lesssim 100$-pc-sized star-forming clumps embedded in a $\sim 2$ kpc disk, are well reproduced by a lensing galaxy modeled by an isothermal ellipsoid with a 400 pc core. The core radius is consistent with that of the visible stellar light, and the mass-to-light ratio of $\sim 2\,M_{\odot}\,L_{\odot}^{-1}$ is comparable to the locally measured value, suggesting that the inner 1 kpc region is dominated by luminous matter. The position of the predicted mass centroid is consistent to within $\simeq 30$ mas with that of a non-thermal source detected with ALMA, which likely traces an active galactic nucleus of the foreground elliptical galaxy. A point source mass of $> 3 \times 10^8\,M_{\odot}$ mimicking a supermassive black hole is required to explain the non-detection of a central image of the background galaxy, although the black hole mass degenerates with the core radius of the elliptical galaxy. The required mass is consistent with that predicted from the well-known correlation between black hole mass and host velocity dispersion. Our analysis demonstrates the power of ALMA imaging of strong gravitational lensing for studying the innermost mass profiles and the central supermassive black hole of distant elliptical galaxies.
We present observations of pseudobulges in S0 and spiral galaxies using imaging data taken with the Spitzer Infra-Red Array Camera. We have used 2-d bulge-disk-bar decomposition to determine structural parameters of 185 S0 galaxies and 31 nearby spiral galaxies. Using the Sersic index and the position on the Kormendy diagram to classify their bulges as either classical or pseudo, we find that 25 S0s (14%) and 24 spirals (77%) host pseudoblges. The fraction of pseudobulges we find in spiral galaxies is consistent with previous results obtained with optical data and show that the evolution of a large fraction of spirals is governed by secular processes rather than by major mergers. We find that the bulge effective radius is correlated with the disk scale length for pseudobulges of S0s and spirals, as expected for secular formation of bulges from disk instabilities, though the disks in S0s are significantly smaller than those in spirals. We show that early-type pseudobulge hosting spirals can transform to pseudobulge hosting S0s by simple gas stripping. However, simple gas stripping mechanism is not sufficient to transform the late-type pseudobulge hosting spirals into pseudobulge hosting S0s.
We report detections of new exoplanets from a radial velocity (RV) survey of metal-rich FGK stars by using three telescopes. By optimizing our RV analysis method to long time-baseline observations, we have succeeded in detecting five new Jovian-planets around three metal-rich stars HD 1605, HD 1666, and HD 67087 with the masses of $1.3 M_{\odot}$, $1.5 M_{\odot}$, and $1.4 M_{\odot}$, respectively. A K1 subgiant star HD 1605 hosts two planetary companions with the minimum masses of $ M_p \sin i = 0.96 M_{\mathrm{JUP}}$ and $3.5 M_{\mathrm{JUP}}$ in circular orbits with the planets' periods $P = 577.9$ days and $2111$ days, respectively. HD 1605 shows a significant linear trend in RVs. Such a system consisting of Jovian planets in circular orbits has rarely been found and thus HD 1605 should be an important example of a multi-planetary system that is likely unperturbed by planet-planet interactions. HD 1666 is a F7 main sequence star which hosts an eccentric and massive planet of $ M_p \sin i = 6.4 M_{\mathrm{JUP}}$ in the orbit with $a_{\rm p} = 0.94$ AU and an eccentricity $e=0.63$. Such an eccentric and massive planet can be explained as a result of planet-planet interactions among Jovian planets. While we have found the large residuals of $\mathrm{rms} = 35.6\ \mathrm{m\ s^{-1}}$, the periodogram analysis does not support any additional periodicities. Finally, HD 67087 hosts two planets of $ M_p \sin i = 3.1 M_{\mathrm{JUP}}$ and $4.9 M_{\mathrm{JUP}}$ in orbits with $P=352.2$ days and $2374$ days, and $e=0.17$ and $0.76$, respectively. Although the current RVs do not lead to accurate determinations of its orbit and mass, HD 67087 c can be one of the most eccentric planets ever discovered in multiple systems.
The binary system eta Carinae has completed its first 5.54y orbit since the beginning of science operation of the Fermi Large Area Telescope (LAT). We are now able to investigate the high-energy gamma-ray source at the position of eta Carinae over its full orbital period. By this, we can address and confirm earlier predictions for temporal and spectral variability. Newer versions of the LAT datasets, instrument response functions and background models allow for a more accurate analysis. Therefore it is important to re-evaluate the previously analyzed time period along with the new data to further constrain location, spectral shape, and flux time history of the gamma-ray source. We confirm earlier predictions of increasing flux values above 10 GeV toward the next periastron passage. For the most recent part of the data sample, flux values as high as those before the first periastron passage in 2008 are recorded. A comparison of spectral energy distributions around periastron and apastron passages reveals strong variation in the high-energy band. This is due to a second spectral component that is present only around periastron. Improved spatial consistency with the gamma-ray source at the position of eta Carinae along with the confirmation of temporal variability above 10 GeV in conjunction with the orbital period strengthens the argument for unambiguous source identification. Spectral variability provides additional constraints for future modeling of the particle acceleration and gamma-ray emission in colliding-wind binary systems.
We investigate the power spectra of outflow-driven turbulence through high-resolution three-dimensional isothermal numerical simulations where the turbulence is driven locally in real-space by a simple spherical outflow model. The resulting turbulent flow saturates at an average Mach number of ~2.5 and is analysed through density and velocity power spectra, including an investigation of the evolution of the solenoidal and compressional components. We obtain a shallow density power spectrum with a slope of ~-1.2 attributed to the presence of a network of localised dense filamentary structures formed by strong shock interactions. The total velocity power spectrum slope is found to be ~-2.0, representative of Burgers shock dominated turbulence model. The density weighted velocity power spectrum slope is measured as ~-1.6, slightly less than the expected Kolmogorov scaling value (slope of -5/3) found in previous works. The discrepancy may be caused by the nature of our real space driving model and we suggest there is no universal scaling law for supersonic compressible turbulence. We find that on average, solenoidal modes slightly dominate in our turbulence model as the interaction between strong curved compressible shocks generates solenoidal modes, and compressible modes decay faster.
We present the size, shape and flux densities at millimeter continuum wavelengths, based on ALMA science verification observations in Band 3 (~94.6 GHz) and Band 6 (~228.7 GHz), from the binary Mira A (o Ceti) and Mira B. The Mira AB system has been observed with ALMA at a spatial resolution of down to ~25 mas. The extended atmosphere of Mira A and the wind around Mira B sources are resolved and we derive the size of Mira A and of the ionized region around Mira B. The spectral indices within Band 3 (between 89-100 GHz) and between Band 3 and Band 6 are also derived. The spectral index of Mira A is found to change from 1.71+-0.05 within Band 3 to 1.54+-0.04 between Band 3 and 6. The spectral index of Mira B is 1.3+-0.2 in Band 3, in good agreement with measurements at longer wavelengths. However it rises to 1.72+-0.11 between the bands. For the first time the extended atmosphere of a star is resolved at these frequencies and for Mira A the diameter is ~3.7x2.9 AU in Band 3 (with brightness temperature Tb~6000 K) and ~4.0x3.6 AU in Band 6 (Tb~2500 K). Additionally, a bright hotspot of ~0.4 AU and with Tb~10000 K is found on the stellar disc of Mira A. The size of the ionized region around the accretion disk of Mira B is found to be ~2.4 AU. The emission around Mira B is consistent with that from a partially ionized wind of gravitationally bound material from Mira A close to the accretion disk of Mira B. The Mira A atmosphere does not fully match predictions, with brightness temperatures in Band 3 significantly higher than expected, potentially due to shock heating. The hotspot is likely due to magnetic activity and could be related to the previously observed X-ray flare of Mira A.
Galactic stars belonging to the Of?p category are all strongly magnetic objects exhibiting rotationally modulated spectral and photometric changes on timescales of weeks to years. Five candidate Of?p stars in the Magellanic Clouds have been discovered, notably in the context of ongoing surveys of their massive star populations. Here we describe an investigation of their photometric behaviour, revealing significant variability in all studied objects on timescales of one week to more than four years, including clearly periodic variations for three of them. Their spectral characteristics along with these photometric changes provide further support for the hypothesis that these are strongly magnetized O stars, analogous to the Of?p stars in the Galaxy.
The role of disk instabilities, such as bars and spiral arms, and the associated resonances, in growing bulges in the inner regions of disk galaxies have long been studied in the low-redshift nearby Universe. There it has long been probed observationally, in particular through peanut-shaped bulges. This secular growth of bulges in modern disk galaxies is driven by weak, non-axisymmetric instabilities: it mostly produces pseudo-bulges at slow rates and with long star-formation timescales. Disk instabilities at high redshift (z>1) in moderate-mass to massive galaxies (10^10 to a few 10^11 Msun of stars) are very different from those found in modern spiral galaxies. High-redshift disks are globally unstable and fragment into giant clumps containing 10^8-10^9 Msun of gas and stars each, which results in highly irregular galaxy morphologies. The clumps and other features associated to the violent instability drive disk evolution and bulge growth through various mechanisms, on short timescales. The giant clumps can migrate inward and coalesce into the bulge in a few 10^8 yr. The instability in the very turbulent media drives intense gas inflows toward the bulge and nuclear region. Thick disks and supermassive black holes can grow concurrently as a result of the violent instability. This chapter reviews the properties of high-redshift disk instabilities, the evolution of giant clumps and other features associated to the instability, and the resulting growth of bulges and associated sub-galactic components.
The Maunder minimum (MM) of greatly reduced solar activity took place in
1645-1715, but the exact level of sunspot activity is uncertain as based, to a
large extent, on historical generic statements of the absence of spots on the
Sun. Here we aim, using a conservative approach, to assess the level and length
of solar cycle during the Maunder minimum, on the basis of direct historical
records by astronomers of that time. A database of the active and inactive days
(days with and without recorded sunspots on the solar disc respectively) is
constructed for three models of different levels of conservatism (loose ML,
optimum MO and strict MS models) regarding generic no-spot records. We have
used the active day fraction to estimate the group sunspot number during the
MM. A clear cyclic variability is found throughout the MM with peaks at around
1655--1657, 1675, 1684 and 1705, and possibly 1666, with the active day
fraction not exceeding 0.2, 0.3 or 0.4 during the core MM, for the three
models. Estimated sunspot numbers are found very low in accordance with a grand
minimum of solar activity.
We have found, for the core MM (1650-1700), that: (1) A large fraction of
no-spot records, corresponding to the solar meridian observations, may be
unreliable in the conventional database. (2) The active day fraction remained
low (below 0.3-0.4) throughout the MM, indicating the low level of sunspot
activity. (3) The solar cycle appears clearly during the core MM. (4) The
length of the solar cycle during the core MM appears $9\pm 1$ years, but there
is an uncertainty in that. (5) The magnitude of the sunspot cycle during MM is
assessed to be below 5-10 in sunspot numbers;
A hypothesis of the high solar cycles during the MM is not confirmed.
Inversion codes are numerical tools used for the inference of physical properties from the observations. Despite their success, the quality of current spectropolarimetric observations and those expected in the near future presents a challenge to current inversion codes. The pixel-by-pixel strategy of inverting spectropolarimetric data that we currently utilize needs to be surpassed and improved. The inverted physical parameters have to take into account the spatial correlation that is present in the data and that contains valuable physical information. We utilize the concept of sparsity or compressibility to develop an new generation of inversion codes for the Stokes parameters. The inversion code uses numerical optimization techniques based on the idea of proximal algorithms to impose sparsity. In so doing, we allow for the first time to exploit the presence of spatial correlation on the maps of physical parameters. Sparsity also regularizes the solution by reducing the number of unknowns. We compare the results of the new inversion code with pixel-by-pixel inversions, demonstrating the increase in robustness of the solution. We also show how the method can easily compensate for the effect of the telescope point spread function, producing solutions with an enhanced contrast.
In 2016, the NASA Juno spacecraft will initiate its one-year mission around Jupiter and become the first probe to explore the polar regions of Jupiter. The HST UV instruments (STIS and ACS) can greatly contribute to the success of the Juno mission by providing key complementary views of Jupiter's UV aurora from Earth orbit. Juno carries an ultraviolet Spectrograph (UVS) and an infrared spectral mapper (JIRAM) that will obtain high-resolution spectral images providing the auroral counterpart to Juno's in situ particles and fields measurements with the plasma JADE and JEDI particle detectors. The Juno mission will be the first opportunity to measure simultaneously the energetic particles at high latitude and the auroral emissions they produce. Following programmatic and technical limitations, the amount of UVS data transmitted to Earth will be severely restricted. Therefore, it is of extreme importance that HST captures as much additional information as possible on Jupiter's UV aurora during the one-year life of the Juno mission. This white paper is a plea for a "Juno initiative" that will ensure that a sufficient number of orbits is allocated to this unique solar system mission.
A properly calibrated longitudinal magnetograph is an instrument that measures circular polarization and gives an estimation of the magnetic flux density in each observed resolution element. This usually constitutes a lower bound of the field strength in the resolution element, given that it can be made arbitrarily large as long as it occupies a proportionally smaller area of the resolution element and/or becomes more transversal to the observer and still produce the same magnetic signal. Yet, we know that arbitrarily stronger fields are less likely --hG fields are more probable than kG fields, with fields above several kG virtually absent-- and we may even have partial information about its angular distribution. Based on a set of sensible considerations, we derive simple formulae based on a Bayesian analysis to give an improved estimation of the magnetic field strength for magnetographs.
We report energy levels, radiative rates (A-values) and lifetimes for the astrophysically-important Be-like ion C III. For the calculations, 166 levels belonging to the $n \le$ 5 configurations are considered and the {\sc grasp} (General-purpose Relativistic Atomic Structure Package) is adopted. Einstein A-coefficients are provided for all E1, E2, M1 and M2 transitions, while lifetimes are compared with available measurements as well as theoretical results, and no large discrepancies noted. Our energy levels are assessed to be accurate to better than 1\% for a majority of levels, and A-values to better than 20\% for most transitions. Collision strengths are also calculated, for which the Dirac Atomic R-matrix Code ({\sc darc}) is used. A wide energy range, up to 21 Ryd, is considered and resonances resolved in a fine energy mesh in the thresholds region. The collision strengths are subsequently averaged over a Maxwellian velocity distribution to determine effective collision strengths up to a temperature of 8.0$\times$10$^5$ K, sufficient for most astrophysical applications. Our data are compared with the recent $R$-matrix calculations of Fern{\'a}ndez-Menchero, Del Zanna \& Badnell [A\&A 566 (2014) A104], and significant differences (up to over an order of magnitude) are noted for several transitions over the complete temperature range of the results.
Collisions between galaxy clusters provide a test of the non-gravitational forces acting on dark matter. Dark matter's lack of deceleration in the `bullet cluster collision' constrained its self-interaction cross-section \sigma_DM/m < 1.25cm2/g (68% confidence limit) for long-ranged forces. Using the Chandra and Hubble Space Telescopes we have now observed 72 collisions, including both `major' and `minor' mergers. Combining these measurements statistically, we detect the existence of dark mass at 7.6\sigma significance. The position of the dark mass has remained closely aligned within 5.8+/-8.2 kpc of associated stars: implying a self-interaction cross-section \sigma_DM/m < 0.47 cm2/g (95% CL) and disfavoring some proposed extensions to the standard model.
We show that the number of observed voids in galaxy redshift surveys is a sensitive function of the equation of state of dark energy. Using the Fisher matrix formalism we find the error ellipses in the $w_0-w_a$ plane when the equation of state of dark energy is assumed to be of the form $w_{CPL}(z)=w_0 +w_a z/(1+z)$. We forecast the number of voids to be observed with the ESA Euclid satellite and the NASA WFIRST mission, taking into account updated details of the surveys to reach accurate estimates of their power. The theoretical model for the forecast of the number of voids is based on matches between abundances in simulations and the analytical prediction. To take into account the uncertainties within the model, we marginalize over its free parameters when calculating the Fisher matrices. The addition of the void abundance constraints to the data from Planck, HST and supernova survey data noticeably tighten the $w_0-w_a$ parameter space. We thus quantify the improvement in the constraints due to the use of voids and demonstrate that the void abundance is a sensitive new probe for the dark energy equation of state.
We present Chandra and XMM-Newton observations of PLCK G036.7+14.9 from the Chandra-Planck Legacy Program. The high resolution X-ray observations reveal two close subclusters, G036N and G036S, which were not resolved by previous ROSAT, optical, or recent Planck observations. We perform detailed imaging and spectral analyses and use a simplified model to study the kinematics of this system. The basic picture is that PLCK G036.7+14.9 is undergoing a major merger (mass ratio close to unity) between the two massive subclusters, with the merger largely along the line-of-sight and probably at an early stage. G036N hosts a small, moderate cool-core, while G036S has at most a very weak cool-core in the central 40 kpc region. The difference in core cooling times is unlikely to be caused by the ongoing merger disrupting a pre-existing cool-core in G036S. G036N also hosts an unresolved radio source in the center, which may be heating the gas if the radio source is extended. The Planck derived mass is higher than the X-ray measured mass of either subcluster, but is lower than the X-ray measured mass of the whole cluster, due to the fact that Planck does not resolve PLCK G036.7+14.9 into subclusters and interprets it as a single cluster. This mass discrepancy could induce significant bias to the mass function if such previously unresolved systems are common in the Planck cluster sample. High resolution X-ray observations are necessary to identify the fraction of such systems and correct such a bias for the purpose of precision cosmological studies.
The low-background, VUV-sensitive 3-inch diameter photomultiplier tube R11410 has been developed by Hamamatsu for dark matter direct detection experiments using liquid xenon as the target material. We present the results from the joint effort between the XENON collaboration and the Hamamatsu company to produce a highly radio-pure photosensor (version R11410-21) for the XENON1T dark matter experiment. After introducing the photosensor and its components, we show the methods and results of the radioactive contamination measurements of the individual materials employed in the photomultiplier production. We then discuss the adopted strategies to reduce the radioactivity of the various PMT versions. Finally, we detail the results from screening 216 tubes with ultra-low background germanium detectors, as well as their implications for the expected electronic and nuclear recoil background of the XENON1T experiment.
We consider a hypothetical observatory of ultra-high energy cosmic rays consisting of two surface detector arrays that measure independently electromagnetic and muon signals induced by air showers. Using the constant intensity cut method, sets of events ordered according to each of both signal sizes are compared giving the number of matched events. Based on its dependence on the zenith angle, a parameter sensitive to the dispersion of the distribution of the logarithmic mass of cosmic rays is introduced. The results obtained using two post-LHC models of hadronic interactions are very similar and indicate a weak dependence on details of these interactions.
We present AMADA, an interactive web application to analyse multidimensional datasets. The user uploads a simple ASCII file and AMADA performs a number of exploratory analysis together with contemporary visualizations diagnostics. The package performs a hierarchical clustering in the parameter space, and the user can choose among linear, monotonic or non-linear correlation analysis. AMADA provides a number of clustering visualization diagnostics such as heatmaps, dendrograms, chord diagrams, and graphs. In addition, AMADA has the option to run a standard or robust principal components analysis, displaying the results as polar bar plots. The code is written in R and the web interface was created using the Shiny framework. AMADA is freely available at this http URL
This paper summarizes the limits on propagation of ultra high energy particles in the Universe, set up by their interactions with cosmic background of photons and neutrinos. By taking into account cosmic evolution of these backgrounds and considering appropriate interactions we derive the mean free path for ultra high energy photons, protons and neutrinos. For photons the relevant processes are the Breit-Wheeler process as well as the double pair production process. For protons the relevant reactions are the photopion production and the Bethe-Heitler process. We discuss the interplay between the energy loss length and mean free path for the Bethe-Heitler process. Neutrino opacity is determined by its scattering off the cosmic background neutrino. We compute for the first time the high energy neutrino horizon as a function of its energy.
Radio observations are an effective tool to discover particle acceleration regions in colliding-wind binaries, through detection of synchrotron radiation; these regions are natural laboratories for the study of relativistic particles. Wind-collision region (WCR) models can reproduce the radio continuum spectra of massive binaries that contain both thermal and non-thermal radio emission; however, key constraints for models come from high-resolution imaging. Only five WCRs have been resolved to date at radio frequencies at milliarcsec (mas) angular scales. The source HD 93129A, prototype of the very few known O2 I stars, is a promising target for study: recently, a second massive, early-type star about 50 mas away was discovered, and a non-thermal radio source detected in the region. Preliminary long-baseline array data suggest that a significant fraction of the radio emission from the system comes from a putative WCR. We sought evidence that HD 93129A is a massive binary system with colliding stellar winds that produce non-thermal radiation, through spatially resolved images of the radio emitting regions. We completed observations with the Australian Long Baseline Array (LBA) to resolve the system at mas angular resolutions and reduced archival Australia Telescope Compact Array (ATCA) data to derive the total radio emission. We also compiled optical astrometric data of the system in a homogeneous way. We reduced historical Hubble Space Telescope data and obtained absolute and relative astrometry with milliarcsec accuracy. The astrometric analysis leads us to conclude that the two stars in HD 93129A form a gravitationally bound system. The LBA data reveal an extended arc-shaped non-thermal source between the two stars, indicative of a WCR. The wind momentum-rate ratio of the two stellar winds is estimated. The ATCA data show a point source with a change in flux level ...
It was recently demonstrated that contact binaries occur in globular clusters (GCs) only immediately below turn-off point and in the region of blue straggler stars (BSs). In addition, observations indicate that at least a significant fraction of BSs in these clusters was formed by the binary mass-transfer mechanism. The aim of our present investigation is to obtain and analyze a set of evolutionary models of cool, close detached binaries with a low metal abundance, which are characteristic of GC. We computed the evolution of 975 models of initially detached, cool close binaries with different initial parameters. The models include mass exchange between components as well as mass and angular momentum loss due to the magnetized winds for very low-metallicity binaries with Z = 0.001. The models are interpreted in the context of existing data on contact binary and blue straggler members of GCs. The model parameters agree well with the observed positions of the GC contact binaries in the Hertzsprung-Russell diagram. Contact binaries in the lower part of the cluster main sequence are absent because there are no binaries with initial orbital periods shorter than 1.5 d. Contact binaries end their evolution as mergers that appear in the BS region. Binary-formed BSs populate the whole observed BS region in a GC, but a gap is visible between low-mass mergers that are concentrated along the zero-age main sequence and binary BSs occupying the red part of the BS region. Very few binary mergers are expected to rotate rapidly and/or possess chemical peculiarities resulting from the exposure of the layers processed by CNO nuclear reactions. All other binary mergers are indistinguishable from the collisionally formed mergers. The results show that binary-formed BSs may constitute at least a substantial fraction of all BSs in a GC.
GRB 060614 was a unique burst straddling both long and short duration gamma-ray bursts and its physical origin is still debated. Here we re-examine the afterglow data of GRB 060614 and find a significant F814W-band excess at $t\sim 13.6$ day after the burst observed by the {\it Hubble Space Telescope (HST)}. This corresponds to an extremely-low luminosity $\sim 3\times 10^{40}~{\rm erg~s^{-1}}$. The excess component has a very red spectrum and a rapid decline, both unexpected within the present theoretical framework of afterglow. We examine two possible sources of this signal$-$a very weak supernova and a Li-Paczynski Macronova/kilonova. We find that the observed signal is incompatible with a weak supernova. However, it is compatible with the ejection of $\sim 0.03-0.1~M_\odot$ of $r-$process material. If this interpretation is correct GRB 060614 arose from a compact binary (most likely a black hole$-$neutron star but also possibly a double neutron star) merger.
In radial velocity observations, a pair of extrasolar planets near a 2:1 orbital resonance can be misinterpreted as a single eccentric planet, if data are sparse and measurement precision insufficient to distinguish between these models. We determine the fraction of alleged single-planet RV detected systems for which a 2:1 resonant pair of planets is also a viable model and address the question of how the models can be disentangled. By simulation we quantified the mismatch arising from applying the wrong model. Model alternatives are illustrated using the supposed single-planet system HD 27894 for which we also study the dynamical stability of near-2:1 resonant solutions. From the data scatter around the fitted single-planet Keplerians, we find that for $74\% $ of the $254$ putative single-planet systems, a 2:1 resonant pair cannot be excluded as a viable model, since the error due to the wrong model is smaller than the scatter. For $187$ stars $\chi ^2$-probabilities can be used to reject the Keplerian models with a confidence of $95\% $ for $54\% $ of the stars and with $99.9\% $ for $39\% $ of the stars. For HD 27894 a considerable fit improvement is obtained when adding a low-mass planet near half the orbital period of the known Jovian planet. Dynamical analysis demonstrates that this system is stable when both planets are initially placed on circular orbits. For fully Keplerian orbits a stable system is only obtained if the eccentricity of the inner planet is constrained to $<0.3$. A large part of the allegedly RV detected single-planet systems should be scrutinized in order to determine the fraction of systems containing near-2:1 resonant pairs of planets. Knowing the abundance of such systems will allow us to revise the eccentricity distribution for extrasolar planets and provide direct constraints for planetary system formation.
In view of new constraints put forth by recent observations and measurements in the realm of astrophysics and nuclear physics, we update the non-linear realization of the sigma model as to reflect such constraints. By doing this, we obtain new equations of state that may be used to describe neutron stars. Such equations of state are obtained by investigating different ways by which the vector mesons self-interact. Furthermore, we also investigate the role played by the delta mesons in the model. As a result, we are able to develop equations of state that are in better agreement with data, such as nuclear compressibility and slope of the symmetry energy at saturation, star masses, radii, and cooling profiles.
A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.
The CRESST experiment (Cryogenic Rare Event Search with Superconducting
Thermometers) searches for dark matter via the phonon and light signal of
elastic scattering processes in scintillating crystals. The discrimination
between a possible dark matter signal and background requires good energy
resolution of the light detector, therefore a high light yield is important.
In this article, we present a method for understanding the light yield
measured with entire detector modules in terms of the efficiencies of light
production and detection. Based on data taken during a dark matter search
phase, it considers the entire process of conversion of deposited energy into
scintillation light as well as transport and collection of the light that occur
in a detector module. We can confirm the results by using a cross-check method
with different systematic uncertainties.
We found that with the detectors operated in CRESST-II phase 1, about 20% of
the produced scintillation light is detected. A part of the light loss is
likely caused by light absorption creating meta-stable excitations in the
scintillating crystals. We also found that, consistent with the relatively low
detection efficiency, an additional light detector increases the amount of
detected light within an otherwise unmodified detector module.
This is a brief rebuke to arXiv:1502.03821, which falls in the category "not even wrong."
It has recently been shown that second-order corrections to the background distance-redshift relation can build up significantly at large redshifts, due to an aggregation of gravitational lensing events. This shifts the expectation value of the distance to the CMB by 1%. In this paper we show that this shift is already properly accounted for in standard CMB analyses. We clarify the role that the distance to the CMB plays in the presence of second-order lensing corrections.
The unique capabilities of Swift that make it ideal for discovery and follow-up of Gamma-Ray bursts also makes it the idea mission for discovery and monitoring of X-ray Transients in the Milky Way and the Large and Small Magellanic Clouds. The Burst Alert Telescope allows for detection of new transient outbursts, the automated follow-up capabilities of Swift allow for rapid observation and localization of the new transient in X-rays and optical/UV bands, and Swift's rapid slewing capabilities allows for low-overhead short observations to be obtained, opening up the possibility of regular, sensitive, long term monitoring of transient outbursts that are not possible with other currently operational X-ray missions. In this paper I describe the methods of discovery of X-ray transients utilizing Swift's BAT and also collaboration with the MAXI telescope. I also detail two examples of X-ray transient science enabled by Swift: Swift discovery and monitoring observations of MAXI J1659-152, a Black Hole candidate Low Mass X-ray Binary in the Galactic Halo, which has the shortest known orbital period of any such system; and Swift monitoring of IGR J00569-7226, an edge on Be/X-ray binary that displayed a outburst in 2013 and 2014, and which monitoring by Swift allowed for detection of dips, eclipses and the determination of the orbital parameters, utilizing a measurement of doppler shifts in the pulsar period.
Bimetric theory describes gravitational interactions in the presence of an extra spin-2 field. Previous work has suggested that its cosmological solutions are generically plagued by instabilities. We show that by taking the Planck mass for the second metric, $M_f$, to be small, these instabilities can be pushed back to unobservably early times. In this limit, the theory approaches general relativity with an effective cosmological constant which is, remarkably, determined by the spin-2 interaction scale. This provides a late-time expansion history which is extremely close to $\Lambda$CDM, but with a technically-natural value for the cosmological constant. We find $M_f$ should be no larger than the electroweak scale in order for cosmological perturbations to be stable by big-bang nucleosynthesis.
We revisit the modeling of the properties of the black-hole remnant resulting the merger of a black-hole binary as a function of the parameters of the binary. We provide a set of empirical formulas for the final mass, spin and recoil velocity of the final black hole as a function of the mass ratio and individual spins of the progenitor. In order to determine the fitting coefficients for these formulas, we perform a set of 126 new numerical evolutions of precessing, unequal-mass black-hole binaries, and fit to the resulting remnant mass, spin, and recoil. In order to reduce the complexity of the analysis, we chose configurations that have one of the black holes spinning, with dimensionless spin alpha=0.8, at different angles with respect to the orbital angular momentum, and the other non-spinning. In addition to evolving families of binaries with different spin-inclination angles, we also evolved binaries with mass ratios as small as q=1/6. We use the resulting empirical formulas to predict the probabilities of black hole mergers leading to a given recoil velocity, total radiated gravitational energy, and final black hole spin.
In this paper we present a model for accelerated expansion of the universe, both during inflation and the present stage of the expansion, from four dimensional $\mathcal{N}=1$ supergravity. We evaluate the tensor-to-scalar ratio ($r\approx 0.00034$), the scalar spectral index ($n_s\approx 0.970$) and the running spetral index ($dn_s/dk\approx -6\times10^{-5}$), and we notice that these parameters are in agreement with Planck+WP+lensing data and with BICEP2/Keck and Planck joint analysis, at $95\%$ CL. The number of e-folds is $50$ or higher. The reheating period has an associated temperature $T_R\sim10^{12}$ Gev, which agrees with the one required by thermal leptogenesis. Regarding the scalar field as dark energy, the autonomous system for the scalar field in the presence of a barotropic fluid provides a stable fixed point that leads to a late-time accelerated expansion of the universe, with an equation of state that mimics the cosmological constant ($w_\Phi\approx -0.997$).
Scalar-tensor (ST) gravity is considered in the case where the scalar is an external field. We show that General Relativity (GR) and standard ST gravity are particular cases of the External Scalar-Tensor (EST) gravity. It is shown with a particular cosmological example that it is possible to join a part of a GR solution to a part of an ST one such that the complete solution neither belongs to GR nor to ST, but fully satisfies the EST field equations. We argue that external fields may effectively work as a type of screening mechanism for ST theories.
The technology of Silicon Drift Detectors (SDDs) has been selected for the two instruments aboard the Large Observatory For X-ray Timing (LOFT) space mission. LOFT underwent a three year long assessment phase as candidate for the M3 launch opportunity within the "Cosmic Vision 2015 -- 2025" long-term science plan of the European Space Agency. During the LOFT assessment phase, we studied the displacement damage produced in the SDDs by the protons trapped in the Earth's magnetosphere. In a previous paper we discussed the effects of the Non Ionising Energy Losses from protons on the SDD leakage current. In this paper we report the measurement of the variation of Charge Collection Efficiency produced by displacement damage caused by protons and the comparison with the expected damage in orbit.
A study of binary systems composed of two point particles with different masses in the linear regime of the characteristic formulation of general relativity is provided. The boundary conditions at the world tubes generated by the particle's orbits are explored, when the metric variables are decomposed in spin-weighted spherical harmonics. The power lost by the emission of gravitational waves is computed using the News Bondi's functions, and the contribution to the gravitational radiation of several multipole terms is shown.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)