We present near-IR spectroscopy of red supergiant (RSG) stars in NGC 6822, obtained with the new VLT-KMOS instrument. From comparisons with model spectra in the J-band we determine the metallicity of 11 RSGs, finding a mean value of [Z] = -0.52 $\pm$ 0.21 which agrees well with previous abundance studies of young stars and HII regions. We also find an indication for a low-significance abundance gradient within the central 1 kpc. We compare our results with those derived from older stellar populations and investigate the difference using a simple chemical evolution model. By comparing the physical properties determined for RSGs in NGC 6822 with those derived using the same technique in the Galaxy and the Magellanic Clouds, we show that there appears to be no significant temperature variation of RSGs with respect to metallicity, in contrast with recent evolutionary models.
We study the relation of AGN accretion, star formation rate (SFR), and stellar mass (M$_*$) using a sample of $\approx$ 8600 star-forming galaxies up to z=2.5 selected with \textit{Herschel} imaging in the GOODS and COSMOS fields. For each of them we derive SFR and M$_*$, both corrected, when necessary, for emission from an active galactic nucleus (AGN), through the decomposition of their spectral energy distributions (SEDs). About 10 per cent of the sample are detected individually in \textit{Chandra} observations of the fields. For the rest of the sample we stack the X-ray maps to get average X-ray properties. After subtracting the X-ray luminosity expected from star formation and correcting for nuclear obscuration, we derive the average AGN accretion rate for both detected sources and stacks, as a function of M$_{*}$, SFR and redshift. The average accretion rate correlates with SFR and with M$_*$. The dependence on SFR becomes progressively more significant at z$>$0.8. This may suggest that SFR is the original driver of these correlations. We find that average AGN accretion and star formation increase in a similar fashion with offset from the star-forming "main-sequence". Our interpretation is that accretion onto the central black hole and star formation broadly trace each other, irrespective of whether the galaxy is evolving steadily on the main-sequence or bursting.
New spectroscopic surveys offer the promise of consistent stellar parameters and abundances ('stellar labels') for hundreds of thousands of stars in the Milky Way: this poses a formidable spectral modeling challenge. In many cases, there is a sub-set of reference objects for which the stellar labels are known with high(er) fidelity. We take advantage of this with The Cannon, a new data-driven approach for determining stellar labels from spectroscopic data. The Cannon learns from the 'known' labels of reference stars how the continuum-normalized spectra depend on these labels by fitting a flexible model at each wavelength; then, The Cannon uses this model to derive labels for the remaining survey stars. We illustrate The Cannon by training the model on only 543 stars in 19 clusters as reference objects, with Teff, log g and [Fe/H] as the labels, and then applying it to the spectra of 56,000 stars from APOGEE DR10. The Cannon is very accurate. Its stellar labels compare well to the stars for which APOGEE pipeline (ASPCAP) labels are provided in DR10, with rms differences that are basically identical to the stated ASPCAP uncertainties. Beyond the reference labels, The Cannon makes no use of stellar models nor any line-list, but needs a set of reference objects that span label-space. The Cannon performs well at lower signal-to-noise, as it delivers comparably good labels even at one ninth the APOGEE observing time. We discuss the limitations of The Cannon and its future potential, particularly, to bring different spectroscopic surveys onto a consistent scale of stellar labels.
Long Gamma-Ray Bursts (GRBs) are produced by ultra-relativistic jets launched from core collapse of massive stars. Most massive stars form in binaries and/or in star clusters, which means that there may be a significant external photon field (EPF) around the GRB progenitor. We calculate the inverse-Compton scattering of EPF by the hot electrons in the GRB jet. Three possible cases of EPFs are considered: the progenitor is (I) in a massive binary system, (II) surrounded by a Wolf-Rayet-star wind, and (III) in a dense star cluster. Typical luminosities of 10^47 - 10^50 erg/s in the 10 - 100 GeV band are expected, depending on the stellar luminosity, binary separation (I), wind mass loss rate (II), stellar number density (III), etc. We calculate the lightcurve and spectrum in each case, taking fully into account the equal-arrival time surfaces and possible pair-production absorption with the prompt gamma-rays. Observations can put constraints on the existence of such EPFs (and hence on the nature of GRB progenitors) and on the radius where the jet internal dissipation process accelerates electrons.
Calcium-rich supernovae (Ca-rich SNe) are peculiar low-luminosity SNe Ib with relatively strong Ca spectral lines at ~2 months after peak brightness. This class also has an extended projected offset distribution, with several members of the class offset from their host galaxies by 30 - 150 kpc. There is no indication of any stellar population at the SN positions. Using a sample of 13 Ca-rich SNe, we present kinematic evidence that the progenitors of Ca-rich SNe originate near the centers of their host galaxies and are kicked to the locations of the SN explosions. Specifically, SNe with small projected offsets have large line-of-sight velocity shifts as determined by nebular lines, while those with large projected offsets have no significant velocity shifts. Therefore, the velocity shifts must not be primarily the result of the SN explosion. There is an excess of SNe with blueshifted velocity shifts within two isophotal radii (5/6 SNe), indicating that the SNe are moving away from their host galaxies and redshifted SNe on the far sides of their galaxies are selectively missed in SN surveys. Additionally, nearly every Ca-rich SN is hosted by a galaxy with indications of a recent merger and/or is in a dense environment. We propose a progenitor model which fits all current data: The progenitor system for a Ca-rich SN is a double white dwarf (WD) system where at least one WD has a significant He abundance. This system, through an interaction with a super-massive black hole (SMBH) is ejected from its host galaxy and the binary is hardened, significantly reducing the merger time. After 10 - 100 Myr (on average), the system explodes with a large physical offset. The rate for such events is significantly enhanced for galaxies which have undergone recent mergers, potentially making Ca-rich SNe new probes of both the galaxy merger rate and (binary) SMBH population. (abridged)
We introduce a force correction term to better model the dynamical friction (DF) experienced by a supermassive black hole (SMBH) as it orbits within its host galaxy. This new approach accurately follows the orbital decay of a SMBH and drastically improves over commonly used advection methods. The force correction introduced here naturally scales with the force resolution of the simulation and converges as resolution is increased. In controlled experiments we show how the orbital decay of the SMBH closely follows analytical predictions when particle masses are significantly smaller than that of the SMBH. In a cosmological simulation of the assembly of a small galaxy, we show how our method allows for realistic black hole orbits. This approach overcomes the limitations of the advection scheme, where black holes are rapidly and artificially pushed toward the halo center and then forced to merge, regardless of their orbits. We find that SMBHs from merging dwarf galaxies can spend significant time away from the center of the remnant galaxy. Improving the modeling of SMBH orbital decay will help in making robust predictions of the growth, detectability, and merger rates of SMBHs, especially at low galaxy masses or at high redshift.
Protostellar outflows provide a means to probe the accretion process of forming stars and their ability to inject energy into their surroundings. However, conclusions based on outflow observations depend upon the degree of accuracy with which their properties can be estimated. We examine the quality of Atacama Large Millimeter/submillimeter Array (ALMA) observations of protostellar outflows by producing synthetic $^{12}$CO(1-0) and $^{13}$CO(1-0) observations of numerical simulations. We use various ALMA configurations, observational parameters, and outflow inclinations to assess how accurately different assumptions and setups can recover underlying properties. We find that more compact arrays and longer observing times can improve the mass and momentum recovery by a factor of two. During the first $\sim$0.3 Myr of evolution, $^{12}$CO(1-0) is optically thick, even for velocities $|v|\ge 1$ km s$^{-1}$, and outflow mass is severely underestimated without an optical depth correction. Likewise, $^{13}$CO(1-0) is optically thick during the first $\simeq 0.1$ Myr. However, underestimation due to shorter observing time, missing flux, and optical depth are partially offset by the assumption of local thermodynamic equilibrium and higher excitation temperatures. Overall, we expect that full ALMA $^{13}$CO(1-0) observations of protostellar sources within 500 pc with observing times $\gtrsim 1$ hrs and assumed excitation temperatures of $T<20$K will reliably measure mass and line-of-sight momentum to within 20\%.
Single stars in ancient globular clusters (GCs) are believed incapable of producing planetary nebulae (PNe), because their post-asymptotic-giant-branch evolutionary timescales are slower than the dissipation timescales for PNe. Nevertheless, four PNe are known in Galactic GCs. Their existence likely requires more exotic evolutionary channels, including stellar mergers and common-envelope binary interactions. I carried out a snapshot imaging search with the Hubble Space Telescope (HST) for PNe in bright Local Group GCs outside the Milky Way. I used a filter covering the 5007 A nebular emission line of [O III], and another one in the nearby continuum, to image 66 GCs. Inclusion of archival HST frames brought the total number of extragalactic GCs imaged at 5007 A to 75, whose total luminosity slightly exceeds that of the entire Galactic GC system. I found no convincing PNe in these clusters, aside from one PN in a young M31 cluster misclassified as a GC, and two PNe at such large angular separations from an M31 GC that membership is doubtful. In a ground-based spectroscopic survey of 274 old GCs in M31, Jacoby et al. (2013) found three candidate PNe. My HST images of one of them suggest that the [O III] emission actually arises from ambient interstellar medium rather than a PN; for the other two candidates, there are broad-band archival UV HST images that show bright, blue point sources that are probably the PNe. In a literature search, I also identified five further PN candidates lying near old GCs in M31, for which follow-up observations are necessary to confirm their membership. The rates of incidence of PNe are similar, and small but non-zero, throughout the GCs of the Local Group.
The magneto-rotational instability (MRI) is considered to be a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. In contrast to accretion disks, radial buoyancy driven by entropy and lepton fraction gradients is expected to have a dynamical role as important as rotation and shear. We investigate the poorly known impact of buoyancy on the non-linear phase of the MRI, by means of three dimensional numerical simulations of a local model in the equatorial plane of a protoneutron star. The use of the Boussinesq approximation allows us to utilise a shearing box model with clean shearing periodic boundary conditions, while taking into account the buoyancy driven by radial entropy and composition gradients. We find significantly stronger turbulence and magnetic fields in buoyantly unstable flows. On the other hand, buoyancy has only a limited impact on the strength of turbulence and magnetic field amplification for buoyantly stable flows in the presence of a realistic thermal diffusion. The properties of the turbulence are, however, significantly affected in the latter case. In particular, the toroidal components of the magnetic field and of the velocity become even more dominant with respect to the poloidal ones. Furthermore, we observed in the regime of stable buoyancy the formation of long lived coherent structures such as channel flows and zonal flows. Overall, our results support the ability of the MRI to amplify the magnetic field significantly even in stably stratified regions of protoneutron stars.
Expanding X-ray cavities observed in hot gas atmospheres of many galaxy groups and clusters generate shock waves and turbulence that are primary heating mechanisms required to avoid uninhibited radiatively cooling flows which are not observed. However, we show here that the evolution of buoyant cavities also stimulates radiative cooling of observable masses of low-temperature gas. During their early evolution, radiative cooling occurs in the wakes of buoyant cavities in two locations: in thin radial filaments parallel to the buoyant velocity and more broadly in gas compressed beneath rising cavities. Radiation from these sustained compressions removes entropy from the hot gas. Gas experiencing the largest entropy loss cools first, followed by gas with progressively less entropy loss. Most cooling occurs at late times, $\sim 10^8-10^9$ yrs, long after the X-ray cavities have disrupted and are impossible to detect. During these late times, slightly denser low entropy gas sinks slowly toward the centers of the hot atmospheres where it cools intermittently, forming clouds near the cluster center. Single cavities of energy $10^{57}-10^{58}$ ergs in the atmosphere of the NGC 5044 group create $10^8 - 10^9$ $M_{\odot}$ of cooled gas, exceeding the mass of extended molecular gas currently observed in that group. The cooled gas clouds we compute share many attributes with molecular clouds recently observed in NGC 5044 with ALMA: self-gravitationally unbound, dust-free, quasi-randomly distributed within a few kpc around the group center.
The first mechanism invoked to explain the existence of the thick disk in the Milky Way Galaxy, were the spiral arms. Up-to-date work summon several other possibilities that together seem to better explain this component of our Galaxy. All these processes must affect differently in distinct types of galaxies, but the contribution of each one has not been straightforward to quantify. In this work, we present a first comprehensive study of the effect of the spiral arms in the formation of thick disks, as going from early to late type disk galaxies, in an attempt to characterize and quantify this specific mechanism in galactic potentials. To this purpose, we perform numerical simulations of test particles in a three-dimensional spiral galaxy potential of normal spiral galaxies (from early to late types). By varying the parameters of the spiral arms we found that the vertical heating of the stellar disk becomes very important in some cases, and strongly depends on the galaxy morphology, pitch angle, arms mass and its pattern speed. The later the galaxy type, the larger is the effect on the disk heating. This study shows that the physical mechanism causing the vertical heating is different from simple resonant excitation. The spiral pattern induce chaotic behavior not linked necessarily to resonances but to direct scattering of disk stars, which leads to an increase of the velocity dispersion. We applied this study to the specific example of the Milky Way Galaxy, for which we have also added an experiment that includes the Galactic bar. From this study we deduce that the effect of spiral arms of a Milky-Way-like potential, on the dynamical vertical heating of the disk is negligible, unlike later galactic potentials for disks.
Recent work suggests blue ellipticals form in mergers and migrate quickly from the blue cloud of star-forming galaxies to the red sequence of passively evolving galaxies, perhaps as a result of black hole feedback. Such rapid reddening of stellar populations implies that large gas reservoirs in the pre-merger star-forming pair must be depleted on short time scales. Here we present pilot observations of atomic hydrogen gas in four blue early-type galaxies that reveal increasing spatial offsets between the gas reservoirs and the stellar components of the galaxies, with advancing post-starburst age. Emission line spectra show associated nuclear activity in two of the merged galaxies, and in one case radio lobes aligned with the displaced gas reservoir. These early results suggest that a kinetic process (possibly feedback from black hole activity) is driving the quick truncation of star formation in these systems, rather than a simple exhaustion of gas supply.
We clarify a relationship of the dynamics of a solar flare and a growing Coronal Mass Ejection (CME) by investigating the dynamics of magnetic fields during the X2.2-class flare taking place in the solar active region 11158 on 2011 February 15, based on simulation results obtained from Inoue et al. 2014. We found that the strongly twisted lines formed through the tether-cutting reconnection in the twisted lines of a nonlinear force-free field (NLFFF) can break the force balance within the magnetic field, resulting in their launch from the solar surface. We further discover that a large-scale flux tube is formed during the eruption as a result of the tether-cutting reconnection between the eruptive strongly twisted lines and these ambient weakly twisted lines. Then the newly formed large flux tube exceeds the critical height of the torus instability. The tether-cutting reconnection thus plays an important role in the triggering a CME. Furthermore, we found that the tangential fields at the solar surface illustrate different phases in the formation of the flux tube and its ascending phase over the threshold of the torus instability. We will discuss about these dynamics in detail.
Powerful winds driven by active galactic nuclei (AGN) are often invoked to play a fundamental role in the evolution of both supermassive black holes (SMBHs) and their host galaxies, quenching star formation and explaining the tight SMBH-galaxy relations. Recent observations of large-scale molecular outflows in ultra-luminous infrared galaxies (ULIRGs) have provided the evidence to support these studies, as they directly trace the gas out of which stars form. Theoretical models suggest an origin of these outflows as energy-conserving flows driven by fast AGN accretion disk winds. Previous claims of a connection between large-scale molecular outflows and AGN activity in ULIRGs were incomplete because they were lacking the detection of the putative inner wind. Conversely, studies of powerful AGN accretion disk winds to date have focused only on X-ray observations of local Seyferts and a few higher redshift quasars. Here we show the clear detection of a powerful AGN accretion disk wind with a mildly relativistic velocity of 0.25c in the X-ray spectrum of IRAS F11119+3257, a nearby (z = 0.189) optically classified type 1 ULIRG hosting a powerful molecular outflow. The AGN is responsible for ~80% of the emission, with a quasar-like luminosity of L_AGN = 1.5x10^46 erg/s. The energetics of these winds are consistent with the energy-conserving mechanism, which is the basis of the quasar mode feedback in AGN lacking powerful radio jets.
The Suzaku data of the highly variable magnetar 1E 1547.0$-$5408, obtained during the 2009 January activity, were reanalyzed. The 2.07 s pulsation of the 15--40 keV emission detected with the HXD was found to be phase modulated, with a period of $36.0^{+4.5}_{-2.5}$ ks and an amplitude of $0.52 \pm 0.14$ s. The modulation waveform is suggested to be more square-wave like rather than sinusoidal. While the effect was confirmed with the 10--14 keV XIS data, the modulation amplitude decreased towards lower energies, becoming consistent with 0 below 4 keV. After the case of 4U 0142+61, this makes the 2nd example of this kind of behavior detected from magnetars. The effect can be interpreted as a manifestation of free precession of this magnetar, which is suggested to be oblately deformed under the presence of strong toroidal field of $\sim 10^{16}$ G.
I show that the circumstellar matter (CSM) of the type Ia supernova 2014J is too massive and its momentum too large to be accounted for by any but the core-degenerate (CD) scenario for type Ia supernovae. Assuming the absorbing gas is of CSM origin, the several shells responsible of the absorption potassium lines are accounted for by a mass loss episode from a massive asymptotic giant branch star during a common envelope phase with a white dwarf companion. The time-varying potassium lines can be accounted for by ionization of neutral potassium and the Na-from-dust absorption (NaDA) model. Before explosion some of the potassium resides in the gas phase and some in dust. Weakening in absorption strength is caused by potassium-ionizing radiation of the supernova, while release of atomic potassium from dust increases the absorption. I conclude that if the absorbing gas originated from the progenitor of SN 2014J, then a common envelope phase took place about 15,000 years ago, leading to the merging of the core with the white dwarf companion, i.e., the core-degenerate scenario. Else, the absorbing material is of interstellar medium origin.
Measuring of the masses of galactic supermassive black holes (SMBHs) is an important task, since they correlate with the host galaxy properties and play an important role in evolution of galaxies. Here we present a new method for measuring of SMBH masses using the polarization of the broad lines emitted from active galactic nuclei (AGNs). We performed spectropolarometric observations of 9 AGNs and find that this method gives measured masses which are in a good agreement with reverberation measurements. An advantage of this method is that it can be used to measure the masses of SMBHs in a consistent way at different cosmological epochs.
We have constructed two types of analytical models for an isothermal filamentary cloud supported mainly by magnetic tension. The first one describes an isolated cloud while the second considers filamentary clouds spaced periodically. Both the models assume that the filamentary clouds are highly flattened. The former is proved to be the asymptotic limit of the latter in which each filamentary cloud is much thinner than the distance to the neighboring filaments. We show that these models reproduce main features of the 2D equilibrium model of Tomisaka (2014) for filamentary cloud threaded by perpendicular magnetic field. It is also shown that the critical mass to flux ratio is $ M /\Phi = (2 \pi \sqrt{G}) ^{-1} $, where $ M $, $ \Phi $ and $ G $ denote the cloud mass, the total magnetic flux of the cloud, and the gravitational constant, respectively. This upper bound coincides with that for an axisymmetric cloud supported by poloidal magnetic fields. We applied the variational principle for studying the Jeans instability of the first model. Our model cloud is unstable against fragmentation as well as the filamentary clouds threaded by longitudinal magnetic field. The fastest growing mode has a wavelength several times longer than the cloud diameter. The second model describes quasi-static evolution of filamentary molecular cloud by ambipolar diffusion.
Only a dozen confirmed FU Orionis-type young outbursting stars (FUors) are known today; this explains the interest in the recent FUor eruption of 2MASS J06593158-0405277. Its outburst and expected decline will be subject to numerous studies in the future. Almost equally important for the understanding of the eruption mechanism, however, is the physical characterization of the FUor's precursor. Here we analyze unpublished archival data and summarize - and partly revise - all relevant photometry from optical to submillimeter wavelengths. Our analysis implies that the FUor is possibly associated with eight T Tauri star candidates and a strong Class 0 source. Adopting a distance of 450 pc for the FUor, we derive a quiescent bolometric luminosity and temperature of L_bol = 4.8 L_Sun and T_bol = 1190 K, typical for young Class II sources. The central star has a temperature of T_eff = 4000 K, a mass of 0.75 M_Sun, and an age of about 6 x 10^5 yr. The SED implies a circumstellar mass of 0.01 - 0.06 M_Sun, and the system is surrounded by a faint infrared nebulosity. Our results provide an almost complete picture of a FUor progenitor, supporting the interpretation of future post-outburst studies.
Context: The north-west photo-dissociation region (PDR) in the reflection nebula NGC 7023 displays a complex structure. Filament-like condensations at the edge of the cloud can be traced via the emission of the main cooling lines, offering a great opportunity to study the link between the morphology and energetics of these regions. Aims: We study the spatial variation of the far-infrared fine-structure lines of [C II] (158 um) and [O I] (63 and 145 um). These lines trace the local gas conditions across the PDR. Methods: We used observations from the Herschel/PACS instrument to map the spatial distribution of these fine-structure lines. The observed region covers a square area of about 110" x 110" with an angular resolution that varies from 4" to 11". We compared this emission with ground-based and Spitzer observations of H2 lines, Herschel/SPIRE observations of CO lines, and Spitzer/IRAC 3.6 um images that trace the emission of polycyclic aromatic hydrocarbons. Results: The [C II] (158 um) and [O I] (63 and 145 um) lines arise from the warm cloud surface where the PDR is located and the gas is warm, cooling the region. We find that although the relative contribution to the cooling budget over the observed region is dominated by [O I]63 um (>30%), H2 contributes significantly in the PDR (35%), as does [C II]158 um outside the PDR (30%). Other species contribute little to the cooling ([O I]145 um 9%, and CO 4%). The [O I] maps resolve these condensations into two structures and show that the peak of [O I] is slightly displaced from the molecular H2 emission. The size of these structures is about 8" (0.015 pc) and in surface cover about 9% of the PDR emission. Finally, we did not detect emission from [N II]122 um, suggesting that the cavity is mostly filled with non-ionised gas.
The origin of elements made by the rapid neutron-capture process (r-process) is not fully understood. Different sources have been proposed, e.g., core-collapse supernovae and neutron star mergers. Old metal-poor stars carry the signature of the astrophysical r-process source(s). Europium is the most indicative element to trace the r-process production, since it is mostly made by the r-process and it is easy to observe compared to other heavy r-process elements. In this work we simulate the evolution of europium in our Galaxy with the inhomogeneous chemical evolution model ICE, and we compare our results with spectroscopic observations. We test the most important parameters affecting the chemical evolution of the r-process element Eu: (a) for neutron star mergers the coalescence time scale of the merger and the probability to experience a neutron star merger event after two supernova explosions occurred and formed a double neutron star system ) and (b) for the sub-class of magneto-rotationally driven Supernovae (Jet-SNe, their occurrence rate compared to standard supernovae ). The main results can be summarized as follows. The observed [Eu/Fe] pattern in the galaxy can be reproduced by a combination of neutron star mergers and magneto-rotationally driven supernovae as r-process sources, while neutron star mergers alone seem to set in at too high metallicities. Jet-SNe provide a cure for this deficiency at low metallicities. Furthermore, we confirm that local inhomogeneities can explain the observed large spread in the europium abundances at low metallicities. We also predict the evolution of [O/Fe] to test whether the spread in alpha-elements for in- homogeneous models agrees with observations and whether this provides constraints on supernova explosion models and their nucleosynthesis.
For the spectral analysis of high-resolution and high-signal-to-noise (S/N) spectra of hot stars, advanced non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These atmospheres are strongly dependent on the reliability of the atomic data that are used to calculate them. Reliable Ga IV - VI oscillator strengths are used to identify Ga lines in the spectra of the DA-type white dwarf G191-B2B and the DO-type white dwarf RE0503-289 and to determine their photospheric Ga abundances. We newly calculated Ga IV - VI oscillator strengths to consider their radiative and collisional bound-bound transitions in detail in our NLTE stellar-atmosphere models for analyzing of Ga lines exhibited in high-resolution and high-S/N UV observations of G191-B2B and RE0503-289. We unambiguously detected 20 isolated and 6 blended (with lines of other species) Ga V lines in the Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum of RE0503-289. The identification of Ga IV and Ga VI lines is uncertain because they are weak and partly blended by other lines. The determined Ga abundance is 3.5 +/- 0.5 x 10**-5 (mass fraction, about 625 times solar). The Ga IV / GA V ionization equilibrium, which is a very sensitive indicator for the effective temperature, is well reproduced in RE0503-289. We identified the strongest Ga IV lines (1258.801, 1338.129 A) in the HST/STIS (Hubble Space Telescope / Space Telescope Imaging Spectrograph) spectrum of G191-B2B and measured a Ga abundance of 2.0 +/- 0.5 x 10**-6 (about 22 times solar). Reliable measurements and calculations of atomic data are a prerequisite for stellar-atmosphere modeling. Observed Ga IV - V line profiles in two white dwarf (G191-B2B and RE0503-289) ultraviolet spectra were well reproduced with our newly calculated oscillator strengths. For the first time, this allowed us to determine the photospheric Ga abundance in white dwarfs.
Hot subdwarfs (sdBs) are core helium-burning stars, which lost almost their entire hydrogen envelope in the red-giant phase. Since a high fraction of those stars are in close binary systems, common envelope ejection is an important formation channel. We identified a total population of 51 close sdB+WD binaries based on time-resolved spectroscopy and multi-band photometry, derive the WD mass distribution and constrain the future evolution of these systems. Most WDs in those binaries have masses significantly below the average mass of single WDs and a high fraction of them might therefore have helium cores. We found 12 systems that will merge in less than a Hubble time and evolve to become either massive C/O WDs, AM\,CVn systems, RCrB stars or even explode as supernovae type Ia.
We study how halo intrinsic dynamical properties are linked to their formation processes, and how both are affected by the large scale tidal field in which halo resides. Halo merger trees obtained from cosmological N-body simulations are used to identify infall halos that are about to merge with their hosts. We find that local tidal field can significantly increase the tangential component of the infall velocity but on average do not change the radial component significantly. These results can be used to explain how the internal velocity anisotropy and spin of halos depend on environment. The position vectors and velocities of infall halos are aligned with the principal axes of the local tidal field, and the alignment depends on the strength of the tidal field. Opposite accretion patterns are found in weak and strong tidal fields, in the sense that in a weak field the accretion flow is dominated by radial motion within the local structure, while a large tangential component is present in a strong field. These findings can be used to understand the strong alignments we find between the principal axes of the internal velocity ellipsoids of halos and the local tidal field, and their dependence on the strength of tidal field. They also explain why halo spin increases with the strength of local tidal field, but only in weak tidal fields does the spin-tidal field alignment follow the prediction of the tidal torque theory. We discuss how our results may be used to understand the spins of disk galaxies and velocity structures of elliptical galaxies and their correlations with large-scale structure.
We present probability distribution functions (PDFs) of the surface densities of ionized and neutral gas in the nearby spiral galaxies M31 and M51, as well as of dust emission and extinction Av in M31. The PDFs are close to lognormal and those for HI and Av in M31 are nearly identical. However, the PDFs for H2 are wider than the HI PDFs and the M51 PDFs have larger dispersions than those for M31. We use a simple model to determine how the PDFs are changed by variations in the line-of-sight (LOS) pathlength L through the gas, telescope resolution and the volume filling factor of the gas, f_v. In each of these cases the dispersion sigma of the lognormal PDF depends on the variable with a negative power law. We also derive PDFs of mean LOS volume densities of gas components in M31 and M51. Combining these with the volume density PDFs for different components of the ISM in the Milky Way (MW), we find that sigma decreases with increasing length L with an exponent of -0.76 +/- 0.06, which is steeper than expected. We show that the difference is due to variations in f_v. As f_v is similar in M31, M51 and the MW, the density structure in the gas in these galaxies must be similar. Finally, we demonstrate that an increase in f_v with increasing distance to the Galactic plane explains the decrease in sigma with latitude of the PDFs of emission measure and FUV emission observed for the MW.
With the discovery of a high-energy neutrino flux in the 0.1 PeV to PeV range from beyond the Earth's atmosphere with the IceCube detector, neutrino astronomy has achieved a major breakthrough in the exploration of the high-energy universe. One of the main goals is the identification and investigation of the still mysterious sources of the cosmic rays which are observed at Earth with energies up to several $10^5$ PeV. In addition to being smoking-gun evidence for the presence of cosmic rays in a specific object, neutrinos escape even dense environments and can reach us from distant places in the universe, thereby providing us with a unique tool to explore cosmic accelerators. This article summarizes our knowledge about the observed astrophysical neutrino flux and current status of the search for individual cosmic neutrino sources. At the end, it gives an overview of plans for future neutrino telescope projects.
The Atmospheric Imaging Assembly in the Solar Dynamics Observatory provides full Sun images every 1 seconds in each of 7 Extreme Ultraviolet passbands. However, for a significant amount of these images, saturation affects their most intense core, preventing scientists from a full exploitation of their physical meaning. In this paper we describe a mathematical and automatic procedure for the recovery of information in the primary saturation region based on a correlation/inversion analysis of the diffraction pattern associated to the telescope observations. Further, we suggest an interpolation-based method for determining the image background that allows the recovery of information also in the region of secondary saturation (blooming).
Some young, massive stars can be found in the Galactic halo. As star formation is unlikely to occur in the halo, they must have been formed in the disk and been ejected shortly afterwards. One explanation is a supernova in a tight binary system. The companion is ejected and becomes a runaway star. HD\,271791 is the kinematically most extreme runaway star known (Galactic restframe velocity $725 \pm 195\, \rm km\,s^{-1}$, which is even larger than the Galactic escape velocity). Moreover, an analysis of the optical spectrum showed an enhancement of the $\alpha$-process elements. This indicates the capture of supernova ejecta, and therefore an origin in a core-collapse supernova. As such high space velocities are not reached by the runaway stars in classical binary supernova ejection scenarios, a very massive but compact primary, probably of Wolf-Rayet type is required. HD\,271791 is therefore a perfect candidate for studying nucleosynthesis in a supernova of probably type Ibc. The goal of this project is to determine the abundances of a large number of elements from the $\alpha$-process, the iron group, and heavier elements by a quantitative analysis of the optical and UV spectral range. Detailed line-formation calculations are employed that account for deviations from local thermodynamic equilibrium (non-LTE). We intend to verify whether core-collapse supernova are a site of r-process element production. Here, we state the current status of the project.
Surface differential rotation (DR) is one major ingredient of the magnetic field generation process in the Sun and likely in other stars. The term solar-like differential rotation describes the observation that solar equatorial regions rotate faster than polar ones. The opposite effect of polar regions rotating faster than equatorial ones (termed as antisolar DR) has only been observed in a few stars, although there is evidence from theoretical dynamo models. We present a new method to detect the sign of DR (i.e. solar-like or antisolar DR) by analyzing long-term high-precision light curves with the Lomb-Scargle periodogram.We compute the Lomb-Scargle periodogram and identify a set of significant periods $P_k$, which we associate with active regions located at different latitudes on the the stellar surface. If detectable, the first harmonics ($P_k'$) of these periods were identified to compute their peak-height-ratios $r_k:=h(P_k')/h(P_k)$. Spots rotating at lower latitudes generate less sine-shaped light curves, which requires additional power in the harmonics, and results in larger ratios $r_k$. Comparing different ratios $r_k$ and the associated periods $P_k$ yields information about the spot latitudes, and reveals the sign of DR. We tested our method on different sets of synthetic light curves all exhibiting solar-like DR. The number of cases where our method detects antisolar DR is the false-positive rate of our method. Depending on the set of light curves, the noise level, the required minimum peak separation, and the presence or absence of spot evolution, our method fails to detect the correct sign in at most 20%. We applied our method to 50 Kepler G stars and found 21-34 stars with solar-like DR and 5-10 stars with antisolar DR, depending on the minimum peak separation.
Hypervelocity stars are those that have speeds exceeding the escape speed and are hence unbound from the Milky Way. We investigate a sample of low-mass hypervelocity candidates obtained using data from the high-precision SDSS Stripe 82 catalogue, which we have combined with spectroscopy from the 200-inch Hale Telescope at Palomar Observatory. We find four good candidates, but without metallicities it is difficult to pin-down their distances and therefore total velocities. Our best candidate has a significant likelihood that it is escaping the Milky Way for a wide-range of metallicities.
At intervals as short as ten thousand years, each white dwarf (WD) passes within a solar radius of a planetoid, i.e., a comet, asteroid, or planet. Gravitational tidal forces tear the planetoid apart; its metal-rich debris falls onto the WD, enriching the atmosphere. A third of WDs exhibit atmospheric "pollution". For roughly every hundred planetoid disruptions, a planetoid collides with a WD. We simulate a small number of collisions, in which "death-by-dynamics" refers to the fate of the planetoid. We also compute the energies and likely durations of a broad sample of collision events, and identify detection strategies at optical and X-ray wavelengths. Collisions with the most massive planetoids can be detected in external galaxies. Some may trigger nuclear burning. If one in $\sim 10^7-10^8$ of WD-planetoid collisions creates the conditions needed for a Type Ia supernova (SN~Ia), "death-by-dynamics" would also refer to the fate of the WD, and could provide a novel channel for the production of SN~Ia. We consider the circumstances under which the rate of SNe~Ia can be increased by interactions with planetoids.
The source(s) of the neutrino excess reported by the IceCube Collaboration is unknown. The TANAMI Collaboration recently reported on the multiwavelength emission of six bright, variable blazars which are positionally coincident with two of the most energetic IceCube events. Such objects are prime candidates to be the source of the highest-energy cosmic rays, and thus of associated neutrino emission. We present an analysis of neutrino emission from the six blazars using observations with the ANTARES neutrino telescope. The standard methods of the ANTARES candidate list search are applied to six years of data to search for an excess of muons - and hence their neutrino progenitors - from the directions of the six blazars described by the TANAMI Collaboration, and which are possibly associated with two IceCube events. Monte Carlo simulations of the detector response to both signal and background particle fluxes are used to estimate the sensitivity of this analysis for different possible source neutrino spectra. A maximum-likelihood approach, using the reconstructed energies and arrival directions of through-going muons, is used to identify events with properties consistent with a blazar origin.Both blazars predicted to be the most neutrino-bright in the TANAMI sample (1653-329 and 1714-336) have a signal flux fitted by the likelihood analysis corresponding to approximately one event. This observation is consistent with the blazar-origin hypothesis of the IceCube event IC14 for a broad range of blazar spectra, although an atmospheric origin cannot be excluded. No ANTARES events are observed from any of the other four blazars, including the three associated with IceCube event IC20. This excludes at a 90% confidence level the possibility that this event was produced by these blazars unless the neutrino spectrum is flatter than -2.4.
Given the stellar density near the galactic center, close encounters between compact object binaries and the supermassive black hole are a plausible occurrence. We present results from a numerical study of close to 13 million such encounters. Consistent with previous studies, we corroborate that, for binary systems tidally disrupted by the black hole, the component of the binary remaining bound to the hole has eccentricity ~ 0.97 and circularizes dramatically by the time it enters the classical LISA band. Our results also show that the population of surviving binaries merits attention. These binary systems experience perturbations to their internal orbital parameters with potentially interesting observational consequences. We investigated the regions of parameter space for survival and estimated the distribution of orbital parameters post-encounter. We found that surviving binaries harden and their eccentricity increases, thus accelerating their merger due gravitational radiation emission and increasing the predicted merger rates by up to 1%.
The first discoveries of exoplanets around Sun-like stars have fueled efforts to find ever smaller worlds evocative of Earth and other terrestrial planets in the Solar System. While gas-giant planets appear to form preferentially around metal-rich stars, small planets (with radii less than four Earth radii) can form under a wide range of metallicities. This implies that small, including Earth-size, planets may have readily formed at earlier epochs in the Universe's history when metals were far less abundant. We report Kepler spacecraft observations of KOI-3158, a metal-poor Sun-like star from the old population of the Galactic thick disk, which hosts five planets with sizes between Mercury and Venus. We used asteroseismology to directly measure a precise age of 11.2+/-1.0 Gyr for the host star, indicating that KOI-3158 formed when the Universe was less than 20% of its current age and making it the oldest known system of terrestrial-size planets. We thus show that Earth-size planets have formed throughout most of the Universe's 13.8-billion-year history, providing scope for the existence of ancient life in the Galaxy.
Stellar dynamo processes can be explored by measuring the magnetic field. This is usually obtained using the atomic and molecular Zeeman effect in spectral lines. While the atomic Zeeman effect can only access warmer regions, the use of molecular lines is of advantage for studying cool objects. The molecules MgH, TiO, CaH, and FeH are suited to probe stellar magnetic fields, each one for a different range of spectral types, by considering the signal that is obtained from modeling various spectral types. We have analyzed the usefulness of different molecules (MgH, TiO, CaH, and FeH) as diagnostic tools for studying stellar magnetism on active G-K-M dwarfs. We investigate the temperature range in which the selected molecules can serve as indicators for magnetic fields on highly active cool stars and present synthetic Stokes profiles for the modeled spectral type. We modeled a star with a spot size of 10% of the stellar disk and a spot comprising either only longitudinal or only transverse magnetic fields and estimated the strengths of the polarization Stokes V and Q signals for the molecules MgH, TiO, CaH, and FeH. We combined various photosphere and spot models according to realistic scenarios. In G dwarfs, the molecules MgH and FeH show overall the strongest Stokes V and Q signals from the starspot, whereas FeH has a stronger Stokes V signal in all G dwarfs, with a spot temperature of 3800K. In K dwarfs, CaH signals are generally stronger, and the TiO signature is most prominent in M dwarfs. Modeling synthetic polarization signals from starspots for a range of G-K-M dwarfs leads to differences in the prominence of various molecular signatures in different wavelength regions, which helps to efficiently select targets and exposure times for observations.
The focus of this report is on the opportunities enabled by the combination of LSST, Euclid and WFIRST, the optical surveys that will be an essential part of the next decade's astronomy. The sum of these surveys has the potential to be significantly greater than the contributions of the individual parts. As is detailed in this report, the combination of these surveys should give us multi-wavelength high-resolution images of galaxies and broadband data covering much of the stellar energy spectrum. These stellar and galactic data have the potential of yielding new insights into topics ranging from the formation history of the Milky Way to the mass of the neutrino. However, enabling the astronomy community to fully exploit this multi-instrument data set is a challenging technical task: for much of the science, we will need to combine the photometry across multiple wavelengths with varying spectral and spatial resolution. We identify some of the key science enabled by the combined surveys and the key technical challenges in achieving the synergies.
Atmosphere is one of the most important noise sources for ground-based Cosmic Microwave Background (CMB) experiments. By increasing optical loading on the detectors, it amplifies their effective noise, while its fluctuations introduce spatial and temporal correlations between detected signals. We present a physically motivated 3d-model of the atmosphere total intensity emission in the millimeter and sub-millimeter wavelengths. We derive an analytical estimate for the correlation between detectors time-ordered data as a function of the instrument and survey design, as well as several atmospheric parameters such as wind, relative humidity, temperature and turbulence characteristics. Using numerical computation, we examine the effect of each physical parameter on the correlations in the time series of a given experiment. We then use a parametric-likelihood approach to validate the modeling and estimate atmosphere parameters from the POLARBEAR-I project first season data set. We compare our results to previous studies and weather station measurements, and find that the polarization fraction of atmospheric emission is below 1.0 percent. The proposed model can be used for realistic simulations of future ground-based CMB observations.
We show that the magnetic dipole and gravitational radiation emitted by a pulsar can undergo superradiant scattering off a spinning black hole companion. We find that the relative amount of superradiant modes in the radiation depends on the pulsar's angular position relative to the black hole's equatorial plane. In particular, when the pulsar and black hole spins are aligned, superradiant modes are dominant at large angles, leading to an amplification of the pulsar's luminosity, whereas for small angles the radiation is dominantly composed of non-superradiant modes and the signal is attenuated. This results in a characteristic double orbital modulation of the pulsar's luminosity that can be used to search for observational signatures of superradiant scattering in astrophysical black holes, providing an important test of general relativity.
From the dynamics of a brane-world with matter fields present in the bulk, the bulk metric and the black string solution near the brane are generalized, when both the dynamics of inhomogeneous dust/generalized dark radiation on the brane-world and inhomogeneous dark radiation in the bulk as well are considered -- as exact dynamical collapse solutions. Based on the analysis on the inhomogeneous static exterior of a collapsing sphere of homogeneous dark radiation on the brane, the associated black string warped horizon is studied, as well as the 5D bulk metric near the brane. Moreover, the black string and the bulk are shown to be more regular upon time evolution, for suitable values for the dark radiation parameter in the model, by analyzing the physical soft singularities.
Searches for the role of spin in gravitation dated before the firm establishment of the electron spin in 1925. Since mass and spin, or helicity in the case of zero mass, are the Casimir invariants of the Poincar\'e group and mass participates in universal gravitation, these searches are natural steps to pursue. In this update, we report on the progress on this topic in the last five years after our last review. We begin with how is Lorentz/Poincar\'e group in local physics arisen from spacetime structure as seen by photon and matter through experiments/observations. The cosmic verification of the Galileo Equivalence Principle for photons/electromagnetic wave packets (Universality of Propagation in spacetime independent of photon energy and polarization, i.e. nonbirefringence) constrains the spacetime constitutive tensor to high precision to a core metric form with an axion degree and a dilaton degree of freedom. Hughes-Drever-type experiments then constrain this core metric to agree with the matter metric. Thus comes the metric with axion and dilation. In local physics this metric gives the Lorentz/Poincar\'e covariance. Constraints on axion and dilaton from polarized/unpolarized laboratory/astrophysical/cosmic experiments/observations are presented. In the end, we review the theoretical progress on the issue of gyrogravitational ratio for fundamental particles and the experimental progress on the measurements of possible long range/intermediate range spin-spin, spin-monopole and spin-cosmos interactions.
Physical equivalence between different conformal frames in scalar-tensor theory of gravity is a known fact. However, assuming that matter minimally couples to the metric of a particular frame, which we call the matter Jordan frame, the matter point of view of the universe may vary from frame to frame. Thus, there is a clear distinction between gravitational sector (curvature and scalar field) and matter sector. In this paper, focusing on a simple power-law inflation model in the Einstein frame, two examples are considered; a super-inflationary and a bouncing universe Jordan frames. Then we consider a spectator curvaton minimally coupled to a Jordan frame, and compute its contribution to the curvature perturbation power spectrum. In these specific examples, we find a blue tilt at short scales for the super-inflationary case, and a blue tilt at large scales for the bouncing case.
We present Montblanc, a GPU implementation of the Radio interferometer
measurement equation (RIME) in support of the Bayesian inference for radio
observations (BIRO) technique. BIRO uses Bayesian inference to select sky
models that best match the visibilities observed by a radio interferometer. To
accomplish this, BIRO evaluates the RIME multiple times, varying sky model
parameters to produce multiple model visibilities. Chi-squared values computed
from the model and observed visibilities are used as likelihood values to drive
the Bayesian sampling process and select the best sky model.
As most of the elements of the RIME and chi-squared calculation are
independent of one another, they are highly amenable to parallel computation.
Additionally, Montblanc caters for iterative RIME evaluation to produce
multiple chi-squared values. Only modified model parameters are transferred to
the GPU between each iteration.
We implemented Montblanc as a Python package based upon NVIDIA's CUDA
architecture. As such, it is easy to extend and implement different pipelines.
At present, Montblanc supports point and Gaussian morphologies, but is designed
for easy addition of new source profiles. Montblanc's RIME implementation is
performant: On an NVIDIA K40, it is approximately 250 times faster than
MeqTrees on a dual hexacore Intel E5{2620v2 CPU. Compared to the OSKAR
simulator's GPU-implemented RIME components it is 7.7 and 12 times faster on
the same K40 for single and double-precision oating point respectively.
However, OSKAR's RIME implementation is more general than Montblanc's
BIRO-tailored RIME.
Theoretical analysis of Montblanc's dominant CUDA kernel suggests that it is
memory bound. In practice, profiling shows that is balanced between compute and
memory, as much of the data required by the problem is retained in L1 and L2
cache.
This paper gives a complete characterization of resonant orbits in a Kerr spacetime. A resonant orbit is defined as a geodesic for which the longitudinal and radial orbital frequencies are commensurate. Our analysis is based on expressing the resonance condition in its most symmetric form using Carlson's integrals. We provide a number of concise formulae for the dependence of resonances on the system parameters. Resonant effects may be observable during the in-spiral of a compact object into a super-massive black hole. When the slowly evolving orbital frequencies pass through a series of low-order resonances, rapid changes in the orbital parameters could produce measurable phase shifts in the emitted gravitational radiation (GW). Resonant orbits may also capture dust leading to electromagnetic emission. The KAM theorem indicates that, low order resonant orbits demarcate the regions where the onset of chaos could occur around a perturbed black-hole. We find that the 1/2 and 2/3 resonances occur at ~4 and 5.4 Schwarzschild radii (Rs) from the event horizon. For compact object in-spirals around super-massive black holes, this region lies within the sensitivity band of space-based GW detectors. For Sgr A*, length scales of ~41 and 55 microarcseconds and timescales of 50 and 79 min respectively should be associated with resonant effects, if Sgr A* is non-spinning. Spin decreases these values by up to ~32% and ~28%. These length-scales are potentially resolvable with VLBI measurements. We find that all low-order resonances are localized to the strong field region r < 50 Rs. This fact guarantees the validity of using approximations based on averaging to model the frequency evolution of a test object in region 50 Rs <r <1000 Rs. The systematic determination of the multipole moments of the central object by observing the orbit of a pulsar, free of chaotic effects, is thus possible.
We investigate the impact on the inflationary predictions from various reheating histories which are characterized by an e-folding number $N_{\mathrm{reh}}$ and an effective equation-of-state parameter $w_{\mathrm{reh}}$ during reheating process. For Higgs inflation with a non-minimal coupling to gravity, the predictions are obtained on the $N_{\mathrm{reh}}\!\!-\!w_{\mathrm{reh}}$ reheating phase diagram. We find that the predictions are insensitive to reheating phase. Within the $1\sigma$ region of the scalar spectral index $n_s$ reported by Planck 2014 Preliminary, almost all possible reheating histories are allowed on the reheating phase diagram, where Higgs inflation with canonical reheating history $w_{\mathrm{reh}}=0$ lies near the upper edge of the $1\sigma$ range of $n_s$. Future measurements of $n_s$ with high precision will identify the reheating physics of Higgs inflation.
We give the results of a study on the 222Rn decay we performed in the Gran Sasso Laboratory (LNGS) by detecting the gamma rays from the radon progeny. The motivation was to monitor the stability of radioactivity measuring several times per year the half-life of a short lifetime (days) source instead of measuring over a long period the activity of a long lifetime (tens or hundreds of years) source. In particular, we give the reason of the large periodical fluctuations in the count rate of the gamma rays due to radon inside a closed canister which has been described in literature and which has been attributed to a possible influence of a component in the solar irradiation affecting the nuclear decay rates. We then provide the result of four half-life measurements we performed underground at LNGS in the period from May 2014 to January 2015 with radon diffused into olive oil. Briefly, we did not measure any change of the 222Rn half-life with a 8*10^-5 precision. Finally, we provide the most precise value for the 222Rn half-life: 3.82146(16){stat}(4){syst} days.
In this work we study the theory of extended quasidilaton massive gravity together with the presence of matter fields. After discussing the homogeneous and isotropic fully dynamical background equations, which governs the exact expansion history of the universe, we consider small cosmological perturbations around these general FLRW solutions. The stability of tensor, vector and scalar perturbations on top of these general background solutions give rise to slightly different constraints on the parameters of the theory than those obtained in the approximative assumption of the late-time asymptotic form of the expansion history, which does not correspond to our current epoch. This opens up the possibility of stable FLRW solutions to be compared with current data on cosmic expansion with the restricted parameter space based on theoretical ground.
We here argue that the "knee" of the cosmic ray energy distribution at $E_c \sim 1$ PeV represents a second order phase transition of cosmic proportions. The discontinuity of the heat capacity per cosmic ray particle is given by $\Delta c=0.450196\ k_B$. However the idea of a deeper critical point singularity cannot be ruled out by present accuracy in neither theory nor experiment. The quantum phase transition consists of cosmic rays dominated by bosons for the low temperature phase E<E_c and dominated by fermions for high temperature phase $E > E_c$. The low temperature phase arises from those nuclei described by the usual and conventional collective boson models of nuclear physics. The high temperature phase is dominated by protons. The transition energy $E_c$ may be estimated in terms of the photo-disintegration of nuclei.
In a recent letter, the AMS collaboration reported the detailed and extensive data concerning the distribution in energy of electron and positron cosmic rays. A central result of the experimental work resides in the energy regime $30\ {\rm GeV} < E < 1\ {\rm TeV}$ wherein the power law exponent of the energy distribution is measured to be $\alpha ({\rm experiment})=3.17$. In virtue of the Fermi statistics obeyed by electrons and positrons, a theoretical value was predicted as $\alpha ({\rm theory})=3.151374$ in very good agreement with experimental data. The consequences of this agreement between theory and experiment concerning the sources of cosmic ray electrons and positrons are briefly explored.
We discuss the "constant speed of sound" (CSS) parameterization of the
equation of state of high density matter, and its application to the Field
Correlator Method (FCM) model of quark matter. We show how observational
constraints on the maximum mass and typical radius of neutron stars are
expressed as constraints on the CSS parameters. We find that the observation of
a $2\,{\rm M}_\odot$ star already severely constrains the CSS parameters, and
is particularly difficult to accommodate if the squared speed of sound in the
high density phase is assumed to be around $1/3$ or less.
We show that the FCM equation of state can be accurately represented by the
CSS parameterization. We display the mapping between the FCM and CSS
parameters, and see that FCM only allows equations of state in a restricted
subspace of the CSS parameters.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
We present new, late-time Hubble Space Telescope and Spitzer Space Telescope observations of the archetypal SN impostor SN 1997bs. We show that SN 1997bs remains much fainter than its progenitor, ruling out the canonical picture of late-time obscuration by dust forming in a shell ejected during the transient. The possibility that the star survived cloaked behind a dusty, steady wind is also disfavored. The simplest explanation is that SN 1997bs was a subluminous Type IIn SN, although it is impossible to rule out the possibility that the star survived, but with a significantly decreased intrinsic luminosity.
In comparison to gas and dust in star-forming galaxies at the peak epoch of galaxy assembly, which are presently the topic of intense study, little is known about the interstellar medium (ISM) of distant, passively evolving galaxies. We report on a deep 3 mm-band search with IRAM/PdBI for molecular gas in a massive ($M_{\star}{\sim}6{\times}10^{11}M_{\odot}$) elliptical galaxy at z=1.4277, the first observation of this kind ever attempted. We place a 3$\sigma$ upper limit of 0.30 Jy km/s on the flux of the CO($J$=$2\rightarrow$1) line or $L'_{\rm CO}$$<$8.3$\times$10$^{9}$ K km/s pc$^2$, assuming a line width in accordance with the stellar velocity dispersion of $\sigma_{\star}{\sim}330$ km/s. This translates to a molecular gas mass of $<$3.6$\times$10$^{10}$($\alpha_{\rm CO}$/4.4)$M_{\odot}$ or a gas fraction of $\lesssim$5% assuming a Salpeter initial mass function (IMF) and an ISM dominated by molecular gas, as observed in local early-type galaxies (ETGs). This low gas fraction approaches that of local ETGs, suggesting that the low star formation activity in massive, high-z passive galaxies reflects a true dearth of gas and a secondary role for inhibitive mechanisms like morphological quenching.
[Abridged] We have only been able to comprehensively characterize the atmospheres of a handful of transiting planets, because most orbit faint stars. TESS will discover transiting planets orbiting the brightest stars, enabling, in principle, an atmospheric survey of 10^2 to 10^3 bright hot Jupiters and warm sub-Neptunes. Uniform observations of such a statistically significant sample would provide leverage to understand---and learn from---the diversity of short-period planets. We argue that the best way to maximize the scientific returns of TESS is with a follow-up space mission consisting of a ~1 m telescope with an optical--NIR spectrograph: it could measure molecular absorption for non-terrestrial planets, as well as eclipses and phase variations for the hottest jovians. Such a mission could observe up to 10^3 transits per year, thus enabling it to survey a large fraction of the bright (J<11) TESS planets. JWST could be used to perform detailed atmospheric characterization of the most interesting transiting targets (transit, eclipse, and---when possible---phase-resolved spectroscopy). TESS is also expected to discover a few temperate terrestrial planets transiting nearby M-Dwarfs. Characterizing these worlds will be time-intensive: JWST will need months to provide tantalizing constraints on the presence of an atmosphere, planetary rotational state, clouds, and greenhouse gases. Future flagship missions should be designed to provide better constraints on the habitability of M-Dwarf temperate terrestrial planets.
We present kinematical analyses of 22 Galactic globular clusters using the Hubble Space Telescope proper motion (HSTPROMO) catalogues recently presented in Bellini et al. (2014). For most clusters, this is the first proper-motion study ever performed, and, for many, this is the most detailed kinematic study of any kind. We use cleaned samples of bright stars to determine binned velocity-dispersion and velocity-anisotropy radial profiles and two-dimensional velocity-dispersion spatial maps. Using these profiles, we search for correlations between cluster kinematics and structural properties. We find that: (1) more centrally-concentrated clusters have steeper radial velocity-dispersion profiles; (2) on average, at 1\sigma confidence in two dimensions, the photometric and kinematic centres of globular clusters agree to within ~1", with a cluster-to-cluster rms of 4" (including observational uncertainties); (3) on average, the cores of globular clusters have isotropic velocity distributions to within 1% (\sigma_t/\sigma_r = 0.992 +/- 0.005), with a cluster-to-cluster rms of 2% (including observational uncertainties); (4) clusters generally have mildly radially anisotropic velocity distributions (\sigma_t/\sigma_r ~ 0.8-1.0) near the half-mass radius, with bigger deviations from isotropy for clusters with longer relaxation times; (5) there is a relation between \sigma_minor/\sigma_major and ellipticity, such that the more flattened clusters in the sample tend to be more anisotropic, with \sigma_minor/\sigma_major ~ 0.9-1.0. Aside from these general results and correlations, the profiles and maps presented here can provide a basis for detailed dynamical modelling of individual globular clusters. Given the quality of the data, this is likely to provide new insights into a range of topics concerning globular cluster mass profiles, structure, and dynamics.
We develop a method to enable collaborative modelling of gravitational lenses and lens candidates, that could be used by non-professional lens enthusiasts. It uses an existing free-form modelling program (glass), but enables the input to this code to be provided in a novel way, via a user-generated diagram that is essentially a sketch of an arrival-time surface. We report on an implementation of this method, SpaghettiLens, which has been tested in a modelling challenge using 29 simulated lenses drawn from a larger set created for the Space Warps citizen science strong lens search. We find that volunteers from this online community asserted the image parities and time ordering consistently in some lenses, but made errors in other lenses depending on the image morphology. While errors in image parity and time ordering lead to large errors in the mass distribution, the enclosed mass was found to be more robust: the model-derived Einstein radii found by the volunteers were consistent with those produced by one of the professional team, suggesting that given the appropriate tools, gravitational lens modelling is a data analysis activity that can be crowd-sourced to good effect. Ideas for improvement are discussed, these include (a) overcoming the tendency of the models to be shallower than the correct answer in test cases, leading to systematic overestimation of the Einstein radius by 10 per cent at present, and (b) detailed modelling of arcs.
This study confirms MLS110213:022733+130617 as a new eclipsing polar. We performed optical spectroscopic, polarimetric and photometric follow-up of this variable source identified by the Catalina Real Time Transient Survey. Using the mid-eclipse times, we estimated an orbital period of 3.787 h, which is above the orbital period gap of the cataclysmic variables. There are nine other known polars with longer orbital periods, and only four of them are eclipsing. We identified high and low-brightness states and high polarization modulated with the orbital period. The spectra are typical of polars, with strong high ionization emission lines and inverted Balmer decrement. The HeII4686 A line is as strong as Hbeta. We modelled the photometric and polarimetric bright-state light curves using the CYCLOPS code. Our modelling suggests an extended emitting region on the WD surface, with a mean temperature of 9 keV and B in the range 18 -- 33 MG. The WD mass estimated from the shock temperature is 0.67 Msol. The distance was estimated as 406+- 54 pc using the Period-Luminosity-Colours method. MLS110213 populates a rare sub-group of polars, near the upper limit of the period distribution, important to understand the evolution of mCVs.
We explored the role of X-ray binaries composed by a black hole and a massive stellar companion (BHXs) as sources of kinetic feedback by using hydrodynamical cosmological simulations. Following previous results, our BHX model selects low metal-poor stars ($Z = [0,10^{-4}]$) as possible progenitors. The model that better reproduces observations assumes that a $\sim 20\%$ fraction of low-metallicity black holes are in binary systems which produce BHXs. These sources are estimated to deposit $\sim 10^{52}$ erg of kinetic energy per event. With these parameters and in the simulated volume, we find that the energy injected by BHXs represents $\sim 30\%$ of the total energy released by SNII and BHX events at redshift $z\sim7$ and then decreases rapidly as baryons get chemically enriched. Haloes with virial masses smaller than $\sim 10^{10} \,M_{\odot}$ (or $T_{\rm vir} \lesssim 10^5 $ K) are the most directly affected ones by BHX feedback. These haloes host galaxies with stellar masses in the range $10^7 - 10^8$ M$_\odot$. Our results show that BHX feedback is able to keep the interstellar medium warm, without removing a significant gas fraction, in agreement with previous analytical calculations. Consequently, the stellar-to-dark matter mass ratio is better reproduced at high redshift. Our model also predicts a stronger evolution of the number of galaxies as a function of the stellar mass with redshift when BHX feedback is considered. These findings support previous claims that the BHXs could be an effective source of feedback in early stages of galaxy evolution.
We predict Lyman-$\alpha$ (Ly$\alpha$) luminosity functions (LFs) of Ly$\alpha$-selected galaxies (Ly$\alpha$ emitters, or LAEs) at $z=3-6$ using the phenomenological model of Dijkstra & Wyithe (2012). This model combines observed UV-LFs of Lyman-break galaxies (LBGs, or drop out galaxies), with constraints on their distribution of Ly$\alpha$ line strengths as a function of UV-luminosity and redshift. Our analysis shows that while Ly$\alpha$ LFs of LAEs are generally not Schechter functions, these provide a good description over the luminosity range of $\log_{10}( L_{\alpha}/{\rm erg}\,{\rm s}^{-1})=41-44$. Motivated by this result, we predict Schechter function parameters at $z=3-6$. Our analysis further shows that (i) the faint end slope of the Ly$\alpha$ LF is steeper than that of the UV-LF of Lyman-break galaxies, (with a median $\alpha_{Ly\alpha} < -2.0$ at $z\gtrsim 4$), and (ii) a turn-over in the Ly$\alpha$ LF of LAEs at Ly$\alpha$ luminosities $10^{40}$ erg s$^{-1}<L_{\alpha}\lesssim 10^{41}$ erg s$^{-1}$ may signal a flattening of UV-LF of Lyman-break galaxies at $-12>M_{\rm UV}>-14$. We discuss the implications of these results - which can be tested directly with upcoming surveys - for the Epoch of Reionization.
Stellar winds observed in asymptotic giant branch (AGB) stars are usually
attributed to a combination of stellar pulsations and radiation pressure on
dust. Shock waves triggered by pulsations propagate through the atmosphere,
compressing the gas and lifting it to cooler regions, which create favourable
conditions for grain growth. If sufficient radiative acceleration is exerted on
the newly formed grains through absorption or scattering of stellar photons, an
outflow can be triggered. Strong candidates for wind-driving dust species in
M-type AGB stars are magnesium silicates (Mg$_2$SiO$_4$ and MgSiO$_3$). Such
grains can form close to the stellar surface, they consist of abundant
materials and, if they grow to sizes comparable to the wavelength of the
stellar flux maximum, they experience strong acceleration by photon scattering.
We use a frequency-dependent radiation-hydrodynamics code with a detailed
description for the growth of Mg$_2$SiO$_4$ grains to calculate the first
extensive set of time-dependent wind models for M-type AGB stars. The resulting
wind properties, visual and near-IR photometry and mid-IR spectra are compared
with observations.We show that the models can produce outflows for a wide range
of stellar parameters. We also demonstrate that they reproduce observed
mass-loss rates and wind velocities, as well as visual and near-IR photometry.
However, the current models do not show the characteristic silicate features at
10 and 18 $\mu$m as a result of the cool temperature of Mg$_2$SiO$_4$ grains in
the wind. Including a small amount of Fe in the grains further out in the
circumstellar envelope will increase the grain temperature and result in
pronounced silicate features, without significantly affecting the photometry in
the visual and near-IR wavelength regions.
We present the results of our observations of the early stages of the 2012--2013 outburst of the transient black hole X-ray binary (BHXRB), Swift J1745$-$26, with the VLA, SMA, and JCMT (SCUBA--2). Our data mark the first multiple-band mm & sub-mm observations of a BHXRB. During our observations the system was in the hard accretion state producing a steady, compact jet. The unique combination of radio and mm/sub-mm data allows us to directly measure the spectral indices in and between the radio and mm/sub-mm regimes, including the first mm/sub-mm spectral index measured for a BHXRB. Spectral fitting revealed that both the mm (230 GHz) and sub-mm (350 GHz) measurements are consistent with extrapolations of an inverted power-law from contemporaneous radio data (1--30 GHz). This indicates that, as standard jet models predict, a power-law extending up to mm/sub-mm frequencies can adequately describe the spectrum, and suggests that the mechanism driving spectral inversion could be responsible for the high mm/sub-mm fluxes (compared to radio fluxes) observed in outbursting BHXRBs. While this power-law is also consistent with contemporaneous optical data, the optical data could arise from either jet emission with a jet spectral break frequency of $\nu_{{\rm break}}\gtrsim1\times10^{14}\,{\rm Hz}$ or the combination of jet emission with a lower jet spectral break frequency of $\nu_{{\rm break}}\gtrsim2\times10^{11}\,{\rm Hz}$ and accretion disc emission. Our analysis solidifies the importance of the mm/sub-mm regime in bridging the crucial gap between radio and IR frequencies in the jet spectrum, and justifies the need to explore this regime further.
Young Galactic supernova remnants (SNRs) are where we can observe closely the supernova (SN) ejecta and its interaction with circumstellar/interstellar medium. Therefore, they provide an opportunity to explore the explosion and the final stage of the evolution of massive stars. Near-infrared (NIR) emission lines in SNRs mostly originate from shocked dense material. In shocked SN ejecta, forbidden lines from heavy ions are prominent, while in shocked circumstellar/interstellar medium, [Fe II] and H2 lines are prominent. [Fe II] lines are strong in both media, and therefore [Fe II] line images provide a good starting point for the NIR study of SNRs. There are about twenty SNRs detected in [Fe II] lines, some of which have been studied in NIR spectroscopy. We will review the NIR [Fe II] observations of SNRs and introduce our recent NIR spectroscopic study of the young core-collapse SNR Cas A where we detected strong [P II] lines.
The disk instability picture gives a plausible explanation for the behavior of soft X-ray transient systems if self-irradiation of the disk is included. We show that there is a simple relation between the peak luminosity (at the start of an outburst) and the decay timescale. We use this relation to place constraints on systems assumed to undergo disk instabilities. The observable X-ray populations of elliptical galaxies must largely consist of long-lived transients, as deduced on different grounds by Piro and Bildsten (2002). The strongly-varying X-ray source HLX-1 in the galaxy ESO 243-49 can be modeled as a disk instability of a highly super-Eddington stellar-mass binary similar to SS433. A fit to the disk instability picture is not possible for an intermediate-mass black hole model for HLX-1. Other, recently identified, super-Eddington ULXs might be subject to disk instability.
Young open clusters located in the outer Galaxy provide us with an opportunity to study star formation activity in a different environment from the solar neighborhood. We present a UBVI and H alpha photometric study of the young open clusters NGC 1624 and NGC 1931 that are situated toward the Galactic anticenter. Various photometric diagrams are used to select the members of the clusters and to determine the fundamental parameters. NGC 1624 and NGC 1931 are, on average, reddened by <E(B-V)> = 0.92 +/- 0.05 and 0.74 +/- 0.17 mag, respectively. The properties of the reddening toward NGC 1931 indicate an abnormal reddening law (Rv,cl = 5.2 +/- 0.3). Using the zero-age main sequence fitting method we confirm that NGC 1624 is 6.0 +/- 0.6 kpc away from the Sun, whereas NGC 1931 is at a distance of 2.3 +/- 0.2 kpc. The results from isochrone fitting in the Hertzsprung-Russell diagram indicate the ages of NGC 1624 and NGC 1931 to be less than 4 Myr and 1.5 - 2.0 Myr, respectively. We derived the initial mass function (IMF) of the clusters. The slope of the IMF (Gamma_NGC 1624 = -2.0 +/- 0.2 and Gamma_NGC 1931 = -2.0 +/- 0.1) appears to be steeper than that of the Salpeter/Kroupa IMF. We discuss the implication of the derived IMF based on simple Monte-Carlo simulations and conclude that the property of star formation in the clusters seems not to be far different from that in the solar neighborhood.
We show that warm dark matter keV fermions (`inos') can be responsible for both core and halo galactic structure, in agreement with current astrophysical/cosmological constraints. We identify the inos with sterile right-handed neutrinos. The possible mass range of up to a few tens of keV, obtained independently from the galactic structure and dark matter astroparticle physics, points towards an important role of the right-handed neutrinos in the cosmic structure.
In our 2013 Astronomical Review article, we discussed the statistics of variability for 633 faint spectral type A-F stars observed by the Kepler spacecraft during Quarters 6-13. We found six stars that showed no variability with amplitude 20 ppm or greater in the range 0.2 to 24.4 cycles/day, but whose positions in the log g--Teff diagram place them in the delta Sct or gamma Dor pulsation instability regions established from pre-Kepler ground-based observations. Here we present results for 2137 additional stars observed during Quarters 14-17, and find 34 stars that lie within the instability regions. In Paper I, we included a +229 K offset to the Kepler Input Catalog Teff to take into account an average systematic difference between the KIC values and the Teff derived from SDSS color photometry for main-sequence F stars (Pinsonneault et al. 2012). Here we compare the KIC Teff value and the Teff derived from spectroscopy taken by the LAMOST instrument (Molenda-Zakowicz et al. 2013, 2014) for 54 stars common to both samples. We find no trend to support applying the offset, but instead find that a small average temperature decrease relative to the KIC Teff may be more appropriate for the stars in our spectral-type range. If the offset is omitted, only 17 of our 34 `constant' stars fall within the instability regions. For the two `constant' stars also observed by LAMOST, the LAMOST Teff values are cooler than the KIC Teff by several hundred K, and would move these stars out of the instability regions. It is possible that a more accurate determination of their Teff and log g would move some of the other `constant' stars out of the instability regions. However, if average (random) errors in Teff are taken into account, 15 to 52 stars may still persist within the instability regions. Explanations for these `constant' stars, both theoretical and observational, remain to be investigated.
We study the warping and tearing of a geometrically thin, non-self-gravitating disk surrounding binary supermassive black holes on an eccentric orbit. The circumbinary disk is significantly misaligned with the binary orbital plane, and is subject to the time-dependent tidal torques. In principle, such a disk is warped and precesses, and is torn into mutually misaligned rings in the region, where the tidal precession torques are stronger than the local viscous torques. We derive the tidal-warp and tearing radii of the misaligned circumbinary disks around eccentric SMBH binaries. We find that in disks with the viscosity parameter, alpha, larger than a critical value depending on the disk aspect ratio, the disk warping appears outside the tearing radius. This condition is expressed as alpha > sqrt{H/3r} for H/r ~<0.1, where H is the disk scale height. If alpha < sqrt{H/3r}, only the disk tearing occurs because the tidal warp radius is inside the tearing radius, where most of disk material is likely to rapidly accrete onto SMBHs. In warped and torn disks, both the tidal-warp and the tearing radii most strongly depend on the binary semi-major axis, although they also mildly depend on the other orbital and disk parameters. This strong dependence enables us to estimate the semi-major axis, once the tidal warp or tearing radius is determined observationally: For the tidal warp radius of 0.1 pc, the semi-major axis is estimated to be ~10^{-2} pc for 10^7 Msun black hole with typical orbital and disk parameters. We also briefly discuss the possibility that central objects of observed warped maser disks in active galactic nuclei are supermassive black hole binaries.
We study the GMC environments surrounding 10 IRDCs, based on 13CO molecular line emission from the Galactic Ring Survey. Using a range of physical scales, we measure the physical properties of the IRDCs and their surrounding molecular material extending out to radii, R, of 30pc. By comparing different methods for defining cloud boundaries and for deriving mass surface densities, Sigma, and velocity dispersions, sigma, we settled on a preferred "CE,tau,G" method of "Connected Extraction" in position-velocity space along with Gaussian fitting to opacity-corrected line profiles for velocity dispersion and mass estimation. We examine how cloud definition affects measurements of the magnitude and direction of line of sight velocity gradients and velocity dispersions, including the associated dependencies on size scale. CE,tau,G-defined IRDCs and GMCs show velocity gradient versus size relations that scale approximately as dv_0/ds~s^(-1/2) and velocity dispersion versus size relations sigma~s^(1/2), which are consistent with the large scale gradients being caused by turbulence. Interpreting velocity gradients as due to rotation, we find a broad spread in rotation directions with respect to Galactic rotation and rotation to gravitational energy fractions beta~0.1. We examine the dynamical state of the clouds finding mean virial parameters alpha_vir~2, consistent with models of magnetized virialized pressure-confined polytropic clouds. "CE,tau,G" IRDCs and GMCs exhibit a strikingly tight correlation of sigma/R^(1/2)~Sigma^n, with n~0.5, the value expected for virial equilibrium. We conclude this is strong evidence of cloud virialization over a wide range of scales from IRDCs to their parent GMCs, perhaps representing a self-similar hierarchy of self-gravitating virialized structures, which is the initial dynamical state of gas that is likely to form star clusters.
The insoluble organic matter (IOM) of an unequilibrated enstatite chondrite Sahara (SAH) 97096 has been investigated using a battery of analytical techniques. As the enstatite chondrites are thought to have formed in a reduced environment at higher temperatures than carbonaceous chondrites, they constitute an interesting comparative material to test the heterogeneities of the IOM in the solar system and to constrain the processes that could affect IOM during solar system evolution. The SAH 97096 IOM is found in situ: as submicrometer grains in the network of fine-grained matrix occurring mostly around chondrules and as inclusions in metallic nodules, where the carbonaceous matter appears to be more graphitized. IOM in these two settings has very similar $\delta^{15}N$ and $\delta^{13}C$; this supports the idea that graphitized inclusions in metal could be formed by metal catalytic graphitization of matrix IOM. A detailed comparison between the IOM extracted from a fresh part and a terrestrially weathered part of SAH 97096 shows the similarity between both IOM samples in spite of the high degree of mineral alteration in the latter. The isolated IOM exhibits a heterogeneous polyaromatic macromolecular structure, sometimes highly graphitized, without any detectable free radicals and deuterium-heterogeneity and having mean H- and N-isotopic compositions in the range of values observed for carbonaceous chondrites. It contains some submicrometer-sized areas highly enriched in $^{15}N$ ($\delta^{15}N$ up to 1600 permil). These observations reinforce the idea that the IOM found in carbonaceous chondrites is a common component widespread in the solar system. Most of the features of SAH 97096 IOM could be explained by the thermal modification of this main component.
In this paper we present the discussion on the salient points of the computational analysis that are at the basis of the paper \emph{Rotation Curves of Galaxies by Fourth Order Gravity} \citep{StSc}. In fact in this paper any galactic component (bulge, disk and Dark matter component) required an onerous numerical computation since the Gauss theorem is not applicable in the Fourth Order Gravity. The computational and data analysis have been made with the software Mathematica$^\circledR$.
We have developed a fast numerical 2-D model of galaxy disk evolution (resolved along the galaxy radius and azimuth) by adopting a scheme of parameterized stochastic self-propagating star formation. We explore the parameter space of the model and demonstrate its capability to reproduce 1-D radial profiles of the galaxy M33: gas surface density, surface brightness in the i and GALEX FUV passbands, and metallicity.
Observations of stable mainly dipolar magnetic fields at the surface of ~7% of single hot stars indicate that these fields are of fossil origin, i.e. they descend from the seed field in the molecular clouds from which the stars were formed. Recent results confort this theory. First, theoretical work and numerical simulations confirm that the properties of the observed fields correspond to those expected from fossil fields. They also showed that rapid rotation does not modify the surfacic dipolar magnetic configurations, but hinders the stability of fossil fields. This explains the lack of correlation between the magnetic field properties and stellar properties in massive stars. It may also explain the lack of detections of magnetic fields in Be stars, which rotate close to their break-up velocity. In addition, observations by the BinaMIcS collaboration of hot stars in binary systems show that the fraction of those hosting detectable magnetic fields is much smaller than for single hot stars. This could be related to results obtained in simulations of massive star formation, which show that the stronger the magnetic field in the original molecular cloud, the more difficult it is to fragment massive cores to form several stars. Therefore, more and more arguments support the fossil field theory.
Our recent search for the presence of a magnetic field in the bright early A-type supergiant HD92207 using FORS2 in spectropolarimetric mode revealed the presence of a longitudinal magnetic field of the order of a few hundred Gauss. However, the definite confirmation of the magnetic nature of this object remained pending due to the detection of short-term spectral variability probably affecting the position of line profiles in left- and right-hand polarized spectra. We present new magnetic field measurements of HD92207 obtained on three different epochs in 2013 and 2014 using FORS2 in spectropolarimetric mode. A 3sigma detection of the mean longitudinal magnetic field using the entire spectrum, <B_z>_all=104+-34G, was achieved in observations obtained in 2014 January. At this epoch, the position of the spectral lines appeared stable. Our analysis of spectral line shapes recorded in opposite circularly polarized light, i.e. in light with opposite sense of rotation, reveals that line profiles in the light polarized in a certain direction appear slightly split. The mechanism causing such a behaviour in the circularly polarized light is currently unknown. Trying to settle the issue of short-term variability, we searched for changes in the spectral line profiles on a time scale of 8-10min using HARPS polarimetric spectra and on a time scale of 3-4min using time series obtained with the CORALIE spectrograph. No significant variability was detected on these time scales during the epochs studied.
We investigate the effects of stochasticity on the observed galaxy parameters by comparing our stochastic star formation two-dimensional (2-D) galaxy evolution models with the commonly used one-dimensional (1-D) models with smooth star formation. The 2-D stochastic models predict high variability of the star formation rate and the surface photometric parameters across the galactic disks and in time.
Neutral Helium multiplets, HeI*3189,3889,10830 are very useful diagnostics to the geometry and physical conditions of the absorbing gas in quasars. So far only a handful of HeI* detections have been reported. Using a newly developed method, we detected HeI*3889 absorption line in 101 sources of a well-defined sample of 285 MgII BAL quasars selected from the SDSS DR5. This has increased the number of HeI* BAL quasars by more than one order of magnitude. We further detected HeI*3189 in 50% (52/101) quasars in the sample. The detection fraction of HeI* BALs in MgII BAL quasars is about 35% as a whole, and increases dramatically with increasing spectral signal-to-noise ratios, from 18% at S/N <= 10 to 93% at S/N >= 35. This suggests that HeI* BALs could be detected in most MgII LoBAL quasars, provided spectra S/N is high enough. Such a surprisingly high HeI* BAL fraction is actually predicted from photo-ionization calculations based on a simple BAL model. The result indicates that HeI* absorption lines can be used to search for BAL quasars at low-z, which cannot be identified by ground-based optical spectroscopic survey with commonly seen UV absorption lines. Using HeI*3889, we discovered 19 BAL quasars at z<0.3 from available SDSS spectral database. The fraction of HeI* BAL quasars is similar to that of LoBAL objects.
The whole set of INTEGRAL observations of type Ia supernova SN2014J, covering the period 16-162 days after the explosion has being analyzed. For spectral fitting the data are split into early and late periods covering days 16-35 and 50-162, respectively, optimized for Ni-56 and Co-56 lines. As expected for the early period much of the gamma-ray signal is confined to energies below $\sim$200 keV, while for the late period it is most strong above 400 keV. In particular, in the late period Co-56 lines at 847 and 1248 keV are detected at 4.7 and 4.3 sigma respectively. The lightcurves in several representative energy bands are calculated for the entire period. The resulting spectra and lightcurves are compared with a subset of models. We confirm our previous finding that the gamma-ray data are broadly consistent with the expectations for canonical 1D models, such as delayed detonation or deflagration models for a near-Chandrasekhar mass WD. Late optical spectra (day 136 after the explosion) show rather symmetric Co and Fe lines profiles, suggesting that unless the viewing angle is special, the distribution of radioactive elements is symmetric in the ejecta.
The knowledge of the intrinsic three-dimensional (3D) structure of galaxy components provides crucial information about the physical processes driving their formation and evolution. In this paper I discuss the main developments and results in the quest to better understand the 3D shape of galaxy bulges. I start by establishing the basic geometrical description of the problem. Our understanding of the intrinsic shape of elliptical galaxies and galaxy discs is then presented in a historical context, in order to place the role that the 3D structure of bulges play in the broader picture of galaxy evolution. Our current view on the 3D shape of the Milky Way bulge and future prospects in the field are also depicted.
We used 15 GHz VLBA observations of 366 sources having at least 5 epochs within a time interval 1995-2013 from the MOJAVE program and/or its predecessor, the 2 cm VLBA Survey. For each source we produced a corresponding stacked image averaging all available epochs for a better reconstruction of the cross section of the flow. We have analyzed jet profiles transverse to the local jet ridge line and derived both apparent and intrinsic opening angles of the parsec-scale outflows. The sources detected by the Fermi Large Area Telescope (LAT) during the first 24 months of operation show wider apparent jet opening angle and smaller viewing angles on a very high level of significance supporting our early findings. Analyzing transverse shapes of the outflows we found that most sources have conical jet geometry at parsec scales, though there are also sources that exhibit active jet collimation.
We identify some of the most HI massive and fastest rotating disk galaxies in the local universe with the aim of probing the processes that drive the formation of these extreme disk galaxies. By combining data from the Cosmic Flows project, which has consistently reanalyzed archival galaxy HI profiles, and 3.6$\mu$m photometry obtained with the Spitzer Space Telescope, with which we can measure stellar mass, we use the baryonic Tully-Fisher (BTF) relationship to explore whether these massive galaxies are distinct. We discuss several results, but the most striking is the systematic offset of the HI-massive sample above the BTF. These galaxies have both more gas and more stars in their disks than the typical disk galaxy of similar rotational velocity. The "condensed" baryon fraction, $f_C$, the fraction of the baryons in a dark matter halo that settle either as cold gas or stars into the disk, is twice as high in the HI-massive sample than typical, and almost reaches the universal baryon fraction in some cases, suggesting that the most extreme of these galaxies have little in the way of a hot baryonic component or cold baryons distributed well outside the disk. In contrast, the star formation efficiency, measured as the ratio of the mass in stars to that in both stars and gas, shows no difference between the HI-massive sample and the typical disk galaxies. We conclude that the star formation efficiency is driven by an internal, self-regulating process, while $f_C$ is affected by external factors. We also found that the most massive HI detected galaxies are located preferentially in filaments. We present the first evidence of an environmental effect on galaxy evolution using a dynamical definition of a filament.
Observations of the solar chromosphere in the line-core of the \Halpha\ line show dark elongated structures called fibrils that show swaying motion. We performed a 3D radiation-MHD simulation of a network region, and computed synthetic \Halpha\ images from this simulation to investigate the relation between fibrils and the magnetic field lines in the chromosphere. The periods, amplitudes and phase-speeds of the synthetic fibrils are consistent with those observed. We analyse the relation between the synthetic fibrils and the field lines threading through them, and find that some fibrils trace out the same field line along the fibril's length, but there are also fibrils that sample different field lines at different locations along their length. Fibrils sample the same field lines on a time scale of $\sim200$~s. This is shorter than their own lifetime. We analysed the evolution of the atmosphere along a number of field lines that thread through fibrils and find that they carry slow-mode waves that load mass into the field line, as well as transverse waves that propagate with the Alfv\'en speed. Transverse waves propagating in opposite directions cause an interference pattern with complex apparent phase speeds. The relationship between fibrils and field lines is complex. It is governed by constant migration and swaying of the field lines, their mass loading by slow modes and subsequent draining, and their actual visibility in \Halpha. Field lines are visible where they lie close to the optical depth unity surface. The location of the latter is governed by the height at which the column mass in the chromosphere reaches a certain value. We conclude that using the swaying motion of fibrils as a tracer of chromospheric transverse oscillations must be done with caution.
We present the first attempt to use a combination of CMB, LIGO, and PPTA data to constrain both the tilt and the running of primordial tensor power spectrum through constraints on the gravitational wave energy density generated in the early universe. Combining measurements at different cosmological scales highlights how complementary data can be used to test the predictions of early universe models including the inflationary consistency relation. Current data prefers a slightly positive tilt ($n_t = 0.13^{+0.54}_{-0.75}$) and a negative running ($n_{t, {\rm run}} < -0.25$) for the tensor power spectrum spectrum. Interestingly, the addition of direct gravitational wave detector data puts strong bounds on the tensor-to-scalar ratio $r < 0.2 $ since the large positive tensor tilt preferred by the Planck temperature power spectrum is no longer allowed. We comment on possible effects of a large positive tilt on the background expansion and show that depending on the assumptions regarding the UV cutoff ($k_{\rm UV}/k_* = 10^{24}$ Mpc$^{-1}$) of the primordial spectrum of gravitational waves, the strongest (upper) bounds on $n_t = 0.05^{+0.44}_{-0.9}$ are derived from this effect.
The production of energetic particles in the universe remains one of the great mysteries of modern science. The mechanisms of acceleration in astrophysical sources and the details about the propagation through the galactic and extragalactic media are still to be defined. In recent years, the cosmic ray flux has been measured with high precision in the energy range from \energy{10} to \energyEV{20.5} by several experiments using different techniques. In some energy ranges, it has been possible to determine the flux of individual elements (hydrogen to iron nuclei). This paper explores an astrophysical scenario in which only our Galaxy and the radio galaxy Cen A produce all particles measured on Earth in the energy range from \energy{10} to \energyEV{20.5}. Data from AMS-02, CREAM, KASCADE, KASCADE-Grande and the Pierre Auger Observatories are considered. The model developed here is able to describe the total and individual particle flux of all experiments considered. It is shown that the theory used here is able to describe the smooth transition from space-based to ground-based measurements. The flux of each element as determined by KASCADE and KASCADE-Grande and the mass sensitivity parameter \xmax measured by the Pierre Auger Observatory above \energyEV{18} are also explored within the framework of the model. The transition from \energy{16} to \energyEV{18} is carefully analyzed. It is shown that the data measured in this energy range suggest the existence of an extra component of cosmic rays yet to be understood.
We discuss moderate resolution spectra, multicolor photometry, and light curves of thirty-one of the most luminous stars and variables in the giant spiral M101. The majority are intermediate A to F-type supergiants. We present new photometry and light curves for three known "irregular blue variables" V2, V4 and V9) and identify a new candidate. Their spectra and variability confirm that they are LBV candidates and V9 may be in an LBV-like maximum light state or eruption.
The Warkworth Radio Astronomical Observatory is operated by the Institute for Radio Astronomy and Space Research (IRASR), AUT University, Auckland, New Zealand. Here we review the characteristics of the VLBI station facilities and report on a number of activities and technical developments in 2014.
We present a new measurement of the mass-concentration relation and the stellar-to-halo mass ratio over a 5*10^(12) solar mass to 2*10^(14) solar mass range. To achieve this, we use the CFHT Stripe 82 Survey (CS82) weak lensing data combined with a well defined catalog of clusters (the redMaPPer catalogue) and the LOWZ/CMASS galaxies of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release (SDSS-III BOSS DR10). The stacked lensing signals around these samples are modeled as a sum of contributions from the central galaxy, the dark matter halo, and the neighboring halos. We measure the mass-concentration relation: c200(M)=A(M200/M0)^(B) with A=5.25+/-1.67, B=-0.13+/-0.12 for 0.2<z<0.4 and A=6.77+/-1.13, B=-0.15+/-0.06 for 0.4<z<0.6. We conclude that the amplitude A and slope B are both consistent with the simulation predictions by Klypin et al. (2014) within the errors. We also measure the stellar-to-halo mass ratio and find it to be flatter than previous measurement for high stellar masses because of the complex structures and merger history in massive dark matter halos.
This paper presents a solely thermal flare, which we detected in the microwave range from the thermal gyro- and free-free emission it produced. An advantage of analyzing thermal gyro emission is its unique ability to precisely yield the magnetic field in the radiating volume. When combined with observationally-deduced plasma density and temperature, these magnetic field measurements offer a straightforward way of tracking evolution of the magnetic and thermal energies in the flare. For the event described here, the magnetic energy density in the radio-emitting volume declines over the flare rise phase, then stays roughly constant during the extended peak phase, but recovers to the original level over the decay phase. At the stage where the magnetic energy density decreases, the thermal energy density increases; however, this increase is insufficient, by roughly an order of magnitude, to compensate for the magnetic energy decrease. When the magnetic energy release is over, the source parameters come back to nearly their original values. We discuss possible scenarios to explain this behavior.
We study the radio--far infrared (FIR) correlation in "blue cloud" galaxies chosen from the PRism MUltiobject Survey (PRIMUS) up to redshift ($z$) of 1.2 in the XMM-LSS field. We use rest-frame emission at 1.4 GHz in the radio and both monochromatic (at 70$\mu$m) and bolometric (between $8-1000~\mu$m) emission in the FIR. To probe the nature of the correlation up to $z\sim1.2$, where direct detection of blue star-forming galaxies is impossible with current technology, we employ the technique of image stacking at 0.325 and 1.4 GHz in the radio and in six infrared bands, viz. 24, 70, 160, 250, 350 and $500~\mu$m. For comparison, we also study the correlation for more luminous galaxies that are directly detected. The stacking analysis allows us to probe the radio--FIR correlation for galaxies that are up to 2 orders of magnitude fainter than the ones detected directly. The $k-$correction in the infrared wavebands is obtained by fitting the observed spectral energy distribution (SED) with a composite mid-IR power law and a single temperature greybody model. We find that the radio luminosity at 1.4 GHz ($L_{\rm 1.4GHz}$) is strongly correlated with monochromatic FIR luminosity at 70 $\mu$m ($L_{\rm 70\mu m}$) having slope $1.09\pm0.05$ and with bolometric luminosity ($L_{\rm TIR}$) having slope $1.11\pm0.04$. The quantity $q_{\rm TIR} (=\log_{10}[L_{\rm TIR}/(3.75\times 10^{12} L_{\rm 1.4 GHz})])$ is observed to decrease with redshift as $q_{\rm TIR} \propto (1+z)^{-0.16\pm0.03}$ probably caused due to the non-linear slope of the radio--FIR correlation. Within the uncertainties of our measurement and the limitations of our flux-limited and color-selected sample, we do not find any evolution of the radio--FIR correlation with redshift.
We employed field star decontaminated 2MASS photometry to study four nearby optical embedded clusters -- Collinder34, NGC3293, NGC3766 and NGC6231 -- obtaining deep colour-magnitude diagrams and stellar radial density profiles. We found what seems to be pre-main sequences detached in different amounts from main sequences in these diagrams. The structural analysis of each cluster revealed different radial distributions for these two sequences. We argued that the detached evolutionary sequences in our sample may be evidence of sequential star formation. Finally, we compared the sample cluster parameters with those of other young clusters in the literature and point out evidence that NGC3766 and NGC6231 might be evolving to OB associations.
We report on a long-term single-dish and VLBI monitoring for intermittent flare activities of a Dominant Blue-Shifted H$_{2}$O Maser (DBSM) associated with a southern high mass young stellar object, G353.273+0.641. Bi-weekly single-dish monitoring using Hokkaido University Tomakomai 11-m radio telescope has shown that a systematic acceleration continues over four years beyond a lifetime of individual maser features. This fact suggests that the H$_{2}$O maser traces a region where molecular gas is steadily accelerated. There were five maser flares during five-years monitoring, and maser distributions in four of them were densely monitored by the VLBI Exploration of Radio Astrometry (VERA). The overall distribution of the maser features suggests the presence of a bipolar jet, with the 3D kinematics indicating that it is almost face-on (inclination angle of $\sim$ 8$^{\fdg}$--17$^{\fdg}$ from the line-of-sight). Most of maser features were recurrently excited within a region of 100$\times$100 AU$^{2}$ around the radio continuum peak, while their spatial distributions significantly varied between each flare. This confirms that episodic propagations of outflow shocks recurrently invoke intermittent flare activities. We also measured annual parallax, deriving the source distance of 1.70 $^{+0.19}_{-0.16}$ kpc that is consistent with the commonly-used photometric distance.
The condensation of complex silicates with pyroxene and olivine composition at conditions prevailing in molecular clouds has been experimentally studied. For this purpose, molecular species comprising refractory elements were forced to accrete on cold substrates representing the cold surfaces of surviving dust grains in the interstellar medium. The efficient formation of amorphous and homogeneous magnesium iron silicates at temperatures of about 12 K has been monitored by IR spectroscopy. The gaseous precursors of such condensation processes in the interstellar medium are formed by erosion of dust grains in supernova shock waves. In the laboratory, we have evaporated glassy silicate dust analogs and embedded the released species in neon ice matrices that have been studied spectroscopically to identify the molecular precursors of the condensing solid silicates. A sound coincidence between the 10 micron band of the interstellar silicates and the 10 micron band of the low-temperature siliceous condensates can be noted.
We present a semi-analytic, physically motivated model for dark matter halo concentration as a function of halo mass and redshift. The semi-analytic model is intimately based on hierarchical structure formation. It uses an analytic model for the halo mass accretion history, based on extended Press Schechter (EPS) theory, and an empirical relation between concentration and an appropriate definition of formation time obtained through fits to the results of numerical simulations. The resulting concentration-mass relations are tested against the simulations and do not exhibit an upturn at high masses or high redshifts as claimed by recent works. Because our semi-analytic model is based on EPS theory, it can be applied to wide ranges in mass, redshift and cosmology. We predict a change of slope in the z=0 concentration-mass relation at a mass scale of $10^{11}\rm{M}_{\odot}$, that is caused by the varying power in the density perturbations. We provide best-fitting expressions of the $c-M$ relations as well as numerical routines. We investigate how halo mass accretion histories affect the evolution of concentrations, finding that the decrease in the accretion rate during the dark energy epoch, produced by the accelerated expansion of the Universe, allows dark matter halos to virialize, relax and contract, and thus concentrations to grow. We also analyzed how the concentration-mass relation predicted by this work affects the power produced by dark matter annihilation.
We use X-ray and optical microlensing measurements of 47 image pairs in 18 lens systems to study the shape of the dark matter density profile in the lens galaxies and the size of the (soft) X-ray emission region. We show that single epoch X-ray microlensing is sensitive to the source size. Our results, in good agreement with previous estimates, show that the X-ray size scales roughly linearly with the black hole mass, with a half light radius of $R_{1/2}\simeq(20\pm12) r_g$ (for $r_g=GM_{BH}/c^2$). This corresponds to a size of $\sim$ 1 light day for a black hole mass of $M_{BH}=10^9 M_\sun$. We simultaneously estimated the fraction of the local surface mass density in stars, finding that the stellar mass fraction is $\alpha=0.20\pm0.05$ at an average radius of $\sim 1.9 R_{e}$, where $R_e$ is the effective radius of the lens. This stellar mass fraction is insensitive to the X-ray source size and in excellent agreement with our earlier results based on optical data. By combining the X-ray and optical microlensing data, we can divide this larger sample into two radial bins. We find that the surface mass density in the form of stars is $\alpha=0.31\pm0.15$ and $\alpha=0.13\pm0.05$ at $(1.3\pm0.3) R_{e}$ and $(2.3\pm0.3) R_{e}$, respectively, in good agreement with expectations and some previous results.
A striking coincidence of revolution periods of S-stars orbiting a supermassive black hole at the Galactic Center of the Milky Way and oscillation periods of such solar and terrestrial observables as the sunspot number, the geomagnetic field Y-component and the global temperature is established on basis of the corresponding experimental data. Rejecting randomness of this discovered coincidence, we put forward a hypothesis that modulation of dark matter flows in the Milky Way by the S-stars is responsible for such a frequency transfer from the Galactic Center to the Solar System.
Double-barred galaxies account for almost one third of all barred galaxies, suggesting that secondary stellar bars, which are embedded in large-scale primary bars, are long-lived structures. However, up to now it has been hard to self-consistently simulate a disc galaxy that sustains two nested stellar bars for longer than a few rotation periods. N-body/hydrodynamical simulations including star formation recipes have been performed. Their properties have been compared with the most recent observational data in order to prove that they are representative of double-barred galaxies, even SB0. Overlaps in dynamical resonances and bar modes have been looked for using Fourier spectrograms. Double-barred galaxies have been successfully simulated with lifetimes as long as 7 Gyr. The stellar population of the secondary bar is younger on average than for the primary large-scale bar. An important feature of these simulations is the absence of any resonance overlap for several Gyr. In particular, there is no overlap between the primary bar ILR and the secondary bar corotation. Therefore, mode coupling cannot sustain the secondary bar mode. Star formation is identified here as possibly being responsible for bringing energy to the nuclear mode. Star formation is also responsible for limiting the amount of gas in the central region which prevents the orbits sustaining the secondary bar from being destroyed. Therefore, the secondary bar can dissolve but reappear after approx. 1 Gyr. When star formation is switched off the dynamical perturbation associated with the secondary bar needs several Gyr to fully vanish. Double-bars can be long-lived in numerical simulations with a gaseous component, even in the absence of overlap of resonances or mode coupling, provided that star formation remains active in the central region where the nuclear bar lies.
We simulate evolution of cometary H II regions based on several champagne flow models and bow shock models, and calculate the profiles of the [Ne II] fine-structure line at $12.81\mu m$, the $H30\alpha$ recombination line and the [Ne III] fine-structure line at $15.55\mu m$ for these models at different inclinations of $0^o, 30^o \textrm{and} 60^o$. We find that the profiles in the bow shock models are generally different from those in the champagne flow models, but the profiles in the bow shock with lower stellar velocity ($\leq5km s^{-1}$) are similar to those in the champagne flow models. In champagne flow models, both the velocity of peak flux and the flux weighted central velocities of all three lines are pointing outward from molecular clouds. In bow shock models, the directions of these velocities rely on the speed of stars. They have the similar motion in high stellar speed case but opposite directions in low stellar speed case. We notice that the line profiles from the slit along the symmetrical axis of the projected 2D image of these models are useful for distinguishing bow shock models and champagne flow models. It is also confirmed by the calculation that the flux weighted central velocity and the line luminosity of the [Ne III] line can be estimated from the [Ne II] line and the $H30\alpha$ line.
We present the results of a pilot JVLA project aimed at studying the bulk of the radio-emitting AGN population, unveiled by the NVSS/FIRST and SDSS surveys.We obtained A-array observations at the JVLA at 1.4, 4.5, and 7.5 GHz for 12 sources of the SDSS/NVSS sample. The radio maps reveal compact unresolved or slightly resolved radio structures on a scale of 1-3 kpc, with only one exception of a FRI/FRII source extended over $\sim$40 kpc. We isolate the radio core component in most of them. The sample splits into two groups. Four sources have small black hole (BH) masses (mostly $\sim$10$^{7}$ M$_{\odot}$) and are hosted by blue galaxies, often showing evidence of a contamination from star formation to their radio emission and associated with radio-quiet AGN. The second group consists in seven radio-loud AGN, which live in red massive ($\sim10^{11}$ M$_{\odot}$) early-type galaxies, with large BH masses ($\gtrsim$10$^{8}$ M$_{\odot}$), and spectroscopically classified as Low Excitation Galaxies, all characteristics typical of FRI radio galaxies. They also lie on the correlation between radio core power and [O III] line luminosity defined by FRIs. However, they are more core dominated (by a factor of $\sim$30) than FRIs and show a deficit of extended radio emission. We dub these sources 'FR0' to emphasize their lack of prominent extended radio emission, the single distinguishing feature with respect to FRIs. The differences in radio properties between FR0s and FRIs might be ascribed to an evolutionary effect, with the FR0 sources undergoing to rapid intermittency that prevents the growth of large scale structures. In our preferred scenario the lack of extended radio emission in FR0s is due to their smaller jet Lorentz $\Gamma$ factor with respect to FRIs, causing possible instabilities and their premature disruption.[abridged]
The spatial and temporal relationships between stellar age, kinematics, and chemistry are a fundamental tool for uncovering the physics driving galaxy formation and evolution. Observationally, these trends are derived using carefully selected samples isolated via the application of appropriate magnitude, colour, and gravity selection functions of individual stars; conversely, the analysis of chemodynamical simulations of galaxies has traditionally been restricted to the age, metallicity, and kinematics of `composite' stellar particles comprised of open cluster-mass simple stellar populations. As we enter the Gaia era, it is crucial that this approach changes, with simulations confronting data in a manner which better mimics the methodology employed by observers. Here, we use the \textsc{SynCMD} synthetic stellar populations tool to analyse the metallicity distribution function of a Milky Way-like simulated galaxy, employing an apparent magnitude plus gravity selection function similar to that employed by the RAdial Velocity Experiment (RAVE); we compare such an observationally-motivated approach with that traditionally adopted - i.e., spatial cuts alone - in order to illustrate the point that how one analyses a simulation can be, in some cases, just as important as the underlying sub-grid physics employed.
Non-local, time-dependent convection models have been used in the literature to explain the location of double-mode pulsations in Cepheids in the HR diagram as well as the existence and location of the red edge of the Cepheid instability strip. These properties are highly sensitive to model parameters. We use 2D radiation hydrodynamical simulations with realistic microphysics and grey radiative-transfer to model the upper 42 % of a short period Cepheid. The simulations show that the strength of the convection zone varies significantly over the pulsation period and exhibits a phase shift of the convective flux relative to the variations in radius. We evaluate the convective flux and the work performed by volume expansion as predicted by the most commonly used convection models. It turns out that over one pulsation cycle the model parameter $\alpha_{\rm c}$, which is proportional to the convective flux, has to be varied by up to a factor of beyond 2 to match the convective flux obtained from the simulations. To bring convective fluxes integrated over the He II convection zone and the overshoot zone below into agreement, this parameter has to be varied by a factor of $\sim 3$ (Stellingwerf) or even up to $\sim 7.5$ (Kuhfu{\ss}). We then present results on the energetics of the convection and overshoot zone by radially symmetric and fluctuating quantities. To successfully model this scenario by a static, one dimensional or even by a simple time-dependent model appears extremely challenging. We conclude that significant improvements are needed in models of Cepheids to make predictions based on them more robust and to improve the reliability of conclusions drawn on the convection-pulsation coupling from such models. Simulations as the present ones can provide guidelines for developing descriptions of convection then applied in traditional 1D modelling of radially pulsating stars.
Emission from blazar jets in the ultraviolet, optical, and infrared is polarized. If these low-energy photons were inverse-Compton scattered, the upscattered high-energy photons retain a fraction of the polarization. Current and future X-ray and gamma-ray polarimeters such as INTEGRAL-SPI, PoGOLITE, X-Calibur, Gamma-Ray Burst Polarimeter, GEMS-like missions, ASTRO-H, and POLARIX have the potential to discover polarized X-rays and gamma-rays from blazar jets for the first time. Detection of such polarization will open a qualitatively new window into high-energy blazar emission; actual measurements of polarization degree and angle will quantitatively test theories of jet emission mechanisms. We examine the detection prospects of blazars by these polarimetry missions using examples of 3C 279, PKS 1510-089, and 3C 454.3, bright sources with relatively high degrees of low-energy polarization. We conclude that while balloon polarimeters will be challenged to detect blazars within reasonable observational times (with X-Calibur offering the most promising prospects), space-based missions should detect the brightest blazars for polarization fractions down to a few percent. Typical flaring activity of blazars could boost the overall number of polarimetric detections by nearly a factor of five to six purely accounting for flux increase of the brightest of the comprehensive, all-sky, Fermi-LAT blazar distribution. The instantaneous increase in the number of detections is approximately a factor of two, assuming a duty cycle of 20% for every source. The detectability of particular blazars may be reduced if variations in the flux and polarization fraction are anticorrelated. Simultaneous use of variability and polarization trends could guide the selection of blazars for high-energy polarimetric observations.
(Abridged) Gaia aims to make a 3-dimensional map of 1,000 million stars in our Milky Way to unravel its kinematical, dynamical, and chemical structure and evolution. Gaia's on-board detection software discriminates stars from spurious objects like cosmic rays and Solar protons. For this, parametrised point-spread-function-shape criteria are used. This study aims to provide an optimum set of parameters for these filters. We developed an emulation of the on-board detection software, which has 20 free, so-called rejection parameters which govern the boundaries between stars on the one hand and sharp or extended events on the other hand. We evaluate the detection and rejection performance of the algorithm using catalogues of simulated single stars, double stars, cosmic rays, Solar protons, unresolved galaxies, and asteroids. We optimised the rejection parameters, improving - with respect to the functional baseline - the detection performance of single and double stars, while, at the same time, improving the rejection performance of cosmic rays and of Solar protons. We find that the minimum separation to resolve a close, equal-brightness double star is 0.23 arcsec in the along-scan and 0.70 arcsec in the across-scan direction, independent of the brightness of the primary. We find that, whereas the optimised rejection parameters have no significant impact on the detectability of de Vaucouleurs profiles, they do significantly improve the detection of exponential-disk profiles. We also find that the optimised rejection parameters provide detection gains for asteroids fainter than 20 mag and for fast-moving near-Earth objects fainter than 18 mag, albeit this gain comes at the expense of a modest detection-probability loss for bright, fast-moving near-Earth objects. The major side effect of the optimised parameters is that spurious ghosts in the wings of bright stars essentially pass unfiltered.
A black hole accretion is necessarily transonic. In presence of sufficiently high viscosity and cooling effects, a low-angular momentum transonic flow can become a standard Keplerian disc except close to the where hole where it must pass through the inner sonic point. However, if the viscosity is not high everywhere and cooling is not efficient everywhere, the flow cannot completely become a Keplerian disc. In this paper, we show results of rigorous numerical simulations of a transonic flow having vertically varying viscosity parameter (being highest on the equatorial plane) and optical depth dependent cooling processes to show that the flow indeed segregates into two distinct components as it approaches a black hole. The component on the equatorial plane has properties of a standard Keplerian disc, though the flow is not truncated at the inner- most stable circular orbit. This component extends till the horizon as a sub-Keplerian flow. This standard disc is found to be surrounded by a hot, low angular momentum component forming a centrifugal barrier dominated oscillating shock wave, consistent with the Chakrabarti-Titarchuk two component advective flow configuration.
We have carried out high resolution observations with Atacama Large Millimeter/Submillimeter Array (ALMA) of continuum emission from Orion KL region. We identify 11 compact sources at ALMA band 6 (245 GHz) and band 7 (339 GHz), including Hot Core, Compact Ridge, SMA1, IRc4, IRc7, and a radio source I (Source I). Spectral energy distribution (SED) of each source is determined by using previous 3 mm continuum emission data. Physical properties such as size, mass, hydrogen number density and column density are discussed based on the dust graybody SED. Among 11 identified sources, Source I, a massive protostar candidate, is a dominant energy source in Orion KL. We extensively investigate its SED from centimeter to submillimeter wavelengths. The SED of Source I can be fitted with a single power-law index of 1.97 suggesting an optically thick emission. We employ the H$^{-}$ free-free emission as an opacity source of this optically thick emission. The temperature, density, and mass of the circumstellar disk associated with Source I are constrained by the SED of H$^{-}$ free-free emission. Still the fitting result shows a significant deviation from the observed flux densities. Combined with the thermal dust graybody SED to explain excess emission at higher frequency, a smaller power-law index of 1.60 for the H$^{-}$ free-free emission is obtained in the SED fitting. The power-law index smaller than 2 would suggest a compact source size or a clumpy structure unresolved with the present study. Future higher resolution observations with ALMA are essential to reveal more detailed spatial structure and physical properties of Source I.
Asteroid (234) Barbara is the prototype of a category of asteroids that has been shown to be extremely rich in refractory inclusions, the oldest material ever found in the Solar System. It exhibits several peculiar features, most notably its polarimetric behavior. In recent years other objects sharing the same property (collectively known as "Barbarians") have been discovered. Interferometric observations in the mid-infrared with the ESO VLTI suggested that (234) Barbara might have a bi-lobated shape or even a large companion satellite. We use a large set of 57 optical lightcurves acquired between 1979 and 2014, together with the timings of two stellar occultations in 2009, to determine the rotation period, spin-vector coordinates, and 3-D shape of (234) Barbara, using two different shape reconstruction algorithms. By using the lightcurves combined to the results obtained from stellar occultations, we are able to show that the shape of (234) Barbara exhibits large concave areas. Possible links of the shape to the polarimetric properties and the object evolution are discussed. We also show that VLTI data can be modeled without the presence of a satellite.
In this paper a description is given of the SFXC software correlator, developed and maintained at the Joint Institute for VLBI in Europe (JIVE). The software is designed to run on generic Linux-based computing clusters. The correlation algorithm is explained in detail, as are some of the novel modes that software correlation has enabled, such as wide-field VLBI imaging through the use of multiple phase centres and pulsar gating and binning. This is followed by an overview of the software architecture. Finally, the performance of the correlator as a function of number of CPU cores, telescopes and spectral channels is shown.
We present a model for producing tidal streams from disrupting progenitors in arbitrary potentials, utilizing the idea that the majority of stars escape from the progenitor's two Lagrange points. The method involves releasing test particles at the Lagrange points as the satellite orbits the host and dynamically evolving them in the potential of both host and progenitor. The method is sufficiently fast to allow large-dimensional parameter exploration using Monte Carlo methods. We provide the first direct modelling of 6-D stream observations -- assuming a stream rather than an orbit -- by applying our methods to GD-1. This is a kinematically cold stream spanning $60^{\circ}$ of the sky and residing in the outer Galaxy $\approx 15$ kpc distant from the centre. We assume the stream moves in a flattened logarithmic potential characterised by an asymptotic circular velocity $v_0$ and a flattening $q$. We recover values of normalisation $v_0$ = $227.2^{+15.6}_{-18.2}$ kms$^{-1}$ and flattening $q$ = $0.91^{+0.04}_{-0.1}$, if the stream is assumed to leading, and $v_0$ = $226.5^{+17.9}_{-17.0}$ kms$^{-1}$, $q$ = $0.90^{+0.05}_{-0.09}$, if it is assumed to be trailing. This can be compared to the values $v_0 = 224 \pm 13$ kms$^{-1}$ and $q= 0.87^{+0.07}_{-0.04}$ obtained by Koposov et al (2010) using the simpler technique of orbit fitting. Although there are differences between stream and orbit fitting, we conclude that orbit fitting can provide accurate results given the current quality of the data, at least for this kinematically cold stream in this logarithmic model of the Galaxy.
During the formation of a star, material is ejected along powerful jets that impact the ambient material. This outflow regulates star formation by e.g. inducing turbulence and heating the surrounding gas. Understanding the associated shocks is therefore essential to the study of star formation. We present comparisons of shock models with CO, H2, and SiO observations in a 'pure' shock position in the BHR71 bipolar outflow. These comparisons provide an insight into the shock and pre-shock characteristics, and allow us to understand the energetic and chemical feedback of star formation on Galactic scales. New CO (Jup = 16, 11, 7, 6, 4, 3) observations from the shocked regions with the SOFIA and APEX telescopes are presented and combined with earlier H2 and SiO data (from the Spitzer and APEX telescopes). The integrated intensities are compared to a grid of models that were obtained from a magneto-hydrodynamical shock code which calculates the dynamical and chemical structure of these regions combined with a radiative transfer module based on the 'large velocity gradient' approximation. The CO emission leads us to update the conclusions of our previous shock analysis: pre-shock densities of 1e4 cm-3 and shock velocities around 20-25 km s-1 are still constrained, but older ages are inferred ( 4000 years). We evaluate the contribution of shocks to the excitation of CO around forming stars. The SiO observations are compatible with a scenario where less than 4% of the pre-shock SiO belongs to the grain mantles. We infer outflow parameters: a mass of 1.8x1e-2 Msun was measured in our beam, in which a momentum of 0.4 Msun km s-1 is dissipated, for an energy of 4.2x1e43erg. We analyse the energetics of the outflow species by species. Comparing our results with previous studies highlights their dependence on the method: H2 observations only are not sufficient to evaluate the mass of outflows.
Major components of ices on interstellar grains in molecular clouds - water and carbon oxides - occur at various optical depths. This implies that selective desorption mechanisms are at work. An astrochemical model of a contracting low-mass molecular cloud core is presented. Ice was treated as consisting of the surface and three subsurface layers (sublayers). Photodesorption, reactive desorption, and indirect reactive desorption were investigated. The latter manifests itself through desorption from H+H reaction on grains. Desorption of shallow subsurface species was included. Modeling results suggest the existence of a "photon-dominated ice" during the early phases of core contraction. Subsurface ice is chemically processed by interstellar photons, which produces complex organic molecules. Desorption from the subsurface layer results in high COM gas-phase abundances at Av = 2.4...10mag. This may contribute towards an explanation for COM observations in dark cores. It was found that photodesorption mostly governs the onset of ice accumulation onto grains. Reaction-specific reactive desorption is efficient for small molecules that form via highly exothermic atom-addition reactions. Higher reactive desorption efficiency results in lower gas-phase abundances of COMs. Indirect reactive desorption allows to closely reproduce the observed H2O:CO:CO2 ratio towards a number of background stars. Presumably this can be done by any mechanism whose efficiency fits with the sequence CO > CO2 >> H2O. After the freeze-out has ended, the three sublayers represent chemically distinct parts of the mantle. 8...10.5mag is the likely AV threshold for the appearance of CO ice. The lower value is supported by observations.
New integral field spectroscopy has been obtained for IZw18, the nearby lowest-metallicity galaxy considered our best local analog of systems forming at high-z. Here we report the spatially resolved spectral map of the nebular HeII4686 emission in IZw18, from which we derived for the first time its total HeII-ionizing flux. Nebular HeII emission implies the existence of a hard radiation field. HeII-emitters are observed to be more frequent among high-z galaxies than for local objects. So investigating the HeII-ionizing source(s) in IZw18 may reveal the ionization processes at high-z. HeII emission in star-forming galaxies, has been suggested to be mainly associated with Wolf-Rayet stars (WRs), but WRs cannot satisfactorily explain the HeII-ionization at all times, in particular at lowest metallicities. Shocks from supernova remnants, or X-ray binaries, have been proposed as additional potential sources of HeII-ionizing photons. Our data indicate that conventional HeII-ionizing sources (WRs, shocks, X-ray binaries) are not sufficient to explain the observed nebular HeII4686 emission in IZw18. We find that the HeII-ionizing radiation expected from models for either low-metallicity super-massive O stars or rotating metal-free stars could account for the HeII-ionization budget measured, while only the latter models could explain the highest values of HeII4686/Hbeta observed. The presence of such peculiar stars in IZw18 is suggestive and further investigation in this regard is needed. This letter highlights that some of the clues of the early Universe can be found here in our cosmic backyard.
We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.
We investigate the orbital dynamics of a \textit{barred-spiral} model when the system is rotating slowly and corotation is located beyond the end of the spiral arms. In the characteristic of the central family of periodic orbits we find a "bistable region". In the response model we observe a ring surrounding the bar and spiral arms starting tangential to the ring. This is a morphology resembling barred-spiral systems with inner rings. However, the dynamics associated with this structure in the case we study is different from that of a typical bar ending close to corotation. The ring of our model is round, or rather elongated perpendicular to the bar. It is associated with a folding (an "S" shaped feature) of the characteristic of the central family, which is typical in bistable bifurcations. Along the "S" part of the characteristic we have a change in the orientation of the periodic orbits from a x1-type to a x2-type morphology. The orbits populated in the response model change rather abruptly their orientation when reaching the lowest energy of the "S". The spirals of the model follow a standard "precessing ellipses flow" and the orbits building them have energies beyond the "S" region. The bar is structured mainly by sticky orbits from regions around the stability islands of the central family. This leads to the appearance of X-features in the bars \textit{on} the galactic plane. Such a bar morphology appears in the unsharp-masked images of some moderately inclined galaxies.
Type Ia supernovae have been proposed to be much better distance indicators at near-infrared compared to optical wavelengths -- the effect of dust extinction is expected to be lower and it has been shown that SNe Ia behave more like `standard candles' at NIR wavelengths. To better understand the physical processes behind this increased uniformity, we have studied the $Y$, $J$ and $H$-filter light curves of 91 SNe Ia from the literature. We show that the phases and luminosities of the first maximum in the NIR light curves are extremely uniform for our sample. The phase of the second maximum, the late-phase NIR luminosity and the optical light curve shape are found to be strongly correlated, in particular more luminous SNe Ia reach the second maximum in the NIR filters at a later phase compared to fainter objects. We also find a strong correlation between the phase of the second maximum and the epoch at which the SN enters the Lira law phase in its optical colour curve (epochs $\sim$ 15 to 30 days after $B$ band maximum). The decline rate after the second maximum is very uniform in all NIR filters. We suggest that these observational parameters are linked to the nickel and iron mass in the explosion, providing evidence that the amount of nickel synthesised in the explosion is the dominating factor shaping the optical and NIR appearance of SNe Ia.
Giant gas planets in close proximity to their host stars experience strong irradiation. In extreme cases photoevaporation causes a transonic, planetary wind and the persistent mass loss can possibly affect the planetary evolution. We have identified nine hot Jupiter systems in the vicinity of the Sun, in which expanded planetary atmospheres should be detectable through Lyman alpha transit spectroscopy according to predictions. We use X-ray observations with Chandra and XMM-Newton of seven of these targets to derive the high-energy irradiation level of the planetary atmospheres and the resulting mass loss rates. We further derive improved Lyman alpha luminosity estimates for the host stars including interstellar absorption. According to our estimates WASP-80 b, WASP-77 b, and WASP-43 b experience the strongest mass loss rates, exceeding the mass loss rate of HD 209458 b, where an expanded atmosphere has been confirmed. Furthermore, seven out of nine targets might be amenable to Lyman alpha transit spectroscopy. Finally, we check the possibility of angular momentum transfer from the hot Jupiters to the host stars in the three binary systems among our sample, but find only weak indications for increased stellar rotation periods of WASP-77 and HAT-P-20.
A new approach to cosmological perturbation theory has been recently introduced by Bartelmann et al., relying on non-equilibrium statistical theory of classical particles, and treating the gravitational interaction perturbatively. They compute analytic expressions for the non-linear matter power spectrum, to first order in the interaction, and at one-loop order in the linear power spectrum. The resulting power spectrum is well-behaved even at large wavenumbers and seems in good agreement with results from numerical simulations. In this paper, we rederive their results concisely with a different approach, starting from the implicit integral solution to particle trajectories. We derive the matter power spectrum to first order in the interaction, but to arbitrary order in the linear power spectrum, from which the one-loop result follows. We also show that standard linear perturbation theory can only be recovered at infinite order in the gravitational interaction. At finite order in the interaction, we find that the linear power spectrum is systematically and significantly underestimated. A comprehensive study of the convergence of the theory with the order of the interaction for non-linear scales will be the subject of future work.
Non-Gaussianity of the primordial density perturbations provides an important measure to constrain models of inflation. At cubic order the non-Gaussianity is captured by two parameters $\tau_{\rm NL}$ and $g_{\rm NL}$ that determine the amplitude of the density perturbation trispectrum. Here we report measurements of the kurtosis power spectra of the cosmic microwave background (CMB) temperature as mapped by Planck by making use of correlations between square temperature-square temperature and cubic temperature-temperature anisotropies. In combination with noise simulations, we find the best joint estimates to be $\tau_{\rm{NL}}=0.3 \pm 0.9 \times 10^4$ and $g_{\rm{NL}}=-1.2 \pm 2.8 \times 10^5$. If $\tau_{\rm NL}=0$, we find $g_{\rm NL}= -1.3\pm 1.8 \times 10^5$.
Left-over, ablated material from a possible non-degenerate companion can reveal itself after about one year in spectra of Type Ia SNe (SNe Ia). We have searched for such material in spectra of SN 2011fe (at 294 days after the explosion) and for SN 2014J (315 days past explosion). The observations are compared with numerical models simulating the expected line emission. The spectral lines sought for are H-alpha, [O I] 6300 and [Ca II] 7291,7324, and the expected width of these lines is about 1000 km/s. No signs of these lines can be traced in any of the two supernovae. When systematic uncertainties are included, the limits on hydrogen-rich ablated gas in SNe 2011fe and 2014J are 0.003 M_sun and 0.0085 M_sun, respectively, where the limit for SN 2014J is the second lowest ever, and the limit for SN 2011fe is a revision of a previous limit. Limits are also put on helium-rich ablated gas. These limits are used, in conjunction with other data, to argue that these supernovae can stem from double-degenerate systems, or from single-degenerate systems with a spun up/spun down super-Chandrasekhar white dwarf. For SN 2011fe, other types of hydrogen-rich donors can likely be ruled out, whereas for SN 2014J a main-sequence donor system with large intrinsic separation is still possible. Helium-rich donor systems cannot be ruled out for any of the two supernovae, but the expected short delay time for such progenitors makes this possibility less likely, especially for SN 2011fe. The broad [Ni II] 7378 emission in SN 2014J is redshifted by about +1300 km/s, as opposed to the known blueshift of roughly -1100 km/s for SN 2011fe. [Fe II] 7155 is also redshifted in SN 2014J. SN 2014J belongs to a minority of SNe Ia that both have a nebular redshift of [Fe II] 7155 and [Ni II] 7378, and a slow decline of the Si II 6355 absorption trough just after B-band maximum.
BICEP2 and the Keck Array are polarization-sensitive microwave telescopes that observe the cosmic microwave background (CMB) from the South Pole at degree angular scales in search of a signature of inflation imprinted as B-mode polarization in the CMB. BICEP2 was deployed in late 2009, observed for three years until the end of 2012 at 150 GHz with 512 antenna-coupled transition edge sensor bolometers, and has reported a detection of B-mode polarization on degree angular scales. The Keck Array was first deployed in late 2010 and will observe through 2016 with five receivers at several frequencies (95, 150, and 220 GHz). BICEP2 and the Keck Array share a common optical design and employ the field-proven BICEP1 strategy of using small-aperture, cold, on-axis refractive optics, providing excellent control of systematics while maintaining a large field of view. This design allows for full characterization of far-field optical performance using microwave sources on the ground. Here we describe the optical design of both instruments and report a full characterization of the optical performance and beams of BICEP2 and the Keck Array at 150 GHz.
In a companion paper we have reported a $>5\sigma$ detection of degree scale $B $-mode polarization at 150 GHz by the BICEP2 experiment. Here we provide a detailed study of potential instrumental systematic contamination to that measurement. We focus extensively on spurious polarization that can potentially arise from beam imperfections. We present a heuristic classification of beam imperfections according to their symmetries and uniformities, and discuss how resulting contamination adds or cancels in maps that combine observations made at multiple orientations of the telescope about its boresight axis. We introduce a technique, which we call "deprojection", for filtering the leading order beam-induced contamination from time ordered data, and show that it removes power from BICEP2's $BB$ spectrum consistent with predictions using high signal-to-noise beam shape measurements. We detail the simulation pipeline that we use to directly simulate instrumental systematics and the calibration data used as input to that pipeline. Finally, we present the constraints on $BB$ contamination from individual sources of potential systematics. We find that systematics contribute $BB$ power that is a factor $\sim10\times$ below BICEP2's 3-year statistical uncertainty, and negligible compared to the observed $BB$ signal. The contribution to the best-fit tensor/scalar ratio is at a level equivalent to $r=(3-6)\times10^{-3}$.
We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5\deg$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $\mu$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies above 150 GHz to a lensed-$\Lambda$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$). We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,\sigma$ significance.
We have developed antenna-coupled transition-edge sensor (TES) bolometers for a wide range of cosmic microwave background (CMB) polarimetry experiments, including BICEP2, Keck Array, and the balloon borne SPIDER. These detectors have reached maturity and this paper reports on their design principles, overall performance, and key challenges associated with design and production. Our detector arrays repeatedly produce spectral bands with 20%-30% bandwidth at 95, 150, or 220~GHz. The integrated antenna arrays synthesize symmetric co-aligned beams with controlled side-lobe levels. Cross-polarized response on boresight is typically ~0.5%, consistent with cross-talk in our multiplexed readout system. End-to-end optical efficiencies in our cameras are routinely 35% or higher, with per detector sensitivities of NET~300 uKrts. Thanks to the scalability of this design, we have deployed 2560 detectors as 1280 matched pairs in Keck Array with a combined instantaneous sensitivity of ~9 uKrts, as measured directly from CMB maps in the 2013 season. Similar arrays have recently flown in the SPIDER instrument, and development of this technology is ongoing.
In the Lema\^{\i}tre -- Tolman (L--T) models that have nonconstant bang-time function $t_B(r)$, light emitted close to those points of the Big Bang where $\dril {t_B} r \neq 0$ is blueshifted at the receiver position. The blueshifted rays are expected to perturb the temperature of the cosmic microwave background (CMB) radiation along the lines of sight of the present central observer. It is shown here that, in a general L--T model, the $t_B(r)$ can be chosen so that the blueshift-generating region is hidden before the recombination time, where the L--T model does not apply. The rest of the paper is devoted to investigating blueshifts in one specific L--T model, called L--T$(t_B)$ -- the one that duplicates the luminosity distance vs. redshift relation of the $\Lambda$CDM model using nonconstant $t_B(r)$ alone. The location of the blueshift-generating region in the L--T$(t_B)$ spacetime is determined. Profiles of redshift/blueshift along several rays intersecting the past light cone of the present central observer are calculated. The L--T$(t_B)$ model matched to Friedmann is considered, and profiles of redshift/blueshift in such a composite model are calculated. The requirement of invisibility of blueshifts makes the L--T$(t_B)$ model astrophysically unacceptable if it should apply back to the recombination time, but does not "rule out" a general L--T model -- it only puts limits on $\dril {t_B} r$.
The parallelization, design and scalability of the \sky code to search for periodic gravitational waves from rotating neutron stars is discussed. The code is based on an efficient implementation of the F-statistic using the Fast Fourier Transform algorithm. To perform an analysis of data from the advanced LIGO and Virgo gravitational wave detectors' network, which will start operating in 2015, hundreds of millions of CPU hours will be required - the code utilizing the potential of massively parallel supercomputers is therefore mandatory. We have parallelized the code using the Message Passing Interface standard, implemented a mechanism for combining the searches at different sky-positions and frequency bands into one extremely scalable program. The parallel I/O interface is used to escape bottlenecks, when writing the generated data into file system. This allowed to develop a highly scalable computation code, which would enable the data analysis at large scales on acceptable time scales. Benchmarking of the code on a Cray XE6 system was performed to show efficiency of our parallelization concept and to demonstrate scaling up to 50 thousand cores in parallel.
The stability of magnetized strange quark matter (MSQM) is studied in the MIT bag model with the density dependent bag pressure. In the consistent thermodynamic description of MSQM, the quark chemical potentials, the total thermodynamic potential and the anisotropic pressure acquire the corresponding additional term proportional to the density derivative of the bag pressure. The model parameter space is determined, for which MSQM is absolutely stable, i.e., its energy per baryon is less than that of the most stable $^{56}$Fe nucleus under the zero external pressure and vanishing temperature. It is shown that there exists the magnetic field strength $H_{u\,max}$ at which the upper bound $B_\infty^u$ on the asymptotic bag pressure $B_\infty\equiv B(\varrho_B\gg \varrho_0$) ($\varrho_0$ being the nuclear saturation density) from the absolute stability window vanishes. The value of this field, \hbox{$H_{u\,max}\sim$$(1$--$3)\cdot10^{18}$}~G, represents the upper bound on the magnetic field strength, which can be reached in a strongly magnetized strange quark star. It is clarified how the absolute stability window and upper bound on the magnetic field strength are affected by varying the parameters in the Gaussian parametrization for the density dependence of the bag pressure.
The moon-forming impact and the subsequent evolution of the proto-Earth is strongly dependent on the properties of materials at the extreme conditions generated by this violent collision. We examine the high pressure behavior of MgO, one of the dominant constituents in the earth's mantle, using high-precision, plate impact shock compression experiments performed on Sandia National Laboratories Z-Machine and extensive quantum simulations using Density Functional Theory (DFT) and quantum Monte Carlo (QMC). The combined data span from ambient conditions to 1.2 TPa and 42,000 K, showing solid-solid and solid-liquid phase boundaries. Furthermore our results indicate under impact that the solid and liquid phases coexist for more than 100 GPa, pushing complete melting to pressures in excess of 600 GPa. The high pressure required for complete shock melting places a lower bound on the relative velocities required for the moon forming impact.
We introduce the notion of a local shadow for a black hole and determine its shape for the particular case of a distorted Schwarzschild black hole. Considering the lowest-order even and odd multiple moments, we compute the relation between the deformations of the shadow of a Schwarzschild black hole and the distortion multiple moments. For the range of values of multiple moments that we consider, the horizon is deformed much less than its corresponding shadow, suggesting the horizon is more `rigid'. Quite unexpectedly we find that a prolate distortion of the horizon gives rise to an oblate distortion of the shadow, and vice-versa.
Within the framework of shallow-water magnetohydrodynamics, we investigate the linear instability of horizontal shear flows, influenced by an aligned magnetic field and stratification. Various classical instability results, such as H{\o}iland's growth rate bound and Howard's semi-circle theorem, are extended to this shallow-water system for quite general profiles. Two specific piecewise-constant velocity profiles, the vortex sheet and the rectangular jet, are studied analytically and asymptotically; it is found that the magnetic field and stratification (as measured by the Froude number) are generally both stabilising, but weak instabilities can be found at arbitrarily large Froude number. Numerical solutions are computed for corresponding smooth velocity profiles, the hyperbolic-tangent shear layer and the Bickley jet, for a uniform background field. A generalisation of the long-wave asymptotic analysis of Drazin & Howard (1962) is employed in order to understand the instability characteristics for both profiles. For the shear layer, the mechanism underlying the primary instability is interpreted in terms of counter-propagating Rossby waves, thereby allowing an explication of the stabilising effects of the magnetic field and stratification.
The direct detection of gravitational waves with the next generation detectors, like Advanced LIGO, provides the opportunity to measure deviations from the predictions of General Relativity. One such departure would be the existence of alternative polarizations. To measure these, we study a single detector measurement of a continuous gravitational wave from a triaxial pulsar source. We develop methods to detect signals of any polarization content and distinguish between them in a model independent way. We present LIGO S5 sensitivity estimates for 115 pulsars.
X-ray spectrometer is one of the satellite payloads on Chang'E-2 satellite. The soft X-ray detector is one of the device on X-ray spectrometer which is designed to detect the major rock-forming elements within 0.5-10keV range on lunar surface. In this paper, energy linearity and energy resolution calibration is done using a weak Fe55 source, while temperature and time effect is considered not take big error. The total uncertainty is estimated to be within 5% after correction.
We study experimental and cosmological constraints on the extension of the Standard Model by three right handed neutrinos with masses between those of the pion and W-boson. This low scale seesaw scenario allows to simultaneously explain the observed neutrino oscillations and baryon asymmetry of the universe. We combine indirect experimental constraints from neutrinoless double $\beta$-decay, lepton flavour violation and neutrino oscillation data with bounds from past direct searches and big bang nucleosynthesis. For masses of a few GeV the heavy right handed neutrinos can be found in meson decays at LHCb, BELLE II or the proposed SHiP experiment, for larger masses they can be searched for in ATLAS and CMS. The chances for a discovery would be even better at a Future Circular Collider.
A hypothesis is proposed that the gamma-ray bursts (GRBs) may arise by blueshifting the emission radiation of hydrogen and helium generated during the last scattering epoch. The blueshift mechanism is provided by such a Lema\^{\i}tre -- Tolman (L--T) model, in which the bang-time function $t_B(r)$ is not everywhere constant. Blueshift arises on \textit{radial} rays that are emitted over regions where $\dril{t_B} r \neq 0$. The paper presents an L--T model adapted for this purpose and shows how it accounts for the observed properties of the GRBs; some properties are accounted for only qualitatively.
The evolution of the filamentation instability produced by two counter-streaming pair plasmas is studied with particle-in-cell (PIC) simulations in both one (1D) and two (2D) spatial dimensions. Radiation friction effects on particles are taken into account. After an exponential growth of both the magnetic field and the current density, a nonlinear quasi-stationary phase sets up characterized by filaments of opposite currents. During the nonlinear stage, a strong broadening of the particle energy spectrum occurs accompanied by the formation of a peak at twice their initial energy. A simple theory of the peak formation is presented. The presence of radiative losses does not change the dynamics of the instability but affects the structure of the particle spectra.
We report simulations of the inspiral and merger of binary neutron stars performed with \texttt{WhiskyTHC}, the first of a new generation of numerical relativity codes employing higher than second-order methods for both the spacetime and the hydrodynamic evolution. We find that the use of higher-order schemes improves substantially the quality of the gravitational waveforms extracted from the simulations when compared to those computed using traditional second-order schemes. The reduced de-phasing and the faster convergence rate allow us to estimate the phase evolution of the gravitational waves emitted, as well as the magnitude of finite-resolution effects, without the need of phase- or time-alignments or rescalings of the waves, as sometimes done in other works. Furthermore, by using an additional unpublished simulation at very high resolution, we confirm the robustness of our high convergence order of $3.2$.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
We report on high resolution CO(1-0), CS(2-1) and 3mm continuum Combined Array for Research in Millimeter Astronomy (CARMA) observations of the molecular outflow host and nearest quasar Markarian 231. We use the CS(2-1) measurements to derive a dense gas mass within Mrk 231 of $1.8\pm0.3\times10^{10}$ $M_\odot$, quite consistent with previous measurements. The CS(2-1) data also seem to indicate that the molecular disk of Mrk 231 is forming stars at normal efficiency. The high resolution CARMA observations were able to resolve the CO(1-0) outflow into two distinct lobes, allowing for a size estimate to be made and further constraining the molecular outflow dynamical time, further constraining the molecular gas escape rate. We find that 15% of the molecular gas within the Mrk 231 outflow actually exceeds the escape velocity in the central kiloparsec. Assuming that molecular gas is not constantly being accelerated, we find the depletion timescale of molecular gas in Mrk 231 to be 49 Myr, rather than 32 Myr, more consistent with the poststarburst stellar population observed in the system.
The Planck satellite has identified several patches of sky with low polarized dust emission, obvious targets for searches for the cosmic-microwave-background (CMB) B-mode signal from inflationary gravitational waves. Still, given the Planck measurement uncertainties, the polarized dust foregrounds in these different candidate patches may differ by an order of magnitude or more. Here we show that a brief initial experiment to map these candidate patches more deeply at a single high frequency can efficiently zero in on the cleanest patch(es) and thus improve significantly the sensitivity of subsequent B-mode searches. A ground-based experiment with current detector technology operating at >~220 GHz for 3 months can efficiently identify a low-dust-amplitude patch and thus improve by up to 20%-60% on the sensitivity to cosmic B modes of the subsequent lower-frequency deep integration. A balloon experiment with current detector sensitivities covering the set of patches and operating at ~350 GHz can reach a similar result in less than two weeks. This strategy may prove crucial in accessing the smallest gravitational-wave signals possible in large-field inflation. The high-frequency data from this exploratory experiment should also provide valuable foreground templates to subsequent experiments that integrate on any of the candidate patches explored.
We present accurate models of the gravitational potential produced by a radially exponential disk mass distribution. The models are produced by combining three separate Miyamoto-Nagai disks. Such models have been used previously to model the disk of the Milky Way, but here we extend this framework to allow its application to disks of any mass, scalelength, and a wide range of thickness from infinitely thin to near spherical (ellipticities from 0 to 0.9). The models have the advantage of simplicity of implementation, and we expect faster run speeds over a double exponential disk treatment. The potentials are fully analytical, and differentiable at all points. The mass distribution of our models deviates from the radial mass distribution of a pure exponential disk by <0.4% out to 4 disk scalelengths, and <1.9% out to 10 disk scalelengths. We tabulate fitting parameters which facilitate construction of exponential disks for any scalelength, and a wide range of disk thickness (a user-friendly, web-based interface is also available). Our recipe is well suited for numerical modelling of the tidal effects of a giant disk galaxy on star clusters or dwarf galaxies. We consider three worked examples; the Milky Way thin and thick disk, and a disky dwarf galaxy.
Large-scale peculiar motion modulates the observed luminosity distribution of galaxies. Using about half a million SDSS galaxies, this can be harnessed to obtain bounds on peculiar velocity moments, the amplitude of the linear matter power spectrum, $\sigma_{8}$, and the growth rate of density perturbations at z ~ 0.1. Results obtained from this approach agree well with the predictions of the $\Lambda$CDM model and are consistent with the reported ~ 1% zero-point tilt in the SDSS photometry.
After the discovery of the first connection between GRBs and SNe almost two decades ago, tens of SN-like rebrightenings have been discovered and about seven solid associations have been spectroscopically confirmed to date. Using GROND optical/NIR data and Swift X-ray/UV data, we estimate the intrinsic extinction, luminosity, and evolution of three SN rebrightenings in GRB afterglow light curves at z~0.5. The SNe 2008hw, 2009nz, and 2010ma exhibit 0.80, 1.15, and 1.78 times the optical (r band) luminosity of SN 1998bw, respectively. While SN 2009nz evolves similarly to SN 1998bw, SNe 2008hw and 2010ma show earlier peak times. The quasi-bolometric light curves were corrected for the contribution of the NIR bands using data available in the literature and blackbody fits. The large luminosity of SN 2010ma (1.4x10^43 erg/s) is confirmed, while SNe 2008hw and 2009nz reached a peak luminosity closer to SN 1998bw. Physical parameters of the SN explosions, such as synthesised nickel mass, ejecta mass, and kinetic energy, are estimated using Arnett's analytic approach, which resulted in nickel masses of around 0.4-0.5 Msun. By means of the a very comprehensive data set, we found that the luminosity and the nickel mass of SNe 2008hw, 2009nz, and 2010ma resembles those of other known GRB-associated SNe. This findings strengthens previous claims of GRB-SNe being brighter than type-Ic SNe unaccompanied by GRBs.
We study the dependence of protoplanetary disk evolution on stellar mass using a large sample of young stellar objects in nearby young star-forming regions. We update the protoplanetary disk fractions presented in our recent work (paper I of this series) derived for 22 nearby (< 500 pc) associations between 1 and 100 Myr. We use a subsample of 1 428 spectroscopically confirmed members to study the impact of stellar mass on protoplanetary disk evolution. We divide this sample into two stellar mass bins (2 M$_{\odot}$ boundary) and two age bins (3 Myr boundary), and use infrared excesses over the photospheric emission to classify objects in three groups: protoplanetary disks, evolved disks, and diskless. The homogeneous analysis and bias corrections allow for a statistically significant inter-comparison of the obtained results. We find robust statistical evidence of disk evolution dependence with stellar mass. Our results, combined with previous studies on disk evolution, confirm that protoplanetary disks evolve faster and/or earlier around high-mass (> 2 M$_{\odot}$) stars. We also find a roughly constant level of evolved disks throughout the whole age and stellar mass spectra. We conclude that protoplanetary disk evolution depends on stellar mass. Such a dependence could have important implications for gas giant planet formation and migration, and could contribute to explaining the apparent paucity of hot Jupiters around high-mass stars.
We have conducted an optical long-slit spectroscopic survey of 1022 galaxies using the 10m Hobby-Eberly Telescope (HET) at McDonald Observatory. The main goal of the HET Massive Galaxy Survey (HETMGS) is to find nearby galaxies that are suitable for black hole mass measurements. In order to measure accurately the black hole mass, one should kinematically resolve the region where the black hole dominates the gravitational potential. For most galaxies, this region is much less than an arcsecond. Thus, black hole masses are best measured in nearby galaxies with telescopes that obtain high-spatial resolution. The HETMGS focuses on those galaxies predicted to have the largest sphere-of-influence, based on published stellar velocity dispersions or the galaxy fundamental plane. To ensure coverage over galaxy types, the survey targets those galaxies across a face-on projection of the fundamental plane. We present the sample selection and resulting data products from the long-slit observations, including central stellar kinematics and emission line ratios. The full dataset, including spectra and resolved kinematics, is available online. Additionally, we show that the current crop of black hole masses are highly biased towards dense galaxies and that especially large disks and low dispersion galaxies are under-represented. This survey provides the necessary groundwork for future systematic black hole mass measurement campaigns.
We construct dynamical models of the Milky Way's Box/Peanut (B/P) bulge, using the recently measured 3D density of Red Clump Giants (RCGs) as well as kinematic data from the BRAVA survey. We match these data using the NMAGIC Made-to-Measure method, starting with N-body models for barred discs in different dark matter haloes. We determine the total mass in the bulge volume of the RCGs measurement (+-2.2 x +- 1.4 x +- 1.2 kpc) with unprecedented accuracy and robustness to be 1.84 +- 0.07 x10^10 Msun. The stellar mass in this volume varies between 1.25-1.6 x10^10 Msun, depending on the amount of dark matter in the bulge. We evaluate the mass-to-light and mass-to-clump ratios in the bulge and compare them to theoretical predictions from population synthesis models. We find a mass-to-light ratio in the K-band in the range 0.8-1.1. The models are consistent with a Kroupa or Chabrier IMF, but a Salpeter IMF is ruled out for stellar ages of 10 Gyr. To match predictions from the Zoccali IMF derived from the bulge stellar luminosity function requires about 40% or 0.7 x10^10 Msun dark matter in the bulge region. The BRAVA data together with the RCGs 3D density imply a low pattern speed for the Galactic B/P bulge of 25-30 km.s-1.kpc-1. This would place the Galaxy among the slow rotators (R >= 1.5). Finally, we show that the Milky Way's B/P bulge has an off-centred X structure, and that the stellar mass involved in the peanut shape accounts for at least 20% of the stellar mass of the bulge, significantly larger than previously thought.
We present the formation of a Kinematically Decoupled Core (KDC) in an elliptical galaxy, resulting from a major merger simulation of two disk galaxies. We show that although the two progenitor galaxies are initially following a prograde orbit, strong reactive forces during the merger can cause a short-lived change of their orbital spin; the two progenitors follow a retrograde orbit right before their final coalescence. This results in a central kinematic decoupling and the formation of a large-scale (~2 kpc radius) counter-rotating core (CRC) at the center of the final elliptical-like merger remnant (M*=1.3x10^11 Msun), while its outer parts keep the rotation direction of the initial orbital spin. The stellar velocity dispersion distribution of the merger remnant galaxy exhibits two symmetrical off-centered peaks, comparable to the observed "2-sigma galaxies". The KDC/CRC consists mainly of old, pre-merger population stars (older than 5 Gyr), remaining prominent in the center of the galaxy for more than 2 Gyr after the coalescence of its progenitors. Its properties are consistent with KDCs observed in massive elliptical galaxies. This new channel for the formation of KDCs from prograde mergers is in addition to previously known formation scenarios from retrograde mergers and can help towards explaining the substantial fraction of KDCs observed in early-type galaxies.
We develop a general framework for data analysis and phenomenology of the CMB four-point function or trispectrum. To lowest order in the derivative expansion, the inflationary action admits three quartic operators consistent with symmetry: $\dot\sigma^4$, $\dot\sigma^2 (\partial\sigma^2)$, and $(\partial\sigma)^4$. In single field inflation, only the first of these operators can be the leading non-Gaussian signal. A Fisher matrix analysis shows that there is one near-degeneracy among the three CMB trispectra, so we parameterize the trispectrum with two coefficients $g_{NL}^{\dot\sigma^4}$ and $g_{NL}^{(\partial\sigma)^4}$, in addition to the coefficient $g_{NL}^{\rm loc}$ of $\zeta^3$-type local non-Gaussianity. This three-parameter space is analogous to the parameter space $(f_{NL}^{\rm loc}, f_{NL}^{\rm equil}, f_{NL}^{\rm orth})$ commonly used to parameterize the CMB three-point function. We next turn to data analysis and show how to represent these trispectra in a factorizable form which leads to computationally fast operations such as evaluating a CMB estimator or simulating a non-Gaussian CMB. We discuss practical issues in CMB analysis pipelines, and perform an optimal analysis of WMAP data. Our minimum-variance estimates are $g_{NL}^{\rm loc} = (-3.80 \pm 2.19) \times 10^5$, $g_{NL}^{\dot\sigma^4} = (-3.20 \pm 3.09) \times 10^6$, and $g_{NL}^{(\partial\sigma)^4} = (-10.8 \pm 6.33) \times 10^5$ after correcting for the effects of CMB lensing. No evidence of a nonzero inflationary four-point function is seen.
Recent suggestions of a "photon underproduction crisis" (Kollmeier et al 2014) have generated concern over the intensity and spectrum of ionizing photons in the metagalactic ultraviolet background (UVB). The balance of hydrogen photoionization and recombination determines the opacity of the low-redshift intergalactic medium (IGM). We calibrate the hydrogen photoionization rate ($\Gamma_H$) by comparing Hubble Space Telescope spectroscopic surveys of the low-redshift column density distribution of H I absorbers to new cosmological simulations. The distribution, $f(N_{HI}, z) = d^2N/d(log N_{HI}) dz$, is consistent with an increased UVB that includes contributions from both quasars and galaxies. Our recommended fit, $\Gamma_H(z) = (4.6x10^{-14}~s^{-1})(1+z)^{4.4}$ for $0 < z < 0.5$, corresponds to unidirectional LyC photon flux $\Phi_0 = 5700$ cm$^{-2}$ s$^{-1}$ at z = 0. This flux agrees with observed IGM metal ionization ratios (C III / C IV and Si III / Si IV) and suggests a 25-30\% contribution of the Lya absorbers to the cosmic baryon inventory. The primary uncertainties in the low-redshift UVB are the contribution from massive stars in galaxies and the LyC escape fraction ($f_{esc}$), a highly directional quantity that is difficult to constrain statistically. We suggest that low-mass starburst galaxies are important contributors to the ionizing UVB at z < 2. Their additional flux would resolve any crisis in photon underproduction.
We have developed a near-infrared camera called ANIR (Atacama Near-InfraRed camera) for the University of Tokyo Atacama Observatory 1.0m telescope (miniTAO) installed at the summit of Cerro Chajnantor (5640 m above sea level) in northern Chile. The camera provides a field of view of 5'.1 $\times$ 5'.1 with a spatial resolution of 0".298 /pixel in the wavelength range of 0.95 to 2.4 $\mu$m. Taking advantage of the dry site, the camera is capable of hydrogen Paschen-$\alpha$ (Pa$\alpha$, $\lambda=$1.8751 $\mu$m in air) narrow-band imaging observations, at which wavelength ground-based observations have been quite difficult due to deep atmospheric absorption mainly from water vapor. We have been successfully obtaining Pa$\alpha$ images of Galactic objects and nearby galaxies since the first-light observation in 2009 with ANIR. The throughputs at the narrow-band filters ($N1875$, $N191$) including the atmospheric absorption show larger dispersion (~10%) than those at broad-band filters (a few %), indicating that they are affected by temporal fluctuations in Precipitable Water Vapor (PWV) above the site. We evaluate the PWV content via the atmospheric transmittance at the narrow-band filters, and derive that the median and the dispersion of the distribution of the PWV are 0.40+/-0.30 mm for $N1875$ and 0.37+/-0.21 mm for $N191$, which are remarkably smaller (49+/-38% for $N1875$ and 59+/-26% for $N191$) than radiometry measurements at the base of Cerro Chajnantor (5100 m alt.). The decrease in PWV can be explained by the altitude of the site when we assume that the vertical distribution of the water vapor is approximated at an exponential profile with scale heights within 0.3-1.9 km (previously observed values at night). We thus conclude that miniTAO/ANIR at the summit of Cerro Chajnantor indeed provides us an excellent capability for a "ground-based" Pa$\alpha$ observation.
The origin of the Cosmic Infrared Background (CIB) is still poorly understood, and represents a challenge from both the theoretical and observational points of view. We analysed 18 ALMA continuum maps in band 6 and 7, with rms down to 7.8 $\mu$Jy, to estimate differential number counts down to 60 $\mu$Jy and 100 $\mu$Jy at $\lambda =$1.3 mm and $\lambda =$1.1 mm, respectively. We improved the source extraction method to detect sources down to S/N = 3.5. We detected 50 faint sources ($<$ 1 mJy) down to 60 $\mu$Jy. Determining the fraction of CIB resolved by the ALMA observations is hampered by the large uncertainties plaguing the CIB measurements (a factor of four in flux). However, our results provide a solid lower limit to CIB intensity. Moreover, the flattening of the integrated number counts at faint fluxes strongly suggest that we are close to resolving 100% of the CIB. Our data imply that galaxies with $\rm SFR < 40 ~M_{\odot}/yr$ certainly contribute less than 50% to the CIB, and probably a much lower fraction, while the bulk of the CIB must be produced by galaxies with $\rm SFR > 40 ~M_{\odot}/yr$. The differential number counts are in nice agreement with recent semi-analytical models of galaxy formation down to our faint fluxes, therefore supporting the galaxy evolutionary scenarios and assumptions made in these models.
We present an analysis of archival Chandra observations of the mixed-morphology remnant 3C400.2. We analysed spectra of different parts of the remnant to observe if the plasma properties provide hints on the origin of the mixed-morphology class. These remnants often show overionization, which is a sign of rapid cooling of the thermal plasma, and super-solar abundances of elements which is a sign of ejecta emission. Our analysis shows that the thermal emission of 3C400.2 can be well explained by a two component non-equilibrium ionization model, of which one component is underionized, has a high temperature ($kT \approx 3.9$ keV) and super-solar abundances, while the other component has a much lower temperature ($kT \approx 0.14$ keV), solar abundances and shows signs of overionization. The temperature structure, abundance values and density contrast between the different model components suggest that the hot component comes from ejecta plasma, while the cooler component has an interstellar matter origin. This seems to be the first instance of an overionized plasma found in the outer regions of a supernova remnant, whereas the ejecta component of the inner region is underionized. In addition, the non-ionization equilibrium plasma component associated with the ejecta is confined to the central, brighter parts of the remnant, whereas the cooler component is present mostly in the outer regions. Therefore our data can most naturally be explained by an evolutionary scenario in which the outer parts of the remnant are cooling rapidly due to having swept up high density ISM, while the inner parts are very hot and cooling inefficiently due to low density of the plasma. This is also known as the relic X-ray scenario.
The Keck Array is a system of cosmic microwave background (CMB) polarimeters, each similar to the BICEP2 experiment. In this paper we report results from the 2012 and 2013 observing seasons, during which the Keck Array consisted of five receivers all operating in the same (150 GHz) frequency band and observing field as BICEP2. We again find an excess of B-mode power over the lensed-$\Lambda$CDM expectation of $> 5 \sigma$ in the range $30 < \ell < 150$ and confirm that this is not due to systematics using jackknife tests and simulations based on detailed calibration measurements. In map difference and spectral difference tests these new data are shown to be consistent with BICEP2. Finally, we combine the maps from the two experiments to produce final Q and U maps which have a depth of 57 nK deg (3.4 $\mu$K arcmin) over an effective area of 400 deg$^2$ for an equivalent survey weight of 250,000 $\mu$K$^{-2}$. The final BB band powers have noise uncertainty a factor of 2.3 times better than the previous results, and a significance of detection of excess power of $> 6\sigma$.
The debate on the nature of the gamma-ray emission from young supernova remnants is still open. Ascribing such emission to hadronic rather than leptonic processes would provide an evidence for the acceleration of protons and nuclei, and this fact would fit with the very popular (but not proven) paradigm that supernova remnants are the sources of Galactic cosmic rays. Here, we discuss this issue with a particular focus on the best studied gamma-ray-bright supernova remnant: RX~J1713.7-3946.
We present observations by the Fermi Gamma-Ray Space Telescope Gamma-Ray Burst Monitor (GBM) of the nearby (z=0.55) GRB 101219B. This burst is a long GRB, with an associated supernova and with a blackbody component detected in the early afterglow observed by the Swift X-ray Telescope (XRT). Here we show that the prompt gamma-ray emission has a blackbody spectrum, making this the second such burst observed by Fermi GBM. The properties of the blackbody, together with the redshift and our estimate of the radiative efficiency, makes it possible to calculate the absolute values of the properties of the outflow. We obtain an initial Lorentz factor Gamma=138\pm 8, a photospheric radius r_phot=4.4\pm 1.9 \times 10^{11} cm and a launch radius r_0=2.7\pm 1.6 \times 10^{7} cm. The latter value is close to the event horizon for a stellar-mass black hole and suggests that the jet has a relatively unobstructed path through the star. There is no smooth connection between the blackbody components seen by GBM and XRT, ruling out the scenario that the late emission is due to high-latitude effects. In the interpretation that the XRT blackbody is prompt emission due to late central engine activity, the jet either has to be very wide or have a clumpy structure where the emission originates from a small patch. Other explanations for this component, such as emission from a cocoon surrounding the jet, are also possible.
We present a "kiloday" optical nebular spectrum of the nearby Type Ia supernova (SN) 2011fe, obtained +981 days after explosion. SN 2011fe exhibits little evolution since the +593 day optical spectrum, but there are several curious aspects in this new extremely late-time regime. We see new emission features at 5800 and 6400 \AA\ that could be attributed to [Fe I], but also consider (and reject) the possibility that the latter is the high-velocity hydrogen predicted by Mazzali et al. (2014). The nebular feature at 5200 \AA\ exhibits a linear velocity evolution of $\sim350$ $\rm km\ s^{-1}$ per 100 days from at least +220 to +980 days, but the line's shape also changes in this time, suggesting that line blending contributes to the evolution. At $\sim 1000$ days after explosion, flux from the SN has declined to a point where contribution from a hot, luminous secondary could be detected. In this work we make the first observational tests for a post-impact remnant star and constrain its luminosity to $<(3-5) \times 10^3$ $\rm L_{\odot}$. Additionally, we do not see any evidence for narrow H$\alpha$ emission, placing an upper limit of $3 \times 10^{-4}$ $\rm M_{\odot}$ on the mass of hydrogen in the system. We conclude that observations continue to strongly exclude many single-degenerate scenarios for SN 2011fe.
In this chapter we review observations and techniques for measuring both bulk flows in prominences and prominence mass. Measuring these quantities is essential to development and testing of models discussed throughout this book. Prominence flows are complex and various, ranging from the relatively linear flows along prominence spines to the complex, turbulent patterns exhibited by hedgerow prominences. Techniques for measuring flows include time slice and optical flow techniques used for motions in the plane of the sky and the use of spectral line profiles to determine Doppler velocities along the line of sight. Prominence mass measurement is chiefly done via continuum absorption measurements, but mass has also been estimated using cloud modeling and white light measurements.
There is growing evidence that when magnetic reconnection occurs in high Lundquist number plasmas such as in the Solar Corona or the Earth's Magnetosphere it does so within a fragmented, rather than a smooth current layer. Within the extent of these fragmented current regions the associated magnetic flux transfer and energy release occurs simultaneously in many different places. This investigation focusses on how best to quantify the rate at which reconnection occurs in such layers. An analytical theory is developed which describes the manner in which new connections form within fragmented current layers in the absence of magnetic nulls. It is shown that the collective rate at which new connections form can be characterized by two measures; a total rate which measures the true rate at which new connections are formed and a net rate which measures the net change of connection associated with the largest value of the integral of $E_{\|}$ through all of the non-ideal regions. Two simple analytical models are presented which demonstrate how each should be applied and what they quantify.
We report the first results from a survey for 1665, 1667, and 1720 MHz OH emission over a small region of the Outer Galaxy centered at $l \approx 105.0\deg , b \approx +1.0\deg$ . This sparse, high-sensitivity survey ($\Delta Ta \approx \Delta Tmb \approx 3.0 - 3.5$ mK rms in 0.55 km/s channels), was carried out as a pilot project with the Green Bank Telescope (GBT, FWHM $\approx 7.6'$) on a 3 X 9 grid at $0.5\deg$ spacing. The pointings chosen correspond with those of the existing $^{12}$CO(1-0) CfA survey of the Galaxy (FWHM $\approx 8.4'$). With 2-hr integrations, 1667 MHz OH emission was detected with the GBT at $\gtrsim 21$ of the 27 survey positions ($\geq 78\%$ ), confirming the ubiquity of molecular gas in the ISM as traced by this spectral line. With few exceptions, the main OH lines at 1665 and 1667 MHz appear in the ratio of 5:9 characteristic of LTE at our sensitivity levels. No OH absorption features are recorded in the area of the present survey, in agreement with the low levels of continuum background emission in this direction. At each pointing the OH emission appears in several components extending over a range of radial velocity and coinciding with well-known features of Galactic structure such as the Local Arm and the Perseus Arm. In contrast, little CO emission is seen in the survey area; less than half of the $\gtrsim 50$ identified OH components show detectable CO at the CfA sensitivity levels, and these are generally faint. There are no CO profiles without OH emission. With few exceptions, peaks in the OH profiles coincide with peaks in the GBT HI spectra (obtained concurrently, FWHM $8.9'$), although the converse is not true. We conclude that main-line OH emission is a promising tracer for the "dark molecular gas" in the Galaxy discovered earlier in Far-IR and gamma-ray emission. Further work is needed to establish the quantitative details of this connection.
Wide area X-ray and far infrared surveys are a fundamental tool to investigate the link between AGN growth and star formation, especially in the low-redshift universe (z<1). The Herschel Terahertz Large Area survey (H-ATLAS) has covered 550 deg^2 in five far-infrared and sub-mm bands, 16 deg^2 of which have been presented in the Science Demonstration Phase (SDP) catalogue. Here we introduce the XMM-Newton observations in H-ATLAS SDP area, covering 7.1 deg^2 with flux limits of 2e-15, 6e-15 and 9e-15 erg/s/cm^2 in the 0.5--2, 0.5--8 and 2--8 keV bands, respectively. We present the source detection and the catalogue, which includes 1700, 1582 and 814 sources detected by Emldetect in the 0.5--8, 0.5--2 and 2--8 keV bands, respectively; the number of unique sources is 1816. We extract spectra and derive fluxes from power-law fits for 398 sources with more than 40 counts in the 0.5--8 keV band. We compare the best-fit fluxes with the catalogue ones, obtained by assuming a common photon index of Gamma=1.7; we find no bulk difference between the fluxes, and a moderate dispersion of s=0.33 dex. Using wherever possible the fluxes from the spectral fits, we derive the 2--10 keV LogN-LogS, which is consistent with a Euclidean distribution. Finally, we release computer code for the tools developed for this project.
We have observed the pulsar in the Crab Nebula at high radio frequencies and high time resolution. We present continuously sampled data at 640-ns time resolution, and individual bright pulses recorded at down to 0.25-ns time resolution. Combining our new data with previous data from our group and from the literature shows the dramatic changes in the pulsar's radio emission between low and high radio frequencies. Below about 5 GHz the mean profile is dominated by the bright Main Pulse and Low-Frequency Interpulse. Everything changes, however, above about 5 GHz; the Main Pulse disappears, the mean profile of the Crab pulsar is dominated by the High-Frequency Interpulse (which is quite different from its low-frequency counterpart) and the two High-Frequency Components. We present detailed observational characteristics of these different components which future models of the pulsar's magnetosphere must explain.
The knowledge of stellar ages directly impacts the characterization of a planetary system as it puts strong constraints on the moment when the system was born. Unfortunately, the determination of precise stellar ages is a very difficult task. Different methods can be used to do so (based on isochrones or chemical element abundances) but they usually provide large uncertainties. During its evolution a star goes through processes leading to loss of angular momentum but also changes in its magnetic activity. Building rotation, magnetic, age relations would be an asset to infer stellar ages model independently. Several attempts to build empirical relations between rotation and age (namely gyrochronology) were made with a focus on cluster stars where the age determination is easier and for young stars on the main sequence. For field stars, we can now take advantage of high-precision photometric observations where we can perform asteroseismic analyses to improve the accuracy of stellar ages. Furthermore, the variability in the light curves allow us to put strong constraints on the stellar rotation and magnetic activity. By combining these precise measurements, we are on the way of understanding and improving relations between magnetic activity, rotation, and age, in particular at different stages of stellar evolution. I will review the status on gyrochronology relationships based on observations of young cluster stars. Then I will focus on solar-like stars and describe the inferences on stellar ages, rotation, and magnetism that can be provided by high-quality photometric observations such as the ones of the Kepler mission, in particular through asteroseismic analyses.
How can we select a cohort of promising astrophysicists before they have made their discoveries? This is the fundamental challenge of academic planning. I argue that science can only blossom if young researchers are rewarded for acquired skills and growth rather than inherited academic ancestry.
We report the Fermi Large Area Telescope (LAT) detection of a $\gamma$-ray source at the position of SAX J1808.4$-$3658. This transient low-mass X-ray binary contains an accreting millisecond puslar, which is only seen during its month-long outbursts and likely switches to be rotation powered during its quiescent state. Emission from the $\gamma$-ray source can be described by a power law with an exponential cutoff, the characteristic form for pulsar emission. Folding the source's 2.0--300 GeV photons at the binary orbital period, a weak modulation is seen (with an H-test value of $\sim$17). In addition, three sets of archival XMM-Newton data for the source field are analyzed, and we find only one X-ray source with 3--4$\sigma$ flux variations in the 2$\sigma$ error circle of the $\gamma$-ray source. However based on the X-ray properties, this X-ray source is not likely a background AGN, the major class of Fermi sources detected by LAT. These results support the possible association between the $\gamma$-ray source and SAX J1808.4$-$3658 and thus the scenario that the millisecond pulsar is rotation powered in the quiescent state. Considering a source distance of 3.5 kpc for SAX J1808.4$-$3658, the 0.1--300 GeV luminosity is 5.7$\times 10^{33}$ erg s$^{-1}$, implying a $\gamma$-ray conversion efficiency of 63\% for the pulsar in this binary.
In a single optical spectrum, the quasar narrow-line region (NLR) reveals low density, photoionized gas in the host galaxy interstellar medium, while the immediate vicinity of the central engine generates the accretion disk continuum and broad emission lines. To isolate these two components, we construct a library of high S/N optical composite spectra created from the Sloan Digital Sky Survey (SDSS-DR7). We divide the sample into bins of continuum luminosity and Hbeta FWHM that are used to construct median composites at different redshift steps up to 0.75. We measure the luminosities of the narrow-emission lines [NeV]3427, [NeIII]3870, [OIII]5007, and [OII]3728 with ionization potentials (IPs) of 97, 40, 35, and 13.6 eV respectively. The high IP lines' luminosities show no evidence of increase with redshift consistent with no evolution in the AGN SED or the host galaxy ISM illuminated by the continuum. In contrast, we find that the [OII] line becomes stronger at higher redshifts, and we interpret this as a consequence of enhanced star formation contributing to the [OII] emission in host galaxies at higher redshifts. The SFRs estimated from the [OII] luminosities show a flatter increase with z than non-AGN galaxies given our assumed AGN contribution to the [OII] luminosity. Finally, we confirm an inverse correlation between the strength of the FeII4570 complex and both the [OIII] EW (though not the luminosity) and the width of the Hbeta line as known from the eigenvector 1 correlations.
We explore the connection between the central supermassive blackholes (SMBH) in galaxies and the dark matter halo through the relation between the masses of the SMBHs and the maximum circular velocities of the host galaxies, as well as the relationship between stellar velocity dispersion of the spheroidal component and the circular velocity. Our assumption here is that the circular velocity is a proxy for the mass of the dark matter halo. We rely on a heterogeneous sample containing galaxies of all types. The only requirement is that the galaxy has a direct measurement of the mass of its SMBH and a direct measurement of its circular velocity and its velocity dispersion. Previous studies have analyzed the connection between the SMBH and dark matter halo through the relationship between the circular velocity and the bulge velocity dispersion, with the assumption that the bulge velocity dispersion stands in for the mass of the SMBH, via the well{}-established SMBH mass{}-bulge velocity dispersion relation. Using intermediate relations may be misleading when one is studying them to decipher the active ingredients of galaxy formation and evolution. We believe that our approach will provide a more direct probe of the SMBH and the dark matter halo connection. We find that the correlation between the mass of supermassive blackholes and the circular velocities of the host galaxies is extremely weak, leading us to state the dark matter halo may not play a major role in regulating the blackhole growth in the present Universe.
Stellar spectropolarimetry has become extremely popular during the last decade, and has led to major advances in our understanding of stellar magnetic fields and of their impact on stellar structure and evolution. Many important discoveries have been obtained thanks to observations performed with the FORS low-resolution spectropolarimeters of the ESO Very Large Telescope. We first review and summarise the major results of a homogeneous re-reduction and analysis of all single-slit FORS1 spectropolarimetric observations. This work revealed a non-negligible dependence of the results upon the adopted reduction and analysis procedure, as well as the presence of instabilities, revealing that photon noise is not the only source of uncertainty. As a consequence of our new analysis and assessment of the uncertainties, we are not able to confirm a large number of magnetic field detections presented in the past for a variety of stars. We further summarise the results of FORS2 spectropolarimetric observations of the A0 supergiant HD92207 which allowed us to explore further the nature of the instabilities, roughly constraining their maximum impact on the derived Stokes profiles and magnetic field values. We finally present new results obtained with a further independent pipeline on the FORS2 data of HD92207, confirming our previous analysis, and discuss simple quality-check controls which can be performed on the data in order to distinguish between genuine and spurious signals. All together, our results reveal that the FORS spectropolarimeters are indeed reliable instruments, when their capabilities are not pushed beyond the limits of a Cassegrain mounted low-resolution spectrograph.
The standard errors of the end-of-mission Gaia astrometry have been re-assessed after conclusion of the in-orbit commissioning phase of the mission. An analytical relation is provided for the parallax standard error as function of Gaia G magnitude (and V-I colour) which supersedes the pre-launch relation provided in de Bruijne (2012).
An established technique for the measurement of ultra-high-energy-cosmic-rays is the detection of the fluorescence light induced in the atmosphere of the Earth, by means of telescopes equipped with photomultiplier tubes. Silicon photomultipliers (SiPMs) promise an increase in the photon detection efficiency which outperforms conventional photomultiplier tubes. In combination with their compact package, a moderate bias voltage of several ten volt and single photon resolution, the use of SiPMs can improve the energy and spatial resolution of air fluorescence measurements, and lead to a gain in information on the primary particle. Though, drawbacks like a high dark-noise-rate and a strong temperature dependency have to be managed. FAMOUS is a refracting telescope prototype instrumented with 64 SiPMs of which the main optical element is a Fresnel lens of 549.7 mm diameter and 502.1 mm focal length. The sensitive area of the SiPMs is increased by a special light collection system consisting of Winston cones. The total field of view of the telescope is approximately 12 $^\circ$. The frontend electronics automatically compensates for the temperature dependency of the SiPMs and will provide trigger information for the readout. Already for this prototype, the Geant4 detector simulation indicates full detection efficiency of extensive air showers of $E=10^{18}\,\text{eV}$ up to a distance of 6 km. We present the first working version of FAMOUS with a focal plane prototype providing seven active pixels.
On 14 October 2014 the Sun Watcher with Active Pixels and Image Processing (SWAP) EUV solar telescope on-board the Project for On-Board Autonomy 2 (PROBA2) spacecraft observed an eruption that led to the formation of perhaps the largest post-eruptive loop system seen in the solar corona in solar cycle 24. The initial eruption occurred at about 18:30 UT on 14 October, behind the East Solar limb, and was observed as a a coronal mass ejection and an M2.2 solar flare. In the 48 hours following the eruption, the associated post-eruptive loops grew to a height of approximately 400000 km (>0.5 solar-radii) at rates between 2-6 km/s. We conclude from our observations of this event that ordinary post-eruptive loops and so-called post-flare giant arches are fundamentally the same and are formed by the same magnetic reconnection mechanism.
The search for extrasolar planets has developed rapidly and, today, more than 1700 planets have been found orbiting stars. Thanks to Gaia, we will collect high-accuracy astrometric orbits of thousands of new low-mass celestial objects, such as extra-solar planets and brown dwarfs. These measurements in combination with spectroscopy and with present day and future extrasolar planet search programs (like HARPS, ESPRESSO) will have a crucial contribution to several aspects of planetary astrophysics (formation theories, dynamical evolution, etc.). Moreover, Gaia will have a strong contribution on the stellar chemical and kinematic characterisation studies. In this paper we present a short overview of the importance of Gaia in the context of exoplanet research. As preparatory work for Gaia, we will then present a study where we derived stellar parameters for a sample of field giant stars.
We develop a method to estimate photometric metallicities by simultaneously fitting the dereddened colors u-g, g-r, r-i and i-z from the SDSS with those predicted by the metallicity-dependent stellar loci. The method is tested with a spectroscopic sample of main-sequence stars in Stripe 82 selected from the SDSS DR9 and three open clusters. With 1 per cent photometry, the method is capable of delivering photometric metallicities precise to about 0.05, 0.12, and 0.18 dex at metallicities of 0.0, -1.0, and -2.0, respectively, comparable to the precision achievable with low-resolution spectroscopy at a signal-to-noise ratio of 10. We apply this method to the re-calibrated Stripe 82 catalog and derive metallicities for about 0.5 million stars of colors 0.3 < g-i < 1.6 mag and distances between 0.3 -- 18 kpc. Potential systematics in the metallicities thus derived, due to the contamination of giants and binaries, are investigated. Photometric distances are also calculated. About 91, 72, and 53 per cent of the sample stars are brighter than r = 20.5, 19.5, and 18.5 mag, respectively. The median metallicity errors are around 0.19, 0.16, 0.11, and 0.085 dex for the whole sample, and for stars brighter than r = 20.5, 19.5, and 18.5 mag, respectively. The median distance errors are 8.8, 8.4, 7.7, and 7.3 per cent for the aforementioned four groups of stars, respectively. The data are publicly available. Potential applications of the data in studies of the distribution, (sub)structure, and chemistry of the Galactic stellar populations, are briefly discussed. The results will be presented in future papers.
We present a novel approach to test the consistency of the cosmological models with multiband CMB data using a nonparametric approach. In our analysis we calibrate the REACT (Risk Estimation and Adaptation after Coordinate Transformation) confidence levels associated with distances in function space (confidence distances) based on the Monte Carlo simulations in order to test the consistency of an assumed cosmological model with observation. To show the applicability of our algorithm, we confront Planck 2013 temperature data with concordance model of cosmology considering two different Planck spectra combination. In order to have an accurate quantitative statistical measure to compare between the data and the theoretical expectations, we calibrate REACT confidence distances and perform a bias control using many realizations of the data. Our results in this work using Planck 2013 temperature data put the best fit $\Lambda$CDM model at $95\% (\sim 2\sigma)$ confidence distance from the center of the nonparametric confidence set which hints towards considerable inconsistency. Repeating the analysis excluding the Planck $217 \times 217$ GHz spectrum data, the best fit $\Lambda$CDM model shifts to $70\% (\sim 1\sigma)$ confidence distance from the center of the nonparametric confidence set. The most prominent features in the data deviating from the best fit $\Lambda$CDM model seems to be at low multipoles $ 18 < \ell < 26$ at greater than $2\sigma$, $\ell \sim 750$ at $\sim1$ to $2\sigma$ and $\ell \sim 1800$ at greater than $2\sigma$ level. Excluding the $217\times217$ GHz spectrum the feature at $\ell \sim 1800$ becomes substantially less significance at $\sim1$ to $2 \sigma$ confidence level. Results of our analysis based on the new approach we propose in this work are in agreement with other analysis done using alternative methods.
The non-Maxwellian $\kappa$-distributions have been detected in the solar transition region and flares. These distributions are characterized by a high-energy tail and a near-Maxwellian core and are known to have significant impact on the resulting optically thin spectra arising from collisionally dominated astrophysical plasmas. We developed the KAPPA package (this http URL) for synthesis of such line and continuum spectra. The package is based on the freely available CHIANTI database and software, and can be used in a similar manner. Ionization and recombination rates together with the ionization equilibria are provided for a range of $\kappa$ values. Distribution-averaged collision strengths for excitation are obtained by an approximate method for all transitions in all ions available within CHIANTI. The validity of this approximate method is tested by comparison with direct calculations. Typical precisions of better than 5% are found, with all cases being within 10%. Tools for calculation of synthetic line and continuum intensities are provided and described. Examples of the synthetic spectra and SDO/AIA responses to emission for the $\kappa$-distributions are given.
The UVMag consortium proposed the space mission project Arago to ESA at its M4 call. It is dedicated to the study of the dynamic 3D environment of stars and planets. This space mission will be equipped with a high-resolution spectropolarimeter working from 119 to 888 nm. A preliminary optical design of the whole instrument has been prepared and is presented here. The design consists of the telescope, the instrument itself, and the focusing optics. Considering not only the scienti?c requirements, but also the cost and size constraints to ?t a M-size mission, the telescope has a 1.3 m diameter primary mirror and is a classical Cassegrain-type telescope that allows a polarization-free focus. The polarimeter is placed at this Cassegrain focus. This is the key element of the mission and the most challenging to be designed. The main challenge lies in the huge spectral range offered by the instrument; the polarimeter has to deliver the full Stokes vector with a high precision from the FUV (119 nm) to the NIR (888 nm). The polarimeter module is then followed by a high-resolution echelle-spectrometer achieving a resolution of 35000 in the visible range and 25000 in the UV. The two channels are separated after the echelle grating, allowing a speci?c cross-dispersion and focusing optics for the UV and visible ranges. Considering the large ?eld of view and the high numerical aperture, the focusing optic for both the UV and visible channels is a Three-Mirror-Anastigmat (TMA) telescope, in order to focus the various wavelengths and many orders onto the detectors.
Increasing data volumes delivered by a new generation of radio interferometers require computationally efficient and robust calibration algorithms. In this paper, we propose distributed calibration as a way of improving both computational cost as well as robustness in calibration. We exploit the data parallelism across frequency that is inherent in radio astronomical observations that are recorded as multiple channels at different frequencies. Moreover, we also exploit the smoothness of the variation of calibration parameters across frequency. Data parallelism enables us to distribute the computing load across a network of compute agents. Smoothness in frequency enables us reformulate calibration as a consensus optimization problem. With this formulation, we enable flow of information between compute agents calibrating data at different frequencies, without actually passing the data, and thereby improving robustness. We present simulation results to show the feasibility as well as the advantages of distributed calibration as opposed to conventional calibration.
Glycolaldehyde is a key molecule in the formation of biologically relevant molecules such as ribose. We report its detection with the Plateau de Bure interferometer towards the Class 0 young stellar object NGC1333 IRAS2A, which is only the second solar-type protostar for which this prebiotic molecule is detected. Local thermodynamic equilibrium analyses of glycolaldehyde, ethylene glycol (the reduced alcohol of glycolaldehyde) and methyl formate (the most abundant isomer of glycolaldehyde) were carried out. The relative abundance of ethylene glycol to glycolaldehyde is found to be ~5 -higher than in the Class 0 source IRAS 16293-2422 (~1), but comparable to the lower limits derived in comets ($\geq$3-6). The different ethylene glycol-to-glycolaldehyde ratios in the two protostars could be related to different CH3OH:CO compositions of the icy grain mantles. In particular, a more efficient hydrogenation on the grains in NGC 1333 IRAS2A would favor the formation of both methanol and ethylene glycol. In conclusion, it is possible that, like NGC 1333 IRAS2A, other low-mass protostars show high ethylene glycol-to-glycolaldehyde abundance ratios. The cometary ratios could consequently be inherited from earlier stages of star formation, if the young Sun experienced conditions similar to NGC1333 IRAS2A.
The Euclid space mission proposes to survey 15000 square degrees of the extragalactic sky during 6 years, with a step-and-stare technique. The scheduling of observation sequences is driven by the primary scientific objectives, spacecraft constraints, calibration requirements and physical properties of the sky. We present the current reference implementation of the Euclid survey and on-going work on survey optimization.
Photometric observations in V and I bands of six eclipsing binaries at the lower limit of the orbital periods of W UMa stars are presented. Three of them are newly discovered eclipsing systems. The light curve solutions revealed that all short-period targets were contact or overcontact binaries and added new six binaries to the family of short-period systems with estimated parameters. Four binaries have equal in size components and mass ratio near 1. The phase variability of the V-I colors of all targets may be explained by lower temperatures of their back surfaces than those of their side surfaces. Five systems revealed O'Connell effect that was reproduced by cool spots on the side surfaces of their primary components. The light curves of V1067 Her in 2011 and 2012 were fitted by diametrically opposite spots. The applying of the criteria for subdivision of the W UMa stars to our targets led to ambiguous results.
The light curves of 2768 stars with effective temperatures and surface gravities placing them near the gamma Doradus/delta Scuti instability region were observed as part of the Kepler Guest Observer program from Cycles 1 through 5. The light curves were analyzed in a uniform manner to search for gamma Doradus, delta Scuti, and hybrid star pulsations. The gamma Doradus, delta Scuti, and hybrid star pulsations extend asteroseismology to stars slightly more massive (1.4 to 2.5 solar masses) than our Sun. We find 207 gamma Doradus, 84 delta Scuti, and 32 hybrid candidate stars. Many of these stars are cooler than the red edge of the gamma Doradus instability strip as determined from ground-based observations made before Kepler. A few of our gamma Doradus candidate stars lie on the hot side of the ground-based gamma Doradus instability strip. The hybrid candidate stars cover the entire region between 6200 K and the blue edge of the ground-based delta Scuti instability strip. None of our candidate stars are hotter than the hot edge of the ground-based delta Scuti instability strip. Our discoveries, coupled with the work of others, shows that Kepler has discovered over 2000 gamma Doradus, delta Scuti, and hybrid star candidates in the 116 square degree Kepler field of view. We found relatively few variable stars fainter than magnitude 15, which may be because they are far enough away to lie between spiral arms in our Galaxy, where there would be fewer stars.
Supernova generated shock waves are responsible for most of the destruction of dust grains in the interstellar medium (ISM). Calculations of the dust destruction timescale have so far been carried out using plane parallel steady shocks, however that approximation breaks down when the destruction timescale becomes longer than that for the evolution of the supernova remnant (SNR) shock. In this paper we present new calculations of grain destruction in evolving, radiative SNRs. To facilitate comparison with the previous study by Jones et al. (1996), we adopt the same dust properties as in that paper. We find that the efficiencies of grain destruction are most divergent from those for a steady shock when the thermal history of a shocked gas parcel in the SNR differs significantly from that behind a steady shock. This occurs in shocks with velocities >~ 200 km/s for which the remnant is just beginning to go radiative. Assuming SNRs evolve in a warm phase dominated ISM, we find dust destruction timescales are increased by a factor of ~2 compared to those of Jones et al. (1996), who assumed a hot gas dominated ISM. Recent estimates of supernova rates and ISM mass lead to another factor of ~3 increase in the destruction timescales, resulting in a silicate grain destruction timescale of ~2-3 Gyr. These increases, while not able resolve the problem of the discrepant timescales for silicate grain destruction and creation, are an important step towards understanding the origin, and evolution of dust in the ISM.
Currently there is no definitive description for the accelerated expansion of the Universe at both early and late times; we know these two periods as the epochs of inflation and dark energy. Contained within this Thesis are two studies of inflation and one in the context of dark energy. The first study involves two noncanonical kinetic terms each in a two-field scenario, and their effects on the generation of isocurvature modes. As a result, these terms affect the isocurvature perturbations produced, and consequently the Cosmic Microwave Background. In the following study, the impact of a sharp transition upon the effective Planck mass is considered in both a single-field and two-field model. A feature in the primordial power spectrum arising from these transitions is found in single-field models, but not for two-field models. The final model discussed is on the subject of dark energy. A type of nonconformal coupling is examined namely the "disformal" coupling; in this scenario a scalar field is disformally coupled to matter species. Two consistency checks are undertaken, the first to provide a fluid description and the second, a kinetic theory. From this, observables are constructed and used to create constraints on the individual coupling strengths.
Context: The gamma-ray binary LSI+61$^{\circ}$303 shows multiple periodicities. The timing analysis of 6.7 yr of GBI radio data and of 6 yr of Fermi-LAT GeV gamma-ray data both have found two close periodicities $P_{1, {\rm GBI}} = 26.49 \pm 0.07$ d, $P_{2, {\rm GBI}} = 26.92 \pm 0.07$ d and $P_{1, \gamma} = 26.48 \pm 0.08$ d, $P_{2, \gamma} = 26.99 \pm 0.08$ d. Aims: The system LSI+61$^{\circ}$303 is the object of several continuous monitoring programs at low and high energies. The frequency difference between $\nu_1$ and $\nu_2$ of only 0.0006~d$^{-1}$ requires long-term monitoring because the frequency resolution in timing analysis is related to the inverse of the overall time interval. The Owens Valley Radio Observatory (OVRO) 40~m telescope has been monitoring the source at 15~GHz for five~years and overlaps with Fermi-LAT monitoring. The aim of this work is to establish whether the two frequencies are also resolved in the OVRO monitoring. Methods: We analysed OVRO data with the Lomb-Scargle method. We also updated the timing analysis of Fermi-LAT observations. Results: The periodograms of OVRO data confirm the two periodicities $P_{1, {\rm OVRO}} = 26.5 \pm 0.1$ d and $P_{2, {\rm OVRO}} = 26.9 \pm 0.1$ d. Conclusions: The three indipendent measurements of $P_1$ and $P_2$ with GBI, OVRO, and Fermi-LAT observations confirm that the periodicities are permanent features of the system LSI+61$^{\circ}$303. The similar behaviours of the emission at high (GeV) and low (radio) energy when the compact object in LSI+61$^{\circ}$303 is toward apastron suggest that the emission is caused by the same periodically ($P_1$) ejected population of electrons in a precessing ($P_2$) jet.
The nature of particle acceleration at the Sun, whether through flare reconnection processes or through shocks driven by coronal mass ejections (CMEs), is still under scrutiny despite decades of research. The measured properties of solar energetic particles (SEPs) have long been modeled in different particle-acceleration scenarios. The challenge has been to disentangle to the effects of transport from those of acceleration. The Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) instrument, enables unique observations of SEPs including composition and the angular distribution of the particles about the magnetic field, i.e. pitch angle distribution, over a broad energy range (>80 MeV) -- bridging a critical gap between space-based measurements and ground-based. We present high-energy SEP data from PAMELA acquired during the 2012 May 17 SEP event. These data exhibit differential anisotropies and thus transport features over the instrument rigidity range. SEP protons exhibit two distinct pitch angle distributions; a low-energy population that extends to 90{\deg} and a population that is beamed at high energies (>1 GeV), consistent with neutron monitor measurements. To explain a low-energy SEP population that exhibits significant scattering or redistribution accompanied by a high-energy population that reaches the Earth relatively unaffected by dispersive transport effects, we postulate that the scattering or redistribution takes place locally. We believe these are the first comprehensive measurements of the effects of solar energetic particle transport in the Earth's magnetosheath.
We used integrated-light medium-resolution spectra of six Galactic globular clusters and model stellar atmospheres to carry out population synthesis and to derive chemical composition and age of the clusters. We used medium-resolution spectra of globular clusters published by Schiavon et al. (2005), as well as our long-slit observations with the 1.93 m telescope of the Haute Provence Observatory. The observed spectra were fitted to the theoretical ones interactively. As an initial approach, we used masses, radii and log g of stars in the clusters corresponding to the best fitting isochrones in the observed color-magnitude diagrams. The computed synthetic blanketed spectra of stars were summed according to the Chabrier mass function. To improve the determination of age and helium content, the shape and depth of the Balmer absorption lines was analysed. The abundances of Mg, Ca, C and several other elements were derived. A reasonable agreement with the literature data both in chemical composition and in age of the clusters is found. Our method might be useful for the development of stellar population models and for a better understanding of extragalactic star clusters.
An asteroid family is typically formed when a larger parent body undergoes a catastrophic collisional disruption, and as such family members are expected to show physical properties that closely trace the composition and mineralogical evolution of the parent. Recently a number of new datasets have been released that probe the physical properties of a large number of asteroids, many of which are members of identified families. We review these data sets and the composite properties of asteroid families derived from this plethora of new data. We also discuss the limitations of the current data, and the open questions in the field.
The Flexible Image Transport System (FITS) standard has been a great boon to
astronomy, allowing observatories, scientists and the public to exchange
astronomical information easily. The FITS standard, however, is showing its
age. Developed in the late 1970s, the FITS authors made a number of
implementation choices that, while common at the time, are now seen to limit
its utility with modern data. The authors of the FITS standard could not
anticipate the challenges which we are facing today in astronomical computing.
Difficulties we now face include, but are not limited to, addressing the need
to handle an expanded range of specialized data product types (data models),
being more conducive to the networked exchange and storage of data, handling
very large datasets, and capturing significantly more complex metadata and data
relationships.
There are members of the community today who find some or all of these
limitations unworkable, and have decided to move ahead with storing data in
other formats. If this fragmentation continues, we risk abandoning the
advantages of broad interoperability, and ready archivability, that the FITS
format provides for astronomy. In this paper we detail some selected important
problems which exist within the FITS standard today. These problems may provide
insight into deeper underlying issues which reside in the format and we provide
a discussion of some lessons learned. It is not our intention here to prescribe
specific remedies to these issues; rather, it is to call attention of the FITS
and greater astronomical computing communities to these problems in the hope
that it will spur action to address them.
We report on the progress in characterization of the Hamamatsu model R11410-21 Photomultiplier tubes (PMTs) for XENON1T dark matter experiment. The absolute quantum efficiency (QE) of the PMT was measured at low temperatures down to -110 $^0$C (a typical the PMT operation temperature in liquid xenon detectors) in a spectral range from 154.5 nm to 400 nm. At -110 $^0$C the absolute QE increased by 10-15\% at 175 nm compared to that measured at room temperature. A new low power consumption, low radioactivity voltage divider for the PMTs is being developed. The measurement results showed that the PMT with the current version of the divider demonstrated a linear response (within 5\%) down to 5$\cdot$10$^4$ photoelectrons at a rate of 200 Hz. The radioactive contamination induced by the PMT and the PMT voltage divider materials satisfies the requirements for XENON1T detector not to exceed a total radioactive contamination in the detector of 0.5 evts/year/1tonn. Most of the PMTs received from the manufacturer showed a high quantum efficiency exceeding 30\%. In the mass production tests the measurements at room temperature showed clear single photoelectron peaks for all PMTs been under study. The optimal operation conditions were found at a gain of 2$\cdot$10$^6$. The operation stability for most of the PMTs was also demonstrated at a temperature of -100 $^0$C. A dedicated setup was built for testing the PMTs in liquid xenon using the XENON1T signal readout components including voltage dividers, cables and feedthroughs. The PMTs tested in liquid xenon demonstrated a stable operation for a time period of more than 5 months.
We present examples of non-Gaussian statistics that can induce bispectra matching local and non-local (including equilateral) templates in biased sub-volumes. We find cases where the biasing from coupling to long wavelength modes affects only the power spectrum, only the bispectrum or both. Our results suggest that ruling out multi-field scenarios is quite difficult: some measurements of the scalar correlations, including the shape of the bispectrum, can be consistent with single-clock predictions even when cosmic variance from super-horizon modes is at work. Furthermore, if observations of the density perturbations rule out single-clock inflation, we will face a serious cosmic variance obstacle in drawing any further conclusions about the particle physics origin of the scalar fluctuations.
Planetary Nebulae (PNe) are amongst the most spectacular objects produced by stellar evolution, but the exact identity of their progenitors has never been established for a large and homogeneous observational sample. We investigate the relationship between PNe and their stellar progenitors in the Large Magellanic Cloud (LMC) through the statistical comparison between a highly complete spectroscopic catalog of PNe and the spatially resolved age distribution of the underlying stellar populations. We find that most PN progenitors in the LMC have main-sequence lifetimes in a narrow range between 5 and 8 Gyr, which corresponds to masses between 1.2 and 1.0 M$_{\odot}$, and produce PNe that last $26^{+6}_{-7}$~kyr on average. We tentatively detect a second population of PN progenitors, with main-sequence lifetimes between 35 and 800~Myr, i.e., masses between 8.2 and 2.1 M$_{\odot}$, and average PN lifetimes of $11^{+6}_{-7}$ kyr. These two distinct and disjoint populations of progenitors strongly suggest the existence of at least two physically distinct formation channels for PNe. Our determination of PN lifetimes and progenitor masses has implications for the understanding of PNe in the context of stellar evolution models, and for the role that rotation, magnetic fields, and binarity can play in the shaping of PN morphologies.
In collisional fluids, a number of key processes rely on the frequency of binary collisions. Collisions seem necessary to generate a shock wave when two fluids collide fast enough, to fulfill the Rankine-Hugoniot relations, to establish an equation of state or a Maxwellian distribution. Yet, these seemingly collisional features are routinely either observed or assumed, in relation with collision\emph{less} astrophysical plasmas. This article will review our current answers to the following questions: How do colliding collisionless plasmas end-up generating a shock as if they were fluids? To which extent are the Rankine-Hugoniot relations fulfilled in this case? Do collisionless shocks propagate like fluid ones? Can we use an equation of state to describe collisionless plasmas, like MHD codes for astrophysics do? Why are Maxwellian distributions ubiquitous in Particle-In-Cell simulations of collisionless shocks? Time and length scales defining the border between the collisional and the collisionless behavior will be given when relevant. In general, when the time and length scales involved in the collisionless processes responsible for the fluid-like behavior may be neglected, the system may be treated like a fluid.
This paper investigates the manner in which classical universes are obtained in the no-boundary quantum state. In this context, universes can be characterised as classical (in a WKB sense) when the wavefunction is highly oscillatory, i.e. when the ratio of the change in the amplitude of the wavefunction becomes small compared to the change in the phase. In the presence of a positive or negative exponential potential, the WKB condition is satisfied in proportion to a factor $e^{-(\epsilon - 3)N/(\epsilon -1)},$ where $\epsilon$ is the (constant) slow-roll/fast-roll parameter and $N$ designates the number of e-folds. Thus classicality is reached exponentially fast in $N$, but only when $\epsilon < 1$ (inflation) or $\epsilon > 3$ (ekpyrosis). Furthermore, when the potential switches off and the ekpyrotic phase goes over into a phase of kinetic domination, the level of classicality obtained up to that point is preserved. Similar results are obtained in a cyclic potential, where a dark energy plateau is added. Finally, for a potential of the form $-\phi^n$ (with $n=4,6,8$), where the classical solution becomes increasingly kinetic-dominated, there is an initial burst of classicalisation which then quickly levels off. These results demonstrate that inflation and ekpyrosis, which are the only dynamical mechanisms known for smoothing the universe, share the perhaps even more fundamental property of rendering space and time classical in the first place.
In the early seventies, Alan Sandage defined cosmology as the search for two
numbers: Hubble parameter ${{H}_{0}}$ and deceleration parameter ${{q}_{0}}$.
The first of the two basic cosmological parameters (the Hubble parameter)
describes the linear part of the time dependence of the scale factor. Treating
the Universe as a dynamical system it is natural to assume that it is
non-linear: indeed, linearity is nothing more than approximation, while
non-linearity represents the generic case. It is evident that future models of
the Universe must take into account different aspects of its evolution. As soon
as the scale factor is the only dynamical variable, the quantities which
determine its time dependence must be essentially present in all aspects of the
Universe' evolution. Basic characteristics of the cosmological evolution, both
static and dynamical, can be expressed in terms of the parameters ${{H}_{0}}$
and ${{q}_{0}}$. The very parameters (and higher time derivatives of the scale
factor) enable us to construct model-independent kinematics of the cosmological
expansion.
Time dependence of the scale factor reflects main events in history of the
Universe. Moreover it is the deceleration parameter who dictates the expansion
rate of the Hubble sphere and determines the dynamics of the observable galaxy
number variation: depending on the sign of the deceleration parameter this
number either grows (in the case of decelerated expansion), or we are going to
stay absolutely alone in the cosmos (if the expansion is accelerated).
The intended purpose of the report is reflected in its title --- "Cosmology
in terms of the deceleration parameter". We would like to show that practically
any aspect of the cosmological evolution is tightly bound to the deceleration
parameter.
We present a particle-in-cell simulation of the generation of a collisionless turbulent shock in a dense plasma driven by an ultra-high-intensity laser pulse. From the linear analysis, we highlight the crucial role of the laser-heated and return-current electrons in triggering a strong Weibel-like instability, giving rise to a magnetic turbulence able to isotropize the target ions.
The emergence of the seeds of cosmic structure from a perfect isotropic and homogeneous Universe has not been fully explained by the standard version of inflationary models. To handle this shortcoming, D. Sudarsky and collaborators have developed a proposal: "the self-induced collapse hypothesis." In this scheme, the collapse of the inflaton wave function is responsible for the emergence of inhomogeneity and anisotropy at each scale. In previous papers, the proposal was developed with an almost exact de Sitter space-time approximation for the background that lead to a perfect scale-invariant power spectrum. In this paper, we consider a quasi-de Sitter expansion factor and calculate the primordial power spectrum for three different choices of the self-induced collapse. The consideration of a quasi-de Sitter background allow us to distinguish departures from an exact scale-invariant power spectrum that are due to the inclusion of the collapse hypothesis. These deviations are also different from the prediction of standard inflationary models with running spectral index. Comparison with the primordial power spectrum preferred by the latest observational data is also discussed. From the analysis performed in this paper, it follows that the collapse schemes analysed in this paper are viable candidates to explain present observations of the CMB fluctuation spectrum.
We consider the causal structure of generalized uncharged McVittie spacetimes with increasing central mass $m (t)$ and positive Hubble factor $H (t)$. Under physically reasonable conditions, namely, a big bang singularity in the past, a positive cosmological constant and an upper limit to the central mass, we prove that the patch of the spacetime described by the cosmological time and areal radius coordinates is always geodesically incomplete, which implies the presence of event horizons in the spacetime. We also show that, depending on the asymptotic behavior of the $m$ and $H$ functions, the generalized McVittie spacetime can have a single black hole, a black-hole/white-hole pair or, differently from classic fixed-mass McVittie, a single white hole. A simple criterion is given to distinguish the different causal structures.
We perform a detailed and illustrative study of the production of keV sterile neutrino Dark Matter (DM) by decays of singlet scalars in the early Universe. In the current study we focus on providing a clear and general overview of this production mechanism. For the first time we study all regimes possible on the level of momentum distribution functions, which we obtain by solving a system of Boltzmann equations. These quantities contain the full information about the production process, which allows us to not only track the evolution of the DM generation but to also take into account all bounds related to the spectrum, such as constraints from structure formation or from avoiding too much dark radiation. In particular we show that this simple production mechanism can, depending on the regime, lead to strongly non-thermal DM spectra which may even feature more than one peak in the momentum distribution. These cases could have particularly interesting consequences for cosmological structure formation, as their analysis requires more refined tools than the simplistic estimate using the free-streaming horizon. Here we present the mechanism including all concepts and subtleties involved, for now using the assumption that the effective number of relativistic degrees of freedom is constant during DM production, which is applicable in a significant fraction of the parameter space. This allows us to derive analytical results to back up our detailed numerical computations, thus leading to the most comprehensive picture of keV sterile neutrino DM production by singlet scalar decays that exists up to now.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
We present a high spatial resolution ($\approx 20$ pc) of $^{12}$CO($2-1$) observations of the lenticular galaxy NGC4526. We identify 103 resolved Giant Molecular Clouds (GMCs) and measure their properties: size $R$, velocity dispersion $\sigma_v$, and luminosity $L$. This is the first GMC catalog of an early-type galaxy. We find that the GMC population in NGC4526 is gravitationally bound, with a virial parameter $\alpha \sim 1$. The mass distribution, $dN/dM \propto M^{-2.39 \pm 0.03}$, is steeper than that for GMCs in the inner Milky Way, but comparable to that found in some late-type galaxies. We find no size-linewidth correlation for the NGC4526 clouds, in contradiction to the expectation from Larson's relation. In general, the GMCs in NGC4526 are more luminous, denser, and have a higher velocity dispersion than equal size GMCs in the Milky Way and other galaxies in the Local Group. These may be due to higher interstellar radiation field than in the Milky Way disk and weaker external pressure than in the Galactic center. In addition, a kinematic measurement of cloud rotation shows that the rotation is driven by the galactic shear. For the vast majority of the clouds, the rotational energy is less than the turbulent and gravitational energy, while the four innermost clouds are unbound and will likely be torn apart by the strong shear at the galactic center. We combine our data with the archival data of other galaxies to show that the surface density $\Sigma$ of GMCs is not approximately constant as previously believed, but varies by $\sim 3$ orders of magnitude. We also show that the size and velocity dispersion of GMC population across galaxies are related to the surface density, as expected from the gravitational and pressure equilibrium, i.e. $\sigma_v R^{-1/2} \propto \Sigma^{1/2}$.
We present multi-epoch optical spectroscopy of seven southern Fermi-monitored blazars from 2008 - 2013 using the Small and Medium Aperture Research Telescope System (SMARTS), with supplemental spectroscopy and polarization data from the Steward Observatory. We find that the emission lines are much less variable than the continuum; 4 of 7 blazars had no detectable emission line variability over the 5 years. This is consistent with photoionization primarily by an accretion disk, allowing us to use the lines as a probe of disk activity. Comparing optical emission line flux with Fermi $\gamma$-ray flux and optical polarized flux, we investigate whether relativistic jet variability is related to the accretion flow. In general, we see no such dependence, suggesting the jet variability is likely caused by internal processes like turbulence or shock acceleration rather than a variable accretion rate. However, three sources showed statistically significant emission line flares in close temporal proximity to very large Fermi $\gamma$-ray flares. While we do not have sufficient emission line data to quantitatively assess their correlation with the $\gamma$-ray flux, it appears that in some cases, the jet might provide additional photoionizing flux to the broad line region, which implies some gamma-rays are produced within the broad line region, at least for these large flares.
We carry out direct numerical simulations of turbulent astrophysical media that explicitly track ionizations, recombinations, and species-by-species radiative cooling. The simulations assume solar composition and follows the evolution of hydrogen, helium, carbon, oxygen, sodium, and magnesium, but they do not include the presence of an ionizing background. In this case, the medium reaches a global steady state that is purely a function of the one-dimensional turbulent velocity dispersion, $\sigma_{\rm 1D},$ and the product of the mean density and the driving scale of turbulence, $n L.$ Our simulations span a grid of models with $\sigma_{\rm 1D}$ ranging from 6 to 58 km s$^{-1}$ and $n L$ ranging from 10$^{16}$ to 10$^{20}$ cm$^{-2},$ which correspond to turbulent Mach numbers from $M=0.2$ to 10.6. The species abundances are well described by single-temperature estimates whenever $M$ is small, but local equilibrium models can not accurately predict the global equilibrium abundances when $M \gtrsim 1.$ To allow future studies to account for nonequilibrium effects in turbulent media, we gather our results into a series of tables, which we will extend in the future to encompass a wider range of elements, compositions, and ionizing processes.
At a distance of 50 kpc and with a dark matter mass of $\sim10^{10}$ M$_{\odot}$, the Large Magellanic Cloud (LMC) is a natural target for indirect dark matter searches. We use five years of data from the Fermi Large Area Telescope (LAT) and updated models of the gamma-ray emission from standard astrophysical components to search for a dark matter annihilation signal from the LMC. We perform a rotation curve analysis to determine the dark matter distribution, setting a robust minimum on the amount of dark matter in the LMC, which we use to set conservative bounds on the annihilation cross section. The LMC emission is generally very well described by the standard astrophysical sources, with at most a $1-2\sigma$ excess identified near the kinematic center of the LMC once systematic uncertainties are taken into account. We place competitive bounds on the dark matter annihilation cross section as a function of dark matter particle mass and annihilation channel.
Observations of gamma-ray-bursts and jets from active galactic nuclei reveal that the jet flow is characterized by a high radiative efficiency and that the dissipative mechanism must be a powerful accelerator of non-thermal particles. Shocks and magnetic reconnection have long been considered as possible candidates for powering the jet emission. Recent progress via fully-kinetic particle-in-cell simulations allows us to revisit this issue on firm physical grounds. We show that shock models are unlikely to account for the jet emission. In fact, when shocks are efficient at dissipating energy, they typically do not accelerate particles far beyond the thermal energy, and vice versa. In contrast, we show that magnetic reconnection can deposit more than 50% of the dissipated energy into non-thermal leptons as long as the energy density of the magnetic field in the bulk flow is larger than the rest mass energy density. The emitting region, i.e., the reconnection downstream, is characterized by a rough energy equipartition between magnetic fields and radiating particles, which naturally accounts for a commonly observed property of blazar jets.
The radio frequency 1.4 GHz transition of the atomic hydrogen is one of the important tracers of the diffuse neutral interstellar medium. Radio astronomical observations of this transition, using either a single dish telescope or an array interferometer, reveal different properties of the interstellar medium. Such observations are particularly useful to study the multiphase nature and turbulence in the interstellar gas. Observations with multiple radio telescopes have recently been used to study these two closely related aspects in greater detail. Using various observational techniques, the density and the velocity fluctuations in the Galactic interstellar medium was found to have a Kolmogorov-like power law power spectra. The observed power law scaling of the turbulent velocity dispersion with the length scale can be used to derive the true temperature distribution of the medium. Observations from a large ongoing atomic hydrogen absorption line survey have also been used to study the distribution of gas at different temperature. The thermal steady state model predicts that the multiphase neutral gas will exist in cold and warm phase with temperature below 200 K and above 5000 K respectively. However, these observations clearly show the presence of a large fraction of gas in the intermediate unstable phase. These results raise serious doubt about the validity of the standard model, and highlight the necessity of alternative theoretical models. Interestingly, numerical simulations suggest that some of the observational results can be explained consistently by including the effects of turbulence in the models of the multiphase medium. This review article presents a brief outline of some of the basic ideas of radio astronomical observations and data analysis, summarizes the results from recent observations, and discusses possible implications of the results.
One of the most fundamental features of a black hole in general relativity is its event horizon: a boundary from which nothing can escape. There has been a recent surge of interest in the nature of these event horizons and their local neighbourhoods. In an attempt to resolve black hole information paradox(es), and more generally, to better understand the path towards quantum gravity, firewalls have been proposed as an alternative to black hole event horizons. In this letter, we explore the phenomenological implications of black holes possessing a surface or firewall. We predict a potentially detectable signature of these firewalls in the form of a high energy astrophysical neutrino flux. We compute the spectrum of this neutrino flux in different models and show that it is a possible candidate for the source of the PeV neutrinos recently detected by IceCube. We further show that, independent of the generation mechanism, IceCube data can be explained (at $1\sigma$ confidence level) by conversion of accretion onto stellar mass (supermassive) black holes, into neutrinos at efficiencies of $\gtrsim 10^{-4} (10^{-2})$.
The use of large, X-ray selected galaxy cluster catalogs for cosmological analyses requires a thorough understanding of the X-ray mass estimates, including the possibility of biases due to the assumption of hydrostatic equilibrium. Weak gravitational lensing is an ideal method to shed light on such issues, due to its insensitivity to the cluster dynamical state. We perform a weak lensing calibration of 166 galaxy clusters from the RBC cluster catalog and compare our results to the hydrostatic X-ray masses from that catalog. To interpret the weak lensing signal in terms of cluster masses, we compare the lensing signal to simple theoretical Navarro-Frenk-White models and to simulated cluster lensing profiles, including complications such as cluster substructure, projected large-scale structure, and Malmquist bias. We find evidence of underestimation in the X-ray masses, as expected, with $<M_{\mathrm{X}}/M_{\mathrm{WL}}>= 0.66_{-0.12}^{+0.07}$ for our best-fit model, a more than $4\sigma$ detection of a bias between X-ray and weak lensing masses. The biases in cosmological parameters in a typical cluster abundance measurement that ignores this mass bias will typically exceed the statistical errors.
We measure the location and evolutionary vectors of 69 Herschel-detected broad-line active galactic nuclei (BLAGNs) in the M_BH-M_* plane. BLAGNs are selected from the COSMOS and CDF-S fields, and span the redshift range 0.2< z<2.1. Black-hole masses are calculated using archival spectroscopy and single-epoch virial mass estimators, and galaxy total stellar masses are calculated by fitting the spectral energy distribution (subtracting the BLAGN component). The mass-growth rates of both the black hole and galaxy are calculated using Chandra/XMM-Newton X-ray and Herschel far-infrared data, reliable measures of the BLAGN accretion and galaxy star formation rates, respectively. We use Monte Carlo simulations to account for biases in our sample, due to both selection limits and the steep slope of the massive end of the galaxy stellar-mass distribution. We find our sample is consistent with no evolution in the M_BH-M_* relation from z~2 to z~0. BLAGNs and their host galaxies which lie off the black hole mass galaxy total stellar mass relation tend to have evolutionary vectors anti-correlated with their mass ratios: that is, galaxies with over-massive (under-massive) black holes tend to have a low (high) ratio of the specific accretion rate to the specific star formation rate. We also use the measured growth rates to estimate the preferred AGN duty cycle for our galaxies to evolve most consistently onto the local M_BH-M_Bul relation. Under reasonable assumptions of exponentially declining star formation histories, the data suggest a non-evolving (no more than a factor of a few) BLAGN duty cycle among star-forming galaxies of 10% (1sigma range of 1-42% at z<1 and 2-34% at z>1).
NICMOS 2 observations are crucial for constraining distances to most of the existing sample of z > 1 SNe Ia. Unlike the conventional calibration programs, these observations involve long exposure times and low count rates. Reciprocity failure is known to exist in HgCdTe devices and a correction for this effect has already been implemented for high and medium count-rates. However observations at faint count-rates rely on extrapolations. Here instead, we provide a new zeropoint calibration directly applicable to faint sources. This is obtained via inter-calibration of NIC2 F110W/F160W with WFC3 in the low count-rate regime using z ~ 1 elliptical galaxies as tertiary calibrators. These objects have relatively simple near-IR SEDs, uniform colors, and their extended nature gives superior signal-to-noise at the same count rate than would stars. The use of extended objects also allows greater tolerances on PSF profiles. We find ST magnitude zeropoints (after the installation of the NICMOS cooling system, NCS) of 25.296 +- 0.022 for F110W and 25.803 +- 0.023 for F160W, both in agreement with the calibration extrapolated from count-rates 1,000 times larger (25.262 and 25.799). Before the installation of the NCS, we find 24.843 +- 0.025 for F110W and 25.498 +- 0.021 for F160W, also in agreement with the high-count-rate calibration (24.815 and 25.470). We also check the standard bandpasses of WFC3 and NICMOS 2 using a range of stars and galaxies at different colors and find mild tension for WFC3, limiting the accuracy of the zeropoints. To avoid human bias, our cross-calibration was "blinded" in that the fitted zeropoint differences were hidden until the analysis was finalized.
The Main Sequence (MS) of star-forming galaxies plays a fundamental role in driving galaxy evolution and in our efforts to understand it. However, different studies find significant differences in the normalization, slope and shape of the MS. These discrepancies arise mainly from the different selection criteria adopted to isolate star-forming galaxies, that may include or exclude galaxies with specific star formation rate (SFR) substantially below the MS value. To obviate this limitation of all current criteria, we propose an objective definition of the MS that does not rely at all on a pre-selection of star-forming galaxies. Constructing the 3D SFR-Mass-Number plot, the MS is then defined as the ridge line of the star-forming peak, as illustrated with various figures. The advantages of such definition are manifold. If generally adopted it will facilitate the inter-comparison of results from different groups using the same star formation rate (SFR) and stellar mass diagnostics, or to highlight the relative systematics of different diagnostics. All this could help understanding MS galaxies as systems in a quasi-steady state equilibrium and would also provide a more objective criterion for identifying quenching galaxies.
We study the influence of the presence of a strong bar in disc galaxies which host an active galactic nucleus (AGN). Using data from the Sloan Digital Sky Survey and morphological classifications from the Galaxy Zoo 2 project, we create a volume-limited sample of 19,756 disc galaxies at $0.01<z<0.05$ which have been visually examined for the presence of a bar. Within this sample, AGN host galaxies have a higher overall percentage of bars (51.8%) than inactive galaxies exhibiting central star formation (37.1%). This difference is primarily due to known effects; that the presence of both AGN and galactic bars is strongly correlated with both the stellar mass and integrated colour of the host galaxy. We control for this effect by examining the difference in AGN fraction between barred and unbarred galaxies in fixed bins of mass and colour. Once this effect is accounted for, there remains a small but statistically significant increase that represents 16% of the average barred AGN fraction. Using the $L_{\rm{[O III]}}/M_{BH} $ratio as a measure of AGN strength, we show that barred AGN do not exhibit stronger accretion than unbarred AGN at a fixed mass and colour. The data are consistent with a model in which bar-driven fueling does contribute to the probability of an actively growing black hole, but in which other dynamical mechanisms must contribute to the direct AGN fueling via smaller, non-axisymmetric perturbations.
We present host stellar velocity dispersion measurements for a sample of 88 broad-line quasars at 0.1<z<1 (46 at z>0.6) from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. High signal-to-noise ratio coadded spectra (average S/N~30 per 69 km/s pixel) from SDSS-RM allowed decomposition of the host and quasar spectra, and measurement of the host stellar velocity dispersions and black hole (BH) masses using the single-epoch (SE) virial method. The large sample size and dynamic range in luminosity (L5100=10^(43.2-44.7) erg/s) lead to the first clear detection of a correlation between SE virial BH mass and host stellar velocity dispersion far beyond the local universe. However, the observed correlation is significantly flatter than the local relation, suggesting that there are selection biases in high-z luminosity-threshold quasar samples for such studies. Our uniform sample and analysis enable an investigation of the redshift evolution of the M-sigma relation free of caveats by comparing different samples/analyses at disjoint redshifts. We do not observe evolution of the M-sigma relation in our sample, up to z~1, but there is an indication that the relation flattens towards higher redshifts. Coupled with the increasing threshold luminosity with redshift in our sample, this again suggests certain selection biases are at work, and simple simulations demonstrate that a constant M-sigma relation is favored to z~1. Our results highlight the scientific potential of deep coadded spectroscopy from quasar monitoring programs, and offer a new path to probe the co-evolution of BHs and galaxies at earlier times.
With Australia Telescope Compact Array observations, we detect a highly elongated Mpc-scale diffuse radio source on the eastern periphery of the Bullet cluster 1E0657-55.8, which we argue has the positional, spectral and polarimetric characteristics of a radio relic. This powerful relic (2.3+/-0.1 x 10^25 W Hz^-1) consists of a bright northern bulb and a faint linear tail. The bulb emits 94% of the observed radio flux and has the highest surface brightness of any known relic. Exactly coincident with the linear tail we find a sharp X-ray surface brightness edge in the deep Chandra image of the cluster -- a signature of a shock front in the hot intracluster medium (ICM), located on the opposite side of the cluster to the famous bow shock. This new example of an X-ray shock coincident with a relic further supports the hypothesis that shocks in the outer regions of clusters can form relics via diffusive shock (re-)acceleration. Intriguingly, our new relic suggests that seed electrons for reacceleration are coming from a local remnant of a radio galaxy, which we are lucky to catch before its complete disruption. If this scenario, in which a relic forms when a shock crosses a well-defined region of the ICM polluted with aged relativistic plasma -- as opposed to the usual assumption that seeds are uniformly mixed in the ICM -- is also the case for other relics, this may explain a number of peculiar properties of peripheral relics.
Organic matter and hydrous silicates are intimately mixed in the matrix of chondrites and in-situ determination of their individual D/H ratios is therefore challenging. Nevertheless, the D/H ratio of each pure component in this mixture should yield a comprehensible signature of the origin and evolution of water and organic matter in our solar system. We measured hydrogen isotope ratios of organic and hydrous silicates in the matrices of two carbonaceous chondrites (Orgueil CI1 and Renazzo CR2) and one unequilibrated ordinary chondrite (Semarkona, LL3.0). A novel protocol was adopted, involving NanoSIMS imaging of H isotopes of monoatomatic ($H^-$) and molecular ($OH^-$) secondary ions collected at the same location. This allowed the most enriched component with respect to D to be identified in the mixture. Using this protocol, we found that in carbonaceous chondrites the isotopically homogeneous hydrous silicates are mixed with D-rich organic matter. The opposite was observed in Semarkona. Hydrous silicates in Semarkona display highly heterogeneous D/H ratios, ranging from $150$ to $1800$ ${\times}$ $10^{-6}$ (${\delta}D_{SMOW} = -40$ to $10,600$ permil). Organic matter in Semarkona does not show such large isotopic variations. This suggests limited isotopic exchange between the two phases during aqueous alteration. Our study greatly expands the range of water isotopic values measured so far in solar system objects. This D-rich water reservoir was sampled by the LL ordinary chondrite parent body and an estimate (up to 9 %) of its relative contribution to the D/H ratio of water in Oort cloud family comets is proposed.
Prediction of the shocks' arrival times (SATs) at the Earth is very important for space weather forecast. There is a well-known SAT model, STOA, which is widely used in the space weather forecast. However, the shock transit time from STOA model usually has a relative large error compared to the real measurements. In addition, STOA tends to yield too much `yes' prediction, which causes a large number of false alarms. Therefore, in this work, we work on the modification of STOA model. First, we give a new method to calculate the shock transit time by modifying the way to use the solar wind speed in STOA model. Second, we develop new criteria for deciding whether the shock will arrive at the Earth with the help of the sunspot numbers and the angle distances of the flare events. It is shown that our work can improve the SATs prediction significantly, especially the prediction of flare events without shocks arriving at the Earth.
Monitor of All sky X-ray Image (MAXI) discovered a new outburst of an X-ray transient source named MAXI J1421-613. Because of the detection of three X-ray bursts from the source, it was identified as a neutron star low-mass X-ray binary. The results of data analyses of the MAXI GSC and the Swift XRT follow-up observations suggest that the spectral hardness remained unchanged during the first two weeks of the outburst. All the XRT spectra in the 0.5-10 keV band can be well explained by thermal Comptonization of multi-color disk blackbody emission. The photon index of the Comptonized component is $\approx$ 2, which is typical of low-mass X-ray binaries in the low/hard state. Since X-ray bursts have a maximum peak luminosity, it is possible to estimate the (maximum) distance from its observed peak flux. The peak flux of the second X-ray burst, which was observed by the GSC, is about 5 photons cm$^{-2}$ s$^{-1}$. By assuming a blackbody spectrum of 2.5 keV, the maximum distance to the source is estimated as 7 kpc. The position of this source is contained by the large error regions of two bright X-ray sources detected with Orbiting Solar Observatory-7 (OSO-7) in 1970s. Besides this, no past activities at the XRT position are reported in the literature. If MAXI J1421-613 is the same source as (one of) them, the outburst observed with MAXI may have occurred after the quiescence of 30-40 years.
In order to better understand the nature of post-starburst quasars (PSQs) in the context of galaxy evolution, we compare their properties to those of post-starburst galaxies and quasars from appropriately selected samples possessing similar redshift ($z \sim 0.3$), luminosity ($M_r \sim -$23), and data quality. We consider morphologies, spectral features, and derived physical properties of the stellar populations and central supermassive black hole. PSQs themselves come in two types: the more luminous AGNs with more luminous post-starburst stellar populations hosted by elliptical galaxies, some which are clearly merger products, and the less luminous systems existing within relatively undisturbed spiral galaxies and possessing signs of a more extended period of star formation. Post-starburst galaxies (PSQs) have elliptical and disturbed/post-merger morphologies similar to those of the more luminous PSQs, display similar spectral properties, but also can have younger stellar populations for a given starburst mass. Quasars at similar redshifts and luminosities around the Seyfert/quasar transition possess similar AGN characteristics, in terms of black hole mass and accretion rate, compared with those of PSQs, but do not appear to be hosted by galaxies with significant post-starburst populations. Recent studies of more luminous quasars find hosts consistent with those of our luminous PSQs, suggesting that these PSQs may be in transition between post-starburst galaxies and a more luminous quasar stage when the post-starburst stellar population remains dominant. The lower luminosity PSQs appear to differ from lower luminosity quasars (Seyfert galaxies) in terms of more significant star formation in their past.
The UV emission from X-ray binaries, is more likely to be produced by reprocessing of X-rays by the outer regions of an accretion disk. The structure of the outer disk may be altered due to the presence of X-ray irradiation and we discuss the physical regimes where this may occur and point out certain X-ray binaries where this effect may be important. The long term X-ray variability of these sources is believed to be due to stochastic fluctuations in the outer disk, which propagate inwards giving rise to accretion rate variation in the X-ray producing inner regions. The X-ray variability will induce structural variations in the outer disk which in turn may affect the inner accretion rate. To understand the qualitative behaviour of the disk in such a scenario, we adopt simplistic assumptions that the disk is fully ionised and is not warped. We develop and use a time dependent global hydrodynamical code to study the effect of a sinusoidal accretion rate perturbation introduced at a specific radius. The response of the disk, especially the inner accretion rate, to such perturbations at different radii and with different time periods is shown. While we didn't find any oscillatory or limit cycle behaviour, our results show that irradiation enhances the X-ray variability at time-scales corresponding to the viscous time-scales of the irradiated disk.
Primordial black holes (PBHs) are theoretical black holes which may be formed during the radiation dominant era and, basically, caused by the gravitational collapse of radiational overdensities. It has been well known that in the context of the structure formation in our Universe such collapsed objects, e.g., halos/galaxies, could be considered as bias tracers of underlying matter fluctuations and the halo/galaxy bias has been studied well. Employing a peak-background split picture which is known to be a useful tool to discuss the halo bias, we consider the large scale clustering behavior of the PBH and propose an almost mass-independent constraint to the scenario that dark matters (DMs) consist of PBHs. We consider the case where the statistics of the primordial curvature perturbations is almost Gaussian, but with small local-type non-Gaussianity. If PBHs account for the DM abundance, such a large scale clustering of PBHs behaves as nothing but the matter isocurvature perturbation and constrained strictly by the observations of cosmic microwave backgrounds (CMB). From this constraint, we show that, in the case a certain single field causes both CMB temperature perturbations and PBH formations, the PBH-DM scenario is excluded even with quite small local-type non-Gaussianity, $|f_\mathrm{NL}|\sim\mathcal{O}(0.01)$, while we give the constraints to parameters in the case where the source field of PBHs is different from CMB perturbations.
IC 310 has recently been identified as a gamma-ray emitter based on observations at GeV energies with Fermi-LAT and at very high energies (VHE, E > 100 GeV) with the MAGIC telescopes. Despite IC 310 having been classified as a radio galaxy with the jet observed at an angle > 10 degrees, it exhibits a mixture of multiwavelength properties of a radio galaxy and a blazar, possibly making it a transitional object. On the night of 12/13th of November 2012 the MAGIC telescopes observed a series of violent outbursts from the direction of IC 310 with flux-doubling time scales faster than 5 min and a peculiar spectrum spreading over 2 orders of magnitude. Such fast variability constrains the size of the emission region to be smaller than 20% of the gravitational radius of its central black hole, challenging the shock acceleration models, commonly used in explanation of gamma-ray radiation from active galaxies. Here we will show that this emission can be associated with pulsar-like particle acceleration by the electric field across a magnetospheric gap at the base of the jet.
We perform realistic simulations of the Sun's birth cluster in order to predict the current distribution of solar siblings in the Galaxy. We study the possibility of finding the solar siblings in the Gaia catalogue by using only positional and kinematic information. We find that the number of solar siblings predicted to be observed by Gaia will be around 100 in the most optimistic case, and that a phase space only search in the Gaia catalogue will be extremely difficult. It is therefore mandatory to combine the chemical tagging technique with phase space selection criteria in order to have any hope of finding the solar siblings.
In this paper, we constrain the dimensionless Compton wavelength parameter $B_0$ of $f(R)$ gravity as well as the mass of sterile neutrino by using the cosmic microwave background observations, the baryon acoustic oscillation surveys, and the linear growth rate measurements. Since both the $f(R)$ model and the sterile neutrino generally predict scale-dependent growth rates, we utilize the growth rate data measured in different wavenumber bins with the theoretical growth rate approximatively scale-independent in each bin. The employed growth rate data come from the peculiar velocity measurements at $z=0$ in five wavenumber bins, and the redshift space distortions measurements at $z=0.25$ and $z=0.37$ in one wavenumber bin. By constraining the $f(R)$ model alone, we get a tight 95% upper bound of $\log_{10}B_0<-4.1$. This result is slightly weakened to $\log_{10}B_0<-3.8$ (at 2$\sigma$ level) once we simultaneously constrain the $f(R)$ model and the sterile neutrino mass, due to the degeneracy between the parameters of the two. For the massive sterile neutrino parameters, we get the effective sterile neutrino mass $m_{\nu,{\rm{sterile}}}^{\rm{eff}}<0.62$ eV (2$\sigma$) and the effective number of relativistic species $N_{\rm eff}=3.40^{+0.10}_{-0.35}$ in the $f(R)$ model. As a comparison, we also obtain $m_{\nu,{\rm{sterile}}}^{\rm{eff}}<0.56$ eV (2$\sigma$) and $N_{\rm eff}=3.43^{+0.10}_{-0.38}$ in the standard $\Lambda$CDM model.
AIMS: The aim of this study is to investigate the origin and evolution of Sc, V, Mn, and Co for a homogeneous and statistically significant sample of stars probing the different populations of the Milky Way, in particular the thin and thick disks. METHODS: Using high-resolution spectra obtained with MIKE, FEROS, SOFIN, FIES, UVES and HARPS spectrographs, we determine Sc, V, Mn, and Co abundances for a large sample of F and G dwarfs in Solar neighbourhood. The method is based on spectral synthesis and using one-dimensional (1-D), plane-parallel, LTE model stellar atmospheres calculated with the MARCS 2012 code. NLTE corrections from literature were applied to Mn and Co. RESULTS: We find that the abundance trends derived for Sc (594 stars), V (466 stars), and Co (567 stars) are very similar to what has been observed for the alpha-elements in the thin and thick disks. On the contrary Mn (569 stars) is generally underabundant relative to the Sun (i.e. [Mn/Fe]<0) for [Fe/H]<0. Also, for Mn, when NLTE corrections are applied, the trend changes and is almost flat over the entire metallicity range that the stars in our sample span (-2<[Fe/H]<+0.4). The [Sc/Fe]-[Fe/H] abundance trends show a small separation between the thin and thick disks, while for V and Co they completely overlap. For Mn there is a small difference in [Mn/Fe] but only when NLTE corrections are used. Comparisons with Ti as a reference element show flat trends for all the elements except for Mn that show well separated [Mn/Ti]-[Ti/H] trends for the thin and thick disks. CONCLUSIONS: Sc and V present trends compatible with production from SNII events. In addition, Sc clearly shows a metallicity dependence for [Fe/H]<-1. Mn instead is produced in SNII events for [Fe/H]<-0.4 and then SNIa start to produce Mn. Finally, Co appears to be produced mainly in SNII with suggestion of enrichment from hypernovae at low metallicities.
The form of the solar meridional circulation is a very important ingredient for mean field flux transport dynamo models. Yet a shroud of mystery still surrounds this large-scale flow, given that its measurement using current helioseismic techniques is challenging. In this work we use results from 3D global simulations of solar convection to infer the dynamical behavior of the established meridional circulation. We make a direct comparison between the meridional circulation that arises in these simulations and the latest observations. Based on our results we argue that there should be an equatorward flow at the base of the convection zone at mid latitudes, below the current maximum depth helioseismic measures can probe (0.75 R). We also provide physical arguments to justify this behaviour. The simulations indicate that the meridional circulation undergoes substantial changes in morphology as the magnetic cycle unfolds. We close by discussing the importance of these dynamical changes for current methods of observation that involve long averaging periods of helioseismic data. Also noteworthy is the fact that these topological changes indicate a rich interaction between magnetic fields and plasma flows, which challenges the ubiquitous kinematic approach used in the vast majority of mean field dynamo simulations.
Inspired by the multiverse scenario, we study a heterotic dark energy model in which there are two parts, the first being the cosmological constant and the second being the holographic dark energy, thus this model is named the $\Lambda$HDE model. By studying the $\Lambda$HDE model theoretically, we find that the parameters $c$ and $\Omega_{hde}$ are divided into a few domains in which the fate of the universe is quite different. We investigate dynamical behaviors of this model, and especially the future evolution of the universe. We perform fitting analysis on the cosmological parameters in the $\Lambda$HDE model by using the recent observational data. We find the model yields $\chi^2_{\rm min}=426.27$ when constrained by Planck+SNLS3+BAO+HST, comparable to the results of the HDE model (428.20) and the concordant $\Lambda$CDM model (431.35). At 68.3\% CL, we obtain $-0.07<\Omega_{\Lambda0}<0.68$ and correspondingly $0.04<\Omega_{hde0}<0.79$, implying at present there is considerable degeneracy between the holographic dark energy and cosmological constant components in the $\Lambda$HDE model.
Context. Within the core accretion scenario of planetary formation, most
simulations performed so far always assume the accreting envelope to have a
solar composition. From the study of meteorite showers on Earth and numerical
simulations, we know that planetesimals must undergo thermal ablation and
disruption when crossing a protoplanetary envelope. Once the protoplanet has
acquired an atmosphere, the primordial envelope gets enriched in volatiles and
silicates from the planetesimals. This change of envelope composition during
the formation can have a significant effect in the final atmospheric
composition and on the formation timescale of giant planets.
Aims. To investigate the physical implications of considering the envelope
enrichment of protoplanets due to the disruption of icy planetesimals during
their way to the core. Particular focus is placed on the effect on the critical
core mass for envelopes where condensation of water can occur.
Methods. Internal structure models are numerically solved with the
implementation of updated opacities for all ranges of metallicities and the
software CEA to compute the equation of state. CEA computes the chemical
equilibrium for an arbitrary mixture of gases and allows the condensation of
some species, including water. This means that the latent heat of phase
transitions is consistently incorporated in the total energy budget.
Results. The critical core mass is found to decrease significantly when an
enriched envelope composition is considered in the internal structure
equations. A particular strong reduction of the critical core mass is obtained
for planets whose envelope metallicity is larger than Z=0.45 when the outer
boundary conditions are suitable for condensation of water to occur in the top
layers of the atmosphere. We show that this effect is qualitatively preserved
when the atmosphere is out of chemical equilibrium.
We demonstrate the eclipsing binary detection performance of the Gaia variability analysis and processing pipeline using Hipparcos data. The automated pipeline classifies 1,067 (0.9%) of the 118,204 Hipparcos sources as eclipsing binary candidates. The detection rate amounts to 89% (732 sources) in a subset of 819 visually confirmed eclipsing binaries, with the period correctly identified for 80% of them, and double or half periods obtained in 6% of the cases.
Cherenkov telescope experiments, such as H.E.S.S., have been very successful in astronomical observations in the very-high-energy (VHE; E $>$ 100 GeV) regime. As an integral part of the detector, such experiments use Earth's atmosphere as a calorimeter. For the calibration and energy determination, a standard model atmosphere is assumed. Deviations of the real atmosphere from the model may therefore lead to an energy misreconstruction of primary gamma rays. To guarantee satisfactory data quality with respect to difficult atmospheric conditions, several atmospheric data quality criteria are implemented in the H.E.S.S. software. These quantities are sensitive to clouds and aerosols. Here, the Cherenkov transparency coefficient will be presented. It is a new monitoring quantity that is able to measure long-term changes in the atmospheric transparency. The Cherenkov transparency coefficient derives exclusively from Cherenkov data and is quite hardware-independent. Furthermore, its positive correlation with independent satellite measurements, performed by the Multi-angle Imaging SpectroRadiometer (MISR), will be presented.
The Virtual Observatory Registry is a distributed directory of information systems and other resources relevant to astronomy. To make it useful, facilities to query that directory must be provided to humans and machines alike. This article reviews the development and status of such facilities, also considering the lessons learnt from about a decade of experience with Registry interfaces. After a brief outline of the history of the standards development, it describes the use of Registry interfaces in some popular clients as well as dedicated UIs for interrogating the Registry. It continues with a thorough discussion of the design of the two most recent Registry interface standards, RegTAP on the one hand and a full-text-based interface on the other hand. The article finally lays out some of the less obvious conventions that emerged in the interaction between providers of registry records and Registry users as well as remaining challenges and current developments.
We suggest how we can use the mass profile of galaxy clusters beyond their virial radius to measure their mass accretion rate, a key prediction of structure formation models. The mass profile can be estimated by applying the caustic technique to dense redshift surveys of clusters and their outskirts, where dynamical equilibrium does not necessarily hold. An additional probe of the mass growth of clusters is their mass fraction in substructures. We show that the caustic technique, that identifies cluster substructures as a by-product, returns catalogs of substructures with mass larger than a few $10^{13}h^{-1}M_\odot$ that are between 60% and 80% complete, depending on the density of the redshift survey.
We study the outbursts of the black hole X-ray binaries MAXI J1659-152, SWIFT J1753.5--0127 and GX 339-4 with the Swift X-ray Telescope. The bandpass of the X-ray Telescope has access to emission from both components of the accretion flow: the accretion disk and the corona/hot flow. This allows a correlated spectral and variability study, with variability from both components of the accretion flow. We present for the first time, a combined study of the evolution of spectral parameters (disk temperature and radius) and timing parameters (frequency and strength) of all power spectral components in different spectral states. Comparison of the correlations in different spectral states shows that the frequency and strength of the power spectral components exhibit dependencies on the disk temperature that are different in the (low-)hard and the hard-intermediate states; most of these correlations that are clearly observed in the hard-intermediate state (in MAXI J1659-152 and GX 339-4) are not seen in the (low-)hard state (in GX 339-4 and SWIFT J1753.5-0127). Also, the responses of the individual frequency components to changes in the disk temperature are markedly different from one component to the next. Hence, the spectral-timing evolution cannot be explained by a single correlation that spans both these spectral states. We discuss our findings in the context of the existing models proposed to explain the origin of variability.
We present a proper-motion study on models of the dwarf spheroidal galaxy Sculptor, based on the predicted proper-motion accuracy of Gaia measurements. Gaia will measure proper motions of several hundreds of stars for a Sculptor-like system. Even with an uncertainty on the proper motion of order 1.5 times the size of an individual proper-motion value of ~10 mas/century, we find that it is possible to recover Sculptor's systemic proper motion at its distance of 79 kpc.
We study the dependence of the properties of group galaxies on the surrounding large-scale environment, using SDSS-DR7 data. Galaxies are ranked according to their luminosity within each group and classified morphologically by the S\'ersic index. We have considered samples of the host groups in superstructures of galaxies, and elsewhere. We find a significant dependence of the properties of late-type brightest group galaxies on the large-scale environment: they show statistically significant higher luminosities and stellar masses, redder u-r colours, lower star formation activity and longer star-formation time-scale when embedded in superstructures. By contrast, the properties of the early-type brightest group galaxies are remarkably similar regardless of the group global environment. The other group member galaxies exhibit only the local influence of the group they inhabit. Our analysis comprises tests against the dependence on the host group luminosity and we argue that group brightest member properties are not only determined by the host halo, but also by the large-scale structure which can influence the accretion process onto their late-type brightest galaxies.
The supernova remnant G15.4+0.1 is considered to be the possible counterpart of the gamma-ray source HESSJ1818-154. With the goal of getting a complete view of this remnant and understanding the nature of the gamma-ray flux, we conducted a detailed radio study that includes the search for pulsations and a model of the broadband emission for the G15.4+0.1/HESSJ1818-154 system. Low-frequency imaging at 624 MHz and pulsar observations at 624 and 1404 MHz towards G15.4+0.1 were carried out with the Giant Metrewave Radio Telescope (GMRT). We correlated the new radio data with observations of the source at X-ray and infrared wavelengths from XMM-Newton and Herschel observatories, respectively. To characterize the neutral hydrogen medium (HI) towards G15.4+0.1, we used data from the Southern Galactic Plane Survey. We modelled the spectral energy distribution using both hadronic and leptonic scenarios. From the combination of the new GMRT observations with existing data, we derived a continuum spectral index alpha=-0.62+-0.03 for the whole remnant. The local synchrotron spectra of G15.4+0.1, calculated from the combination of the GMRT data with 330 MHz observations from the VLA, tends to be flatter in the central part of the remnant, accompanying the region where the blast wave is impinging molecular gas. No spectral index trace was found indicating the radio counterpart to the pulsar wind nebula proposed from X-ray observations. In addition, the search for radio pulsations yielded negative results. Emission at far-infrared wavelengths is observed in the region where the SNR shock is interacting with dense molecular clumps. We also identified HI features forming a shell that wraps most of the outer border of G15.4+0.1. Characteristic parameters were estimated for the shocked HI gas. We found that either a purely hadronic or leptonic model is compatible with the broadband emission known so far.
The recent discover of intermediate mass black holes (IMBHs) and dense compact objects in the centre of several galactic and extra-galactic stellar clusters opened a series of questions about their formation and evolution. One possibility is that gravitational encounters carry heavy stars to concentrate toward the cluster centre, leading to the formation of a compact sub-system composed mainly by high-mass stars. Correlations between those compact objects and the cluster in which they are contained could help to constraint their formation scenario. Using theoretical and statistical arguments, we show here that gravitational encounters lead to the formation of a dense system placed at the centre of the host stellar cluster. Moreover, we used one-to-one N -body realization of star clusters to follow the core formation process and its evolution. We derive scaling relations connecting the cluster mass and the deposited, central mass, comparing such results with observations. This work is the first showing that mass segregation in a stellar cluster gives rise, within the cluster centre, to a massive core whose mass correlates with the cluster mass as observations suggest.
The complex spectra of high energy ${\rm e^\pm}$ cosmic rays (CRs) observed near Earth are those expected from standard model physics. In particular, the observed hardening of their spectra with increasing energy reported by the AMS-02 collaboration can be produced by the transition of their energy-loss by inverse Compton scattering off Galactic light from the Thomson to the Klein-Nishina regime. The "cut-off" near TeV in the combined ${\rm e^\pm}$ flux observed with H.E.S.S can be due to pair production in ${\rm e^\pm}\gamma$ collisions in source rather than a bump produced by the decay/annihilation of dark matter particles with a mass of $\sim$TeV. Beyond this "cutoff", the $e^\pm$ CRs are mostly produced by the decay of mesons from hadronic collisions of CR protons in/near source with a positron fraction $\sim 0.57$.
The Yarkovsky effect describes a small but significant force that affects the orbital motion of meteoroids and asteroids smaller than $30-40$ kilometers in diameter. It is caused by sunlight; when these bodies heat up in the Sun, they eventually re-radiate the energy away in the thermal waveband, which in turn creates a tiny thrust. This recoil acceleration is much weaker than solar and planetary gravitational forces, but it can produce measurable orbital changes over decades and substantial orbital effects over millions to billions of years. The same physical phenomenon also creates a thermal torque that, complemented by a torque produced by scattered sunlight, can modify the rotation rates and obliquities of small bodies as well. This rotational variant has been coined the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. During the past decade or so, the Yarkovsky and YORP effects have been used to explore and potentially resolve a number of unsolved mysteries in planetary science dealing with small bodies. Here we review the main results to date, and preview the goals for future work.
The standard theory of electromagnetic cascades onto a photon background predicts a quasi-universal shape for the resulting non-thermal photon spectrum. This has been applied to very disparate fields, including non-thermal big bang nucleosynthesis (BBN). However, once the energy of the injected photons falls below the pair-production threshold the spectral shape is very different, a fact that has been overlooked in past literature. This loophole may have important phenomenological consequences, since it generically alters the BBN bounds on non-thermal relics: for instance it allows to re-open the possibility of purely electromagnetic solutions to the so-called "cosmological lithium problem", which were thought to be excluded by other cosmological constraints. We show this with a proof-of-principle example and a simple particle physics model, compared with previous literature.
We show that the ratio between the stellar mass of central galaxy and the mass of its host halo, $f_c \equiv M_{*,c}/M_{\rm h}$, can be used as an observable proxy of halo assembly time, in that galaxy groups with higher $f_c$ assembled their masses earlier. Using SDSS groups of Yang et al., we study how $f_c$ correlates with galaxy properties such as color, star formation rate, metallicity, bulge to disk ratio, and size. Central galaxies of a given stellar mass in groups with $f_c>0.02$ tend to be redder in color, more quenched in star formation, smaller in size, and more bulge dominated, as $f_c$ increases. The trends in color and star formation appear to reverse at $f_c<0.02$, reflecting a down-sizing effect that galaxies in massive halos formed their stars earlier although the host halos themselves assembled later (lower $f_c$). No such reversal is seen in the size of elliptical galaxies, suggesting that their assembly follows halo growth more closely than their star formation. Satellite galaxies of a given stellar mass in groups of a given halo mass tend to be redder in color, more quenched in star formation and smaller in size as $f_c$ increases. For a given stellar mass, satellites also tend to be smaller than centrals. The trends are stronger for lower mass groups. For groups more massive than $\sim 10^{13}{\rm M}_\odot$, a weak reversed trend is seen in color and star formation. The observed trends in star formation are qualitatively reproduced by an empirical model based on halo age abundance matching, but not by a semi-analytical model tested here.
We have measured v sin i with high precision for a sample of dM3 stars (86
targets). We detected rotation in 82 stars (73 dM3 stars, and 9 dM3e stars). We
compare our measurements of v sin i for all the stars in our dM0, dM2, dM3 and
dM4 samples to those from other authors. We find a good agreement down to v sin
i values of less than 1 km/s. The mean of the differences between measurements
is only 0.42 km/s.
We find that the distribution of P/sini for our dM3 stars is different from
the distribution of P/sini among our samples of dM2 and dM4 stars. The mean
rotation rate for the dM3 stars (excluding dM3e and sdM3 stars) is
significantly slower (25.8 days) than for dM2 (14.4 days) and dM4 stars (11.4
days). Analogous behavior also emerges among the faster rotators (dMe stars):
we find that a longer rotation period also occurs at spectral sub-type dM3e.
Our data suggest that, as regards the rotational properties of lower main
sequence stars, spectral sub-type dM3 stands out as exhibiting unusual slow
rotation compared to that of adjoining sub-types. Our data lead us to suggest
that the unusual rotational properties of M3 dwarfs may represent a signature
of the transition to complete convection (TTCC).
We introduce a numerical method to integrate tidal effects on collisional systems, using any definition of the external potential as a function of space and time. Rather than using a linearisation of the tidal field, this new method follows a differential technique to numerically evaluate the tidal acceleration and its time derivative. Theses are then used to integrate the motions of the components of the collisional systems, like stars in star clusters, using a predictor-corrector scheme. The versatility of this approach allows the study of star clusters, including their tidal tails, in complex, multi-components, time-evolving external potentials. The method is implemented in the code nbody6 (Aarseth 2003).
We present a 67--93.6 GHz spectral line survey of Orion-KL with the new 4 mm Receiver on the Green Bank Telescope (GBT). The survey reaches unprecedented depths and covers the low-frequency end of the 3 mm atmospheric window which has been relatively unexplored previously. The entire spectral-line survey is published electronically for general use by the astronomical community. The calibration and performance of 4 mm Receiver on the GBT is also summarized.
We present high-resolution rotation curves and mass models of 26 dwarf galaxies from LITTLE THINGS. LITTLE THINGS is a high-resolution Very Large Array HI survey for nearby dwarf galaxies in the local volume within 11 Mpc. The rotation curves of the sample galaxies derived in a homogeneous and consistent manner are combined with Spitzer archival 3.6 micron and ancillary optical U, B, and V images to construct mass models of the galaxies. We decompose the rotation curves in terms of the dynamical contributions by baryons and dark matter halos, and compare the latter with those of dwarf galaxies from THINGS as well as Lambda CDM SPH simulations in which the effect of baryonic feedback processes is included. Being generally consistent with THINGS and simulated dwarf galaxies, most of the LITTLE THINGS sample galaxies show a linear increase of the rotation curve in their inner regions, which gives shallower logarithmic inner slopes alpha of their dark matter density profiles. The mean value of the slopes of the 26 LITTLE THINGS dwarf galaxies is alpha =-0.32 +/- 0.24 which is in accordance with the previous results found for low surface brightness galaxies (alpha = -0.2 +/- 0.2) as well as the seven THINGS dwarf galaxies (alpha =-0.29 +/- 0.07). However, this significantly deviates from the cusp-like dark matter distribution predicted by dark-matter-only Lambda CDM simulations. Instead our results are more in line with the shallower slopes found in the Lambda CDM SPH simulations of dwarf galaxies in which the effect of baryonic feedback processes is included. In addition, we discuss the central dark matter distribution of DDO 210 whose stellar mass is relatively low in our sample to examine the scenario of inefficient supernova feedback in low mass dwarf galaxies predicted from recent Lambda SPH simulations of dwarf galaxies where central cusps still remain.
We report the detection of giant pulse emission from PSR B0950+08 in 24 hours of observations made at 39.4 MHz, with a bandwidth of 16 MHz, using the first station of the Long Wavelength Array, LWA1. We detected 119 giant pulses from PSR B0950+08 (at its dispersion measure), which we define as having SNRs at least 10 times larger than for the mean pulse in our data set. These 119 pulses are 0.035\% of the total number of pulse periods in the 24 hours of observations. The rate of giant pulses is about 5.0 per hour. The cumulative distribution of pulse strength $S$ is a steep power law, $N(>S)\propto S^{-4.7}$, but much less steep than would be expected if we were observing the tail of a Gaussian distribution of normal pulses. We detected no other transient pulses in a dispersion measure range from 1 to 90 pc cm$^{-3}$, in the beam tracking PSR B0950+08. The giant pulses have a narrower temporal width than the mean pulse (17.8 ms, on average, vs.\ 30.5 ms). The pulse widths are consistent with a previously observed weak dependence on observing frequency, which may be indicative of a deviation from a Kolmogorov spectrum of electron density irregularities along the line of sight. The rate and strength of these giant pulses is less than has been observed at $\sim$100 MHz. Additionally, the mean (normal) pulse flux density we observed is less than at $\sim$100 MHz. These results suggest this pulsar is weaker and produces less frequent giant pulses at 39 MHz than at 100 MHz.
We report the discovery and follow-up observations of a system of three objects identified by the ALFALFA extragalactic HI survey, cataloged as (almost) dark extragalactic sources, i.e., extragalactic HI detections with no discernible counterpart in publicly available, wide-field, imaging surveys. We have obtained deep optical imaging with WIYN pODI and HI synthesis maps with WSRT of the HI1232+20 system. The source with the highest HI flux has a newly discovered ultra-low surface brightness (LSB) optical counterpart associated with it, while the other two sources have no detected optical counterparts in our images. Our optical observations show that the detected LSB optical counterpart has a peak surface brightness of ~26.4 mag/arcsec^2 in g', which is exceptionally faint. This source (AGC 229385) has the largest accurately measured HI mass-to-light ratio of an isolated object: MHI/Lg'=46 Msun/Lsun, and has an HI mass of 7.2*10^8 Msun. The other two HI sources (with HI masses 2.0*10^8 and 1.2*10^8 Msun) without optical counterparts have upper limit surface brightnesses of 27.9 and 27.8 mag/arcsec^2 in g', and lower limits on their gas mass-to-light ratio of MHI/Lg'>57 and >31 Msun/Lsun. This system lies relatively close in projection to the Virgo Cluster, but velocity flow models indicate that it is located at ~25 Mpc, substantially beyond Virgo. The system appears to be quite isolated, with no known object closer than 500 kpc. These HI sources may represent both sides of the threshold between "dark" star-less galaxies and galaxies with stellar populations. We discuss a variety of possible formation scenarios for the HI1232+20 system.
About 4-5% of chondrules are compound: two separate chondrules stuck together. This is commonly believed to be the result of the two component chondrules having collided shortly after forming, while still molten. This allows high velocity impacts to result in sticking. However, at T ~ 1100K, the temperature below which chondrules collide as solids (and hence usually bounce), coalescence times for droplets of appropriate composition are measured in tens of seconds. Even at 1025K, at which temperature theory predicts that the chondrules must have collided extremely slowly to have stuck together, the coalescence time scale is still less than an hour. These coalescence time scales are too short for the collision of molten chondrules to explain the observed frequency of compound chondrules. We suggest instead a scenario where chondrules stuck together in slow collisions while fully solid; and the resulting chondrule pair was subsequently briefly heated to a temperature in the range of 900-1025K. In that temperature window the coalescence time is finite but long, covering a span of hours to a decade. This is particularly interesting because those temperatures are precisely the critical window for thermally ionized MRI activity, so compound chondrules provide a possible probe into that vital regime.
Recent work by Levitan et al has expanded the long-term photometric database for AM CVn stars. In particular, their outburst properties are well-correlated with orbital period, and allow constraints to be placed on the secular mass transfer rate between secondary and primary if one adopts the disk instability model for the outbursts. We use the observed range of outbursting behavior for AM CVn systems as a function of orbital period to place a constraint on mass transfer rate versus orbital period P. We infer a rate ~5 x 10^{-9} Msun/yr (P/1000 s)^{-5.2}. We show the functional form so obtained is consistent with the recurrence time-orbital period relation found by Levitan et al using a simple theory for the recurrence time. Also, we predict their steep dependence of outburst duration on orbital period will flatten considerably once the longer orbital period systems have more complete observations.
Global magnetic field extrapolations are now revealing the huge complexity of the Sun's corona, and in particular the structure of the boundary between open and closed magnetic flux. Moreover, recent developments indicate that magnetic reconnection in the corona likely occurs in highly fragmented current layers, and that this typically leads to a dramatic increase in the topological complexity beyond that of the equilibrium field. In this paper we investigate the consequences of reconnection at the open-closed flux boundary ("interchange reconnection") in a fragmented current layer. We demonstrate that it leads to a situation in which magnetic flux (and therefore plasma) from open and closed field regions is efficiently mixed together. This corresponds to an increase in the length and complexity of the open-closed boundary. Thus, whenever reconnection occurs at a null point or separator of the open-closed boundary, the associated separatrix arc of the so-called "S-web" in the high corona becomes not a single line but a band of finite thickness within which the open-closed flux boundary is highly structured. This has significant implications for the acceleration of the slow solar wind, for which the interaction of open and closed field is thought to be important, and may also explain the coronal origins of certain solar energetic particles. The topological structures examined contain magnetic null points, separatrices and separators, and include a model for a pseudo-streamer. The potential for understanding both the large scale morphology and fine structure observed in flare ribbons associated with coronal nulls is also discussed.
Despite their importance, very few RR Lyrae (RRL) stars have been known to reside in binary systems. We report on a search for binary RRL in the OGLE-III Galactic bulge data. Our approach consists in the search for evidence of the light-travel time effect in so-called observed minus calculated ($O-C$) diagrams. Analysis of 1952 well-observed fundamental-mode RRL in the OGLE-III data revealed an initial sample of 29 candidates. We used the recently released OGLE-IV data to extend the baselines up to 17 years, leading to a final sample of 12 firm binary candidates. We provide $O-C$ diagrams and binary parameters for this final sample, and also discuss the properties of 8 additional candidate binaries whose parameters cannot be firmly determined at present. We also estimate that $\gtrsim 4$ per cent of the RRL reside in binary systems.
The Pierre Auger Observatory, located on a vast, high plain in western Argentina, is the world's largest cosmic ray observatory. The objectives of the Observatory are to probe the origin and characteristics of cosmic rays above $10^{17}$ eV and to study the interactions of these, the most energetic particles observed in nature. The Auger design features an array of 1660 water-Cherenkov particle detector stations spread over 3000 km$^2$ overlooked by 24 air fluorescence telescopes. In addition, three high elevation fluorescence telescopes overlook a 23.5 km$^2$, 61 detector infill array. The Observatory has been in successful operation since completion in 2008 and has recorded data from an exposure exceeding 40,000 km$^2$ sr yr. This paper describes the design and performance of the detectors, related subsystems and infrastructure that make up the Auger Observatory.
A localized perturbation of a magnetic flux tube produces a pair of wave trains that propagate in opposite directions along the tube. These wave packets disperse as they propagate, where the extent of dispersion depends on the physical properties of the magnetic structure, on the length of the initial excitation, and on its nature (e.g., transverse or axisymmetric). In Oliver et al. (2014) we considered a transverse initial perturbation, whereas the temporal evolution of an axisymmetric one is examined here. In both papers we use a method based on Fourier integrals to solve the initial value problem. Previous studies on wave propagation in magnetic wave guides have emphasized that the wave train dispersion is influenced by the particular dependence of the group velocity on the longitudinal wavenumber. Here we also find that long initial perturbations result in low amplitude wave packets and that large values of the magnetic tube to environment density ratio yield longer wave trains. To test the detectability of propagating transverse or axisymmetric wave packets in magnetic tubes of the solar atmosphere (e.g., coronal loops, spicules, or prominence threads) a forward modelling of the perturbations must be carried out. This is left for a future work.
Theories where the Planck scale is dynamically generated from dimensionless interactions provide predictive inflationary potentials and super-Planckian field variations. We first study the minimal single-field realisation in the low-energy effective field theory limit, finding the predictions $n_s \approx 0.96$ for the spectral index and $r \approx 0.13$ for the tensor-to-scalar ratio, close to those of a quadratic potential. Next we consider agravity as a dimensionless quantum gravity theory finding a multi-field inflation that converges towards an attractor trajectory that predicts $n_s\approx 0.96$ and $0.003<r<0.13$, interpolating between the quadratic and Starobinsky inflation. These theories relate the smallness of the weak scale to the smallness of inflationary perturbations: both arise naturally because of small couplings, implying a reheating temperature of $10^{7-9}$ GeV. A measurement of $r$ by Keck/Bicep3 would give us information on quantum gravity in the dimensionless scenario.
We study the Q-ball dark matter in the gauge-mediated supersymmetry breaking, and seek for the detection possibility in the IceCube experiment. We find that the Q balls would be the dark matter in the parameter region different from that for the gravitino dark matter. In particular, the Q ball is a good dark matter candidate for low reheating temperature, which may be suitable for the Affleck-Dine baryogenesis and/or nonthermal leptogenesis. The dark matter Q balls are detectable by the IceCube-like experiments in the future, which is the peculiar feature compared to the case of the gravitino dark matter.
We calculate the energy per particle of symmetric nuclear matter and pure neutron matter using the microscopic many-body Brueckner-Hartree-Fock (BHF) approach and employing the Argonne V18 (AV18) nucleon-nucleon (NN) potential supplemented with two different three-nucleon force models recently constructed to reproduce the binding energy of $^3$H, $^3$He and $^4$He nuclei as well as the neutron-deuteron doublet scattering length. We find that none of these new three-nucleon force models is able to reproduce simultaneously the empirical saturation point of symmetric nuclear matter and the properties of three- and four-nucleon systems.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)
Ark 564 (z=0.0247) is an X-ray bright NLS1. By using advanced X-ray timing techniques, Legg et al. (2012) discovered an excess of "delayed" emission in the hard X-ray band (4-7.5 keV) following about 1000 seconds after "flaring" light in the soft X-ray band (0.4-1 keV). We report on the X-ray spectral analysis of eight XMM-Newton and one Suzaku observation of Ark 564. High-resolution spectroscopy was performed with the RGS in the soft X-ray band, while broad-band spectroscopy was performed with the EPIC-pn and XIS/PIN instruments. We analysed time-averaged, flux-selected, and time-resolved spectra. Despite the large variability in flux, the broad band spectral shape of Ark 564 is not dramatically varying and can be reproduced either by a superposition of a power law and a blackbody emission, or by a Comptonized power law emission model. High resolution spectroscopy revealed the presence of ionised gas along the line of sight at the systemic redshift of the source, with a low column density and a range of ionisation states. Broad band spectroscopy revealed a very steep intrinsic continuum and a rather weak emission feature in the iron K band; modelling this feature with a reflection component requires the presence of highly ionised gas. Either a reflection-dominated or an absorption-dominated model are able to well reproduce the time-averaged data from a statistical point of view, in both cases requiring contrived geometries and/or unlikely physical parameters. Finally, through time-resolved analysis we spectroscopically identified the "delayed" emission discovered by Legg et al. (2012) as a spectral hardening above ~4 keV; the most likely interpretation for this component is reprocessing of the "flaring" light by gas located at 10-100 r_g from the central supermassive black hole and so hot to be able to Compton upscatter the flaring intrinsic continuum emission.
We have constructed merger trees for galaxies in the Illustris Simulation by directly tracking the baryonic content of subhalos. These merger trees are used to calculate the galaxy-galaxy merger rate as a function of descendant stellar mass, progenitor stellar mass ratio, and redshift. We demonstrate that the most appropriate definition for the mass ratio of a galaxy-galaxy merger consists in taking both progenitor masses at the time when the secondary progenitor reaches its maximum stellar mass. Additionally, we avoid effects from `orphaned' galaxies by allowing some objects to `skip' a snapshot when finding a descendant, and by only considering mergers which show a well-defined `infall' moment. Adopting these definitions, we obtain well-converged predictions for the galaxy-galaxy merger rate with the following main features, which are qualitatively similar to the halo-halo merger rate except for the last one: a strong correlation with redshift that evolves as $\sim (1+z)^{2.4-2.8}$, a power law with respect to mass ratio, and an increasing dependence on descendant stellar mass, which steepens significantly for descendant stellar masses greater than $\sim 2 \times 10^{11} \, {\rm M_{\odot}}$. These trends are consistent with observational constraints for medium-sized galaxies ($M_{\ast} \gtrsim 10^{10} \, {\rm M_{\odot}}$), but in tension with some recent observations of the close pair fraction for massive galaxies ($M_{\ast} \gtrsim 10^{11} \, {\rm M_{\odot}}$), which report a nearly constant or decreasing evolution with redshift. Finally, we provide a fitting function for the galaxy-galaxy merger rate which is accurate over a wide range of stellar masses, progenitor mass ratios, and redshifts.
We have serendipitously identified the first lithium-rich giant star located close to the red giant branch bump in a globular cluster. Through intermediate-resolution FLAMES spectra we derived a lithium abundance of A(Li)=2.55 (assuming local thermodynamical equilibrium), which is extremely high considering the star's evolutionary stage. Kinematic and photometric analysis confirm the object as a member of the globular cluster NGC 362. This is the fourth Li-rich giant discovered in a globular cluster but the only one known to exist at a luminosity close to the bump magnitude. The three previous detections are clearly more evolved, located close to, or beyond the tip of their red giant branch. Our observations are able to discard the accretion of planets/brown dwarfs, as well as an enhanced mass-loss mechanism as a formation channel for this rare object. Whilst the star sits just above the cluster bump luminosity, its temperature places it towards the blue side of the giant branch in the colour-magnitude diagram. We require further dedicated observations to unambiguously identify the star as a red giant: we are currently unable to confirm whether Li production has occurred at the bump of the luminosity function or if the star is on the pre zero-age horizontal branch. The latter scenario provides the opportunity for the star to have synthesised Li rapidly during the core helium flash or gradually during its red giant branch ascent via some extra mixing process.
Currently-proposed galaxy quenching mechanisms predict very different behaviours during major halo mergers, ranging from significant quenching enhancement (e.g., clump-induced gravitational heating models) to significant star formation enhancement (e.g., gas starvation models). To test real galaxies' behaviour, we present an observational galaxy pair method for selecting galaxies whose host haloes are preferentially undergoing major mergers. Applying the method to central L* (10^10 Msun < M* < 10^10.5 Msun) galaxies in the Sloan Digital Sky Survey (SDSS) at z<0.06, we find that major halo mergers can at most modestly reduce the star-forming fraction, from 59% to 47%. Consistent with past research, however, mergers accompany enhanced specific star formation rates for star-forming L* centrals: ~10% when a paired galaxy is within 200 kpc (approximately the host halo's virial radius), climbing to ~70% when a paired galaxy is within 30 kpc. No evidence is seen for even extremely close pairs (<30 kpc separation) rejuvenating star formation in quenched galaxies. For galaxy formation models, our results suggest: (1) quenching in L* galaxies likely begins due to decoupling of the galaxy from existing hot and cold gas reservoirs, rather than a lack of available gas or gravitational heating from infalling clumps, (2) state-of-the-art semi-analytic models currently over-predict the effect of major halo mergers on quenching, and (3) major halo mergers can trigger enhanced star formation in non-quenched central galaxies.
This paper introduces the multiband periodogram, a general extension of the well-known Lomb-Scargle approach for detecting periodic signals in time-domain data. In addition to advantages of the Lomb-Scargle method such as treatment of non-uniform sampling and heteroscedastic errors, the multiband periodogram significantly improves period finding for randomly sampled multiband light curves (e.g., Pan-STARRS, DES and LSST). The light curves in each band are modeled as arbitrary truncated Fourier series, with the period and phase shared across all bands. The key aspect is the use of Tikhonov regularization which drives most of the variability into the so-called base model common to all bands, while fits for individual bands describe residuals relative to the base model and typically require lower-order Fourier series. This decrease in the effective model complexity is the main reason for improved performance. We use simulated light curves and randomly subsampled SDSS Stripe 82 data to demonstrate the superiority of this method compared to other methods from the literature, and find that this method will be able to efficiently determine the correct period in the majority of LSST's bright RR Lyrae stars with as little as six months of LSST data. A Python implementation of this method, along with code to fully reproduce the results reported here, is available on GitHub.
${\it Swift}$ J2058.4+0516 (Sw J2058+05, hereafter) has been suggested as the second member (after Sw J1644+57) of the rare class of tidal disruption events accompanied by relativistic ejecta. Here we report a multiwavelength (X-ray, ultraviolet/optical/infrared, radio) analysis of Sw J2058+05 from 3 months to 3 yr post-discovery in order to study its properties and compare its behavior with that of Sw J1644+57. Our main results are as follows. (1) The long-term X-ray light curve of Sw J2058+05 shows a remarkably similar trend to that of Sw J1644+57. After a prolonged power-law decay, the X-ray flux drops off rapidly by a factor of $\gtrsim 160$ within a span of $\Delta$$t$/$t$ $\le$ 0.95. Associating this sudden decline with the transition from super-Eddington to sub-Eddington accretion, we estimate the black hole mass to be in the range of $10^{4-6}$ M$_{\odot}$. (2) We detect rapid ($\lesssim 500$ s) X-ray variability before the dropoff, suggesting that, even at late times, the X-rays originate from close to the black hole (ruling out a forward-shock origin). (3) We confirm using ${\it HST}$ and VLBA astrometry that the location of the source coincides with the galaxy's center to within $\lesssim 400$ pc (in projection). (4) We modeled Sw J2058+05's ultraviolet/optical/infrared spectral energy distribution with a single-temperature blackbody and find that while the radius remains more or less constant at a value of $63.4 \pm 4.5$ AU ($\sim 10^{15}$ cm) at all times during the outburst, the blackbody temperature drops significantly from $\sim$ 30,000 K at early times to a value of $\sim$ 15,000 K at late times (before the X-ray dropoff). Our results strengthen Sw J2058+05's interpretation as a tidal disruption event similar to Sw J1644+57. For such systems, we suggest the rapid X-ray dropoff as a diagnostic for black hole mass.
With the discovery of Y dwarfs by the WISE mission, the population of field brown dwarfs now extends to objects with temperatures comparable to those of Solar System planets. To investigate the atmospheres of these newly identified brown dwarfs, we have conducted a pilot study monitoring an initial sample of three late T-dwarfs (T6.5, T8 and T8.5) and one Y-dwarf (Y0) for infrared photometric variability at multiple epochs. With J-band imaging, each target was observed for a period of 1.0h to 4.5h per epoch, which covers a significant fraction of the expected rotational period. These measurements represent the first photometric monitoring for these targets. For three of the four targets (2M1047, Ross 458C and WISE0458), multi-epoch monitoring was performed, with the time span between epochs ranging from a few hours to ~2 years. During the first epoch, the T8.5 target WISE0458 exhibited variations with a remarkable min-to-max amplitude of 13%, while the second epoch light curve taken ~2 years later did not note any variability to a 3% upper limit. With an effective temperature of ~600 K, WISE0458 is the coldest variable brown dwarf published to-date, and combined with its high and variable amplitude makes it a fascinating target for detailed follow-up. The three remaining targets showed no significant variations, with a photometric precision between 0.8% and 20.0%, depending on the target brightness. Combining the new results with previous multi-epoch observations of brown dwarfs with spectral types of T5 or later, the currently identified variables have locations on the colour-colour diagram better matched by theoretical models incorporating cloud opacities rather than cloud-free atmospheres. This preliminary result requires further study to determine if there is a definitive link between variability among late-T dwarfs and their location on the colour-colour diagram.
Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical parameters such as the attenuation. Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Halpha, 8, 24, 70, and 100 micron maps, we have systematically compared the emission in individual bands with various SFR estimators from a resolution of 33 pc to 2084 pc. Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators and Halpha+24 micron. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential reddening is compatible with the canonical value found for starburst galaxies. Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be poorly known.
We report small-scale clustering measurements from the PRIMUS spectroscopic redshift survey as a function of color and luminosity. We measure the real-space cross-correlations between 62,106 primary galaxies with PRIMUS redshifts and a tracer population of 545,000 photometric galaxies over redshifts from z=0.2 to z=1. We separately fit a power-law model in redshift and luminosity to each of three independent color-selected samples of galaxies. We report clustering amplitudes at fiducial values of z=0.5 and L=1.5 L*. The clustering of the red galaxies is ~3 times as strong as that of the blue galaxies and ~1.5 as strong as that of the green galaxies. We also find that the luminosity dependence of the clustering is strongly dependent on physical scale, with greater luminosity dependence being found between r=0.0625 Mpc/h and r=0.25 Mpc/h, compared to the r=0.5 Mpc/h to r=2 Mpc/h range. Moreover, over a range of two orders of magnitude in luminosity, a single power-law fit to the luminosity dependence is not sufficient to explain the increase in clustering at both the bright and faint ends at the smaller scales. We argue that luminosity-dependent clustering at small scales is a necessary component of galaxy-halo occupation models for blue, star-forming galaxies as well as for red, quenched galaxies.
Ly$\alpha$ photons scattered by neutral hydrogen atoms in the circumgalactic media or produced in the halos of star-forming galaxies are expected to lead to extended Ly$\alpha$ emission around galaxies. Such low surface brightness Ly$\alpha$ halos (LAHs) have been detected by stacking Ly$\alpha$ images of high-redshift star-forming galaxies. We study the origin of LAHs by performing radiative transfer modeling of nine $z=3.1$ Lyman-Alpha Emitters (LAEs) in a high resolution hydrodynamic galaxy formation simulation. We develop a method of computing the mean Ly$\alpha$ surface brightness profile of each LAE by effectively integrating over many different observing directions. Without adjusting any parameters, our model yields an average Ly$\alpha$ surface brightness profile in remarkable agreement with observations. We find that observed LAHs can not be accounted for solely by photons originating from the central LAE and scattered to large radii by hydrogen atoms in the circumgalactic gas. Instead, Ly$\alpha$ emission from regions in the outer halo is primarily responsible for producing the extended LAHs seen in observations, which potentially includes both star-forming and cooling radiation. The contribution from star formation in the outer halo regions can be strongly constrained to be negligible by the observed absence of an extended ultra-violet (UV) halo. Our results therefore suggest that cooling radiation from the outer halo regions of LAEs plays a major role in forming their extended LAHs. We discuss the implications and caveats of such a picture.
The impending era of wide-field radio surveys has the potential to revolutionize our understanding of astrophysical transients. Here we evaluate the prospects of a wide range of planned and hypothetical radio surveys using the properties and volumetric rates of known and hypothetical classes of extragalactic synchrotron radio transients (e.g., on- and off-axis gamma-ray bursts [GRB], supernovae, tidal disruption events [TDE], compact object mergers). Utilizing these sources and physically motivated considerations we assess the allowed phase-space of radio luminosity and peak timescale for extragalactic transients. We also include for the first time effects such as redshift evolution of the rates, K-corrections, and non-Euclidean luminosity distance, which affect the detection rates of the most sensitive surveys. The number of detected events is calculated by means of a Monte Carlo method, using the various survey properties (depth, cadence, area) and realistic detection criteria that include a cut on the minimum variability of the transients during the survey and an assessment of host galaxy contamination. Near-term GHz frequency surveys (ASKAP/VAST, Very Large Array Sky Survey) will detect few events: <~30-50 on- and off-axis long GRBs and off-axis tidal disruption events, and ~10-20 neutron star binary mergers if ~1% of the mergers result in a stable millisecond magnetar. Low-frequency surveys (e.g., LOFAR) are unlikely to detect any transients, while a hypothetical large-scale mm survey may detect ~40 on-axis long GRBs. On the other hand, SKA surveys at ~0.1-1 GHz have the potential to uncover thousands of transients, mainly on- and off-axis long GRBs, on-axis short GRBs, off-axis TDEs, and neutron star binary mergers with magnetar remnants.
Strong lensing time delay cosmography has excellent complementarity with other dark energy probes, and will soon have abundant systems detected. We investigate two issues in the imaging and spectroscopic followup required to obtain the time delay distance. The first is optimization of spectroscopic resources. We develop a code to optimize the cosmological leverage under the constraint of constant spectroscopic time, and find that sculpting the lens system redshift distribution can deliver a 40% improvement in dark energy figure of merit. The second is the role of systematics, correlated between different quantities of a given system or model errors common to all systems. We show how the levels of different systematics affect the cosmological parameter estimation, and derive guidance for the fraction of double image vs quad image systems to follow as a function of differing systematics between them.
We present a detailed analysis of a H2-bearing metal-rich sub-damped Lyman-alpha system at zabs = 0.10115 towards the radio-loud quasar J0441-4313, at a projected separation of ~7.6 kpc from a star-forming galaxy. The H2, C I and Na I absorption are much stronger in the redder of the two components seen in the Hubble Space Telescope / Cosmic Origins Spectrograph spectrum. The best single component fit to the strong H2 component gives log N(H2) = 16.61 +/- 0.05. However, possible hidden saturation in the medium resolution spectrum can allow for log N(H2) to be as high as 18.9. The rotational excitation temperature of H2 in this component is 133 +33/-22 K. Photoionization models suggest 30-80% of the total N(H I) is associated with the strong H2 component, that has a density <= 100 cm^-3 and is subject to a radiation field that is <= 0.5 times the Galactic mean field. The Very Large Baseline Array 1.4 GHz continuum image of the radio source contains only 27% of the arcsecond scale emission. Using a previously published spectrum, no 21-cm absorption is found to be associated with the strong H2 component. This suggests that either the N(H I) associated with this component is <= 50% of the total N(H I) or the gas covering factor is <= 0.27. This is consistent with the results of the photoionization model that uses UV radiation due to stars in the associated galaxy. The 21-cm absorption previously reported from the weaker H2 component suggests a spin temperature of <= 90 K, at odds with the weakness of H2, C I and Na I absorption in this component. From the inferred physical and chemical conditions, we suggest that the gas may be tracing a recent metal-rich outflow from the host-galaxy.
A large set of 2.5D MHD simulations has been carried out for axisymmetric viscous/diffusive disc accretion to rotating magnetized stars for the purpose of assessing the conditions where the outflows or jets are asymmetric relative to the equatorial plane. Observations of jets from young stellar objects reveal the asymmetric outflows from some sources. The considered initial magnetic fields are symmetric about the equatorial plane and consist of a radially distributed field threading the disc (disc-field) and a stellar dipole field. (1). For pure disc-fields the symmetry or asymmetry of the outflows is affected by the midplane plasma $\beta$ of the disc. For the low density discs with small plasma $\beta$ values, outflows are observed to be symmetric about the equatorial plane to within 10% over timescales of hundreds of inner disc orbits. For the denser higher $\beta$ discs, the coupling of the upper and lower coronal plasmas is broken, and quasi-periodic field motion in the two hemispheres becomes different. This asymmetry leads to asymmetric episodic outflows. (2.) Accreting stars with a stellar dipole field and no disc-field exhibit episodic, two component outflows - a magnetospheric wind and an inner disc wind from somewhat larger radial distances. Both are characterized by similar velocity profiles but the magnetospheric wind has densities 10 times that of the disc wind. The winds are highly asymmetric with outflow from one hemisphere and funnel flow accretion in the opposite hemisphere. (3.) Adding a disc-field which is anti-parallel to the stellar dipole field in the disc acts to suppress the magnetospheric and disc winds. In contrast, adding a disc-field parallel to the stellar dipole field acts to enhance the episodic magnetospheric and disc winds. The winds are again highly asymmetric about the equatorial plane. They are not influenced by the initial plasma $\beta$ of the disc.
Distant regions close to the plane of our Galaxy are largely unexplored by optical surveys as they are hidden by dust. We have used near-infrared data (that minimizes dust obscuration) from the ESO Public survey VISTA Variables of the Via Lactea (VVV) (Minniti et al. 2011; Saito et al. 2012; henceforth S12) to search for distant stars at low latitudes. We have discovered four Cepheid variables within an angular extent of one degree centered at Galactic longitude of $l = -27.4^\circ$ and Galactic latitude of $b = -1.08 ^\circ$. We use the tightly constrained period-luminosity relationship that these pulsating stars obey (Persson et al. 2004; Matsunaga et al. 2011) to derive distances. We infer an average distance to these Cepheid variables of 90 kpc. The Cepheid variables are highly clustered in angle (within one degree) and in distance (the standard deviation of the distances is 12 kpc). They are at an average distance of $\sim 2~\rm kpc$ from the plane and their maximum projected separation is $\sim 1~ \rm kpc$. These young ($\sim$ 100 Myr old), pulsating stars (Bono et al. 2005) are unexpected at such large distances from the Galactic disk, which terminates at $\sim$ 15 kpc (Minniti et al. 2011). The highly clustered nature in distance and angle of the Cepheid variables suggests that the stars may be associated with a dwarf galaxy, one that was earlier predicted by a dynamical analysis (Chakrabarti \& Blitz 2009).
Yellowballs are a collection of approximately 900 compact, infrared sources identified and named by volunteers participating in the Milky Way Project (MWP), a citizen-science project that uses GLIMPSE/MIPSGAL images from Spitzer to explore topics related to Galactic star formation. In this paper, through a combination of catalog cross-matching and infrared color analysis, we show that yellowballs are a mix of compact star-forming regions, including ultra-compact and compact HII regions, as well as analogous regions for less massive B-type stars. The resulting MWP yellowball catalog provides a useful complement to the Red MSX Source (RMS) survey. It similarly highlights regions of massive star formation, but the selection of objects purely on the basis of their infrared morphology and color in Spitzer images identifies a signature of compact star-forming regions shared across a broad range of luminosities, and by inference, masses. We discuss the origin of their striking mid-infrared appearance, and suggest that future studies of the yellowball sample will improve our understanding of how massive and intermediate-mass star-forming regions transition from compact to more extended bubble-like structures.
The origin of Galactic cosmic rays remains unconfirmed, but promising candidates for their sources are found in star-forming regions. We report a series of X-ray observations, with Suzaku, toward the nearby star-forming region of Cygnus X. They aim at comparing diffuse X-ray emissions on and off the $\gamma$-ray cocoon of hard cosmic rays revealed by Fermi LAT. After excluding point sources and small-scale structures and subtracting the non-X-ray and cosmic X-ray backgrounds, the 2--10~keV X-ray intensity distribution is found to monotonically decrease with increasing Galactic latitude. This indicates that most of the extended emission detected by Suzaku originates from the Galactic ridge. In two observations, we derive upper limits of $3.4 \times 10^{-8}~{\rm erg~s^{-1}~cm^{-2}~sr^{-1}}$ and $1.3 \times 10^{-8}~{\rm erg~s^{-1}~cm^{-2}~sr^{-1}}$ to X-ray emission in the 2--10 keV range from the gamma-ray cocoon. These limits exclude the presence of cosmic-ray electrons with energies above about 50 TeV at a flux level capable of explaining the gamma-ray spectrum. They are consistent with the emission cut-off observed near a TeV in gamma rays. The properties of Galactic ridge and local diffuse X-rays are also discussed.
Magnetic reconnection is a process of magnetic field topology change, which is one of the most fundamental processes in magnetized plasmas. In most astrophysical environments the Reynolds numbers are large and therefore the transition to turbulence is inevitable. This turbulence must be taken into account for any theory of magnetic reconnection, since the initially laminar configurations can transit to the turbulence state, what is demonstrated by 3D high resolution numerical simulations. We discuss ideas of how turbulence can modify reconnection with the focus on the Lazarian & Vishniac (1999) reconnection model and present numerical evidence supporting the model and demonstrate that it is closely connected to the concept of Richardson diffusion and compatible with the Lagrangian dynamics of magnetized fluids. We point out that the Generalized Ohm's Law, that accounts for turbulent motion, predicts the subdominance of the microphysical plasma effects for a realistically turbulent media. We show that on of the most dramatic consequences of turbulence is the violation of the generally accepted notion of magnetic flux freezing. This notion is a corner stone of most theories dealing with magnetized plasmas and therefore its change induces fundamental shifts in accepted paradigms like turbulent reconnection entailing the diffusion process that is essential for understanding star formation. We argue, that at sufficiently high Reynolds numbers the process of tearing reconnection should transfer to turbulent reconnection. We discuss flares predicted by turbulent reconnection and relate them to solar flares and gamma ray bursts. We analyze solar observations, measurements in the solar wind or heliospheric current sheet, and show their correspondence with turbulent reconnection predictions. Finally, we discuss 1st Order Fermi acceleration as a natural consequence of the turbulent reconnection.
High-magnification gravitational microlensing events provide an important channel of detecting planetary systems with multiple giants located at their birth places. In order to investigate the potential existence of additional planets, we reanalyze the light curves of the eight high-magnification microlensing events for each of which a single planet was previously detected. The analyzed events include OGLE-2005-BLG-071, OGLE-2005-BLG-169, MOA-2007-BLG-400, MOA-2008-BLG-310, MOA-2009-BLG-319, MOA-2009-BLG-387, MOA-2010-BLG-477, and MOA-2011-BLG-293. We find that including an additional planet improves fits with $\Delta\chi^2 < 80$ for seven out of eight analyzed events. For MOA-2009-BLG-319, the improvement is relatively big with $\Delta\chi^2 \sim 143$. From inspection of the fits, we find that the improvement of the fits is attributed to systematics in data. Although no clear evidence of additional planets is found, it is still possible to constrain the existence of additional planets in the parameter space. For this purpose, we construct exclusion diagrams showing the confidence levels excluding the existence of an additional planet as a function of its separation and mass ratio. We also present the exclusion ranges of additional planets with 90\% confidence level for Jupiter, Saturn, and Uranus-mass planets.
Among the transient black hole binary systems, 4U 1630-47 is one of the most active sources exhibiting outbursts every few hundred days, with every outburst lasting typically around hundred days. During the 2002-2004 outburst the appearance of quasi-periodic oscillations (QPOs) coincide with the onset of anomalous state. There are two distinct QPO states: namely the single QPO state with one QPO and the twin QPO state with two QPOs not related harmonically. The spectral features of this state are corroborated by a previous outburst in 1998. The evolution of the inner disc temperature and the inner disc radius suggest the possible onset of geometrically thicker, slim disc as a possible explanation of the energy spectral features. The other possibilities involving the Comptonizing cloud also exist that may explain the anomalous state. The two different QPO states exhibit different spectral features, and we provide a possible empirical physical scenario from the observation of the evolution of the spectral features.
Based on high-spectral-resolution observations (R=60000) performed with the 6-m BTA telescope in combination with the echelle spectrograph NES, we have studied the optical spectra of three A-type supergiants: a peculiar supergiant 3 Pup, a post-AGB star BD+48 1220, and a massive $\alpha$ Cyg, which belong to essentially different stages of evolution. A spectral atlas for these stars is prepared in the wavelength interval of 3920 to 6720 \AA.
In the framework of relativistic mean field theory, the relativistic neutrino emissivity of the nucleonic and hyperonic direct Urca processes in the degenerate baryon matter of neutron stars are studied. We investigate particularly the influence of the isovector scalar interaction which is considered by exchanging $\delta$ meson on the nucleonic and hyperonic direct Urca processes. The results indicate that $\delta$ mesons lead to obvious enhancement of the total neutrino emissivity, which must result in more rapid cooling rate of neutron star matter.
PSR J2032+4127 is a gamma-ray and radio-emitting pulsar which has been regarded as a young luminous isolated neutron star. However, its recent spin-down rate has extraordinarily increased by a factor of two. We present evidence that this is due to its motion as a member of a highly-eccentric binary system with a 15-solar-mass Be star, MT91~213. Timing observations show that, not only are the positions of the two stars coincident within 0.4 arcsec, but timing models of binary motion of the pulsar fit the data much better than a model of a young isolated pulsar. MT91~213, and hence the pulsar, lie in the Cyg~OB2 stellar association, which is at a distance of only 1.4-1.7 kpc. The pulsar is currently on the near side of, and accelerating towards, the Be star, with an orbital period of 20-30 years. The next periastron is well-constrained to occur in early 2018, providing an opportunity to observe enhanced high-energy emission as seen in other Be-star binary systems.
This HDR-thesis is devoted to the study of the rotation of the natural satellites of the giant planets and of Mercury. These bodies have a resonant rotation. Most of the natural satellites rotate synchronously, showing the same hemisphere to their parent planet (1:1 spin-orbit resonance). The case of Mercury is unique since its spin rate is exactly 1.5 its mean motion (3:2 spin-orbit resonance). These two configurations are dynamical equilibria, reached after damping of the initial rotation of the relevant bodies. Thus, the rotation quantities are a signature of the interior, in particular of a putative global ocean. This manuscript divides into 3 parts. The first part is devoted to the synchronous resonance. It presents different models of rotation from a fully rigid body to a one with a global subsurfacic ocean. We always consider all the degrees of freedom simultaneously, using analytical and numerical resolutions. These models are applied on Titan, Callisto, Janus, Epimetheus, Mimas, Hyperion, and Io. The second part presents the resonant rotation of Mercury, target of the two space missions MESSENGER and BepiColombo. We reveal in particular how it got trapped into its 3:2 resonance. The final part presents an algorithm I have elaborated to tackle the rotational problems.
Context: MWC 656 has recently been established as the first observationally detected high-mass X-ray binary system containing a Be star and a black hole (BH). The system has been associated with a gamma-ray flaring event detected by the AGILE satellite in July 2010. Aims: Our aim is to evaluate if the MWC 656 gamma-ray emission extends to very high energy (VHE > 100 GeV) gamma rays. Methods. We have observed MWC 656 with the MAGIC telescopes for $\sim$23 hours during two observation periods: between May and June 2012 and June 2013. During the last period, observations were performed contemporaneously with X-ray (XMM-Newton) and optical (STELLA) instruments. Results: We have not detected the MWC 656 binary system at TeV energies with the MAGIC Telescopes in either of the two campaigns carried out. Upper limits (ULs) to the integral flux above 300 GeV have been set, as well as differential ULs at a level of $\sim$5\% of the Crab Nebula flux. The results obtained from the MAGIC observations do not support persistent emission of very high energy gamma rays from this system at a level of 2.4\% the Crab flux.
The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.
We present an updated description of the Planck Low Frequency (LFI) data processing pipeline, associated with the 2015 data release. We point out the places in which our results and methods have remained unchanged since the 2013 paper and we highlight the changes made for the 2015 release, describing the products (especially timelines) and the ways in which they were obtained. We demonstrate that the pipeline is self-consistent (principally based on simulations) and report all null tests. We refer to other related papers where more detailed descriptions on the LFI data processing pipeline may be found if needed.
This paper presents the characterization of the in-flight beams, the beam window functions, and the associated uncertainties for the Planck Low Frequency Instrument (LFI). The structure of the paper is similar to that presented in the 2013 Planck release; the main differences concern the beam normalization and the delivery of the window functions to be used for polarization analysis. The in-flight assessment of the LFI main beams relies on measurements performed during observations of Jupiter. By stacking data from seven Jupiter transits, the main beam profiles are measured down to -25 dB at 30 and 44 GHz, and down to -30 dB at 70 GHz. The agreement between the simulated beams and the measured beams is confirmed to be better than 1% at each LFI frequency band (within the 20 dB contour from the peak, the rms values are: 0.1% at 30 and 70 GHz; 0.2% at 44 GHz). Simulated polarized beams are used for the computation of the effective beam window functions. The error budget for the window functions is estimated from both main beam and sidelobe contributions, and accounts for the radiometer band shapes. The total uncertainties in the effective beam window functions are: 0.7% and 1% at 30 and 44 GHz, respectively (at $\ell \approx 600$); and 0.5% at 70 GHz (at $\ell \approx 1000$).
This paper describes the mapmaking procedure applied to Planck LFI (Low Frequency Instrument) data. The mapmaking step takes as input the calibrated timelines and pointing information. The main products are sky maps of $I,Q$, and $U$ Stokes components. For the first time, we present polarization maps at LFI frequencies. The mapmaking algorithm is based on a destriping technique, enhanced with a noise prior. The Galactic region is masked to reduce errors arising from bandpass mismatch and high signal gradients. We apply horn-uniform radiometer weights to reduce effects of beam shape mismatch. The algorithm is the same as used for the 2013 release, apart from small changes in parameter settings. We validate the procedure through simulations. Special emphasis is put on the control of systematics, which is particularly important for accurate polarization analysis. We also produce low-resolution versions of the maps, and corresponding noise covariance matrices. These serve as input in later analysis steps and parameter estimation. The noise covariance matrices are validated through noise Monte Carlo simulations. The residual noise in the map products is characterized through analysis of half-ring maps, noise covariance matrices, and simulations.
The Planck High Frequency Instrument (HFI) has observed the full sky at six frequencies (100, 143, 217, 353, 545, and 857 GHz) in intensity and at four frequencies in linear polarization (100, 143, 217, and 353 GHz). In order to obtain sky maps, the time-ordered information (TOI) containing the detector and pointing samples must be processed and the angular response must be assessed. The full mission TOI is included in the Planck 2015 release. This paper describes the HFI TOI and beam processing for the 2015 release. HFI calibration and map-making are described in a companion paper. The main pipeline has been modified since the last release (2013 nominal mission in intensity only), by including a correction for the non-linearity of the warm readout and by improving the model of the bolometer time response. The beam processing is an essential tool that derives the angular response used in all the Planck science papers and we report an improvement in the effective beam window function uncertainty of more than a factor 10 relative to the 2013 release. Noise correlations introduced by pipeline filtering function are assessed using dedicated simulations. Angular cross-power spectra using datasets which are decorrelated in time are immune to the main systematic effects.
This paper describes the processing applied to the Planck High Frequency Instrument (HFI) cleaned, time-ordered information to produce photometrically calibrated maps in temperature and (for the first time) in polarization. The data from the 2.5 year full mission include almost five independent full-sky surveys. HFI observes the sky over a broad range of frequencies, from 100 to 857 GHz. To get the best accuracy on the calibration over such a large range, two different photometric calibration schemes have been used. The 545 and 857 GHz data are calibrated using models of planetary atmospheric emission. The lower frequencies (from 100 to 353 GHz) are calibrated using the time-variable cosmological microwave background dipole which we call the orbital dipole. This source of calibration only depends on the satellite velocity with respect to the solar system and permits an independent measurement of the amplitude of the CMB solar dipole (3364.5 +/- 0.8 \mu K) which is 1\sigma\ higher than the WMAP measurement with a direction that is consistent between both experiments. We describe the pipeline used to produce the maps of intensity and linear polarization from the HFI timelines, and the scheme used to set the zero level of the maps a posteriori. We also summarize the noise characteristics of the HFI maps in the 2015 Planck data release and present some null tests to assess their quality. Finally, we discuss the major systematic effects and in particular the leakage induced by flux mismatch between the detectors leading to spurious polarization signal.
Planck has mapped the microwave sky in nine frequency bands between 30 and 857 GHz in temperature and seven bands between 30 and 353 GHz in polarization. In this paper we consider the problem of diffuse astrophysical component separation, and process these maps within a Bayesian framework to derive a consistent set of full-sky astrophysical component maps. For the temperature analysis, we combine the Planck observations with the 9-year WMAP sky maps and the Haslam et al. 408 MHz map to derive a joint model of CMB, synchrotron, free-free, spinning dust, CO, line emission in the 94 and 100 GHz channels, and thermal dust emission. Full-sky maps are provided with angular resolutions varying between 7.5 arcmin and 1 deg. Global parameters (monopoles, dipoles, relative calibration, and bandpass errors) are fitted jointly with the sky model, and best-fit values are tabulated. For polarization, the model includes CMB, synchrotron, and thermal dust emission. These models provide excellent fits to the observed data, with rms temperature residuals smaller than 4 uK over 93% of the sky for all Planck frequencies up to 353 GHz, and fractional errors smaller than 1% in the remaining 7% of the sky. The main limitations of the temperature model at the lower frequencies are degeneracies among the spinning dust, free-free, and synchrotron components; additional observations from external low-frequency experiments will be essential to break these. The main limitations of the temperature model at the higher frequencies are uncertainties in the 545 and 857 GHz calibration and zero-points. For polarization, the main outstanding issues are instrumental systematics in the 100-353 GHz bands on large angular scales in the form of temperature-to-polarization leakage, uncertainties in the analog-to-digital conversion, and very long time constant corrections, all of which are expected to improve in the near future.
We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.
We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the cosmological constant scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of state, principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When estimating the density of DE at early times, we significantly improve present constraints. We then move to general parameterizations of the DE or MG perturbations that encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use baryonic acoustic oscillations, type-Ia supernovae and local measurements of the Hubble constant. We further show the impact of measurements of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. All results that include only background parameterizations are in agreement with LCDM. When testing models that also change perturbations (even when the background is fixed to LCDM), some tensions appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck TT+lowP when parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 sigma when external data sets are included. It however disappears when including CMB lensing.
We present the most significant measurement of the cosmic microwave background (CMB) lensing potential to date (at a level of 40 sigma), using temperature and polarization data from the Planck 2015 full-mission release. Using a polarization-only estimator we detect lensing at a significance of 5 sigma. We cross-check the accuracy of our measurement using the wide frequency coverage and complementarity of the temperature and polarization measurements. Public products based on this measurement include an estimate of the lensing potential over approximately 70% of the sky, an estimate of the lensing potential power spectrum in bandpowers for the multipole range 40<L<400 and an associated likelihood for cosmological parameter constraints. We find good agreement between our measurement of the lensing potential power spectrum and that found in the best-fitting LCDM model based on the Planck temperature and polarization power spectra. Using the lensing likelihood alone we obtain a percent-level measurement of the parameter combination Sigma_8 Omega_m^{0.25} = 0.591+-0.021. We combine our determination of the lensing potential with the E-mode polarization also measured by Planck to generate an estimate of the lensing B-mode. We show that this lensing B-mode estimate is correlated with the B-modes observed directly by Planck at the expected level and with a statistical significance of 10 sigma, confirming Planck's sensitivity to this known sky signal. We also correlate our lensing potential estimate with the large-scale temperature anisotropies, detecting a cross-correlation at the 3 sigma level, as expected due to dark energy in the concordance LCDM model.
The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, due to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally-dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]
Full-sky CMB maps from the 2015 Planck release allow us to detect departures from global isotropy on the largest scales. We present the first searches using CMB polarization for correlations induced by a non-trivial topology with a fundamental domain intersecting, or nearly intersecting, the last scattering surface (at comoving distance $\chi_{rec}$). We specialize to flat spaces with toroidal and slab topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection of a compact topology at a scale below the diameter of the last scattering surface. The limits on the radius $R_i$ of the largest sphere inscribed in the topological domain (at log-likelihood-ratio $\Delta\ln{L}>-5$ relative to a simply-connected flat Planck best-fit model) are $R_i>0.97\chi_{rec}$ for the cubic torus and $R_i>0.56\chi_{rec}$ for the slab. The limit for the cubic torus from the matched-circles search is numerically equivalent, $R_i>0.97\chi_{rec}$ (99% CL) from polarisation data alone. We also perform a Bayesian search for a Bianchi VII$_h$ geometry. In the non-physical setting where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour the inclusion of a Bianchi component. However, the cosmological parameters generating this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to Planck data requires an amplitude of $-0.1\pm0.04$ compared to +1 if the model were to be correct. In the physical setting where the Bianchi parameters are fit simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models to $(\omega/H)_0<7.6\times10^{-10}$ (95% CL). [Abridged]
We predict and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG when we consider a maximally helical field. For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.1$ nG and $B_{1\,\mathrm{Mpc}}<0.7$ nG if the impact of PMFs on the ionization history of the Universe is included in the analysis. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfv\'en waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.
This paper presents a study of the ISW effect from the Planck 2015 temperature and polarization data release. The CMB is cross-correlated with different LSS tracers: the NVSS, SDSS and WISE catalogues, and the Planck 2015 convergence lensing map. This cross-correlation yields a detection at $4\,\sigma$, where most of the signal-to-noise is due to the Planck lensing and NVSS. In fact, the ISW effect is detected only from the Planck data (through the ISW-lensing bispectrum) at $\approx 3\,\sigma$, which is similar to the detection level achieved by combining the cross-correlation signal coming from all the catalogues. This cross-correlation analysis is performed only with the Planck temperature data, since the polarization scales available in the 2015 release do not permit significant improvement of the CMB-LSS cross-correlation detectability. Nevertheless, polarization data is used to study the anomalously large ISW signal previously reported through the aperture photometry on stacked CMB features at the locations of known superstructures, which is in conflict with $\Lambda$CDM expectations. We find that the current Planck polarization data do not reject that this signal could be caused by the ISW effect. In addition, the stacking of the Planck lensing map on the superstructures' locations exhibits a positive cross-correlation with these large-scale structures. Finally, we have improved our previous reconstruction of the ISW temperature fluctuations by combining the information encoded in all the previously mentioned LSS tracers. In particular, we construct a map of the ISW secondary anisotropies and the corresponding uncertainties map, obtained from simulations. We also explore the reconstruction of the ISW anisotropies caused by the LSS traced by the 2MPZ survey, by directly inverting the density field into the gravitational potential field.
We have constructed all-sky y-maps of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 30 to 857 GHz frequency channel maps from the Planck satellite survey. These reconstructed y-maps are delivered as part of the Planck 2015 release. The y-maps are characterised in terms of noise properties and residual foreground contamination, mainly thermal dust emission at large angular scales and CIB and extragalactic point sources at small angular scales. Specific masks are defined to minimize foreground residuals and systematics. Using these masks we compute the y-map angular power spectrum and higher order statistics. From these we conclude that the y-map is dominated by tSZ signal in the multipole range, 20-600. We compare the measured tSZ power spectrum and higher order statistics to various physically motivated models and discuss the implications of our results in terms of cluster physics and cosmology.
We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise of six, and is more than a factor of two larger than the 2013 Planck cluster cosmology sample. The counts are consistent with those from 2013 and yield compatible constraints under the same modelling assumptions. Taking advantage of the larger catalogue, we extend our analysis to the two-dimensional distribution in redshift and signal-to-noise. We use mass estimates from two recent studies of gravitational lensing of background galaxies by Planck clusters to provide priors on the hydrostatic bias parameter, $1-b$. In addition, we use lensing of cosmic microwave background (CMB) temperature fluctuations by Planck clusters as a third independent constraint on this parameter. These various calibrations imply constraints on the present-day amplitude of matter fluctuations in varying degrees of tension with those coming from Planck analysis of primary fluctuations in the CMB; for the lowest estimated values of $1-b$ the tension is mild, only a little over one standard deviation, while for the largest estimated value it remains substantial. We also examine constraints on extensions to the base flat $\Lambda CDM$ model by combining the cluster and CMB constraints. The combination appears to favour non-minimal neutrino masses, but this possibility does little to relieve the overall tension because it simultaneously lowers the implied value of the Hubble parameter, thereby exacerbating the discrepancy with most current astrophysical estimates. Improving the precision of cluster mass calibrations from the current 10%-level to 1% would significantly strengthen these combined analyses and provide a stringent test of the base $\Lambda CDM$ model.
We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data-sets, and is the first SZ-selected cluster survey containing > 103 confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the Y5R500 estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires. the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.
We present the Planck Catalogue of Galactic Cold Clumps (PGCC), an all-sky catalogue of Galactic cold clump candidates detected by Planck. This catalogue is the full version of the Early Cold Core (ECC) catalogue, which was made available in 2011 with the Early Release Compact Source Catalogue (ERCSC) and contained 915 high S/N sources. It is based on the Planck 48 months mission data that are currently being released to the astronomical community. The PGCC catalogue is an observational catalogue consisting exclusively of Galactic cold sources. The three highest Planck bands (857, 545, 353 GHz) have been combined with IRAS data at 3 THz to perform a multi-frequency detection of sources colder than their local environment. After rejection of possible extragalactic contaminants, the PGCC catalogue contains 13188 Galactic sources spread across the whole sky, i.e., from the Galactic plane to high latitudes, following the spatial distribution of the main molecular cloud complexes. The median temperature of PGCC sources lies between 13 and 14.5 K, depending on the quality of the flux density measurements, with a temperature ranging from 5.8 to 20 K after removing sources with the 1% largest temperature estimates. Using seven independent methods, reliable distance estimates have been obtained for 5574 sources, which allows us to derive their physical properties such as their mass, physical size, mean density and luminosity. The PGCC sources are located mainly in the solar neighbourhood, up to a distance of 10.5 kpc towards the Galactic centre, and range from low-mass cores to large molecular clouds. Because of this diversity and because the PGCC catalogue contains sources in very different environments, the catalogue is useful to investigate the evolution from molecular clouds to cores. Finally, the catalogue also includes 54 additional sources located in the SMC and LMC.
So far, only one rotating disk has been clearly identified and studied in AGB
or post-AGB objects (in the Red Rectangle), by means of observations with high
spectral and spatial resolution. However, disks are thought to play a key role
in the late stellar evolution and are suspected to surround many evolved stars.
We aim to extend our knowledge on these structures.
We present interferometric observations of CO J=2-1 emission from the nebula
surrounding the post-AGB star AC Her, a source belonging to a class of objects
that share properties with the Red Rectangle and show hints of Keplerian disks.
We clearly detect the Keplerian dynamics of a second disk orbiting an evolved
star. Its main properties (size, temperature, central mass) are derived from
direct interpretation of the data and model fitting. With this we confirm that
there are disks orbiting the stars of this relatively wide class of post-AGB
objects
We show that the cosmic ray (CR) knee can be entirely explained by energy-dependent CR leakage from the Milky Way, with an excellent fit to all existing data. We test this hypothesis calculating the trajectories of individual CRs in the Galactic magnetic field. We find that the CR escape time $\tau_{\rm esc}(E)$ exhibits a knee-like structure around $E/Z={\rm few}\times 10^{15}$ eV for small coherence lengths and strengths of the turbulent magnetic field. The resulting intensities for different groups of nuclei are consistent with the ones determined by KASCADE and KASCADE-Grande, using simple power-laws as injection spectra. The transition from Galactic to extragalactic CRs is terminated at $\approx 2\times 10^{18}$ eV, while extragalactic CRs contribute sizeable to the subdominant proton flux already for $\gtrsim 2\times 10^{16}$ eV. The natural source of extragalactic CRs in the intermediate energy region up to the ankle are in this model normal and starburst galaxies. The escape model provides a good fit to $\ln(A)$ data; it predicts that the phase of the CR dipole varies strongly in the energy range between $1\times 10^{17}$ and $3\times 10^{18}$ eV, while our estimate for the dipole magnitude is consistent with observations.
When the predecoupling plasma is thermodinamically reversible its fluctuations are classified in terms of the adiabatic and entropic modes. A different category of physical solutions, so far unexplored, arises when the inhomogeneities of the viscosity coefficients induce computable curvature perturbations. The viscous modes are explicitly illustrated and compared with the conventional isocurvature solutions.
We present the first high spectral resolution abundance analysis of two newly discovered Galactic globular clusters, namely Mercer 5 and 2MASS GC02 residing in regions of high interstellar reddening in the direction of the Galactic center. The data were acquired with the Phoenix high-resolution near-infrared echelle spectrograph at Gemini South (R~50000) in the 15500.0 A - 15575.0 A spectral region. Iron, Oxygen, Silicon, Titanium and Nickel abundances were derived for two red giant stars, in each cluster, by comparing the entire observed spectrum with a grid of synthetic spectra generated with MOOG. We found [Fe/H] values of -0.86 +/- 0.12 and -1.08 +/- 0.13 for Mercer 5 and 2MASS GC02 respectively. The [O/Fe], [Si/Fe] and [Ti/Fe] ratios of the measured stars of Mercer 5 follow the general trend of both bulge field and cluster stars at this metallicity, and are enhanced by > +0.3. The 2MASS GC02 stars have relatively lower ratios, but still compatible with other bulge clusters. Based on metallicity and abundance patterns of both objects we conclude that these are typical bulge globular clusters.
Asteroids formed in a dynamically quiescent disk but their orbits became gravitationally stirred enough by Jupiter to lead to high-speed collisions. As a result, many dozen large asteroids have been disrupted by impacts over the age of the Solar System, producing groups of fragments known as asteroid families. Here we explain how the asteroid families are identified, review their current inventory, and discuss how they can be used to get insights into long-term dynamics of main belt asteroids. Electronic tables of the membership for 122 notable families are reported on the Planetary Data System node.
OTELO is an emission-line object survey carried out with the red tunable filter of the instrument OSIRIS at the GTC, whose aim is to become the deepest emission-line object survey to date. With 100% of the data of the first pointing finally obtained in June 2014, we present here some aspects of the processing of the data and the very first results of the OTELO survey. We also explain the next steps to be followed in the near future.
Massive and intermediate mass stars play a crucial role in astrophysics. Indeed, massive stars are the main producers of heavy elements, explode in supernovae at the end of their short lifetimes, and may be the progenitors of gamma ray bursts. Intermediate mass stars, although not destined to explode in supernovae, display similar phenomena, are much more numerous, and have some of the richest pulsation spectra. A key to understanding these stars is understanding the effects of rapid rotation on their structure and evolution. These effects include centrifugal deformation and gravity darkening which can be observed immediately, and long terms effects such as rotational mixing due to shear turbulence, which prolong stellar lifetime, modify chemical yields, and impact the stellar remnant at the end of their lifetime. In order to understand these effects, a number of models have been and are being developed over the past few years. These models lead to increasingly sophisticated predictions which need to be tested through observations. A particularly promising source of constraints is seismic observations as these may potentially lead to detailed information on their internal structure. However, before extracting such information, a number of theoretical and observational hurdles need to be overcome, not least of which is mode identification. The present proceedings describe recent progress in modelling these stars and show how an improved understanding of their pulsations, namely frequency patterns, mode visibilities, line profile variations, and mode excitation, may help with deciphering seismic observations.
The ultra-high contrast capability required to form images of other solar systems is arguably the highest-profile challenge in astronomy today. The current high-contrast imaging efforts all require background subtraction to separate the planetary image from the image of the host star. Background estimation is difficult due to the presence of non-common path aberrations (NCPAs) that change with time. The only major source of information that is not being utilized by current efforts is the random encoding of the planetary image and the NCPAs by the atmosphere on millisecond time-scales. Here, a method that utilizes this information in order to avoid background subtraction altogether is proposed. This new paradigm will allow simultaneous estimation of the time-dependent NCPAs and the planetary image via rigorous statistical inference procedures. These procedures are fully compatible with other information sources, such as diurnal field rotation and spectral diversity. Given the open-ended nature of the background subtraction issues, the ideas explained herein may well the key to imaging habitable planets with Extremely Large Telescopes (ELTs). Fully exploiting the information content of millisecond exposures will require significant design modifications of the ELT wavefront sensors and science camera systems, if ultra-high contrast imaging is to be priority.
A novel idea is presented for removing mag- netic flux from the galactic disc in order to satisfy the boundary conditions of the alpha- omega dynamo. The idea involves making use of superbubbles that break out of the galactic disc. When this happens, their shells break up into many fragments, and from these fragments, spikes can arise that can move small pieces of the galactic flux lines in the superbubble shell and move them so high enough into the halo that they do not return. As a result the flux lines remaining in the disc are effectively cut into short pieces of finite length that, although weakly connected to the halo, have no tensile strength at their ends. They are thus free to randomly rotate until no net flux is left. In the first section, the boundary condition problem is discussed and then this solution is discussed.. In the second section the proper- ties of the superbubble that are necessary for the model are laid out. It is shown by a rough numerical estimate, that enough lines can be cut to resolve the dynamo problem. Lastly, conditions under which the spike instability exists are discussed.
A covariant action principle for ideal relativistic magnetohydrodynamics (MHD) in terms of natural Eulerian field variables is given. This is done by generalizing the covariant Poisson bracket theory of Marsden et al., which uses a noncanonical bracket to effect constrained variations of an action functional. Various implications and extensions of this action principle are also discussed. Two significant by-products of this formalism are the introduction of a new divergence-free 4-vector variable for the magnetic field, and a new Lie-dragged form for the theory.
We discuss the dynamics of multiple scalar fields and the possibility of realistic inflation in the maximal gauged supergravity. In this paper, we address this problem in the framework of recently discovered 1-parameter deformation of ${\rm SO}(4,4)$ and ${\rm SO}(5,3)$ dyonic gaugings, for which the base point of the scalar manifold corresponds to an unstable de Sitter vacuum. In the gauge-field frame where the embedding tensor takes the value in the sum of the {\bf 36} and {\bf 36'} representations of ${\rm SL}(8)$, we present a scheme that allows us to derive an analytic expression for the scalar potential. With the help of this formalism, we derive the full potential and gauge coupling functions in analytic forms for the ${\rm SO}(3)\times {\rm SO}(3)$-invariant subsectors of ${\rm SO}(4,4)$ and ${\rm SO}(5,3)$ gaugings, and argue that there exist no new critical points in addition to those discovered so far. For the ${\rm SO}(4,4)$ gauging, we also study the behavior of 6-dimensional scalar fields in this sector near the Dall'Agata-Inverso de Sitter critical point at which the negative eigenvalue of the scalar mass square with the largest modulus goes to zero as the deformation parameter $s$ approaches a critical value $s_{\rm c}$. We find that when the deformation parameter $s$ is taken sufficiently close to the critical value, arbitrarily long inflation can be realized without a fine tuning of the initial point. It turns out that the spectral index $n_s$ of the curvature perturbation at the time of the 60 e-folding number is always about $0.96$ and within the $1\sigma$ range $n_s=0.9639\pm0.0047$ obtained by Planck, irrespective of the value of the $\eta$ parameter at the critical saddle point. The tensor-scalar ratio predicted by this model is around $10^{-3}$ and is close to the value in the Starobinsky model.
The idea that the vacuum energy density $\rho_{\Lambda}$ could be time dependent is a most reasonable one in the expanding Universe; in fact, much more reasonable than just a rigid cosmological constant for the entire cosmic history. Being $\rho_{\Lambda}=\rho_{\Lambda}(t)$ dynamical, it offers a possibility to tackle the cosmological constant problem in its various facets. Furthermore, for a long time (most prominently since Dirac's first proposal on a time variable gravitational coupling) the possibility that the fundamental "constants" of Nature are slowly drifting with the cosmic expansion has been continuously investigated. In the last two decades, and specially in recent times, mounting experimental evidence attests that this could be the case. In this paper, we consider the possibility that these two groups of facts might be intimately connected, namely that the observed acceleration of the Universe and the possible time variation of the fundamental constants are two manifestations of the same underlying dynamics. We call it: the "micro and macro connection", and on its basis we expect that the cosmological term in Einstein's equations, Newton's coupling and the masses of all the particles in the Universe, both the dark matter particles and the ordinary baryons and leptons, should all drift with the cosmic expansion. Here we discuss specific cosmological models realizing such possibility in a way that preserves the principle of covariance of General Relativity.
Electron density profiles of the Venus' ionosphere are inverted from the Venus Express (VEX) one-way open-loop radio occultation experiments carried out by Shanghai 25 m antenna from November 2011 to January 2012 at solar maximum conditions and by New Norcia 35 m antenna from August 2006 to June 2008 at solar intermediate conditions. The electron density profile (from 110 km to 400 km) retrieved from the X-band egress observation at Shanghai station, shows a single peak near 147 km with a peak density of about $2 \times 10^4 \rm{cm}^{-3}$ at a solar zenith angle of 94$^{\circ}$. As a comparison, the VEX radio science (VeRa) observations at New Norcia station were also examined, including S-, X-band and dual-frequency data in the ingress mode. The results show that the electron density profiles retrieved from the S-band data are more analogous to the dual-frequency data in the profile shape, compared with the X-band data. Generally, the S-band results slightly underestimate the magnitude of the peak density, while the X-band results overestimate that. The discrepancy in the X-band profile is probably due to the relatively larger unmodeled orbital errors. It is also expected that the ionopause height is sensitive to the solar wind dynamical pressure in high and intermediate solar activities, usually in the range of 200 km - 1000 km on the dayside and much higher on the nightside. Structural variations ("bulges" and fluctuations) can be found in the electron density profiles in intermediate solar activity, which may be caused by the interaction of the solar wind with the ionosphere. Considerable ionizations can be observed in the Venus' nightside ionosphere, which are unexpected for the Martian nightside ionosphere in most cases.
The main goal of the present contribution is a pedagogical introduction to the fascinating world of neutron stars by relying on relativistic density functional theory. Density functional theory provides a powerful--and perhaps unique--framework for the calculation of both the properties of finite nuclei and neutron stars. Given the enormous densities that may be reached in the core of neutron stars, it is essential that such theoretical framework incorporates from the outset the basic principles of Lorentz covariance and special relativity. After a brief historical perspective, we present the necessary details required to compute the equation of state of dense, neutron-rich matter. As the equation of state is all that is needed to compute the structure of neutron stars, we discuss how nuclear physics--particularly certain kind of laboratory experiments--can provide significant constrains on the behavior of neutron-rich matter.
The value of the safety factor on the magnetic axis of a finite-beta spheromak is shown to be a function of beta in contrast to what was used in P. M. Bellan, Phys. Plasmas 9, 3050 (2002); this dependence on beta substantially reduces the gradient of the safety factor compared to the previous calculation. The method for generating finite-beta spheromak equilibria is extended to generate equilibria describing toroidal magnetic "bubbles" where the hydrodynamic pressure on the magnetic axis is less than on the toroid surface. This "anti-confinement" configuration can be considered an equilibrium with an inverted beta profile and is relevant to interplanetary magnetic clouds as these clouds have lower hydrodynamic pressure in their interior than on their surface.
The shadows cast by non-rotating and rotating modified gravity (MOG) black holes are determined by the two parameters mass $M$ and angular momentum $J=Ma$. The sizes of the shadows cast by the spherically symmetric static Schwarzschild-MOG and Kerr-MOG rotating black holes increase significantly as the free parameter $\alpha$ is increased from zero. The Event Horizon Telescope (EHT) shadow image measurements can determine whether Einstein's general relativity is correct or whether it should be modified in the presence of strong gravitational fields.
Links to: arXiv, form interface, find, astro-ph, recent, 1502, contact, help (Access key information)