Dark matter numerical simulations and the N-body method are essential for understanding how structure forms and evolves in the Universe. However, the discrete nature of N-body simulations can affect its accuracy when modelling collisionless systems. We introduce a new approach to simulate the gravitational evolution of cold collisionless fluids by solving the Vlasov-Poisson equations in terms of adaptively refineable "Lagrangian phase space elements". These geometrical elements are piecewise smooth maps between three-dimensional Lagrangian space and six-dimensional Eulerian phase space and approximate the continuum structure of the distribution function. They allow for dynamical adaptive splitting to follow the evolution even in regions of very strong mixing. We discuss various test problems which demonstrate the correctness and performance of our method. We show that it has several advantages compared to standard N-body algorithms by i) explicitly tracking the fine-grained distribution function, ii) naturally representing caustics, iii) providing an arbitrarily regular density field defined everywhere in space, iv) giving a smooth and regular gravitational potential field, thus eliminating the need for any type of ad-hoc force softening. Finally, we illustrate the feasibility of using our method for cosmological studies by simulating structure formation in a warm dark matter cosmology. We show that spurious collisionality and discreteness noise of N-body methods are both strongly suppressed, which eliminates artificial fragmentation of filaments while providing access to the full deterministic evolution of the fluid in phase space. Therefore, we argue that our new approach improves on the N-body method when simulating self-gravitating cold and collisionless fluids, and is the first method that allows to explicitly follow the fine-grained evolution in six-dimensional phase space.
The radial density profiles of stellar galaxy discs can be well approximated as an exponential. Compared to this canonical form, however, the profiles in the majority of disc galaxies show downward or upward breaks at large radii. Currently, there is no coherent explanation in a galaxy formation context of the radial profile per se, along with the two types of profile breaks. Using a set of controlled hydrodynamic simulations of disc galaxy formation, we find a correlation between the host halo's initial angular momentum and the resulting radial profile of the stellar disc: galaxies that live in haloes with a low spin parameter {\lambda} <~ 0.03 show an up-bending break in their disc density profiles, while galaxies in haloes of higher angular momentum show a down-bending break. We find that the case of pure exponential profiles ({\lambda} ~ 0.035) coincides with the peak of the spin parameter distribution from cosmological simulations. Our simulations not only imply an explanation of the observed behaviours, but also suggest that the physical origin of this effect is related to the amount of radial redistribution of stellar mass, which is anti-correlated with {\lambda}.
We present the first successful simulation of a neutrino-driven supernova explosion in three dimensions (3D), using the Prometheus-Vertex code with an axis-free Yin-Yang grid and a sophisticated treatment of three-flavor, energy-dependent neutrino transport. The progenitor is a non-rotating, zero-metallicity 9.6 Msun star with an iron core. While in spherical symmetry outward shock acceleration sets in later than 300 ms after bounce, a successful explosion starts at ~130 ms post-bounce in two dimensions (2D). The 3D model explodes at about the same time but with faster shock expansion than in 2D and a more quickly increasing and roughly 10 percent higher explosion energy. The more favorable explosion conditions in 3D are explained by lower temperatures and thus reduced neutrino emission in the cooling layer below the gain radius. This moves the gain radius inward and leads to a bigger mass in the gain layer, whose larger recombination energy boosts the explosion energy in 3D. These differences are caused by less coherent, less massive, and less rapid convective downdrafts associated with post-shock convection in 3D. The less violent impact of these accretion downflows in the cooling layer produces less dissipative heating and therefore diminishes energy losses by neutrino emission. We thus have, for the first time, identified a reduced mass accretion rate, lower infall velocities, and a smaller surface filling factor of convective downdrafts as consequences of 3D postshock turbulence that facilitate neutrino-driven explosions and strengthen them compared to the 2D case.
The "green valley" is a wide region separating the blue and the red peaks in the ultraviolet-optical color magnitude diagram, first revealed using GALEX UV photometry. The term was coined by Christopher Martin in 2005. Green valley highlights the discriminating power of UV to very low relative levels of ongoing star formation, to which the optical colors, including u-r, are insensitive. It corresponds to massive galaxies below the star-forming "main" sequence, and therefore represents a critical tool for the study of the quenching of star formation and its possible resurgence in otherwise quiescent galaxies. This article reviews the results pertaining to morphology, structure, environment, dust content and gas properties of green valley galaxies in the local universe. Their relationship to AGN is also discussed. Attention is given to biases emerging from defining the "green valley" using optical colors. We review various evolutionary scenarios and we present evidence for a new, quasi-static view of the green valley, in which the majority of galaxies currently in the green valley were only partially quenched in the distant past and now participate in a slow cosmic decline of star formation, which also drives down the activity on the main sequence, presumably as a result of the dwindling accretion/cooling onto galaxy disks.
We present multi-epoch non-redundant masking observations of the T Cha transition disk, taken at the VLT and Magellan in H, Ks, and L' bands. T Cha is one of a small number of transition disks that host companion candidates discovered by high-resolution imaging techniques, with a putative companion at a position angle of 78 degrees, separation of 62 mas, and contrast at L' of 5.1 mag. We find comparable binary parameters in our re-reduction of the initial detection images, and similar parameters in the 2011 L', 2013 NaCo L', and 2013 NaCo Ks data sets. We find a close-in companion signal in the 2012 NaCo L' dataset that cannot be explained by orbital motion, and a non-detection in the 2013 MagAO/Clio2 L' data. However, Monte-carlo simulations show that the best fits to the 2012 NaCo and 2013 MagAO/Clio2 followup data may be consistent with noise. There is also a significant probability of false non-detections in both of these data sets. We discuss physical scenarios that could cause the best fits, and argue that previous companion and scattering explanations are inconsistent with the results of the much larger dataset presented here.
Since z~1, the stellar mass density locked in low mass groups and clusters has grown by a factor of ~8. Here we make the first statistical measurements of the stellar mass content of low mass X-ray groups at 0.5<z<1, enabling the calibration of stellar-to-halo mass scales for wide-field optical and infrared surveys. Groups are selected from combined Chandra and XMM-Newton X-ray observations in the Chandra Deep Field South (CDFS). These ultra-deep observations allow us to identify bona fide low mass groups at high redshift and enable measurements of their total halo masses. We compute aggregate stellar masses for these halos using galaxies from the Carnegie-Spitzer-IMACS (CSI) spectroscopic redshift survey. Stars comprise ~3-4% of the total mass of group halos with masses 10^{12.8}<M200/Msun<10^{13.5} (about the mass of Fornax and 1/50th the mass of Virgo). Complementing our sample with higher mass halos at these redshifts, we find that the stellar-to-halo mass ratio decreases toward higher halo masses, consistent with other work in the local and high redshift universe. The observed scatter about the stellar-halo mass relation is ~0.25 dex, which is relatively small and suggests that total group stellar mass can serve as a rough proxy for halo mass. We find no evidence for any significant evolution in the stellar-halo mass relation since z<1. Quantifying the stellar content in groups since this epoch is critical given that hierarchical assembly leads to such halos growing in number density and hosting increasing shares of quiescent galaxies.
Proto-galaxies forming in low-mass dark matter haloes are thought to provide the majority of ionising photons needed to reionise the Universe, due to their high escape fractions of ionising photons. We study how the escape fraction in high-redshift galaxies relates to the physical properties of the halo in which the galaxies form by computing escape fractions for 75801 haloes between redshifts 27 and 6 that were extracted from the FiBY project, high-resolution cosmological hydrodynamics simulations of galaxy formation. We find that the main constraint on the escape fraction is the presence of dense gas within 10 pc of the young sources that emit the majority of the ionising photons produced over the lifetime of the stellar population. This results in a strong mass dependence of the escape fraction. The lower potential well in haloes with virial mass below 10^8 solar mass results in lower column densities close to the sources that can be penetrated by the radiation from young, massive stars. In general only a single stellar population forms in these haloes, so supernova feedback sets in too late to strongly affect the escape fraction. In haloes with higher virial mass supernova feedback plays an important role, but only 30% of the haloes in this mass range has an escape fraction higher than 1%. We find a large range of escape fractions in haloes with similar properties, caused by different distributions of the dense gas in the halo. On average the escape fraction of HeI-ionising photons is higher than HI-ionising photons, but almost no HeII-ionising photons escape. Due to the inhomogeneous distribution of the dense gas the escape fraction is highly anisotropic. The strong mass dependence, the large spread and the large anisotropy of the escape fraction may strongly affect the topology of reionisation and is something current models of reionisation should strive to take into account.
NGC 3201 is a globular cluster suspected to have an intrinsic spread in the iron content. We re-analysed a sample of 21 cluster stars observed with UVES-FLAMES at the Very Large Telescope and for which Simmerer et al. found a 0.4 dex wide [Fe/H] distribution with a metal-poor tail. We confirmed that when spectroscopic gravities are adopted, the derived [Fe/H] distribution spans ~0.4 dex. On the other hand, when photometric gravities are used, the metallicity distribution from Fe I lines remains large, while that derived from Fe II lines is narrow and compatible with no iron spread. We demonstrate that the metal-poor component claimed by Simmerer et al. is composed by asymptotic giant branch stars that could be affected by non local thermodynamical equilibrium effects driven by iron overionization. This leads to a decrease of the Fe I abundance, while leaving the Fe II abundance unaltered. A similar finding has been already found in asymptotic giant branch stars of the globular clusters M5 and 47 Tucanae. We conclude that NGC 3201 is a normal cluster, with no evidence of intrinsic iron spread.
We have observed an N-body/Smoothed Particle Hydrodynamics simulation of a Milky Way like barred spiral galaxy. We present a simple method that samples N-body model particles into mock Gaia stellar observations and takes into account stellar populations, dust extinction and Gaia's science performance estimates. We examine the kinematics around a nearby spiral arm at a similar position to the Perseus arm at three lines of sight in the disc plane; (l,b)=(90,0), (120,0) and (150,0) degrees. We find that the structure of the peculiar kinematics around the co-rotating spiral arm, which is found in Kawata et al. (2014b), is still visible in the observational data expected to be produced by Gaia despite the dust extinction and expected observational errors of Gaia. These observable kinematic signatures will enable testing whether the Perseus arm of the Milky Way is similar to the co-rotating spiral arms commonly seen in N-body simulations.
We investigate the impact of sinks of ionizing radiation on the reionization-era 21-cm signal, focusing on 1-point statistics. We consider sinks in both the intergalactic medium and inside galaxies. At a fixed filling factor of HII regions, sinks will have two main effects on the 21-cm morphology: (i) as inhomogeneous absorbers of ionizing photons they result in smaller and more widespread cosmic HII patches; and (ii) as reservoirs of neutral gas they contribute a non-zero 21-cm signal in otherwise ionized regions. Both effects damp the contrast between neutral and ionized patches during reionization, making detection of the epoch of reionization with 21-cm interferometry more challenging. Here we systematically investigate these effects using the latest semi-numerical simulations. We find that sinks dramatically suppress the peak in the redshift evolution of the variance, corresponding to the midpoint of reionization. As previously predicted, skewness changes sign at midpoint, but the fluctuations in the residual HI suppress a late-time rise. Furthermore, large levels of residual HI dramatically alter the evolution of the variance, skewness and power spectrum from that seen at lower levels. In general, the evolution of the large-scale modes provides a better, cleaner, higher signal-to-noise probe of reionization.
In the Local Group, almost all satellite dwarf galaxies that are within the virial radius of the Milky Way (MW) and M31 exhibit strong environmental influence. The orbital histories of these satellites provide the key to understanding the role of the MW/M31 halo, lower-mass groups, and cosmic reionization on the evolution of dwarf galaxies. We examine the virial-infall histories of satellites with M_star = 10 ^ {3 - 9} M_sun using the ELVIS suite of cosmological zoom-in dissipationless simulations of 48 MW/M31-like halos. Satellites at z = 0 fell into the MW/M31 halos typically 5 - 8 Gyr ago at z = 0.5 - 1. However, they first fell into any host halo typically 7 - 10 Gyr ago at z = 0.7 - 1.5. This difference arises because many satellites experienced "group preprocessing" in another host halo, typically of M_vir ~ 10 ^ {10 - 12} M_sun, before falling into the MW/M31 halos. Satellites with lower-mass and/or those closer to the MW/M31 fell in earlier and are more likely to have experienced group preprocessing; half of all satellites with M_star < 10 ^ 6 M_sun were preprocessed in a group. Infalling groups also drive most satellite-satellite mergers within the MW/M31 halos. Finally, none of the surviving satellites at z = 0 were within the virial radius of their MW/M31 halo during reionization (z > 6), and only < 4% were satellites of any other host halo during reionization. Thus, effects of cosmic reionization versus host-halo environment on the formation histories of surviving dwarf galaxies in the Local Group occurred at distinct epochs and are separable in time.
High energy observations of extreme BL Lac objects, such as 1ES 0229+200 or 1ES 0347-121, recently focused interest both for blazar and jet physics and for the implication on the extragalactic background light and intergalactic magnetic field estimate. However, the number of these extreme highly peaked BL Lac objects (EHBL) is still rather small. Aiming at increase their number, we selected a group of EHBL candidates starting from the BL Lac sample of Plotkin et al. (2011), considering those undetected (or only barely detected) by the Large Area Telescope onboard Fermi and characterized by a high X-ray vs. radio flux ratio. We assembled the multi-wavelength spectral energy distribution of the resulting 9 sources, profiting of publicly available archival observations performed by the Swift, Galex and Fermi satellites, confirming their nature. Through a simple one-zone synchrotron self-Compton model we estimate the expected VHE flux, finding that in the majority of cases it is within the reach of present generation of Cherenkov arrays or of the forthcoming Cherenkov Telescope Array (CTA).
We compute the production rate of photons in the ionizing Lyman continua (LyC) of H I (lambda < 912 A), He I (lambda < 504 A), and He II (lambda < 228 A) using recent stellar evolutionary tracks coupled to a grid of non-LTE, line-blanketed (WM-basic) model atmospheres. The median LyC production efficiency is Q_LyC = (6+/-2)x10^60 LyC photons per Msun of star formation (range [3.1-9.4]x10^60) corresponding to a revised calibration of 10^{53.3+/-0.2} photons/s per Msun/yr. Efficiencies in the helium continua are Q_HeI ~ 10^60 photons/Msun and Q_HeII ~ 10^56 photons/Msun at solar metallicity and larger at low metallicity. The critical star formation rate needed to maintain reionization against recombinations at z = 7 is rho_SFR = (0.012 Msun/yr/Mpc^3) [(1+z)/8]^3 [(C_H /3) (0.2/ f_esc) for fiducial values of IGM clumping factor C_H = 3 and LyC escape fraction f_esc = 0.2. The boost in LyC production efficiency is an important ingredient, together with metallicity, C_H, and f_esc, in assessing whether IGM reionization was complete by z ~ 7. Monte-Carlo sampled spectra of coeval starbursts during the first 5 Myr have intrinsic flux ratios of F(1500)/F(900) = 0.4-0.5 and F(912^-)/F(912^+) = 0.4-0.7 in the far-UV (1500A), the LyC (900A), and at the Lyman edge (912A). These ratios can be used to calibrate the LyC escape fractions in starbursts.
We explore the possibility of measuring the mass accretion rate of galaxy clusters by using dense galaxy redshift surveys of their outer regions. By approximating the accretion with the infall of a spherical shell, the mass accretion rate only depends on the mass profile of the cluster in a thin shell at radii larger than $R_{200}$. This approximation is rather crude in hierarchical clustering scenarios, where both smooth accretion and aggregation of smaller dark matter haloes contribute to the mass accretion of clusters. Nevertheless, in the redshift range $z=[0,1]$, our prescription returns an average mass accretion rate within $20 \%$ of the average rate derived with the more realistic merger trees of dark matter haloes extracted from $N$-body simulations. The mass accretion rate of galaxy clusters has been the topic of numerous detailed numerical and theoretical investigations, but so far it has remained inaccessible to measurements in the real Universe. Our result suggests that measuring the mass accretion rate of galaxy clusters is actually feasible, thus providing a potential new observational test of the cosmological and structure formation models.
The M31 globular cluster X-ray binary XB158 (a.k.a. Bo 158) exhibits intensity dips on a 2.78 hr period in some observations, but not others. The short period suggests a low mass ratio, and an asymmetric, precessing disk due to additional tidal torques from the donor star since the disk crosses the 3:1 resonance. Previous theoretical 3D smoothed particle hydrodynamical modeling suggested a super-orbital disk precession period 29$\pm$1 times the orbital period, i.e. $\sim$81$\pm$3 hr. We conducted a Swift monitoring campaign of 30 observations over ~1 month in order to search for evidence of such a super-orbital period. Fitting the 0.3--10 keV Swift XRT luminosity lightcurve with a sinusoid yielded a period of 5.65+/-0.05 days, and a >5$\sigma$ improvement in $\chi^2$ over the best fit constant intensity model. A Lomb-Scargle periodogram revealed that periods 5.4--5.8 days were detected at a >3$\sigma$ level, with a peak at 5.6 days. We consider this strong evidence for a 5.65 day super-orbital period, ~70\% longer than the predicted period. The 0.3--10 keV luminosity varied by a factor ~5, consistent with variations seen in long-term monitoring from Chandra. We conclude that other X-ray binaries exhibiting similar long-term behaviour are likely to also be X-ray binaries with low mass ratios and super-orbital periods.
Several observations of transition discs show lopsided dust-distributions. A potential explanation is the formation of a large-scale vortex acting as a dust-trap at the edge of a gap opened by a giant planet. Numerical models of gap-edge vortices have thus far employed locally isothermal discs, but the theory of this vortex-forming or `Rossby wave' instability was originally developed for adiabatic discs. We generalise the study of planetary gap stability to non-isothermal discs using customised numerical simulations of disc-planet systems where the planet opens an unstable gap. We include in the energy equation a simple cooling function with cooling timescale $t_c=\beta\Omega_k^{-1}$, where $\Omega_k$ is the Keplerian frequency, and examine the effect of $\beta$ on the stability of gap edges and vortex lifetimes. We find increasing $\beta$ lowers the growth rate of non-axisymmetric perturbations, and the dominant azimuthal wavenumber $m$ decreases. We find a quasi-steady state consisting of one large-scale, over-dense vortex circulating the outer gap edge, typically lasting $O(10^3)$ orbits. Vortex lifetimes were found to generally increase with cooling times up to an optimal value, beyond which vortex lifetimes decrease. This non-monotonic dependence is qualitatively consistent with recent studies using strictly isothermal discs that vary the disc aspect ratio.
We present a suite of cosmological radiation-hydrodynamical simulations of the assembly of galaxies driving the reionization of the intergalactic medium (IGM) at z >~ 6. The simulations account for the hydrodynamical feedback from photoionization heating and the explosion of massive stars as supernovae (SNe). Our reference simulation, which was carried out in a box of size 25 comoving Mpc/h using 2 x 512^3 particles, produces a reasonable reionization history and matches the observed UV luminosity function of galaxies. Simulations with different box sizes and resolutions are used to investigate numerical convergence, and simulations in which either SNe or photoionization heating or both are turned off, are used to investigate the role of feedback from star formation. Ionizing radiation is treated using accurate radiative transfer at the high spatially adaptive resolution at which the hydrodynamics is carried out. SN feedback strongly reduces the star formation rates (SFRs) over nearly the full mass range of simulated galaxies and is required to yield SFRs in agreement with observations. Photoheating helps to suppress star formation in low-mass galaxies, but its impact on the cosmic SFR is small. Because the effect of photoheating is masked by the strong SN feedback, it does not imprint a signature on the UV galaxy luminosity function. Photoheating does provide a strong positive feedback on reionization because it smooths density fluctuations in the IGM, which lowers the IGM recombination rate substantially. Our simulations demonstrate a tight non-linear coupling of galaxy formation and reionization, motivating the need for the accurate and simultaneous inclusion of photoheating and SN feedback in models of the early Universe.
We review results from cosmic X-ray surveys of active galactic nuclei (AGNs) over the past ~ 15 yr that have dramatically improved our understanding of growing supermassive black holes (SMBHs) in the distant universe. First, we discuss the utility of such surveys for AGN investigations and the capabilities of the missions making these surveys, emphasizing Chandra, XMM-Newton, and NuSTAR. Second, we briefly describe the main cosmic X-ray surveys, the essential roles of complementary multiwavelength data, and how AGNs are selected from these surveys. We then review key results from these surveys on the AGN population and its evolution ("demographics"), the physical processes operating in AGNs ("physics"), and the interactions between AGNs and their environments ("ecology"). We conclude by describing some significant unresolved questions and prospects for advancing the field.
In addition to producing a strong gravitational signal, a short gamma-ray burst (GRB), and a compact remnant, neutron star mergers eject significant masses at significant kinetic energies. This mass ejection takes place via dynamical mass ejection and a GRB jet but other processes have also been suggested: a shock-breakout material, a cocoon resulting from the interaction of the jet with other ejecta, and viscous and neutrino driven winds from the central remnant or the accretion disk. The different components of the ejected masses include up to a few percent of a solar mass, some of which is ejected at relativistic velocities. The interaction of these ejecta with the surrounding interstellar medium will produce a long lasting radio flare, in a similar way to GRB afterglows or to radio supernovae. The relative strength of the different signals depends strongly on the viewing angle. An observer along the jet axis or close to it will detect a strong signal at a few dozen days from the radio afterglow (or the orphan radio afterglow) produced by the highly relativistic GRB jet. For a generic observer at larger viewing angles, the dynamical ejecta, whose contribution peaks a year or so after the event, will generally dominate. Depending on the observed frequency and the external density, other components may also give rise to a significant contribution. We also compare these estimates with the radio signature of the short GRB 130603B. The radio flare from the dynamical ejecta might be detectable with the EVLA and the LOFAR for the higher range of external densities $n\gtrsim 0.5{\rm cm^{-3}}$.
The starbusting, nearby (D = 32.9 Mpc) spiral (Sc) galaxy NGC2276 belongs to the sparse group dominated by the elliptical galaxy NGC2300. NGC2276 is a remarkable galaxy, as it displays a disturbed morphology at many wavelengths. This is possibly due to gravitational interaction with the central elliptical galaxy of the group. Previous ROSAT and XMM-Newton observations resulted in the detection of extended hot gas emission and of a single very bright (~1.e41 erg/s) ultraluminous X-ray source (ULX) candidate. Here we report on a study of the X-ray sources of NGC2276 based on CHANDRA data taken in 2004. CHANDRA was able to resolve 16 sources, 8 of which are ULXs, and to reveal that the previous ULX candidate is actually composed of a few distinct objects. We construct the luminosity function of NGC2276, which can be interpreted as dominated by high mass X-ray binaries, and estimate the star formation rate (SFR) to be ~5-15 Msun/yr, consistent with the values derived from optical and infrared observations. By means of numerical simulations, we show that both ram pressure and viscous transfer effects are necessary to produce the distorted morphology and the high SFR observed in NGC2276, while tidal interaction have a marginal effect.
The first part of this article is a historical and physical introduction to quasars and their close cousins, called Active Galactic Nuclei (AGN). In the second part, I argue that our progress in understanding them has been unsatisfactory and in fact somewhat illusory since their discovery fifty years ago, and that much of the reason is a pervasive lack of critical thinking in the research community. It would be very surprising if other fields do not suffer similar failings.
We present the third Fermi Large Area Telescope source catalog (3FGL) of sources in the 100~MeV--300~GeV range. Based on the first four years of science data from the Fermi Gamma-ray Space Telescope mission, it is the deepest yet in this energy range. Relative to the 2FGL catalog, the 3FGL catalog incorporates twice as much data as well as a number of analysis improvements, including improved calibrations at the event reconstruction level, an updated model for Galactic diffuse gamma-ray emission, a refined procedure for source detection, and improved methods for associating LAT sources with potential counterparts at other wavelengths. The 3FGL catalog includes 3033 sources above 4 sigma significance, with source location regions, spectral properties, and monthly light curves for each. Of these, 78 are flagged as potentially being due to imperfections in the model for Galactic diffuse emission. Twenty-five sources are modeled explicitly as spatially extended, and overall 232 sources are considered as identified based on angular extent or correlated variability (periodic or otherwise) observed at other wavelengths. For 1009 sources we have not found plausible counterparts at other wavelengths. More than 1100 of the identified or associated sources are active galaxies of the blazar class; several other classes of non-blazar active galaxies are also represented in the 3FGL. Pulsars represent the largest Galactic source class. From source counts of Galactic sources we estimate the contribution of unresolved sources to the Galactic diffuse emission is ~3% at 1 GeV.
We have carried out 200,000 N-body simulations of three identical stellar embryos with masses from a Chabrier IMF and embedded in a molecular core. The bodies are initially non-hierarchical and undergo chaotic motions, while accreting using Bondi-Hoyle accretion. The coupling of dynamics and accretion often leads to one or two dominant bodies controlling the center of the cloud core, while banishing the other(s) to the lower-density outskirts, leading to stunted growth. Eventually each system transforms either to a bound hierarchical configuration or breaks apart into separate single and binary components. The orbital motion is followed for 100 Myr. To illustrate the simulations we introduce the 'triple diagnostic diagram', which plots two dimensionless numbers against each other, representing the binary mass ratio and the mass ratio of the third body to the total system mass. Numerous freefloating BD binaries are formed in these simulations. The separation distribution function is in good correspondence with observations, showing a steep rise at close separations, peaking around 13 AU and then declining more gently. Unresolved BD triple systems may appear as wider BD binaries. Mass ratios are strongly peaked towards unity, as observed, but this is partially due to the initial assumptions. Eccentricities gradually increase towards higher values, due to the lack of viscous interactions in the simulations, which would both shrink the orbits and decrease their eccentricities. The main threat to newly born triple systems is internal instabilities, not external perturbations. Dynamical interactions in newborn triple systems of stellar embryos embedded in and accreting from a cloud core naturally form a population of freefloating BD binaries, and this mechanism may constitute a significant pathway for the formation of BD binaries. (Abstract abbreviated).
Tidal disruption events (TDE) in which a star is devoured by a massive black hole at a galac- tic center pose a challenge to our understanding of accretion processes. Within a month the accretion rate reaches super-Eddington levels. It then drops gradually over a time scale of a year to sub-Eddington regimes. The initially geometrically thick disk becomes a thin one and eventually an ADAF at very low accretion rates. As such, TDEs explore the whole range of accretion rates and configurations. A challenging question is what the corresponding light curves of these events are. We explore numerically the disk luminosity and the conditions within the inner region of the disk using a fully general relativistic slim disk model. Those conditions determine the magnitude of the magnetic field that engulfs the black hole and this, in turn, determines the Blandford-Znajek jet power. We estimate this power in two different ways and show that they are self-consistent. We find, as expected earlier from analytic argu- ments (Krolik & Piran 2012), that neither the disk luminosity nor the jet power follows the accretion rate throughout the disruption event. The disk luminosity varies only logarithmi- cally with the accretion rate at super-Eddington luminosities. The jet power follows initially the accretion rate but remains a constant after the transition from super- to sub- Eddington. At lower accretion rates at the end of the MAD phase the disk becomes thin and the jet may stop altogether. These new estimates of the jet power and disk luminosity that do not simply follow the mass fallback rate should be taken into account when searching for TDEs and analysing light curves of TDE candidates. Identification of some of the above mentioned transitions may enable us to estimate better TDE parameters.
We present Gemini spectroscopy for 21 candidate optical counterparts to X-ray sources discovered in the Galactic Bulge Survey (GBS). For the majority of the 21 sources, the optical spectroscopy establishes that they are indeed the likely counterparts. One of the criteria we used for the identification was the presence of an Ha emission line. The spectra of several sources revealed an Ha emission line only after careful subtraction of the F or G stellar spectral absorption lines. In a sub-class of three of these sources the residual Halpha emission line is broad (> 400 km/s) which suggests that it is formed in an accretion disk, whereas in other cases the line width is such that we currently cannot determine whether the line emission is formed in an active star/binary or in an accretion disk. GBS source CX377 shows this hidden accretion behaviour most dramatically. The previously-identified broad Ha emission of this source is not present in its Gemini spectra taken about 1 year later. However, broad emission is revealed after subtracting an F6 template star spectrum. The Gemini spectra of three sources (CX446, CX1004, and CXB2) as well as the presence of possible eclipses in light curves of these sources suggest that these sources are accreting binaries viewed under a high inclination.
Here we describe the \emph{Siding Spring Southern Seyfert Spectroscopic Snapshot Survey} (S7) and present results on 64 galaxies drawn from the first data release. The S7 uses the Wide Field Spectrograph (WiFeS) mounted on the ANU 2.3m telescope located at the Siding Spring Observatory to deliver an integral field of $38\times25$~ arcsec at a spectral resolution of $R=7000$ in the red ($530-710$nm), and $R=3000$ in the blue ($340-560$nm). {From these data cubes we have extracted the Narrow Line Region (NLR) spectra from a 4 arc sec aperture centred on the nucleus. We also determine the H$\beta$ and [OIII]~$\lambda$5007 fluxes in the narrow lines, the nuclear reddening, the reddening-corrected relative intensities of the observed emission lines, and the H$\beta$ and \lOIII\ luminosities {determined from spectra for which the stellar continuum has been removed.} We present a set of images of the galaxies in [OIII]~$\lambda$5007, [NII]~$\lambda$6584 and H$\alpha$ which serve to delineate the spatial extent of the extended narrow line region (ENLR) and {\bf also to} reveal the structure and morphology of the surrounding \HII\ regions. Finally, we provide a preliminary discussion of those Seyfert~1 and Seyfert~2 galaxies which display coronal emission lines in order to explore the origin of these lines.
We use observations from the Compact Reconnaissance Imaging Spectrometer for
Mars (CRISM) of the north polar cap during late summer for two Martian years,
to monitor the complete summer cycle of albedo and water ice grain size in
order to place quantitative limits of the amount of water ice deposited in late
summer.
We establish here for the first time the complete spring to summer cycle of
water ice grain sizes on the north polar cap. The apparent grain sizes grow
until Ls=132, when they appear to shrink again, until they are obscured at the
end of summer by the north polar hood.
Under the assumption that the shrinking of grain sizes is due to the
deposition of find grained ice, we quantify the amount of water ice deposited
per Martian boreal summer, and estimate the amount of water ice that must be
transported equatorward.
Interestingly, we find that the relative amount of water ice deposited in the
north cap during boreal summer (0.7-7 microns) is roughly equivalent to the
average amount of water ice deposited on the south polar cap during austral
summer (0.6-6 microns).
Recently, an excess of GeV gamma ray near the Galactic Centre has beenvreported. The spectrum obtained can be best fitted with the annihilationvof $30-40$ GeV dark matter particles through $b \bar{b}$ channel. In this letter, I show that this annihilation model can also solve the mysteries of heating source in x-ray plasma and the unexpected high gamma-ray luminosity. The cross section constrained by these observations give excellent agreements with both the predicted range by using Fermi-LAT data and the canonical thermal relic abundance cross section.
We present a proof-of-concept analysis of photometric redshifts with Bayesian priors on physical properties of galaxies. This concept is particularly suited for upcoming/on-going large imaging surveys, in which only several broad-band filters are available and it is hard to break some of the degeneracies in the multi-color space. We construct model templates of galaxies using a stellar population synthesis code and apply Bayesian priors on physical properties such as stellar mass and star formation rate. These priors are a function of redshift and they effectively evolve the templates with time in an observationally motivated way. We demonstrate that the priors help reduce the degeneracy and deliver significantly improved photometric redshifts. Furthermore, we show that a template error function, which corrects for systematic flux errors in the model templates as a function of rest-frame wavelength, delivers further improvements. One great advantage of our technique is that we simultaneously measure redshifts and physical properties of galaxies in a fully self-consistent manner, unlike the two-step measurements with different templates often performed in the literature. One may rightly worry that the physical priors bias the inferred galaxy properties, but we show that the bias is smaller than systematic uncertainties inherent in physical properties inferred from the SED fitting and hence is not a major issue. We will extensively test and tune the priors in the on-going Hyper Suprime-Cam survey and will make the code publicly available in the future.
We study the impact of baryonic physics on cosmological parameter estimation with weak lensing surveys. We run a set of cosmological hydrodynamics simulations with different galaxy formation models. We then perform ray-tracing simulations through the total matter density field to generate 100 independent convergence maps of 25 deg$^2$ field-of-view, and use them to examine the ability of the following three lensing statistics as cosmological probes; power spectrum, peak counts, and Minkowski Functionals. For the upcoming wide-field observations such as Subaru Hyper Suprime-Cam (HSC) survey with a sky coverage of 1400 deg$^2$, the higher-order statistics provide tight constraints on the matter density, density fluctuation amplitude, and dark energy equation of state, but appreciable parameter bias is induced by the baryonic processes such as gas cooling and stellar feedback. When we use power spectrum, peak counts, and Minkowski Functionals, the relative bias in the dark energy equation of state parameter $w$ is at a level of, respectively, $\sim0.06\sigma$, $0.5-0.6\sigma$, and $0.01-0.1\sigma$ where $\sigma$ is the overall error derived from Fisher analysis. We find the bias is induced in different directions in the parameter space depending on the statistics employed. While the two-point statistics, i.e. power spectrum, yield robust results against baryonic effects, the overall constraining power is weak compared with the other higher-order statistics. On the other hand, using higher-order statistics alone results in significantly biased parameter estimate. We suggest to use an optimized combination of, for example, power spectrum and higher-order statistics so that the baryonic effects on parameter estimation are mitigated. Such `calibrated' combination can place stringent and robust constraints on cosmological parameters.
The standard approach for time-resolved X-ray spectral analysis of
thermonuclear bursts involves subtraction of the pre-burst emission as
background. This approach implicitly assumes that the persistent flux remains
constant throughout the burst. We reanalyzed 332 photospheric radius expansion
bursts observed from 40 sources by the Rossi X-ray Timing Explorer, introducing
a multiplicative factor $f_a$ to the persistent emission contribution in our
spectral fits. We found that for the majority of spectra the best-fit value of
$f_a$ is significantly greater than 1, suggesting that the persistent emission
typically increases during a burst. Elevated $f_a$ values were not found solely
during the radius expansion interval of the burst, but were also measured in
the cooling tail. The modified model results in a lower average value of the
$\chi^2$ fit statistic, indicating superior spectral fits, but not yet to the
level of formal statistical consistency for all the spectra.
We interpret the elevated $f_a$ values as an increase of the mass accretion
rate onto the neutron star during the burst, likely arising from the effects of
Poynting-Robertson drag on the disk material. We measured an inverse
correlation of $f_a$ with the persistent flux, consistent with theoretical
models of the disc response. We suggest that this modified approach may provide
more accurate burst spectral parameters, as well as offering a probe of the
accretion disk structure.
In the last decade, enormous progress has been achieved in the understanding of the various facets of coalescing double neutron star and neutron black hole binary systems. One hopes that the mergers of such compact binaries can be routinely detected with the advanced versions of the ground-based gravitational wave detector facilities, maybe as early as in 2016. From the theoretical side, there has also been mounting evidence that compact binary mergers could be major sources of heavy elements and these ideas have gained recent observational support from the detection of an event that has been interpreted as a "macronova", an electromagnetic transient powered by freshly produced, radioactively decaying heavy elements. In addition, compact binaries are the most plausible triggers of short gamma-ray bursts (sGRBs) and the last decade has witnessed the first detection of a sGRB afterglow and subsequent observations have delivered a wealth of information on the environments in which such bursts occur. To date, compact binary mergers can naturally explain most --though not all-- of the observed sGRB properties. This article reviews major recent developments in various areas related to compact binary mergers.
Some extrasolar giant planets in close orbits---"hot Jupiters"---exhibit larger radii than that of a passively cooling planet. The extreme irradiation $L_{\rm eq}$ these hot Jupiters receive from their close in stars creates a thick isothermal layer in their envelopes, which slows down their convective cooling, allowing them to retain their inflated size for longer. This is yet insufficient to explain the observed sizes of the most inflated planets. Some models invoke an additional power source, deposited deep in the planet's envelope. Here we present an analytical model for the cooling of such irradiated, and internally heated gas giants. We show that a power source $L_{\rm dep}$, deposited at an optical depth $\tau_{\rm dep}$, creates an exterior convective region, between optical depths $L_{\rm eq}/L_{\rm dep}$ and $\tau_{\rm dep}$, beyond which a thicker isothermal layer exists, which in extreme cases may extend to the center of the planet. This convective layer, which occurs only for $L_{\rm dep}\tau_{\rm dep}>L_{\rm eq}$, further delays the cooling of the planet. Such a planet is equivalent to a planet irradiated with $L_{\rm eq}\left(1+L_{\rm dep}\tau_{\rm dep}/L_{\rm eq}\right)^\beta$, where $\beta\approx 0.35$ is an effective power-law index describing the radiative energy density as function of the optical depth for a convective planet $U\propto\tau^\beta$. Our simple analytical model reproduces the main trends found in previous numerical works, and provides an intuitive understanding. We derive scaling laws for the cooling rate of the planet, its central temperature, and radius. These scaling laws can be used to estimate the effects of tidal or Ohmic dissipation, wind shocks, or any other mechanism involving energy deposition, on sizes of hot Jupiters.
We have assembled a sample of 5 X-ray-absorbed and submm-luminous type 1 QSOs at $z \sim 2$ which are simultaneously growing their central black holes through accretion and forming stars copiously. We present here the analysis of their rest-frame UV to submm Spectral Energy Distributions (SEDs), including new Herschel data. Both AGN (direct and reprocessed) and Star Formation (SF) emission are needed to model their SEDs. From the SEDs and their UV-optical spectra we have estimated the masses of their black holes $M_{BH}\sim 10^{9}-10^{10}\,M_{\odot}$, their intrinsic AGN bolometric luminosities $L_{BOL}\sim(0.8 - 20)\times 10^{13} L_{\odot}$, Eddington ratios $L_{BOL}/L_{Edd}\sim 0.1 - 1.1$ and bolometric corrections $L_{BOL}/L_{X,2-10}\sim 30 - 500$. These values are common among optically and X-ray-selected type 1 QSOs (except for RX~J1249), except for the bolometric corrections, which are higher. These objects show very high far-infrared luminosities $L_{FIR}\sim$ (2 - 8)$\times10^{12}\,M_{\odot}$ and Star Formation Rates SFR$\sim 1000 M_{\odot}/$y. From their $L_{FIR}$ and the shape of their FIR-submm emission we have estimated star-forming dust masses of $M_{DUST}\sim 10^9\,M_\odot$. We have found evidence of a tentative correlation between the gas column densities of the ionized absorbers detected in X-ray (N$_{H_{ion}}$) and $SFR$. Our computed black hole masses are amongst the most massive known.
Multi-Conjugate Adaptive Optics systems based on sodium Laser Guide Stars may exploit Natural Guide Stars to solve intrinsic limitations of artificial beacons (tip-tilt indetermination and anisoplanatism) and to mitigate the impact of the sodium layer structure and variability. The sodium layer may also have transverse structures leading to differential effects among Laser Guide Stars. Starting from the analysis of the input perturbations related to the Sodium Layer variability, modeled directly on measured sodium layer profiles, we analyze, through a simplified end-to-end simulation code, the impact of the low/medium orders induced on global performance of the European Extremely Large Telescope Multi-Conjugate Adaptive Optics module MAORY.
This is continuation of our programme to search for the elusive radio-quiet BL Lacs, by carrying out a systematic search for intranight optical variability (INOV) in a subset of `weak-line quasars' which are already designated as `high-confidence BL Lac candidate' and are also known to be radio-quiet. For 6 such radio-quiet weak-line quasars (RQWLQs), we present here new INOV observations taken in 11 sessions of duration >3 hours each. Combining these data with our previously published INOV monitoring of RQWLQs in 19 sessions yields INOV observations for a set of 15 RQWLQs monitored in 30 sessions, each lasting more than 3 hours. The 30 differential light curves, thus obtained for the 15 RQWLQs, were subjected to a statistical analysis using the F-test, and the deduced INOV characteristics of the RQWLQs then compared with those published recently for several prominent AGN classes, also applying the F-test. From our existing INOV observations, there is a hint that RQWLQs in our sample show a significantly higher INOV duty cycle than radio-quiet quasars and radio lobe-dominated quasars. Two sessions when we have detected strong (blazar-like) INOV for RQWLQs are pointed out, and these two RQWLQs are therefore the best known candidates for radio-quiet BL Lacs, deserving to be pursued. For a proper comparison with the INOV properties already established for (brighter) members of several prominent classes of AGN, a factor of 2-3 improvement in the INOV detection threshold for the RQWLQs is needed and it would be very interesting to check if that would yield a significantly higher estimate for INOV duty cycle than is found here.
Measurement of Cosmic Microwave Background (CMB) anisotropies has been playing a lead role in precision cosmology by providing some of the tightest constrains on cosmological models and parameters. However, precision can only be meaningful when all major systematic effects are taken into account. Non-circular beams in CMB experiments can cause large systematic deviation in the angular power spectrum, not only by modifying the measurement at a given multipole, but also introducing coupling between different multipoles through a deterministic bias matrix. Here we add a mechanism for emulating the effect of a full bias matrix to the Planck likelihood code through the parameter estimation code SCoPE. We show that if the angular power spectrum was measured with a non-circular beam, the assumption of circular Gaussian beam or considering only the diagonal part of the bias matrix can lead to huge error in parameter estimation. We demonstrate that, at least for elliptical Gaussian beams, use of scalar beam window functions obtained via Monte Carlo simulations starting from a fiducial spectrum, as implemented in Planck analyses for example, leads to em only few percent of sigma deviation of the best-fit parameters. However, we notice more significant differences in the posterior distributions for some of the parameters, which would in turn lead to incorrect errorbars. These differences can be reduced, so that the errorbars match within few percent, by adding an iterative reanalysis step, where the beam window function would be recomputed using the best-fit spectrum estimated in the first step.
Detection of supernovae (SNe) and, more generally, of transient events in large surveys can provide numerous false detections. In the case of a deferred processing of survey images, this implies reconstructing complete light curves for all detections, requiring sizable processing time and resources. Optimizing the detection of transient events is thus an important issue for both present and future surveys. We present here the optimization done in the SuperNova Legacy Survey (SNLS) for the 5-year data differed photometric analysis. In this analysis, detections are derived from stacks of subtracted images with one stack per lunation. The 3-year analysis provided 300,000 detections dominated by signals of bright objects that were not perfectly subtracted. We developed a subtracted image stack treatment to reduce the number of non SN-like events using morphological component analysis. This technique exploits the morphological diversity of objects to be detected to extract the signal of interest. At the level of our subtraction stacks, SN-like events are rather circular objects while most spurious detections exhibit different shapes. A two-step procedure was necessary to have a proper evaluation of the noise in the subtracted image stacks and thus a reliable signal extraction. We also set up a new detection strategy to obtain coordinates with good resolution for the extracted signal. SNIa Monte-Carlo (MC) generated images were used to study detection efficiency and coordinate resolution. When tested on SNLS 3-year data this procedure decreases the number of detections by two, losing only 5% of SN-like events, all faint ones. MC results show that SNIa detection efficiency is equivalent to that of the original method for bright events, while the coordinate resolution is improved.
General relativistic magnetohydrodynamic (GRMHD) flows along magnetic fields threading a black hole can be divided into inflow and outflow part, according to the result of the competition between the black hole gravity and magneto-centrifugal forces along the field line. Here we present the first self-consistent, semi-analytical solution for a cold, Poynting flux-dominated (PFD) GRMHD flow, which pass all four critical (inner and outer, Alfv\'en and fast-magnetosonic) points along a parabolic streamline. By assuming that the dominating (electromagnetic) component of the energy flux per flux tube is conserved at the surface where the inflow and outflow are separated, the outflow part of the solution can be constraint by the inflow part of the solution.The semi-analytical method can provide fiducial and complementary solutions for GRMHD simulations around the rotating black hole, given that the black hole spin, global streamline, and magnetizaion (i.e., a mass-loading at the inflow/outflow separation) are prescribed. For reference, we demonstrate a self-consistent result with the work by McKinney in a quantitative level.
Multiwavelength variability data, combined with spectral-timing analysis techniques, provides information about the causal relationship between different physical components in accreting black holes. Using fast-timing data and long-term monitoring, we can probe the behaviour of the same components across the black hole mass scale. In this chapter we review the observational status of multiwavelength variability in accreting black holes, from black hole X-ray binaries to AGN, and consider the implications for models of accretion and ejection, primarily considering the evidence for accretion disc and jet variability in these systems. We end with a consideration of future prospects in this quickly-developing field.
Inspired on the well known dynamical dichotomy predicted in voids, where some underdense regions expand whereas others collapse due to overdense surrounding regions, we explored the interplay between the void inner dynamics and its large scale environment. The environment is classified depending on its density as in previous works. We analyse the dynamical properties of void-centered spherical shells at different void-centric distances depending on this classification. The above dynamical properties are given by the angular distribution of the radial velocity field, its smoothness, the field dependence on the tracer density and shape, and the field departures from linear theory. We found that the velocity field in expanding voids follows more closely the linear prediction, with a more smooth velocity field. However when using velocity tracers with large densities such deviations increase. Voids with sizes around $18\,h^{-1}\,Mpc$ are in a transition regime between regions with expansion overpredicted and underpredicted from linear theory. We also found that velocity smoothness increases as the void radius, indicating the laminar flow dominates the expansion of larger voids (more than $18\,h^{-1}\,Mpc$). The correlations observed suggest that nonlinear dynamics of the inner regions of voids could be dependent on the evolution of the surrounding structures. These also indicate possible scale couplings between the void inner expansion and the large scale regions where voids are embedded. These results shed some light to the origin of nonlinearities in voids, going beyond the fact that voids just quickly becomes nonlinear as they become emptier.
Strange quark matter (SQM) may be the true ground state of hadronic matter, indicating that the observed pulsars may actually be strange stars, but not neutron stars. According to this SQM hypothesis, the existence of a hydrostatically stable sequence of strange quark matter stars has been predicted, ranging from 1 --- 2 solar mass strange stars, to smaller strange dwarfs and even strange planets. While gravitational wave (GW) astronomy is expected to open a new window to the universe, it will shed light on the searching for SQM stars. Here we show that due to their extreme compactness, strange planets can spiral very close to their host strange stars, without being tidally disrupted. Like inspiraling neutron stars or black holes, these systems would serve as a new kind of sources for GW bursts, producing strong gravitational waves at the final stage. The events occurring in our local Universe can be detected by the upcoming gravitational wave detectors, such as Advanced LIGO and the Einstein Telescope. This effect provides a unique probe to SQM objects and is hopefully a powerful tool for testing the SQM hypothesis.
We study the Integrated Sachs-Wolfe (ISW) effect in ghost-free, massive bigravity, where only one metric couples to matter. We focus on the infinite-branch bigravity (IBB) model which exhibits viable cosmic expansion histories and stable linear perturbations, while the cosmological constant is set to zero and the late-time accelerated expansion of the Universe is due solely to the gravitational interaction terms. The ISW contribution to the CMB auto-correlation power spectrum is predicted, as well as the cross-correlation between the CMB temperature anisotropies and the large-scale structure. We use ISW amplitudes as observed in the WMAP 9-year temperature data together with galaxy and AGN data provided by the WISE mission, in order to compare the theoretical predictions to the observations. The ISW amplitudes in IBB are found to be larger than the corresponding ones in the standard LCDM model by roughly a factor of 1.5, but are still consistent with the observations.
Many atmospheric and climatic criteria have to be taken into account for the selection of a suitable site for the next generation of imaging air-shower Cherenkov telescopes, the "Cherenkov Telescope Array" CTA. Such data are not available with sufficient precision or the comparability to allow for a comprehensive characterization of the proposed sites to be made. Identical cross-calibrated instruments have been developed which allow for precise comparison between sites, the cross-validation of existing data, and the ground-validation of satellite data. The site characterization work package of the CTA consortium opted to construct and deploy 9 copies of an autonomous multi-purpose weather sensor, incorporating an infrared cloud sensor a newly developed sensor for measuring the light of the night sky, and an All-Sky-Camera, the whole referred to as Autonomous Tool for Measuring Observatory Site COnditions PrEcisely (ATMOSCOPE). We present here the hardware that was combined into the ATMOSCOPE and characterize its performance.
We studied the radio emission from four radio-loud and gamma-ray-loud narrow-line Seyfert 1 galaxies. The goal was to investigate whether a relativistic jet is operating at the source, and quantify its characteristics. We relied on the most systematic monitoring of such system in the cm and mm radio bands which is conducted with the Effelsberg 100 m and IRAM 30 m telescopes and covers the longest time-baselines and the most radio frequencies to date. We extract variability parameters and compute variability brightness temperatures and Doppler factors. The jet powers were computed from the light curves to estimate the energy output. The dynamics of radio spectral energy distributions were examined to understand the mechanism causing the variability. All the sources display intensive variability that occurs at a pace faster than what is commonly seen in blazars. The flaring events show intensive spectral evolution indicative of shock evolution. The brightness temperatures and Doppler factors are moderate, implying a mildly relativistic jet. The computed jet powers show very energetic flows. The radio polarisation in one case clearly implies a quiescent jet underlying the recursive flaring activity. Despite the generally lower flux densities, the sources appear to show all typical characteristics seen in blazars that are powered by relativistic jets.
We summarize basic observational results on Sagittarius~A* obtained from the radio, infrared and X-ray domain. Infrared observations have revealed that a dusty S-cluster object (DSO/G2) passes by SgrA*, the central super-massive black hole of the Milky Way. It is still expected that this event will give rise to exceptionally intense activity in the entire electromagnetic spectrum. Based on February to September 2014 SINFONI observations. The detection of spatially compact and red-shifted hydrogen recombination line emission allows a new estimate of the orbital parameters of the DSO. We have not detected strong pre-pericenter blue-shifted or post-pericenter red-shifted emission above the noise level at the position of SgrA* or upstream the orbit. The periapse position was reached in May 2014. Our 2004-2012 infrared polarization statistics shows that SgrA* must be a very stable system - both in terms of geometrical orientation of a jet or accretion disk and in terms of the variability spectrum which must be linked to the accretion rate. Hence polarization and variability measurements are the ideal tool to probe for any change in the system as a function of the DSO/G2 fly-by. Due to the 2014 fly-by of the DSO, increased accretion activity of SgrA* may still be upcoming. Future observations of bright flares will improve the derivation of the spin and the inclination of the SMBH from NIR/sub-mm observations.
The chemical composition of exoplanet host stars is an important factor in understanding the formation and characteristics of their orbiting planets. The best example of this to date is the planet-metallicity correlation. Other proposed correlations are thus far less robust, in part due to uncertainty in the chemical history of stars pre- and post-planet formation. Binary host stars of similar type present an opportunity to isolate the effects of planets on host star abundances. Here we present a differential elemental abundance analysis of the XO-2 stellar binary, in which both G9 stars host giant planets, one of which is transiting. Building on our previous work, we report 16 elemental abundances and compare the $\Delta$(XO-2N-XO-S) values to elemental condensation temperatures. The $\Delta$(N-S) values and slopes with condensation temperature resulting from four different pairs of stellar parameters are compared to explore the effects of changing the relative temperature and gravity of the stars. We find that most of the abundance differences between the stars depend on the chosen stellar parameters, but that Fe, Si, and potentially Ni are consistently enhanced in XO-2N regardless of the chosen stellar parameters. This study emphasizes the power of binary host star abundance analysis for probing the effects of giant planet formation, but also illustrates the potentially large uncertainties in abundance differences and slopes induced by changes in stellar temperature and gravity.
We discuss the potential of a next generation space-borne CMB experiment for studies of extragalactic sources with reference to COrE+, a project submitted to ESA in response to the M4 call. We consider three possible options for the telescope size: 1m, 1.5m and 2m (although the last option is probably impractical, given the M4 boundary conditions). The proposed instrument will be far more sensitive than Planck and will have a diffraction-limited angular resolution. These properties imply that even the 1m telescope option will perform substantially better than Planck for studies of extragalactic sources. The source detection limits as a function of frequency have been estimated by means of realistic simulations. The most significant improvements over Planck results are presented for each option. COrE+ will provide much larger samples of truly local star-forming galaxies, making possible analyses of the properties of galaxies (luminosity functions, dust mass functions, star formation rate functions, dust temperature distributions, etc.) across the Hubble sequence. Even more interestingly, COrE+ will detect, at |b|> 30 deg, thousands of strongly gravitationally lensed galaxies. Such large samples are of extraordinary astrophysical and cosmological value in many fields. Moreover, COrE+ high frequency maps will be optimally suited to pick up proto-clusters of dusty galaxies, i.e. to investigate the evolution of large scale structure at larger redshifts than can be reached by other means. Thanks to its high sensitivity COrE+ will also yield a spectacular advance in the blind detection of extragalactic sources in polarization. This will open a new window for studies of radio source polarization and of the global properties of magnetic fields in star forming galaxies and of their relationships with SFRs.
The super-massive 4 million solar mass black hole Sagittarius~A* (SgrA*) shows flare emission from the millimeter to the X-ray domain. A detailed analysis of the infrared light curves allows us to address the accretion phenomenon in a statistical way. The analysis shows that the near-infrared flare amplitudes are dominated by a single state power law, with the low states in SgrA* limited by confusion through the unresolved stellar background. There are several dusty objects in the immediate vicinity of SgrA*. The source G2/DSO is one of them. Its nature is unclear. It may be comparable to similar stellar dusty sources in the region or may consist predominantly of gas and dust. In this case a particularly enhanced accretion activity onto SgrA* may be expected in the near future. Here the interpretation of recent data and ongoing observations are discussed.
We have observed the HII region RCW175 with the 64m Parkes telescope at 8.4GHz and 13.5GHz in total intensity, and at 21.5GHz in both total intensity and polarization. High angular resolution, high sensitivity, and polarization capability enable us to perform a detailed study of the different constituents of the HII region. For the first time, we resolve three distinct regions at microwave frequencies, two of which are part of the same annular diffuse structure. Our observations enable us to confirm the presence of anomalous microwave emission (AME) from RCW175. Fitting the integrated flux density across the entire region with the currently available spinning dust models, using physically motivated assumptions, indicates the presence of at least two spinning dust components: a warm component with a relatively large hydrogen number density n_H=26.3/cm^3 and a cold component with a hydrogen number density of n_H=150/cm^3. The present study is an example highlighting the potential of using high angular-resolution microwave data to break model parameter degeneracies. Thanks to our spectral coverage and angular resolution, we have been able to derive one of the first AME maps, at 13.5GHz, showing clear evidence that the bulk of the AME arises in particular from one of the source components, with some additional contribution from the diffuse structure. A cross-correlation analysis with thermal dust emission has shown a high degree of correlation with one of the regions within RCW175. In the center of RCW175, we find an average polarized emission at 21.5GHz of 2.2\pm0.2(rand.)\pm0.3(sys.)% of the total emission, where we have included both systematic and statistical uncertainties at 68% CL. This polarized emission could be due to sub-dominant synchrotron emission from the region and is thus consistent with very faint or non-polarized emission associated with AME.
The gamma-ray binary LS I +61 303 is composed of a Be star and a compact companion orbiting in an eccentric orbit. Variable flux modulated with the orbital period of ~26.5 d has been detected from radio to very high-energy gamma rays. In addition, the system presents a superorbital variability of the phase and amplitude of the radio outburst with a period of ~4.6 yr. We present optical photometric observations of LS I +61 303 spanning ~1.5 yr and contemporaneous Halpha equivalent width (EW Halpha) data. The optical photometry shows, for the first time, that the known orbital modulation suffers a positive orbital phase shift and an increase in flux for data obtained 1-yr apart. This behavior is similar to that already known at radio wavelengths, indicating that the optical flux follows the superorbital variability as well. The orbital modulation of the EW Halpha presents the already known superorbital flux variability but shows, also for the first time, a positive orbital phase shift. In addition, the optical photometry exhibits a lag of ~0.1-0.2 in orbital phase with respect to the EW Halpha measurements at similar superorbital phases, and presents a lag of ~0.1 and ~0.3 orbital phases with respect noncontemperaneous radio and X-ray outbursts, respectively. The phase shifts detected in the orbital modulation of thermal indicators, such as the optical flux and the EW Halpha, are in line with the observed behavior for nonthermal indicators, such as X-ray or radio emission. This shows that there is a strong coupling between the thermal and nonthermal emission processes in the gamma-ray binary LS I +61 303. The orbital phase lag between the optical flux and the EW Halpha is naturally explained considering different emitting regions in the circumstellar disk, whereas the secular evolution might be caused by the presence of a moving one-armed spiral density wave in the disk.
We explore the feasibility of estimating primary cosmic ray composition at high energies from the study of two parameters of Extensive Air Showers (EAS) at ground and underground level with Monte Carlo simulations using the new EPOS and QGSJETII hadronic models tuned with LHC data. Namely, the slope and density at a given distance of the muon lateral distribution function are analysed in this work. The power to discriminate primary masses is quantified in terms of merit factor for each parameter. The analysis considers three different primary particles (proton, iron and gamma), four different zenith angles (0$^{\circ}$, 15$^{\circ}$, 30$^{\circ}$ and 45$^{\circ}$) and primary energies of $10^{17.25}$ eV, $10^{17.50}$ eV and $10^{17.75}$ eV.
The light curves observed in microlensing events due to binary lenses span an extremely wide variety of forms, characterised by U-shaped caustic crossings and/or additional smoother peaks. However, all peaks of the binary-lens light curve can be traced back to features of caustics of the lens system. Moreover, all peaks can be categorised as one of only four types (cusp-grazing, cusp-crossing, fold-crossing or fold-grazing). This enables us to present the first complete map of the parameter space of the equal-mass case by identifying regions in which light curves feature the same number and nature of peaks. We find that the total number of morphologies that can be obtained is 73 out of 232 different regions. The partition of the parameter space so-obtained provides a new key to optimise modelling of observed events through a clever choice of initial conditions for fitting algorithms.
The polarization of an X-ray beam that produces electrons with velocity components perpendicular to the beam generates an azimuthal distribution of the ejected electrons. We present methods for simulating and for analyzing the angular dependence of electron detections which enable us to derive simple analytical expressions for useful statistical properties of observable data. The derivations are verified by simulations. While we confirm the results of previous work on this topic, we provide an extension needed for analytical treatment of the full range of possible polarization amplitudes.
We extend models of our Galaxy based on distribution functions (DFs) that are
analytic functions of the action integrals to extended distribution functions
(EDFs), which have an analytic dependence on metallicity as well. We use a
simple, but physically-motivated, functional forms for the metallicity of the
interstellar medium as a function of radius and time and for the star-formation
rate, and a model for the diffusion of stars through phase space to suggest the
required functional form of an EDF. We introduce a simple prescription for
radial migration that preserves the overall profile of the disc while allowing
individual stars to migrate throughout the disc. Our models explicitly consider
the thin and thick discs as two distinct components separated in age.
We show how an EDF can be used to incorporate realistic selection functions
in models, and to construct mock catalogues of observed samples. We show that
the selection function of the Geneva-Copenhagen Survey (GCS) biases in favour
of young stars, which have atypically small random velocities. With the
selection function taken into account our models produce good fits of the GCS
data in chemo-dynamical space and the Gilmore and Reid (1983) density data.
From our EDF, we predict the structure of the SEGUE G-dwarf sample. The
kinematics are successfully predicted. The predicted metallicity distribution
has too few stars with [Fe/H]$\simeq-0.5$ dex and too many metal-rich stars. A
significant problem may be the lack of any chemical-kinematic correlations in
our thick disc. We argue that EDFs will prove essential tools for the analysis
of both observational data and sophisticated models of Galaxy formation and
evolution.
We present a general expression for the values of the average kinetic energy and of the temperature of kinetic decoupling of a WIMP, valid for any cosmological model. We show an example of the usage of our solution when the Hubble rate has a power-law dependence on temperature, and we show results for the specific cases of kination cosmology and low- temperature reheating cosmology.
The interstellar dust content in galaxies can be traced in extinction at optical wavelengths, or in emission in the far-infrared. Several studies have found that radiative transfer models that successfully explain the optical extinction in edge-on spiral galaxies generally underestimate the observed FIR/submm fluxes by a factor of about three. In order to investigate this so-called dust energy balance problem, we use two Milky Way-like galaxies produced by high-resolution hydrodynamical simulations. We create mock optical edge-on views of these simulated galaxies (using the radiative transfer code SKIRT), and we then fit the parameters of a basic spiral galaxy model to these images (using the fitting code FitSKIRT). The basic model includes smooth axisymmetric distributions along a S\'ersic bulge and exponential disc for the stars, and a second exponential disc for the dust. We find that the dust mass recovered by the fitted models is about three times smaller than the known dust mass of the hydrodynamical input models. This factor is in agreement with previous energy balance studies of real edge-on spiral galaxies. On the other hand, fitting the same basic model to less complex input models (e.g. a smooth exponential disc with a spiral perturbation or with random clumps), does recover the dust mass of the input model almost perfectly. Thus it seems that the complex asymmetries and the inhomogeneous structure of real and hydrodynamically simulated galaxies are a lot more efficient at hiding dust than the rather contrived geometries in typical quasi-analytical models. This effect may help explain the discrepancy between the dust emission predicted by radiative transfer models and the observed emission in energy balance studies for edge-on spiral galaxies.
Tracking the spectral evolution of transiently accreting neutron stars between outburst and quiescence probes relatively poorly understood accretion regimes. Such studies are challenging because they require frequent monitoring of sources with luminosities below the thresholds of current all-sky X-ray monitors. We present the analysis of over 30 observations of the neutron star low-mass X-ray binary SAX J1750.8-2900 taken across four years with the X-ray telescope aboard Swift. We find spectral softening with decreasing luminosity both on long ($\sim$1 year) and short ($\sim$days to week) timescales. As the luminosity decreases from $4\times10^{36}$ erg s$^{-1}$ to $ \sim1\times10^{35} $ erg s$^{-1}$ (0.5-10 keV), the power law photon index increases from from 1.4 to 2.9. Although not statistically required, our spectral fits allow an additional soft component that displays a decreasing temperature as the luminosity decreases from $4 \times 10^{36} $ to $6 \times 10^{34}$ erg s$^{-1}$. Spectral softening exhibited by SAX J1750.8-2900 is consistent both with accretion emission whose spectral shape steepens with decreasing luminosity and also with being dominated by a changing soft component, possibly associated with accretion onto the neutron star surface, as the luminosity declines.
We present a pilot study on the origin and assembly history of the ICL for four galaxy clusters at 0.44<z<0.57 observed with the Hubble Space Telescope from the Cluster Lensing and Supernova Survey with Hubble (CLASH) sample. Using this sample of clusters we set an empirical limit on the amount of scatter in ICL surface brightness profiles of such clusters at z=0.5 and constrain the progenitor population and formation mechanism of the ICL by measuring the ICL surface brightness profile, the ICL color and color gradient, and the total ICL luminosity within 10<r<110 kpc. The observed scatter is physical, which we associate with differences in ICL assembly process, formation epoch, and/or ICL content. Using stellar population synthesis models we transform the observed colors to metallicity. For three of the four clusters we find clear negative gradients that, on average, decrease from super solar in the central regions of the BCG to sub-solar in the ICL. Such negative color/metallicity gradients can arise from tidal stripping of L* galaxies and/or the disruption of dwarf galaxies, but not major mergers with the BCG. We also find that the ICL at 110 kpc has a color comparable to m*+2 red sequence galaxies and a total luminosity between 10<r<110 kpc of 4-8 L*. This suggests that the ICL is dominated by stars liberated from galaxies with L>0.2 L* and that neither dwarf disruption nor major mergers with the BCG alone can explain the observed level of luminosity and remain consistent with either the observed evolution in the faint end slope of the luminosity function or predictions for the number of BCG major mergers since z=1. Taken together, the results of this pilot study are suggestive of a formation history for these clusters in which the ICL is built-up by the stripping of >0.2 L* galaxies, and disfavor significant contribution to the ICL by dwarf disruption or major mergers with the BCG.
Recent studies of dust in the interstellar medium have challenged the capabilities and validity of current dust models, indicating that the properties of dust evolve as it transits between different phases of the interstellar medium. We conduct a multi-wavelength study of the dust emission from the ionized gas of the IC 434 emission nebula, and combine this with modeling, from large scales that provide insight into the history of the IC 434/L1630 region, to small scales that allow us to infer quantitative properties of the dust content inside the H II region. The dust enters the H II region through momentum transfer with a champagne flow of ionized gas, set up by a chance encounter between the L1630 molecular cloud and the star cluster of $\sigma$ Ori. We observe two clearly separated dust populations inside the ionized gas, that show different observational properties, as well as contrasting optical properties. Population A is colder ($\sim$ 25 K) than predicted by widely-used dust models, its temperature is insensitive to an increase of the impinging radiation field, is momentum-coupled to the gas, and efficiently absorbs radiation pressure to form a dust wave at 1.0 pc ahead of $\sigma$ Ori AB. Population B is characterized by a constant [20/30] flux ratio throughout the HII region, heats up to $\sim$ 75 K close to the star, and is less efficient in absorbing radiation pressure, forming a dust wave at 0.1 pc from the star. We conclude that the dust inside IC 434 is bimodal. The characteristics of population A are remarkable and can not be explained by current dust models. Population B are grains that match the classical description of spherical, compact dust. Our results confirm recent work that stress the importance of variations in the dust properties between different regions of the interstellar medium.
The symmetric peak observed in linear polarization in the core of the solar sodium D$_1$ line at 5896 \AA\ has remained enigmatic since its discovery nearly two decades ago. One reason is that the theory of polarized scattering has not been experimentally tested for multi-level atomic systems in the relevant parameter domains, although the theory is continually being used for the interpretation of astrophysical observations. A laboratory experiment that was set up a decade ago to find out whether the D$_1$ enigma is a problem of solar physics or quantum physics revealed that the D$_1$ system has a rich polarization structure in situations where standard scattering theory predicts zero polarization, even when optical pumping of the $m$ state populations of the hyperfine-split ground state is accounted for. Here we show that the laboratory results can be modeled in great quantitative detail if the theory is extended to include the coherences in both the initial and final states of the scattering process. Radiative couplings between the allowed dipole transitions generate coherences in the initial state. Corresponding coherences in the final state are then demanded by a phase closure selection rule. The experimental results for the well understood D$_2$ line are used to constrain the two free parameters of the experiment, collision rate and optical depth, to suppress the need for free parameters when fitting the D$_1$ results.
Light degrees of freedom that modify gravity on cosmological scales must be "screened" on solar system scales in order to be compatible with data. The Vainshtein mechanism achieves this through a breakdown of classical perturbation theory, as large interactions involving new degrees of freedom become important below the so-called Vainshtein radius. We begin to develop an extension of the Parameterized Post-Newtonian (PPN) formalism that is able to handle Vainshteinian corrections. We argue that theories with a unique Vainshtein scale must be expanded using two small parameters. In this Parameterized Post-Newtonian-Vainshteinian (PPNV) expansion, the primary expansion parameter that controls the PPN order is, as usual, the velocity $v$. The secondary expansion parameter, $\alpha$, controls the strength of the Vainshteinian correction and is a theory-specific combination of the Schwarzschild radius and the Vainshtein radius of the source that is independent of its mass. We present the general framework and apply it to Cubic Galileon theory both inside and outside the Vainshtein radius. The PPNV framework can be used to determine the compatibility of such theories with solar system and other strong-field data.
We investigate the cosmological evolution of mimetic matter model with arbitrary scalar potential. The cosmological reconstruction is explicitly done for different choices of potential. The cases that mimetic matter model shows the evolution as Cold Dark Matter(CDM), wCDM model, dark matter and dark energy with dynamical $Om(z)$ or phantom dark energy with phantom-non-phantom crossing are presented in detail. The cosmological perturbations for such evolution are studied in mimetic matter model. For instance, the evolution behavior of the matter density contrast which is different from usual one, i.e. $\ddot \delta + 2 H \dot \delta - \kappa ^2 \rho \delta /2 = 0$ is investigated. The possibility of peculiar evolution of $\delta$ in the model under consideration is shown. Special attention is paid to the behavior of matter density contrast near to future singularity where decay of perturbations may occur much earlier the singularity.
Links to: arXiv, form interface, find, astro-ph, recent, 1501, contact, help (Access key information)
Galactic feedback strongly affects the way galactic environments are enriched. We examine this connection by performing a suite of cosmological hydrodynamic simulations, exploring a range of parameters based on the galaxy formation model developed in Vogelsberger et al. 2013 (henceforth V13). We examine the effects of AGN feedback, wind mass loading, wind specific energy, and wind metal-loading on the properties of the circumgalactic medium (CGM) of galaxies with $M_\text{halo} > 10^{11} M_\odot$. Note that while the V13 model was tuned to match observations including the stellar mass function, no explicit tuning was done for the CGM. The wind energy per unit outflow mass has the most significant effect on the CGM enrichment. High energy winds launch metals far beyond the virial radius. AGN feedback also has a significant effect, but only at $z < 3$. We compare to high redshift HI and CIV observations. All our simulations produce the observed number of Damped Lyman-$\alpha$ Absorbers. At lower column density, several of our simulations produce enough Lyman Limit Systems (LLS) $100$ kpc from the galaxy, but in all cases the LLS abundance falls off with distance faster than observations, with too few LLS at $200$ kpc. Further, in all models the CIV abundance drops off too sharply with distance, with too little CIV $100$-$200$ kpc from the galaxy. Higher energy wind models produce more extended CIV but also produce less stars, in tension with star-formation rate density observations. This highlights the fact that circumgalactic observations are a strong constraint on galactic feedback models.
More than 50 years after the dawn of X-ray astronomy, the dynamical parameters of the prototypical X-ray binary Sco X-1 are still unknown. We combine a Monte Carlo analysis, which includes all the previously known orbital parameters of the system, along with the K-correction to set dynamical constraints to the masses of the compact object (M1<1.73 Msun) and the companion star (0.28 Msun<M2<0.70 Msun). For the case of a canonical neutron star mass of M1= 1.4 Msun, the orbital inclination is found to be lower than 40 degree. We also present the best near-infrared spectrum of the source to date. There is no evidence of donor star features on it, but we are able to constrain the veiling factor as a function of the spectral type of the secondary star. The combination of both techniques restricts the spectral type of the donor to be later than K4 and luminosity class IV. It also constrains the contribution of the companion light to the infrared emission of Sco X-1 to be lower than 33 percent. This implies that the accretion related luminosity of the system in the K band is larger than 4 x 10E35 erg/s.
The CHemical Abundances of Spirals (CHAOS) project leverages the combined power of the Large Binocular Telescope (LBT) with the broad spectral range and sensitivity of the Multi Object Double Spectrograph (MODS) to measure "direct" abundances in large samples of HII regions in spiral galaxies. We present LBT MODS observations of 62 HII regions in the nearby NGC628. We measure one or more auroral lines ([OIII] 4363, [NII] 5755, [SIII] 6312, or [OII] 7320,7330) in a large number of HII regions (40). Comparing derived temperatures from multiple auroral line measurements, we find: a strong correlation between temperatures based on [SIII] and [NII]; and large discrepancies for some temperatures based on [OII] and [OIII]. These trends are consistent with other observations in the literature, yet, given the widespread use and acceptance of [OIII] as a temperature determinant, the magnitude of the T[OIII] discrepancies still came as a surprise. Based on these results, we conduct a uniform abundance analysis using the temperatures derived from [SIII] and [NII], and report the gas-phase abundance gradients for NGC628. Relative abundances of S/O, Ne/O, and Ar/O are constant across the galaxy, consistent with no systematic change in the upper IMF over the sampled range in metallicity. These alpha-element ratios, along with N/O, all show small dispersions consistent with no intrinsic dispersion (0.05) over 70% of the azimuthally averaged radius. We interpret these results as an indication that, at a given radius, the interstellar medium in NGC628 is chemically well-mixed. Unlike the nearly temperature-independent gradients, O/H abundances have a larger intrinsic dispersion of ~0.13 dex. We posit that this dispersion represents an upper limit to the true dispersion in O/H at a given radius and that some of that dispersion is due to systematic uncertainties arising from temperature measurements.
We have observed NGC5194 (M51a) as part of the CHemical Abundances of Spirals (CHAOS) project. Using the Multi Object Double Spectrographs (MODS) on the Large Binocular Telescope (LBT) we are able to measure one or more of the temperature-sensitive auroral lines ([O III] 4363, [N II] 5755, [S III] 6312) and thus measure "direct" gas-phase abundances in 29 individual HII regions. [O III] 4363 is only detected in two HII regions both of which show indications of excitation by shocks. We compare our data to previous direct abundances measured in NGC5194 and find excellent agreement for all but one region (Delta[log(O/H)] ~ 0.04). We find no evidence of trends in Ar/O, Ne/O, or S/O within NGC5194 or compared to other galaxies. We find modest negative gradients in both O/H and N/O with very little scatter (sigma < 0.08 dex), most of which can be attributed to random error and not to intrinsic dispersion. The gas-phase abundance gradient is consistent with the gradients observed in other interacting galaxies, which tend to be shallower than gradients measured in isolated galaxies. The N/O ratio (<log(N/O)> = -0.62) suggests secondary nitrogen production is responsible for a significantly larger fraction of nitrogen (e.g., factor of 8-10) relative to primary production mechanisms than predicted by theoretical models.
Recent discoveries have put the picture of stellar clusters being simple stellar populations into question. In particular, the color-magnitude diagrams of intermediate age (1-2 Gyr) massive clusters in the Large Magellanic Cloud (LMC) show features that could be interpreted as age spreads of 100-500 Myr. If multiple generations of stars are present in these clusters then, as a consequence, young (<1 Gyr) clusters with similar properties should have age spreads of the same order. In this paper we use archival Hubble Space Telescope (HST) data of eight young massive LMC clusters (NGC 1831, NGC 1847, NGC 1850, NGC 2004, NGC 2100, NGC 2136, NGC 2157 and NGC 2249) to test this hypothesis. We analyzed the color-magnitude diagrams of these clusters and fitted their star formation history to derive upper limits of potential age spreads. We find that none of the clusters analyzed in this work shows evidence for an extended star formation history that would be consistent with the age spreads proposed for intermediate age LMC clusters. Tests with artificial single age clusters show that the fitted age dispersion of the youngest clusters is consistent with spreads that are purely induced by photometric errors. As an additional result we determined a new age of NGC 1850 of ~100 Myr, significantly higher than the commonly used value of about 30 Myr, although consistent with early HST estimates.
[ABRIDGED] $\omega$ Centauri (NGC 5139) contains large numbers of variable
stars of different types and, in particular, more than a hundred RR Lyrae
stars. We have conducted a variability survey of $\omega$ Cen in the NIR, using
ESO's 4.1m Visible and Infrared Survey Telescope for Astronomy (VISTA). This is
the first paper of a series describing our results.
$\omega$ Cen was observed using VIRCAM mounted on VISTA. A total of 42 and
100 epochs in $J$ and $K_{\rm S}$, respectively, were obtained, distributed
over a total timespan of 352 days. PSF photometry was performed, and periods of
the known variable stars were improved when necessary using an ANOVA analysis.
An unprecedented homogeneous and complete NIR catalogue of RR Lyrae stars in
the field of $\omega$ Cen was collected, allowing us to study, for the first
time, all the RR Lyrae stars associated to the cluster, except 4 located far
away from the cluster center. Membership status, subclassifications between
RRab and RRc subtypes, periods, amplitudes, and mean magnitudes were derived
for all the stars in our sample. Additionally, 4 new RR Lyrae stars were
discovered, 2 of them with high probability of being cluster members. The
distribution of $\omega$ Cen stars in the Bailey (period-amplitude) diagram is
also discussed. Reference lines in this plane, for both Oosterhoff type I (OoI)
and II (OoII) components, are provided.
In the present paper, we clarify the status of many (candidate) RR Lyrae
stars that had been unclear in previous studies. This includes stars with
anomalous positions in the color-magnitude diagram, uncertain periods or/and
variability types, and possible field interlopers. We conclude that $\omega$
Cen hosts a total of 88 RRab and 101 RRc stars, for a grand total of 189 likely
members. We confirm that most RRab stars in the cluster belong to an OoII
component, as previously found using visual data.
SuperSpec is a novel on-chip spectrometer we are developing for multi-object, moderate resolution (R = 100 - 500), large bandwidth (~1.65:1) submillimeter and millimeter survey spectroscopy of high-redshift galaxies. The spectrometer employs a filter bank architecture, and consists of a series of half-wave resonators formed by lithographically-patterned superconducting transmission lines. The signal power admitted by each resonator is detected by a lumped element titanium nitride (TiN) kinetic inductance detector (KID) operating at 100-200 MHz. We have tested a new prototype device that is more sensitive than previous devices, and easier to fabricate. We present a characterization of a representative R=282 channel at f = 236 GHz, including measurements of the spectrometer detection efficiency, the detector responsivity over a large range of optical loading, and the full system optical efficiency. We outline future improvements to the current system that we expect will enable construction of a photon-noise-limited R=100 filter bank, appropriate for a line intensity mapping experiment targeting the [CII] 158 micron transition during the Epoch of Reionization
Deep surveys with the SKA1-MID array offer for the first time the opportunity to systematically explore the polarization properties of the microJy source population. Our knowledge of the polarized sky approaching these levels is still very limited. In total intensity the population will be dominated by star-forming and normal galaxies to intermediate redshifts ($z \sim1-2$), and low-luminosity AGN to high redshift. The polarized emission from these objects is a powerful probe of their intrinsic magnetic fields and of their magnetic environments. For redshift of order 1 and above the broad bandwidth of the mid-bands span the Faraday thick and thin regimes allowing study of the intrinsic polarization properties of these objects as well as depolarization from embedded and foreground plasmas. The deep field polarization images will provide Rotation Measures data with very high solid angle density allowing a sensitive statistical analysis of the angular variation of RM on critical arc-minute scales from a magnetic component of Large Scale Structure of the Universe.
The atmospheres of exoplanets are commonly studied by observing the transit of the planet passing in front of its parent star. The obscuration of part of the stellar disk during a transit will reveal aspects of its surface structure resulting from general centre-to-limb variations (CLVs). These become apparent when forming the ratio between the stellar light in and out of transit. These phenomena can be seen particularly clearly during the progress of a penumbral lunar eclipse, where the Earth transits the solar disk and masks different regions of the solar disk as the eclipse progresses. When inferring the properties of the planetary atmosphere, it is essential that this effect originating at the star is properly accounted for. Using the data observed from the 2014-April-15 lunar eclipse with the ESPaDOnS spectrograph mounted on the Canada France Hawaii Telescope (CFHT), we have obtained for the first time a time sequence of the penumbral spectra. These penumbral spectra enable us to study the centre-to-limb variations of solar Fraunhofer lines when the Earth is transiting Sun. The Na i and Ca ii absorption features reported from previous lunar eclipse observations are demonstrated to be CLV features, which dominate the corresponding line profiles and mask possible planetary signal. Detecting atmospheric species in exoplanets via transit spectroscopy must account for the CLV effect.
We investigate the possibility of detecting the 3D cross correlation power spectrum of the Ly-$\alpha$ forest and HI 21 cm signal from the post reionization epoch. Other than the direct dependence on the dark matter power spectrum, the cross-correlation signal is found to be sensitive to the respective bias parameters which dictates the strength of anisotropy in redshift space. We find that the cross-correlation power spectrum can be detected using $400 ~ \, \rm hrs$ observation with SKA-mid (phase 1) and a futuristic BOSS like experiment with a quasar density of $30 ~ \rm deg^{-2}$ at a peak SNR of $15$ for a single field experiment at redshift $z = 2.5$. We model the clustering of HI distribution for the Ly-$\alpha$ forest and 21 cm signal on large scales using the linear bias model. We probe the possibility of independently constraining these parameters using the cross power spectrum. We find that with the same experiments $1 \sigma$ marginalized errors on the 21-cm linear redshift space distortion parameter $\beta_T$ and $\beta_F$ corresponding to the Ly-$\alpha $ forest are $\sim 2.7 \%$ and $\sim 1.4 \%$ respectively for $10$ independent pointings of the SKA. This prediction indicates a significant improvement over existing measurements. We claim that the cross correlation of the Ly-$\alpha$ and 21 cm observation in 3D not only ascertains the cosmological origin of the signal in presence of astrophysical foregrounds but also provides stringent constraints on large scale HI bias. This provides an independent probe towards understanding cosmological structure formation.
We investigate the gravitational wave (GW) signal generated by a population of double neutron-star binaries (DNS) with eccentric orbits caused by kicks during supernova collapse and binary evolution. The DNS population of a standard Milky-Way type galaxy has been studied as a function of star formation history, initial mass function (IMF) and metallicity and of the binary-star common-envelope ejection process. The model provides birth rates, merger rates and total numbers of DNS as a function of time. The GW signal produced by this population has been computed and expressed in terms of a hypothetical space GW detector (eLISA) by calculating the number of discrete GW signals at different confidence levels, where `signal' refers to detectable GW strain in a given frequency-resolution element. In terms of the parameter space explored, the number of DNS-originating GW signals is greatest in regions of recent star formation, and is significantly increased if metallicity is reduced from 0.02 to 0.001, consistent with Belczynski10a. Increasing the IMF power-law index (from --2.5 to --1.5) increases the number of GW signals by a large factor. This number is also much higher for models where the common-envelope ejection is treated using the $\alpha-$mechanism (energy conservation) than when using the $\gamma-$mechanism (angular-momentum conservation). We have estimated the total number of detectable DNS GW signals from the Galaxy by combining contributions from thin disc, thick disc, bulge and halo. The most probable numbers for an eLISA-type experiment are 0-1600 signals per year at S/N$\geqslant$1, 0-900 signals per year at S/N$\geqslant$3, and 0-570 at S/N$\geqslant$5, coming from about 0-65, 0-60 and 0-50 resolved DNS respectively.
We report the discovery of the merging cluster, RXCJ2359.3-6042, from the REFLEX II cluster survey and present our results from all three detectors combined in the imaging and spectral analysis of the XMM-Newton data. Also known as Abell 4067, this is a unique system, where a compact bullet penetrates an extended, low density cluster at redshift z=0.099 clearly seen from our follow-up XMM-Newton observation. The bullet goes right through the central region of the cluster without being disrupted and we can clearly watch the process how the bullet component is stripped of its layers outside the core. There is an indication of a shock heated region in the East of the cluster with a higher temperature. The bulk temperature of the cluster is about 3.12 keV implying a lower mass system. Spearheading the bullet is a cool core centred by a massive early type galaxy. The temperatures and metallicities of a few regions in the cluster derived from the spectral analysis supports our conjecture based on the surface brightness image that a much colder compact component at 1.55 keV with large metallicity (0.75 Zsol) penetrates the main cluster, where the core of the infalling component survived the merger leaving stripped gas behind at the centre of the main cluster. We also give an estimate of the total mass within r500, which is about 2e14Msol from the deprojected spherical-beta modelling of the cluster in good agreement with other mass estimates from the M--Tx and M-sigma_v relations.
We measured the sizes and morphological parameters of LARS galaxies in the continuum, Lya, and Ha images. We studied morphology by using the Gini coefficient vs M20 and asymmetry vs concentration diagrams. We then simulated LARS galaxies at z~2 and 5.7, performing the same morphological measurements. We also investigated the detectability of LARS galaxies in current deep field observations. The subsample of LAEs within LARS (LARS-LAEs) was stacked to provide a comparison to stacking studies performed at high redshift. LARS galaxies have continuum size, stellar mass, and rest-frame absolute magnitude typical of Lyman break analogues in the local Universe and also similar to 2<z<3 star-forming galaxies and massive LAEs. LARS optical morphology is consistent with the one of merging systems, and irregular or starburst galaxies. For the first time we quantify the morphology in Lya images: even if a variety of intrinsic conditions of the interstellar medium can favour the escape of Lya photons, LARS-LAEs appear small in the continuum, and their Lya is compact. LARS galaxies tend to be more extended in Lya than in the rest-frame UV. It means that Lya photons escape by forming haloes around HII regions of LARS galaxies. The stack of LARS-LAE Lya images is peaked in the centre, indicating that the conditions, which make a galaxy an LAE, tend to produce a concentrated surface brightness profile. On the other hand, the stack of all LARS galaxies is shallower and more extended. This can be caused by the variety of dust and HI amount and distribution, which produces a more complex, patchy, and extended profile, like the one observed for Lyman break galaxies that can contribute to the stack. We cannot identify a single morphological property that controls whether a galaxy emits a net positive Lya flux. However, the LARS-LAEs have continuum properties consistent with merging systems.
Dark Stars (DS) are stellar objects made (almost entirely) of ordinary atomic material but powered by the heat from Dark Matter (DM) annihilation (rather than by fusion). Weakly Interacting Massive Particles (WIMPs), among the best candidates for DM, can be their own antimatter and can accumulate inside the star, with their annihilation products thermalizing with and heating the DS. The resulting DSs are in hydrostatic and thermal equilibrium. The first phase of stellar evolution in the history of the Universe may have been dark stars. Though DM constituted only $<0.1\%$ of the mass of the star, this amount was sufficient to power the star for millions to billions of years. Depending on their DM environment, early DSs can become very massive ($>10^6 M_\odot$), very bright ($>10^9 L_\odot$), and potentially detectable with the James Webb Space Telescope (JWST). Once the DM runs out and the dark star dies, it may collapse to a black hole; thus DSs can provide seeds for the supermassive black holes observed throughout the Universe and at early times. Other sites for dark star formation exist in the Universe today in regions of high dark matter density such as the centers of galaxies. The current review briefly discusses DSs existing today but focuses on the early generation of dark stars.
We present the observational results of the Gamma-ray blazar, 3C 66A, at 2.3, 8.4, and 22 GHz at 4 epochs during 2004-05 with the VLBA. The resulting images show an overall core-jet structure extending roughly to the south with two intermediate breaks occurring in the region near the core. By model-fitting to the visibility data, the northmost component, which is also the brightest, is identified as the core according to its relatively flat spectrum and its compactness. As combined with some previous results to investigate the proper motions of the jet components, it is found the kinematics of 3C 66A is quite complicated with components of inward and outward, subluminal and superluminal motions all detected in the radio structure. The superluminal motions indicate strong Doppler boosting exists in the jet. The apparent inward motions of the innermost components last for at least 10 years and could not be caused by new-born components. The possible reason could be non-stationarity of the core due to opacity change.
We present the calibration of the Swift UVOT grisms, of which there are two, providing low-resolution field spectroscopy in the ultraviolet and optical bands respectively. The UV grism covers the range 1700-5000 Angstrom with a spectral resolution of 75 at 2600 Angstrom for source magnitudes of u=10-16 mag, while the visible grism covers the range 2850-6600 Angstrom with a spectral resolution of 100 at 4000 Angstrom for source magnitudes of b=12-17 mag. This calibration extends over all detector positions, for all modes used during operations. The wavelength accuracy (1-sigma) is 9 Angstrom in the UV grism clocked mode, 17 Angstrom in the UV grism nominal mode and 22 Angstrom in the visible grism. The range below 2740 Angstrom in the UV grism and 5200 Angstrom in the visible grism never suffers from overlapping by higher spectral orders. The flux calibration of the grisms includes a correction we developed for coincidence loss in the detector. The error in the coincidence loss correction is less than 20%. The position of the spectrum on the detector only affects the effective area (sensitivity) by a few percent in the nominal modes, but varies substantially in the clocked modes. The error in the effective area is from 9% in the UV grism clocked mode to 15% in the visible grism clocked mode .
The ultra-long period Cepheids (ULPCs) are classical Cepheids with pulsation periods exceeding $\approx 80$ days. The intrinsic brightness of ULPCs are ~1 to ~3 mag brighter than their shorter period counterparts. This makes them attractive in future distance scale work to derive distances beyond the limit set by the shorter period Cepheids. We have initiated a program to search for ULPCs in M31, using the single-band data taken from the Palomar Transient Factory, and identified eight possible candidates. In this work, we presented the VI-band follow-up observations of these eight candidates. Based on our VI-band light curves of these candidates and their locations in the color-magnitude diagram and the Period-Wesenheit diagram, we verify two candidates as being truly ULPCs. The six other candidates are most likely other kinds of long-period variables. With the two confirmed M31 ULPCs, we tested the applicability of ULPCs in distance scale work by deriving the distance modulus of M31. It was found to be $\mu_{M31,ULPC}=24.30\pm0.76$ mag. The large error in the derived distance modulus, together with the large intrinsic dispersion of the Period-Wesenheit (PW) relation and the small number of ULPCs in a given host galaxy, means that the question of the suitability of ULPCs as standard candles is still open. Further work is needed to enlarge the sample of calibrating ULPCs and reduce the intrinsic dispersion of the PW relation before re-considering ULPCs as suitable distance indicators.
Hot molecular cores (HMCs) are intermediate stages of high-mass star formation and are also known for their rich emission line spectra at (sub-)mm wavebands. The observed spectral feature of HMCs such as total number of emission lines and associated line intensities are also found to vary with evolutionary stages. We developed various 3D models for HMCs guided by the evolutionary scenarios proposed by recent empirical and modeling studies. We then investigated the spatio-temporal variation of temperature and molecular abundances in HMCs by consistently coupling gas-grain chemical evolution with radiative transfer calculations. We explored the effects of varying physical conditions on molecular abundances including density distribution and luminosity evolution of the central protostar(s). The time-dependent temperature structure of the hot core models provides a realistic framework for investigating the spatial variation of ice mantle evaporation as a function of evolutionary timescales. With increasing protostellar luminosity, the water ice evaporation font ($\sim$100K) expands and the spatial distribution of gas phase abundances of these COMs also spreads out. We simulated the synthetic spectra for these models at different evolutionary timescales to compare with observations. A qualitative comparison of the simulated and observed spectra suggests that these self-consistent hot core models can reproduce the notable trends in hot core spectral variation within the typical hot core timescales of 10$^{5}$ year. These models predict that the spatial distribution of various emission line maps will also expand with evolutionary time. The model predictions can be compared with high resolution observation that can probe scales of a few thousand AU in high-mass star forming regions such as from ALMA.[Abridged]
We present Wide Field Spectrograph (WiFeS) integral field spectroscopy and HST FOS spectroscopy for the LINER galaxy NGC 1052. We infer the presence of a turbulent accretion flow forming a small-scale accretion disk. We find a large-scale outflow and ionisation cone along the minor axis of the galaxy. Part of this outflow region is photoionised by the AGN, and shares properties with the ENLR of Seyfert galaxies, but the inner ($R \lesssim 1.0$~arcsec) accretion disk and the region around the radio jet appear shock excited. The emission line properties can be modelled by a "double shock" model in which the accretion flow first passes through an accretion shock in the presence of a hard X-ray radiation, and the accretion disk is then processed through a cocoon shock driven by the overpressure of the radio jets. This model explains the observation of two distinct densities ($\sim10^4$ and $\sim10^6$ cm$^{-3}$), and provides a good fit to the observed emission line spectrum. We derive estimates for the velocities of the two shock components and their mixing fractions, the black hole mass, the accretion rate needed to sustain the LINER emission and derive an estimate for the jet power. Our emission line model is remarkably robust against variation of input parameters, and so offers a generic explanation for the excitation of LINER galaxies, including those of spiral type such as NGC 3031 (M81).
Solid material in protoplanetary discs will suffer one of two fates after the epoch of planet formation; either being bound up into planetary bodies, or remaining in smaller planetesimals to be ground into dust. These end states are identified through detection of sub-stellar companions by periodic radial velocity (or transit) variations of the star, and excess emission at mid- and far-infrared wavelengths, respectively. Since the material that goes into producing the observable outcomes of planet formation is the same, we might expect these components to be related both to each other and their host star. Heretofore, our knowledge of planetary systems around other stars has been strongly limited by instrumental sensitivity. In this work, we combine observations at far-infrared wavelengths by IRAS, Spitzer, and Herschel with limits on planetary companions derived from non-detections in the 16-year Anglo-Australian Planet Search to clarify the architectures of these (potential) planetary systems and search for evidence of correlations between their constituent parts. We find no convincing evidence of such correlations, possibly owing to the dynamical history of the disk systems, or the greater distance of the planet-search targets. Our results place robust limits on the presence of Jupiter analogs which, in concert with the debris disk observations, provides insights on the small-body dynamics of these nearby systems.
We describe a successful effort to produce a laser comb around 1.55 $\mu$m in the astronomical H band using a method based on a line-referenced, electro-optical-modulation frequency comb. We discuss the experimental setup, laboratory results, and proof of concept demonstrations at the NASA Infrared Telescope Facility (IRTF) and the Keck-II telescope. The laser comb has a demonstrated stability of $<$ 200 kHz, corresponding to a Doppler precision of ~0.3 m/s. This technology, when coupled with a high spectral resolution spectrograph, offers the promise of $<$1 m/s radial velocity precision suitable for the detection of Earth-sized planets in the habitable zones of cool M-type stars.
Based on high-spectral resolution observations performed with the echelle spectrograph NES of the 6-meter telescope, we have studied the peculiarities of the spectrum and the velocity field in the atmosphere and envelope of the optical counterpart of the infrared source IRAS 23304+6347. We have concluded about the absence of significant variations in the radial velocity Vr inferred from atmospheric absorptions and about its coincidence with the systemic velocity deduced from radio data. The envelope expansion velocity Vexp=15.5 km/s has been determined from the positions of rotational band lines of the C$_2$ Swan (0; 0) band. A complex emission-absorption profile of the Swan (0; 1) 5635 \AA{} has been recorded. Analysis of the multicomponent NaI~D doublet line profile has revealed interstellar components V(IS)=$-61.6$ and $-13.2$ km/s as well as a circumstellar component with V(CS)=$-41.0$ km/s whose position corresponds to the velocity inferred from C$_2$ features. The presence of the interstellar component with Vr=$-61.6$ km/s in the spectrum allows d=2.5 kpc to be considered as a lower limit for the distance to the star. A splitting of the profiles for strong absorptions of ionized metals (YII, BaII, LaII, SiII) attributable to the presence of a short-wavelength component originating in the circumstellar envelope has been detected in the optical spectrum of IRAS23304+6347 for the first time.
The most probable initial magnetic configuration of a CME is a flux rope consisting of twisted field lines which fill the whole volume of a dark coronal cavity. The flux ropes can be in stable equilibrium in the coronal magnetic field for weeks and even months, but suddenly they loose their stability and erupt with high speed. Their transition to the unstable phase depends on the parameters of the flux rope (i.e., total electric current, twist, mass loading etc.), as well as on the properties of the ambient coronal magnetic field. One of the major governing factors is the vertical gradient of the coronal magnetic field which is estimated as decay index (n). Cold dense prominence material can be collected in the lower parts of the helical flux tubes. Filaments are therefore good tracers of the flux ropes in the corona, which become visible long before the beginning of the eruption. The perspectives of the filament eruptions and following CMEs can be estimated by the comparison of observed filament heights with calculated decay index distributions. The present paper reviews the formation of magnetic flux ropes, their stable and unstable phases, eruption conditions, and also discusses their physical implications in the solar corona.
We present the structural and star formation properties of 59 void galaxies as part of the Void Galaxy Survey (VGS). Our aim is to study in detail the physical properties of these void galaxies and study the effect of the void environment on galaxy properties. We use Spitzer 3.6 $\rm{\mu m}$ and B-band imaging to study the morphology and color of the VGS galaxies. For their star formation properties, we use Halpha and GALEX near-UV imaging. We compare our results to a range of galaxies of different morphologies in higher density environments. We find that the VGS galaxies are in general disk dominated and star forming galaxies. Their star formation rates are, however, often less than 1 $\rm{M_{\odot}}$ $\rm{yr^{-1}}$. There are two early-type galaxies in our sample as well. In $\rm{r_{e}}$ versus $\rm{M_{B}}$ parameter space, VGS galaxies occupy the same space as dwarf irregulars and spirals.
We investigate possible sky survey modes with the Middle Sized Telescopes (MST, aimed at covering the energy range from $\sim$100 GeV to 10 TeV) subsystem of the Cherenkov Telescope Array (CTA). We use the standard CTA tools, CORSIKA and sim_telarray, to simulate the development of gamma-ray showers, proton background and the telescope response. We perform simulations for the H.E.S.S.-site in Namibia, which is one of the candidate sites for the CTA experiment. We study two previously considered modes, parallel and divergent, and we propose a new, convergent mode with telescopes tilted toward the array center. For each mode we provide performance parameters crucial for choosing the most efficient survey strategy. For the non-parallel modes we study the dependence on the telescope offset angle. We show that use of both the divergent and convergent modes results in potential advantages in comparison with use of the parallel mode. The fastest source detection can be achieved in the divergent mode with larger offset angles ($\sim 6^{\circ}$ from the Field of View centre for the outermost telescopes), for which the time needed to perform a scan at a given sensitivity level is shorter by a factor of $\sim$2.3 than for the parallel mode. We note, however, the direction and energy reconstruction accuracy for the divergent mode is even by a factor of $\sim 2$ worse than for other modes. Furthermore, we find that at high energies and for observation directions close to the center of the array field of view, the best performance parameters are achieved with the convergent mode, which favors this mode for deep observations of sources with hard energy spectra.
One favored progenitor model for short duration gamma-ray bursts (SGRBs) is the coalescence of two neutron stars (NS-NS). One possible outcome of such a merger would be a rapidly spinning, strongly magnetized neutron star (known as a millisecond magnetar). These magnetars may be "supra-massive", implying they would collapse to black holes after losing centrifugal support due to magnetic dipole spindown. By systematically analyzing the BAT-XRT light curves of all short GRBs detected by {\em swift}, we test how well the data are consistent with this central engine model of short GRBs. We find that the so-called "extended emission" observed with BAT in some short GRBs are fundamentally the same component as the "internal X-ray plateau" as observed in many short GRBs, which is defined as a plateau in the lightcurve followed by a very rapid drop. Based on how likely a short GRB hosts a magnetar, we characterize the entire {\em Swift} short GRB sample into three categories: the "internal plateau" sample, the "external plateau" sample, and the "no plateau" sample. The derived magnetar surface magnetic field $B_{\rm p}$ and the initial spin period $P_0$ fall into the reasonable range. No GRBs in the internal plateau sample have the total energy exceeding the maximum energy budget of a millisecond magnetar. Assuming that the rapid fall time at the end of the internal plateau is the collapse time of a supra-massive magnetar to a black hole, and applying the measured mass distribution of NS-NS systems in our Galaxy, we constrain the neutron star equation of state (EOS). The data suggest that the NS equation of state is close to the GM1 model, which has a maximum non-rotating NS mass $M_{\rm TOV} \sim 2.37 M_\odot$.
We report observations of the flickering variability of the symbiotic recurrent nova RS Oph at quiescence in five bands (UBVRI). We find evidence of correlation between the peak-to-peak flickering amplitude (\Delta F) and the average flux of the hot component ($F_{av}$). The correlation is highly significant with correlation coefficient 0.85 and p-value $\sim 10^{-15}$. Combining the data from all wavebands, we find a dependence of the type $\Delta F \propto F_{av}^k$, with power-law index $k = 1.02 \pm 0.04$ for the UBVRI flickering of RS Oph. Thus, the rms amplitude of variability scale almost linearly with the average flux of the hot component, $< \sigma_{rms} / F_{av} > = 0.08 \pm 0.02$. The detected correlation is similar to that found in some X-ray binaries. The possible reasons are briefly discussed. The data are available upon request from the authors.
We use our vertical Milky Way disc model together with Galaxia to create mock observations of stellar samples in the solar neighbourhood. We compare these to the corresponding volume complete observational samples of dereddened and binary accounted data from Hipparcos and the Catalogue of Nearby Stars. Sampling the likelihood in the parameter space we determine a new fiducial IMF considering constraints from dwarf and giant stars. The resulting IMF observationally backed in the range from 0.5 to 10 Msun is a two slope broken power law with -1.49 +- 0.08 for the low mass slope, a break at 1.39 +- 0.05 Msun and a high mass slope of -3.02 +- 0.06. The Besancon group also converging to a similar IMF even though their observational sample being quite different to ours shows that the forward modelling technique is a powerful diagnostic to test theoretical concepts like the local field star IMF.
New measurements using radio and plasma-wave instruments in interplanetary space have shown that nanometer-scale dust, or nanodust, is a significant contributor to the total mass in interplanetary space. Better measurements of nanodust will allow us to determine where it comes from and the extent to which it interacts with the solar wind. When one of these nanodust grains impacts a spacecraft, it creates an expanding plasma cloud, which perturbs the photoelectron currents. This leads to a voltage pulse between the spacecraft body and the antenna. Nanodust has a high charge/mass ratio, and therefore can be accelerated by the interplanetary magnetic field to speeds up to the speed of the solar wind: significantly faster than the Keplerian orbital speeds of heavier dust. The amplitude of the signal induced by a dust grain grows much more strongly with speed than with mass of the dust particle. As a result, nanodust can produce a strong signal, despite their low mass. The WAVES instruments on the twin Solar TErrestrial RElations Observatory spacecraft have observed interplanetary nanodust particles since shortly after their launch in 2006. After describing a new and improved analysis of the last five years of STEREO/WAVES Low Frequency Receiver data, a statistical survey of the nanodust characteristics, namely the rise time of the pulse voltage and the flux of nanodust, is presented. Agreement with previous measurements and interplanetary dust models is shown. The temporal variations of the nanodust flux are also discussed.
The remarkable progress in cosmic microwave background (CMB) studies over past decade has led to the era of precision cosmology in striking agreement with the $\Lambda$CDM model. However, the lack of power in the CMB temperature anisotropies at large angular scales (low-$\ell$), as has been confirmed by the recent Planck data also (up to $\ell=40$), is still an open problem. One can avoid to seek an explanation for this problem by attributing the lack of power to cosmic variance or can look for explanations i.e., different inflationary potentials or initial conditions for inflation to begin with, non-trivial topology, ISW effect etc. Features in the primordial power spectrum (PPS) motivated by the early universe physics has been the most common solution to address this problem. In the present work we also follow this approach and consider a set of PPS which have features and constrain the parameters of those using WMAP 9 year and Planck data employing Markov-Chain Monte Carlo (MCMC) analysis. The prominent feature of all the models of PPS that we consider is an infra-red cut off which leads to suppression of power at large angular scales. We consider models of PPS with maximum three extra parameters and use Akaike information criterion ($AIC$) of model selection to compare the models. We find that inflationary models with cut off features lead to a better fit of the observed data compared to simple power law model. For most models we find good constraints for the cut off scale $k_c$, however, for other parameters our constraints are not that good. We find that model with sharp cut off in PPS best-fit the WMAP 9 year data and Starabinsky models is the preferred model for the joint WMAP 9 year + Planck data set, which is also able to produce CMB power suppression up to $\ell\leq30$ to some extent.
Radial-velocity observations of Kepler candidates obtained with the SOPHIE and HARPS-N spectrographs have permitted unveiling the nature of the five giant planets Kepler-41b, Kepler-43b, Kepler-44b, Kepler-74b, and Kepler-75b, the massive companion Kepler-39b, and the brown dwarf KOI-205b. These companions were previously characterized with long-cadence (LC) Kepler data. Here we aim at refining the parameters of these transiting systems by i) modelling the published radial velocities (RV) and Kepler short-cadence (SC) data that provide a much better sampling of the transits, ii) performing new spectral analyses of the SOPHIE and ESPaDOnS spectra, and iii) improving stellar rotation periods hence stellar age estimates through gyrochronology, when possible. Posterior distributions of the system parameters were derived with a differential evolution Markov chain Monte Carlo approach. Our main results are as follows: a) Kepler-41b is significantly larger and less dense than previously found because a lower orbital inclination is favoured by SC data. This also affects the determination of the geometric albedo that is lower than previously derived: Ag < 0.135; b) Kepler-44b is moderately smaller and denser than reported in the discovery paper; c) good agreement was achieved with published Kepler-43, Kepler-75, and KOI-205 system parameters, although the host stars Kepler-75 and KOI-205 were found to be slightly richer in metals and hotter, respectively; d) the previously reported non-zero eccentricities of Kepler-39b and Kepler-74b might be spurious. If their orbits were circular, the two companions would be smaller and denser than in the eccentric case. The radius of Kepler-39b is still larger than predicted by theoretical isochrones. Its parent star is hotter and richer in metals than previously determined. [ABRIDGED]
We present metallicity gradients in 49 local field star-forming galaxies. We derive gas-phase oxygen abundances using two widely adopted metallicity calibrations based on the [OIII]/Hbeta, [NII]/Halpha and [NII]/[OII] line ratios. The two derived metallicity gradients are usually in good agreement within +/-0.14 dex/R25 (R25 is the B-band iso-photoal radius), but the metallicity gradients can differ significantly when the ionisation parameters change systematically with radius. We investigate the metallicity gradients as a function of stellar mass (8<log(M*/Msun)<11) and absolute B-band luminosity (-16 > MB > -22). When the metallicity gradients are expressed in dex/kpc, we show that galaxies with lower mass and luminosity, on average, have steeper metallicity gradients. When the metallicity gradients are expressed in dex/R25, we find no correlation between the metallicity gradients, and stellar mass and luminosity. We provide a local benchmark metallicity gradient of field star-forming galaxies useful for comparison with studies at high redshifts. We investigate the origin of the local benchmark gradient using simple chemical evolution models and observed gas and stellar surface density profiles in nearby field spiral galaxies. Our models suggest that the local benchmark gradient is a direct result of the coevolution of gas and stellar disk under virtually closed-box chemical evolution when the stellar-to-gas mass ratio becomes high (>>0.3). These models imply low current mass accretion rates (<0.3xSFR), and low mass outflow rates (<3xSFR) in local field star-forming galaxies.
We explore the science power of space telescopes used to estimate the true masses of known radial-velocity exoplanets by means of astrometry on direct images. We translate a desired mass accuracy (+/10% in our example) into a minimum goal for the signal-to-noise ratio, which implies a minimum exposure time. When the planet is near a node, the mass measurement becomes difficult if not impossible, because the apparent separation becomes decoupled from the inclination angle of the orbit. The combination of this nodal effect with considerations of solar and anti-solar pointing restrictions, photometric and obscurational completeness, and image blurring due to orbital motion, severely limits the observing opportunities, often to only brief intervals in a five-year mission. We compare the science power of four missions, two with external star shades, EXO-S and WFIRST-S, and two with internal coronagraphs, EXO-C and WFIRST-C. The star shades out-perform the coronagraph in this science program by about a factor of three. For both coronagraphs, the input catalog includes 16 RV planets, of which EXO-C could possibly observe 10, of which 6 would have a 90% guarantee of success. Of the same 16 planets, WFIRST-C could possibly observe 12, of which 9 are guaranteed. For both star-shade missions, the input catalog includes 55 planets, of which EXO-S could possibly observe 37, of which 20 are guaranteed. Of the same 55, WFIRST-S could possibly observe 45, of which 30 are guaranteed. The longer spectroscopic exposure times should be easily accommodated for the RV planets with guaranteed success.
In this work we model the quintessence potential in a Taylor series expansion, up to second order, around the present-day value of the scalar field. The field is evolved in a thawing regime assuming zero initial velocity. We use the latest data from the Planck satellite, baryonic acoustic oscillations observations from the Sloan Digital Sky Survey, and Supernovae luminosity distance information from Union$2.1$ to constrain our models parameters, and also include perturbation growth data from WiggleZ. We show explicitly that the growth data does not perform as well as the other datasets in constraining the dark energy parameters we introduce. We also show that the constraints we obtain for our model parameters, when compared to previous works of nearly a decade ago, have not improved significantly. This is indicative of how little dark energy constraints, overall, have improved in the last decade, even when we add new growth of structure data to previous existent types of data.
The recent discovery of a neutron star accretor in the ultra-luminous X-ray source M82 X-2 challenges our understanding of high-mass X-ray binary formation and evolution. By combining binary population synthesis and detailed mass-transfer models, however, we show that the binary parameters of M82 X-2 are not surprising provided non-conservative mass transfer is allowed. Specifically, the donor-mass lower limit and orbital period measured for M82 X-2 lie near the most probable values predicted by population synthesis models, and systems such as M82 X-2 should be observed in approximately 13% of the galaxies with a star-formation history similar to M82. We conclude that the binary system that formed M82 X-2 is most likely younger than 50 Myr old and contain a donor star which had an initial mass of approximately 8-10 M$_\odot$, while the NS's progenitor star had an initial mass in the $8-25\,\rm M_{\odot}$ range. The donor star still currently resides on the main sequence, and is capable of continued MT on the thermal timescale, while in the ultra-luminous X-ray regime, for as long as 400,000 years.
The standing slow magneto-acoustic oscillations in cooling coronal loops are investigated. There are two damping mechanisms which are considered to generate the standing acoustic modes in coronal magnetic loops namely thermal conduction and radiation. The background temperature is assumed to change temporally due to optically thin radiation. In particular, the background plasma is assumed to be radiatively cooling. The effects of cooling on longitudinal slow MHD modes is analytically evaluated by choosing a simple form of radiative function that ensures the temperature evolution of the background plasma due to radiation coincides with the observed cooling profile of coronal loops. The assumption of low-beta plasma leads to neglect the magnetic field perturbation and eventually reduces the MHD equations to a 1D system modelling longitudinal MHD oscillations in a cooling coronal loop. The cooling is assumed to occur on a characteristic time scale much larger than the oscillation period that subsequently enables using the WKB theory to study the properties of standing wave. The governing equation describing the time-dependent amplitude of waves is obtained and solved analytically. The analytically derived solutions are numerically evaluated to give further insight into the evolution of the standing acoustic waves. We find that the plasma cooling gives rise to a decrease in the amplitude of oscillations. In spite of the reduction in damping rate caused by rising the cooling, the damping scenario of slow standing MHD waves strongly increases in hot coronal loops.
In this paper the next attempt is made to clarify the nature of the Euclidean behavior of the boundary in the angular size-redshift cosmological test. It is shown experimentally that this can be explained by the selection determined by anisotropic morphology and anisotropic radiation of extended radio sources. A catalogue of extended radio sources with minimal flux densities of about 0.01 Jy at 1.4 GHz was compiled for conducting the test. Without the assumption of their size evolution, the agreement between the experiment and calculation was obtained both in the Lambda CDM model (Omega_m=0.27 , Omega_v=0.73.) and the Friedman model (Omega = 0.1 ).
The concept of Space Manifold Dynamics is a new method of space research. We have applied it along with the basic idea of the method of Ott, Grebogi and York (OGY method) to stabilize the motion of a spacecraft around the triangular Lagrange point L5 of the Earth--Moon system. We have determined the escape rate of the trajectories in the general three- and four-body problem and estimated the average lifetime of the particles. Integrating the two models we mapped in detail the phase space around the L5 point of the Earth--Moon system. Using the phase space portrait our next goal was to apply a modified OGY method to keep a spacecraft close to the vicinity of L5. We modified the equation of motions with the addition of a time dependent force to the motion of the spacecraft. In our orbit--keeping procedure there are three free parameters: (i) the magnitude of the thrust, (ii) the start time and (iii) the length of the control. Based on our numerical experiments we were able to determine possible values for these parameters and successfully apply a control phase to a spacecraft to keep it on orbit around L5.
Despite the success of the combination of conservative schemes and staggered
constrained transport algorithms in the last fifteen years, the accurate
description of highly magnetized, relativistic flows with strong shocks
represents still a challenge in numerical RMHD. The present paper focusses in
the accuracy and robustness of several correction algorithms for the conserved
variables, which has become a crucial ingredient in the numerical simulation of
problems where the magnetic pressure dominates over the thermal pressure by
more than two orders of magnitude.
Two versions of non-relativistic and fully relativistic corrections have been
tested and compared using a magnetized cylindrical explosion with high
magnetization ($ \ge 10^4$) as test. In the non-relativistic corrections, the
total energy is corrected for the difference in the classical magnetic energy
term between the average of the staggered fields and the conservative ones,
before (CA1) and after (CA1') recovering the primitive variables. These
corrections are unable to pass the test at any numerical resolution. The two
relativistic approaches (CA2 and CA2'), correcting also the magnetic terms
depending on the flow speed in both the momentum and the total energy, reveal
as much more robust. These algorithms pass the test succesfully and with very
small deviations of the energy conservation ($\le 10^{-4}$), and very low
values of the total momentum ($\le 10^{-8}$). In particular, the algorithm CA2'
(that corrects the conserved variables after recovering the primitive
variables) passes the test at all resolutions.
The numerical code used to run all the test cases is briefly described.
Context. The relation between solar surface rotation and sunspot activity
still remains open. Sunspot activity has dramatically reduced in solar cycle 24
and several solar activity indices and flux measurements experienced
unprecedentedly low levels during the last solar minimum.
Aims. We aim to reveal the momentary variation of solar surface rotation,
especially during the recent years of reducing solar activity. Methods. We used
a dynamic, differentially rotating reference system to determine the best-fit
annual values of the differential rotation parameters of active longitudes of
solar X-ray flares and sunspots in 1977-2012.
Results. The evolution of rotation of solar active longitudes obtained with
X-ray flares and with sunspots is very similar. Both hemispheres speed up since
the late 1990s, with the southern hemisphere rotating slightly faster than the
north. Earlier, in 1980s, rotation in the northern hemisphere was considerably
faster, but experienced a major decrease in the early 1990s. On the other hand,
little change was found in the southern rotation during these decades. This led
to a positive asymmetry in north-south rotation rate in the early part of the
time interval studied.
Conclusions. The rotation of both hemispheres has been speeding up at roughly
the same rate since late 1990s, with the southern hemisphere rotating slightly
faster than the north. This period coincides with the start of dramatic
weakening of solar activity, as observed in sunspots and several other solar,
interplanetary and geomagnetic parameters.
We introduce Coronal-Line Forest Active Galactic Nuclei (CLiF AGN), AGN which have a rich spectrum of forbidden high-ionization lines (FHILs, e.g. [FeVII], [FeX] and [NeV]), as well as relatively strong narrow ($\sim$300 km s$^{-1}$) H$\alpha$ emission when compared to the other Balmer transition lines. We find that the kinematics of the CLiF emitting region are similar to those of the forbidden low-ionization emission-line (FLIL) region. We compare emission line strengths of both FHILs and FLILs to CLOUDY photoionization results and find that the CLiF emitting region has higher densities (10$^{4.5}$ $<$ n$_H$ $<$ 10$^{7.5}$ cm$^{-3}$) when compared to the FLIL emitting region (10$^{3.0}$ $<$ n$_H$ $<$ 10$^{4.5}$ cm$^{-3}$). We use the photoionization results to calculate the CLiF regions radial distances (0.04 $<$ R$_{CLiF}$ $<$ 32.5 pc) and find that they are comparable to the dust grain sublimation distances (0.10 $<$ R$_{SUB}$ $<$ 4.3 pc). As a result we suggest that the inner torus wall is the most likely location of the CLiF region, and the unusual strength of the FHILs is due to a specific viewing angle giving a maximal view of the far wall of the torus without the continuum being revealed.
Prompt gamma-ray and early X-ray afterglow emission in gamma-ray bursts (GRBs) are characterized by a bursty behavior and are often interspersed with long quiescent times. There is compelling evidence that X-ray flares are linked to prompt gamma-rays. However, the physical mechanism that leads to the complex temporal distribution of gamma-ray pulses and X-ray flares is not understood. Here we show that the waiting time distribution (WTD) of pulses and flares exhibits a power-law tail extending over 4 decades with index ~2 and can be the manifestation of a common time-dependent Poisson process. This result is robust and is obtained on different catalogs. Surprisingly, GRBs with many (>=8) gamma-ray pulses are very unlikely to be accompanied by X-ray flares after the end of the prompt emission (3.1 sigma Gaussian confidence). These results are consistent with a simple interpretation: an hyperaccreting disk breaks up into one or a few groups of fragments, each of which is independently accreted with the same probability per unit time. Prompt gamma-rays and late X-ray flares are nothing but different fragments being accreted at the beginning and at the end, respectively, following the very same stochastic process and likely the same mechanism.
We study the effects of a large-scale, ordered magnetic field in protoplanetary disks on Type I planet migration using a linear perturbation analysis in the ideal-MHD limit. We focus on wind-driving disks, in which a magnetic torque $\propto B_{0z} \partial B_{0\varphi}/\partial z$ (where $B_{0z}$ and $B_{0\varphi}$ are the equilibrium vertical and azimuthal field components) induces vertical angular momentum transport. We derive the governing differential equation for the disk response and identify its resonances and turning points. For a disk containing a slightly subthermal, pure-$B_{0z}$ field, the total 3D torque is close to its value in the 2D limit but remains lower than the hydrodynamic torque. In contrast with the 2D pure-$B_{0\varphi}$ field model considered by Terquem (2003), inward migration is not reduced in this case when the field amplitude decreases with radius. The presence of a subdominant $B_{0\varphi}$ component whose amplitude increases from zero at $z=0$ has little effect on the torque when acting alone, but in conjunction with a $B_{0z}$ component it gives rise to a strong torque that speeds up the inward migration by a factor $\gtrsim 200$. This factor could, however, be reduced in a real disk by dissipation and magnetic diffusivity effects. Unlike all previously studied disk migration models, in the $B_{0z}+\partial B_{0\varphi}/\partial z$ case the dominant contributions to the torque add with the same sign from the two sides of the planet. We attribute this behavior to a new mode of interaction wherein a planet moves inward by plugging into the disk's underlying angular momentum transport mechanism.
The GJ 436 system contains a transiting planet GJ 436 b which is a hot analogue of Neptune on an eccentric orbit. Recently, two additional transiting sub-Earth planets have been postulated in the literature. We observed three transits of GJ 436 b over the course of 3 years using two-meter class telescopes, each with a photometric precision better than one millimagnitude. We studied system dynamics based on the existence of the additional planets. We redetermined system parameters, which were in agreement with those found in the literature. We refined the orbital period of GJ 436 b and found no evidence of transit timing variations. The orbital motion of the GJ 436 c planet candidate was found to be significantly affected by the planet b with variations in transit times at a level of 20 minutes. As the orbital period of the GJ 436 d planet candidate remains unknown, our numerical experiments rule out orbits in low-order resonances with GJ 436 b. The GJ 436 system with the hot Neptune and additional two Earth-like planets, if confirmed, would be an important laboratory for studies of formation and evolution of planetary systems.
This paper shows that the inclinations of bodies captured into mean motion resonances in the Kuiper belt have remained very nearly unchanged, being only slightly increased from initial lower values by migration and/or by long-term planetary perturbations. Thus the observed maximum as high as ~ 30 deg of the i's of bodies in resonance must reflect either a broad initial range at least to that level for capturable bodies or an elevating process possibly exemplified by the sweeping of secular resonances. We have obtained capture probabilities for 2 well-populated resonances, showing reduced but finite values for i's up to 35 deg. Whatever led to the present distribution must have produced increases in i for some, but not for all, resonant bodies in the belt.
Stellar black hole (BH) binaries are one of the most promising GW sources for GW detection by the ground-based detectors. Nuclear star clusters (NCs) located at the center of galaxies are known to harbor massive black holes (MBHs) and to be bounded by a gravitational potential by other galactic components such as the galactic bulge. Such an environment of NCs provides a favorable conditions for the BH-BH binary formation by the gravitational radiation (GR) capture due to the high BH number density and velocity dispersion. We carried out detailed numerical study of the formation of BH binaries in the NCs using a series of N-body simulations for equal-mass cases. There is no mass segregation introduced. We have derived scaling relations of the binary formation rate with the velocity dispersion of the stellar system beyond the radius of influence and made estimates of the rate of formation of black hole binaries per unit comoving volume and thus expected detection rate by integrating the binary formation rate over galaxy population within the detection distance of the advanced detectors. We find that the overall formation rates for BH-BH binaries per NC is 10^(-10)/yr for the Milky-Way-like galaxies and weakly dependent on the mass of MBH as M^(3/28). We estimate the detection rate of 0.02-14/yr for advanced LIGO/Virgo considering several factors such as the dynamical evolution of NCs, the variance of the number density of stars and the mass range of MBH giving uncertainties.
The interpretation of gravitational microlensing events caused by planetary systems or multiple stars is based on the n-point mass lens model. The first planets detected by microlensing were well described by the two-point-mass model of a star with one planet. By the end of 2014, four events involving three-point-mass lenses had been announced. Two of the lenses were stars with two planetary companions each; two were binary stars with a planet orbiting one component. While the two-point-mass model is well understood, the same cannot be said for lenses with three or more components. Even the range of possible critical-curve topologies and caustic geometries of the three-point-mass lens remains unknown. In this paper we provide new tools for mapping the critical-curve topology and caustic cusp number in the parameter space of n-point-mass lenses. We perform our analysis in the image plane of the lens. We show that all contours of the Jacobian are critical curves of re-scaled versions of the lens configuration. Utilizing this property further, we introduce the cusp curve to identify cusp-image positions on all contours simultaneously. In order to track cusp-number changes in caustic metamorphoses, we define the morph curve, which pinpoints the positions of metamorphosis-point images along the cusp curve. We demonstrate the usage of both curves on simple two- and three-point-mass lens examples. For the three simplest caustic metamorphoses we illustrate the local structure of the image and source planes.
We report the results of a comprehensive reanalysis of Ulysses observations of interstellar He atoms flowing through the solar system, the goal being to reassess the interstellar He flow vector and to search for evidence of variability in this vector. We find no evidence that the He beam seen by Ulysses changes at all from 1994-2007. The direction of flow changes by no more than ~0.3 deg and the speed by no more than ~0.3 km/s. A global fit to all acceptable He beam maps from 1994-2007 yields the following He flow parameters: V_ISM=26.08+/-0.21 km/s, lambda=75.54+/-0.19 deg, beta=-5.44+/-0.24 deg, and T=7260+/-270 K; where lambda and beta are the ecliptic longitude and latitude direction in J2000 coordinates. The flow vector is consistent with the original analysis of the Ulysses team, but our temperature is significantly higher. The higher temperature somewhat mitigates a discrepancy that exists in the He flow parameters measured by Ulysses and the Interstellar Boundary Explorer, but does not resolve it entirely. Using a novel technique to infer photoionization loss rates directly from Ulysses data, we estimate a density of n_He=0.0196+/-0.0033 cm^-3 in the interstellar medium.
We probe the column densities and masses traced by the ionized and neutral atomic carbon with spectrally resolved maps, and compare them to the diffuse and dense molecular gas traced by [C I] and low-$J$ CO lines toward the star-forming region M17SW. We mapped a 4.1pc x 4.7pc region in the [C I] 609 m$\mu$ line using the APEX telescope, as well as the CO isotopologues with the IRAM 30m telescope. We analyze the data based on velocity channel maps that are 1 km/s wide. We correlate their spatial distribution with that of the [C II] map obtained with SOFIA/GREAT. Optically thin approximations were used to estimate the column densities of [C I] and [C II] in each velocity channel. The spatial distribution of the [C I] and all CO isotopologues emission was found to be associated with that of [C II] in about 20%-80% of the mapped region, with the high correlation found in the central (15-23 km/s ) velocity channels. The excitation temperature of [C I] ranges between 40 K and 100 K in the inner molecular region of M17 SW. Column densities in 1 km/s channels between ~10$^{15}$ and ~10$^{17}$ cm$^{-2}$ were found for [C I]. Just ~20% of the velocity range (~40 km/s) that the [C II] line spans is associated with the star-forming material traced by [C I] and CO. The total gas mass estimated from the [C II] emission gives a lower limit of ~4.4x10$^3$ $M_{\odot}$. At least 64% of this mass is not associated with the star-forming material in M17SW. We also found that about 36%, 17%, and 47% of the [C II] emission is associated with the HII, HI, and H_2 regimes, respectively. Comparisons with the H41$\alpha$ line shows an ionization region mixed with the neutral and part of the molecular gas, in agreement with the clumped structure and dynamical processes at play in M17SW. These results are also relevant to extra-galactic studies in which [C II] is often used as a tracer of star-forming material.
We study the effects of a large-scale, ordered magnetic field in protoplanetary disks on Type I planet migration using a combination of numerical simulations in 2D and 3D and a linear perturbation analysis. Steady-state models of such disks require the inclusion of magnetic diffusivity. To make progress using ideal MHD, we focus on simplified field configurations, involving purely vertical ($B_z$) and azimuthal ($B_\varphi$) field components and a combination of the two. For each of the models we calculate the locations of the relevant resonances and of the turning points, which delineate the propagation regions of the MHD waves that transport angular momentum from the planet to the disk. We use both numerical and semianalytic methods to evaluate the cumulative back torque acting on the planet, and explore the effect of spatial gradients in the disk's physical variables on the results. We conclude that, under realistic (3D) circumstances, a large-scale magnetic field can slow down the inward migration that characterizes the underlying unmagnetized disk --- by up to a factor of $\sim 2$ when the magnetic pressure approaches the thermal pressure --- but it cannot reverse it. A previous inference that a pure-$B_\phi$ field whose amplitude decreases fast enough with radius leads to outward migration applies only in 2D. In fact, we find that, in 3D, a pure-$B_\phi$ disk undergoes a rapid transition to turbulence on account of a magnetorotational instability that is triggered by the planet-induced appearance of a weak $B_z$ component.
We present new Spitzer/IRAC observations of 55 dusty Long Period Variables (LPVs, 48 AGB and 6 RSG stars) in the Galaxy that have different chemistry, variability type, and mass-loss rate. O-rich AGB stars (including intrinsic S-type) tend to have redder [3.6]-[8.0] colors than carbon stars for a given [3.6]-[4.5] color due to silicate features increasing the flux in the 8.0 {\mu}m IRAC band. For colors including the 5.8 {\mu}m band, carbon stars separate into two distinct sequences, likely due to a variable photospheric C$_3$ feature that is only visible in relatively unobscured, low mass-loss rate sources. Semiregular variables tend to have smaller IR excess in [3.6]-[8.0] color than Miras, consistent with the hypothesis that semiregular variables lose mass discontinuously. Miras have redder colors for longer periods while semiregular variables do not. Galactic AGB stars follow the period-luminosity sequences found for the Magellanic Clouds. Mira variables fall along the fundamental pulsation sequence, while semiregular variables are mostly on overtone sequences. We also derive a relationship between mass-loss rate and [3.6]-[8.0] color. The fits are similar in shape to those found by other authors for AGBs in the LMC, but discrepant in overall normalization, likely due to different assumptions in the models used to derive mass-loss rates. We find that IR colors are not unique discriminators of chemical type, suggesting caution when using color selection techniques to infer the chemical composition of AGB dust returned to the ISM.
A new solar feature termed a dark jet is identified from observations of an extended solar coronal hole that was continuously monitored for over 44 hours by the EUV Imaging Spectrometer on board the Hinode spacecraft in 2011 February 8-10. Line-of-sight velocity maps derived from the coronal Fe XII $\lambda$195.12 emission line, formed at 1.5 MK, revealed a number of large-scale, jet-like structures that showed significant blueshifts. The structures had either weak or no intensity signal in 193 A filter images from the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory, suggesting that the jets are essentially invisible to imaging instruments. The dark jets are rooted in bright points and occur both within the coronal hole and at the quiet Sun-coronal hole boundary. They exhibit a wide range of shapes, from narrow columns to fan-shaped structures, and sometimes multiple jets are seen close together. A detailed study of one dark jet showed line-of-sight speeds increasing along the jet axis from 52 to 107 km s$^{-1}$ and a temperature of 1.2-1.3 MK. The low intensity of the jet was due either to a small filling factor of 2% or to a curtain-like morphology. From the HOP 177 sample, dark jets are as common as regular coronal hole jets, but their low intensity suggests a mass flux around two orders of magnitude lower.
We perform a blind search for the variability of the gamma-ray sky in the energy range E>1 GeV using 308 weeks of the Fermi-LAT data. We use the technique based on the comparison of the weekly photon counts and exposures in sky pixels by means of the Kolmogorov-Smirnov test. We consider the flux variations in the region significant if statistical probability of uniformity is less than $4\times10^{-6}$, which corresponds to 0.05 false detections in the whole set of 12288 pixels. Close inspection of the detected variable regions result in identification of 8 sources without previous known variability. Two of them are included in the second Fermi LAT source catalogue (FBQS J122424.1+243623 and GB6 J0043+3426) and one (3EG J1424+3734) was reported by EGRET and also was included in the First Fermi LAT source catalogue (1FGL), but is missing in the 2FGL. Possible identifications of five other sources are obtained using NED and SIMBAD databases (1RXS J161939.9+765515, PMN J2320-6447, PKS 0226-559, PKS J0030-0211, PMN J0225-2603). These new variable gamma-ray sources demonstrate recurring flaring activity with time scale ~weeks and have hard spectra. Their spectral energy distributions deviate significantly from a simple power-law shape and often peak around ~GeV. These properties of activity are typical for flaring blazars.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of the physics of accretion and ejection around compact objects. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of accreting white dwarfs. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of active galactic nuclei. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of very faint X-ray binaries, orbital period distribution of black hole X-ray binaries and neutron star spin up. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of radio-loud Active Galactic Nuclei. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of stellar flares. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of gamma-ray bursts. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of magnetospheres of isolated neutron stars. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of jetted tidal disruption events. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of high-energy radiation from thunderstorms and lightning. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper.
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of high-mass X-ray binaries and ultra-luminous X-ray sources. For a summary, we refer to the paper.
The minimal sub-Planckian axion inflation model accounts for a large scalar-to-tensor ratio via a spiralling trajectory in the field space of a complex field $\Phi$. Here we consider how the predictions of the model are modified by Planck scale-suppressed corrections. In the absence of Planck corrections the model is equivalent to a $\phi^{4/3}$ chaotic inflation model. Planck corrections become important when the dimensionless coupling $\xi$ of $|\Phi|^{2}$ to the topological charge of the strongly-coupled gauge sector $F \tilde{F}$ satisfies $\xi \sim 1$. For values of $|\Phi|$ which allow the Planck corrections to be understood via an expansion in powers of $|\Phi|^{2}/M_{Pl}^{2}$, we show that their effect is produce a significant modification of the tensor-to-scalar ratio from its $\phi^{4/3}$ chaotic inflation value without strongly modifying the spectral index. In addition, to leading order in $|\Phi|^2/M_{Pl}^{2}$, the Planck modifications of $n_{s}$ and $r$ satisfy a consistency relation, $\Delta n_{s} = - \Delta r/16$. Observation of these shifts and their correlation would allow the model to be distinguished from a simple $\phi^{4/3}$ chaotic inflation model and would also provide a signature for the influence of leading-order Planck corrections.
In a paper posted on the arXiv a few weeks ago Berti, Brito and Cardoso \cite{Berti+14} suggest that ultra-high-energy particles can emerge from collisions in a black hole's ergosphere. This can happen if the process involves a particle on an outgoing trajectory very close to the black hole. Clearly such a particle cannot emerge from the black hole. It is argued \cite{Berti+14} that this particle can arise in another collision. Thus the process involves two collisions: one in which an outgoing particle is produced extremely close to the horizon, and a second one in which energy is gained. The real efficiency of this process should take into account, therefore, the energy needed to produce the first particle. We show here that while this process is kinematically possible, it requires a deposition of energy that is divergently large compared with the energy of the escaping particle. Thus, in contradiction to claims of infinitely high efficiencies, the efficiency of the combined process is in fact extremely small, approaching zero for very high output energies. Even under more general conditions than those considered in \cite{Berti+14} the total energy gain never diverges, and is larger only by a factor of a few than the energy gain of the original collisional Penrose process that takes place between two infalling particles \cite{Piran+75,PiranShaham77,Bejger+12}
Using a framework based on the 1+3 formalism we carry out a study on axially and reflection symmetric perfect and geodesic fluids, looking for possible models of sources radiating gravitational waves. Therefore, the fluid should be necessarily shearing, for otherwise the magnetic part of the Weyl tensor vanishes, leading to a vanishing of the super-Poynting vector. However, for the family of perfect, geodesic fluids considered here, it appears that all possible cases reduce to conformally flat, shear--free, vorticity-free, fluids, i.e Friedmann-Roberston-Walker. The super-Poynting vector vanishes and therefore no gravitational radiation is expected to be produced. The physical meaning of the obtained result is discussed.
We propose an extension of natural inflation, where the inflaton potential is a general periodic function. Specifically, we study elliptic inflation where the inflaton potential is given by Jacobi elliptic functions or Jacobi theta functions, which appear in gauge and Yukawa couplings in the string theories compactified on toroidal backgrounds. We show that the predicted values of the spectral index and the tensor-to-scalar ratio interpolate from natural inflation to exponential inflation such as $R^2$- and Higgs inflation or brane inflation, where the spectral index asymptotes to $n_s = 1-2/N \simeq 0.967$ for the e-folding number $N = 60$. Such elliptic inflation can be thought of as a specific realization of multi-natural inflation, where the inflaton potential consists of multiple sinusoidal functions. We also discuss examples in string theory where the Jacobi theta function appears in the inflaton potential.
The lore paradigm for solving so-called horizon and flatness problems in cosmology is the primordial inflation. Plethora of inflationary models have been built in last decades and first experimental probes seem to appear in favor of the inflationary paradigm. We will focus here on one of them, the Higgs inflation, and show the combined constraint required for such a model at cosmological as well as gravitational scales, i.e. for compact objects. We will show that Higgs inflation model gives rise to particlelike solutions around compact objects, dubbed Higgs monopoles, characterized by the nonminimal coupling parameter as well as the mass and the compactness of the object. For large values of the nonminimal coupling constant and specific compactness, the amplitude of the Higgs field inside the matter distribution can be arbitrarily large.
The galilean genesis scenario is an alternative to inflation in which the universe starts expanding from Minkowski in the asymptotic past by violating the null energy condition stably. Several concrete models of galilean genesis have been constructed so far within the context of galileon-type scalar-field theories. We give a generic, unified description of the galilean genesis scenario in terms of the Horndeski theory, i.e., the most general scalar-tensor theory with second-order field equations. In doing so we generalize the previous models to have a new parameter (denoted by {\alpha}) which results in controlling the evolution of the Hubble rate. The background dynamics is investigated to show that the generalized galilean genesis solution is an attractor, similarly to the original model. We also study the nature of primordial perturbations in the generalized galilean genesis scenario. In all the models described by our generalized genesis Lagrangian, amplification of tensor perturbations does not occur as opposed to what happens in quasi-de Sitter inflation. We show that the spectral index of curvature perturbations is determined solely from the parameter {\alpha} and does not depend on the other details of the model. In contrast to the original model, a nearly scale-invariant spectrum of curvature perturbations is obtained for a specific choice of {\alpha}.
IceCube has measured a diffuse astrophysical flux of TeV-PeV neutrinos. The most plausible sources are unique high energy cosmic ray accelerators like hypernova remnants (HNRs) and remnants from gamma ray bursts in star-burst galaxies, which can produce primary cosmic rays with the required energies and abundance. In this case, however, ordinary supernova remnants (SNRs), which are far more abundant than HNRs, produce a comparable or larger neutrino flux in the ranges up to 100-150 TeV energies, implying a spectral break in the IceCube signal around these energies. The SNRs contribution in the diffuse flux up to these hundred TeV energies provides a natural baseline and then constrains the expected PeV flux.
Recent work in the literature has shown that the leading long distance quantum corrections to the Newtonian potential imply tiny but observable effects in the restricted three-body problem of celestial mechanics, i.e., at the Lagrangian libration points of stable equilibrium the planetoid is not exactly at equal distance from the two bodies of large mass, but the Newtonian values of its coordinates are changed by a few millimeters in the Earth-Moon system. First, we assess such a theoretical calculation by exploiting the full theory of the quintic equation, i.e., its reduction to Bring-Jerrard form and the resulting expression of roots in terms of generalized hypergeometric functions. By performing the numerical analysis of the exact formulas for the roots, we confirm and slightly improve the theoretical evaluation of quantum corrected coordinates of Lagrangian libration points of stable equilibrium. Second, we discuss the prospects to measure, with the help of laser ranging, the above departure from the equilateral triangle picture, which is a challenging task. On the other hand, a modern version of the planetoid is the solar sail, and much progress has been made, in recent years, on the displaced periodic orbits of solar sails at all libration points, both stable and unstable. The present paper investigates therefore, eventually, a restricted three-body problem involving Earth, Moon and a solar sail. By taking into account the quantum corrections to the Newtonian potential, displaced periodic orbits of the solar sail at libration points are again found to exist.
We study the cosmology of bimetric theory with a composite matter coupling. We find two possible branches of background evolution. We investigate the question of stability of cosmological perturbations. For the tensor and vector perturbations, we derive conditions on the absence of ghost and gradient instabilities. For the scalar modes, we obtain conditions for avoiding ghost degrees. In the first branch, we find that one of the scalar modes becomes a ghost at the late stages of the evolution. Conversely, this problem can be avoided in the second branch. However, we also find that the constraint for the second branch prevents the doubly coupled matter fields from being the standard ingredients of cosmology. We thus conclude that a realistic and stable cosmological model requires additional minimally coupled matter fields.
Links to: arXiv, form interface, find, astro-ph, recent, 1501, contact, help (Access key information)
The redshift-space distortion (RSD) of galaxies surrounding massive clusters is emerging as a promising testbed for theories of modified gravity. Conventional applications of this method rely upon the assumption that the velocity field in the cluster environment is uniquely determined by the cluster mass profile. Yet, real dark matter halos in N-body simulations are known to violate the assumption that virial mass determines the configuration space distribution, an effect known as assembly bias. In this Letter, I show that assembly bias in simulated dark matter halos also manifests in velocity space. In the 1-10 Mpc environment surrounding a cluster, high-concentration "tracer" halos exhibit a 10-20% larger pairwise-velocity dispersion profile relative to low-concentration tracer halos of the same mass. This difference is comparable to the size of the RSD signal predicted by f(R) models designed to account for the cosmic acceleration. I use the age matching technique to study how color-selection effects may influence the cluster RSD signal, finding a ~10% effect due to redder satellites preferentially occupying higher mass halos, and a ~5% effect due to assembly-biased colors of centrals. In order to use cluster RSD measurements to robustly constrain modified gravity, we likely will need to develop empirical galaxy formation models more sophisticated than any in the current literature.
We model the chemical evolution of six UFDs: Bootes I, Canes Venatici II, Coma Berenices, Hercules, Leo IV and Ursa Major I, based on their recently determined star formation histories. We show that two single-age bursts cannot explain the observed [$\alpha$/Fe] vs [Fe/H] distribution in these galaxies and that some self-enrichment is required within the first burst. An alternative scenario is modelled, in which star formation is continuous except for short interruptions when one or more supernovae temporarily blow the dense gas out from the centre of the system. This model allows for self-enrichment and can reproduce the chemical abundances of the UFDs in which the second burst is only a trace population. We conclude that the most likely star formation history is one or two extended periods of star formation, with the first burst lasting for at least 100~Myr. As found in earlier work, the observed properties of UFDs can be explained by formation at a low mass ($M_{\rm{vir}}\sim10^7$~M$_\odot$), rather than being stripped remnants of much larger systems.
The mass and structural evolution of massive galaxies is one of the hottest topics in galaxy formation. This is because it may reveal invaluable insights into the still debated evolutionary processes governing the growth and assembly of spheroids. However, direct comparison between models and observations is usually prevented by the so-called "progenitor bias", i.e., new galaxies entering the observational selection at later epochs, thus eluding a precise study of how pre-existing galaxies actually evolve in size. To limit this effect, we here gather data on high-redshift brightest group and cluster galaxies, evolve their (mean) host halo masses down to z=0 along their main progenitors, and assign as their "descendants" local SDSS central galaxies matched in host halo mass. At face value, the comparison between high redshift and local data suggests a noticeable increase in stellar mass of a factor of >2 since z~1, and of >2.5 in mean effective radius. We then compare the inferred stellar mass and size growth with those predicted by hierarchical models for central galaxies, selected at high redshifts to closely match the halo and stellar mass bins as in the data. Only hierarchical models characterized by very limited satellite stellar stripping and parabolic orbits are capable of broadly reproducing the stellar mass and size increase of a factor ~2-4 observed in cluster galaxies since z ~1. The predicted, average (major) merger rate since z~1 is in good agreement with the latest observational estimates.
We use 317,000 emission-line galaxies from the Sloan Digital Sky Survey to investigate line-ratio selection of active galactic nuclei (AGNs). In particular, we demonstrate that "star formation dilution" by HII regions causes a significant bias against AGN selection in low-mass, blue, star-forming, disk-dominated galaxies. This bias is responsible for the observed preference of AGNs among high-mass, green, moderately star-forming, bulge-dominated hosts. We account for the bias and simulate the intrinsic population of emission-line AGNs using a physically-motivated Eddington ratio distribution, intrinsic AGN narrow line region line ratios, a luminosity-dependent Lbol/L[OIII] bolometric correction, and the observed Mbh-sigma relation. These simulations indicate that, in massive (log(M*/Msun) > 10) galaxies, AGN accretion is correlated with specific star formation rate but is otherwise uniform with stellar mass. There is some hint of lower black hole occupation in low-mass (log(M*/Msun) < 10) hosts, although our modeling is limited by uncertainties in measuring and interpreting the velocity dispersions of low-mass galaxies. The presence of star formation dilution means that AGNs contribute little to the observed strong optical emission lines (e.g., [OIII] and Ha) in low-mass and star-forming hosts. However the AGN population recovered by our modeling indicates that AGN feedback has nearly uniform efficiency at all stellar masses, star formation rates, and morphologies. Taken together, our characterization of the observational bias and resultant AGN occupation function suggest that AGNs are unlikely to be the dominant source of star formation quenching in galaxies, but instead are fueled by the same gas which drives star formation activity.
Interaction of charges in CCDs with the already accumulated charge distribution causes both a flux dependence of the point-spread function (an increase of observed size with flux, also known as the brighter/fatter effect) and pixel-to-pixel correlations of the noise in flat fields. We describe these effects in the Dark Energy Camera (DECam) with charge dependent shifts of effective pixel borders, i.e. the Antilogus et al. (2014) model, which we fit to measurements of flat-field noise correlations. The latter fall off approximately as a power-law r^-2.5 with pixel separation r, are isotropic except for an asymmetry in the direct neighbors along rows and columns, are stable in time, and are weakly dependent on wavelength. They show variations from chip to chip at the 20% level that correlate with the silicon resistivity. The charge shifts predicted by the model cause biased shape measurements, primarily due to their effect on bright stars, at levels exceeding weak lensing science requirements. We measure the flux dependence of star images and show that the effect can be mitigated by applying the reverse charge shifts at the pixel level during image processing. Differences in stellar size, however, remain significant due to residuals at larger distance from the centroid.
We aim to constrain the evolution of AGN as a function of obscuration using an X-ray selected sample of $\sim2000$ AGN from a multi-tiered survey including the CDFS, AEGIS-XD, COSMOS and XMM-XXL fields. The spectra of individual X-ray sources are analysed using a Bayesian methodology with a physically realistic model to infer the posterior distribution of the hydrogen column density and intrinsic X-ray luminosity. We develop a novel non-parametric method which allows us to robustly infer the distribution of the AGN population in X-ray luminosity, redshift and obscuring column density, relying only on minimal smoothness assumptions. Our analysis properly incorporates uncertainties from low count spectra, photometric redshift measurements, association incompleteness and the limited sample size. We find that obscured AGN with $N_{H}>{\rm 10^{22}\, cm^{-2}}$ account for ${77}^{+4}_{-5}\%$ of the number density and luminosity density of the accretion SMBH population with $L_{{\rm X}}>10^{43}\text{ erg/s}$, averaged over cosmic time. Compton-thick AGN account for approximately half the number and luminosity density of the obscured population, and ${38}^{+8}_{-7}\%$ of the total. We also find evidence that the evolution is obscuration-dependent, with the strongest evolution around $N_{H}\thickapprox10^{23}\text{ cm}^{-2}$. We highlight this by measuring the obscured fraction in Compton-thin AGN, which increases towards $z\sim3$, where it is $25\%$ higher than the local value. In contrast the fraction of Compton-thick AGN is consistent with being constant at $\approx35\%$, independent of redshift and accretion luminosity. We discuss our findings in the context of existing models and conclude that the observed evolution is to first order a side-effect of anti-hierarchical growth.
We present a census of blue horizontal branch (BHB) and blue straggler (BS) stars belonging to dwarf galaxies and globular clusters, and compare these counts to that of the Milky Way stellar halo. We find, in agreement with earlier studies, that the ratio of BS-to-BHB stars in these satellite populations is dependent on stellar mass. Dwarf galaxies show an increasing BS-to-BHB ratio with luminosity. In contrast, globular clusters display the reverse trend, with N_BS/N_BHB (< 1) decreasing with luminosity. The faintest (L < 10^5 L_Sun) dwarfs have similar numbers of BS and BHB stars (N_BS/N_BHB ~ 1), whereas more massive dwarfs tend to be dominated by BS stars (N_BS/N_BHB ~ 2-40). We find that the BS-to-BHB ratio in the stellar halo is relatively high (N_BS/N_BHB ~ 5-6), and thus inconsistent with the low ratios found in both ultra-faint dwarfs and globular clusters. Our results favour more massive dwarfs as the dominant "building blocks" of the stellar halo, in good agreement with current predictions from LambdaCDM models.
We present an in-depth analysis of stellar activity and its effects on radial velocity (RV) for the M2 dwarf GJ 176 based on spectra taken over 10 years from the High Resolution Spectrograph on the Hobby-Eberly Telescope. These data are supplemented with spectra from previous observations with the HIRES and HARPS spectrographs, and V- and R-band photometry taken over 6 years at the Dyer and Fairborn observatories. Previous studies of GJ 176 revealed a super-Earth exoplanet in an 8.8-day orbit. However, the velocities of this star are also known to be contaminated by activity, particularly at the 39-day stellar rotation period. We have examined the magnetic activity of GJ 176 using the sodium I D lines, which have been shown to be a sensitive activity tracer in cool stars. In addition to rotational modulation, we see evidence of a long-term trend in our Na I D index, which may be part of a long-period activity cycle. The sodium index is well correlated with our RVs, and we show that this activity trend drives a corresponding slope in RV. Interestingly, the rotation signal remains in phase in photometry, but not in the spectral activity indicators. We interpret this phenomenon as the result of one or more large spot complexes or active regions which dominate the photometric variability, while the spectral indices are driven by the overall magnetic activity across the stellar surface. In light of these results, we discuss the potential for correcting activity signals in the RVs of M dwarfs.
(abridged) Spatially resolved observations of molecular line emission provide unique constraints on protoplanetary disk turbulence. Using local non-ideal MHD simulations and radiative transfer calculations, we assess the ability of ALMA observations to robustly detect and characterize disk turbulence. We specifically predict the outcome of the magnetorotational instability (MRI) in the disk around HD 163296, a promising observational target. We find that the MRI can support the observed level of accretion if the outer disk surface is ionized by far-UV photons and threaded by a weak net vertical magnetic field. We identify two classes of MRI solution - dynamo solutions in which the surface magnetic field reverses periodically, and non-dynamo solutions in which much of the Maxwell stress is steady and large scale. In both classes the small-scale turbulence increases in strength with height above the mid-plane, and can be represented as a microturbulent component. Using vertical profiles of the turbulent velocity from simulations at different radii, we use radiative transfer calculations to quantify the observational signatures. We show that the peak to line center flux ratio is a robust diagnostic of turbulence that is only mildly degenerate with uncertainties in disk temperature. For the CO(3-2) line variations in the predicted peak-to-trough ratio between our most and least turbulent models are ~15%. We develop predictions for other molecular lines and for channel maps whose morphology allows for independent constraints on turbulence.
We investigate the composition of interstellar grains along the line of sight toward Zeta Ophiuchi, a well-studied environment near the diffuse-dense cloud transition. A spectral decomposition analysis of the solid-state absorbers is performed using archival spectroscopic observations from the Spitzer Space Telescope and Infrared Space Observatory. We find strong evidence for the presence of sub-micron-sized amorphous silicate grains, principally comprised of olivine-like composition, with no convincing evidence of H2O ice mantles. However, tentative evidence for thick H2O ice mantles on large (a ~ 2.8 microns) grains is presented. Solid-state abundances of elemental Mg, Si, Fe, and O are inferred from our analysis and compared to standard reference abundances. We find that nearly all of elemental Mg and Si along the line of sight are present in amorphous silicate grains, while a substantial fraction of elemental Fe resides in compounds other than silicates. Moreover, we find that the total abundance of elemental O is largely inconsistent with the adopted reference abundances, indicating that as much as ~156 ppm of interstellar O is missing along the line of sight. After taking into account additional limits on the abundance of elemental O in other O-bearing solids, we conclude that any missing reservoir of elemental O must reside on large grains that are nearly opaque to infrared radiation.
Very-long baseline interferometric observations have resolved structure on scales of only a few Schwarzschild radii around the supermassive black holes at the centers of our Galaxy and M87. In the near future, such observations are expected to image the shadows of these black holes together with a bright and narrow ring surrounding their shadows. For a Kerr black hole, the shape of this photon ring is nearly circular unless the black hole spins very rapidly. Whether or not, however, astrophysical black holes are truly described by the Kerr metric as encapsulated in the no-hair theorem still remains an untested assumption. For black holes that differ from Kerr black holes, photon rings have been shown numerically to be asymmetric for small to intermediate spins. In this paper, I calculate semi-analytic expressions of the shapes of photon rings around black holes described by a new Kerr-like metric which is valid for all spins. I show that photon rings in this spacetime are affected by two types of deviations from the Kerr metric which can cause the ring shape to be highly asymmetric. I argue that the ring asymmetry is a direct measure of a potential violation of the no-hair theorem and that both types of deviations can be detected independently if the mass and distance of the black hole are known. In addition, I obtain approximate expressions of the diameters, displacements, and asymmetries of photon rings around Kerr and Kerr-like black holes.
The spins of a number of supermassive and stellar-mass black holes have been measured based on detections of thermal continuum emission and relativistically broadened iron lines in their x-ray spectra. Likewise, quasiperiodic variability has been observed in several sources. Such measurements commonly make the assumption that black holes are described by the Kerr metric, which according to the no-hair theorem characterizes black holes uniquely in terms of their masses and spins. This fundamental property of black holes can be tested observationally by measuring potential deviations from the Kerr metric introduced by a parametrically deformed Kerr-like spacetime. Thermal spectra, iron lines, and variability have already been studied extensively in several such metrics, which usually depend on only one particular type of deviation or contain unphysical regions outside of the compact object. In this paper, I study these x-ray probes in the background of a new Kerr-like metric which depends on four independent deviation functions and is free of pathological regions outside of the event horizon. I show that the observed signals depend significantly on primarily two types of deviations and that the strong correlation between the spin and the deviation parameters found previously in other Kerr-like metrics is partially broken for rapidly spinning black holes. This suggests that high-spin sources are the best candidates for tests of the no-hair theorem with x-rays and I obtain first constraints on such deviations from the stellar-mass black hole Cygnus X-1.
There is mounting observational evidence that most galactic nuclei host both supermassive black holes (SMBHs) and young populations of stars. With an abundance of massive stars, core-collapse supernovae are expected in SMBH spheres of influence. We develop a novel numerical method, based on the Kompaneets approximation, to trace supernova remnant (SNR) evolution in these hostile environments, where radial gas gradients and SMBH tides are present. We trace the adiabatic evolution of the SNR shock until 50% of the remnant is either in the radiative phase or is slowed down below the SMBH Keplerian velocity and is sheared apart. In this way, we obtain shapes and lifetimes of SNRs as a function of the explosion distance from the SMBH, the gas density profile and the SMBH mass. As an application, we focus here exclusively on quiescent SMBHs, because their light may not hamper detections of SNRs and because we can take advantage of the unsurpassed detailed observations of our Galactic Centre. Assuming that properties such as gas and stellar content scale appropriately with the SMBH mass, we study SNR evolution around other quiescent SMBHs. We find that, for SMBH masses over ~10^7 Msun, tidal disruption of SNRs can occur at less than 10^4 yr, leading to a shortened X-ray emitting adiabatic phase, and to no radiative phase. On the other hand, only modest disruption is expected in our Galactic Centre for SNRs in their X-ray stage. This is in accordance with estimates of the lifetime of the Sgr A East SNR, which leads us to expect one supernova per 10^4 yr in the sphere of influence of Sgr A*.
PTF11iqb was initially classified as a TypeIIn event caught very early after explosion. It showed narrow Wolf-Rayet (WR) spectral features on day 2, but the narrow emission weakened quickly and the spectrum morphed to resemble those of Types II-L and II-P. At late times, Halpha emission exhibited a complex, multipeaked profile reminiscent of SN1998S. In terms of spectroscopic evolution, we find that PTF11iqb was a near twin of SN~1998S, although with weaker interaction with circumstellar material (CSM) at early times, and stronger CSM interaction at late times. We interpret the spectral changes as caused by early interaction with asymmetric CSM that is quickly (by day 20) enveloped by the expanding SN ejecta photosphere, but then revealed again after the end of the plateau when the photosphere recedes. The light curve can be matched with a simple model for weak CSM interaction added to the light curve of a normal SN~II-P. This plateau requires that the progenitor had an extended H envelope like a red supergiant, consistent with the slow progenitor wind speed indicated by narrow emission. The cool supergiant progenitor is significant because PTF11iqb showed WR features in its early spectrum --- meaning that the presence of such WR features in an early SN spectrum does not necessarily indicate a WR-like progenitor. [abridged] Overall, PTF11iqb bridges SNe~IIn with weaker pre-SN mass loss seen in SNe II-L and II-P, implying a continuum between these types.
We compare model results from our semi-analytic merger tree based framework for high-redshift ($z \simeq 5-20$) galaxy formation against reionization indicators including the Planck electron scattering optical depth ($\tau_{es}$) and the ionizing photon emissivity ($\dot n_{ion}$) to constrain the particle mass of Warm Dark Matter (WDM). Our framework traces the Dark Matter (DM) and baryonic assembly of galaxies in 4 DM cosmologies: Cold Dark Matter (CDM) and WDM with a particle mass of $m_x = 2.25,3$ and 5 keV. It includes all the key processes of star formation, supernova feedback, the merger/accretion/ejection driven evolution of gas and stellar mass, and the effect of the ultra-violet background (UVB) created during reionization in photo-evaporating the gas content of galaxies in halos with $M_h \leq 10^9 M_\odot$. We show that current Planck $\tau_{es}$ values rule out $m_x \leq 2.5$ keV WDM, even in the physically unlikely scenario that all ionizing photons produced by these galaxies escape and contribute to reionization (i.e. $f_{esc}=1$). With the largest number of UVB-suppressed galaxies, CDM faces a "stalling" of the reionization process with this effect decreasing with the disappearance of small-scale structure with decreasing $m_x$. Finally, we find the bulk of the reionization photons come from galaxies with a halo mass $M_h \leq 10^9M_\odot$, stellar mass $M_* \leq 10^7M_\odot$ and UV magnitude $ -18 \leq M_{UV} \leq -13$ in CDM. The progressive suppression of low-mass halos with decreasing $m_x$ leads to a shift in the "reionization" population to larger (halo and stellar) masses of $M_h \geq 10^9M_\odot$ and $M_* \geq 10^7M_\odot$ for $m_x \geq 3$ keV WDM, although the UV limits effectively remain unchanged.
We model the pulse profiles and the phase resolved spectra of the anomalous X-ray pulsar 1E 1048.1-5937 obtained with XMM-Newton to map its surface temperature distribution during an active and a quiescent epoch. We develop and apply a model that takes into account the relevant physical and geometrical effects on the neutron star surface, magnetosphere, and spacetime. Using this model, we determine the observables at infinity as a function of pulse phase for different numbers and sizes of hot spots on the surface. We show that the pulse profiles extracted from both observations can be modeled with a single hot spot and an antipodal cool component. The size of the hot spot changes from $\approx 80^{\circ}$ in 2007, 3 months after the onset of a dramatic flux increase, to $\approx 30^{\circ}$ during the quiescent observation in 2011, when the pulsed fraction returned to the pre-outburst $\approx$ 65\% level. For the 2007 observation, we also find that a model consisting of a single 0.4 keV hot spot with a magnetic field strength of $1.8 \times 10^{14}$ G accounts for the spectra obtained at three different pulse phases but under predicts the flux at the pulse minimum, where the contribution to the emission from the cooler component is non-negligible. The inferred temperature of the spot stays approximately constant between different pulse phases, in agreement with a uniform temperature, single hot spot model. These results suggest that the emitting area grows significantly during outbursts but returns to its persistent and significantly smaller size within a few year timescale.
The use of Type~Ia SNe has thus far produced the most reliable measurement of the expansion history of the Universe, suggesting that $\Lambda$CDM offers the best explanation for the redshift--luminosity distribution observed in these events. But the analysis of other kinds of source, such as cosmic chronometers, gamma ray bursts, and high-$z$ quasars, conflicts with this conclusion, indicating instead that the constant expansion rate implied by the $R_{\rm h}=ct$ Universe is a better fit to the data. The central difficulty with the use of Type~Ia SNe as standard candles is that one must optimize three or four nuisance parameters characterizing supernova luminosities simultaneously with the parameters of an expansion model. Hence in comparing competing models, one must reduce the data independently for each. We carry~out such a comparison of $\Lambda$CDM and the $R_{\rm h}=ct$ Universe, using the Supernova Legacy Survey (SNLS) sample of 252 SN~events, and show that each model fits its individually reduced data very well. But since $R_{\rm h}=ct$ has only one free parameter (the Hubble constant), it follows from a standard model selection technique that it is to be preferred over $\Lambda$CDM, the minimalist version of which has three (the Hubble constant, the scaled matter density and either the spatial curvature constant or the dark-energy equation-of-state parameter). We estimate by the Bayes Information Criterion that in a pairwise comparison, the likelihood of $R_{\rm h}=ct$ is $\sim 90\%$, compared with only $\sim 10\%$ for a minimalist form of $\Lambda$CDM, in which dark energy is simply a cosmological constant. Compared to $R_{\rm h}=ct$, versions of the standard model with more elaborate parametrizations of dark energy are judged to be even less likely.
We report the blind detection of 12CO emission from a Distant Red Galaxy, HS1700.DRG55. We have used the IRAM PdBI-WIDEX, with its 3.6GHz of instantaneous dual-polarization bandwidth, to target 12CO(3--2) from galaxies lying in the proto-cluster at z=2.300 in the field HS1700+64. If indeed this line in DRG55 is 12CO(3--2), it's detection at 104.9GHz indicates a z_CO=2.296. None of the other eight known z~2.30 proto-cluster galaxies lying within the primary beam (PB) are detected in 12CO, although the limits are ~2x worse towards the edge of the PB where several lie. The optical/near-IR magnitudes of DRG55 (R_AB>27, K_AB=22.3) mean that optical spectroscopic redshifts are difficult with 10m-class telescopes, but near-IR redshifts would be feasible. The 24um-implied SFR (210 M_odot yr-1), stellar mass (~10^11 M-odot) and 12CO line luminosity (3.6x10^10 K km s-1 pc^2) are comparable to other normal 12CO-detected star forming galaxies in the literature, although the galaxy is some ~2 mag (~6x) fainter in the rest-frame UV than 12CO-detected galaxies at z>2. The detection of DRG55 in 12CO complements three other 12CO detected UV-bright galaxies in this proto-cluster from previous studies, and suggests that many optically faint galaxies in the proto-cluster may host substantial molecular gas reservoirs, and a full blind census of 12CO in this overdense environment is warranted.
We investigate whether the large scale structure environment of galaxy clusters imprints a selection bias on Sunyaev Zel'dovich (SZ) catalogs. Such a selection effect might be caused by line of sight (LoS) structures that add to the SZ signal or contain point sources that disturb the signal extraction in the SZ survey. We use the Planck PSZ1 union catalog (Planck Collab- oration et al. 2013a) in the SDSS region as our sample of SZ selected clusters. We calculate the angular two-point correlation function (2pcf) for physically correlated, foreground and background structure in the RedMaPPer SDSS DR8 catalog with respect to each cluster. We compare our results with an optically selected comparison cluster sample and with theoretical predictions. In contrast to the hypothesis of no environment-based selection, we find a mean 2pcf for background structures of -0.049 on scales of $\lesssim 40'$, significantly non-zero at $\sim 4 \sigma$, which means that Planck clusters are more likely to be detected in regions of low background density. We hypothesize this effect arises either from background estimation in the SZ survey or from radio sources in the background. We estimate the defect in SZ signal caused by this effect to be negligibly small, of the order of $\sim 10^{-4}$ of the signal of a typical Planck detection. Analogously, there are no implications on X-ray mass measurements. However, the environ- mental dependence has important consequences for weak lensing follow up of Planck galaxy clusters: we predict that projection effects account for half of the mass contained within a 15' radius of Planck galaxy clusters. We did not detect a background underdensity of CMASS LRGs, which also leaves a spatially varying redshift dependence of the Planck SZ selection function as a possible cause for our findings.
We report on a method, PUSH, for triggering core-collapse supernova explosions of massive stars in spherical symmetry. We explore basic explosion properties and calibrate PUSH such that the observables of SN1987A are reproduced. Our simulations are based on the general relativistic hydrodynamics code AGILE combined with the detailed neutrino transport scheme IDSA for electron neutrinos and ALS for the muon and tau neutrinos. To trigger explosions in the otherwise non-exploding simulations, we rely on the neutrino-driven mechanism. The PUSH method locally increases the energy deposition in the gain region through energy deposition by the heavy neutrino flavors. Our setup allows us to model the explosion for several seconds after core bounce. We explore the progenitor range 18-21M$_{\odot}$. Our studies reveal a distinction between high compactness (HC) and low compactness (LC) progenitor models, where LC models tend to explore earlier, with a lower explosion energy, and with a lower remnant mass. HC models are needed to obtain explosion energies around 1 Bethe, as observed for SN1987A. However, all the models with sufficiently high explosion energy overproduce $^{56}$Ni. We conclude that fallback is needed to reproduce the observed nucleosynthesis yields. The nucleosynthesis yields of $^{57-58}$Ni depend sensitively on the electron fraction and on the location of the mass cut with respect to the initial shell structure of the progenitor star. We identify a progenitor and a suitable set of PUSH parameters that fit the explosion properties of SN1987A when assuming 0.1M$_{\odot}$ of fallback. We predict a neutron star with a gravitational mass of 1.50M$_{\odot}$. We find correlations between explosion properties and the compactness of the progenitor model in the explored progenitors. However, a more complete analysis will require the exploration of a larger set of progenitors with PUSH.
Results are presented from a 3-D Pluto general circulation model (PGCM) that includes a subsurface model and volatile cycle. Conductive heating and cooling are present, as is non-local thermodynamic equilibrium (non-LTE) heating by methane at 2.3 and 3.3 microns, non-LTE cooling by heating by methane at 7.6 microns, and LTE CO rotational line cooling. This model is novel in that it has both detailed subsurface and atmospheric model components. Yet, there is little dependence of the model results on surface albedo, emissivity, or conductivity. Predictions are also provided for the Alice and REX instruments on New Horizons and for ground-based stellar occultations. Due to the weak temperature gradients, Alice (both solar and background star consultations) and REX are predicted to observe nearly the same temperature profiles on immersion and emersion. In the stratosphere, differences of up to 20 K are possible, while at higher altitudes (100-350 km), the differences are as large as 10 K. For both methane concentration and initial surface pressure, it should be possible to distinguish between the 0.2 and 1.0 methane concentrations and 8 and 24 microbar initial surface concentrations used here. For the ground-based stellar occultation, there is a detectable difference between light curves with the different methane concentrations used here, but not for the initial surface pressures.
PSR B1821$-$24 is a solitary millisecond pulsar (MSP) which radiates multi-wavelength pulsed photons. It has complex radio, X-ray and $\gamma$-ray pulse profiles with distinct peak phase-separations that challenge the traditional caustic emission models. Using the single-pole annular gap model with suitable magnetic inclination angle ($\alpha=40^\circ$) and viewing angle ($\zeta=75^\circ$), we managed to reproduce its pulse profiles of three wavebands. It is found that the middle radio peak is originated from the core gap region at high altitudes, and the other two radio peaks are originated from the annular gap region at relatively low altitudes. Two peaks of both X-ray and $\gamma$-ray wavebands are fundamentally originated from annular gap region, while the $\gamma$-ray emission generated from the core gap region contributes somewhat to the first $\gamma$-ray peak. Precisely reproducing the multi-wavelength pulse profiles of PSR B1821$-$24 enables us to understand emission regions of distinct wavebands and justify pulsar emission models.
Aims. We use advanced 3D NLTE radiative magnetohydrodynamic simulations of the solar atmosphere to carry out detailed tests of chromospheric diagnostics at millimeter and submillimeter wavelengths. Methods. We focused on the diagnostics of the thermal structure of the chromosphere in the wavelength bands from 0.4 mm up to 9.6 mm that can be accessed with the Atacama Large Millimeter/Submillimeter Array (ALMA) and investigated how these diagnostics are affected by the instrumental resolution. Results. We find that the formation height range of the millimeter radiation depends on the location in the simulation domain and is related to the underlying magnetic structure. Nonetheless, the brightness temperature is a reasonable measure of the gas temperature at the effective formation height at a given location on the solar surface. There is considerable scatter in this relationship, but this is significantly reduced when very weak magnetic fields are avoided. Our results indicate that although instrumental smearing reduces the correlation between brightness and temperature, millimeter brightness can still be used to reliably diagnose electron temperature up to a resolution of 1". If the resolution is more degraded, then the value of the diagnostic diminishes rapidly. Conclusions. We conclude that millimeter brightness can image the chromospheric thermal structure at the height at which the radiation is formed. Thus multiwavelength observations with ALMA with a narrow step in wavelength should provide sufficient information for a tomographic imaging of the chromosphere.
We performed 12CO(J=1-0) (hereafter, CO) observations towards 12 normal star-forming galaxies with stellar mass of Mstar=10^10.6-10^11.3 Msun at z=0.1-0.2 with the 45-m telescope at the Nobeyama Radio Observatory (NRO). The samples are selected with Dn(4000) that is a strength of the 4000 \AA break, instead of commonly used far-infrared (FIR) flux. We successfully detect the CO emissions from eight galaxies with signal-to-noise ratio (S/N) larger than three, demonstrating the effectiveness of the Dn(4000)-based sample selection. For the first time, we find a tight anti-correlation between Dn(4000) and molecular gas fraction (fmol) using literature data of nearby galaxies in which the galaxies with more fuel for star formation have younger stellar populations. We find that our CO-detected galaxies at z~0.1-0.2 also follow the same relation of nearby galaxies. This implies that the galaxies evolve along this Dn(4000)-fmol relation, and that Dn(4000) seems to be used as a proxy for fmol which requires many time-consuming observations. Based on the comparison with the model calculation with a population synthesis code, we find that star formation from metal enriched gas and its quenching in the early time are necessary to reproduce galaxies with large Dn(4000) and non-zero gas fraction.
The redshifted 21cm line of neutral hydrogen (HI), potentially observable at low radio frequencies (~50-200 MHz), should be a powerful probe of the physical conditions of the inter-galactic medium during Cosmic Dawn and the Epoch of Reionisation (EoR). The sky-averaged HI signal is expected to be extremely weak (~100 mK) in comparison to the foreground of up to 10000 K at the lowest frequencies of interest. The detection of such a weak signal requires an extremely stable, well characterised system and a good understanding of the foregrounds. Development of a nearly perfectly (~mK accuracy) calibrated total power radiometer system is essential for this type of experiment. We present the BIGHORNS (Broadband Instrument for Global HydrOgen ReioNisation Signal) experiment which was designed and built to detect the sky-averaged HI signal from the EoR at low radio frequencies. The BIGHORNS system is a mobile total power radiometer, which can be deployed in any remote location in order to collect radio-interference (RFI) free data. The system was deployed in remote, radio quiet locations in Western Australia and low RFI sky data have been collected. We present a description of the system, its characteristics, details of data analysis and calibration. We have identified multiple challenges to achieving the required measurement precision, which triggered two major improvements for the future system.
We present the first extensive study of the coronal line variability in an active galaxy. Our data set for the nearby source NGC 4151 consists of six epochs of quasi-simultaneous optical and near-infrared spectroscopy spanning a period of about eight years and five epochs of X-ray spectroscopy overlapping in time with it. None of the coronal lines showed the variability behaviour observed for the broad emission lines and hot dust emission. In general, the coronal lines varied only weakly, if at all. Using the optical [Fe VII] and X-ray O VII emission lines we estimate that the coronal line gas has a relatively low density of n~10^3 cm^-3 and a relatively high ionisation parameter of log U~1. The resultant distance of the coronal line gas from the ionising source is about two light years, which puts this region well beyond the hot inner face of the obscuring dusty torus. The high ionisation parameter implies that the coronal line region is an independent entity rather than part of a continuous gas distribution connecting the broad and narrow emission line regions. We present tentative evidence for the X-ray heated wind scenario of Pier & Voit. We find that the increased ionising radiation that heats the dusty torus also increases the cooling efficiency of the coronal line gas, most likely due to a stronger adiabatic expansion.
With references to both key and oft-forgotten pioneering works, this article starts by presenting a review into how we came to believe in the existence of massive black holes at the centres of galaxies. It then presents the historical development of the near-linear (black hole)-(host spheroid) mass relation, before explaining why this has recently been dramatically revised. Past disagreement over the slope of the (black hole)-(velocity dispersion) relation is also explained, and the discovery of sub-structure within the (black hole)-(velocity dispersion) diagram is discussed. As the search for the fundamental connection between massive black holes and their host galaxies continues, the competing array of additional black hole mass scaling relations for samples of predominantly inactive galaxies are presented.
We use Hubble Space Telescope (HST) to reach the end of the white dwarf (WD) cooling sequence (CS) in the solar-metallicity open cluster NGC 6819. Our photometry and completeness tests show a sharp drop in the number of WDs along the CS at magnitudes fainter than mF606W = 26.050+/- 0.075. This implies an age of 2.25+/-0.20 Gyr, consistent with the age of 2.25+/-0.30 Gyr obtained from fits to the main-sequence turn-off. The use of different WD cooling models and initial-final-mass relations have a minor impact the WD age estimate, at the level of ~0.1 Gyr. As an important by-product of this investigation we also release, in electronic format, both the catalogue of all the detected sources and the atlases of the region (in two filters). Indeed, this patch of sky studied by HST (of size ~70 arcmin sq.) is entirely within the main Kepler-mission field, so the high-resolution images and deep catalogues will be particularly useful.
The spatial and temporal invariance in the spectra of energetic particles in the gradual solar events is reproduced in the simulations. Based on a numerical solution of the focused transport equation, we obtain the intensity time profiles of solar energetic particles (SEPs) accelerated by an interplanetary shock in the three-dimensional interplanetary space. The shock is treated as a moving source of energetic particles with a distribution function. The time profiles of particle flux with different energies are calculated in the ecliptic at $1$ AU. We find that the spatial and temporal invariance in SEP spectra are the results of the effects of perpendicular diffusion and adiabatic cooling in the interplanetary space in our model. Furthermore, a spectra invariant region, which agrees with observations but is different than the one suggested by Reames and co-workers, is proposed based on our simulations.
We explored the impact of the synergy between the Euclid near-infrared
photometric surveys and the SKA radio continuum surveys on the studies of the
cosmic star formation. The Euclid satellite is expected to perform a Wide and
Deep photometric surveys to an infrared limit of H ~ 24 and H ~ 26 respectively
and a spectroscopy survey with a flux limit of $\sim 3 \times 10^{-16}$ erg
cm$^{-2}$ s$^{-1}$ in the Halpha line. Combining the H band Euclid selected
samples with the ground based ancillary data (fundamental for the SFR
estimation) we will be able to detect the star forming galaxies down to SFRs of
order of unit to z ~ 2 and down to SFR ~ 10 to z ~ 3, sampling the majority of
the star forming galaxies up to z ~3 and beyond and placing definitive
constraints on the star formation history of the universe at z<4-5 (is there a
peak a z ~2 or a plateau at 1 <z <5 ?) and on the galaxies evolution models.
The only tools able to provide a accurate dust-free calculation of their SFR
are the SKA continuum surveys.
The observational parameters of the Deep Tier SKA1 reference survey (a 0.2-
0.5 arcsec resolution and a 5 sigma detection limit of 1 microJy over 30 deg2
at Band 1/2 ) are the perfect complement of the Euclid survey. We showed, in
fact, that with this kind of SKA survey we will be able to determine a dust
unbiased SFR for a huge fraction (~85 %) of the Euclid SFG providing strong
constraints on the star formation history of the Universe.
By employing a simple semi-analytical star formation model where the formation rates of Population (Pop) I/II and III stars can be calculated, respectively, we account for the number distribution of gamma-ray bursts (GRBs) with high pseudo-redshifts that was derived from an empirical luminosity-indictor relationship. It is suggested that a considerable number of Pop III GRBs could exist in the present sample of Swift GRBs. By further combining the implication for the star formation history from the optical depth of the CMB photons, it is also suggested that only a very small fraction 0.6% of Pop III GRBs could have triggered the Swift BAT. These results could provide an useful basis for estimating future detectability of Pop III stars and their produced transient phenomena.
We represent a method to reconstruct the equation of state for dark energy directly from observational Hubble parameter data in a nonparametric way. We use principal component analysis (PCA) to extract the signal from data with noise. In addition, we modify Akaike information criteria (AIC) to guarantee the quality of reconstruction and avoid over-fitting simultaneously. The results show that our method is robust in reconstruction of dark energy equation of state. Although current observational Hubble parameter data alone can not give a strong constraint yet, future observations with more accurate data can help to improve the quality of reconstruction significantly, which is consistent with the results of H.-R. Yu et al.
(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $\mu$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $\mu$m, [OI] 63 $\mu$m, and [OI] 145 $\mu$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $\mu$m, and [CI] 609 $\mu$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $\mu$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $\mu$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.
We present the new code ALCAR developed to model multidimensional, multi energy-group neutrino transport in the context of supernovae and neutron-star mergers. The algorithm solves the evolution equations of the 0th- and 1st-order angular moments of the specific intensity, supplemented by an algebraic relation for the 2nd-moment tensor to close the system. The scheme takes into account frame-dependent effects of order O(v/c) as well as the most important types of neutrino interactions. The transport scheme is significantly more efficient than a multidimensional solver of the Boltzmann equation, while it is more accurate and consistent than the flux-limited diffusion method. The finite-volume discretization of the essentially hyperbolic system of moment equations employs methods well-known from hydrodynamics. For the time integration of the potentially stiff moment equations we employ a scheme in which only the local source terms are treated implicitly, while the advection terms are kept explicit, thereby allowing for an efficient computational parallelization of the algorithm. We investigate various problem setups in one and two dimensions to verify the implementation and to test the quality of the algebraic closure scheme. In our most detailed test, we compare a fully dynamic, one-dimensional core-collapse simulation with two published calculations performed with well-known Boltzmann-type neutrino-hydrodynamics codes and we find very satisfactory agreement.
Galaxy intrinsic alignment can be a severe source of error in weak-lensing studies. The problem has been widely studied by numerical simulations and with heuristic models, but without a clear theoretical justification of its origin and amplitude. In particular, it is still unclear whether intrinsic alignment of galaxies is dominated by formation and accretion processes or by the effects of the instantaneous tidal field acting upon them. We investigate this question by developing a simple model of intrinsic alignment for elliptical galaxies, based on the instantaneous tidal field. Making use of the galaxy stellar distribution function, we estimate the intrinsic alignment signal and find that although it has the expected dependence on the tidal field, it is too weak to account for the observed signal. This is an indirect validation of the standard view that intrinsic alignment is caused by formation and/or accretion processes.
We present new optical observations of the HH 111 Herbig-Haro jet using the Gemini Multi Object Spectrograph in its Integral Field Unit mode. Eight fields of 5" x 3.5" have been positioned along and across the HH 111 jet, covering the spatial region from knot E to L in HH 111 (namely, knots E, F, G, H, J, K and L). We present images and velocity channel maps for the [O I] 6300+6360, Halpha, [N II] 6548+6583 and [S II] 6716+6730 lines, as well as for the [S II]6716/6730 line ratio. We find that the HH 111 jet has an inner region with lower excitation and higher radial velocity, surrounded by a broader region of higher excitation and lower radial velocity. Also, we find higher electron densities at lower radial velocities. These results imply that the HH 111 jet has a fast, axial region with lower velocity shocks surrounded by a lower velocity sheath with higher velocity shocks.
Clathrate hydrates are believed to play a significant role in various solar system environments, e.g. comets, and the surfaces and interiors of icy satellites, however the structural factors governing their formation and dissociation are poorly understood. We demonstrate the use of a high pressure gas cell, combined with variable temperature cooling and time-resolved data collection, to the in situ study of clathrate hydrates under conditions relevant to solar system environments. Clathrates formed and processed within the cell are monitored in situ using synchrotron X-ray powder diffraction and Raman spectroscopy. X-ray diffraction allows the formation of clathrate hydrates to be observed as CO2 gas is applied to ice formed within the cell. Complete conversion is obtained by annealing at temperatures just below the ice melting point. A subsequent rise in the quantity of clathrate is observed as the cell is thermally cycled. Four regions between 100-5000cm-1 are present in the Raman spectra that carry features characteristic of both ice and clathrate formation. This novel experimental arrangement is well suited to studying clathrate hydrates over a range of temperature (80-500K) and pressure (1-100bar) conditions and can be used with a variety of different gases and starting aqueous compositions. We propose the increase in clathrate formation observed during thermal cycling may be due to the formation of a quasi liquid-like phase that forms at temperatures below the ice melting point, but which allows easier formation of new clathrate cages, or the retention and delocalisation of previously formed clathrate structures, possibly as amorphous clathrate. The structural similarities between hexagonal ice, the quasi liquid-like phase, and crystalline CO2 hydrate mean that differences in the Raman spectrum are subtle; however, all features out to 5000cm-1 are diagnostic of clathrate structure.
We report on early results from a pilot program searching for metal-poor stars with LAMOST and follow-up high-resolution observation acquired with the MIKE spectrograph attached to the Magellan~II telescope. We performed detailed abundance analysis for eight objects with iron abundances [Fe/H] < -2.0, including five extremely metal-poor (EMP; [Fe/H] < -3.0) stars with two having [Fe/H] < -3.5. Among these objects, three are newly discovered EMP stars, one of which is confirmed for the first time with high-resolution spectral observations. Three program stars are regarded as carbon-enhanced metal-poor (CEMP) stars, including two stars with no enhancement in their neutron-capture elements, which thus possibly belong to the class of CEMP-no stars; one of these objects also exhibits significant enhancement in nitrogen, and is thus a potential carbon and nitrogen-enhanced metal-poor star. The [X/Fe] ratios of the sample stars generally agree with those reported in the literature for other metal-poor stars in the same [Fe/H] range. We also compared the abundance patterns of individual program stars with the average abundance pattern of metal-poor stars, and find only one chemically peculiar object with abundances of at least two elements (other than C and N) showing deviations larger than 0.5dex. The distribution of [Sr/Ba] versus [Ba/H] agrees that an additional nucleosynthesis mechanism is needed aside from a single r-process. Two program stars with extremely low abundances of Sr and Ba support the prospect that both main and weak r-process may have operated during the early phase of Galactic chemical evolution. The distribution of [C/N] shows that there are two groups of carbon-normal giants with different degrees of mixing. However, it is difficult to explain the observed behavior of the [C/N] of the nitrogen-enhanced unevolved stars based on current data.
Context. Galaxy clusters can be used as cosmological probes, but to this end, they need to be thoroughly understood. Combining all cluster observables in a consistent way will help us to understand their global properties and their internal structure. Aims. We provide proof of the concept that the projected gravitational potential of galaxy clusters can directly be reconstructed from X-ray observations. We also show that this joint analysis can be used to locally test the validity of the equilibrium assumptions in galaxy clusters. Methods. We used a newly developed reconstruction method, based on Richardson-Lucy deprojection, that allows reconstructing projected gravitational potentials of galaxy clusters directly from X-ray observations. We applied this algorithm to the well-studied cluster Abell 1689 and compared the gravitational potential reconstructed from X-ray observables to the potential obtained from gravitational lensing measurements. [...] Results. Assuming spherical symmetry and hydrostatic equilibrium, the potentials recovered from gravitational lensing and from X-ray emission agree very well beyond 500 kpc. Owing to the fact that the Richardson-Lucy deprojection algorithm allows deprojecting each line of sight independently, this result may indicate that non-gravitational effects and/or asphericity are strong in the central regions of the clusters. Conclusions. We demonstrate the robustness of the potential reconstruction method based on the Richardson-Lucy deprojection algorithm and show that gravitational lensing and X-ray emission lead to consistent gravitational potentials. Our results illustrate the power of combining galaxy-cluster observables in a single, non-parametric, joint reconstruction of consistent cluster potentials that can be used to locally constrain the physical state of the gas.
In this letter we present, for the first time, a study of star formation rate, gas fraction and galaxy morphology of a constrained simulation of the Milky Way (MW) and Andromeda (M31) galaxies, compared to other MW-mass galaxies. By combining with unconstrained simulations we cover a sufficient volume to compare these galaxies environmental densities ranging from the field to that of the Local Group (LG). This is particularly relevant as it has been shown that, quite generally, galaxy properties depend intimately upon their environment, most prominently when galaxies in clusters are compared to those in the field. For galaxies in loose groups such as the LG, however, environmental effects have been less clear. We consider the galaxy's environmental density in spheres of 1200 kpc (comoving) and find that whilst environment does not appear to directly affect morphology, there is a positive trend with star formation rates. This enhancement in star formation occurs systematically for galaxies in higher density environments, regardless whether they are part of the LG or in filaments. Our simulations suggest that the richer environment at Mpc-scales may help replenish the star-forming gas, allowing higher specific star formation rates in galaxies such as the MW.
The solar system started to form about 4.56 Gyr ago and despite the long intervening time span, there still exist several clues about its formation. The three major sources for this information are meteorites, the present solar system structure and the planet-forming systems around young stars. In this introduction we give an overview of the current understanding of the solar system formation from all these different research fields. This includes the question of the lifetime of the solar protoplanetary disc, the different stages of planet formation, their duration, and their relative importance. We consider whether meteorite evidence and observations of protoplanetary discs point in the same direction. This will tell us whether our solar system had a typical formation history or an exceptional one. There are also many indications that the solar system formed as part of a star cluster. Here we examine the types of cluster the Sun could have formed in, especially whether its stellar density was at any stage high enough to influence the properties of today's solar system. The likelihood of identifying siblings of the Sun is discussed. Finally, the possible dynamical evolution of the solar system since its formation and its future are considered.
The extremely high sensitivity and resolution of the Square Kilometre Array (SKA) will be useful for addressing a wide set of themes relevant for cosmology, in synergy with current and future cosmic microwave background (CMB) projects. Many of these themes also have a link with future optical-IR and X-ray observations. We discuss the scientific perspectives for these goals, the instrumental requirements and the observational and data analysis approaches, and identify several topics that are important for cosmology and astrophysics at different cosmic epochs.
Measurements of the neutral hydrogen gas content of a sample of 93 post-merger galaxies are presented, from a combination of matches to the ALFALFA.40 data release and new Arecibo observations. By imposing completeness thresholds identical to that of the ALFALFA survey, and by compiling a mass-, redshift- and environment-matched control sample from the public ALFALFA.40 data release, we calculate gas fraction offsets (Delta f_gas) for the post-mergers, relative to the control sample. We find that the post-mergers have HI gas fractions that are consistent with undisturbed galaxies. However, due to the relative gas richness of the ALFALFA.40 sample, from which we draw our control sample, our measurements of gas fraction enhancements are likely to be conservative lower limits. Combined with comparable gas fraction measurements by Fertig et al. in a sample of galaxy pairs, who also determine gas fraction offsets consistent with zero, we conclude that there is no evidence for significant neutral gas consumption throughout the merger sequence. From a suite of 75 binary merger simulations we confirm that star formation is expected to decrease the post-merger gas fraction by only 0.06 dex, even several Gyr after the merger. Moreover, in addition to the lack of evidence for gas consumption from gas fraction offsets, the observed HI detection fraction in the complete sample of post-mergers is twice as high as the controls, which suggests that the post-merger gas fractions may actually be enhanced. We demonstrate that a gas fraction enhancement in post-mergers, relative to a stellar mass-matched control sample, would indeed be the natural result of merging randomly drawn pairs from a parent population which exhibits a declining gas fraction with increasing stellar mass.
Galactic and extragalactic relativistic jets have rich environments that are full of moving objects, such as stars and dense clumps. These objects can enter into the jets and generate shocks and non-thermal emission. We characterize the emitting properties of the downstream region of a standing shock formed due to the interaction of a relativistic jet with an obstacle. We focus on the case of red giants interacting with an extragalactic jet. We perform relativistic axisymmetric hydrodynamical simulations of a relativistic jet meeting an obstacle of very large inertia. The results are interpreted in the framework of a red giant whose dense and slow wind interacts with the jet of an active galactic nucleus. Assuming that particles are accelerated in the standing shock generated in the jet as it impacts the red giant wind, we compute the non-thermal particle distribution, the Doppler boosting enhancement, and the non-thermal luminosity in gamma rays. The available non-thermal energy from jet-obstacle interactions is potentially enhanced by a factor of $\sim 100$ when accounting for the whole surface of the shock induced by the obstacle, instead of just the obstacle section. The observer gamma-ray luminosity, including the flow velocity and Doppler boosting effects, can be ~300(g/10)^2 times higher than when the emitting flow is assumed at rest and only the obstacle section is considered, where g is the jet Lorentz factor. For a whole population of red giants inside the jet of an AGN, the predicted persistent gamma-ray luminosities may be potentially detectable for a jet pointing to the observer. Obstacles interacting with relativistic outflows, for instance clouds and populations of stars for extragalactic jets, or stellar wind inhomogeneities in microquasar jets and in winds of pulsars in binaries, should be taken into account when investigating the non-thermal emission from these sources.
We present constraints on testing general relativity (GR) at cosmological scales using recent data sets and the impact of galaxy intrinsic alignment (IA) in the CFHTLenS lensing data on those constraints. We consider CMB temperature data from Planck, the galaxy power spectrum from WiggleZ, weak lensing tomography from the CFHTLenS, ISW-galaxy cross correlations, and BAO data from 6dF, SDSS DR7, and BOSS DR9. We use a parameterization of the modified gravity (MG) that is binned in redshift and scale, a parameterization that evolves monotonically in scale but is binned in redshift, and a functional parameterization that evolves only in redshift. We present the results in terms of the MG parameters $Q$ and $\Sigma$. We employ an IA model with an amplitude $A_{CFHTLenS}$ that is included in the parameter analysis. We find an improvement in the constraints on the MG parameters corresponding to $40-53\%$ increase on the figure of merit compared to previous studies, and GR is found consistent with the data at the $95\%$ CL. The bounds found on $A_{CFHTLenS}$ are sensitive to whether the MG parameterization is scale dependent, and the correlations between $A_{CFHTLenS}$ and MG parameters are found weak to moderate. $A_{CFHTLenS}$ is found consistent with zero for the 3 MG parameterizations and when the whole lensing sample is used. A significantly non-zero $A_{CFHTLenS}$ for GR and the scale-independent MG parameterization is found when we use the optimized early-type galaxy sample of (Heymans et al. 2013). We find that the tensions observed in previous studies persist, and there is an indication that CMB data and lensing data prefer different values for MG parameters, particularly for the parameter $\Sigma$. The analysis of the confidence contours and probability distributions suggest that the bimodality found follows that of the known tension in the $\sigma_8$ parameter. (Abridged)
We formulate within a generalized distributional approach the treatment of the stability against radial perturbations for both neutral and charged stratified stars in Newtonian and Einstein's gravity. We obtain from this approach the boundary conditions connecting two any phases within a star and underline its relevance for realistic models of compact stars with phase transitions, owing to the modification of the star's set of eigenmodes with respect to the continuous case.
According to the no-hair theorem, astrophysical black holes are uniquely characterized by their masses and spins and are described by the Kerr metric. Several parametric spacetimes which deviate from the Kerr metric have been proposed in order to test this theorem with observations of black holes in both the electromagnetic and gravitational-wave spectra. Such metrics often contain naked singularities or closed timelike curves in the vicinity of the compact objects that can limit the applicability of the metrics to compact objects that do not spin rapidly, and generally admit only two constants of motion. The existence of a third constant, however, can facilitate the calculation of observables, because the equations of motion can be written in first-order form. In this paper, I design a Kerr-like black hole metric which is regular everywhere outside of the event horizon, possesses three independent constants of motion, and depends nonlinearly on four free functions that parameterize potential deviations from the Kerr metric. This metric is generally not a solution to the field equations of any particular gravity theory, but can be mapped to known four-dimensional black hole solutions of modified theories of gravity for suitable choices of the deviation functions. I derive expressions for the energy, angular momentum, and epicyclic frequencies of a particle on a circular equatorial orbit around the black hole and compute the location of the innermost stable circular orbit. In addition, I write the metric in a Kerr-Schild-like form, which allows for a straightforward implementation of fully relativistic magnetohydrodynamic simulations of accretion flows in this metric. The properties of this metric make it a well-suited spacetime for strong-field tests of the no-hair theorem in the electromagnetic spectrum with black holes of arbitrary spin.
Higgs G-inflation takes advantage of a Galileon-like ghost-free derivative coupling. It is a nonrenormalizable operator and is strongly coupled at high energy scales. Perturbative analysis has no longer predictive power there. In general, when the Lagrangian is expanded around the vacuum, the strong coupling scale is identified as the mass scale that appears in nonrenormalizable operators. In inflationary models, however, the identification of the strong coupling scale is subtle, since the structures of the kinetic term as well as the interaction itself are modified by the background inflationary dynamics. As a result, the strong coupling scale is back ground field dependent. In this letter, we evaluate the strong coupling scale of the fluctuations around the inflationary background including the Nambu Goldstone mode associated with the symmetry breaking in the Higgs G-inflation. We find that the system is weakly coupled when the scales which we now observe exit the horizon during inflation, and the observational predictions with the semiclassical treatment are valid. However, we also find that the inflaton field value where the strong coupling scale and the Hubble scale meet is less than the Planck scale. Therefore, we cannot describe the model from the Planck scale, or the chaotic initial condition.
We study the $f$-mode frequencies and damping times of nonrotating neutron stars (NS) in general relativity (GR) by solving the linearized perturbation equations, with the aim to establish "universal" relations that depend only weakly on the equations of state (EOS). Using a more comprehensive set of EOSs, we re-examine some proposed linearizations that describe the $f$-mode parameters in terms of mass and radius of the neutron star (NS), and we test a more recent proposal for expressing the $f$-mode parameters as quadratic functions of the effective compactness. Our extensive results for each equation of state considered allow us to study the accuracy of each proposal. In particular, we find that the damping time deviates quite considerably from the proposed linearization. We introduce a new universal relation for the product of the $f$-mode frequency and damping time as a function of the (ordinary) compactness, which proved to be more accurate. The relations using the effective compactness on the other hand also fit our data accurately. Our results show that the maximum oscillation frequency depends strongly on the EOS, such that the measurement of a high oscillation frequency would rule out several EOSs. Lastly, we compare the exact mode frequencies to those obtained in the Cowling approximation, and also to results obtained with a nonlinear evolution code, validating the implementations of the different approaches.
We consider the evolution of electromagnetic fields coupled to conduction currents during the reheating era, under the assumption that the currents may be described by second order casual hydrodynamics. The resulting theory is not conformally invariant. The expansion of the Universe produces temperature gradients which couple to the current and generally oppose Ohmic dissipation. Although the effect is not strong, it suggests that the unfolding of hydrodynamic instabilities in these models may follow a different pattern than in first order theories, and even than in second order theories on non expanding backgrounds.
Interactions between dark matter and dark energy, allowing both conformal and and disformal couplings, are studied in detail. We discuss the background evolution, anisotropies in the cosmic microwave background and large scale structures. One of our main findings is that a large conformal coupling is not necessarily disallowed in the presence of a general disformal term. On the other hand, we find that negative disformal couplings very often lead to instabilities in the scalar field. Studying the background evolution and linear perturbations only, our results show that it is observationally challenging to disentangle disformal from purely conformal couplings.
The electromagnetic bremsstrahlung spectrum for the dipole which falling by a spiral orbit into the Schwarzschild black hole was found. The characteristic features in this electromagnetic spectrum can be used for determine of the black hole mass by the new way. This new way (if implemented) provides higher accuracy in determining of the black hole mass. Also these features in the spectrum can be used for determine of the certain characteristics in the black hole magnetosphere or in the accretion disk characteristics around the black hole. It is also shown that the asymptotic behavior of this spectrum (at high frequencies) is practically independent from the impact parameter of the falling dipole.
Links to: arXiv, form interface, find, astro-ph, recent, 1501, contact, help (Access key information)
We analyze 221 eclipsing binaries (EBs) in the Large Magellanic Cloud with B-type main-sequence (MS) primaries ($M_1$ $\approx$ 4 - 14 M$_{\odot}$) and orbital periods $P$ = 20 - 50 days that were photometrically monitored by the Optical Gravitational Lensing Experiment. We utilize our three-stage automated pipeline to (1) classify all 221 EBs, (2) fit physical models to the light curves of 130 detached well-defined EBs from which unique parameters can be determined, and (3) recover the intrinsic binary statistics by correcting for selection effects. We uncover two statistically significant trends with age. First, younger EBs tend to reside in dustier environments with larger photometric extinctions, an empirical relation that can be implemented when modeling stellar populations. Second, younger EBs generally have large eccentricities. This demonstrates that massive binaries at moderate orbital periods are born with a Maxwellian "thermal" orbital velocity distribution, which indicates they formed via dynamical capture. In addition, the age-eccentricity anticorrelation provides a direct constraint for tidal evolution in highly eccentric binaries containing hot MS stars with radiative envelopes. The intrinsic fraction of B-type MS stars with stellar companions $q$ = $M_2$/$M_1$ $>$ 0.2 and orbital periods $P$ = 20 - 50 days is (8 $\pm$ 2)%. We find early-type binaries at $P$ = 20 - 50 days are weighted significantly toward small mass ratios $q$ $\approx$ 0.2 - 0.3, which is different than previous observations of close binaries with $P$ $<$ 20 days. This indicates that early-type binaries at slightly wider orbital separations have experienced substantially less coevolution and competitive accretion during their formation in the circumbinary disk.
Age determination is undertaken for nearby early-type (BAF) stars, which constitute attractive targets for high-contrast debris disk and planet imaging surveys. Our analysis sequence consists of: acquisition of uvby{\beta} photometry from catalogs, correction for the effects of extinction, interpolation of the photometry onto model atmosphere grids from which atmospheric parameters are determined, and finally, comparison to the theoretical isochrones from pre-main sequence through post-main sequence stellar evolution models, accounting for the effects of stellar rotation. We calibrate and validate our methods at the atmospheric parameter stage by comparing our results to fundamentally determined Teff and log g values. We validate and test our methods at the evolutionary model stage by comparing our results on ages to the accepted ages of several benchmark open clusters (IC 2602, {\alpha} Persei, Pleiades, Hyades). Finally, we apply our methods to estimate stellar ages for 3493 field stars, including several with directly imaged exoplanet candidates.
We present an analysis of the galaxy-scale gaseous outflows from the FIRE (Feedback in Realistic Environments) simulations. This suite of hydrodynamic cosmological zoom simulations provides a sample of halos where star-forming giant molecular clouds are resolved to z=0, and features an explicit stellar feedback model on small scales. In this work, we focus on quantifying the gas mass ejected out of galaxies in winds and how this material travels through the halo. We correlate these quantities to star formation in galaxies throughout cosmic history. Our simulations reveal that a significant portion of every galaxy's evolution, particularly at high redshift, is dominated by bursts of star formation, which are followed by powerful gusts of galactic outflow that sweep up a large fraction of gas in the interstellar medium and send it through the circumgalactic medium. The dynamical effect of these outflows can significantly limit the amount of star formation within the affected galaxy. At low redshift, however, sufficiently massive galaxies corresponding to L*-progenitors develop stable disks and switch into a continuous and quiescent mode of star formation that does not drive outflows into the halo. We find inflow to be more continuous than outflow, although filamentary accretion onto the galaxy can be temporarily disrupted by recently ejected outflows. Using a variety of techniques, we measure outflow rates and use them to derive mass-loading factors, and their dependence on circular velocity, halo mass, and stellar mass for a large sample of galaxies in the FIRE simulation suite, spanning four decades in halo mass, six decades in stellar mass, and a redshift range of 4.0 > z > 0. Mass-loading factors for L*-progenitors are eta ~= 10 at high redshift, but decrease to eta << 1 at low redshift. [continued in text]
We derive [K/Fe] abundance ratios for 119 stars in the globular cluster NGC 2808, all of them having O, Na, Mg and Al abundances homogeneously measured in previous works. We detect an intrinsic star-to-star spread in the Potassium abundance. Moreover [K/Fe] abundance ratios display statistically significant correlations with [Na/Fe] and [Al/Fe], and anti-correlations with [O/Fe] and [Mg/Fe]. All the four Mg deficient stars ([Mg/Fe]<0.0) discovered so far in NGC 2808 are enriched in K by ~0.3 dex with respect to those with normal [Mg/Fe]. NGC 2808 is the second globular cluster, after NGC 2419, where a clear Mg-K anti-correlation is detected, albeit of weaker amplitude. The simultaneous correlation/anti-correlation of [K/Fe] with all the light elements usually involved in the chemical anomalies observed in globular cluster stars, strongly support the idea that these abundance patterns are due to the same self-enrichment mechanism that produces Na-O and Mg-Al anti-correlations. This finding suggests that detectable spreads in K abundances may be typical in the massive globular clusters where the self-enrichment processes are observed to produce their most extreme manifestations.
Low-cost mass-produced sensors and optics have recently made it feasible to build telescope arrays which observe the entire accessible sky simultaneously. In this article we discuss the scientific motivation for these telescopes, including exoplanets, stellar variability and extragalactic transients. To provide a concrete example we detail the goals and expectations for the Evryscope, an under-construction 780 MPix telescope which covers 8,660 square degrees in each two-minute exposure; each night, 18,400 square degrees will be continuously observed for an average of approximately 6 hours. Despite its small 61mm aperture, the system's large field of view provides an etendue which is ~10% of LSST. The Evryscope, which places 27 separate individual telescopes into a common mount which tracks the entire accessible sky with only one moving part, will return 1%-precision, many-year-length, high-cadence light curves for every accessible star brighter than mV=16.5, with brighter stars having few-millimagnitude photometric precision in long-term light curves. It will be capable of searching for transiting giant planets around the brightest and most nearby stars, where the planets are much easier to characterize; it will also search for small planets nearby M-dwarfs, for planetary occultations of white dwarfs, and will perform comprehensive nearby microlensing and eclipse-timing searches for exoplanets inaccessible to other planet-finding methods. The Evryscope will also monitor outbursting young stars, white dwarf activity, and stellar activity of all types, along with finding a large sample of very-long-period M-dwarf eclipsing binaries. When relatively rare transients events occur, such as gamma-ray bursts (GRBs), nearby supernovae, or even gravitational wave detections, the array will return minute-by-minute light curves without needing pointing towards the event as it occurs. (abridged)
We discuss the question of gauge choice when analysing second-order relativistic perturbations and when defining the galaxy bias at second order. Some misconceptions in the recent literature on the comoving-synchronous gauge are addressed and we show that this gauge is appropriate at second order to describe the matter overdensity and to define local Lagrangian bias.
The origin of ultra-compact dwarfs (UCDs)--objects larger and more massive than typical globular clusters (GCs), but more compact than typical dwarf galaxies--has been hotly debated in the 15 years since their discovery. Even whether UCDs should be considered galactic in origin, or simply the most extreme GCs, is not yet settled. We present the dynamical properties of 97 spectroscopically confirmed UCDs (rh >~10 pc) and 911 GCs associated with central cD galaxy of the Virgo cluster, M87. Our UCDs, of which 89% have M_star > ~2X10^6 M_sun and 92% are as blue as the classic blue GCs, nearly triple the sample of previous confirmed Virgo UCDs, providing by far the best opportunity for studying the global dynamics of a UCD system. We found that (1) UCDs have a surface number density profile that is shallower than that of the blue GCs in the inner ~ 70 kpc and as steep as that of the red GCs at larger radii; (2) UCDs exhibit a significantly stronger rotation than the GCs, and the blue GCs seem to have a velocity field that is more consistent with that of the surrounding dwarf ellipticals than with that of UCDs; (3) UCDs have a radially increasing orbital anisotropy profile, and are tangentially-biased at radii < ~ 40 kpc and radially-biased further out. In contrast, the blue GCs become more tangentially-biased at larger radii beyond ~ 40 kpc; (4) GCs with M_star > 2X10^6 M_sun have rotational properties indistinguishable from the less massive ones, suggesting that it is the size, instead of mass, that differentiates UCDs from GCs as kinematically distinct populations. We conclude that most UCDs in M87 are not consistent with being merely the most luminous and extended examples of otherwise normal GCs. The radially-biased orbital structure of UCDs at large radii is in general agreement with the "tidally threshed dwarf galaxy" scenario.
Massive young stellar objects (MYSOs) with hot cores are classic sources of
complex organic molecules. The origins of these molecules in such sources, as
well as the small- and large-scale differentiation between nitrogen- and
oxygen-bearing complex species, are poorly understood.
We aim to use complex molecule abundances toward a chemically less explored
class of MYSOs with weak hot organic emission lines to constrain the impact of
hot molecular cores and initial ice conditions on the chemical composition
toward MYSOs.
We use the IRAM 30m and the Submillimeter Array to search for complex organic
molecules over 8-16 GHz in the 1~mm atmospheric window toward three MYSOs with
known ice abundances, but without luminous molecular hot cores.
Complex molecules are detected toward all three sources at comparable
abundances with respect to CH$_3$OH to classical hot core sources. The relative
importance of CH$_3$CHO, CH$_3$CCH, CH$_3$OCH$_3$, CH$_3$CN, and HNCO differ
between the organic-poor MYSOs and hot cores, however. Furthermore, the
N-bearing molecules are generally concentrated toward the source centers, while
most O- and C-bearing molecules are present both in the center and in the
colder envelope. Gas-phase HNCO/CH$_3$OH ratios are tentatively correlated with
the ratios of NH$_3$ ice over CH$_3$OH ice in the same lines of sight, which is
consistent with new gas-grain model predictions.
Hot cores are not required to form complex organic molecules, and source
temperature and initial ice composition both seem to affect complex organic
distributions toward MYSOs. To quantify the relative impact of temperature and
initial conditions requires, however, a larger spatially resolved survey of
MYSOs with ice detections.
(Abridged) In the implicit large eddy simulation (ILES) paradigm, the dissipative nature of high-resolution shock-capturing schemes is exploited to provide an implicit model of turbulence. Recent 3D simulations suggest that turbulence might play a crucial role in core-collapse supernova explosions, however the fidelity with which turbulence is simulated in these studies is unclear. Especially considering that the accuracy of ILES for the regime of interest in CCSN, weakly compressible and strongly anisotropic, has not been systematically assessed before. In this paper we assess the accuracy of ILES using numerical methods most commonly employed in computational astrophysics by means of a number of local simulations of driven, weakly compressible, anisotropic turbulence. We report a detailed analysis of the way in which the turbulent cascade is influenced by the numerics. Our results suggest that anisotropy and compressibility in CCSN turbulence have little effect on the turbulent kinetic energy spectrum and a Kolmogorov $k^{-5/3}$ scaling is obtained in the inertial range. We find that, on the one hand, the kinetic energy dissipation rate at large scales is correctly captured even at relatively low resolutions, suggesting that very high effective Reynolds number can be achieved at the largest scales of the simulation. On the other hand, the dynamics at intermediate scales appears to be completely dominated by the so-called bottleneck effect, \ie the pile up of kinetic energy close to the dissipation range due to the partial suppression of the energy cascade by numerical viscosity. An inertial range is not recovered until the point where relatively high resolution $\sim 512^3$, which would be difficult to realize in global simulations, is reached. We discuss the consequences for CCSN simulations.
We make use of two suites of ultra-high resolution N-body simulations of individual dark matter haloes, the Phoenix and the Aquarius Projects, to investigate the systematics of subhalo assembly histories in host haloes differing by a factor of 1000 in mass. We find that the progenitors of the present day subhalo population are relatively more abundant in high mass haloes, in contrast to previous studies claiming a universal abundance independent of the mass of the host halo. This is mainly because these studies count progenitors that pass through the halo and are later re-accreted more than once. The fraction of these 'wavering' progenitors is larger in less massive haloes. The typical accretion time for all progenitors varies strongly with host halo mass: $z \sim 5$ for the Galactic-scale Aquarius haloes and $z \sim 2.5$ for the cluster-scale Phoenix haloes. Once progenitors start to orbit their parent haloes, they rapidly lose their original mass, but nevertheless more than 80 (70) percent of them survive to present day in the Phoenix (Aquarius) haloes. At given redshift, the fraction of subhaloes that survive is independent of the host halo mass, whilst the fraction of mass lost by subhaloes is larger in higher mass haloes. These systematics explain many similarities and differences between subhalo populations in haloes of different masses at the present day.
We present an expanded kinematic study of the young cluster NGC 2264 based
upon optical radial velocities measured using multi-fiber echelle spectroscopy
at the 6.5 meter MMT and Magellan telescopes. We report radial velocities for
695 stars, of which approximately 407 stars are confirmed or very likely
members. Our results more than double the number of members with radial
velocities from F{\H u}r{\'e}sz et al., resulting in a much better defined
kinematic relationship between the stellar population and the associated
molecular gas.
In particular, we find that there is a significant subset of stars that are
systematically blueshifted with respect to the molecular ($^{13}$CO) gas. The
detection of Lithium absorption and/or infrared excesses in this blue-shifted
population suggests that at least some of these stars are cluster members; we
suggest some speculative scenarios to explain their kinematics. Our results
also more clearly define the redshifted population of stars in the northern end
of the cluster; we suggest that the stellar and gas kinematics of this region
are the result of a bubble driven by the wind from O7 star S Mon. Our results
emphasize the complexity of the spatial and kinematic structure of NGC 2264,
important for eventually building up a comprehensive picture of cluster
formation.
Almost all galaxies along the Hubble sequence host a compact massive object (CMO) in their center. The CMO can be either a supermassive black hole (SMBH) or a very dense stellar cluster, also known as nuclear star cluster (NSC). Generally, heavier galaxies (mass >~ 10^{11} solar masses) host a central SMBH while lighter show a central NSC. Intermediate mass hosts, instead, contain both a NSC and a SMBH. One possible formation mechanisms of a NSC relies on the dry-merger (migratory) scenario, in which globular clusters (GCs) decay toward the center of the host galaxy and merge. In this framework, the absence of NSCs in high-mass galaxies can be imputed to destruction of the infalling GCs by the intense tidal field of the central SMBH. In this work, we report preliminary results of N-body simulations performed using our high-resolution, direct, code HiGPUs, to investigate the effects of a central SMBH on a single GC orbiting around it. By varying either the mass of the SMBH and the mass of the host galaxy, we derived an upper limit to the mass of the central SMBH, and thus to the mass of the host, above which the formation of a NSC is suppressed.
Constraining the properties of Population III (Pop III) stars will be very challenging because they reside in small galaxies at high redshift which will be difficult to detect even with future instruments such as the James Webb Space Telescope (JWST). In this paper, we suggest that intensity mapping may be a promising method to study Pop III stars. Intensity mapping is a technique proposed to measure large-scale fluctuations of galaxy line emission in three dimensions without resolving individual sources. This technique is well suited for observing many faint galaxies because it can measure their cumulative emission even if they cannot be directly detected. We focus on intensity mapping of He II recombination lines, and in particular He II 1640 \AA{}. These lines are much stronger in Pop III stars than Pop II stars because the harder spectra of Pop III stars are expected to produce many He II ionizing photons. Measuring the He II 1640 \AA{} intensity mapping signal, along with the signals from other lines such as Ly$\alpha$, H$\alpha$, and metal lines, could give constraints on the initial mass function (IMF) and star formation rate density (SFRD) of Pop III stars as a function of redshift. To demonstrate the feasibility of these observations, we estimate the strength of the Pop III He II 1640 \AA{} intensity mapping signal from $z=10-20$. We show that at $z\approx10$, the signal could be measured accurately by two different hypothetical future instruments, one which cross-correlates He II 1640 \AA{} with CO(1-0) line emission from galaxies and the other with 21 cm emission from the intergalactic medium (IGM).
M82 X-1, the brightest X-ray source in the galaxy M82, has been thought to be an intermediate-mass black hole (100 to 10,000 solar masses) because of its extremely high luminosity and variability characteristics, although some models suggest that its mass may be only about 20 solar masses. The previous mass estimates were based on scaling relations that use low-frequency characteristic timescales which have large intrinsic uncertainties. For stellar-mass black holes, we know that the high-frequency quasi-periodic oscillations (100-450 hertz) in the X-ray emission that occur in a 3:2 frequency ratio are stable and scale in frequency inversely with black hole mass with a reasonably small dispersion. The discovery of such stable oscillations thus potentially offers an alternative and less ambiguous means of mass determination for intermediate-mass black holes, but has hitherto not been realized. Here we report stable, twin-peak (3:2 frequency ratio) X-ray quasi-periodic oscillations from M82 X-1 at frequencies of 3.32$\pm$0.06 hertz and 5.07$\pm$0.06 hertz. Assuming that we can extrapolate the inverse-mass scaling that holds for stellar-mass black holes, we estimate the black hole mass of M82 X-1 to be 428$\pm$105 solar masses. In addition, we can estimate the mass using the relativistic precession model, from which we get a value of 415$\pm$63 solar masses.
We present deep HST/STIS coronagraphic images of the Beta Pic debris disk obtained at two epochs separated by 15 years. The new images and the re-reduction of the 1997 data provide the most sensitive and detailed views of the disk at optical wavelengths as well as the yet smallest inner working angle optical coronagraphic image of the disk. Our observations characterize the large-scale and inner-disk asymmetries and we identify multiple breaks in the disk radial surface brightness profile. We study in detail the radial and vertical disk structure and show that the disk is warped. We explore the disk at the location of the Beta Pic b super-jupiter and find that the disk surface brightness slope is continuous between 0.5 and 2.0 arcsec, arguing for no change at the separations where Beta Pic b orbits. The two epoch images constrain the disk surface brightness evolution on orbital and radiation pressure blow-out timescales. We place an upper limit of 3% on the disk surface brightness change between 3-5 arcsec, including the locations of the disk warp, and the CO and dust clumps. We discuss the new observations in the context of high-resolution multi-wavelength images and divide the disk asymmetries in two groups: axisymmetric and non-axisymmetric. The axisymmetric structures (warp, large-scale butterfly, etc.) are consistent with disk structure models that include interactions of a planetesimal belt and a non-coplanar giant planet. The non-axisymmetric features, however, require a different explanation.
Aluminum is the heaviest light element displaying large star--to--star variations in Galactic globular clusters (GCs). This element may provide additional insight into the origin of the multiple populations now known to be a common place in GCs, and also the nature of the first-generation stars responsible for a cluster's chemical inhomogeneities. In a previous analysis we found that, unlike more metal-poor GCs, 47 Tuc did not exhibit a strong Na-Al correlation, which motivates a careful study of the similar metallicity but less massive GC M71. We present chemical abundances of O, Na, Al, and Fe for 33 giants in M71 using spectra obtained with the WIYN-Hydra spectrograph. Our spectroscopic analysis finds that, similar to 47 Tuc and in contrast with more metal-poor GCs, M71 stars do not exhibit a strong Na-Al correlation and span a relatively narrow range in [Al/Fe], which are characteristics that GC formation models must reproduce.
We present a new method for engineering the artificial capture of asteroids. Based on theories of the chaos-assisted capture of natural satellites of the giant planets, we show how an unbound asteroid that passes close to a regular region of phase space can be easily moved onto the nearby KAM tori and essentially permanently captured with the Earth's Hill sphere without closing the zero velocity curves. The method has the advantages of a relatively low delta-v requirement and no need for control strategies. An illustration of the method is given for an example asteroid trajectory, demonstrating that it is a viable strategy for the final capture stage of asteroids in the Earth's neighbourhood.
We simulate the formation of a low metallicity (0.01 Zsun) stellar cluster in a dwarf galaxy at redshift z~14. Beginning with cosmological initial conditions, the simulation utilizes adaptive mesh refinement and sink particles to follow the collapse and evolution of gas past the opacity limit for fragmentation, thus resolving the formation of individual protostellar cores. A time- and location-dependent protostellar radiation field, which heats the gas by absorption on dust, is computed by integration of protostellar evolutionary tracks with the MESA code. The simulation also includes a robust non-equilibrium chemical network that self-consistently treats gas thermodynamics and dust-gas coupling. The system is evolved for 18 kyr after the first protostellar source has formed. In this time span, 30 sink particles representing protostellar cores form with a total mass of 81 Msun. Their masses range from ~0.1 Msun to 14.4 Msun with a median mass ~0.5-1 Msun. Massive protostars grow by competitive accretion while lower-mass protostars are stunted in growth by close encounters and many-body ejections. In the regime explored here, the characteristic mass scale is determined by the temperature floor set by the cosmic microwave background and by the onset of efficient dust-gas coupling. It seems unlikely that host galaxies of the first bursts of metal-enriched star formation will be detectable with the James Webb Space Telescope or other next-generation infrared observatories. Instead, the most promising access route to the dawn of cosmic star formation may lie in the scrutiny of metal-poor, ancient stellar populations in the Galactic neighborhood. The observable targets that correspond to the system simulated here are ultra-faint dwarf satellite galaxies such as Bootes II, Segue I and II, and Willman I.
Lieu et al. (2015) have recently claimed that it is possible to substantially improve the sensitivity of radio astronomical observations. In essence, their proposal is to make use of the intensity of the photon shot noise as a measure of the photon arrival rate. Lieu et al. (2015) provide a detailed quantum-mechanical calculation of a proposed measurement scheme that uses two detectors and conclude that this scheme avoids the sensitivity degradation that is associated with photon bunching. If correct, this result could have a profound impact on radio astronomy. Here I present a detailed analysis of the sensitivity attainable using shot-noise measurement schemes that use either one or two detectors, and demonstrate that neither scheme can avoid the photon bunching penalty. I perform both semiclassical and fully quantum calculations of the sensitivity, obtaining consistent results, and provide a formal proof of the equivalence of these two approaches. These direct calculations are furthermore shown to be consistent with an indirect argument based on a correlation method that establishes an independent limit to the sensitivity of shot-noise measurement schemes. Collectively, these results conclusively demonstrate that the photon bunching sensitivity penalty applies to shot noise measurement schemes just as it does to ordinary photon counting, in contradiction to the fundamental claim made by Lieu et al. (2015). The source of this contradiction is traced to a logical fallacy in their argument.
We present moderately-high resolution echelle observations of the nucleus of NGC 5195, the line of sight to which samples intervening interstellar material associated with the outer spiral arm of M51. Our spectra reveal absorption from interstellar Na I, K I, Ca II, and CH+, and from a number of diffuse interstellar bands (DIBs), at a velocity close to that exhibited by H I 21 cm emission from M51 at the position of NGC 5195. The H I column density implied by the equivalent width of the 5780.5 DIB, based on the relationship between W(5780.5) and N(H I) derived for sight lines in the local Galactic interstellar medium, is consistent with the column density obtained from the integrated H I emission. The H2 column density predicted from the observed column density of K I, using the Galactic relationship between N(K I) and N(H2), is comparable to N(H I), suggesting a high molecular fraction (~0.65) for the M51 gas toward NGC 5195. The DIBs toward NGC 5195 are, on average, ~40% weaker than would be expected based on the K I column density, a further indication that the gas in this direction has a high molecular content. The M51 material is characterized also by a high N(Na I)/N(Ca II) ratio (>11), indicative of a high degree of Ca depletion, and a high W(5797.1)/W(5780.5) ratio (~1.6), suggestive of either a very weak ambient radiation field or a significantly shielded environment. A high N(CH+)/N(CH) ratio (>2.3) for the M51 material toward NGC 5195 may be the result of enhanced turbulence due to interactions between M51 and its companion.
We obtained six observations of PSR J1741-2054 using the $Chandra$ ACIS-S detector totaling $\sim$300 ks. By registering this new epoch of observations to an archival observation taken 3.2 years earlier using X-ray point sources in the field of view, we have measured the pulsar proper motion at $\mu =109 \pm 10$ mas/yr. The spectrum of the pulsar can be described by an absorbed power law with photon index $\Gamma$=2.68$\pm$0.04, plus a blackbody with an emission radius of (4.5$^{+3.2}_{-2.5})d_{0.38}$ km, for a DM-estimated distance of $0.38d_{0.38}$ kpc and a temperature of $61.7\pm3.0$ eV. Emission from the compact nebula is well described by an absorbed power law model with a photon index of $\Gamma$ = 1.67$\pm$0.06, while the diffuse emission seen as a trail extending northeast of the pulsar shows no evidence of synchrotron cooling. We also looked for extended features that might represent a jet or torus-like structure using image deconvolution and PSF-subtraction but we find no conclusive evidence of such features immediately surrounding the pulsar.
We present N-body simulations of resonant planets with inclined orbits that show chaotically evolving eccentricities and inclinations that can persist for at least 10 Gyr. A wide range of behavior is possible, from fast, low amplitude variations to systems in which eccentricities reach 0.9999 and inclinations 179.9 degrees. While the orbital elements evolve chaotically, at least one resonant argument always librates. We show that the HD 73526, HD 45364 and HD 60532 systems may be in chaotically-evolving resonances. Chaotic evolution is apparent in the 2:1, 3:1 and 3:2 resonances, and for planetary masses from lunar- to Jupiter-mass. In some cases, orbital disruption occurs after several Gyr, implying the mechanism is not rigorously stable, just long-lived relative to the main sequence lifetimes of solar-type stars. Planet-planet scattering appears to yield planets in inclined resonances that evolve chaotically in about 0.5% of cases. These results suggest that 1) approximate methods for identifying unstable orbital architectures may have limited applicability, 2) the observed close-in exoplanets may be produced during the high eccentricity phases induced by inclined resonances, 3) those exoplanets' orbital planes may be misaligned with the host star's spin axis, 4) systems with resonances may be systematically younger than those without, 5) the distribution of period ratios of adjacent planets detected via transit may be skewed due to inclined resonances, and 6) potentially habitable planets in resonances may have dramatically different climatic evolution than the Earth. The GAIA spacecraft is capable of discovering giant planets in these types of orbits.
It is now widely accepted that dense filaments of molecular gas are integral to the process of stellar birth. While numerical simulations have succeeded in reproducing filamentary structure in turbulent gas and analytic calculations have predicted the formation of dense gas filaments via radial collapse, the exact process(es) that generate/s such filaments which then form prestellar cores within them, is unclear. In this work we therefore study numerically the formation of a dense filament using a relatively simple set-up of a uniform-density cylinder in pressure equilibrium with its confining medium. In particular, we examine if its propensity to form a dense filament and further, to the formation of prestellar cores within this filament bears on the gravitational state of the initial volume of gas. We report a radial collapse leading to the formation of a dense filamentary cloud is likely when the initial volume of gas is at least critically stable (characterised by the approximate equality between the mass line-density for this volume and its maximum value). Though self-gravitating, this volume of gas, however, is not seen to be in free-fall. This post-collapse filament then fragments along its length due to the growth of a Jeans-like instability to form prestellar cores like \emph{beads on a string}. We suggest, dense filaments in typical star-forming clouds classified as gravitationally super-critical under the assumption of : (i) isothermality when in fact, they are not, and (ii) extended radial profiles as against one that is pressure-truncated, thereby causing significant over-estimation of their mass line-density, are unlikely to experience gravitational free-fall. The radial density and temperature profile derived for this post-collapse filament is consistent with that deduced for typical filamentary clouds mapped in recent surveys of nearby star-forming regions.
We present an exposure-time calculator (ETC) for the Immersion Grating Infrared Spectrograph (IGRINS). The signal and noise values are calculated by taking into account the telluric background emission and absorption, the emission and transmission of the telescope and instrument optics, and the dark current and read noise of the infrared detector arrays. For the atmospheric transmission, we apply models based on the amount of precipitable water vapor along the line of sight to the target. The ETC produces the expected signal-to-noise ratio (S/N) for each resolution element, given the exposure-time and number of exposures. In this paper, we compare the simulated continuum S/N for the early-type star HD 124683 and the late-type star GSS 32, and the simulated emission line S/N for the H2 rovibrational transitions from the Iris Nebula NGC 7023 with the observed IGRINS spectra. The simulated S/N from the ETC is overestimate by 10 - 15 % for the sample continuum targets.
During the embedded stage of star formation, bipolar molecular outflows and UV radiation from the protostar are important feedback processes. Our aim is to quantify the feedback, mechanical and radiative, for a large sample of low-mass sources. The outflow activity is compared to radiative feedback in the form of UV heating by the accreting protostar to search for correlations and evolutionary trends. Large-scale maps of 26 young stellar objects, which are part of the Herschel WISH key program are obtained using the CHAMP+ instrument on the APEX (12CO and 13CO 6-5), and the HARP-B instrument on the JCMT (12CO and 13CO 3-2). Maps are used to determine outflow parameters and envelope models are used to quantify the amount of UV-heated gas and its temperature from 13CO 6-5 observations. All sources in our sample show outflow activity and the outflow force, F_CO, is larger for Class 0 sources than for Class I sources, even if their luminosities are comparable. The outflowing gas typically extends to much greater distances than the power-law envelope and therefore influences the surrounding cloud material directly. Comparison of the CO 6-5 results with Herschel-HIFI H2O and PACS high-J CO lines, both tracing currently shocked gas, shows that the two components are linked, even though the transitions do not probe the same gas. The link does not extend down to CO 3-2. The conclusion is that CO 6-5 depends on the shock characteristics (density and velocity), whereas CO 3-2 is more sensitive to conditions in the surrounding environment (density). The radiative feedback is responsible for increasing the gas temperature by a factor of two, up to 30-50 K, on scales of a few thousand AU, particularly along the direction of the outflow. The mass of the UV heated gas exceeds the mass contained in the entrained outflow in the inner ~3000 AU and is therefore at least as important on small scales.
We present long slit spectrophotometric emission line fluxes of bright and extended (<5 arcsec in diameter) Planetary Nebulae (PNe) selected from Acker et al. 1992 catalog with suitable equitorial coordinates for Northern hemisphere. In total, 17 PNe have been choosen and observed in 2008--2010. To measure absolute fluxes, broad slit sizes, ranging from 3.5\arcsec to 7.5\arcsec were used and thus equivalent widths of all observable emission line fluxes were also calculated. Among 17 PNe's observed, line flux measurements of 12 of them were made for the first time. This work also aims to extend the sky coverage of emission line flux standards in Northern hemisphere (Dopita & Hua 1997 - 52 PNe in Southern hemisphere; Wright et al. 2005 - 6 PNe in Northern hemisphere). Electron temperatures and densities, and chemical abundances of these PNe were also calculated in this work. These data is expected to lead the photometric or spectrometric further work for absolute emission line flux measurements needed for \hii regions, supernova remnants etc.
Investigating the spin parameter distribution of subhaloes in two high resolution isolated halo simulations, re- cent work by Onions et al. suggested that typical subhalo spins are consistently lower than the spin distribution found for field haloes. To further examine this puzzle, we have analyzed simulations of a cosmological volume with sufficient resolution to resolve a significant subhalo population. We confirm the result of Onions et al. and show that the typical spin of a subhalo decreases with decreasing mass and increasing proximity to the host halo center. We interpret this as the growing influence of tidal stripping in removing the outer layers, and hence the higher angular momentum particles, of the subhaloes as they move within the host potential. Investigating the redshift dependence of this effect, we find that the typical subhalo spin is smaller with decreasing redshift. This indicates a temporal evolution as expected in the tidal stripping scenario.
Solar prominences are clouds of cool plasma levitating above the solar surface and insulated from the million-degree corona by magnetic fields. They form in regions of complex magnetic topology, characterized by non-potential fields, which can evolve abruptly, disintegrating the prominence and ejecting magnetized material into the heliosphere. However, their physics is not yet fully understood because mapping such complex magnetic configurations and their evolution is extremely challenging, and must often be guessed by proxy from photometric observations.Using state-of-the-art spectro-polarimetric data, we reconstruct the structure of the magnetic field in a prominence. We find that prominence feet harbor helical magnetic fields connecting the prominence to the solar surface below.
Kepler-78b is one of a growing sample of planets similar, in composition and size, to the Earth. It was first detected with NASA's \emph{Kepler} spacecraft and then characterised in more detail using radial velocity follow-up observations. Not only is its size very similar to that of the Earth ($1.2 R_\oplus$), it also has a very similar density ($5.6$ g cm$^{-2}$). What makes this planet particularly interesting is that it orbits its host star every $8.5$ hours, giving it an orbital distance of only $0.0089$ au. What we investigate here is whether or not such a planet could have been perturbed into this orbit by an outer companion on an inclined orbit. In this scenario, the outer perturber causes the inner orbit to undergo Kozai-Lidov cycles which, if the periapse comes sufficiently close to the host star, can then lead to the planet being tidally circularised into a close orbit. We find that this process can indeed produce such very-close-in planets within the age of the host star ($\sim 600 - 900$ Myr), but it is more likely to find such ultra-short-period planets around slightly older stars ($> 1$ Gyr). However, given the size of the Kepler sample and the likely binarity, our results suggest that Kepler-78b may indeed have been perturbed into its current orbit by an outer stellar companion. The likelihood of this happening, however, is low enough that other processes - such as planet-planet scattering - could also be responsible.
We present the highest-quality polarisation profiles to date of 16
non-recycled pulsars and four millisecond pulsars, observed below 200 MHz with
the LOFAR high-band antennas. Based on the observed profiles, we perform an
initial investigation of expected observational effects resulting from the
propagation of polarised emission in the pulsar magnetosphere and the
interstellar medium.
The predictions of magnetospheric birefringence in pulsars have been tested
using spectra of the pulse width and fractional polarisation from
multifrequency data. The derived spectra offer only partial support for the
expected effects of birefringence on the polarisation properties, with only
about half of our sample being consistent with the model's predictions. It is
noted that for some pulsars these measurements are contaminated by the effects
of interstellar scattering. For a number of pulsars in our sample, we have
observed significant variations in the amount of Faraday rotation as a function
of pulse phase, which is possibly an artefact of scattering. These variations
are typically two orders of magnitude smaller than that observed at 1400 MHz by
Noutsos et al. (2009), for a different sample of southern pulsars. In this
paper we present a possible explanation for the difference in magnitude of this
effect between the two frequencies, based on scattering. Finally, we have
estimated the magnetospheric emission heights of low-frequency radiation from
four pulsars, based on the phase lags between the flux-density and the PA
profiles, and the theoretical framework of Blaskiewicz, Cordes & Wasserman
(1991). These estimates yielded heights of a few hundred km; at least for PSR
B1133+16, this is consistent with emission heights derived based on
radius-to-frequency mapping, but is up to a few times larger than the recent
upper limit based on pulsar timing.
Strong evidence exists for a highly significant correlation between the radio
flux density and gamma-ray energy flux in blazars revealed by Fermi. However,
there are central issues that need to be clarified in this field: what are the
counterparts of the about 30% of gamma-ray sources that are as yet
unidentified? Are they just blazars in disguise or they are something more
exotic, possibly associated with dark matter? How would they fit in the
radio-gamma ray connection studied so far?
With their superb sensitivity, SKA1-MID and SKA1-SUR will help to resolve all
of these questions. Even more, while the radio-MeV/GeV connection has been
firmly established, a radio-VHE connection has been entirely elusive so far.
The advent of CTA in the next few years and the expected CTA-SKA1 synergy will
offer the chance to explore this connection, even more intriguing as it
involves the opposite ends of the electromagnetic spectrum and the acceleration
of particles up to the highest energies.
We are already preparing to address these questions by exploiting data from
the various SKA pathfinders and precursors. We have obtained 18 cm European
VLBI Network observations of E>10 GeV sources, with a detection rate of 83%.
Moreover, we are cross correlating the Fermi catalogs with the MWA
commissioning survey: when faint gamma-ray sources are considered, pure
positional coincidence is not significant enough for selecting counterparts and
we need an additional physical criterion to pinpoint the right object. It can
be radio spectral index, variability, polarization, or compactness, needing
high angular resolution in SKA1-MID; timing studies can also reveal pulsars,
which are often found from dedicated searches of unidentified gamma-ray
sources. SKA will be the ideal instrument for investigating these
characteristics in conjunction with CTA.
(abridged)
The South Galactic Cap $u$-band Sky Survey (SCUSS) provides a deep $u$-band imaging of about 5000 deg$^2$ in south Galactic cap. It is about 1.5 mag deeper than the SDSS $u$-band. In this paper we evaluate the capability of quasar selection using both SCUSS and SDSS data, based on considerations of the deep SCUSS $u$-band imaging and two-epoch $u$-band variability. We find that the combination of the SCUSS $u$-band and the SDSS $griz$ band allows us to select more faint quasars and more quasars at redshift around 2.2 than the selection only with the SDSS $ugriz$ data. Quasars have significant $u$-band variabilities. The fraction of quasars with large two-epoch variability is much higher than that of stars. The selection by variability can select both low-redshift quasars with ultraviolet excess and mid-redshift ($2 < z <3.5$) quasars where quasar selection by optical colors is inefficient. The above two selections are complementary and make full use of the SCUSS u-band advantages.
Through CO mm-line and optical spectroscopy, we investigate the properties of the Fried Egg nebula IRAS 17163-3907, which has recently been proposed to be one of the rare members of the yellow hypergiant class. The CO J=2-1 and J=3-2 emission arises from a region within 20" of the star and is clearly associated with the circumstellar material. The CO lines show a multi-component asymmetrical profile, and an unexpected velocity gradient is resolved in the east-west direction, suggesting a bipolar outflow. This is in contrast with the apparent symmetry of the dust envelope as observed in the infrared. The optical spectrum of IRAS 17163-3907 between 5100 and 9000 {\AA} was compared with that of the archetypal yellow hypergiant IRC+10420 and was found to be very similar. These results build on previous evidence that IRAS 17163-3907 is a yellow hypergiant.
The heat capacity of neutron matter is studied over the range of densities and temperatures prevailing in neutron-star crusts, allowing for the transition to a superfluid phase at temperatures below some critical temperature $T_{sf}$ and including the transition to the classical limit. Finite temperature Hartree-Fock-Bogoliubov equations (FTHFB) are solved and compared to existing approximate expressions. In particular, the formula given by Levenfish and Yakovlev is found to reproduce the numerical results with a high degree of accuracy for temperatures $T\leq T_{sf}$. In the non-superfluid phase, $T\geq T_{sf}$, the linear approximation is valid only at temperature $T\ll T_{{\rm F} n}$ ($T_{{\rm F} n}$ being the Fermi temperature of the neutron gas) which is rarely the case in the shallow layers of the neutron star's crust. A non-perturbative interpolation between the quantal and the classical regimes is proposed here. The heat capacity, conveniently parametrized solely in terms of $T_{sf}$, $T_{{\rm F} n}$, and the neutron number density $n_n$, can be easily implemented in neutron-star cooling simulations.
The formation of Wolf-Rayet central stars of planetary nebulae ([WR] CSPNe) whose spectroscopic appearance mimics massive WR stars remains poorly understood. Least understood is the nature and frequency of binary companions to [WR] CSPNe that may explain their H-deficiency. We have conducted a systematic radial velocity (RV) study of 6 [WR] CSPNe to search for post-common-envelope (post-CE) [WR] binaries. We used a cross-correlation method to construct the RV time-series as successfully done for massive close binary WR stars. No significant RV variability was detected for the late-[WC] type nuclei of Hen 2-113, Hen 3-1333, PMR~2 and Hen 2-99. Significant, large-amplitude variability was found in the [WC4] nucleus of NGC 5315. In the [WO1] nucleus of NGC 5189 we discovered significant periodic variability that reveals a close binary with $P_\mathrm{orb}=4.04\pm0.1$ d. We measured a semi-amplitude of $62.3\pm1.3$ km s$^{-1}$ that gives a companion mass of $m_2\ge0.5$ $M_\odot$ or $m_2=0.84$ $M_\odot$ (assuming $i=45^\circ$). The most plausible companion type is a massive WD as found in Fleming 1. The spectacular nebular morphology of NGC 5189 fits the pattern of recently discovered post-CE PNe extremely well with its dominant low-ionisation structures (e.g. as in NGC 6326) and collimated outflows (e.g. as in Fleming 1). The anomalously long 4.04 d orbital period is either a once-off (e.g. NGC 2346) or it may indicate there is a sizeable population of [WR] binaries with massive WD companions in relatively wide orbits, perhaps influenced by interactions with the strong [WR] wind.
We analyze publicly available void catalogs of the Baryon Oscillation Spectroscopic Survey Data Release 10 at redshifts $0.4<z<0.7$. The first goal of this paper is to extend the Cosmic Microwave Background stacking analysis of previous spectroscopic void samples at $z<0.4$. In addition, the DR10 void catalog provides the first chance to spectroscopically probe the volume of the Granett et al. (2008) supervoid catalog that constitutes the only set of voids which has shown a significant detection of a cross-correlation signal between void locations and average CMB chill. We found that the positions of voids identified in the spectroscopic DR10 CMASS galaxy catalog typically do not coincide with the locations of the Granett et al. (2008) supervoids in the overlapping volume, in spite of the presence of large underdense regions of high void-density in DR10. The stacking of filtered CMB temperatures at these different void locations shows a $\Delta T = - 6.2 \pm 3.4 ~\mu K$ signal for the 120 largest voids, otherwise the correlation is washed out by statistical uncertainties. This correlation is, however, significantly lower than the $\Delta T = - 11.5 \pm 3.7 ~\mu K$ we found by stacking 35 of the 50 Granett et al. (2008) supervoids available in the DR10 volume. This failure to reproduce the signal with a different void catalog may be due to systematic differences in the detection of voids in photometric and spectroscopic samples.
WINERED is a newly built high-efficiency (throughput$ > 25-30\%$) and high-resolution spectrograph customized for short NIR bands at 0.9-1.35 ${\rm \mu}$m. WINERED is equipped with ambient temperature optics and a cryogenic camera using a 1.7 ${\rm \mu}$m cut-off HgCdTe HAWAII-2RG array detector. WINERED has two grating modes: one with a conventional reflective echelle grating (R$\sim$28,300), which covers 0.9-1.35 $\mu$m simultaneously, the other with ZnSe or ZnS immersion grating (R$\sim$100,000). We have completed the development of WINERED except for the immersion grating, and started engineering and science observations at the Nasmyth platform of the 1.3 m Araki Telescope at Koyama Astronomical Observatory of Kyoto-Sangyo University in Japan. We confirmed that the spectral resolution ($R\sim$ 28,300) and the throughput ($>$ 40\% w/o telescope/atmosphere/array QE) meet our specifications. We measured ambient thermal backgrounds (e.g., 0.06 ${\rm [e^{-}/sec/pixel]}$ at 287 K), which are roughly consistent with that we expected. WINERED is a portable instrument that can be installed at any telescope with Nasmyth focus as a PI-type instrument. If WINERED is installed on a 10 meter telescope, the limiting magnitude is expected to be J=18-19, which can provide high-resolution spectra with high quality even for faint distant objects.
The black hole system 1E${\thinspace}$1740.7$-$2942 is usually the brightest hard X-ray source (above 20 keV) near the Galactic Center, but presents some epochs of low emission (below the INTEGRAL detection limit, for example). In this work, we present the results of studies on 1E${\thinspace}$1740.7$-$2942 over 10 years, using the instruments ISGRI/IBIS and JEM-X, both on board the INTEGRAL observatory. We fit the spectra with both the compTT and cutoffpl models. According to the fits and taking the mean value over the 10 years, we have obtained a plasma temperature in the range $\sim$20$\,-\,$90$\thinspace$ keV, and an average powerlaw index of 1.41 ($\sigma$=0.25). We have also made a Lomb$\,-\,$Scargle periodogram of the flux in the 50$\,-\,$20${\thinspace}$keV band and found two tentative periods at 2.90 and 3.99${\thinspace}$days. We present here the preliminary results of this ongoing work.
IceTop, the surface component of the IceCube detector, has been used to measure the energy spectrum of cosmic ray primaries in the range between 1.58 PeV and 1.26 EeV. It can also be used to study the low energy muons in air showers by looking at large distances (> 300m) from the shower axis. We will show the muon lateral distribution function at large lateral distances as measured with IceTop and discuss the implications of this measurement. We also discuss the prospects for low energy muon studies with IceTop.
We describe the analysis of the seven broad-band X-ray continuum observations of the archetypal Seyfert 1 galaxy NGC 5548 that were obtained with XMM-Newton or Chandra, simultaneously with high-energy (> 10 keV) observations with NuSTAR and INTEGRAL. These data were obtained as part of a multiwavelength campaign undertaken from the summer of 2013 till early 2014. We find evidence of a high-energy cut-off in at least one observation, which we attribute to thermal Comptonization, and a constant reflected component that is likely due to neutral material at least a few light months away from the continuum source. We confirm the presence of strong, partial covering X-ray absorption as the explanation for the sharp decrease in flux through the soft X-ray band. The obscurers appear to be variable in column density and covering fraction on time scales as short as weeks. A fit of the average spectrum over the range 0.3-400 keV with a realistic Comptonization model indicates the presence of a hot corona with a temperature of 40(+40,-10) keV and an optical depth of 2.7(+0.7,-1.2) if a spherical geometry is assumed.
We study a sample of 11 Type II supernovae (SNe) discovered by the OGLE-IV survey. All objects have well sampled I-band light curves, and at least one spectrum. We find that 3 or 4 of the 11 SNe have a declining light curve, making them SNe II-L, while the rest have plateaus that can be as short as 70d, unlike the 100d typically found in nearby galaxies. These SNe are also brighter than found in the local Universe, and show that magnitude limited surveys find SNe that are different than found in nearby galaxies. We discuss this sample in the context of understanding Type II SNe as a class and their suggested use as standard candles.
Theories of galaxy formation predict the existence of extended gas halo around spiral galaxies. If there are 10-100 nG magnetic fields at several ten kpc distances from the galaxies, extended galactic cosmic ray (CR) haloes could also exist. Galactic CRs could interact with the tenuous hot halo gas to produce observable $\gamma$-rays. In this paper we have performed search for such a halo around the M31 galaxy -- the closest large spiral galaxy. Our analysis of 5.5 years of the Fermi LAT data revealed the presence of a spatially extended emission excess around M31. The data can be fitted using the simplest morphology of a uniformly bright circle. The best fit gave a 4.4$\sigma$ significance for a $3^{\circ}$ (40 kpc) halo with photon flux of $\sim (1.9\pm1.1)\times 10^{-9} ~\mathrm{cm^{-2}s^{-1}}$ and luminosity $(8.4\pm4.6)\times 10^{38} ~\mathrm{erg~s^{-1}}$ in the energy range 0.3--100 GeV. The presence of such a halo compellingly shows that a 10-100 nG magnetic field should extend around M31 up to a 40 kpc distance.
We present results demonstrating the time resolution and $\mu$/e separation capabilities with a new concept of an EAS detector capable for measurements of cosmic rays arriving with large zenith angles. This kind of detector has been designed to be a part of a large area (several square kilometers) surface array designed to measure Ultra High Energy (10-200 PeV) $\tau$ neutrinos using the Earth-skimming technique. A criteria to identify electron-gammas is also shown and the particle identification capability is tested by measurements in coincidence with the KASKADE-GRANDE experiment in Karlsruhe, Germany.
We address how to construct an infinitely cyclic universe model. A major consideration is to make the entropy cyclic which requires the entropy to be reset to zero in each cycle expansion to turnaround, to contraction, to bounce, etc. Here we reset entropy at the turnaround by selecting the visible universe from the multiverse which is generated by the accelerated expansion. In the model, the observed homogeneity is explained by the low entropy at the bounce, The observed flatness arises from the contraction together with the reduction in size between the expanding and contracting universe. The present flatness is predicted to be very precise.
We investigate a simplified model of dark matter where a Majorana fermion $\chi$ coannihilates with a colored scalar top partner $\tilde{t}$. We explore the cosmological history, with particular emphasis on the most relevant low-energy parameters: the mass splitting between the dark matter and the coannihilator, and the Yukawa coupling $y_\chi$ that connects these fields to the Standard Model top quarks. We also allow a free quartic coupling $\lambda_h$ between a pair of Higgs bosons and $\tilde{t}$ pairs. We pay special attention to the case where the values take on those expected where $\tilde{t}$ corresponds to the superpartner of the right-handed top, and $\chi$ is a bino. Direct detection, indirect detection, and colliders are complementary probes of this simple model.
We investigate the direct detection phenomenology of a class of dark matter (DM) models in which DM does not directly interact with nuclei, {but rather} the products of its annihilation do. When these annihilation products are very light compared to the DM mass, the scattering in direct detection experiments is controlled by relativistic kinematics. This results in a distinctive recoil spectrum, a non-standard and or even {\it absent} annual modulation, and the ability to probe DM masses as low as a $\sim$10 MeV. We use current LUX data to show that experimental sensitivity to thermal relic annihilation cross sections has already been reached in a class of models. Moreover, the compatibility of dark matter direct detection experiments can be compared directly in $E_{min}$ space without making assumptions about DM astrophysics. Lastly, when DM has direct couplings to nuclei, the limit from annihilation to relativistic particles in the Sun can be stronger than that of conventional non-relativistic direct detection by more than three orders of magnitude for masses in a 2-8 GeV window.
Motivated by the mathematic theory of split-complex numbers (or hyperbolic numbers, also perplex numbers) and the split-quaternion numbers (or coquaternion numbers), we define the notion of split-complex scalar field and the split-quaternion scalar field. Then we explore the cosmic evolution of these scalar fields in the background of spatially flat Friedmann-Robertson-Walker Universe. We find that both the quintessence field and the phantom field could naturally emerge in these scalar fields. Introducing the metric of field space, these theories fall into a subclass of the multi-field theories which have been extensively studied in inflationary cosmology. Using the brane world model, the split-complex Dirac-Born-Infeld Lagrangian is constructed and analyzed.
Among the prominent low-mass dark matter candidates is the QCD axion but also other light and weakly interacting particles beyond the Standard Model. We review briefly the case for such dark matter and give an overview on most recent experimental efforts within laboratory searches, where we focus on experiments exploiting a potential electromagnetic coupling of such particles.
Cadmium Zinc Telluride and Cadmium Telluride are the detector materials of choice for the detection of X-rays in the X-ray energy band E >= 5keV with excellent spatial and spectral resolution and without cryogenic cooling. Owing to recent breakthroughs in grazing incidence mirror technology, next-generation hard X-ray telescopes will achieve angular resolution between 5 and 10 arc seconds - about an order of magnitude better than that of the NuSTAR hard X-ray telescope. As a consequence, the next generation of X-ray telescopes will require pixelated X-ray detectors with pixels on a grid with a lattice constant of <= 250um. Additional detector requirements include a low energy threshold of less than 5keV and an energy resolution of less than one keV. The science drivers for a high angular-resolution X-ray mission include studies and measurements of black hole spins, the cosmic evolution of super-massive black holes, active galactic nuclei feedback, and the behaviour of matter at very high densities. In this contribution, we report on our R&D studies with the goal to optimise small-pixel Cadmium Zinc Telluride and Cadmium Telluride detectors.
Links to: arXiv, form interface, find, astro-ph, recent, 1501, contact, help (Access key information)
The XENON100 experiment is the second phase of the XENON program for the direct detection of the dark matter in the universe. The XENON100 detector is a two-phase Time Projection Chamber filled with 161 kg of ultra pure liquid xenon. The results from 224.6 live days of dark matter search with XENON100 are presented. No evidence for dark matter in the form of WIMPs is found, excluding spin-independent WIMP-nucleon scattering cross sections above 2 $\times$ 10$^{-45}$ cm$^2$ for a 55 GeV/c$^2$ WIMP at 90% confidence level (C.L.). The most stringent limit is established on the spin-dependent WIMP-neutron interaction for WIMP masses above 6 GeV/c$^2$, with a minimum cross section of 3.5 $\times$ 10$^{-40}$ cm$^2$ (90% C.L.) for a 45 GeV/c$^2$ WIMP. The same dataset is used to search for axions and axion-like-particles. The best limits to date are set on the axion-electron coupling constant for solar axions, $g_{Ae}$ < 7.7 $\times$ 10$^{-12}$ (90% C.L.), and for axion-like-particles, $g_{Ae}$ < 1 $\times$ 10$^{-12}$ (90% C.L.) for masses between 5 and 10 keV/c$^2$.
Extracting Times of Arrival from pulsar radio signals depends on the knowledge of the pulsars pulse profile and how this template is generated. We examine pulsar template generation with Bayesian methods. We will contrast the classical generation mechanism of averaging intensity profiles with a new approach based on Bayesian inference. We introduce the Bayesian measurement model imposed and derive the algorithm to reconstruct a "statistical template" out of noisy data. The properties of these "statistical templates" are analysed with simulated and real measurement data from PSR B1133+16. We explain how to put this new form of template to use in analysing secondary parameters of interest and give various examples: We implement a nonlinear filter for determining ToAs of pulsars. Applying this method to data from PSR J1713+0747 we derive ToAs self consistently, meaning all epochs were timed and we used the same epochs for template generation. While the average template contains fluctuations and noise as unavoidable artifacts, we find that the "statistical template" derived by Bayesian inference quantifies fluctuations and remaining uncertainty. This is why the algorithm suggested turns out to reconstruct templates of statistical significance from ten to fifty single pulses. A moving data window of fifty pulses, taking out one single pulse at the beginning and adding one at the end of the window unravels the characteristics of the methods to be compared. It shows that the change induced in the classical reconstruction is dominated by random fluctuations for the average template, while statistically significant changes drive the dynamics of the proposed method's reconstruction. The analysis of phase shifts with simulated data reveals that the proposed nonlinear algorithm is able to reconstruct correct phase information along with an acceptable estimation of the remaining uncertainty.
We investigate the position of the radio core in a blazar by multi-epoch astrometric observations at 43 GHz. Using the VLBI Exploration of Radio Astrometry (VERA), we have conducted four adjacent observations in February 2011 and another four in October 2011, and succeeded in measuring the position of the radio core in the TeV blazar Mrk 501 relative to a distant compact quasar NRAO 512. During our observations, we find that (1) there is no positional change within ~0.2 mas or ~2.0 pc de-projected with $\pm1\sigma$ error for the weighted-mean phase-referenced positions of Mrk 501 core relative to NRAO 512 over four adjacent days, and (2) there is an indication of position change for 3C 345 core relative to NRAO 512. By applying our results to the standard internal shock model for blazars, we constrain the bulk Lorenz factors of the ejecta.
We study the impact of relatively weak AGN feedback on the interstellar medium of intermediate and massive elliptical galaxies. We find that the AGN activity, while globally heating the ISM, naturally stimulates some degree of hot gas cooling on scales of several kpc. This process generates the persistent presence of a cold ISM phase, with mass ranging between 10$^4$ and $\gtrsim$ 5 $\times$ 10$^7$ M$_\odot$, where the latter value is appropriate for group centered, massive galaxies. Widespread cooling occurs where the ratio of cooling to free-fall time before the activation of the AGN feedback satisfies $t_{cool}/t_{ff} \lesssim 70$, that is we find a less restrictive threshold than commonly quoted in the literature. This process helps explaining the body of observations of cold gas (both ionized and neutral/molecular) in Ellipticals and, perhaps, the residual star formation detected in many early-type galaxies. The amount and distribution of the off-center cold gas vary irregularly with time. The cold ISM velocity field is irregular, initially sharing the (outflowing) turbulent hot gas motion. Typical velocity dispersions of the cold gas lie in the range 100-200 km/s. Freshly generated cold gas often forms a cold outflow and can appear kinematically misaligned with respect to the stars. We also follow the dust evolution in the hot and cold gas. We find that the internally generated cold ISM has a very low dust content, with representative values of the dust-to-gas ratio of 10$^{-4}$- 10$^{-5}$. Therefore, this cold gas can escape detection in the traditional dust-absorption maps.
Several problems regarding the process of galaxy formation are still open. One of them is the role played by AGN phenomena in contributing to galaxy build-up and to Star Formation (SF) quenching. On the other hand, the theory of AGN formation predicts these phenomena to be correlated with the host-galaxy environment, thus opening for links between SF quenching, environment and AGN phenomena in the galaxy formation paradigm. This work is focused on the study of the correlation between environmental density and radio AGN presence. Using data from the photometric COSMOS survey and the radio 1.4GHz VLA-COSMOS, a sample of radio AGNs has been defined. The environment has been studied throughout the use of the richness distributions inside a parallelepipedon with base side of 1 Mpc and height proportional to the photometric redshift precision. Richness distributions have been compared as a function of both the redshift and the relative evolution of the stellar masses of galaxies and AGN hosts up to $z=2$. Radio AGNs are always located in environments that are significantly denser than those around galaxies in which radio emission is absent. Therefore, the environment seems to enhance the probability of a galaxy to host a radio AGN. Moreover, a distinction between high power and low power radio AGNs leads to the conclusion that the significance in the environmental effect is maintained only for low power radio sources. By studying the evolution of stellar masses it is possible to conclude that radio AGN presence is a phenomenon related to quiescent galaxies up to $z=2$, with a significant increase of the fraction of quiescent galaxies hosting a radio AGN with decreasing redshift. Hints of an environmental effect are present. The results found with this work lead to conclude that denser environments play a significant role in enhancing the probability of a galaxy to host a low power radio AGN.
We study the broad emission line blazars detected in the gamma-ray band by the Large Area Telescope onboard the Fermi satellite and with the optical spectrum studied by Shaw et al. (2012, 2013). The observed broad line strength provides a measure of the ionizing luminosity of the accretion disk, while the gamma-luminosity is a proxy for the bolometric non-thermal beamed jet emission. The resulting sample, composed by 217 blazars, is the best suited to study the connection between accretion and jet properties. We compare the broad emission line properties of these blazars with those of radio-quiet and radio-loud quasars present in the Sloan Digital Sky Survey, to asses differences and similarities of the disc luminosity and the virial black hole mass. For most sources, we could derive the black hole mass by reproducing the IR-optical-UV data with a standard accretion disc spectrum, and we compared the black hole masses derived with the two methods. The distributions of the masses estimated in the two ways agree satisfactorily. We then apply a simple, one-zone, leptonic model to all the 217 objects of our sample. The knowledge of the black hole mass and disc luminosity helps to constrain the jet parameters. On average they are similar to what found by previous studies of smaller samples of sources.
Ultra-deep observations of ECDF-S with Chandra and XMM-Newton enable a search for extended X-ray emission down to an unprecedented flux of $2\times10^{-16}$ ergs s$^{-1}$ cm$^{-2}$. We present the search for the extended emission on spatial scales of 32$^{\prime\prime}$ in both Chandra and XMM data, covering 0.3 square degrees and model the extended emission on scales of arcminutes. We present a catalog of 46 spectroscopically identified groups, reaching a redshift of 1.6. We show that the statistical properties of ECDF-S, such as logN-logS and X-ray luminosity function are broadly consistent with LCDM, with the exception that dn/dz/d$\Omega$ test reveals that a redshift range of $0.2<z<0.5$ in ECDF-S is sparsely populated. The lack of nearby structure, however, makes studies of high-redshift groups particularly easier both in X-rays and lensing, due to a lower level of clustered foreground. We present one and two point statistics of the galaxy groups as well as weak-lensing analysis to show that the detected low-luminosity systems are indeed low-mass systems. We verify the applicability of the scaling relations between the X-ray luminosity and the total mass of the group, derived for the COSMOS survey to lower masses and higher redshifts probed by ECDF-S by means of stacked weak lensing and clustering analysis, constraining any possible departures to be within 30\% in mass. Abridged.
By considering linear-order departures from general relativity, we compute a novel expression for the weak lensing convergence power spectrum under alternative theories of gravity. This comprises an integral over a 'kernel' of general relativistic quantities multiplied by a theory-dependent 'source' term. The clear separation between theory-independent and -dependent terms allows for an explicit understanding of each physical effect introduced by altering the theory of gravity. We take advantage of this to explore the degeneracies between gravitational parameters in weak lensing observations.
Ever since WMAP announced its first results, different analyses have shown that there is weak evidence for several large-scale anomalies in the CMB data. While the evidence for each anomaly appears to be weak, the fact that there are multiple seemingly unrelated anomalies makes it difficult to account for them via a single statistical fluke. So, one is led to considering a combination of these anomalies. But, if we "hand-pick" the anomalies (test statistics) to consider, we are making an \textit{a posteriori} choice. In this article, we propose two statistics that do not suffer from this problem. The statistics are linear and quadratic combinations of the $a_{\ell m}$'s with random co-efficients, and they test the null hypothesis that the $a_{\ell m}$'s are independent, normally-distributed, zero-mean random variables with an $m$-independent variance. The motivation for such statistics is generality; equivalently, it is a non \textit{a posteriori} choice. But, a very useful by-product of considering such statistics is this: Because most physical models that lead to large-scale anomalies result in coupling multiple $\ell$ and $m$ modes, the "coherence" of this coupling should get enhanced if a combination of different modes is considered. Using fiducial data, we demonstrate that the method works and discuss how it can be used with actual CMB data to make quite general statements about how incompatible the data are with the null hypothesis.
As of 2014 August, the Nuclear Spectroscopic Telescope Array (NuSTAR) had observed ~30 X-ray binaries either as part of the planned program, as targets of opportunity, or for instrument calibration. The main science goals for the observations include probing the inner part of the accretion disk and constraining black hole spins via reflection components, providing the first observations of hard X-ray emission from quiescent Low Mass X-ray Binaries (LMXBs), measuring cyclotron lines from accreting pulsars, and studying type I X-ray bursts from neutron stars. Here, we describe the science objectives in more depth and give an overview of the NuSTAR observations that have been carried out to achieve the objectives. These include observation of four "IGR" High Mass X-ray Binaries (HMXBs) discovered by INTEGRAL. We also summarize the results that have been obtained and their implications. Among the IGR HMXBs, we focus on the discovery of a cyclotron line in the spectrum of IGR J17544-2619.
VERITAS is an array of four 12-m atmospheric Cherenkov telescopes, designed to observe the very-high-energy (VHE; E $\geq$100 GeV) sky. Since 2007, it has detected more than 20 extra-galactic sources, the majority of which are active-galactic-nuclei of the blazar class. In this paper we present a selection of the most recent results from the VERITAS blazar observing program, in particular in the context of broad-band multi-wavelength campaigns. Four results are highlighted: the detection of the flat-spectrum-radio-quasar PKS 1222+216 (4C +21.35) during March 2014; the $\gamma$-ray flare from the BL Lac object 1ES 1727+502, observed under bright-moonlight conditions during May 2013; the bright $\gamma$-ray flare from the BL Lac object 1ES 1011+496 during February 2014; and the long-term campaign on the BL Lac object PKS 1424+240, currently the farthest (z > 0.60), persistent VHE emitter
We present an analysis of a new, large sample of field blue-straggler stars (BSSs) in the thick disk and halo system of the Galaxy, based on stellar spectra obtained during the Sloan Digital Sky Survey (SDSS) and the Sloan Extension for Galactic Understanding and Exploration (SEGUE). Using estimates of stellar atmospheric parameters obtained from application of the SEGUE Stellar Parameter Pipeline, we obtain a sample of some 8000 BSSs, which are considered along with a previously selected sample of some 4800 blue horizontal-branch (BHB) stars. We derive the ratio of BSSs to BHB stars, F$_{\rm BSS/BHB}$, as a function of Galactocentric distance and distance from the Galactic plane. The maximum value found for F$_{\rm BSS/BHB}$ is $\sim~$4.0 in the thick disk (at 3 kpc $<$ $|$Z$|$ $<$ 4 kpc), declining to on the order of $\sim~1.5-2.0$ in the inner-halo region; this ratio continues to decline to $\sim~$1.0 in the outer-halo region. We associate a minority of field BSSs with a likely extragalactic origin; at least 5$\%$ of the BSS sample exhibit radial velocities, positions, and distances commensurate with membership in the Sagittarius Stream.
Employing the data from orbital periods and masses of extra-solar planets in 166 multiple planetary systems, the period-ratio and mass-ratio of adjacent planet pairs are studied. The correlation between the period-ratio and mass-ratio is confirmed and found to have a correlation coefficient of 0.5303 with a 99% confidence interval (0.3807, 0.6528). A comparison with the distribution of synthetic samples from a Monte Carlo simulation reveals the imprint of planet-planet interactions on the formation of adjacent planet pairs in multiple planetary systems.
The observed tight radio/X-ray correlation in the low spectral state of some black hole X-ray binaries implies the strong coupling of the accretion and jet. The correlation of $L_{\rm R} \propto L_{\rm X}^{\sim 0.5-0.7}$ was well explained by the coupling of a radiatively inefficient accretion flow and a jet. Recently, however, a growing number of sources show more complicated radio/X-ray correlations, e.g., $L_{\rm R} \propto L_{\rm X}^{\sim 1.4}$ for $L_{\rm X}/L_{\rm Edd} \gtrsim 10^{-3}$, which is suggested to be explained by the coupling of a radiatively efficient accretion flow and a jet. In this work, we interpret the deviation from the initial radio/X-ray correlation for $L_{\rm X}/L_{\rm Edd} \gtrsim 10^{-3}$ with a detailed disc corona-jet model. In this model, the disc and corona are radiatively and dynamically coupled. Assuming a fraction of the matter in the accretion flow, $\eta\equiv \dot M_{\rm jet}/\dot M$, is ejected to form the jet, we can calculate the emergent spectrum of the disc corona-jet system. We calculate $L_{\rm R}$ and $L_{\rm X}$ at different $\dot M$, adjusting $\eta$ to fit the observed radio/X-ray correlation of the black hole X-ray transient H1743-322 for $L_{\rm X}/L_{\rm Edd}> 10^{-3}$. It is found that always the X-ray emission is dominated by the disc corona and the radio emission is dominated by the jet. We noted that the value of $\eta$ for the deviated radio/X-ray correlation for $L_{\rm X}/L_{\rm Edd} > 10^{-3}$, is systematically less than that of the case for $L_{\rm X}/L_{\rm Edd} < 10^{-3}$, which is consistent with the general idea that the jet is often relatively suppressed at the high luminosity phase in black hole X-ray binaries.
We study a nonsingular bounce inflation model, which can drive the early universe from a contracting phase, bounce into an ordinary inflationary phase, followed by the reheating process. Besides the bounce that avoided the Big-Bang singularity which appears in the standard cosmological scenario, we make use of the Horndesky theory and design the kinetic and potential forms of the lagrangian, so that neither of the two big problems in bouncing cosmology, namely the ghost and the anisotropy problems, will appear. The cosmological perturbations can be generated either in the contracting phase or in the inflationary phase, where in the latter the power spectrum will be scale-invariant and fit the observational data, while in the former the perturbations will have nontrivial features that will be tested by the large scale structure experiments. We also fit our model to the CMB TT power spectrum.
We employ a suite of 75 simulations of galaxies in idealised major mergers (stellar mass ratio ~2.5:1), with a wide range of orbital parameters, to investigate the spatial extent of interaction-induced star formation. Although the total star formation in galaxy encounters is generally elevated relative to isolated galaxies, we find that this elevation is a combination of intense enhancements within the central kpc and moderately suppressed activity at large galacto-centric radii. The radial dependence of the star formation enhancement is stronger in the less massive galaxy than in the primary, and is also more pronounced in mergers of more closely aligned disc spin orientations. Conversely, these trends are almost entirely independent of the encounter's impact parameter and orbital eccentricity. Our predictions of the radial dependence of triggered star formation, and specifically the suppression of star formation beyond kph-scales, will be testable with the next generation of integral-field spectroscopic surveys.
We present the results of a search for potential transit signals in the full 17-quarter data set collected during Kepler's primary mission that ended on May 11, 2013, due to the on-board failure of a second reaction wheel needed to maintain high precision, fixed, pointing. The search includes a total of 198,646 targets, of which 112,001 were observed in every quarter and 86,645 were observed in a subset of the 17 quarters. We find a total of 12,669 targets that contain at least one signal that meets our detection criteria: periodicity of the signal, a minimum of three transit events, an acceptable signal-to-noise ratio, and four consistency tests that suppress false positives. Each target containing at least one transit-like pulse sequence is searched repeatedly for other signals that meet the detection criteria, indicating a multiple planet system. This multiple planet search adds an additional 7,698 transit-like signatures for a total of 20,367. Comparison of this set of detected signals with a set of known and vetted transiting planet signatures in the Kepler field of view shows that the recovery rate of the search is 90.3%. We review ensemble properties of the detected signals and present various metrics useful in validating these potential planetary signals. We highlight previously undetected planetary candidates, including several small potential planets in the habitable zone of their host stars.
In the cold dark matter (CDM) paradigm, bulges easily form through galaxy mergers, either major or minor, or through clumpy disks in the early universe, where clumps are driven to the center by dynamical friction. Also pseudo-bulges, with a more disky morphology and kinematics, can form more slowly through secular evolution of a bar, where resonant stars are elevated out of the plane, in a peanut/box shape. As a result, in CDM cosmological simulations, it is very difficult to find a bulgeless galaxy, while they are observed very frequently in the local universe. A different picture emerges in alternative models of the missing mass problem. In MOND (MOdified Newtonian Dynamics), galaxy mergers are much less frequent, since the absence of dark matter halos reduces the dynamical friction between two galaxies. Also, while clumpy galaxies lead to rapid classical bulge formation in CDM, the inefficient dynamical friction with MOND in the early-universe galaxies prevents the clumps to coalesce together in the center to form spheroids. This leads to less frequent and less massive classical bulges. Bars in MOND are more frequent and stronger, and have a more constant pattern speed, which modifies significantly the pseudo-bulge morphology. The fraction of pseudo-bulges is expected to be dominant in MOND.
We investigate the correlation between extinction and H~{\sc i} and CO emission at intermediate and high Galactic latitudes ($|b|>10\degr$) within the footprint of the Xuyi Schmidt Telescope Photometric Survey of the Galactic anticentre (XSTPS-GAC) on small and large scales. In Paper I (Chen et al. 2014), we present a three-dimensional dust extinction map within the footprint of XSTPS-GAC, covering a sky area of over 6,000\,deg$^2$ at a spatial angular resolution of 6\,arcmin. In the current work, the map is combined with data from gas tracers, including H~{\sc i} data from the Galactic Arecibo L-band Feed Array H~{\sc i} survey and CO data from the Planck mission, to constrain the values of dust-to-gas ratio $DGR=A_V/N({\rm H})$ and CO-to-$\rm H_2$ conversion factor $X_{\rm CO}=N({\rm H_2})/W_{\rm CO}$ for the entire GAC footprint excluding the Galactic plane, as well as for selected star-forming regions (such as the Orion, Taurus and Perseus clouds) and a region of diffuse gas in the northern Galactic hemisphere. For the whole GAC footprint, we find $DGR=(4.15\pm0.01) \times 10^{-22}$\,$\rm mag\,cm^{2}$ and $X_{\rm CO}=(1.72 \pm 0.03) \times 10^{20}$\,$\rm cm^{-2}\,(K\,km\,s^{-1})^{-1}$. We have also investigated the distribution of "CO-dark" gas (DG) within the footprint of GAC and found a linear correlation between the DG column density and the $V$-band extinction: $N({\rm DG}) \simeq 2.2 \times 10^{21} (A_V - A^{c}_{V})\,\rm cm^{-2}$. The mass fraction of DG is found to be $f_{\rm DG}\sim 0.55$ toward the Galactic anticentre, which is respectively about 23 and 124 per cent of the atomic and CO-traced molecular gas in the same region. This result is consistent with the theoretical work of Papadopoulos et al. but much larger than that expected in the $\rm H_2$ cloud models by Wolfire et al.
Based on archive photographic photometry and recent CCD photometric data for red novae V4332 Sgr and V838 Mon, we established their stellar composition, exploded components, and the nature of explosions. Low temperature in the outburst maximum is due to quasi-adiabatic expansion of a massive stellar envelope after the central energy surge preceded the outburst.
A recent revision of black hole scaling relations (Kormendy & Ho 2013), indicates that the local mass density in black holes should be increased by up to a factor of five with respect to previously determined values. The local black hole mass density is connected to the mean radiative efficiency of accretion through the time integral of the AGN volume density and a significant increase of the local black holes mass density would have interesting consequences on AGN accretion properties and demography. One possibility to explain a large black hole mass density is that most of the Black Hole growth is via radiatively inefficient channels such as super Eddington accretion, however, given the intrinsic degeneracies in the Soltan argument, this solution is not unique. Here we show how it is possible to accommodate a larger fraction of heavily buried, Compton thick AGN, without violating the limit imposed by the hard X-ray and mid-infrared backgrounds spectral energy density.
Transits on single stars are rare. The probability rarely exceeds a few per cent. Furthermore, this probability rapidly approaches zero at increasing orbital period. Therefore transit surveys have been predominantly limited to the inner parts of exoplanetary systems. Here we demonstrate how circumbinary planets allow us to beat these unfavourable odds. By incorporating the geometry and the three-body dynamics of circumbinary systems, we analytically derive the probability of transitability, a configuration where the binary and planet orbits overlap on the sky. We later show that this is equivalent to the transit probability, but at an unspecified point in time. This probability, at its minimum, is always higher than for single star cases. In addition, it is an increasing function with mutual inclination. By applying our analytical development to eclipsing binaries, we deduce that transits are highly probable, and in some case guaranteed. For example, a circumbinary planet revolving at 1 AU around a 0.3 AU eclipsing binary is certain to eventually transit - a 100% probability - if its mutual inclination is greater than 0.6 deg. We show that the transit probability is generally only a weak function of the planet's orbital period; circumbinary planets may be used as practical tools for probing the outer regions of exoplanetary systems to search for and detect warm to cold transiting planets.
The rate at which interstellar gas is converted into stars, and its dependence on environment, is one of the pillars on which our understanding of the visible Universe is build. We present a comparison of the surface density of young stars (Sigma_*) and dust surface density (Sigma_d) across NGC346 (N66) in 115 independent pixels of 6x6 pc^2. We find a correlation between Sigma_* and Sigma_d with a considerable scatter. A power law fit to the data yields a steep relation with an exponent of 2.6+-0.2. We convert Sigma_d to gas surface density (Sigma_g) and Sigma_* to star formation rate (SFR) surface densities (Sigma_SFR), using simple assumptions for the gas-to-dust mass ratio and the duration of star formation. The derived total SFR (4+-1 10^-3 M_sun/yr) is consistent with SFR estimated from the Ha emission integrated over the Ha nebula. On small scales the Sigma_SFR derived using Ha systematically underestimates the count-based Sigma_SFR, by up to a factor of 10. This is due to ionizing photons escaping the area, where the stars are counted. We find that individual 36 pc^2 pixels fall systematically above integrated disc-galaxies in the Schmidt-Kennicutt diagram by on average a factor of ~7. The NGC346 average SFR over a larger area (90 pc radius) lies closer to the relation but remains high by a factor of ~3. The fraction of the total mass (gas plus young stars) locked in young stars is systematically high (~10 per cent) within the central 15 pc and systematically lower outside (2 per cent), which we interpret as variations in star formation efficiency. The inner 15 pc is dominated by young stars belonging to a centrally condensed cluster, while the outer parts are dominated by a dispersed population. Therefore, the observed trend could reflect a change of star formation efficiency between clustered and non-clustered star-formation.
We present a comprehensive analysis of the whole sample of available XMM-Newton observations of High Mass X-ray Binaries (HMXBs) until August, 2013, focusing on the FeK{\alpha} emission line. This line is a key tool to better understand the physical properties of the material surrounding the X-ray source within a few stellar radii (the circumstellar medium). We have collected observations from 46 HMXBs, detecting FeK{\alpha} in 21 of them. We have used the standard classification of HMXBs to divide the sample in different groups. We find that: (1) FeK{\alpha} is centred at a mean value of 6.42 keV. Considering the instrumental and fits uncertainties, this value is compatible with ionization states lower than FeXVIII. (2) The flux of the continuum is well correlated with the flux of the line, as expected. Eclipse observations show that the Fe fluorescence emission comes from an extended region surrounding the X-ray source. (3) FeK{\alpha} is narrow (width lower than 0.15keV), reflecting that the reprocessing material does not move at high speeds. We attempt to explain the broadness of the line in terms of three possible broadening phenomena: line blending, Compton scattering and Doppler shifts (with velocities of the reprocessing material V=1000-2000 km/s). (4) The equivalent hydrogen column (NH) directly correlates with the EW of FeK{\alpha}, displaying clear similarities to numerical simulations. It highlights the strong link between the absorbing and the fluorescent matter. The obtained results clearly point to a very important contribution of the donors wind in the FeK{\alpha} emission and the absorption when the donor is a supergiant massive star.
While astrochemical models are successful in reproducing many of the observed interstellar species, they have been struggling to explain the observed abundances of complex organic molecules. Current models tend to privilege grain surface over gas phase chemistry in their formation. One key assumption of those models is that radicals trapped in the grain mantles gain mobility and react on lukewarm (>30 K) dust grains. Thus, the recent detections of methyl formate (MF) and dimethyl ether (DME) in cold objects represent a challenge and may clarify the respective role of grain surface and gas phase chemistry. We propose here a new model to form DME and MF with gas phase reactions in cold environments, where DME is the precursor of MF via an efficient reaction overlooked by previous models. Furthermore, methoxy, a precursor of DME, is also synthetized in the gas phase from methanol, which is desorbed by a non-thermal process from the ices. Our new model reproduces fairy well the observations towards L1544. It also explains, in a natural way, the observed correlation between DME and MF. We conclude that gas phase reactions are major actors in the formation of MF, DME and methoxy in cold gas. This challenges the exclusive role of grain-surface chemistry and favours a combined grain-gas chemistry.
In the context of count-in-cells statistics, the joint probability distribution of the density in two concentric spherical shells is predicted from first first principle for sigmas of the order of one. The agreement with simulation is found to be excellent. This statistics allows us to deduce the conditional one dimensional probability distribution function of the slope within under dense (resp. overdense) regions, or of the density for positive or negative slopes. The former conditional distribution is likely to be more robust in constraining the cosmological parameters as the underlying dynamics is less evolved in such regions. A fiducial dark energy experiment is implemented on such counts derived from Lambda-CDM simulations.
[abridged] We use the latest release of CIGALE, a galaxy SED fitting model relying on energy balance, to study the influence of an AGN in estimating both the SFR and stellar mass in galaxies, as well as the contribution of the AGN to the power output of the host. Using the galaxy formation SAM GALFORM, we create mock galaxy SEDs using realistic star formation histories (SFH) and add an AGN of Type 1, Type 2, or intermediate type whose contribution to the bolometric luminosity can be variable. We perform an SED fitting of these catalogues with CIGALE assuming three different SFHs: a single- and double-exponentially-decreasing, and a delayed SFH. Constraining thecontribution of an AGN to the LIR (fracAGN) is very challenging for fracAGN<20%, with uncertainties of ~5-30% for higher fractions depending on the AGN type, while FIR and sub-mm are essential. The AGN power has an impact on the estimation of $M_*$ in Type 1 and intermediate type AGNs but has no effect for galaxies hosting Type 2 AGNs. We find that in the absence of AGN emission, the best estimates of $M_*$ are obtained using the double-exponentially-decreasing model but at the expense of realistic ages of the stellar population. The delayed SFH model provides good estimates of $M_*$ and SFR, with a maximum offset of 10% as well as better estimates of the age. Our analysis shows that the underestimation of the SFR increases with fracAGN for Type 1 systems, as well as for low contributions of an intermediate AGN type, but it is quite insensitive to the emission of Type 2 AGNs up to fracAGN~45%. Similarly the UV emission is critical in accurately retrieving the $M_*$ for Type 1 and intermediate type AGN, and the SFR of all of the three AGN types. We show that the presence of AGN emission introduces a scatter to the SFR-$M_*$ main sequence relation derived from SED fitting, which is driven by the uncertainties on $M_*$.
We present Herschel PACS mapping observations of the [OI]63 micron line towards protostellar outflows in the L1448, NGC1333-IRAS4, HH46, BHR71 and VLA1623 star forming regions. We detect emission spatially resolved along the outflow direction, which can be associated with a low excitation atomic jet. In the L1448-C, HH46 IRS and BHR71 IRS1 outflows this emission is kinematically resolved into blue- and red-shifted jet lobes, having radial velocities up to 200 km/s. In the L1448-C atomic jet the velocity increases with the distance from the protostar, similarly to what observed in the SiO jet associated with this source. This suggests that [OI] and molecular gas are kinematically connected and that this latter could represent the colder cocoon of a jet at higher excitation. Mass flux rates (\.M$_{jet}$(OI)) have been measured from the [OI]63micron luminosity adopting two independent methods. We find values in the range 1-4 10$^{-7}$ Mo/yr for all sources but HH46, for which an order of magnitude higher value is estimated. \.M$_{jet}$(OI) are compared with mass accretion rates (\.M$_{acc}$) onto the protostar and with \.M$_{jet}$ derived from ground-based CO observations. \.M$_{jet}$(OI)/\.M$_{acc}$ ratios are in the range 0.05-0.5, similar to the values for more evolved sources. \.M$_{jet}$(OI) in HH46 IRS and IRAS4A are comparable to \.M$_{jet}$(CO), while those of the remaining sources are significantly lower than the corresponding \.M$_{jet}$(CO). We speculate that for these three sources most of the mass flux is carried out by a molecular jet, while the warm atomic gas does not significantly contribute to the dynamics of the system.
The project Massive Unseen Companions to Hot Faint Underluminous Stars from SDSS (MUCHFUSS) aims at finding hot subdwarf stars with massive compact companions like massive white dwarfs (M>1.0 M$_\odot$), neutron stars, or stellar-mass black holes. We present orbital and atmospheric parameters and put constraints on the nature of the companions of 12 close hot subdwarf B star (sdB) binaries found in the course of the MUCHFUSS project. The systems show periods between 0.14 and 7.4 days. Three systems most likely have white dwarf companions. SDSS J083006.17+475150.3 is likely to be a rare example of a low-mass helium-core white dwarf. SDSS J095101.28+034757.0 shows an excess in the infrared that probably originates from a third companion in a wide orbit. SDSS J113241.58-063652.8 is the first helium deficient sdO star with a confirmed close companion. This study brings to 142 the number of sdB binaries with orbital periods of less than 30 days and with measured mass functions. We present an analysis of the minimum companion mass distribution and show that it is bimodal. One peak around 0.1 M$_\odot$ corresponds to the low-mass main sequence and substellar companions. The other peak around 0.4 M$_\odot$ corresponds to the white dwarf companions. The derived masses for the white dwarf companions are significantly lower than the average mass for single carbon-oxygen white dwarfs. In a T$_{\rm eff}$-log(g) diagram of sdB+dM companions, we find signs that the sdB components are more massive than the rest of the sample. The full sample was compared to the known population of extremely low-mass white dwarf binaries as well as short-period white dwarfs with main sequence companions. Both samples show a significantly different companion mass distribution. We calculate merger timescales and timescales when the companion will fill its Roche Lobe and the system evolves into a cataclysmic variable.
Molecular clouds at very high latitude ($b>60^{\circ}$) away from the Galactic plane are rare and in general are expected to be non-star-forming. However, we report the discovery of two embedded clusters (Camargo 438 and Camargo 439) within the high-latitude molecular cloud HRK 81.4-77.8 using WISE. Camargo 439 with Galactic coordinates $\ell=81.11^{\circ}$ and $b=-77.84^{\circ}$ is an $\sim2$ Myr embedded cluster (EC) located at a distance from the Sun of $d_{\odot}=5.09\pm0.47$ kpc. Adopting the distance of the Sun to the Galactic centre $R_{\odot}=7.2$ kpc we derive for Camargo 439 a Galactocentric distance of $R_{GC}=8.70\pm0.26$ kpc and a vertical distance from the plane of $-4.97\pm0.46$ kpc. Camargo 438 at $\ell=79.66^{\circ}$ and $b=-78.86^{\circ}$ presents similar values. The derived parameters for these two ECs put HRK 81.4-77.8 in the halo at a distance from the Galactic centre of $\sim8.7$ kpc and $\sim5.0$ kpc from the disc. Star clusters provide the only direct means to determine the high latitude molecular cloud distances. The present study shows that the molecular cloud HRK 81.4-77.8 is currently forming stars, apparently an unprecedented event detected so far among high latitude clouds. We carried out a preliminary orbit analysis. It shows that this ECs are the most distant known embedded clusters from the plane and both cloud and clusters are probably falling ballistically from the halo onto the Galactic disc, or performing a flyby.
In this paper, we study two sets of local geomagnetic indices from 26 stations using the principal component (PC) and the independent component (IC) analysis methods. We demonstrate that the annually averaged indices can be accurately represented as linear combinations of two first components with weights systematically depending on latitude. We show that the annual contributions of coronal mass ejections (CMEs) and high speed streams (HSSs) to geomagnetic activity are highly correlated with the first and second IC. The first and second ICs are also found to be very highly correlated with the strength of the interplanetary magnetic field (IMF) and the solar wind speed, respectively, because solar wind speed is the most important parameter driving geomagnetic activity during HSSs while IMF strength dominates during CMEs. These results help in better understanding the long-term driving of geomagnetic activity and in gaining information about the long-term evolution of solar wind parameters and the different solar wind structures.
We present the stellar population content of early-type galaxies from the
Atlas3D survey. Using spectra integrated within apertures covering up to one
effective radius, we apply two methods: one based on measuring line-strength
indices and applying single stellar population (SSP) models to derive
SSP-equivalent values of stellar age, metallicity, and alpha enhancement; and
one based on spectral fitting to derive non-parametric star-formation
histories, mass-weighted average values of age, metallicity, and half-mass
formation timescales. Using homogeneously derived effective radii and
dynamically-determined galaxy masses, we present the distribution of stellar
population parameters on the Mass Plane (M_JAM, Sigma_e, R_maj), showing that
at fixed mass, compact early-type galaxies are on average older, more
metal-rich, and more alpha-enhanced than their larger counterparts.
From non-parametric star-formation histories, we find that the duration of
star formation is systematically more extended in lower mass objects. Assuming
that our sample represents most of the stellar content of today's local
Universe, approximately 50% of all stars formed within the first 2 Gyr
following the big bang. Most of these stars reside today in the most massive
galaxies (>10^10.5 M_sun), which themselves formed 90% of their stars by z~2.
The lower-mass objects, in contrast, have formed barely half their stars in
this time interval. Stellar population properties are independent of
environment over two orders of magnitude in local density, varying only with
galaxy mass. In the highest-density regions of our volume (dominated by the
Virgo cluster), galaxies are older, alpha-enhanced and have shorter
star-formation histories with respect to lower density regions.
The Apache Point Observatory Galactic Evolution Experiment (APOGEE), part of the Sloan Digital Sky Survey III, explores the stellar populations of the Milky Way using the Sloan 2.5-m telescope linked to a high resolution (R~22,500), near-infrared (1.51-1.70 microns) spectrograph with 300 optical fibers. For over 100,000 predominantly red giant branch stars that APOGEE targeted across the Galactic bulge, disks and halo, the collected high S/N (>100 per half-resolution element) spectra provide accurate (~0.1 km/s) radial velocities, stellar atmospheric parameters, and precise (~0.1 dex) chemical abundances for about 15 chemical species. Here we describe the basic APOGEE data reduction software that reduces multiple 3D raw data cubes into calibrated, well-sampled, combined 1D spectra, as implemented for the SDSS-III/APOGEE data releases (DR10 and DR12). The processing of the near-IR spectral data of APOGEE presents some challenges for reduction, including automated sky subtraction and telluric correction over a 3 degree diameter field and the combination of spectrally dithered spectra. We also discuss areas for future improvement.
We present the discovery of five millisecond pulsars (MSPs) from the PALFA Galactic plane survey using Arecibo. Four of these (PSRs J0557+1551, J1850+0244, J1902+0300, and J1943+2210) are binary pulsars whose companions are likely white dwarfs, and one (PSR J1905+0453) is isolated. Phase-coherent timing solutions, ranging from $\sim$1 to $\sim$3 years in length, and based on observations from the Jodrell Bank and Arecibo telescopes, provide precise determinations of spin, orbital, and astrometric parameters. All five pulsars have large dispersion measures ($>100$ pc cm$^{-3}$, within the top 20% of all known Galactic field MSPs) and are faint (1.4 GHz flux density < 0.1 mJy, within the faintest 5% of all known Galactic field MSPs), illustrating PALFA's ability to find increasingly faint, distant MSPs in the Galactic plane. In particular, PSR J1850+0244 has a dispersion measure of 540 pc cm$^{-3}$, the highest of all known MSPs. Such distant, faint MSPs are important input for accurately modeling the total Galactic MSP population.
We investigate the dynamics of the FLRW flat cosmological models in which the vacuum energy varies with redshift. A particularly well motivated model of this type is the so-called quantum field vacuum, in which both kind of terms $H^{2}$ and constant appear in the effective dark energy density affecting the evolution of the main cosmological functions at the background and perturbation levels. Specifically, it turns out that the functional form of the quantum vacuum endows the vacuum energy of a mild dynamical evolution which could be observed nowadays and appears as dynamical dark energy. Interestingly, the low-energy behaviour is very close to the usual $\Lambda$CDM model, but it is by no means identical. Finally, within the framework of the quantum field vacuum we generalize the large scale structure properties, namely growth of matter perturbations, cluster number counts and spherical collapse model.
The young massive OB association Cygnus OB2, in the Cygnus X complex, is the closest (1400 pc) star forming region to the Sun hosting thousands of young low mass stars and up to 1000 OB stars, among which are some of the most massive stars known in our Galaxy. This region holds great importance for several fields of modern astrophysics, such as the study of the physical properties of massive and young low-mass stars and the feedback provided by massive stars on star and planet formation process. Cygnus OB2 has been recently observed with Chandra/ACIS-I as part of the 1.08Msec Chandra Cygnus OB2 Legacy Project. This survey detected 7924 X-ray sources in a square degree area centered on Cyg OB2. Since a proper classification and study of the observed X-ray sources also requires the analysis of their optical and infrared counterparts, we combined a large and deep set of optical and infrared catalogs available for this region with our new X-ray catalog. In this paper we describe the matching procedure and present the combined catalog containing 5703 sources. We also briefly discuss the nature of the X-ray sources with optical and infrared counterparts using their position in the color-magnitude and color-color diagrams.
We studied a sample of stars associated with the Sh 2-296 nebula, part of the
reflection nebulae complex in the region of Canis Major (CMa R1). Our sample
corresponds to optical counterparts of X-ray sources detected from observations
with the XMM-Newton satellite, which revealed dozens of possible low-mass young
stars not yet known in this region.
A sample of 58 young star candidates were selected based on optical spectral
features, mainly H{\alpha} and lithium lines, observed with multi-objects
spectroscopy performed by the Gemini South telescope. Among the candidates, we
find 41 confirmed T Tauri and 15 very likely young stars. Based on the
H{\alpha} emission, the T Tauri stars were distinguished between classical
(17%) and weak-lined (83%), but no significant difference was found in the age
and mass distribution of these two classes.
The characterization of the sample was complemented by near- and mid-infrared
data, providing an estimate of ages and masses from the comparison with
pre-main-sequence evolutionary models. While half of the young stars have an
age of 1-2 Myrs or less, only a small fraction (~25%) shows evidence of IR
excess revealing the presence of circumstellar discs. This low fraction is
quite rare compared to most young star-forming regions, suggesting that some
external factor has accelerated the disc dissipation.
We present a dense model grid with tailored input chemical composition appropriate for the Large Magellanic Cloud. We use a one-dimensional hydrodynamic stellar evolution code, which accounts for rotation, transport of angular momentum by magnetic fields, and stellar wind mass loss to compute our detailed models. We calculate stellar evolution models with initial masses of 70-500 Msun and with initial surface rotational velocities of 0-550 km/s, covering the core-hydrogen burning phase of evolution. We find our rapid rotators to be strongly influenced by rotationally induced mixing of helium, with quasi-chemically homogeneous evolution occurring for the fastest rotating models. Above 160 Msun, homogeneous evolution is also established through mass loss, producing pure helium stars at core hydrogen exhaustion independent of the initial rotation rate. Surface nitrogen enrichment is also found for slower rotators, even for stars that lose only a small fraction of their initial mass. For models above 150 MZAMS, and for models in the whole considered mass range later on, we find a considerable envelope inflation due to the proximity of these models to their Eddington limit. This leads to a maximum zero-age main sequence surface temperature of 56000 K, at 180 Msun, and to an evolution of stars in the mass range 50-100 Msun to the regime of luminous blue variables in the HR diagram with high internal Eddington factors. Inflation also leads to decreasing surface temperatures during the chemically homogeneous evolution of stars above 180 Msun. The cool surface temperatures due to the envelope inflation in our models lead to an enhanced mass loss, which prevents stars at LMC metallicity from evolving into pair-instability supernovae. The corresponding spin-down will also prevent very massive LMC stars to produce long-duration gamma-ray bursts, which might, however, originate from lower masses.
Small, cool planets represent the typical end-products of planetary formation. Studying the archi- tectures of these systems, measuring planet masses and radii, and observing these planets' atmospheres during transit directly informs theories of planet assembly, migration, and evolution. Here we report the discovery of three small planets orbiting a bright (Ks = 8.6 mag) M0 dwarf using data collected as part of K2, the new transit survey using the re-purposed Kepler spacecraft. Stellar spectroscopy and K2 photometry indicate that the system hosts three transiting planets with radii 1.5-2.1 R_Earth, straddling the transition region between rocky and increasingly volatile-dominated compositions. With orbital periods of 10-45 days the planets receive just 1.5-10x the flux incident on Earth, making these some of the coolest small planets known orbiting a nearby star; planet d is located near the inner edge of the system's habitable zone. The bright, low-mass star makes this system an excellent laboratory to determine the planets' masses via Doppler spectroscopy and to constrain their atmospheric compositions via transit spectroscopy. This discovery demonstrates the power of K2 and future space-based transit searches to find many fascinating objects of interest.
We report on the interannual variability of the atmospheric ice/dust cycle in the Martian polar regions for Mars Years 28-30. We used CRISM emission phase function measurements to derive atmospheric dust optical depths and data from the MARCI instrument to derive atmospheric water ice optical depths. We have used autocorrelation and cross correlation functions in order to quantify the degree to which dust and ice are correlated throughout both polar regions during Mars Years 28-29. We find that in the south polar region, dust has the tendency to "self clear", demonstrated by negative autocorrelation around the central peak. This does not occur in the north polar region. In the south polar region, dust and ice are temporally and spatially anti correlated. In the north polar region, this relationship is reversed, however temporal correlation of northern dust and ice clouds is weak - 6 times weaker than the anticorrelation in the south polar region. Our latitudinal autocorrelation functions allow us to put average spatial sizes of event cores and halos. Dust events in the south are largest, affecting almost the entire pole, whereas dust storms are smaller in the north. Ice clouds in north are similar in latitudinal extent to those in the south (both have halos < 10{\deg}). Using cross-correlation functions of water ice and dust, we find that dust events temporally lag ice events by 35-80 degrees of solar longitude in the north and south poles, which is likely due to seasonality of dust and ice events.
We propose a method to measure difference between positions of centroids of polarized and total flux of astronomical object in visible light, which we call for brevity polaroastrometry. Method efficiency is being demonstrated by the example of reduction of observations carried out with an instrument combining features of two-beam polarimeter with rotating half-wave plate and speckle interferometer, at 70-cm telescope in $V$ band. For the total number of accumulated photoelectrons $N_e=10^9$, which corresponds to series duration of 500 sec and object magnitude $V=6^m$, the method precision is 60-70 $\mu$as in terms of 1$\sigma$. At smaller $N_e$ precision decays as $\sim 1.7^{\prime\prime}/\sqrt{N_e}$, while at larger $N_e$ it remains the same due to imperfections of the half-wave plate. For the main sequence stars, non-polarized and polarized by interstellar dust, we didn't detect significant polaroastrometric signal. For the Mira variable star $\chi$ Cyg total amplitudes of the polaricenters were found to be $310\pm70$ $\mu$as and $300\pm70$ $\mu$as, for Stokes parameters $Q$ and $U$, respectively. For $o$ Cet these values are $490\pm100$ $\mu$as and $1160\pm100$ $\mu$as, for R Tri polaricenters coincide with the photocenter to the level of error.
In this work we consider the effect that the appearance of pseudoscalar condensates in a neutron star can have on its cooling rate. We make no particular assumption on the origin and characteristics of these possible condensates and only assume that in regions where the pseudoscalar density varies the propagation of photons is governed by modified Maxwell-Chern-Simons electrodynamics. We find that this gives non-trivial reflection coefficients between regions of different pseudoscalar density and may affect very substantially the star cooling rate. While quantitative results do depend on precise details that can only be answered once a proper equation of state is determined, the general trend is quite universal and serious consideration should be given to this possibility.
In this paper we investigate the limits imposed by thermodynamics to a dark energy fluid. We obtain the heat capacities and the compressibilities for a dark energy fluid. These thermodynamical variables are easily accessible experimentally for any terrestrial fluid. The thermal and mechanical stabilities require these quantities to be positive. We show that such requirements forbid the existence of a cosmic fluid with negative constant EoS parameter which excludes vacuum energy as a candidate to explain the cosmic acceleration. We also show that the current observational data from SN Ia, BAO and $H(z)$ are in conflict with the physical constraints that a general dark energy fluid with a time-dependent EoS parameter must obey which can be interpreted as an evidence against the dark energy hypothesis. Although our result excludes the vacuum energy, a geometrical cosmological term as originally introduced by Einstein in the field equations remains untouched.
Observations by the Fermi-LAT telescope have uncovered a significant $\gamma$-ray excess toward the Milky Way Galactic Center. There has been no detection of a similar signal in the direction of the Milky Way dwarf spheroidal galaxies. Additionally, astronomical observations indicate that dwarf galaxies and other faint galaxies are less dense than predicted by the simplest cold dark matter models. We show that a self-interacting dark matter model with a particle mass of roughly 50 GeV annihilating to the mediator responsible for the strong self-interaction can simultaneously explain all three observations. The mediator is necessarily unstable and its mass must be below about 100 MeV in order to lower densities in faint galaxies. If the mediator decays to electron-positron pairs with a cross section on the order of the thermal relic value, then we find that these pairs can up-scatter the interstellar radiation field and produce the observed $\gamma$-ray excess. We show that this model is compatible with all current constraints and highlight detectable signatures unique to self-interacting dark matter models.
During inflation, the geometry of spacetime is described by a (quasi-)de Sitter phase. Inflationary observables are determined by the underlying (softly broken) de Sitter isometry group SO(1, 4) which acts like a conformal group on R^3: when the fluctuations are on super-Hubble scales, the correlators of the scalar fields are constrained by conformal invariance. Heavy fields with mass m larger than the Hubble rate H correspond to operators with imaginary dimensions in the dual Euclidean three-dimensional conformal field theory. By making use of the dS/CFT correspondence we show that, besides the Boltzmann suppression expected from the thermal properties of de Sitter space, the generic effect of heavy fields in the inflationary correlators of the light fields is to introduce power-law suppressed corrections of the form O(H^2/m^2). This can be seen, for instance, at the level of the four-point correlator for which we provide the correction due to a massive scalar field exchange.
We study the impacts of reheating temperature on the inflationary predictions of the spectral index and tensor-to-scalar ratio. Assuming that reheating process is very fast, the reheating temperature can be constrained for sinusoidal oscillation within a factor of 10 - 100 or even better with the prospect of future observations. Beyond this, we find that the predictions can also be insensitive to the reheating temperature in certain models, including the Higgs inflation.
In the effective theory of isoscalar and isovector dark matter-nucleon interactions mediated by a heavy spin-1 or spin-0 particle, 8 isotope-dependent nuclear response functions can be generated in the dark matter scattering by nuclei. We compute the 8 nuclear response functions for the 16 most abundant elements in the Sun, i.e. H, $^{3}$He, $^{4}$He, $^{12}$C, $^{14}$N, $^{16}$O, $^{20}$Ne, $^{23}$Na, $^{24}$Mg, $^{27}$Al, $^{28}$Si, $^{32}$S, $^{40}$Ar, $^{40}$Ca, $^{56}$Fe, and $^{59}$Ni, through detailed numerical shell model calculations. We use our response functions to compute the rate of dark matter capture by the Sun for all isoscalar and isovector dark matter-nucleon effective interactions, including several operators previously considered for dark matter direct detection only. We study in detail the dependence of the capture rate on specific dark matter-nucleon interaction operators, and on the different elements in the Sun. We find that a so far neglected momentum dependent dark matter coupling to the nuclear vector charge gives a larger contribution to the capture rate than the constant spin-dependent interaction commonly included in experimental searches. Our investigation lays the foundations for model independent analyses of dark matter induced neutrino signals from the Sun. The nuclear response functions obtained in this study are listed in analytic form in an appendix, ready to be used in other projects.
Inhomogeneous plasmas and fluids contain energy stored in inhomogeneity and they naturally tend to relax into lower energy states by developing instabilities or by diffusion. But the actual amount of energy in such inhomogeneities has remained unknown. In the present work the amount of energy stored in a density gradient is calculated for several specific density profiles in a cylindric configuration. This is of practical importance for drift wave instability in various plasmas, and in particular in its application in models dealing with the heating of solar corona because the instability is accompanied with stochastic heating, so the energy contained in inhomogeneity is effectively transformed into heat. It is shown that even for a rather moderate increase of the density at the axis in magnetic structures in the corona by a factor 1.5 or 3, the amount of excess energy per unit volume stored in such a density gradient becomes several orders of magnitude greater than the amount of total energy losses per unit volume (per second) in quiet regions in the corona. Consequently, within the life-time of a magnetic structure such energy losses can easily be compensated by the stochastic drift wave heating.
The Advanced LIGO and Virgo experiments are poised to detect gravitational waves (GWs) directly for the first time this decade. The ultimate prize will be joint observation of a compact binary merger in both gravitational and electromagnetic channels. However, GW sky locations that are uncertain by hundreds of square degrees will pose a challenge. I describe a real-time detection pipeline and a rapid Bayesian parameter estimation code that will make it possible to search promptly for optical counterparts in Advanced LIGO. Having analyzed a comprehensive population of simulated GW sources, we describe the sky localization accuracy that the GW detector network will achieve as each detector comes online and progresses toward design sensitivity. Next, in preparation for the optical search with the intermediate Palomar Transient Factory (iPTF), we have developed a unique capability to detect optical afterglows of gamma-ray bursts (GRBs) detected by the Fermi Gamma-ray Burst Monitor (GBM). Its comparable error regions offer a close parallel to the Advanced LIGO problem, but Fermi's unique access to MeV-GeV photons and its near all-sky coverage may allow us to look at optical afterglows in a relatively unexplored part of the GRB parameter space. We present the discovery and broadband follow-up observations of eight GBM-iPTF afterglows. Two of the bursts are at low redshift, are sub-luminous with respect to "standard" cosmological bursts, and have spectroscopically confirmed broad-line type Ic supernovae. These two bursts are possibly consistent with mildly relativistic shocks breaking out from the progenitor envelopes rather than the standard mechanism of internal shocks within an ultra-relativistic jet. On a technical level, the GBM-iPTF effort is a prototype for locating and observing optical counterparts of GW events in Advanced LIGO with the Zwicky Transient Facility.
We present effective collision strengths for electron excitation and de-excitation of the ten forbidden transitions between the five lowest energy levels of the astronomically abundant doubly-ionised oxygen ion, O^{2+}. The raw collision strength data were obtained from an R-matrix intermediate coupling calculation using the Breit-Pauli relativistic approximation published previously by the authors. The effective collision strengths were calculated with kappa-distributed electron energies and are tabulated as a function of the electron temperature and kappa.
Links to: arXiv, form interface, find, astro-ph, recent, 1501, contact, help (Access key information)
