The chromospheric activity index logR'HK of stars hosting transiting hot Jupiters appears to be correlated with the planets' surface gravity. One of the possible explanations is based on the presence of condensations of planetary evaporated material located in a circumstellar cloud that absorbs the CaII H&K and MgII h&k resonance line emission flux, used to measure chromospheric activity. A larger column density in the condensations, or equivalently a stronger absorption in the chromospheric lines, is obtained when the evaporation rate of the planet is larger, which occurs for a lower gravity of the planet. We analyze here a sample of stars hosting transiting hot Jupiters tuned in order to minimize systematic effects (e.g., interstellar medium absorption). Using a mixture model, we find that the data are best fit by a two-linear-regression model. We interpret this result in terms of the Vaughan-Preston gap. We use a Monte Carlo approach to best take into account the uncertainties, finding that the two intercepts fit the observed peaks of the distribution of logR'HK for main-sequence solar-like stars. We also find that the intercepts are correlated with the slopes, as predicted by the model based on the condensations of planetary evaporated material. Our findings bring further support to this model, although we cannot firmly exclude different explanations. A precise determination of the slopes of the two linear components would allow one to estimate the average effective stellar flux powering planetary evaporation, which can then be used for theoretical population and evolution studies of close-in planets.
Massive black hole binaries (MBHBs) are a natural byproduct of galaxy mergers. Previous studies have shown that flares from stellar tidal disruption events (TDEs) are modified by the presence of a secondary perturber, causing interruptions in the light curve. We study the dynamics of TDE debris in the presence of a milliparsec-separated MBHB by integrating ballistic particle orbits in the time-varying potential of the binary. We find that gaps in the light curve appear when material misses the accretion radius on its first return to pericentre. Subsequent recurrences can be decomposed into "continuous" and "delayed" components, which exhibit different behaviour. We find that this potential can substantially alter the locations of stream self-intersections. When debris is confined to the plane, we find that close encounters with the secondary BH leave noticeable signatures on the fallback rate and can result in significant accretion onto the secondary BH. Tight, equal-mass MBHBs accrete equally, periodically trading the infalling stream.
We present a detailed analysis of a red quasar at z=2.32 with an intervening damped Lyman-alpha absorber (DLA) at z=2.13. Using high quality data from the X-shooter spectrograph at ESO Very Large Telescope we find that the absorber has a metallicity consistent with Solar. We observe strong C I and H$_2$ absorption indicating a cold, dense absorbing medium. Partial coverage effects are observed in the C I lines, from which we infer a covering fraction of $27 \pm 6$ % and a physical diameter of the cloud of 0.1 pc. From the covering fraction and size, we estimate the size of the background quasar's broad line region. We search for emission from the DLA counterpart in optical and near-infrared imaging. No emission is observed in the optical data. However, we see tentative evidence for a counterpart in the H and K' band images. The DLA shows high depletion (as probed by [Fe/Zn]=-1.22) indicating that significant amounts of dust must be present in the DLA. By fitting the spectrum with various dust reddened quasar templates we find a best-fitting amount of dust in the DLA of $A(V)_{\rm DLA}=0.28 \pm 0.01|_{\rm stat} \pm 0.07|_{\rm sys}$. We conclude that dust in the DLA is causing the colours of this intrinsically very luminous background quasar to appear much redder than average quasars, thereby not fulfilling the criteria for quasar identification in the Sloan Digital Sky Survey. Such chemically enriched and dusty absorbers are thus underrepresented in current samples of DLAs.
The detection of optical surface brightness structures in the sky with magnitudes fainter than 30 mag/arcsec^2 (3sigma in 10x10 arcsec boxes; r-band) has remained elusive in current photometric deep surveys. Here we show how present-day 10 meter class telescopes can provide broadband imaging 1.5-2 mag deeper than most previous results within a reasonable amount of time (i.e. <10h on source integration). In particular, we illustrate the ability of the 10.4 m Gran Telescopio de Canarias (GTC) telescope to produce imaging with a limiting surface brightness of 31.5 mag/arcsec^2 (3sigma in 10x10 arcsec boxes; r-band) using 8.1 hours on source. We apply this power to explore the stellar halo of the galaxy UGC00180, a galaxy analogous to M31 located at ~150 Mpc, by obtaining a surface brightness radial profile down to mu_r~33 mag/arcsec^2. This depth is similar to that obtained using star counts techniques of Local Group galaxies, but is achieved at a distance where this technique is unfeasible. We find that the mass of the stellar halo of this galaxy is ~4x10^9 Msun, i.e. 3+-1% of the total stellar mass of the whole system. This amount of mass in the stellar halo is in agreement with current theoretical expectations for galaxies of this kind.
At present, all physical models of diffuse Galactic gamma-ray emission assume that the distribution of cosmic-ray sources traces the observed populations of either OB stars, pulsars, or supernova remnants. However, since H2-rich regions host significant star formation and numerous supernova remnants, the morphology of observed H2 gas should also provide a physically motivated, high-resolution tracer for cosmic-ray injection. We assess the impact of utilizing H2 as a tracer for cosmic-ray injection on models of diffuse Galactic gamma-ray emission. We employ state-of-the-art 3D particle diffusion and gas density models, along with a physical model for the star-formation rate based on global Schmidt laws. Allowing a fraction, f_H2, of cosmic-ray sources to trace the observed H2 density, we find that a theoretically well-motivated value f_H2 ~ 0.20 -- 0.25 (i) provides a significantly better global fit to the diffuse Galactic gamma-ray sky and (ii) highly suppresses the intensity of the residual gamma-ray emission from the Galactic center region. Specifically, in models utilizing our best global fit values of f_H2 ~ 0.20 -- 0.25, the spectrum of the galactic center gamma-ray excess is drastically affected, and the morphology of the excess becomes inconsistent with predictions for dark matter annihilation.
Using our cosmological radiative transfer code, we study the implications of the updated QSO emissivity and star formation history for the escape fraction (f_esc) of hydrogen ionizing photons from galaxies. We estimate the f_esc that is required to reionize the Universe and to maintain the ionization state of the intergalactic medium in the post-reionization era. At z>5.5, we show that a constant f_esc of 0.14 to 0.22 is sufficient to reionize the Universe. At z<3.5, consistent with various observations, we find that f_esc can have values from 0 to 0.05. However, a steep rise in f_esc, of at least a factor of ~3, is required between z=3.5 to 5.5. It results from a rapidly decreasing QSO emissivity at z>3 together with a nearly constant measured H I photoionization rates at 3<z<5. We show that, this requirement of a steep rise in f_esc over a very short time can be relaxed if we consider the contribution from a recently found large number density of faint QSOs at z>4. In addition, a simple extrapolation of the contribution of such QSOs to high-z suggests that QSOs alone can reionize the Universe. This implies, at z>3.5, that either the properties of galaxies should evolve rapidly to increase the f_esc or most of the low mass galaxies should host massive blackholes and sustain accretion over a prolonged period. These results motivate a careful investigation of theoretical predictions of these alternate scenarios that can be distinguished using future observations. Moreover, it is also very important to revisit the measurements of H I photoionization rates that are crucial to the analysis presented here.
We examine the evolution of the relation between stellar mass surface density, velocity dispersion and half-light radius-the stellar mass fundamental plane-for quiescent galaxies at z<0.7. We measure the local relation from galaxies in the Sloan Digital Sky Survey and the intermediate redshift relation from ~500 quiescent galaxies with stellar masses 10 < log(M_stellar/M_solar) < 11.5. Nearly half of the quiescent galaxies in our intermediate redshift sample are compact. After accounting for important selection and systematic effects, the size and velocity dispersion distributions of galaxies at intermediate redshifts are similar to galaxies in the local universe. The orientation and zero-point of the stellar mass fundamental plane is independent of redshift for massive quiescent galaxies at z<0.7. Compact quiescent galaxies fall on the same relation as the extended objects. We confirm that compact quiescent galaxies are the tail of the size and mass distribution of the normal quiescent galaxy population.
A method for the automatic mapping of coronal holes (CH) using simultaneous
multi-instrument EUV imaging data is described. Synchronized EUV images from
STEREO/EUVI A\&B 195\AA\ and SDO/AIA 193\AA\ undergo preprocessing steps that
include PSF-deconvolution and the application of data-derived intensity
corrections that account for center-to-limb variations (limb-brightening) and
inter-instrument intensity normalization. A systematic approach is taken to
derive a robust limb-brightening correction technique that takes advantage of
unbiased long-term averages of data and respects the physical nature of the
problem. The new preprocessing greatly assists in CH detection, allowing for
the use of a simplified variable-connectivity two-threshold region growing
image segmentation algorithm to obtain consistent detection results. Some
examples of the generated synchronic EUV and CH maps are shown, as well as
preliminary analysis of CH evolutions.
Several data and code products are made available to the community ({\tt
www.predsci.com/chd}): For the period of this study (06/10/2010 to 08/18/14) we
provide synchronic EUV and coronal hole map data at 6-hour cadence, data
derived limb brightening corrections for STEREO/EUVI A\&B 195\AA\ and SDO/AIA
193\AA, and inter-instrument correction factors to equate their intensities. We
also provide the coronal hole image segmentation code modules ({\tt ezseg})
which are implemented in both FORTRAN OpenMP and GPU-accelerated C-CUDA. A
complete implementation of our coronal hole detection pipeline in the form of a
ready-to-use MATLAB driver script {\tt euv2chm} utilizing {\tt ezseg} is also
made available.
Dramatically improved data from observatories like the CoRoT and Kepler spacecraft have recently facilitated nonlinear time series analysis and phenomenological modeling of variable stars, including the search for strange (aka fractal) or chaotic dynamics. We recently argued [Lindner et al., Phys. Rev. Lett. 114 (2015) 054101] that the Kepler data includes "golden" stars, whose luminosities vary quasiperiodically with two frequencies nearly in the golden ratio, and whose secondary frequencies exhibit power-law scaling with exponent near -1.5, suggesting strange nonchaotic dynamics and singular spectra. Here we use a series of phenomenological models to make plausible the connection between golden stars and fractal spectra. We thereby suggest that at least some features of variable star dynamics reflect universal nonlinear phenomena common to even simple systems.
We apply the Jeans Axisymmetric Multi-Gaussian Expansion method to the stellar kinematic maps of 40 Sa-Sd EDGE-CALIFA galaxies and derive their circular velocity curves (CVCs). The CVCs are classified using the Dynamical Classification method developed in Kalinova et al. (2015) . We also calculate the observational baryon efficiency, OBE, where $M_*/M_b=M_*/(M_*+M_{HI}+M_{H_2})$ of the galaxies using their stellar mass, total neutral hydrogen mass and total molecular gas from CO luminosities. Slow-rising, Flat and Round-peaked CVC types correspond to specific OBEs, stellar and dark matter (DM) halo mass values, while the Sharp-peaked CVCs span in the whole DM halo mass range of $10^{11}-10^{14} M_{\odot}$.
The compressibility of molecular cloud (MC) turbulence plays a crucial role in star formation models, because it controls the amplitude and distribution of density fluctuations. The relation between the compressive ratio (the ratio of powers in compressive and solenoidal motions) and the statistics of turbulence has been studied systematically only in idealized simulations with random external forces. In this work, we analyze a simulation of large-scale turbulence(250 pc) driven by supernova (SN) explosions that has been shown to yield realistic MC properties. We demonstrate that SN driving results in MC turbulence that is only mildly compressive, with the turbulent ratio of compressive to solenoidal modes ~0.3 on average, lower than the equilibrium value of 0.5 found in the inertial range of isothermal simulations with random solenoidal driving. We also find that the compressibility of the turbulence is not noticeably affected by gravity, nor is the mean cloud expansion or contraction velocity (MCs do not collapse as a whole even if their own prestellar cores collapse to form stars). Furthermore, the clouds follow the same relation between the rms density and the rms velocity as in isothermal turbulence and their average gas density PDF is described well by a lognormal distribution, with the addition of a high-density power-law tail when self-gravity is included.
Motivated by recent discussions, both in private and in the literature, we use a Monte Carlo simulation of planetary systems to investigate sources of bias in determining the mass-radius distribution of exoplanets for the two primary techniques used to measure planetary masses---Radial Velocities (RVs) and Transit Timing Variations (TTVs). We assert that mass measurements derived from these two methods are comparably reliable---as the physics underlying their respective signals is well understood. Nevertheless, their sensitivity to planet mass varies with the properties of the planets themselves. We find that for a given planet size, the RV method tends to find planets with higher mass while the sensitivity of TTVs is more uniform. This ``sensitivity bias'' implies that a complete census of TTV systems is likely to yield a more robust estimate of the mass-radius distribution provided there are not important physical differences between planets near and far from mean-motion resonance. We discuss differences in the sensitivity of the two methods with orbital period and system architecture, which may compound the discrepancies between them (e.g., short period planets detectable by RVs may be more dense due to atmospheric loss). We advocate for continued mass measurements using both approaches as a means both to measure the masses of more planets and to identify potential differences in planet structure that may result from their dynamical and environmental histories.
We present predictions of the Si iv ions in turbulent mixing layers (TMLs) between hot and cool gas and in cool high-velocity clouds (HVCs) that travel through a hot halo, complementing the C iv, N v, and O vi predictions in Kwak & Shelton, Kwak et al., and Henley et al. We find that the Si iv ions are most abundant in regions where the hot and cool gases first begin to mix or where the mixed gas has cooled significantly. The predicted column densities of high velocity Si iv and the predicted ratios of Si iv to C iv and O vi found on individual sightlines in our HVC simulations are in good agreement with observations of high velocity gas. Low velocity Si iv is also seen in the simulations, as a result of decelerated gas in the case of the HVC simulations and when looking along directions that pass perpendicular to the direction of motion in the TML simulations. The ratios of low velocity Si iv to C iv and O vi in the TML simulations are in good agreement with those recorded for Milky Way halo gas, while the ratio of Si iv to O vi from the decelerated gas in the HVC simulations is lower than that observed at normal velocity in the Milky Way halo. We attribute the shortfall of normal velocity Si iv to not having modeled the effects of photoionization and, following Henley et al., consider a composite model that includes decelerated HVC gas, supernova remnants, galactic fountain gas, and the effect of photoionization.
In this second paper in the series, we investigate occurrence frequencies of apparent unipolar processes, cancellation, and emergence of patch structures in quiet regions. Apparent unipolar events are considerably more frequent than cancellation and emergence as per our definition, which is consistent with Lamb et al. (2013). Furthermore, we investigate the frequency distributions of changes in flux during apparent unipolar processes are and found that they concentrate around the detection limit of the analysis. Combining these findings with the results of our previous paper, Iida et al. (2012), that merging and splitting are more dominant than emergence and cancellation, these results support the understanding that apparent unipolar processes are actually interactions with and among patches below the detection limit and that there still are numerous flux interactions between the flux range in this analysis and below the detection limit. We also investigate occurrence frequency distributions of flux decrease during cancellation. We found a relatively strong dependence, 2.48$\pm$0:26 as a power-law index. This strong dependence on flux is consistent with the model, which is suggested in the previous paper.
Photodetachment cross-section $\sigma_{ph}(p_e)$ of the negatively charged hydrogen ion H$^{-}$ is determined with the use of highly accurate variational wave functions constructed for this ion. Photodetachment cross-sections of the H$^{-}$ ion are also studied for very small and very large values of the photo-electron momentum $p_e$. Maximum of this cross-section has been evaluated to very high accuracy and we have found that $[\sigma_{ph}(p_e)]_{\max} \approx$ 3.8627035742 $\cdot 10^{-17}$ $cm^2$ at $p_e \approx$ 0.113206(1) $a.u.$ Photodetachment of the H$^{-}$ ion at very small and very large $p_e$ values is also considered. Our method is based upon the Rayleigh's formula for spherical Bessel functions.
Since NGC300 is a bulge-less, isolated low-mass galaxy and has not experienced radial migration during its evolution history, it can be treated as an ideal laboratory to test simple galactic chemical evolution models. By assuming its disk forms gradually from continuous accretion of primordial gas and including the gas-outflow process, we construct a simple chemical evolution model for NGC300 to build a bridge between its SFH and its observed data, especially the present-day radial profiles and global observed properties (e.g., cold gas mass, star-formation rate and metallicity). By means of comparing the model predictions with the corresponding observations, we adopt the classical $\chi^{2}$ methodology to find out the best combination of free parameters $a$, $b$ and $b_{\rm out}$. Our results show that, by assuming an inside-out formation scenario and an appropriate outflow rate, our model reproduces well most of the present-day observational values, not only the radial profiles but also the global observational data for the NGC300 disk. Our results suggest that NGC300 may experience a rapid growth of its disk. Through comparing the best-fitting model predicted SFH of NGC300 with that of M33, we find that the mean stellar age of NGC300 is older than that of M33 and there is a lack of primordial gas infall onto the disk of NGC300 recently. Our results also imply that the local environment may paly a key role in the secular evolution of NGC300.
The integrated Sachs-Wolfe (ISW) effect and its non-linear extension Rees-Sciama (RS) effect provide us the information of the time evolution of gravitational potential. The cross-correlation between the cosmic microwave background (CMB) and the large scale structure (LSS) is known as a promising way to extract the ISW (RS) effect. It is known that the RS effect shows the unique behavior by changing the anti-correlated cross correlation between the CMB and the mass tracer into the positively correlated cross correlation compared to the linear ISW effect. We show that the dependence of this flipping scale of the cross-correlation between RS and weak lensing on dark energy models. However, there exists the degeneracy between DE and $\Omega_{\rm{m}0}$ which might be broken by redshift dependent observables. The cross-correlation between the momentum field and the density field might be served as the better observable to be used for this purpose.
The anisotropic 2-point correlation function (2PCF) of galaxies measures pairwise clustering as a function of the pair separation's angle to the line of sight. The latter is often defined as either the angle bisector of the observer-galaxy-pair triangle or the vector from the observer to the separation midpoint. Here we show how to accelerate either of these measurements with Fourier Transforms, using a slight generalization of the Yamamoto et al. (2006) estimator in which each member of the pair is used successively as the line of sight. We also present perturbation theory predictions for our generalized estimator including wide-angle corrections.
We report the discovery of an extremely dense group of massive galaxies at the centre of the protocluster at $z=3.09$ in the SSA22 field from near-infrared spectroscopy conducted with the Multi-Object InfraRed Camera and Spectrograph (MOIRCS) equipped on the Subaru Telecope. The newly discovered group comprises seven galaxies confirmed at $z_{\rm spec}\approx3.09$ within 180 kpc including five massive objects with the stellar masses larger than $10^{10.5}~M_{\odot}$ and is associated with a bright sub-mm source SSA22-AzTEC14. The dynamical mass of the group estimated from the line-of-sight velocity dispersion of the members is $M_{\rm dyn}\sim1.6\pm0.3\times10^{13}~M_{\odot}$. Such a dense group is expected to be very rare at high redshift as we found only a few comparable systems in large-volume cosmological simulations. Such rare groups in the simulations are hosted in collapsed halos with $M_{\rm vir}=10^{13.4}-10^{14.0}~M_{\odot}$ and evolve into the brightest cluster galaxies (BCGs) of the most massive clusters at present. The observed AzTEC14 group at $z=3.09$ is therefore very likely to be a proto-BCG in the multiple merger phase. The observed total stellar mass of the group is $5.8^{+5.1}_{-2.0}\times10^{11}~M_{\odot}$. It suggests that over half the stellar mass of its descendant had been formed by $z=3$. Moreover, we identified over two members for each of the four Ly$\alpha$ blobs (LABs) using our new spectroscopic data. This verifies our previous argument that many of the LABs in the SSA22 protocluster associated with multiple developed stellar components.
We analyzed observations of interstellar neutral helium (ISN~He) obtained from the Interstellar Boundary Explorer (IBEX) satellite during its first six years of operation. We used a refined version of the ISN~He simulation model, presented in the companion paper by Sokol_et al. 2015, and a sophisticated data correlation and uncertainty system and parameter fitting method, described in the companion paper by Swaczyna et al 2015. We analyzed the entire data set together and the yearly subsets, and found the temperature and velocity vector of ISN~He in front of the heliosphere. As seen in the previous studies, the allowable parameters are highly correlated and form a four-dimensional tube in the parameter space. The inflow longitudes obtained from the yearly data subsets show a spread of ~6 degree, with the other parameters varying accordingly along the parameter tube, and the minimum chi-square value is larger than expected. We found, however, that the Mach number of the ISN~He flow shows very little scatter and is thus very tightly constrained. It is in excellent agreement with the original analysis of ISN~He observations from IBEX and recent reanalyses of observations from Ulysses. We identify a possible inaccuracy in the Warm Breeze parameters as the likely cause of the scatter in the ISN~He parameters obtained from the yearly subsets, and we suppose that another component may exist in the signal, or a process that is not accounted for in the current physical model of ISN~He in front of the heliosphere. From our analysis, the inflow velocity vector, temperature, and Mach number of the flow are equal to lambda_ISNHe = 255.8 +/- 0.5 degree, beta_ISNHe = 5.16 +/- 0.10 degree, T_ISNHe = 7440 +/- 260 K, v_ISNHe = 25.8 +/- 0.4$ km/s, and M_ISNHe = 5.079 +/- 0.028, with uncertainties strongly correlated along the parameter tube.
The C II 133.5 nm multiplet has been observed by NASA's Interface Region Imaging Spectrograph (IRIS) in unprecedented spatial resolution. The aims of this work are to characterize these new observations of the C II lines, place them in context with previous work, and to identify any additional value the C II lines bring when compared with other spectral lines. We make use of wide, long exposure IRIS rasters covering the quiet Sun and an active region. Line properties such as velocity shift and width are extracted from individual spectra and analyzed. The lines have a variety of shapes (mostly single-peak or double-peak), are strongest in active regions and weaker in the quiet Sun. The ratio between the 133.4 nm and 133.5 nm components is always less than 1.8, indicating that their radiation is optically thick in all locations. Maps of the C II line widths are a powerful new diagnostic of chromospheric structures, and their line shifts are a robust velocity diagnostic. Compared with earlier quiet Sun observations, we find similar absolute intensities and mean line widths, but smaller red shifts; this difference can perhaps be attributed to differences in spectral resolution and spatial coverage. The C II intensity maps are somewhat similar to those of transition region lines, but also share some features with chromospheric maps such as those from the Mg II k line, indicating that they are formed between the upper chromosphere and transition region. C II intensity, width, and velocity maps can therefore be used to gather additional information about the upper chromosphere.
We have carried out a study of temperature evolution during thermonuclear bursts in LMXBs using broad band data from two instruments onboard BeppoSAX, the MECS and the PDS. However, instead of applying the standard technique of time resolved spectroscopy, we have determined the temperature in small time intervals using the ratio of count rates in the two instruments assuming a blackbody nature of burst emission and different interstellar absorption for different sources. Data from a total of twelve observations of six sources were analysed during which 22 bursts were detected. We have obtained temperatures as high as ~3.0 keV, even when there is no evidence of photospheric radius expansion. These high temperatures were observed in the sources within different broadband spectral states (soft and hard).
We study preheating after hilltop inflation where the inflaton couples to another scalar field, e.g. a right-handed sneutrino, which provides a mechanism for generating the correct initial conditions for inflation and also a decay channel for the inflaton that allows for reheating and non-thermal leptogenesis. In the presence of such a coupling, we find that after the known phases of preheating during which the inflaton field becomes fully inhomogeneous, there can be a subsequent preheating phase where the fluctuations of the other field get resonantly enhanced, from initial vacuum fluctuations up to amplitudes of the same order (and even larger) as the ones of the inflaton field. This resonant enhancement differs from the usual parametric resonance as the inflaton field is highly inhomogeneous at the time the enhancement takes place. We study this effect using lattice simulations as well as semi-analytically with a generalized Floquet analysis for inhomogeneous background fields.
We have developed a refined and optimized version of the Warsaw Test Particle Model of interstellar neutral gas in the heliosphere, specially tailored for analysis of IBEX-Lo observations. The former version of the model was used in the analysis of neutral He observed by IBEX that resulted in an unexpected conclusion that the interstellar neutral He flow vector was different than previously thought and that a new population of neutral He, dubbed the Warm Breeze, exists in the heliosphere. It was also used in the reanalysis of Ulysses observations that confirmed the original findings on the flow vector, but suggested a significantly higher temperature. The present version model has two strains targeted for different applications, based on an identical paradigm, but differing in the implementation and in the treatment of ionization losses. We present the model in detail and discuss numerous effects related to the measurement process that potentially modify the resulting flux of ISN~He observed by IBEX, and identify those of them that should not be omitted in the simulations to avoid biasing the results. This paper is part of a coordinated series of papers presenting the current state of analysis of IBEX-Lo observations of ISN~He. Details of the analysis method are presented by Swaczyna et al. 2015, and results of the analysis are presented by Bzowski et al. 2015.
Neutral interstellar helium has been observed by the Interstellar Boundary Explorer (IBEX) since 2009 with a signal-to-noise ratio well above 1000. Because of the geometry of the observations, the signal observed from January to March each year is the easiest to identify. However, as we show via simulations, the portion of the signal in the range of intensities from 10^{-3} to 10^{-2} of the peak value, previously mostly left out from the analysis, may bring important information about the details of the distribution function of interstellar He gas in front of the heliosphere. In particular, these observations may inform us about possible departures of the parent interstellar He population from equilibrium. We compare the expected distribution of the signal for the canonical assumption of a single Maxwell-Boltzmann population with the distributions for a superposition of the Maxwell-Boltzmann primary population and the recently discovered Warm Breeze, and for a single primary population given by a kappa function. We identify the regions on the sky where the differences between those cases are expected to be the most visible against the background. We discuss the diagnostic potential of the fall peak of the interstellar signal, reduced by a factor of 50 due to the Compton-Getting effect but still above the detection limit of IBEX. We point out the strong energy dependence of the fall signal and suggest that searching for this signal in the data could bring an independent assessment of the low-energy measurement threshold of the IBEX-Lo sensor.
A tidal disruption event (TDE) occurs when a star wanders close enough to a black hole to be disrupted by its tidal force. The debris of a tidally disrupted star are expected to form an accretion disc around the supermassive black hole. The light curves of these events sometimes show a quasi-periodic modulation of the flux that can be associated with the precession of the accretion disc due to the Lense-Thirring ("frame-dragging") effect. Since the initial star orbit is in general inclined with respect to the black hole spin, this misalignment combined with the Lense-Thirring effect leads to a warp in the disc. In this paper we provide a simple model of the system composed by a thick and narrow accretion disc surrounding a spinning supermassive black hole, with the aim to: (a) compute the expected precession period as a function of the system parameters, (b) discuss the conditions that have to be satisfied in order to have rigid precession, (c) investigate the alignment process, highlighting how different mechanisms play a role leading the disc and the black hole angular momenta into alignment.
A code used previously to predict O-star mass fluxes as a function of metallicity is used to compute a grid of models with the metallicity of the Small Magellanic Cloud (SMC). These models allow mass-loss rates to be derived by interpolation for all O-type supergiants in the SMC, with the possible exception of extremely massive stars close to the Eddington limit.
This paper combines observational datasets and cosmological simulations to generate realistic numerical replicas of the nearby Universe. These latter are excellent laboratories for studies of the non-linear process of structure formation in our neighborhood. With measurements of radial peculiar velocities in the Local Universe (cosmicflows-2) and a newly developed technique, we produce Constrained Local UniversE Simulations (CLUES). To assess the quality of these constrained simulations, we compare them with random simulations as well as with local observations. The cosmic variance, defined as the mean one-sigma scatter of cell-to-cell comparison between two fields, is significantly smaller for the constrained simulations than for the random simulations. Within the inner part of the box where most of the constraints are, the scatter is smaller by a factor 2 to 3 on a 5 Mpc/h scale with respect to that found for random simulations. This one-sigma scatter obtained when comparing the simulated and the observation-reconstructed velocity fields is only 104 +/- 4 km/s i.e. the linear theory threshold. These two results demonstrate that these simulations are in agreement with each other and with the observations of our neighborhood. For the first time, simulations constrained with observational radial peculiar velocities resemble the Local Universe up to a distance of 150 Mpc/h on a scale of a few tens of megaparsecs. When focusing on the inner part of the box, the resemblance with our cosmic neighborhood extends to a few megaparsecs (< 5 Mpc/h). The simulations provide a proper Large Scale environment for studies of the formation of nearby objects.
We report the discovery of a population of deeply embedded protostellar candidates in the 20 km s$^{-1}$ cloud, one of the massive molecular clouds in the Central Molecular Zone (CMZ) of the Milky Way, using interferometric submillimeter continuum and H$_2$O maser observations. The submillimeter continuum emission shows five 1-pc scale clumps, each of which further fragments into several 0.1-pc scale cores. We identify 17 dense cores, among which 12 are gravitationally bound. Among the 18 H$_2$O masers detected, 13 coincide with the cores and probably trace outflows emanating from the protostars. There are also 5 gravitationally bound dense cores without H$_2$O maser detection. In total the 13 masers and 5 cores may represent 18 protostars with spectral types later than B1 or potential growing more massive stars at earlier evolutionary stage, given the non-detection in the centimeter radio continuum. In combination with previous studies of CH$_3$OH masers, we conclude that the star formation in this cloud is at an early evolutionary phase, before the presence of any significant ionizing or heating sources. Our findings indicate that star formation in this cloud may be triggered by a tidal compression as it approaches pericenter, similar to the case of G0.253+0.016 but with a higher star formation rate, and demonstrate that high angular resolution, high sensitivity maser and submillimeter observations are a promising technique to unveil deeply embedded star formation in the CMZ.
We present a detailed study of the kinematics of M dwarfs in the CARMENES (Calar Alto high-Resolution search for M dwarfs with Exoearths with Near-infrared and optical Echelle Spectrographs) input catalog. We have selected all M dwarfs with known parallactic distance or a good photometric distance estimation, precise proper motion in the literature or as determined by us, and radial velocity measurements. Using these parameters, we computed the M dwarfs galactic space motions (U, V, W). For the stars with U and V velocity components inside or near the boundaries that determine the young disk population, we have analyzed the possible membership in the classical moving groups and nearby loose associations with ages between 10 and 600 Ma. For the candidate members, we have compiled information available in the literature in order to constrain their membership by applying other age-dating methods.
Planet detection through microlensing is usually limited by a well-known degeneracy in the Einstein timescale $t_E$, which prevents mass and distance of the lens to be univocally determined. Gould (2013) has shown that a satellite in geosynchronous orbit could provide masses and distances for most standard planetary events ($t_E \approx 20$ days) via a microlens parallax measurement. This paper extends the analysis to shorter Einstein timescales, $t_E \approx 1$ day, dealing with the case of Jupiter-mass lenses. We then study the capabilities of a low Earth orbit satellite at even shorter timescales, $t_E \approx 0.1$ days. A Fisher matrix analysis is employed to predict how the 1-sigma error on parallax depends on $t_E$ and the peak magnification of the microlensing event. It is shown that a geosynchronous satellite could detect parallaxes for Jupiter-mass free-floaters and discover planetary systems around very low-mass brown dwarfs. Moreover, a low Earth orbit satellite could lead to the discovery of Earth-mass free-floating planets. Limitations to these results can be the strong requirements on the photometry, the effects of blending and, in case of the low orbit, the Earth umbra.
We present a qualitative analysis of the variability of quasar broad absorption lines using the large multi-epoch spectroscopic dataset of the Sloan Digital Sky Survey Data Release 10. We confirm that variations of absorption lines are highly coordinated among different components of the same ion or the same absorption component of different ions for C IV, Si IV and N V. Furthermore, we show that the equivalent widths of the lines decrease or increase statistically when the continuum brightens or dims. This is further supported by the synchronized variations of emission and absorption line equivalent width, when the well established intrinsic Baldwin effect for emission lines is taken into account. We find that the emergence of an absorption component is usually accompanying with dimming of the continuum while the disappearance of an absorption line component with brightening of the continuum. This suggests that the emergence or disappearance of a C IV absorption component is only the extreme case, when the ionic column density is very sensitive to continuum variations or the continuum variability amplitude is larger. These results support the idea that absorption line variability is driven mainly by changes in the gas ionization in response to continuum variations, that the line-absorbing gas is highly ionized, and in some extreme cases, too highly ionized to be detected in UV absorption lines. Due to uncertainties in the spectroscopic flux calibration, we cannot quantify the fraction of quasars with asynchronized continuum and absorption line variations.
We report on an experimental and theoretical investigation of the importance of anharmonicity in the 3 micron CH stretching region of Polycyclic Aromatic Hydrocarbon (PAH) molecules. We present mass-resolved, high-resolution spectra of the gas-phase cold (~4K) linear PAH molecules naphthalene, anthracene, and tetracene. The measured IR spectra show a surprisingly high number of strong vibrational bands. For naphthalene, the observed bands are well separated and limited by the rotational contour, revealing the band symmetries. Comparisons are made to the harmonic and anharmonic approaches of the widely used Gaussian software. We also present calculated spectra of these acenes using the computational program SPECTRO, providing anharmonic predictions enhanced with a Fermi-resonance treatment that utilises intensity redistribution. We demonstrate that the anharmonicity of the investigated acenes is strong, dominated by Fermi resonances between the fundamental and double combination modes, with triple combination bands as possible candidates to resolve remaining discrepancies. The anharmonic spectra as calculated with SPECTRO lead to predictions of the main modes that fall within 0.5% of the experimental frequencies. The implications for the Aromatic Infrared Bands, specifically the 3 micron band are discussed.
We investigate the impact of coulomb screening on primordial nucleosynthesis in a universe having scale factor that evolves linearly with time. Coulomb screening affects primordial nucleosynthesis via enhancement of thermonuclear reaction rates. This enhancement is determined by the solving Poisson equation within the context of mean field theory (under appropriate conditions during the primordial nucleosynthesis). Using these results, we claim that the mean field estimates of coulomb screening hardly affect the predicted element abundances and nucleosynthesis parameters$, \{\eta_9,\xi_e\}$. The deviations from mean field estimates are also studied in detail by boosting genuine screening results with the screening parameter ($\omega_s$). These deviations show negligible effect on the element abundances and on nucleosynthesis parameters. This work thus rules out the coulomb screening effects on primordial nucleosynthesis in slow evolving models and confirms that constraints in ref.[7] on nucleosynthesis parameters remain unaltered.
A number of simulators have argued that major mergers can sometimes preserve discs (e.g. Springel & Hernquist 2005), but the possibility that they could explain the emergence of lenticular galaxies (S0s) has been generally neglected. In fact, observations of S0s reveal a strong structural coupling between their bulges and discs, which seems difficult to reconcile with the idea that they come from major mergers. However, in Querejeta et al. (2015a) we have used N-body simulations of binary mergers to show that, under favourable conditions, discs are first destroyed but soon regrow out of the leftover debris, matching observational photometric scaling relations (e.g. Laurikainen et al. 2010). Additionally, in Querejeta et al. (2015b) we have shown how the merger scenario agrees with the recent discovery that S0s and most spirals are not compatible in an angular momentum--concentration plane. This important result from CALIFA constitutes a serious objection to the idea that spirals transform into S0s mainly by fading (e.g. via ram-pressure stripping, as that would not explain the observed simultaneous change in $\lambda_\mathrm{Re}$ and concentration), but our simulations of major mergers do explain that mismatch. From such a 3D comparison we conclude that mergers must be a relevant process in the build-up of the current population of S0s.
Diffuse radio emission in galaxy clusters is known to be related to cluster mass and cluster dynamical state. We collect the observed fluxes of radio halos, relics and mini-halos for a sample of galaxy clusters from literature, and calculate their radio powers. We then obtain the values of cluster mass or mass proxies from previous observations, and also obtain the various dynamical parameters of these galaxy clusters from optical and X-ray data. The radio power of relics, halos and mini-halos are correlated with the cluster masses or mass proxies, as found by previous authors, with the correlations concerning giant radio halos being, in general, the strongest ones. We found that the inclusion of dynamical parameters as the third dimension can significantly reduce the data scatter for the scaling relations, especially for radio halos. We therefore conclude that the substructures in X-ray images of galaxy clusters and the irregular distributions of optical brightness of member galaxies can be used to quantitatively characterize the shock waves and turbulence in intracluster medium responsible for reaccelerating particles to generate the observed diffuse radio emission. The power of radio halos and relics are correlated with cluster mass proxies and dynamical parameters in the form of a fundamental plane.
Pulsed sodium laser guide stars (LGS) are useful because they allow for Rayleigh blanking and fratricide avoidance in multiple-LGS systems. Bloch-equation simulations of sodium-light interactions show that these may be able to achieve photon returns nearly equal to, and in some cases greater than, what is seen from continuous-wave (CW) excitation. In this work, we study the time-dependent characteristics of sodium fluorescence, and investigate the optimal format for the new fiber laser LGS that will be part of the upgraded adaptive optics (AO) system on the Shane telescope at Mt. Hamilton. Results of this analysis are examined in the context of their general applicability to other LGS systems and the potential benefits of uplink correction are considered. Comparisons of simulation predictions with measurements from existing LGS are also presented and discussed.
We set to search for Rayleigh scattering and K and Na absorption signatures
from the atmosphere of TrES-3b using ground-based transmission spectroscopy
covering the wavelength range from 530 to 950 nm as observed with OSIRIS@GTC.
Our analysis is based on a Bayesian approach where the light curves covering
a set of given passbands are fitted jointly with PHOENIX-calculated stellar
limb darkening profiles. The analysis is carried out assuming both white and
red -- temporally correlated -- noise, with two approaches (Gaussian processes
and divide-by-white) to account for the red noise.
An initial analysis reveals a transmission spectrum that shows a strong
Rayleigh-like increase in extinction towards the blue end of the spectrum, and
enhanced extinction around the K I resonance doublet near 767 nm. However, the
signal amplitudes are significantly larger than expected from theoretical
considerations. A detailed analysis reveals that the K I-like feature is
entirely due to variability in the telluric O$_2$ absorption, but the
Rayleigh-like feature remains unexplained.
Nucleosynthesis of heavy elements requires the use of different experimental and theoretical methods to determine astrophysical reaction rates than light element nucleosynthesis. Additionally, there are also larger uncertainties involved in the astrophysical models, both because the sites are not well known and because of differing numerical treatments in different models. As an example for the latter, the production of p-nuclei is compared in two different stellar models, demonstrating that a model widely used for postproduction calculations may have a zone grid too coarse to follow the synthesis of p-nuclei in detail.
In the last decade or so there has been debate over the possibility that the fuzzy quantum nature of spacetime might decohere wavefronts emanating from very distant sources. Consequences of that could be "blurred" or "faded" images of compact structures in galaxies, primarily at z>1 for their emitted X-rays and gamma-rays, but perhaps even in ultraviolet through optical light at higher redshift. So far there are only inconclusive hints of this from z~4 active-galactic nucleii and gamma-ray bursts viewed with Fermi and Hubble Space Telescope. If correct though, that would impose a significant, fundamental resolution limit for galaxies out to z~8 in the era of the James Webb Space Telescope and the next generation of ground-based telescopes using adaptive optics.
We explore the hypothesis of having an approximate lepton number conservation as a way to achieve a successful leptogenesis in low-scale seesaw mechanisms. The smallness of the active neutrino masses, as well as a strong degeneracy in the mass spectrum of the heavy sterile states, are both consequence of the assumed approximate symmetry. We propose a minimal extension of the Standard Model in order to implement the idea, and perform an analytical and numerical study to determine the viable solutions in the model and the testability of this leptogenesis scenario in future experiments.
In this paper, we consider rolling tachyon, with steep run-away type of potentials non-minimally coupled to massive neutrino matter. The coupling dynamically builds up at late times as neutrino matter turns non-relativistic. In case of scaling and string inspired potentials, we have shown that non-minimal coupling leads to minimum in the field potential. Given a suitable choice of model parameters, it is shown to give rise to late-time acceleration with the desired equation of state.
We investigate the pre-inflationary dynamics of inflation with the Starobinsky potential, favored by recent data from the Planck mission, using techniques developed to study cosmological perturbations on quantum spacetimes in the framework of loop quantum gravity. We find that for a large part of the initial data, inflation compatible with observations occurs. There exists a subset of this initial data that leads to quantum gravity signatures that are potentially observable. Interestingly, despite the different inflationary dynamics, these quantum gravity corrections to the powerspectra are similar to those obtained for inflation with a quadratic potential, including suppression of power at large scales. Furthermore, for super horizon modes the tensor modes show deviations from the standard inflationary paradigm that are unique to the Starobinsky potential and could be important for non-Gaussian modulation and tensor fossils.
Ultra high energy cosmic rays (UHECRs) probably originate in extreme conditions in which extra dimension effects might be important. In this paper we calculate the correction in black hole accretion mechanisms due to extra dimension effects in the static and rotating cases. A parametrization of the external Kerr horizons in both cases is presented and analysed. We use previous calculations of upper limits on the UHECR flux to set limits on the UHECR production efficiency of nine sources. The upper limit on the UHECR luminosity calculation is based on GeV-TeV gamma-ray measurements. The total luminosity due to the accretion mechanism is compared to the upper limit on UHECRs. The dependence of the UHECR production efficiency upper limit on black hole mass is also presented and discussed.
Links to: arXiv, form interface, find, astro-ph, recent, 1510, contact, help (Access key information)
We study how the hot Universe evolves and acquires the prevailing vacuum state, demonstrating that in specific conditions which are believed to apply, the Universe becomes frozen into the state with the smallest value of Higgs vacuum field $v=\langle h\rangle$, even if this is not the state of lowest energy. This supports the false vacuum dark energy $\Lambda$-model. Under several likely hypotheses we determine the temperature in the evolution of the Universe at which two vacuua $v_1, v_2$ can swap between being true and false. We evaluate the dynamical surface pressure on domain walls between low and high mass vaccua due to the presence of matter and show that the low mass state remains the preferred vacuum of the Universe.
We report on Atacama Large Millimeter/submillimeter Array (ALMA) detections of molecular absorption lines in Bands 3, 6 and 7 toward four radio-loud quasars, which were observed as the bandpass and complex gain calibrators. The absorption systems, three of which are newly detected, are found to be Galactic origin. Moreover, HCO absorption lines toward two objects are detected, which almost doubles the number of HCO absorption samples in the Galactic diffuse medium. In addition, high HCO to H13CO+ column density ratios are found, suggesting that the interstellar media (ISM) observed toward the two calibrators are in photodissociation regions, which observationally illustrates the chemistry of diffuse ISM driven by ultraviolet (UV) radiation. These results demonstrate that calibrators in the ALMA Archive are potential sources for the quest for new absorption systems and for detailed investigation of the nature of the ISM.
Type I X-ray bursts are thermonuclear flashes observed from the surface of accreting neutron stars in Low Mass X-ray Binaries. Oscillations have been observed during the rise and/or decay of some of these X-ray bursts. Those seen during the rise can be well explained by a spreading hot spot model, but to date there haven't been any quantitative studies that consistently track the oscillation amplitude both during the rise and decay (tail) of bursts. Here we compute the light curves and amplitudes of oscillations in X-ray burst models that realistically account for both flame spreading and subsequent cooling. We present results for two such "cooling wake" models, a "canonical" cooling model where each patch on the neutron star surface heats and cools identically, and an "asymmetric" model where parts of the star cool at different rates. We show that while canonical cooling models can generate oscillations in the tails of bursts, they cannot easily produce the highest observed modulation amplitudes. Alternatively, a simple model with asymmetric cooling can achieve higher amplitudes consistent with the observations.
We report on the energy-momentum correlators obtained with recent numerical simulations of the Abelian Higgs model, essential for the computation of cosmic microwave background and matter perturbations of cosmic strings. Due to significant improvements both in raw computing power and in our parallel simulation framework, the dynamical range of the simulations has increased four-fold both in space and time, and for the first time we are able to simulate strings with a constant physical width in both the radiation and matter eras. The new simulations improve the accuracy of the measurements of the correlation functions at the horizon scale and confirm the shape around the peak. The normalization is slightly higher in the high wave-number tails, due to a small increase in the string density. We study for the first time the behaviour of the correlators across cosmological transitions, and discover that the correlation functions evolve adiabatically, ie the network adapts quickly to changes in the expansion rate. We propose a new method for constructing source functions for Einstein-Boltzmann integrators, comparing it with two other methods previously used. The new method is more consistent, easier to implement, and significantly more accurate.
We study the physical properties of a spectroscopic sample of 28 star-forming galaxies in a large filamentary structure in the COSMOS field at $z\sim$0.53, with spectroscopic data taken with the Keck/DEIMOS spectrograph, and compare them with a control sample of 30 field galaxies. We spectroscopically confirm the presence of a large galaxy filament ($\sim$ 8 Mpc), along which five confirmed X-ray groups exist. We show that within the uncertainties, the ionization parameter, equivalent width (EW), EW versus specific star-formation rate (sSFR) relation, EW versus stellar mass relation, line-of-sight velocity dispersion, dynamical mass, and stellar-to-dynamical mass ratio are similar for filament and field star-forming galaxies. However, we show that on average, filament star-forming galaxies are more metal-enriched ($\sim$ 0.1$-$0.15 dex), possibly due to the inflow of the already enriched intrafilamentary gas into filament galaxies. Moreover, we show that electron densities are significantly lower (a factor of $\sim$17) in filament star-forming systems compared to those in the field, possibly because of a longer star-formation timescale for filament star-forming galaxies. Our results highlight the potential pre-processing role of galaxy filaments and intermediate-density environments on the evolution of galaxies, which has been highly underestimated.
We present a novel approach to constrain the formation channels of
Ultra-Compact Dwarf Galaxies (UCDs). This inhomogeneous class of objects of
remnants of tidally stripped dwarf elliptical galaxies and high mass globular
clusters. We use three methods to unravel their nature: 1) we analysed their
surface brightness profiles, 2) we carried out a direct search for tidal
features around UCDs and 3) we compared the spatial distribution of GCs and
UCDs in the halo of their host galaxy.
Based on FORS2 observations, we have studied the detailed structural
composition of a large sample of 97 UCDs in the halo of NGC1399, the central
galaxy of the Fornax cluster, by analysing theirsurface brightness profiles. We
derived the structural parameters of 13 extended UCDs modelling them with a
single Sersic function and decomposing them into composite King and Sersic
profiles. We find evidence for faint stellar envelopes at mu=~26 mag\arcsec^-2
surrounding the UCDs up to an extension of 90pc in radius.
We also show new evidence for faint asymmetric structures and tidal tail-like
features surrounding several of these UCDs, a possible tracer of their origin
and assembly history within their host galaxy halos. In particular, we present
evidence for the first discovery of a significant tidal tail with an extension
of ~350pc around UCD-FORS2.
We searched for local overdensities in the spatial distribution of globular
clusters within the halo of NGC1399, to see if they are related to the
positions of the UCDs. We found a local overabundance of globular clusters on a
scale of <1kpc around UCDs, when we compare it to the distribution of globulars
from the host galaxy. This effect is strongest for the metal-poor blue GCs. We
discuss how likely it is that these clustered globulars were originally
associated with the UCD.
We present predictions for the surface density of ultracool dwarfs (with spectral types M8-T8) for a host of deep fields that are likely to be observed with the James Webb Space Telescope. Based on simple thin and thick/thin disk (exponential) models, we show the typical distance modulus is mu~9.8 mag, which at high Galactic latitude is 5log(2 z_scl)-5. Since this is a property of the density distribution of an exponential disk, it is independent of spectral type or stellar sample. Using the published estimates of the ultracool dwarf luminosity function, we show that their number counts typically peak around J~24 mag with a total surface density of Sigma ~ 0.3 arcmin^-2, but with a strong dependence on galactic coordinate and spectral type. Owing to the exponential shape of the disk, the ultracool dwarfs are very rare at faint magnitudes (J>~27 mag), with typical densities of Sigma~0.005 arcmin^-2 (or ~20% of the total contribution within the field). Therefore in the very narrow and deep fields, we predict there are only a few ultracool dwarfs, and hence these stars are likely not a severe contaminant in searches for high-redshift galaxies. Furthermore the ultracool dwarfs are expected to be considerably brighter than the high-redshift galaxies, so samples near the faint-end of the high-redshift galaxy population will be the purest. We present the star-count formalism in a simplified way so that observers may easily predict the number of stars for their conditions (field, depth, wavelength, etc.).
We present results from high-resolution semi-global simulations of neutrino-driven convection in core-collapse supernovae. We employ an idealized setup with parametrized neutrino heating/cooling and nuclear dissociation at the shock front. We study the internal dynamics of neutrino-driven convection and its role in re-distributing energy and momentum through the gain region. We find that even if buoyant plumes are able to locally transfer heat up to the shock, convection is not able to create a net positive energy flux and overcome the downwards transport of energy from the accretion flow. Turbulent convection does, however, provide a significant effective pressure support to the accretion flow as it favors the accumulation of energy, mass and momentum in the gain region. We derive an approximate equation that is able to explain and predict the shock evolution in terms of integrals of quantities such as the turbulent pressure in the gain region or the effects of non-radial motion of the fluid. We use this relation as a way to quantify the role of turbulence in the dynamics of the accretion shock. Finally, we investigate the effects of grid resolution, which we change by a factor 20 between the lowest and highest resolution. Our results show that the shallow slopes of the turbulent kinetic energy spectra reported in previous studies are a numerical artefact. Kolmogorov scaling is progressively recovered as the resolution is increased.
We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi-Dirac cutoff power law model along with a narrow Fe K line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~$5 \times 10^{37} $erg s$^{-1}$) and lowest (~$5 \times 10^{36} $erg s$^{-1}$) observed 3-60 keV luminosities.
A popular theory of star formation is gravito-turbulent fragmentation, in which self-gravitating structures are created by turbulence-driven density fluctuations. Simple theories of isothermal fragmentation successfully reproduce the core mass function (CMF) which has a very similar shape to the initial mass function (IMF) of stars. However, numerical simulations of isothermal turbulent fragmentation thus far have not succeeded in identifying a fragment mass scale that is independent of the simulation resolution. Moreover, the fluid equations for magnetized, self-gravitating, isothermal turbulence are scale-free, and do not predict any characteristic mass. In this paper we show that, although an isothermal self-gravitating flow does produce a CMF with a mass scale imposed by the initial conditions, this scale changes as the parent cloud evolves. In addition, the cores that form undergo further fragmentation and after sufficient time forget about their initial conditions, yielding a scale-free pure power-law distribution $\mathrm{d} N/\mathrm{d} M\propto M^{-2}$ for the stellar IMF. We show that this problem can be alleviated by introducing a simple model for stellar radiation feedback. Radiative heating, powered by accretion onto forming stars, arrests the fragmentation cascade and imposes a characteristic mass scale that is nearly independent of the time-evolution or initial conditions in the star-forming cloud, and that agrees well with the peak of the observed IMF. In contrast, models that introduce a stiff equation of state for denser clouds but that do not explicitly include the effects of feedback do not yield an invariant IMF.
Our understanding of stars and their fates is based on coupling observations to theoretical models. Unlike laboratory physicists, we cannot perform experiments on stars, but rather must patiently take what nature allows us to observe. Simulation offers a means of virtual experimentation, enabling a detailed understanding of the most violent ongoing explosions in the Universe---the deaths of stars.
An eruptive nova-like event took place in 1994 in the stellar-merger candidate V4332 Sgr. Following the eruption, dust consisting of refractory silicate rich dust grains containing a significant component of AlO bonding was formed sometime between 1998 and 2003. Observations using Spitzer between 2005 and 2009 show significant changes in the 10 micron silicate stretch feature. There is a deepening of the 10 micron silicate stretch as well as the development of a feature between about 13 and 20 microns consistent with a blend of the MgO and FeO stretching features and the O-Si-O bending mode of increasingly ordered silicate dust. Near-infrared observations show the presence of AlO and water vapor in the outflow in 2003, 2004 and 2005: the AlO has significantly decreased in spectra obtained in 2014 while the water vapor remains largely unchanged. An attempt is made to correlate these observations and understand the significance of these changes using DUSTY modeling. The observations appear consistent with the kinetically-controlled, condensation of highly under-oxidized SiO/AlO/Fe/Mg dust grains in the outflow followed by the continuous evolution of the initial condensate due to thermal annealing and oxidation of the dust via reaction with ambient O, OH and H2O in the expanding, cooling shell. Periodic monitoring of this dust shell over the mid-infrared spectral range could yield useful information on the evolution of under-oxidized silicate condensates exposed to hot water vapor in more conventional circumstellar environments.
We analyse color-magnitude diagrams of eight Globular Clusters (GCs) in the outer Galactic Halo. Images were taken with the Wide Field Channel of the Advanced Camera for Survey and the Ultraviolet and Visual Channel of the Wide Field Camera 3 on board of the Hubble Space Telescope. We have determined the fraction of binary stars along the main sequence and combined results with those of a recent paper where some of us have performed a similar analysis on 59 Galactic GCs. In total, binaries have been now studied homogeneously in 67 GCs. We studied the radial and luminosity distributions of the binary systems, the distribution of their mass-ratios and investigated univariate relations with several parameters of the host GCs. We confirm the anti-correlation between the binary fraction and the luminosity of the host cluster, and find that low-luminosity clusters can host a large population in excess of ~40% in the cluster core. However, our results do not support a significant correlation with the cluster age as suggested in the literature. In most GCs, binaries are more centrally concentrated than single stars. If the fraction of binaries is normalised to the core binary fraction the radial density profiles follow a common trend. It has a maximum in the center and declines by a factor of two at a distance of about two core radii from the cluster center. After dropping to its minimum at a radial distance of $\sim$5 core radii it stays approximately constant at larger radii. We also find that the mass-ratio and the distribution of binaries as a function of the mass of the primary star is almost flat.
We have constructed timing solutions for 81 gamma-ray pulsars covering more than five years of Fermi data. The sample includes 37 radio-quiet or radio-faint pulsars which cannot be timed with other telescopes. These timing solutions and the corresponding pulse times of arrival are prerequisites for further study, e.g. phase-resolved spectroscopy or searches for mode switches. Many gamma-ray pulsars are strongly affected by timing noise, and we present a new method for characterizing the noise process and mitigating its effects on other facets of the timing model. We present an analysis of timing noise over the population using a new metric for characterizing its strength and spectral shape, namely its time-domain correlation. The dependence of the strength on spin frequency and spin-down rate is in good agreement with previous studies. We find that noise process power spectra $S(f)$ for unrecycled pulsars are steep, with strong correlations over our entire data set and spectral indices $S(f)\propto f^{-\alpha}$ of 5 to 9. One possible explanation for these results is the occurrence of unmodelled, episodic 'microglitches'. Finally, we show that our treatment of timing noise results in robust parameter estimation, and in particular we measure a precise timing position for each pulsar. We extensively validate our results with multi-wavelength astrometry, and using our updated position, we firmly identify the X-ray counterpart of PSR J1418-6058.
{\bf R}-matrix calculations combined with the adiabatic nuclei approximation are used to compute electron-impact rotiational excitation rates for three closed-shell diatomic cations, HeH$^+$, CH$^+$, ArH$^+$. Comparisons with previous studies show that an improved treatment of threshold effects leads to significant changes in the low temperature rates, furthermore the new calculations suggest that excitation of CH$^+$ is dominated by $\Delta J =1$ transitions as is expected for cations with a large dipole moment. A model for ArH$^+$ excitation in the Crab Nebula is presented which gives results consistent with the observations for electron densities in the range $2-3\times 10^3$~cm$^{-3}$.
We conducted an optical and near-infrared polarimetric observation of the highly dormant Jupiter-Family Comet, 209P/LINEAR. Because of its low activity, we were able to determine the linear polarization degrees of the coma dust particles and nucleus independently, that is $P_n$=30.3$^{+1.3}_{-0.9}$% at $\alpha$=92.2$^\circ$ and $P_n$=31.0$^{+1.0}_{-0.7}$% at $\alpha$=99.5$^\circ$ for the nucleus, and $P_c$=28.8$^{+0.4}_{-0.4}$% at $\alpha$=92.2$^\circ$ and 29.6$^{+0.3}_{-0.3}$% at $\alpha$=99.5$^\circ$ for the coma. We detected no significant variation in $P$ at the phase angle coverage of 92.2$^\circ$-99.5$^\circ$, which may imply that the obtained polarization degrees are nearly at maximum in the phase-polarization curves. By fitting with an empirical function, we obtained the maximum values of linear polarization degrees $P_\mathrm{max}$=30.8% for the nucleus and $P_\mathrm{max}$=29.6% for the dust coma. The $P_\mathrm{max}$ of the dust coma is consistent with those of dust-rich comets. The low geometric albedo of $P_v$=0.05 was derived from the slope-albedo relationship and was associated with high $P_\mathrm{max}$. We examined $P_\mathrm{max}$-albedo relations between asteroids and 209P, and found that the so-called Umov law seems to be applicable on this cometary surface.
By using subsets of the HATNet and K2 (Kepler two-wheel) Campaign 1 databases, we examine the effectiveness of filtering out systematics from photometric time series while simultaneously searching for periodic signals. We carry out tests to recover simulated sinusoidal and transit signals added to time series with both real and artificial noise. We find that the simple (and more traditional) method that performs correction for systematics first and signal search thereafter, produces higher signal recovery rates on the average, while also being substantially faster than the simultaneous method. Independently of the method of search, once the signal is found, a far less time consuming full-fledged model, incorporating both the signal and systematics, must be employed to recover the correct signal shape. As a by-product of the tests on the K2 data, we find that for longer period sinusoidal signals the detection rate decreases (after an optimum value is reached) as the number of light curves used for systematics filtering increases. The decline of the detection rate is observable in both methods of filtering, albeit the simultaneous method performs better in the regime of relative high template number. We suspect that the observed phenomenon is linked to the increased role of low amplitude intrinsic stellar variability in the space-based data. This assumption is also supported by the substantially higher stability of the detection rates for transit signals against the increase of the template number.
We investigate the water deuteration ratio and ortho-to-para nuclear spin ratio of H2 (OPR(H2)) during the formation and early evolution of a molecular cloud, following the scenario that accretion flows sweep and accumulate HI gas to form molecular clouds. We follow the physical evolution of post-shock materials using a one-dimensional shock model, with post-processing gas-ice chemistry simulations. This approach allows us to study the evolution of the OPR(H2) and water deuteration ratio without an arbitrary assumption concerning the initial molecular abundances, including the initial OPR(H2). When the conversion of hydrogen into H2 is almost complete, the OPR(H2) is already much smaller than the statistical value of three due to the spin conversion in the gas phase. As the gas accumulates, the OPR(H2) decreases in a non-equilibrium manner. We find that water ice can be deuterium-poor at the end of its main formation stage in the cloud, compared to water vapor observed in the vicinity of low-mass protostars where water ice is likely sublimated. If this is the case, the enrichment of deuterium in water should mostly occur at somewhat later evolutionary stages of star formation, i.e., cold prestellar/protostellar cores. The main mechanism to suppress water ice deuteration in the cloud is the cycle of photodissociation and reformation of water ice, which efficiently removes deuterium from water ice chemistry. The removal efficiency depends on the main formation pathway of water ice. The OPR(H2) plays a minor role in water ice deuteration at the main formation stage of water ice.
I review some steps in the conversion of molecular cloud gas into stars and planets, with an emphasis in this presentation on the early stage molecular cloud fragmentation that leads to elongated filaments/ribbons. Magnetic fields can play a crucial role in all stages and need to be invoked particularly for early stage fragmentation as well as in late core collapse where it may control disk formation. I also review some elements of hydrodynamic modeling of disk evolution.
The DiskMass survey recently provided measurements of the vertical velocity dispersions of disk stars in a sample of nearly face-on galaxies. By setting the disk scale-heights to be equal to those of edge-on galaxies with similar scale-lengths, it was found that these disks must be sub-maximal, with surprisingly low K-band mass-to-light ratios of the order of $M_\star/L_K \simeq 0.3 M_\odot/L_\odot$. This study made use of a simple relation between the disk surface density and the measured velocity dispersion and scale height of the disk, neglecting the shape of the rotation curve and the dark matter contribution to the vertical force, which can be especially important in the case of sub-maximal disks. Here, we point out that these simplifying assumptions led to an overestimation of the stellar mass-to-light ratios. Relaxing these assumptions, we compute even lower values than previously reported for the mass-to-light ratios, with a median $M_\star/L_K \simeq 0.18 M_\odot/L_\odot$, where 14 galaxies have $M_\star/L_K < 0.11$. Invoking prolate dark matter halos made only a small difference to the derived $M_\star/L_K$, although extreme prolate halos ($q>1.5$ for the axis ratios of the potential) might help. The cross-terms in the Jeans equation are also generally negligible. These deduced K-band stellar mass-to-light ratios are even more difficult to reconcile with stellar population synthesis models than the previously reported ones.
Radioactive nuclei play an important role in planetary evolution by providing an internal heat source, which affects planetary structure and helps facilitate plate tectonics. A minimum level of nuclear activity is thought to be necessary --- but not sufficient --- for planets to be habitable. Extending previous work that focused on short-lived nuclei, this paper considers the delivery of long-lived radioactive nuclei to circumstellar disks in star forming regions. Although the long-lived nuclear species are always present, their abundances can be enhanced through multiple mechanisms. Most stars form in embedded cluster environments, so that disks can be enriched directly by intercepting ejecta from supernovae within the birth clusters. In addition, molecular clouds often provide multiple episodes of star formation, so that nuclear abundances can accumulate within the cloud; subsequent generations of stars can thus receive elevated levels of radioactive nuclei through this distributed enrichment scenario. This paper calculates the distribution of additional enrichment for $^{40}$K, the most abundant of the long-lived radioactive nuclei. We find that distributed enrichment is more effective than direct enrichment. For the latter mechanism, ideal conditions lead to about 1 in 200 solar systems being directly enriched in $^{40}$K at the level inferred for the early solar nebula (thereby doubling the abundance). For distributed enrichment from adjacent clusters, about 1 in 80 solar systems are enriched at the same level. Distributed enrichment over the entire molecular cloud is more uncertain, but can be even more effective.
This paper presents a series of 95 new measurements of the longitudinal (effective) magnetic field $B_e$ of the Ap star $\gamma$ Equ (HD 201601). Observations were obtained at the coud\'e focus of the 1-m reflector at the Special Astrophysical Observatory (SAO RAS) in Russia over a time period of 4190 days (more than 11 years). We compiled a long record of $B_e$ points, adding our measurements to all published data. The time series of magnetic data consists of 395 $B_e$ points extending for 24488 days, or over 67 years. Various methods of period determination were examined for the case in which the length of the observed time series is rather short and amounts only to ~69 percent of the period. We argue that the fitting of a sine wave to the observed $B_e$ points by least squares yields the most reliable period in the case of $\gamma$ Equ. Therefore, the best period for long-term magnetic variations of $\gamma$ Equ, and hence the rotational period, is $P_{\rm rot}=35462.5 \pm 1149$ days $= 97.16 \pm 3.15$ years.
Observations of the 21 cm line radiation coming from the epoch of reionization have a great capacity to study the cosmological growth of the Universe. Also, CMB polarization produced by gravitational lensing has a large amount of information about the growth of matter fluctuations at late time. In this thesis, we investigate their sensitivities to the impact of neutrino property on the growth of density fluctuations, such as the total neutrino mass, the neutrino mass hierarchy, the effective number of neutrino species (extra radiation), and the lepton asymmetry of our Universe. We will show that by combining the precise CMB polarization observations with Square Kilometer Array (SKA) we can measure the impact of non-zero neutrino mass on the growth of density fluctuation, and determine the neutrino mass hierarchy at 2 sigma level if the total neutrino mass is smaller than 0.1 eV. Additionally, we will show that by using these combinations we can constrain the lepton asymmetry better than big-bang nucleosynthesis (BBN). Besides we discuss constraints on that in the presence of some extra radiation, and show that the 21 cm line observations can substantially improve the constraints obtained by CMB alone, and allow us to distinguish the effects of the lepton asymmetry from those of extra radiation.
Mg is the alpha element of choice for Galactic population and chemical evolution studies as it is easily detectable in all late-type stars. Such studies require precise elemental abundances, and thus departures from LTE need to be accounted for. Our goal is to provide reliable departure coefficients and equivalent widths in non-LTE, and for reference in LTE, for diagnostic lines of Mg studied in late-type stars. These can be used e.g., to correct LTE spectra and abundances. Using the model atom built and tested in the preceding paper in this series, we performed non-LTE radiative transfer calculations in a grid of 3945 stellar 1D atmospheric models. We used a sub-grid of 86 models to explore the propagation of errors in the recent atomic collision calculations to the radiative transfer results. We obtained departure coefficients for all the levels and equivalent widths (in LTE and non-LTE) for all the radiative transitions included in the "final" model atom of Osorio et al.. We present and describe our results and show some examples of applications of the data. The errors due to uncertainties in the collisional data are investigated and tabulated. The results for equivalent widths and departure coefficients are made freely available. Giants tend to have negative abundance corrections while dwarfs have positive, though small, corrections. Error analysis results show that uncertainties related to the atomic collision data are typically of order 0.01 dex or less, although for few stellar models in specific lines uncertainties can be as large as 0.03 dex. As these errors are less than or of the same order as typical corrections, we expect that we can use these results to extract Mg abundances from high quality spectra more reliably than from classical LTE analysis.
Papers on point source searches submitted to the 34th International Cosmic Ray Conference (ICRC 2015, The Hague) by the IceCube Collaboration.
Papers on atmospheric and astrophysical diffuse neutrino searches of all flavors submitted to the 34th International Cosmic Ray Conference (ICRC 2015, The Hague) by the IceCube Collaboration.
Papers on cosmic rays submitted to the 34th International Cosmic Ray Conference (ICRC 2015, The Hague) by the IceCube Collaboration.
Papers on searches for dark matter and exotic particles submitted to the 34th International Cosmic Ray Conference (ICRC 2015, The Hague) by the IceCube Collaboration.
Papers on neutrino oscillations and supernova searches submitted to the 34th International Cosmic Ray Conference (ICRC 2015, The Hague) by the IceCube Collaboration.
Papers submitted to the 34th International Cosmic Ray Conference (ICRC 2015, The Hague) by the IceCube-Gen2 Collaboration.
Several thousand solar masses of molecular, atomic and ionized gas lie in the innermost ~10 pc of our Galaxy. The most relevant structure of molecular gas is the circumnuclear ring (CNR), a dense and clumpy ring surrounding the supermassive black hole (SMBH), with a radius of ~2 pc. We propose that the CNR formed through the tidal disruption of a molecular cloud, and we investigate this scenario by means of N-body smoothed-particle hydrodynamics simulations. We ran a grid of simulations with different cloud mass (4X10^4, 1.3X10^5 solar masses), different initial orbital velocity (v_in=0.2-0.5 v_esc, where v_esc is the escape velocity from the SMBH), and different impact parameter (b=8, 26 pc). The disruption of the molecular cloud leads to the formation of very dense and clumpy gas rings, containing most of the initial cloud mass. If the initial orbital velocity of the cloud is sufficiently low (v_in<0.4 v_esc, for b=26 pc) or the impact parameter is sufficiently small (b<10 pc, for v_in>0.5 v_esc), at least two rings form around the SMBH: an inner ring (with radius ~0.4 pc) and an outer ring (with radius ~2-4 pc). The inner ring forms from low-angular momentum material that engulfs the SMBH during the first periapsis passage, while the outer ring forms later, during the subsequent periapsis passages of the disrupted cloud. The inner and outer rings are misaligned by ~24 degrees, because they form from different gas streamers, which are affected by the SMBH gravitational focusing in different ways. The outer ring matches several properties (mass, rotation velocity, temperature, clumpiness) of the CNR in our Galactic centre. We speculate that the inner ring might account for the neutral gas observed in the central cavity.
We report re-analyses of the Suzaku observations of the Galactic supernova remnant (SNR), G272.2$-$3.2, for which the previous studies were limited below 3 keV. With careful data reduction and background subtraction, we discover the K-shell lines of Ar, Ca, and Fe above 3 keV. The X-ray spectrum of G272.2$-$3.2 consists of two components, a low-temperature collisional ionization equilibrium (CIE) plasma ($kT_{\rm e} \sim 0.2$ keV) and a high-temperature non-equilibrium ionization (NEI) plasma ($kT_{\rm e} = 0.6$-$3$ keV). The CIE plasma has solar abundances over the entire area, hence it would originate from the interstellar medium. On the other hand, the abundances of the NEI plasma increase toward the inner region, suggesting the ejecta origin. The line center energy of the Fe K-shell emission ($\sim 6.4$ keV) suggests that the ejecta are recently heated by the reverse shock, a common feature in Type Ia SNRs.
We describe the technique of spectropolarimetric observations allowing for the measurements of the Stokes parameters in one of the observational modes of the SCORPIO focal reducer of the 6-m BTA telescope of the SAO RAS. The characteristics of the instrument in the spectropolarimetric mode of observations are given. We present the algorithm of observational data reduction. The capabilities of the SCORPIO spectropolarimetric mode are demonstrated on the examples of observations of various astronomical objects.
Using 2D smoothed particle hydrodynamics, we investigate the distribution of wait times between strong shocks in a turbulent, self-gravitating accretion disc. We show the resulting distributions do not depend strongly on the cooling time or resolution of the disc and that they are consistent with the predictions of earlier work (Young & Clarke 2015; Cossins et al. 2009, 2010). We use the distribution of wait times between shocks to estimate the likelihood of stochastic fragmentation by gradual contraction of shear-resistant clumps on the cooling time scale. We conclude that the stochastic fragmentation mechanism (Paardekooper 2012) cannot change the radius at which fragmentation is possible by more than ~20%, restricting direct gravitational collapse as a mechanism for giant planet formation to the outer regions of protoplanetary discs.
Rubble piles are a common feature of solar system bodies. They are composed of monolithic elements of ice or rock bound by gravity. Voids occupy a significant fraction of the volume of a rubble pile. They can exist up to pressure $P\approx \epsy\mu$, where $\epsy$ is the monolithic material's yield strain and $\mu$ its rigidity. At low $P$, contacts between neighboring elements are confined to a small fraction of their surface areas. As a result, the effective thermal conductivity of a rubble pile, $\kcon\approx k(P/(\epsy\mu))^{1/2}$, can be orders of magnitude smaller than, $k$, the thermal conductivity of its monolithic elements. In a fluid-free environment, only radiation can transfer energy across voids. It contributes an additional component, $\krad=16\ell\sigma T^3/3$, to the total effective conductivity, $\keff=\kcon +\krad$. Here $\ell$, the inverse of the opacity per unit volume, is of order the size of the elements and voids. An important distinction between $\kcon$ and $\krad$ is that the former is independent of the size of the elements whereas the latter is proportional to it. Our expression for $\keff$ provides a good fit to the depth dependence of thermal conductivity in the top $140\,\mathrm{cm}$ of the lunar regolith. It also offers a good starting point for detailed modeling of thermal inertias for asteroids and satellites. Measurement of the response of surface temperature to variable insolation is a valuable diagnostic of a regolith. There is an opportunity for careful experiments under controlled laboratory conditions to test models of thermal conductivity such as the one we outline.
We use dynamical models that include bulk rotation, velocity dispersion anisotropy and both stars and dark matter to explore the conditions that give rise to the early-type galaxy scaling relations referred to as the Fundamental Plane (FP) and Manifold (FM). The modelled scaling relations generally match the observed relations and are remarkably robust to all changes allowed within these models. The empirical relationships can fail beyond the parameter ranges where they were calibrated and we discuss the nature of those failures. Because the location of individual models relative to the FP and FM is sensitive to the adopted physical scaling of the models, unconstrained rescaling produces a much larger scatter about the scaling relations than that observed. We conclude that only certain combinations of scaling values, which define the physical radial and kinematic scale of the model, produce low scatter versions of the FP and FM. These combinations further result in reproducing a condition observed previously for galaxies, $r_c \rho_0 = $ constant, where $r_c$ is the scaling radius and $\rho_0$ is the central density. As such, we conclude that this empirical finding and global galaxy scaling relations are not independent and that finding the physical cause of one should lead to the solution to the other. Although our models are strictly for pressure supported galaxies, these results may well hold generally because the central density constraint was first identified in dwarf spheroidals but later extended to rotating giant galaxies and the FM applies to galaxies of any morphological type and luminosity class.
We report Keck/ESI and VLT/UVES observations of three super-damped Lyman-alpha quasar absorbers with H I column densities log N(HI) >= 21.7 at redshifts z=2-2.5. All three absorbers show similar metallicities (-1.3 to -1.5 dex), and dust depletion of Fe, Ni, and Mn. Two of the absorbers show supersolar [S/Zn] and [Si/Zn]. We combine our results with those for other DLAs to examine trends between N(HI), metallicity, dust depletion. A larger fraction of the super-DLAs lie close to or above the line [X/H]=20.59-log N(HI) in the metallicity vs. N(HI) plot, compared to the less gas-rich DLAs, suggesting that super-DLAs are more likely to be rich in molecules. Unfortunately, our data for Q0230-0334 and Q0743+1421 do not cover H2 absorption lines. For Q1418+0718, some H2 lines are covered, but not detected. CO is not detected in any of our absorbers. For DLAs with log N(HI) < 21.7, we confirm strong correlation between metallicity and Fe depletion, and find a correlation between metallicity and Si depletion. For super-DLAs, these correlations are weaker or absent. The absorbers toward Q0230-0334 and Q1418+0718 show potential detections of weak Ly-alpha emission, implying star formation rates of about 1.6 and 0.7 solar masses per year, respectively (ignoring dust extinction). Upper limits on the electron densities from C II*/C II or Si II*/Si II are low, but are higher than the median values in less gas-rich DLAs. Finally, systems with log N(HI) > 21.7 may have somewhat narrower velocity dispersions delta v_90 than the less gas-rich DLAs, and may arise in cooler and/or less turbulent gas.
Observations show that spiral galaxies in galaxy clusters tend to have on average less neutral hydrogen (HI) than galaxies of the same type and size in the field. There is accumulating evidence that such HI-deficient galaxies are also relatively frequent in galaxy groups. An important question is, which mechanisms are responsible for the gas deficiency in galaxy groups. To gain a better understanding of how environment affects the gas content of galaxies, we identified a sample of six HI-deficient galaxies from the HI Parkes All Sky Survey (HIPASS) using HI-optical scaling relations. One of the galaxies is located in the outskirts of the Fornax cluster, four are in loose galaxy groups and one is in a galaxy triplet. We present new high resolution HI observations with the Australia Telescope Compact Array (ATCA) of these galaxies. We discuss the possible cause of HI-deficiency in these galaxies based on HI observations and various multi-wavelength data. We find that the galaxies have truncated HI disks, lopsided gas distribution and some show asymmetries in their stellar disks. We conclude that both ram pressure stripping and tidal interactions are important gas removal mechanisms in low density environments.
The total magnification due to a point lens has been of particular interest as the theorem that gravitational lensing results in light amplification for all observers appears to contradict the conservation of photon number. This has been discussed several times, and various resolutions have been offered. In this note, we use a kinematic approach to provide a formula for the magnification factor for the primary image accurate to first order and valid for rays leaving the source at any trajectory. We thus determine the magnification over a sphere surrounding the system. A new result found is that while the magnification dips below unity far from the optical axis as noted by others, it returns to unity directly behind the source.
Filamentary structures are ubiquitous in the interstellar medium, yet their formation, internal structure, and longevity have not been studied in detail. We report the results from a comprehensive numerical study that investigates the characteristics, formation, and evolution of filaments arising from magnetohydrodynamic interactions between supersonic winds and dense clouds. Here we improve on previous simulations by utilising sharper density contrasts and higher numerical resolutions. By following multiple density tracers, we find that material in the envelopes of the clouds is removed and deposited downstream to form filamentary tails, while the cores of the clouds serve as footpoints and late-stage outer layers of these tails. Aspect ratios >12, subsonic velocity dispersions ~0.1-0.3 of the wind sound speed, and magnetic field amplifications ~100 are found to be characteristic of these filaments. We also report the effects of different magnetic field strengths and orientations. The magnetic field strength regulates vorticity production: sinuous filamentary towers arise in non-magnetic environments, while strong magnetic fields inhibit small-scale perturbations at boundary layers making tails less turbulent. Magnetic field components aligned with the direction of the flow favour the formation of pressure-confined flux ropes inside the tails, whilst transverse components tend to form current sheets. Softening the equation of state to nearly isothermal leads to suppression of dynamical instabilities and further collimation of the tail. Towards the final stages of the evolution, we find that small cloudlets and distorted filaments survive the break-up of the clouds and become entrained in the winds, reaching velocities ~0.1 of the wind speed.
Galactic black-hole X-ray binaries emit a compact, optically thick, mildy relativistic radio jet when they are in the hard and hard-intermediate states. In a series of papers, we have developed a jet model and have shown, through Monte Carlo simulations, that our model can explain many observational results. In this work, we investigate one more constraining relationship between the cutoff energy and the phase lag during the early stages of an X-ray outburst of the black-hole X-ray binary GX 339-4: the cutoff energy decreases while the phase lag increases during the brightening of the hard state. We demonstrate that our jet model naturally explains the above correlation, with a minor modification consisting of introducing an acceleration zone at the base of the jet. The observed correlation between the cutoff energy and the phase lag suggests that the lags are produced by the hard component. Here we show that this correlation arises naturally if Comptonization in the jet produces these two quantities.
We present here the first abundance analysis of 44 late B, A and F-type members of the young open cluster M6 (NGC 6405, age about 75 Myrs). Spectra, covering the 4500 to 5800 \AA{} wavelength range, were obtained using the FLAMES/GIRAFFE spectrograph attached to the ESO Very Large Telescopes (VLT). We determined the atmospheric parameters using calibrations of the Geneva photometry and by adjusting the $H_{\beta}$ profiles to synthetic ones. The abundances of up to 20 chemical elements, were derived for 19 late B, 16 A and 9 F stars by iteratively adjusting synthetic spectra to the observations. We also derived a mean cluster metallicity of $\mathrm{[Fe/H]=0.07\pm0.03}$ dex from the iron abundances of the F-type stars. We find that, for most chemical elements, the normal late B and A-type stars exhibit larger star-to-star abundance variations than the F-type stars do probably because of the faster rotation of the B and A stars. The abundances of C, O, Mg, Si and Sc appear to be anticorrelated to that of Fe, while the opposite holds for the abundances of Ca, Ti, Cr, Mn, Ni, Y, and Ba about as expected if radiative diffusion is efficient in the envelopes of these stars. In the course of this analysis, we discovered five new peculiar stars: one mild-Am, one Am, and one Fm star (HD 318091, CD-32 13109, GSC 07380-01211), one HgMn star (HD 318126), and one He-weak P-rich (HD 318101) star. We also discovered a new spectroscopic binary, most likely a SB2. We performed a detailed modelling of HD 318101,the new He-weak P-rich CP star, using the Montr\'eal stellar evolution code XEVOL which treats self-consistently all particle transport processes. Although the overall abundance pattern of this star is properly reproduced, we find that detailed abundances (in particular the high P excess) resisted modelling attempts even when a range of turbulence profiles and mass loss rates were considered.
We utilize cosmological simulations of 16 galaxy clusters at redshifts $z=0$ and $z=0.6$ to study the effect of inflowing streams on the properties of the inner Intra-Cluster Medium (ICM). We find that the mass accretion occurs predominantly along streams that originate from the cosmic web and consist of heated gas. Clusters that are unrelaxed in terms of their X-ray morphology are characterized by higher mass inflow rates and deeper penetration of the streams, typically into the inner third of the virial radius. The penetrating streams generate elevated random motions, bulk flows, cold fronts and metal mixing, thus producing Non-Cool-Core clusters. The degree of penetration of the streams may change over time such that clusters can switch from being unrelaxed to relaxed over a time-scale of several Gyrs. The stream properties thus help us understand the distinction between cool-core and non-cool-core clusters.
We use photometry in the F220W, F250W, F330W, F435W filters from the High Resolution Channel of the Advanced Camera for Surveys and photometry in the F555W, F675W, and F814W filters from the Wide Field and Planetary Camera 2 aboard the Hubble Space Telescope to derive individual stellar reddenings and extinctions for stars in the HD 97950 cluster in the giant HII region NGC 3603. The mean line-of-sight reddening for about a hundred main-sequence member stars inside the cluster is $E(F435W-F555W)=1.33\pm0.12$ mag. After correcting for foreground reddening, the total to selective extinction ratio is $R_{F555W}=3.75\pm0.87$ in the cluster. Within the standard deviation associated with $E(\rm \lambda-F555W)/E(F435W-F555W)$ in each filter, the cluster extinction curve at ultraviolet wavelengths tends to be greyer than the average Galactic extinction laws from Cardelli et al. (1989) and Fitzpatrick et al. (1999). It is closer to the extinction law derived by Calzetti et al. (2000) for starburst galaxies, where the 0.2175 $\rm \mu m$ bump is absent. This indicates an anomalous extinction in the HD 97950 cluster, which may due to the clumpy dust distribution within the cluster, and the size of dust grains being larger than the average Galactic ISM.
I review a new rapidly growing area of high-energy plasma astrophysics --- radiative magnetic reconnection, i.e., a reconnection regime where radiation reaction influences reconnection dynamics, energetics, and nonthermal particle acceleration. This influence be may be manifested via a number of astrophysically important radiative effects, such as radiation-reaction limits on particle acceleration, radiative cooling, radiative resistivity, braking of reconnection outflows by radiation drag, radiation pressure, viscosity, and even pair creation at highest energy densities. Self-consistent inclusion of these effects in magnetic reconnection theory and modeling calls for serious modifications to our overall theoretical approach to the problem. In addition, prompt reconnection-powered radiation often represents our only observational diagnostic tool for studying remote astrophysical systems; this underscores the importance of developing predictive modeling capabilities to connect the underlying physical conditions in a reconnecting system to observable radiative signatures. This Chapter gives an overview of recent theoretical progress in developing basic physical understanding of radiative reconnection, with a special emphasis on astrophysically important radiation mechanisms like synchrotron, curvature, and inverse-Compton. It also offers a broad review of key high-energy astrophysical applications of radiative reconnection, such as: pulsar wind nebulae and magnetospheres, accreting black-hole coronae and jets in XRBs and AGN, magnetospheres of magnetars, and Gamma-Ray Bursts. Finally, this Chapter discusses the most critical open questions and outlines the directions for future research in this exciting new frontier of plasma astrophysics.
(Abridged) The chemical behaviour of an ample sample of PNe in NGC6822 is analyzed. Spectrophotometric data of 11 PNe and two H II regions were obtained with the OSIRIS spectrograph attached to the Gran Telescopio Canarias. Data for other 13 PNe and three H II regions were retrieved from the literature. Physical conditions and chemical abundances of O, N, Ne, Ar and S were derived for 19 PNe and 4 H II regions. Abundances in the PNe sample are widely distributed showing 12+log(O/H) from 7.4 to 8.2 and 12+log(Ar/H) from 4.97 to 5.80. Two groups of PNe can be differentiated: one old, with low metallicity (12+log(O/H)<8.0 and 12+log(Ar/H)<5.7) and another younger with metallicities similar to the values of H II regions. The old objects are distributed in a larger volume than the young ones. An important fraction of PNe (>30%) was found to be highly N-rich (Type I PNe). Such PNe occur at any metallicity. In addition, about 60% of the sample presents high ionization (He++/He >= 0.1), possessing a central star with effective temperature larger than 10^6 K. Possible bias in the sample are discussed. From comparison with stellar evolution models by A. Karakas's group of the observed N/O abundance ratios, our PNe should have had initial masses lower than 4 M_sun, although if the comparison is made with Ne vs. O abundances, the initial masses should have been lower than 2 M_sun. It appears that these models of stars of 2-3 M_sun are producing too much 22Ne in the stellar surface at the end of the AGB. On the other hand, the comparison with another set of stellar evolution models by P. Ventura's group with a different treatment of convection and on the assumptions concerning the overshoot of the convective core during the core H-burning phase, provided a reasonable agreement between N/O and Ne/H observed and predicted ratios if initial masses of more massive stars are of about 4 M_sun.
Low mass dwarf spheroidal galaxies are key objects for our understanding of
the chemical evolution of the pristine Universe and the Local Group of
galaxies. Abundance ratios in stars of these objects can be used to better
understand their star formation and chemical evolution. We report on the
analysis of a sample of 11 stars belonging to 5 different ultra faint dwarf
spheroidal galaxies (UfDSph) based on X-Shooter spectra obtained at the VLT.
Medium resolution spectra have been used to determine the detailed chemical
composition of their atmosphere. We performed a standard 1D LTE analysis to
compute the abundances.
Considering all the stars as representative of the same population of low
mass galaxies, we found that the [alpha/Fe] ratios vs [Fe/H] decreases as the
metallicity of the star increases in a way similar to what is found for the
population of stars belonging to dwarf spheroidal galaxies. The main difference
is that the solar [alpha/Fe] is reached at a much lower metallicity for the
UfDSph than the dwarf spheroidal galaxies.
We report for the first time the abundance of strontium in CVnI. The star we
analyzed in this galaxy has a very high [Sr/Fe] and a very low upper limit of
barium which makes it a star with an exceptionally high [Sr/Ba] ratio.
Our results seem to indicate that the galaxies which have produced the bulk
of their stars before the reionization (fossil galaxies) have lower [X/Fe]
ratios at a given metallicity than the galaxies that have experienced a
discontinuity in their star formation rate (quenching).
Our aim is to find and classify OB stars in Sextans A, to later determine
accurate stellar parameters of these blue massive stars in this low metallicity
region $(Z \sim 0.1 \rm Z_{\odot})$.
Using UBV photometry, the reddening-free index Q and GALEX imaging, we built
a list of blue massive star candidates in Sextans A. We obtained low resolution
(R $\sim$ 1000) GTC-OSIRIS spectra for a fraction of them and carried out
spectral classification. For the confirmed O-stars we derive preliminary
stellar parameters.
The target selection criteria and observations were successful and have
produced the first spectroscopic atlas of OB-type stars in Sextans A. From the
whole sample of 18 observed stars, 12 were classified as early OB-types,
including 5 O-stars. The radial velocities of all target stars are in agreement
with their Sextans A membership, although three of them show significant
deviations. We determined the stellar parameters of the O-type stars using the
stellar atmosphere code FASTWIND, and revisited the sub-SMC temperature scale.
Two of the O-stars are consistent with relatively strong winds and enhanced
helium abundances, although results are not conclusive. We discuss the position
of the OB stars in the HRD. Initial stellar masses run from slightly below 20
up to 40 solar masses.
The target selection method worked well for Sextans A, confirming the
procedure developed in Garcia \& Herrero (2013). The stellar temperatures are
consistent with findings in other galaxies. Some of the targets deserve
follow-up spectroscopy because of indications of a runaway nature, an enhanced
helium abundance or a relatively strong wind. We observe a correlation between
HI and OB associations similar to the irregular galaxy IC1613, confirming the
previous result that the most recent star formation of Sextans A is currently
on-going near the rim of the H\,{\sc I} cavity.
Ice mantles that formed on top of dust grains are photoprocessed by the secondary ultraviolet (UV) field in cold and dense molecular clouds. UV photons induce photochemistry and desorption of ice molecules. Experimental simulations dedicated to ice analogs under astrophysically relevant conditions are needed to understand these processes. We present UV-irradiation experiments of a pure CO2 ice analog. Calibration of the QMS allowed us to quantify the photodesorption of molecules to the gas phase. This information was added to the data provided by the FTIR on the solid phase to obtain a complete quantitative study of the UV photoprocessing of an ice analog. Experimental simulations were performed in an ultra-high vacuum chamber. Ice samples were deposited onto an infrared transparent window at 8K and were subsequently irradiated with a microwave-discharged hydrogen flow lamp. After irradiation, ice samples were warmed up until complete sublimation was attained. Photolysis of CO2 molecules initiates a network of photon-induced chemical reactions leading to the formation of CO, CO3 ,O2 , and O3 . During irradiation, photon-induced desorption of CO and, to a lesser extent, O2 and CO2 took place through a process called indirect desorption induced by electronic transitions (DIET), with maximum photodesorption yields (Ypd) of 1.2 x 10-2 molecules/incident photon , 9.3 x 10-4 molecules/incident photon , and 1.1 x 10-4 molecules/incident photon , respectively. Calibration of mass spectrometers allows a direct quantification of photodesorption yields instead of the indirect values that were obtained from infrared spectra in most previous works. Supplementary information provided by infrared spectroscopy leads to a complete quantification, and therefore a better understanding, of the processes taking place in UV-irradiated ice mantles.
In protoplanetary disks micron-size dust grains coagulate to form larger structures with complex shapes and compositions. The coagulation process changes the absorption and scattering properties of particles in the disk in significant ways. To properly interpret observations of protoplanetary disks and to place these observations in the context of the first steps of planet formation, it is crucial to understand the optical properties of these complex structures. We derive the optical properties of dust aggregates using detailed computations of aggregate structures and compare these computa- tionally demanding results with approximate methods that are cheaper to compute in practice. In this way we wish to understand the merits and problems of approximate methods and define the context in which they can or cannot be used to analyze observations of objects where significant grain growth is taking place. For the detailed computations we used the discrete dipole approximation (DDA), a method able to compute the interaction of light with a complexly shaped, inhomogeneous particle. We compared the results to those obtained using spherical and irregular, homogeneous and inhomogeneous particles. While no approximate method properly reproduces all characteristics of large dust aggregates, the thermal properties of dust can be analyzed using irregularly shaped, porous, inhomogeneous grains. The asymmetry of the scattering phase function is a good indicator of aggregate size, while the degree of polarization is probably determined by the size of the constituent particles. Optical properties derived from aggregates significantly differ from the most frequently used standard ("astronomical silicate" in spherical grains). We outline a computationally fast and relatively accurate method that can be used for a multiwavelength analysis of aggregate dust in protoplanetary disks.
We present the first application of a new foreground removal pipeline to the current leading HI intensity mapping dataset, obtained by the Green Bank Telescope (GBT). We study the 15hr and 1hr field data of the GBT observations previously presented in Masui et al. (2013) and Switzer et al. (2013) covering about 41 square degrees at 0.6 < z < 1.0 which overlaps with the WiggleZ galaxy survey employed for the cross-correlation with the maps. In the presented pipeline, we subtract the Galactic foreground continuum and the point source contaminations using an independent component analysis technique (fastica) and develop a description for a Fourier-based optimal weighting estimator to compute the temperature power spectrum of the intensity maps and cross-correlation with the galaxy survey data. We show that fastica is a reliable tool to subtract diffuse and point-source emission by using the non-Gaussian nature of their probability functions. The power spectra of the intensity maps and the cross-correlation with WiggleZ is typically an order of magnitude higher than the previous findings by the GBT team. fastica is a very conservative subtraction technique and is not able to remove anisotropic noise contaminations caused by instrumental systematics unlike the singular value decomposition method which does not discriminate components according to their statistical properties. We confirm that foreground subtraction with fastica is robust against 21cm signal loss as seen by the converged amplitude of the cross-correlation of the intensity maps with the WiggleZ data.
We present an analysis of the photometric data of comet C/2011 L4 (PANSTARRS) observed at heliocentric distance of 4.4 - 4.2 AU. The comet C/2011 L4 (PANSTARRS) shows one significant activity, despite of its quite large heliocentric distance. The color indexes, dust mass-loss rates and radius of the comet are measured.
We report the discovery of four short-period eclipsing systems in the Kepler light curves, consisting of an A-star primary and a low-mass white-dwarf (WD) secondary (dA+WD) - KIC 4169521, KOI-3818, KIC 2851474 and KIC 9285587. The systems show BEaming, Ellipsoidal and Reflection (BEER) phase modulations together with primary and secondary eclipses. These add to the 6 Kepler, and 18 WASP, previously known short-period eclipsing dA+WD binaries. The light curves together with follow-up spectroscopic observations allow us to derive the masses, radii and effective temperatures of the two components of the four systems. The orbital periods, of 1.17-3.82 d, and WD masses, of 0.19-0.22 Msun, are similar to those of the previously known systems. The WD radii of KOI-3818, KIC 2851474, and KIC 9285587 are 0.026, 0.035 and 0.026 Rsun, respectively, the smallest WD radii derived so far for short-period eclipsing dA+WD binaries. These three binaries extend the previously known population to older systems with cooler and smaller WD secondaries. KOI-3818 displays evidence for a fast-rotating primary and a minute but significant eccentricity of ~0.0015. These features are probably the outcome of the mass-transfer process.
The arrangement of the fiber positioning units in LAMOST focal plane may lead to the collisions during the fiber allocation. To avoid these collisions, the soft protection system has to abandon some targets located in the overlapped field of the adjacent fiber units. In this paper, we firstly analyzed the probability of the collisions between fibers and inferred their possible reasons. It is useful to solve the problem of the fiber-positioning units collisions so as to improve LAMOST efficiency. Based on it, a collision handling system is designed by using the master-slave control structure between the micro control unit (MCU) and the microcomputer. The simulated experiments validate that the system can provide real-time inspection and swap the information between the fiber unit controllers and the main controller.
We present 3D hydrodynamic simulations of the adiabatic interaction of a
shock with a dense, spherical cloud. We compare how the nature of the
interaction changes with the Mach number of the shock, M, and the density
contrast of the cloud, chi. We also examine the differences with 2D
axisymmetric calculations, perform detailed resolution tests, and compare
``inviscid'' results to those obtained with the inclusion of a k-epsilon
subgrid turbulence model.
We find that resolutions of 32-64 cells per cloud radius are the minimum
necessary to capture the dominant dynamical processes in 3D simulations. In
contrast to our earlier 2D work, we find that 3D inviscid and k-epsilon
simulations typically show very good agreement. As such, there does not appear
to be any compelling reason for using the k-epsilon subgrid model in 3D
calculations, though it remains very useful for 2D calculations. Clouds
accelerate and mix up to 5 times faster when they are poorly resolved. This has
implications for numerical simulations of multi-phase flows where a fast, low
density medium interacts with slower, higher density clouds (e.g., galactic
winds).
The interaction proceeds very similarly in 2D and 3D - although non-azimuthal
modes lead to different behaviour, there is very little effect on key global
quantities such as the lifetime of the cloud and its acceleration. We do not
find significant differences in the hollowing or ``voiding'' of the cloud
between 2D and 3D simulations with M=10 and chi=10, in contradiction to
expectations. This may be due to the softer edges used for our clouds.
The biggest differences between our 2D and 3D calculations are found when
M=1.5 and chi=10 - the cloud is destroyed more rapidly in 2D simulations,
perhaps because secondary vortices form earlier and are more prevelant in the
higher resolution 2D simulations.
Hot subdwarf B stars (sdBs) are the stripped cores of red giants located at the bluest extension of the horizontal branch. Several different kinds of pulsators are found among those stars. The mechanism that drives those pulsations is well known and the theoretically predicted instability regions for both the short-period p-mode and the long-period g-mode pulsators match the observed distributions fairly well. However, it remains unclear why only a fraction of the sdB stars pulsate, while stars with otherwise very similar parameters do not show pulsations. From an observers perspective I review possible candidates for the missing parameter that makes sdB stars pulsate or not.
Quasi--stellar objects (quasars) located behind nearby galaxies provide an excellent absolute reference system for astrometric studies, but they are difficult to identify because of fore- and background contamination. Deep wide--field, high angular resolution surveys spanning the entire area of nearby galaxies are needed to obtain a complete census of such quasars. We embarked on a program to expand the quasar reference system behind the Large and the Small Magellanic Clouds, the Magellanic Bridge, and the Magellanic Stream, connecting the Clouds with the Milky Way. Hundreds of quasar candidates were selected based on their near--infrared colors and variability properties from the ongoing public ESO VISTA Magellanic Clouds survey. A subset of 49 objects was followed up with optical spectroscopy. We confirmed the quasar nature of 37 objects (34 new identifications), four are low redshift objects, three are probably stars, and the remaining three lack prominent spectral features for a secure classification; bona fide quasars, judging from their broad absorption lines are located, as follows: 10 behind the LMC, 13 behind the SMC, and 14 behind the Bridge. The quasars span a redshift range from z~0.5 to z~4.1. Upon completion the VMC survey is expected to yield a total of ~1500 quasars with Y<19.32 mag, J<19.09 mag, and Ks<18.04 mag.
An evolution of luminosity of galaxies in emission lines or wavelength ranges
in which they are sensitive to the star formation process is caused by burning
out of the most massive O-class stars during a few million years after a
starburst. We study the impact of this effect on the luminosity function (LF)
of a sample of star-forming galaxies.
We introduce several types of LFs: an initial LF after a starburst, current,
time-averaged and sample ones. We find the relations between them in general
and specify them in the case of the luminosity evolution law proposed for the
luminous compact galaxies. We obtain the sample LF for the cases the initial
one is described by the pure Schechter function or the log-normal distribution
and analyze the properties of these LFs. As a result we get two new types of
LFs to fit the LF of a sample of star-forming galaxies.
We study the effect of long gradient modes on large scale observables. When defined correctly, genuine observables should not only be gauge invariant but also devoid of any gauge artifacts. One such gauge artifact is a pure gradient mode. Using the relativistic formulation of large scale observables, we confirm that a long gradient mode which is still outside observer's horizon leaves no imprint on the large scale observables at first order. These include the cosmic rulers and the number counts. This confirms the existing method for relativistically defined observables. The general relativistic bias relation for the halos and galaxies is also invariant under the presence of a long gradient mode perturbation. The observed power spectrum is not affected by this long mode.
By solving analytically the various types of Lane-Emden equations with rotation, we have discovered two new coupled fundamental properties of rotating, self-gravitating, gaseous disks in equilibrium: Isothermal disks must, on average, exhibit strict power-law density profiles in radius $x$ on their equatorial planes of the form $A x^{k-1}$, where $A$ and $k-1$ are the integration constants; and ``flat'' rotation curves precisely such as those observed in spiral galaxy disks. Polytropic disks must, on average, exhibit strict density profiles of the form $\left[\ln(A x^k)\right]^n$, where $n$ is the polytropic index; and ``flat'' rotation curves described by square roots of upper incomplete gamma functions. By ``on average,'' we mean that, irrespective of the chosen boundary conditions, the actual profiles must oscillate around and remain close to the strict mean profiles of the analytic singular equilibrium solutions. We call such singular solutions the ``intrinsic'' solutions of the differential equations because they are demanded by the second-order equations themselves with no regard to the Cauchy problem. The results are directly applicable to gaseous galaxy disks that have long been known to be isothermal and to protoplanetary disks during the extended isothermal and adiabatic phases of their evolution. In galactic gas dynamics, they have the potential to resolve the dark matter--modified gravity controversy in a sweeping manner, as they render both of these hypotheses unnecessary. In protoplanetary disk research, they provide observers with powerful new probing tool, as they predict a clear and simple connection between the radial density profiles and the rotation curves of self-gravitating disks in their very early (pre-Class 0 and perhaps the youngest Class Young Stellar Objects) phases of evolution.
We report the discovery of one RR Lyrae star in the ultra--faint satellite galaxy Hydra II based on time series photometry in the g, r and i bands obtained with the Dark Energy Camera at Cerro Tololo Interamerican Observatory, Chile. The RR Lyrae star has a mean magnitude of $i = 21.30\pm 0.04$ which translates to a heliocentric distance of $151\pm 8$ kpc for Hydra II; this value is $\sim 13\%$ larger than the estimate from the discovery paper based on the average magnitude of several blue horizontal branch star candidates. The new distance implies a slightly larger half-light radius of $76^{+12}_{-10}$ pc and a brighter absolute magnitude of $M_V = -5.1 \pm 0.3$, which keeps this object within the realm of the dwarf galaxies. The pulsational properties of the RR Lyrae star ($P=0.645$ d, $\Delta g = 0.68$ mag) suggest Hydra II may be a member of the intermediate Oosterhoff or Oosterhoff II group. A comparison with other RR Lyrae stars in ultra--faint systems indicates similar pulsational properties among them, which are different to those found among halo field stars and those in the largest of the Milky Way satellites. We also report the discovery of 31 additional short period variables in the field of view (RR Lyrae, SX Phe, eclipsing binaries, and a likely anomalous cepheid). However, given their magnitudes and large angular separation from Hydra II, they must be field stars not related to Hydra II.
Despite the lack of a shear-rich tachocline region low-mass fully convective stars are capable of generating strong magnetic fields, indicating that a dynamo mechanism fundamentally different from the solar dynamo is at work in these objects. We present a self-consistent three dimensional model of magnetic field generation in low-mass fully convective stars. The model utilizes the anelastic magnetohydrodynamic equations to simulate compressible convection in a rotating sphere. A distributed dynamo working in the model spontaneously produces a dipole-dominated surface magnetic field of the observed strength. The interaction of this field with the turbulent convection in outer layers shreds it, producing small-scale fields that carry most of the magnetic flux. The Zeeman-Doppler-Imaging technique applied to synthetic spectropolarimetric data based on our model recovers most of the large-scale field. Our model simultaneously reproduces the morphology and magnitude of the large-scale field as well as the magnitude of the small-scale field observed on low-mass fully convective stars.
We evaluate the local variance of the Hubble Constant $H_0$ with low-z Type Ia Supernovae (SNe). Our analyses are performed using a hemispherical comparison procedure to test whether the bulk flow motion can reconcile the measurement of the Hubble Constant $H_0$ from standard candles ($H_0 = 73.8 \pm 2.4 \; \mathrm{km \; s}^{-1}\; \mathrm{Mpc}^{-1}$) with that of the Planck's Cosmic Microwave Background data ($67.8 \pm 0.9 \; \mathrm{km \; s}^{-1} \mathrm{Mpc}^{-1}$). We obtain that $H_0$ ranges from $68.9 \pm 0.5 \; \mathrm{km \; s}^{-1} \mathrm{Mpc}^{-1}$ to $71.2 \pm 0.7 \; \mathrm{km \; s}^{-1} \mathrm{Mpc}^{-1}$ through the celestial sphere, with maximal dipolar anisotropy towards the $(l,b) = (315^{\circ},27^{\circ})$ direction. Interestingly, this result is in good agreement with both $H_0$ estimations, as well as the bulk flow direction reported in the literature. In addition, we assess the statistical significance of this variance with different prescriptions of Monte Carlo simulations, finding a good concordance for its amplitude given the limitation of the SNe data set in terms of celestial coverage and their distance uncertainties. Additionally, we test the null hypothesis of this anisotropic direction being a random effect, albeit we can reject such hypothesis with good statistical significance. We then conclude that the conflict between these different Hubble Constant determinations can be unified by correctly taking into account the bulk flow motion of the local Universe.
The Erigone family is a C-type group in the inner main belt. Its age has been
estimated by several researchers to be less then 300 My, so it is a relatively
young cluster. Yarko-YORP Monte Carlo methods to study the chronology of the
Erigone family confirm results obtained by other groups. The Erigone family,
however, is also characterized by its interaction with the $z_2$ secular
resonance. While less than 15% of its members are currently in librating states
of this resonance, the number of objects, members of the dynamical group, in
resonant states is high enough to allow to use the study of dynamics inside the
$z_2$ resonance to set constraints on the family age.
Like the ${\nu}_{6}$ and $z_1$ secular resonances, the $z_2$ resonance is
characterized by one stable equilibrium point at $\sigma = 180^{\circ}$ in the
$z_2$ resonance plane $(\sigma, \frac{d\sigma}{dt})$, where $\sigma$ is the
resonant angle of the $z_2$ resonance. Diffusion in this plane occurs on
timescales of $\simeq 12$ My, which sets a lower limit on the Erigone family
age. Finally, the minimum time needed to reach a steady-state population of
$z_2$ librators is about 90 My, which allows to impose another, independent
constraint on the group age.
We present a publicly-available spectral model for thermal X-ray emission from a baryonic jet in an X-ray binary system, inspired by the microquasar SS 433. The jet is assumed to be strongly collimated (half-opening angle $\Theta\sim 1\deg$) and mildly relativistic (bulk velocity $\beta=V_{b}/c\sim 0.03-0.3$). Its X-ray spectrum is found by integrating over thin slices of constant temperature, radiating in optically thin coronal regime. The temperature profile along the jet and corresponding differential emission measure distribution are calculated with full account for gas cooling due to expansion and radiative losses. Since the model predicts both the spectral shape and luminosity of the jet's emission, its normalisation is not a free parameter if the source distance is known. We also explore the possibility of using simple X-ray observables (such as flux ratios in different energy bands) to constrain physical parameters of the jet (e.g. gas temperature and density at its base) without broad-band fitting of high-resolution spectra. We demonstrate this approach in application to Chandra HETGS spectra of SS 433 in its 'edge-on' precession phase, when the contribution from non-jet spectral components is expected to be low. Our model provides a reasonable fit to the 1-3 keV data, while some residuals remain at higher energies, which may be partially attributed to a putative reflection component. Besides SS 433, the model might be used for describing jet components in spectra of other Galactic XRBs (e.g. 4U 1630-47), ULXs (e.g. Holmberg II X-1), and candidate SS 433 analogues like S26 in NGC7793 and the radio transient in M82.
The extreme extragalactic sources known as Ultraluminous X-ray Sources (ULX) represent a unique testing environment for compact objects population studies and the accretion process. Their nature has long been disputed. Their luminosity, well above the Eddington luminosity for a stellar-mass black hole, can imply the presence of an intermediate-mass black hole or a stellar black hole accreting above the Eddington limit. Both these interpretations are important to understand better the accretion process and the evolution of massive black holes. The last few years have seen a dramatic improvement of our knowledge of these sources. In particular, the super-Eddington interpretation for the bulk of the ULX population has gained a strong consensus. Nonetheless, exceptions to this general trend do exist, and in particular one ULX was shown to be a neutron star, and another was shown to be a very likely IMBH candidate. In this paper, I will review the progress done in the last few years.
We propose a model for short duration gamma-ray bursts (sGRBs) based on the formation of a quark star after the merger of two neutron stars. We assume that the sGRB central engine is a proto-magnetar, which has been previously invoked to explain the plateau-like X-ray emission observed following both long and short GRBs. Here, we show that: i) a few milliseconds after the merger it is possible to form a stable and massive star made in part of quarks; ii) during the early cooling phase of the incompletely formed quark star, the flux of baryons ablated from the surface by neutrinos is large and it does not allow the outflow to achieve a bulk Lorentz factor high enough to produce a GRB; iii) after the quark burning front reaches the stellar surface, baryon ablation ceases and the jet becomes too baryon poor to produce a GRB; iv) however, between these two phases a GRB can be produced over the finite timescale required for the baryon pollution to cease; a characteristic timescale of the order of $\sim 0.1 $ s naturally results from the time the conversion front needs to cover the distance between the rotational pole and the latitude of the last closed magnetic field line; v) we predict a correlation between the luminosity of the sGRB and its duration, consistent with the data; vi) our model also predicts a delay of the order of ten seconds between the time of the merger event and the sGRB, allowing for the possibility of precursor emission and implying that the jet will encounter the dense cocoon formed immediately after the merger.
Narrow line Seyfert 1 galaxies (NLS1s) usually do not host relativistic jet and the $\gamma$-ray NLS1s are expected to be rare. All $\gamma$-ray NLS1s reported to date have flat radio spectra and the jets are found to be closely aligned. No $\gamma$-ray mis-aligned NLS1 has been predicted before. In this work we analyze the first seven-year $Fermi$/Large Area Telescope (LAT) data of a steep radio spectrum NLS1 B3 1441+476 and report the {\it first} detection of $\gamma$-rays in such a kind of objects. No rapid variability is observed from radio to $\gamma$ rays and additionally low core dominance ($\lesssim$ 0.7) and Compton dominance ($\lesssim$ 1) are found. B3 1441+476 has a compact radio morphology and a radio spectrum turnover at $\sim$ 100MHz. A radiation model successfully reproducing some steep-spectrum radio quasars can reasonably fit the spectral energy distribution of B3 1441+476. All these facts strongly suggest that B3 1441+476 hosts a mis-aligned and plausibly underdeveloped relativistic jet, which provides a valuable target to reveal the formation and evolution of relativistic jets in NLS1s.
The \textit{Spitzer} SAGE survey has allowed the identification and analysis of significant samples of Young Stellar Object (YSO) candidates in the Large Magellanic Cloud (LMC). However the angular resolution of \textit{Spitzer} is relatively poor meaning that at the distance of the LMC, it is likely that many of the \textit{Spitzer} YSO candidates in fact contain multiple components. We present high resolution \textit{K}-band integral field spectroscopic observations of the three most prominent massive YSO candidates in the N113 H\,{\sc ii} region using VLT/SINFONI. We have identified six \textit{K}-band continuum sources within the three \textit{Spitzer} sources and we have mapped the morphology and velocity fields of extended line emission around these sources. Br$\gamma$, He\,{\sc i} and H$_2$ emission is found at the position of all six \textit{K}-band sources; we discuss whether the emission is associated with the continuum sources or whether it is ambient emission. H$_2$ emission appears to be mostly ambient emission and no evidence of CO emission arising in the discs of YSOs has been found. We have mapped the centroid velocities of extended Br$\gamma$ emission and He {\sc i} emission and found evidence of two expanding compact H\,{\sc ii} regions. One source shows compact and strong H$_2$ emission suggestive of a molecular outflow. The diversity of spectroscopic properties observed is interpreted in the context of a range of evolutionary stages associated with massive star formation.
We re-analyse the publicly available radial velocity (RV) measurements for Alpha Cen B, a star hosting an Earth-mass planet candidate, Alpha Cen Bb, with 3.24 day orbital period. We demonstrate that the 3.24 d signal observed in the Alpha Cen B data almost certainly arises from the window function (time sampling) of the original data. We show that when stellar activity signals are removed from the RV variations, other significant peaks in the power spectrum of the window function are coincidentally suppressed, leaving behind a spurious yet apparently-significant 'ghost' of a signal that was present in the window function's power spectrum to begin with. Even when fitting synthetic data with time sampling identical to the original data, but devoid of any genuine periodicities close to that of the planet candidate, the original model used to infer the presence of Alpha Cen Bb leads to identical conclusions: viz., the 3$\sigma$ detection of a half-a-metre-per-second signal with 3.236 day period. Our analysis underscores the difficulty of detecting weak planetary signals in RV data, and the importance of understanding in detail how every component of an RV data set, including its time sampling, influences final statistical inference.
Stellar seismology appears more and more as a powerful tool for a better determination of the fundamental properties of solar-type stars. However the particular case of Sun is still challenging. The helioseismic sound speed determination continues to disagree with the Standard Solar Model (SSM) prediction for about a decade, questioning the reliability of this model. One of the sources of uncertainty could be in the treatment of the transport of radiation from the solar core to the surface. In this letter, we use the new OPAS opacity tables, recently available for solar modelling, to address this issue. We discuss first the peculiarities of these tables, then we quantify their impact on the solar sound speed and density profiles using the reduced OPAS tables taken on the grids of the OPAL ones. We use the two evolution codes MESA and CLES that led to similar conclusions in the solar radiative zone. In comparison to commonly used OPAL opacity tables, the new solar models computed, for the most recent photospheric composition, with OPAS tables present improvements in the location of the base of the convective zone and in the description of the solar radiative zone in comparison to the helioseismic observations, even if the differences in the Rosseland mean opacity do not exceed 6 %. We finally carry out a comparison to a solar model computed with the OP opacity tables.
We explore the radial velocity performance benefits of coupling starlight to a fast-scanning interferometer and a fast-readout spectrometer with zero readout noise. By rapidly scanning an interferometer we can decouple wavelength calibration errors from precise radial velocity measurements, exploiting the advantages of lock-in amplification. In a Bayesian framework, we investigate the correlation between wavelength calibration errors and resulting radial velocity errors. We construct an end-to-end simulation of this approach to address the feasibility of achieving 10 cm/s radial velocity precision on a typical Sun-like star using existing, 5-meter-class telescopes. We find that such a precision can be reached in a single night, opening up possibilities for ground-based detections of Earth-Sun analog systems.
We investigate the Hierarchical Gravitational Fragmentation scenario through numerical simulations of the prestellar stages of the collapse of a marginally gravitationally unstable isothermal sphere immersed in a strongly gravitationally unstable, uniform background medium. The core developes a Bonnor-Ebert (BE)-like density profile, while at the time of singularity (the protostar) formation the envelope approaches a singular-isothermal-sphere (SIS)-like $r^-2$ density profile. However, these structures are never hydrostatic. In this case, the central flat region is characterized by an infall speed, while the envelope is characterized by a uniform speed. This implies that the hydrostatic SIS initial condition leading to Shu's classical inside-out solution is not expected to occur, and therefore neither should the inside-out solution. Instead, the solution collapses from the outside-in, naturally explaining the observation of extended infall velocities. The core, defined by the radius at which it merges with the background, has a time-variable mass, and evolves along the locus of the ensemble of observed prestellar cores in a plot of $M/M_{BE}$ vs. $M$, where $M$ is the core's mass and $M_{BE}$ is the critical Bonnor-Ebert mass, spanning the range from the "stable" to the "unstable" regimes, even though it is collapsing at all times. We conclude that the presence of an unstable background allows a core to evolve dynamically from the time when it first appears, even when it resembles a pressure-confined, stable BE-sphere. The core can be thought of as a ram-pressure confined BE-sphere, with an increasing mass due to the accretion from the unstable background.
W Serpentids and Double Periodic Variables (DPVs) are candidates for close interacting binaries in a non-conservative evolutionary stage; while W Serpentids are defined by high-excitation ultraviolet emission lines present during most orbital phases, and by usually showing variable orbital periods, DPVs are characterized by a long photometric cycle lasting roughly 33 times the (practically constant) orbital period. We report the discovery of 7 new Galactic DPVs, increasing the number of known DPVs in our Galaxy by 50%. We find that DPVs are tangential-impact systems, i.e. their primaries have radii barely larger than the critical Lubow-Shu radius. These systems are expected to show transient discs, but we find that they host stable discs with radii smaller than the tidal radius. Among tangential-impact systems including DPVs and semi-detached Algols, only DPVs have primaries with masses between 7 and 10 $M_{\odot}$. We find that DPVs are in a Case-B mass transfer stage with donor masses between 1 and 2 M$_{\odot}$ and with primaries resembling Be stars. W Serpentids are impact and non-impact systems, their discs extend until the last non-intersecting orbit and show a larger range of stellar mass and mass ratio than DPVs. Infrared photometry reveals significant color excesses in many DPVs and W Serpentids, usually larger for the latter ones, suggesting variable amounts of circumstellar matter.
Ultra-high energy cosmic rays (UHECRs) are the highest energy messengers of the present universe, with energies up to $10^{20}$ eV. Studies of astrophysical particles (nuclei, electrons, neutrinos and photons) at their highest observed energies have implications for fundamental physics as well as astrophysics. The primary particles interact in the atmosphere and generate extensive air showers. Analysis of those showers enables one not only to estimate the energy, direction and most probable mass of the primary cosmic particles, but also to obtain information about the properties of their hadronic interactions at an energy more than one order of magnitude above that accessible with the current highest energy human-made accelerator. In this contribution we will review the state-of-the-art in UHECRs detection. We will present the leading experiments Pierre Auger Observatory and Telescope Array and discuss the cosmic ray energy spectrum, searches for directional anisotropy, studies of mass composition, the determination of the number of shower muons (which is sensitive to the shower hadronic interactions) and the proton-air cross section.
We present a simple formula for the average dark energy equation of state at redshifts between those of two observations of baryon acoustic oscillations (BAO). The formula is independent of any parametrization or basis of the dark energy equation of state and essentially independent of the cosmological model. We use this formula to study the well-known tension between Lyman alpha forest BAO and other cosmological probes. Using only the line of sight Lyman alpha forest BAO and BOSS CMASS dataset, there is already more than 2 sigma tension with the standard LambdaCDM cosmological model which implies that either (i) The BOSS Lyman alpha forest measurement of the Hubble parameter was too low as a result of a statistical fluctuation or systematic error or else (ii) the dark energy equation of state falls steeply at high redshift.
The James Webb Space Telescope (JWST) has the potential to enhance our understanding of near-Earth objects (NEOs). We present results of investigations into the observability of NEOs given the nominal observing requirements of JWST on elongation (85-135 degrees) and non-sidereal rates ($<$30mas/s). We find that approximately 75% of NEOs can be observed in a given year. However, observers will need to wait for appropriate observing windows. We find that JWST can easily execute photometric observations of meter-sized NEOs which will enhance our understanding of the small NEO population.
Giant planets are believed to host central dense rocky/icy cores that are key actors in the core-accretion scenario for their formation. In the same time, some of their components are unstable in the temperature and pressure regimes of central regions of giant planets and only ab-initio EOS computations can address the question of the state of matter. In this framework, several works demonstrated that erosion and redistribution of core materials in the envelope must be taken into account. These complex mechanisms thus may deeply modify giant planet interiors for which signatures of strong tidal dissipation have been obtained for Jupiter and Saturn. The best candidates to explain this dissipation are the viscoelastic dissipation in the central dense core and turbulent friction acting on tidal inertial waves in their fluid convective envelope. In this work, we study the consequences of the possible melting of central regions for the efficiency of each of these mechanisms.
A scaling theory of long-wavelength electrostatic turbulence in a magnetised, weakly collisional plasma (e.g., drift-wave turbulence driven by temperature gradients) is proposed, with account taken both of the nonlinear advection of the perturbed particle distribution by fluctuating ExB flows and of its phase mixing, which is caused by the streaming of the particles along the mean magnetic field and, in a linear problem, would lead to Landau damping. A consistent theory is constructed in which very little free energy leaks into high velocity moments of the distribution, rendering the turbulent cascade in the energetically relevant part of the wave-number space essentially fluid-like. The velocity-space spectra of free energy expressed in terms of Hermite-moment orders are steep power laws and so the free-energy content of the phase space does not diverge at infinitesimal collisionality (while it does for a linear problem); collisional heating due to long-wavelength perturbations vanishes in this limit (also in contrast with the linear problem, in which it occurs at the finite rate equal to the Landau-damping rate). The ability of the free energy to stay in the low velocity moments of the distribution is facilitated by the "anti-phase-mixing" effect, whose presence in the nonlinear system is due to the stochastic version of the plasma echo (the advecting velocity couples the phase-mixing and anti-phase-mixing perturbations). The partitioning of the wave-number space between the (energetically dominant) region where this is the case and the region where linear phase mixing wins is governed by the "critical balance" between linear and nonlinear timescales (which for high Hermite moments splits into two thresholds, one demarcating the wave-number region where phase mixing predominates, the other where plasma echo does).
Starting with the relativistic Boltzmann equation for a system of particles defined by a distribution function, we have derived the virial relation for a spherical structure within an expanding background in the context of general relativity. This generalized form of the virial relation is then applied to the static case of a spherically symmetric structure to see the difference in the simplest case to the Newtonian relation. A relativistic Mass-Temperature relation for this simple case is also derived which can be applied to compact objects in astrophysics. Our general virial relation is then applied to the non-static case of a structure within an expanding universe where an extra term, usually missed in studies of structures in the presence of the dark energy, appears.
In the energy range from few TeV to 25 TeV, upper bounds on the dark matter decay rate into high energy monochromatic neutrinos have recently become comparable to those on monochromatic gamma-ray lines. This implies clear possibilities of a future double "smoking-gun" evidence for the dark matter particle, from the observation of both a gamma and a neutrino line at the same energy. In particular, we show that a scenario where both lines are induced from the same dark matter particle decay leads to correlations that can already be tested. We study this "double monochromatic" scenario by considering the complete list of lowest dimensional effective operators that could induce such a decay. Furthermore, we argue that, on top of lines from decays into two-body final states, three-body final states can also be highly relevant. In addition to producing a distinct hard photon spectrum, three-body final states also produce a line-like feature in the neutrino spectrum that can be searched for by neutrino telescopes.
A model of dark sector where $O({\rm few~GeV})$ mass dark matter particles $\chi$ are supplied by a lighter dark force mediator $V$, $m_V \ll m_\chi$, is motivated by the recently discovered mismatch between simulated and observed shapes of galactic haloes. Such models, in general, provide a challenge for direct detection efforts and collider searches. We show that for a large range of coupling constants and masses, the production and decay of the bound states of $\chi$, such as $0^{-+}$ and $1^{--}$ states, $\eta_D$ and $ \Upsilon_D$, is an important search channel. We show that $e^+e^-\to \eta_D +V$ or $\Upsilon_D +\gamma$ production at $B$-factories for $\alpha_D > 0.1$ is sufficiently strong to result in multiple pairs of charged leptons and pions via $\eta_D\to 2V \to 2(l^+l^-)$ and $\Upsilon_D\to 3V \to 3(l^+l^-)$ $(l=e,\mu,\pi)$. The absence of such final states in the existing searches performed at BaBar and Belle sets new constraints on the parameter space of the model. We also show that a search for multiple bremsstrahlung of dark force mediators, $e^+e^-\to \chi\bar\chi+nV$, resulting in missing energy and multiple leptons, will further improve the sensitivity to self-interacting dark matter.
We give a review on the SIMP paradigm and discuss a consistent model for SIMP dark mesons in the context of a dark QCD with flavor symmetry. The $Z'$-portal interaction is introduced being compatible with stable dark mesons and is responsible for making the SIMP dark mesons remain in kinetic equilibrium with the SM during the freeze-out process. The SIMP parameter space of the $Z'$ gauge boson can be probed by future collider and direct detection experiments.
In general-relativistic perturbation theory, a point mass accelerates away
from geodesic motion due to its gravitational self-force. Because the
self-force is small, one can often approximate the motion as geodesic. However,
it is well known that self-force effects accumulate over time, making the
geodesic approximation fail on long timescales. It is less well known that this
failure at large times translates to a failure at large distances as well. At
second perturbative order, two large-distance pathologies arise: spurious
secular growth and infrared-divergent retarded integrals. Both stand in the way
of practical computations of second-order self-force effects.
Utilizing a simple flat-space scalar toy model, I develop methods to overcome
these obstacles. The secular growth is tamed with a multiscale expansion that
captures the system's slow evolution. The divergent integrals are eliminated by
matching to the correct retarded solution at large distances. I also show how
to extract conservative self-force effects by taking local-in-time "snapshots"
of the global solution. These methods are readily adaptable to the physically
relevant case of a point mass orbiting a black hole.
We propose a novel leptogenesis scenario at the reheating era. Our setup is minimal in the sense that, in addition to the standard model Lagrangian, we only consider an inflaton and higher dimensional operators. The lepton number asymmetry is produced not by the decay of a heavy particle, but by the scattering between the standard model particles. After the decay of an inflaton, the model is described within the standard model with higher dimensional operators. The Sakharov's three conditions are satisfied by the following way. The violation of the lepton number is realized by the dimension-5 operator. The complex phase comes from the dimension-6 four lepton operator. The universe is out of equilibrium before the reheating is completed. It is found that the successful baryogenesis is realized for the wide range of parameters, the inflaton mass and reheating temperature, depending on the cutoff scale. Since we only rely on the effective Lagrangian, our scenario can be applicable to all mechanisms to generate neutrino Majorana masses.
We report constraints on models of composite inflation based on the slow-roll approximation by using the recent Planck measurement. In so doing, we compare the spectral index of curvature perturbation and the tensor-to-scalar ratio predicted by such models with 2015 Planck data. We find that the results predicted by the models present in this work are still consistent with the Planck analysis.
Using a three-dimensional nonhydrostatic general circulation model, we investigate the response of the thermosphere-ionosphere system to the 5-6 August 2011 major geomagnetic storm. The model is driven by measured storm-time input data of the Interplanetary Magnetic Field (IMF), solar activity, and auroral activity. Simulations for quiet steady conditions over the same period are performed as well in order to assess the response of the neutral and plasma parameters to the storm. During the storm, the high-latitude mean ion flows are enhanced by up to 150-180%. Largest ion flows are found in the main phase of the storm. Overall, the global mean neutral temperature increases by up to 15%, while the maximum thermal response is higher in the winter Southern Hemisphere at high-latitudes than the summer Northern Hemisphere: 40% vs. 20%increase in high-latitude mean temperature, respectively. The global mean Joule heating increases by more than a factor of three. There are distinct hemispheric differences in the magnitude and morphology of the horizontal ion flows and thermospheric flows during the different phases of the storm. The largest hemispheric difference in the thermospheric circulation is found during the main and recovery phases of the storm, demonstrating appreciable geographical variations. The advective forcing is found to contribute to the modeled hemispheric differences.
Oks proposes the existence of stable planetary orbits around binary stars, in the shape of a helix on a conical surface whose axis of symmetry coincides with the interstellar axis. We show that planetary orbits initially meeting this description will not continue to do so as the binary pair rotates.
We study the metric perturbations in the context of restricted $f(R)$ gravity, in which a parameter for deviation from the full diffeomorphisms of space-time is introduced. We demonstrate that one can choose the parameter to remove the induced anisotropic stress, which is present in the usual $f(R)$ gravity. Moreover, to prevent instability for the scalar and tensor metric perturbations, some constraints on the model are obtained.
We argue that the exact degeneracy of vacua in N=1 supergravity can shed light on the smallness of the cosmological constant. The presence of such vacua, which are degenerate to very high accuracy, may also result in small values of the quartic Higgs coupling and its beta function at the Planck scale in the phase in which we live.
This paper describes several applications in astronomy and cosmology that are addressed using probabilistic modelling and statistical inference.
It has recently been shown by Yang. et. al. [Phys. Rev. D {\bf 87}, 041502(R) (2013)] that rotating Kerr black holes are characterized by two distinct sets of quasinormal resonances. These two families of quasinormal resonances display qualitatively different asymptotic behaviors in the extremal ($a/M\to 1$) black-hole limit: The zero-damping modes (ZDMs) are characterized by relaxation times which tend to infinity in the extremal black-hole limit ($\Im\omega\to 0$ as $a/M\to 1$), whereas the damped modes (DMs) are characterized by non-zero damping rates ($\Im\omega\to$ finite-values as $a/M\to 1$). In this paper we refute the claim made by Yang et. al. that co-rotating DMs of near-extremal black holes are restricted to the limited range $0\leq \mu\lesssim\mu_{\text{c}}\approx 0.74$, where $\mu\equiv m/l$ is the dimensionless ratio between the azimuthal harmonic index $m$ and the spheroidal harmonic index $l$ of the perturbation mode. In particular, we use an analytical formula originally derived by Detweiler in order to prove the existence of DMs (damped quasinormal resonances which are characterized by finite $\Im\omega$ values in the $a/M\to 1$ limit) of near-extremal black holes in the $\mu>\mu_{\text{c}}$ regime, the regime which was claimed by Yang et. al. not to contain damped modes. We show that these co-rotating DMs (in the regime $\mu>\mu_{\text{c}}$) are expected to characterize the resonance spectra of rapidly-rotating (near-extremal) black holes with $a/M\gtrsim 1-10^{-9}$.
Links to: arXiv, form interface, find, astro-ph, recent, 1510, contact, help (Access key information)
We present Submillimeter Array (SMA) molecular line observations in two 2
GHz-wide bands centered at 217.5 and 227.5 GHz, toward the massive star forming
region W51 North. We identified 84 molecular line transitions from 17 species
and their isotopologues. The molecular gas distribution of these lines mainly
peaks in the continuum position of W51 North, and has a small tail extending to
the west, probably associated with W51 d2. In addition to the commonly detected
nitrogen and oxygen-bearing species, we detected a large amount of transitions
of the Acetone (CH$_3$COCH$_3$) and Methyl Formate (CH$_3$OCHO), which may
suggest that these molecules are present in an early evolutionary stage of the
massive stars. We also found that W51 North is an ethanol-rich source. There is
no obvious difference in the molecular gas distributions between the
oxygen-bearing and nitrogen-bearing molecules. Under the assumption of Local
Thermodynamic Equilibrium (LTE), with the XCLASS tool, the molecular column
densities, and rotation temperatures are estimated.
We found that the oxygen-bearing molecules have considerable higher column
densities and fractional abundances than the nitrogen-bearing molecules. The
rotation temperatures range from 100 to 200 K, suggesting that the molecular
emission could be originated from a warm environment.
Finally, based on the gas distributions, fractional abundances and the
rotation temperatures, we conclude that CH$_3$OH, C$_2$H$_5$OH, CH$_3$COCH$_3$
and CH$_3$CH$_2$CN might be synthesized on the grain surface, while gas phase
chemistry is responsible for the production of CH$_3$OCH$_3$, CH$_3$OCHO and
CH$_2$CHCN.
We have determined the detailed star formation history and total mass of the globular clusters in the Fornax dwarf spheroidal using archival HST WFPC2 data. Colour magnitude diagrams are constructed in the F555W and F814W bands and corrected for the effect of Fornax field star contamination, after which we use the routine Talos to derive the quantitative star formation history as a function of age and metallicity. The star formation history of the Fornax globular clusters shows that Fornax 1, 2, 3 and 5 are all dominated by ancient~(>10 Gyr) populations. Cluster Fornax 1 and 3 display metallicities as low as [Fe/H]=-2.5 while Fornax 2 and 5 are slightly more metal-rich at [Fe/H]=-2.0, consistent with resolved and unresolved metallicity tracers. Conversely, Fornax 4 displays a more extended star formation history dominated by metal-rich~([Fe/H]=-1.4 dex) stars, with an age of ~10 Gyr, inconsistent with the other clusters. Its central location and complex population mix favours the proposed model that it might be the nucleus of the Fornax dwarf spheroidal. The combined stellar mass in globular clusters as derived from the SFH is (9.73$\pm$0.79)$\times$10$^{5}$ M$_{\odot}$ which corresponds to 2.5$\pm$0.2 percent of the total stellar mass in Fornax. This mass can be further subdivided into metal-poor stars to yield a mass fraction of 9.4$\pm$1.3 percent of the metal-poor Fornax field, or 19.1$\pm$2.8 percent when considering just the full mass of the four most metal-poor clusters. Therefore, the SFH results provide separate supporting evidence for the unusually high mass fraction of the GCs compared to the Fornax field population.
We study the connection of star formation to atomic (HI) and molecular hydrogen (H$_2$) in isolated, low metallicity dwarf galaxies with high-resolution ($m_{\rm gas}$ = 4 M$_\odot$, $N_{\rm ngb}$ = 100) SPH simulations. The model includes self-gravity, non-equilibrium cooling, shielding from an interstellar radiation field, the chemistry of H$_2$ formation, H$_2$-independent star formation, supernova feedback and metal enrichment. We find that the H$_2$ mass fraction is sensitive to the adopted dust-to-gas ratio and the strength of the interstellar radiation field, while the star formation rate is not. Star formation is regulated by stellar feedback, keeping the gas out of thermal equilibrium for densities $n <$ 1 cm$^{-3}$. Because of the long chemical timescales, the H$_2$ mass remains out of chemical equilibrium throughout the simulation. Star formation is well-correlated with cold ( T $\leqslant$ 100 K ) gas, but this dense and cold gas - the reservoir for star formation - is dominated by HI, not H$_2$. In addition, a significant fraction of H$_2$ resides in a diffuse, warm phase, which is not star-forming. The ISM is dominated by warm gas (100 K $<$ T $\leqslant 3\times 10^4$ K) both in mass and in volume. The scale height of the gaseous disc increases with radius while the cold gas is always confined to a thin layer in the mid-plane. The cold gas fraction is regulated by feedback at small radii and by the assumed radiation field at large radii. The decreasing cold gas fractions result in a rapid increase in depletion time (up to 100 Gyrs) for total gas surface densities $\Sigma_{\rm HI+H_2} \lesssim$ 10 M$_\odot$pc$^{-2}$, in agreement with observations of dwarf galaxies in the Kennicutt-Schmidt plane.
We present the evolution of galaxy sizes, from redshift 2 to 0, for actively star forming and passive galaxies in the cosmological hydrodynamical 1003 cMpc3 simulation of the EAGLE project. We find that the sizes increase with stellar mass , but that the relation weakens with increasing redshift. Separating galaxies by their star formation activity, we find that passive galaxies are typically smaller than active galaxies at fixed stellar mass. These trends are consistent with those found in observations and the level of agreement between the predicted and observed size - mass relation is of order 0.1 dex for z < 1 and 0.2-0.3 dex from redshift 1 to 2. We use the simulation to compare the evolution of individual galaxies to that of the population as a whole. While the evolution of the size-stellar mass relation for active galaxies provides a good proxy for the evolution of individual galaxies, the evolution of individual passive galaxies is not well represented by the observed size - mass relation due to the evolving number density of passive galaxies. Observations of z \approx 2 galaxies have revealed an abundance of massive red compact galaxies, that depletes below z \approx 1. We find that a similar population forms naturally in the simulation. Comparing these galaxies to their z = 0 descendants, we find that all compact galaxies grow in size due to the high-redshift stars migrating outwards. Approximately 60% of the compact galaxies increase in size further due to renewed star formation and/or mergers.
The visibility of LyA emitting galaxies during the Epoch of Reionization is controlled by both diffuse HI patches in large-scale bubble morphology and small-scale absorbers. To investigate the impact on LyA photons, we apply a novel combination of analytic and numerical calculations to three scenarios: (i) the `bubble' model, where only diffuse HI outside ionized bubbles is present; (ii) the `web' model, where HI exists only in overdense self-shielded gas; and (iii) the more realistic 'web-bubble' model, which contains both. Our analysis confirms that there is a degeneracy between the ionization structure of the intergalactic medium (IGM) and the HI fraction inferred from LyA surveys, as the three models suppress LyA flux equally with very different HI fractions. We argue that a joint analysis of the LyA luminosity function and the rest-frame equivalent width distribution/LyA fraction can break this degeneracy and provide constraints on the reionization history and its topology. We further show that constraints can improve if we consider the full shape of the M_UV-dependent redshift evolution of the LyA fraction of Lyman break galaxies. Contrary to conventional wisdom, we find that (i) a drop of LyA fraction larger for UV-faint than for UV-bright galaxies can be reproduced with web and web-bubble models and therefore does not provide exclusive evidence of patchy reionization, and (ii) the IGM-transmission PDF is unimodal for bubble models and bimodal in web models. We further highlight the importance of galaxy-absorber cross-correlation. Comparing our models to observations, the neutral fraction at z~7 is likely to be of order of tens of per cent when interpreted with bubble or web-bubble models. Alternatively, we obtain a conservative lower limit ~1% in the web models, if we allow for a drop in the photoionization rate by a factor of ~100 from the post-reionized universe. [abridged]
We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters ($BVRI$ and $ugriz$). Combined with UV data from the $Hubble$ $Space$ $Telescope$ and $Swift$, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from $1158\,{\rm \AA}$ to the $z$-band ($\sim\! 9160\,{\rm \AA}$). We find that the lags at wavelengths longer than the $V$ band are equal to or greater than the lags of high ionization-state emission lines (such as HeII$\lambda 1640$ and $\lambda 4686$), suggesting that the continuum emitting source is of a physical size comparable to the inner broad line region. The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $\tau \propto \lambda^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the inter-band continuum lags due to broad line region emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the $u$ and $U$ bands due to Balmer continuum emission.
We present an analytic model for how momentum deposition from stellar feedback simultaneously regulates star formation and drives outflows in a turbulent interstellar medium (ISM). Because the ISM is turbulent, a given patch of ISM exhibits sub-patches with a range of surface densities. The high-density patches are 'pushed' by feedback, thereby driving turbulence and self-regulating local star formation. Sufficiently low-density patches, however, are accelerated to above the escape velocity before the region can self-adjust and are thus vented as outflows. In the turbulent-pressure-supported regime, when the gas fraction is $\gtrsim 0.3$, the ratio of the turbulent velocity dispersion to the circular velocity is sufficiently high that at any given time, of order half of the ISM has surface density less than the critical value and thus can be blown out on a dynamical time. The resulting outflows have a mass-loading factor ($\eta \equiv M_{\rm out}/M_{\star}$) that is inversely proportional to the gas fraction times the circular velocity. At low gas fractions, the star formation rate needed for local self-regulation, and corresponding turbulent Mach number, decline rapidly; the ISM is 'smoother', and it is actually more difficult to drive winds with large mass-loading factors. Crucially, our model predicts that stellar-feedback-driven outflows should be suppressed at $z \lesssim 1$ in $M_{\star} \gtrsim 10^{10} M_{\odot}$ galaxies. This mechanism allows massive galaxies to exhibit violent outflows at high redshifts and then 'shut down' those outflows at late times, thereby enabling the formation of a smooth, extended thin stellar disk. We provide simple fitting functions for $\eta$ that should be useful for sub-resolution and semi-analytic models. [abridged]
Empirical methods for connecting galaxies to their dark matter halos have become essential in interpreting measurements of the spatial statistics of galaxies. Among the most successful of these methods is the technique of subhalo abundance matching, which has to date been used to associate galaxy properties with a small set of halo properties. We generalize this set of halo properties to allow variable dependence on halo concentration, and parameterize the degree of concentration dependence with a single parameter. This parameter provides a smooth interpolation between abundance matching to peak halo mass and to peak halo circular velocity. We characterize the influence of this parameter on two-point clustering, the satellite fraction, and the degree of galaxy assembly bias. We also evaluate the degeneracies between the concentration dependence and the scatter in the abundance matching relation. We show that low redshift clustering measurements from SDSS prefer a moderate amount of concentration dependence --- more than would be indicated by matching galaxy luminosity to the peak halo mass, and less than would be indicated by matching to the peak halo circular velocity. We also show that these results are robust to moderate changes in cosmological parameters, and that the best-fit model from two-point clustering agrees with previous measurements of the satellite fraction. We note that statistical constraints on these models have been (and still are, in most regimes) limited primarily by sample variance in the limited-size simulations, and not in the data. We discuss physical interpretations of these results and their implications for the galaxy-halo connection.
Filaments are ubiquitous in the universe. They are seen in cosmological structures, in the Milky Way centre and in dense interstellar gas. Recent observations have revealed that stars and star clusters form preferentially at the intersection of dense filaments. Understanding the formation and properties of filaments is therefore a crucial step in understanding star formation. Here we perform three-dimensional high-resolution magnetohydrodynamical simulations that follow the evolution of molecular clouds and the formation of filaments and stars within them. We apply a filament detection algorithm and compare simulations with different combinations of physical ingredients: gravity, turbulence, magnetic fields and jet/outflow feedback. We find that gravity-only simulations produce significantly narrower filament profiles than observed, while simulations that at least include turbulence produce realistic filament properties. For these turbulence simulations, we find a remarkably universal filament width of (0.10+/-0.02) pc, which is independent of the evolutionary stage or the star formation history of the clouds. We derive a theoretical model that provides a physical explanation for this characteristic filament width, based on the sonic scale (lambda_sonic) of molecular cloud turbulence. Our derivation provides lambda_sonic as a function of the cloud diameter L, the velocity dispersion sigma_v, the gas sound speed c_s and the strength of the magnetic field parameterised by plasma beta. For typical cloud conditions in the Milky Way spiral arms, we find theoretically that lambda_sonic = 0.04-0.16 pc, in excellent agreement with the filament width of 0.05-0.15 pc found in observations.
Using a radio-quiet subsample of the Sloan Digital Sky Survey spectroscopic quasar catalog, spanning redshifts 0.5-3.5, we derive the mean millimetre and far-infrared quasar spectral energy densities via a stacking analysis of Atacama Cosmology Telescope and Herschel-SPIRE data. We constrain the form and evolution of the far-infrared emission finding 3-4$\sigma$ evidence for the presence of the thermal Sunyaev-Zel'dovich (SZ) effect in the millimetre bands. We find this signal to be characteristic of a hot ionized gas component with thermal energy $(6.2 \pm 1.7) \times 10^{60}$erg. This amount of thermal energy is an order of magnitude greater than would be expected assuming only hot gas in virial equilibrium with the dark matter haloes of $(1-5)\times 10^{12}h^{-1}$M$_\odot$ that these systems are expected to occupy, though the highest quasar mass estimates found in the literature could explain a large fraction of this energy. We find that our measurements are consistent with a scenario in which quasars deposit up to $(14.5 \pm 3.3)~\tau_8^{-1}$ per cent of their radiative energy into their circumgalactic environment if their typical period of quasar activity is $\tau_8\times 10^8$ years. If quasar host masses are high ($\sim10^{13}h^{-1}$M$_\odot$), then this percentage will be reduced significantly. Furthermore, the uncertainty quoted for this percentage is only statistical and additional systematic uncertainties (e.g., on quasar bolometric luminosity) enter at the 40 per cent level. Finally, emission from thermal dust is significant in these systems, with infrared luminosities of $\log_{10}(L_{\rm ir}/{\rm L}_\odot)=11.4-12.2$, increasing to higher redshift. We consider various models for dust emission. While sufficiently complex dust models can obviate the SZ effect, the SZ interpretation remains favoured at the 3-4$\sigma$ level for most models.
Studying star formation beyond the optical radius of galaxies allows us to test empirical relations in extreme conditions with low average gas density and low molecular fraction. Previous studies discovered galaxies with extended ultraviolet (XUV) disks, which often contain star forming regions with lower Halpha-to-far-UV (FUV) flux ratios compared to inner disk star forming regions. However, most previous studies lack measurements of molecular gas, which is presumably the component of the interstellar medium out of which stars form. We analyzed published CO measurements and upper limits for fifteen star forming regions in the XUV or outer disk of three nearby spiral galaxies and a new CO upper limit from the IRAM 30 m telescope in one star forming region at r = 3.4 r_25 in the XUV disk of NGC 4625. We found that the star forming regions are in general consistent with the same molecular-hydrogen Kennicutt-Schmidt law that applies within the optical radius, independent of whether we used Halpha or FUV as the star formation rate (SFR) tracer. However, a number of the CO detections are significantly offset towards higher SFR surface density for their molecular hydrogen surface density. Deeper CO data may enable us to use the presence or absence of molecular gas as an evolutionary probe to break the degeneracy between age and stochastic sampling of the initial mass function as the explanation for the low Halpha-to-FUV flux ratios in XUV disks.
The analysis of galaxy properties and the relations among them and the environment, can be used to investigate the physical processes driving galaxy evolution. We study the cluster A209 by using the CLASH-VLT spectroscopic data combined with Subaru photometry, yielding to 1916 cluster members down to a stellar mass of 10^{8.6} Msun. We determine: i) the stellar mass function of star-forming and passive galaxies; ii) the intra-cluster light and its properties; iii) the orbits of low- and high-mass passive galaxies; and iv) the mass-size relation of ETGs. The stellar mass function of the star-forming galaxies does not depend on the environment, while the slope found for passive galaxies becomes flatter in the densest region. The color distribution of the intra-cluster light is consistent with the color of passive members. The analysis of the dynamical orbits shows that low-mass passive galaxies have tangential orbits, avoiding small pericenters around the BCG. The mass-size relation of low-mass passive ETGs is flatter than that of high mass galaxies, and its slope is consistent with that of field star-forming galaxies. Low-mass galaxies are also more compact within the scale radius of 0.65 Mpc. The ratio between stellar and number density profiles shows a mass segregation in the center. The comparative analysis of the stellar and total density profiles indicates that this effect is due to dynamical friction. Our results are consistent with a scenario in which the "environmental quenching" of low-mass galaxies is due to mechanisms such as harassment out to R200, starvation and ram-pressure stripping at smaller radii, as supported by the analysis of the mass function, of the dynamical orbits and of the mass-size relation of passive early-types in different regions. Our analyses support the idea that the intra-cluster light is formed through the tidal disruption of subgiant galaxies.
Recent imaging of protoplanetary disks with high resolution and contrast have revealed a striking variety of substructure. Of particular interest are cases where near-infrared scattered light images show evidence for low-intensity annular "gaps". The origins of such structures are still uncertain, but the interaction of the gas disk with planets is a common interpretation. We study the impact that the evolution of the solid material can have on the observable properties of disks in a simple scenario without any gravitational or hydrodynamical disturbances to the gas disk structure. Even with a smooth and continuous gas density profile, we find that the scattered light emission produced by small dust grains can exhibit ring-like depressions similar to those presented in recent observations. The physical mechanisms responsible for these features rely on the inefficient fragmentation of dust particles. The occurrence and position of the proposed "gap" features depend most strongly on the dust-to-gas ratio, the fragmentation threshold velocity, the strength of the turbulence, and the age of the disk, and should be generic (at some radius) for typically adopted disk parameters. The same physical processes can affect the thermal emission at optically thin wavelengths ($\sim$1 mm), although the behavior can be more complex; unlike for disk-planet interactions, a "gap" should not be present at these longer wavelengths.
Locating planets in circumstellar Habitable Zones is a priority for many exoplanet surveys. Space-based and ground-based surveys alike require robust toolsets to aid in target selection and mission planning. We present the Catalog of Earth-Like Exoplanet Survey Targets (CELESTA), a database of Habitable Zones around 36,000 nearby stars. We calculated stellar parameters, including effective temperatures, masses, and radii, and we quantified the orbital distances and periods corresponding to the circumstellar Habitable Zones. We gauged the accuracy of our predictions by contrasting CELESTA's computed parameters to observational data. We ascertain a potential return on investment by computing the number of Habitable Zones probed for a given survey duration. A versatile framework for extending the functionality of CELESTA into the future enables ongoing comparisons to new observations, and recalculations when updates to Habitable Zone models, stellar temperatures, or parallax data become available. We expect to upgrade and expand CELESTA using data from the Gaia mission as the data becomes available.
The dust-content of damped Lyman-alpha systems (DLAs) is an important observable for understanding their origin and the neutral gas reservoirs of galaxies. While the average colour-excess of DLAs, E(B-V), is known to be <15 milli-magnitudes (mmag), both detections and non-detections with ~2 mmag precision have been reported. Here we find 3.2-sigma statistical evidence for DLA dust-reddening of 774 Sloan Digital Sky Survey (SDSS) quasars by comparing their fitted spectral slopes to those of ~7000 control quasars. The corresponding E(B-V) is 3.0 +/- 1.0 mmag, assuming a Small Magellanic Cloud (SMC) dust extinction law, and it correlates strongly (3.5-sigma) with the metal content, characterised by the SiII1526 absorption-line equivalent width, providing additional confidence that the detection is due to dust in the DLAs. Evolution of E(B-V) over the redshift range 2.1 < z < 4.0 is limited to <2.5 mmag per unit redshift (1-sigma), consistent with the known, mild DLA metallicity evolution. There is also no apparent relationship with neutral hydrogen column density, N(HI), though the data are consistent with a mean E(B-V)/N(HI) = (3.5 +/- 1.0) x 10^{-24} mag cm^2, approximately the ratio expected from the SMC scaled to the lower metallicities typical of DLAs. We implement the SDSS selection algorithm in a portable code to assess the potential for systematic, redshift-dependent biases stemming from its magnitude and colour-selection criteria. The effect on the mean E(B-V) is negligible (<5 per cent) over the entire redshift range of interest. Given the broad potential usefulness of this implementation, we make it publicly available.
The formation of the Milky Way stellar halo is thought to be the result of merging and accretion of building blocks such as dwarf galaxies and massive globular clusters. Recently, Deason et al. (2015) suggested that the Milky Way outer halo formed mostly from big building blocks, such as dwarf spheroidal galaxies, based on the similar number ratio of blue straggler (BS) stars to blue horizontal-branch (BHB) stars. Here we demonstrate, however, that this result is seriously biased by not taking into detailed consideration on the formation mechanism of BHB stars from helium enhanced second-generation population. In particular, the high BS-to-BHB ratio observed in the outer halo fields is most likely due to a small number of BHB stars provided by GCs rather than to a large number of BS stars. This is supported by our dynamical evolution model of GCs which shows preferential removal of first generation stars in GCs. Moreover, there are a sufficient number of outer halo GCs which show very high BS-to-BHB ratio. Therefore, the BS-to-BHB number ratio is not a good indicator to use in arguing that more massive dwarf galaxies are the main building blocks of the Milky Way outer halo. Several lines of evidence still suggest that GCs can contribute a significant fraction of the outer halo stars.
We investigate the formation of metal-poor globular clusters (GCs) at the center of two dark matter halos with $M_{halo}\sim4\times10^7\,M_\odot$ at $z>10$ using cosmological radiation-hydrodynamics simulations. We find that very compact ($\lesssim$ 1 pc) and massive ($\sim6\times10^5\,M_\odot$) clusters form rapidly when pristine gas collapses isothermally with the aid of efficient Ly$\alpha$ emission during the transition from molecular-cooling halos to atomic-cooling halos. Because the local free-fall time of dense star-forming gas is very short ($\ll 1\,{\rm Myr}$), a large fraction of the collapsed gas is turned into stars before stellar feedback processes blow out the gas and shut down star formation. Although the early stage of star formation is limited to a small region of the central star-forming disk, we find that the disk quickly fragments due to metal enrichment from supernovae. Sub-clusters formed in the fragmented clouds eventually merge with the main cluster at the center. We estimate using a simple analytic calculation that, if 20 percent of these halos form a GC, they can account for the number of metal-poor GCs observed in the local Universe, making this scenario appealing. However, despite the similarities in mass, star formation histories, and size to the local GCs, we find that there is a substantial spread in metallicities within each simulated GC candidate, as metals are enriched inhomogeneously in star-forming clouds. We discuss a possible solution, involving Pop III stars, to the metal enrichment problem.
From a comprehensive statistical analysis of $Swift$ X-ray light-curves of gamma-ray bursts (GRBs) collected from December 2004 to the end of 2010, we found a three-parameter correlation between the isotropic energy emitted in the rest frame 1-10$^4$ keV energy band during the prompt emission (E$_{\rm{\gamma,iso}}$), the rest frame peak of the prompt emission energy spectrum (E$_{\rm{pk}}$), and the X-ray energy emitted in the rest frame 0.3-30 keV observed energy band (E$_{\rm{X,iso}}$), computed excluding the contribution of the flares. In this paper, we update this correlation with the data collected until June 2014, expanding the sample size with $\sim$35% more objects, where the number of short GRBs doubled. With this larger sample we confirm the existence of a universal correlation that connects the prompt and afterglow properties of long and short GRBs. We show that this correlation does not depend on the X-ray light-curve morphology and that further analysis is necessary to firmly exclude possible biases derived by redshift measurements. In addition we discuss about the behavior of the peculiar objects as ultra-long GRBs and we propose the existence of an intermediate group between long and short GRBs. Interestingly, two GRBs with uncertain classification fall into this category. Finally, we discuss the physics underlying this correlation, in the contest of the efficiency of conversion of the prompt $\gamma$-ray emission energy into the kinetic energy of the afterglow, the photosferic model, and the cannonball model.
The stereoscopic imaging atmospheric Cherenkov technique, developed in the 1980s and 1990s, is now used by a number of existing and planned gamma-ray observatories around the world. It provides the most sensitive view of the very high energy gamma-ray sky (above 30 GeV), coupled with relatively good angular and spectral resolution over a wide field-of-view. This Chapter summarizes the details of the technique, including descriptions of the telescope optical systems and cameras, as well as the most common approaches to data analysis and gamma-ray reconstruction.
We use the CARMA millimeter interferometer to map the Antennae Galaxies (NGC4038/39), tracing the bulk of the molecular gas via the 12CO(1-0) line and denser molecular gas via the high density transitions HCN(1-0), HCO+(1-0), CS(2-1), and HNC(1-0). We detect bright emission from all tracers in both the two nuclei and three locales in the overlap region between the two nuclei. These three overlap region peaks correspond to previously identified "supergiant molecular clouds". We combine the CARMA data with Herschel infrared (IR) data to compare observational indicators of the star formation efficiency (SFR/H2~IR/CO), dense gas fraction (HCN/CO), and dense gas star formation efficiency (IR/HCN). Regions within the Antennae show ratios consistent with those seen for entire galaxies, but these ratios vary by up to a factor of 6 within the galaxy. The five detected regions vary strongly in both their integrated intensities and these ratios. The northern nucleus is the brightest region in mm-wave line emission, while the overlap region is the brightest part of the system in the IR. We combine the CARMA and Herschel data with ALMA CO data to report line ratio patterns for each bright point. CO shows a declining spectral line energy distribution, consistent with previous studies. HCO+(1-0) emission is stronger than HCN(1-0) emission, perhaps indicating either more gas at moderate densities or higher optical depth than is commonly seen in more advanced mergers.
The longevity of Cassini's exploration of Saturn's atmosphere (a third of a Saturnian year) means that we have been able to track the seasonal evolution of atmospheric temperatures, chemistry and cloud opacity over almost every season, from solstice to solstice and from perihelion to aphelion. Cassini has built upon the decades-long ground-based record to observe seasonal shifts in atmospheric temperature, finding a thermal response that lags behind the seasonal insolation with a lag time that increases with depth into the atmosphere, in agreement with radiative climate models. Seasonal hemispheric contrasts are perturbed at smaller scales by atmospheric circulation, such as belt/zone dynamics, the equatorial oscillations and the polar vortices. Temperature asymmetries are largest in the middle stratosphere and become insignificant near the radiative-convective boundary. Cassini has also measured southern-summertime asymmetries in atmospheric composition, including ammonia (the key species for the topmost clouds), phosphine and para-hydrogen (both disequilibrium species) in the upper troposphere; and hydrocarbons deriving from the UV photolysis of methane in the stratosphere (principally ethane and acetylene). These chemical asymmetries are now altering in subtle ways due to (i) the changing chemical efficiencies with temperature and insolation; and (ii) vertical motions associated with large-scale overturning in response to the seasonal temperature contrasts. Similarly, hemispheric contrasts in tropospheric aerosol opacity and coloration that were identified during the earliest phases of Cassini's exploration have now reversed, suggesting an intricate link between the clouds and the temperatures. [Abridged]
In this paper we continue the release of high-level data products from the multiyear photometric survey collected at the 67/92 cm Schmidt Telescope in Asiago. The primary goal of the survey is to discover and to characterise variable objects and exoplanetary transits in four fields containing five nearby open clusters spanning a broad range of ages. This second paper releases a photometric catalogue, in five photometric bands, of the Solar-age, Solar-metallicity open cluster M 67 (NGC 2682). Proper motions are derived comparing the positions observed in 2013 at the Asiago's Schmidt Telescope with those extracted from WFI@2.2m MPG/ESO images in 2000. We also analyse the variable sources within M 67. We detected 68 variables, 43 of which are new detection. Variable periods and proper-motion memberships of a large majority of sources in our catalogue are improved with respect to previous releases. The entire catalogue will be available in electronic format. Besides the general interest on an improved catalogue, this work will be particularly useful because of: (1) the imminent release of Kepler/K2 Campaign-5 data of this cluster, for which our catalogue will provide an excellent, high spatial resolution input list, and (2) characterisation of the M 67 stars which are targets of intense HARPS and HARPS-N radial-velocity surveys for planet search.
Ultra-precise light curves from Kepler provide the best opportunity to determine rates and statistical properties of stellar flares. From 11 months of data on the active M4 dwarf, GJ 1243, we have built the largest catalog of flares for a single star: over 6100 events. Combining 885 of our most pristine flares, we generated an empirical white-light flare template. This high-fidelity template shows a rapid initial rise, and two distinct exponential cooling phases. This template is useful in constraining flare energies and for improved flare detection in many surveys. Complex, multi-peaked events are more common for higher energy flares in this sample. Using our flare template we characterize the structure of complex events. In this contributed talk, I presented results from our boutique study of GJ 1243, and described an expanded investigation of the structure of complex flares and their connection to solar events.
Nuclear star clusters (NCs) are found to exist in the centres of many galaxies and appear to follow scaling relations similar to those of super-massive black holes. Previous analytical work has suggested that such relations are a consequence of feedback regulated growth. We explore this idea using high resolution hydrodynamical simulations, focusing on the validity of the simplifying assumptions made in analytical models. In particular, we investigate feedback emanating from multiple stellar sources rather than from a single source, as is usually assumed, and show that collisions betweens shells of gas swept up by feedback leads to momentum cancellation and the formation of high density clumps and filaments. This high density material is resistant both to expulsion from the galaxy potential and to disruption by feedback; if it falls back onto the NC, we expect the gas to be available for further star formation or for feeding a central black hole. We also note our results may have implications for the evolution of globular clusters and stellar clusters in high redshift dark matter halos.
We present the first large scale high angular resolution survey of ionized nitrogen in the Galactic Plane through emission of its two fine structure transitions ([NII]) at 122 $\mu$m and 205 $\mu$m. The observations were largely obtained with the PACS instrument onboard the Herschel Space Observatory. The lines-of-sight were in the Galactic plane, following those of the Herschel OTKP project GOT C+. Both lines are reliably detected at the 10$^{-8}$ - 10$^{-7}$ $W$m$^{-2}$sr$^{-1}$ level over the range -60$^{o}$ $\leq$ $l$ $\leq$ 60$^{o}$. The $rms$ of the intensity among the 25 PACS spaxels of a given pointing is typically less than one third of the mean intensity, showing that the emission is extended. [NII] is produced in gas in which hydrogen is ionized, and collisional excitation is by electrons. The ratio of the two fine structure transitions provides a direct measurement of the electron density, yielding $n(e)$ largely in the range 10 to 50 cm$^{-3}$ with an average value of 29 cm$^{-3}$ and N$^+$ column densities 10$^{16}$ to 10$^{17}$ cm$^{-2}$. [NII] emission is highly correlated with that of [CII], and we calculate that between 1/3 and 1/2 of the [CII] emission is associated with the ionized gas. The relatively high electron densities indicate that the source of the [NII] emission is not the Warm Ionized Medium (WIM), which has electron densities more than 100 times smaller. Possible origins of the observed [NII] include the ionized surfaces of dense atomic and molecular clouds, the extended low density envelopes of HII regions, and low-filling factor high-density fluctuations of the WIM.
We present new near-infrared (NIR) observations of M63 from the Extended Disk Galaxy Exploration Science (EDGES) Survey. The extremely deep 3.6 $\mu$m mosaic reaches 29 AB mag arcsec$^{-2}$ at the outer reaches of the azimuthally-averaged surface brightness profile. At this depth the consequences of galactic accretion are found within a nearby tidal stream and an up-bending break in the slope of the surface brightness profile. This break occurs at a semi-major axis length of $\sim$8', and is evidence of either an enhanced outer disc or an inner stellar halo. Simulations of galaxy evolution, along with our observations, support an inner halo as the explanation for the up-bending break. The mass of this halo component is the largest found in an individual galaxy thus far. Additionally, our observations detect a nearby tidal stream. The mass of the stream suggests that a handful of such accretion events are necessary to populate the inner stellar halo. We also find that the accretion rate of the galaxy from the stream alone underestimates the accretion rate required to build M63's inner stellar halo.
The statistical description of Giant Molecular Cloud (GMC) properties relies heavily on the performance of automatic identification algorithms, which are often seriously affected by the survey design. The algorithm we designed, SCIMES (Spectral Clustering for Interstellar Molecular Emission Segmentation), is able to overcome some of these limitations by considering the cloud segmentation problem in the broad framework of the graph theory. The application of the code on the CO(3-2) High Resolution Survey (COHRS) data allowed for a robust decomposition of more than 12,000 objects in the Galactic Plane. Together with the wealth of Galactic Plane surveys of the recent years, this approach will help to open the door to a future, systematic cataloging of all discrete molecular features of our own Galaxy.
We examine the influence of a continuous density structuring transverse to coronal slabs on the dispersive properties of fundamental standing kink and sausage modes supported therein. We derive generic dispersion relations (DRs) governing linear fast waves in pressureless straight slabs with general transverse density distributions, and focus on the cases where the density inhomogeneity takes place in a layer of arbitrary width and in arbitrary form. The physical relevance of the solutions to the DRs is demonstrated by the corresponding time-dependent computations. For all profiles examined, the lowest-order kink modes are trapped regardless of longitudinal wavenumber $k$. A continuous density distribution introduces a difference to their periods of $\lesssim 13\%$ when $k$ is the observed range, relative to the case where the density profile takes a step-function form. Sausage modes and other branches of kink modes are leaky at small $k$, and their periods and damping times are heavily influenced by how the transverse density profile is prescribed, the lengthscale in particular. These modes have sufficiently high quality to be observable only for physical parameters representative of flare loops. We conclude that while the simpler DR pertinent to a step-function profile can be used for the lowest-order kink modes, the detailed information on the transverse density structuring needs to be incorporated into studies of sausage modes and higher-order kink modes.
Supernova `Refsdal', multiply imaged by cluster MACS1149.5+2223, represents a rare opportunity to make a true blind test of model predictions in extragalactic astronomy, on a time scale that is short compared to a human lifetime. In order to take advantage of this event, we produced seven gravitational lens models with five independent methods, based on Hubble Space Telescope (HST) Hubble Frontier Field images, along with extensive spectroscopic follow-up from HST and from the Very Large Telescope. We compare the model predictions and show that they agree reasonably well with the measured time delays and magnification ratios between the known images, even though these quantities were not used as input. This agreement is encouraging, considering that the models only provide statistical uncertainties, and do not include additional sources of uncertainties such as structure along the line of sight, cosmology, and the mass sheet degeneracy. We then present the model predictions for the other appearances of SN `Refsdal'. A future image will reach its peak in the first half of 2016, while another image appeared between 1994 and 2004. The past image would have been too faint to be detected in archival images. The future image should be approximately one third as bright as the brightest known images and thus detectable in HST images, as soon as the cluster can be targeted again (beginning 2015 October 30). We will find out soon whether our predictions are correct.
Rotation curves of more than one hundred spiral galaxies were compiled from the literature, and deconvolved into bulge, disk, and dark halo using $\chi^2$ fitting in order to determine their scale radii and masses. Correlation analyses were obtained of the fitting parameters for galaxies that satisfied selection and accuracy criteria. Size-mass relations indicate that the sizes and masses are positively correlated among different components in such a way that the larger or more massive is the dark halo, the larger or more massive are the disk and bulge. Empirical size-mass relations were obtained for bulge, disk and dark halo by the least-squares fitting. The disk-to-halo mass ratio was found to be systematically greater by a factor of three than that predicted by cosmological simulations combined with photometry. A preliminary mass function for dark halo was obtained, which is represented by the Schechter function followed by a power law.
We report the discovery of HATS-17b, the first transiting warm Jupiter of the HATSouth network. HATS-17b transits its bright (V=12.4) G-type (M$_{\star}$=1.131 $\pm$ 0.030 M$_{\odot}$, R$_{\star}$=1.091$^{+0.070}_{-0.046}$ R$_{\star}$) metal-rich ([Fe/H]=+0.3 dex) host star in a circular orbit with a period of P=16.2546 days. HATS-17b has a very compact radius of 0.777 $\pm$ 0.056 R$_J$ given its Jupiter-like mass of 1.338 $\pm$ 0.065 M$_J$. Up to 50% of the mass of HATS-17b may be composed of heavy elements in order to explain its high density with current models of planetary structure. HATS-17b is the longest period transiting planet discovered to date by a ground-based photometric survey, and is one of the brightest transiting warm Jupiter systems known. The brightness of HATS-17b will allow detailed follow-up observations to characterize the orbital geometry of the system and the atmosphere of the planet.
Galaxy mergers produce binaries of supermassive black holes, which emit gravitational waves prior to their coalescence. We perform three-dimensional hydrodynamic simulations to study the tidal disruption of stars by such a binary in the final centuries of its life. We find that the gas stream of the stellar debris moves chaotically in the binary potential and forms accretion disks around both black holes. The accretion light curve is modulated over the binary orbital period owing to relativistic beaming. This periodic signal allows to detect the decay of the binary orbit due to gravitational wave emission by observing two tidal disruption events that are separated by more than a decade.
Accreting, steadily nuclear-burning white dwarfs are associated with so-called close-binary supersoft X-ray sources (SSSs), observed to have temperatures of a few$\times 10^{5}$K and luminosities on the order of $10^{38}$erg/s. These and other types of SSSs are expected to be capable of ionizing their surrounding circumstellar medium, however, to date only one such nebula was detected in the Large Magellanic Cloud (of its 6 known close-binary SSSs), surrounding the accreting, nuclear-burning WD CAL 83. This has led to the conclusion that most SSSs cannot have been both luminous ($\gtrsim 10^{37}$erg/s) and hot ($\gtrsim$ few $\times 10^{4}$K) for the majority of their past accretion history, unless the density of the ISM surrounding most sources is much less than that inferred for the CAL 83 nebula (4--10$\rm{cm}^{-3}$). Here we demonstrate that most SSSs must lie in much lower density media than CAL 83. Past efforts to detect such nebulae have not accounted for the structure of the ISM in star-forming galaxies and, in particular, for the fact that most of the volume is occupied by low density warm \& hot ISM. CAL 83 appears to lie in a region of ISM which is at least $\sim 40$-fold overdense. We compute the probability of such an event to be $\approx 18\%$, in good agreement with observed statistics. We provide a revised model for the "typical" SSS nebula, and outline the requirements of a survey of the Magellanic clouds which could detect the majority of such objects. We then briefly discuss some of the possible implications, should there prove to be a large population of previously undiscovered ionizing sources.
We have estimated a metallicity map of the Large Magellanic Cloud (LMC) using the Magellanic Cloud Photometric Survey (MCPS) and Optical Gravitational Lensing Experiment (OGLE III) photometric data. This is a first of its kind map of metallicity up to a radius of 4 - 5 degrees, derived using photometric data and calibrated using spectroscopic data of Red Giant Branch (RGB) stars. We identify the RGB in the V, (V$-$I) colour magnitude diagrams of small subregions of varying sizes in both data sets. We use the slope of the RGB as an indicator of the average metallicity of a subregion, and calibrate the RGB slope to metallicity using spectroscopic data for field and cluster red giants in selected subregions. The average metallicity of the LMC is found to be [Fe/H] = $-$0.37 dex ($\sigma$[Fe/H] = 0.12) from MCPS data, and [Fe/H] = $-$0.39 dex ($\sigma$[Fe/H] = 0.10) from OGLE III data. The bar is found be the most metal-rich region of the LMC. Both the data sets suggest a shallow radial metallicity gradient up to a radius of 4 kpc ($-$0.049$\pm$0.002 dex kpc$^{-1}$ to $-$0.066$\pm$0.006 dex kpc$^{-1}$). Subregions in which the mean metallicity differs from the surrounding areas do not appear to correlate with previously known features; spectroscopic studies are required in order to assess their physical significance.
There is accumulating evidence for the presence of complex molecules, including carbon-bearing and organic molecules, in the interstellar medium. Much of this evidence comes to us from studies of chemical composition, photo- and mass-spectroscopy in cometary, meteoritic and asteroid samples, indicating a need to better understand the surface chemistry of astrophysical objects. There is also considerable interest in the origins of life-forming and life-sustaining molecules on Earth. Here, we perform reactive molecular dynamics simulations to probe the formation of carbon-rich molecules and clusters on carbonaceous surfaces resembling dust grains and meteoroids. Our results show that large chains form on graphitic surfaces at low temperatures (100K - 500K) and smaller fullerene-like molecules form at higher temperatures (2000K - 3000K). The formation is faster on the surface than in the gas at low temperatures but slower at high temperatures as surface interactions prevent small clusters from coagulation. We find that for efficient formation of molecular complexity, mobility about the surface is important and helps to build larger carbon chains on the surface than in the gas phase at low temperatures. Finally, we show that the temperature of the surface strongly determines what kind of structures forms and that low turbulent environments are needed for efficient formation.
In prospect of its application in cryogenic rare-event searches, we have investigated the low-temperature scintillation properties of CaWO4 crystals down to 3.4 K under {\alpha} and {\gamma} excitation. Concerning the scintillation decay times, we observe a long component in the ms range which significantly contributes to the light yield below 40K. For the first time we have measured the temperature dependence of the {\alpha}/{\gamma}- ratio of the light yield. This parameter which can be used to discriminate {\alpha} and {\gamma} events in scintillating bolometers is found to be about 8-15% smaller at low temperatures compared to room temperature.
The Interferometric Monitoring of Gamma-ray Bright Active galactic nuclei (iMOGABA) program provides not only simultaneous multifrequency observations of bright gamma-ray detected active galactic nuclei (AGN), but also covers the highest Very Large Baseline Interferometry (VLBI) frequencies ever being systematically monitored, up to 129 GHz. However, observation and imaging of weak sources at the highest observed frequencies is very challenging. In the second paper in this series, we evaluate the viability of the frequency phase transfer technique to iMOGABA in order to obtain larger coherence time at the higher frequencies of this program (86 and 129 GHz) and image additional sources that were not detected using standard techniques. We find that this method is applicable to the iMOGABA program even under non-optimal weather conditions.
Tidally-excited inertial waves in stellar convective regions are a key mechanism for tidal dissipation in stars and therefore the evolution of close-in binary or planetary systems. As a first step, we explore here the impact of latitudinal differential rotation on the properties of free inertial modes and identify the different families of modes. We show that they differ from the case of solid-body rotation. Using an analytical approach as well as numerical calculations, we conclude that critical layers (where the Doppler-shifted frequency vanishes) could play a very important role for tidal dissipation.
The thermal desorption of ammonia (NH$_3$) from single crystal forsterite (010) has been investigated using temperature-programmed desorption. The effect of defects on the desorption process has been probed by the use of a rough cut forsterite surface prepared from the cleaved forsterite sample. Several approaches have been used to extract the desorption energy and pre-exponential factor describing the desorption kinetics. In the sub-monolayer coverage regime, the NH$_3$ desorption shows a broad distribution of desorption energies, indicating the presence of different adsorption sites, which results in an apparent coverage-dependent desorption energy. This distribution is sensitive to the surface roughness with the cut forsterite surface displaying a significantly broader distribution of desorption energies compared to the cleaved forsterite surface. The cut forsterite surface exhibits sites with desorption energies up to 62.5 kJ mol$^{-1} $ in comparison to a desorption energy of up to 58.0 kJ mol$^{-1} $ for the cleaved surface. Multilayer desorption is independent of the nature of the forsterite surface used, with a desorption energy of ($25.8\pm0.9$) kJ mol$^{-1} $. On astrophysically relevant heating time-scales, the presence of a coverage dependent desorption energy distribution results in a lengthening of the NH$_3$ desorption time-scale by $5.9\times 10^4$ yr compared to that expected for a single desorption energy. In addition, the presence of a larger number of high-energy adsorption sites on the rougher cut forsterite surface leads to a further lengthening of ca. 7000 yr.
The HERMES spectrograph installed on the 1.2-m Mercator telescope has been used to monitor the radial velocity of 13 low-metallicity carbon stars, among which 7 Carbon-Enhanced Metal-Poor (CEMP) stars and 6 CH stars. All stars but one show clear evidence for binarity. New orbits are obtained for 8 systems. The sample covers an extended range in orbital periods, extending from 3.4 d (for the dwarf carbon star HE 0024-2523) to about 54 yr (for the CH star HD 26, the longest known among barium, CH and extrinsic S stars). Three systems exhibit low-amplitude velocity variations with periods close to 1 yr superimposed on a long-term trend. In the absence of an accurate photometric monitoring of these systems, it is not clear yet whether these variations are the signature of a very low-mass companion, or of regular envelope pulsations. The period - eccentricity (P - e) diagram for the 40 low-metallicity carbon stars with orbits now available shows no difference between CH and CEMP-s stars (the latter corresponding to those CEMP stars enriched in s-process elements, as are CH stars). We suggest that they must be considered as one and the same family and that their different names only stem from historical reasons. Indeed, these two families have as well very similar mass-function distributions, corresponding to companions with masses in the range 0.5 - 0.7 Msun, indicative of white-dwarf companions, adopting 0.8 - 0.9 Msun for the primary component. This result confirms that CH and CEMP-s stars obey the same mass-transfer scenario as their higher-metallicity analogs, the barium stars. The P - e diagrams of barium, CH and CEMP-s stars are indeed very similar. They reveal two different groups of systems: one with short orbital periods (P < 1000 d) and mostly circular or almost circular orbits, and another with longer-period and eccentric (e > 0.1) orbits.
Using astrometric observations spanning more than a century and including a large set of Cassini data, we determine Saturn's tidal parameters through their current effects on the orbits of the eight main and four coorbital moons. We have used the latter to make the first determination of Saturn's Love number, $k_2=0.390 \pm 0.024$, a value larger than the commonly used theoretical value of 0.341 (Gavrilov & Zharkov, 1977), but compatible with more recent models (Helled & Guillot, 2013) for which $k_2$ ranges from 0.355 to 0.382. Depending on the assumed spin for Saturn's interior, the new constraint can lead to a reduction of up to 80% in the number of potential models, offering great opportunities to probe the planet's interior. In addition, significant tidal dissipation within Saturn is confirmed (Lainey et al., 2012) corresponding to a high present-day tidal ratio $k_2/Q=(1.59 \pm 0.74) \times 10^{-4}$ and implying fast orbital expansions of the moons. This high dissipation, with no obvious variations for tidal frequencies corresponding to those of Enceladus and Dione, may be explained by viscous friction in a solid core, implying a core viscosity typically ranging between $10^{14}$ and $10^{16}$ Pa.s (Remus et al., 2012). However, a dissipation increase by one order of magnitude at Rhea's frequency could suggest the existence of an additional, frequency-dependent, dissipation process, possibly from turbulent friction acting on tidal waves in the fluid envelope of Saturn (Ogilvie & Li, 2004). Alternatively, a few of Saturn's moons might themselves experience large tidal dissipation.
The James Webb Space Telescope (JWST), as the largest space-based astronomical observatory with near- and mid-infrared instrumentation, will elucidate many mysterious aspects of comets. We summarize four cometary science themes especially suited for this telescope and its instrumentation: the drivers of cometary activity, comet nucleus heterogeneity, water ice in comae and on surfaces, and activity in faint comets and main-belt asteroids. With JWST, we can expect the most distant detections of gas, especially CO2, in what we now consider to be only moderately bright comets. For nearby comets, coma dust properties can be studied with their driving gases, measured simultaneously with the same instrument or contemporaneously with another. Studies of water ice and gas in the distant Solar System will help us test our understanding of cometary interiors and coma evolution. The question of cometary activity in main-belt comets will be further explored with the possibility of a direct detection of coma gas. We explore the technical approaches to these science cases and provide simple tools for estimating comet dust and gas brightness. Finally, we consider the effects of the observatory's non-sidereal tracking limits, and provide a list of potential comet targets during the first 5 years of the mission.
While searches for young stellar objects (YSOs) with the Spitzer Space Telescope focused on known molecular clouds, photometry from the Wide-field Infrared Survey Explorer (WISE) can be used to extend the search to the entire sky. As a precursor to more expansive searches, we present results for a 100 square degree region centered on the Canis Major clouds.
The enrichment in s-process elements of barium stars is known to be due to pollution by mass transfer from a companion formerly on the thermally-pulsing asymptotic giant branch (AGB), now a carbon-oxygen white-dwarf (WD). This paper investigates the relationship between the s-process enrichment in the barium star and the mass of its WD companion. It is expected that helium WDs, which have masses smaller than about 0.5 Msun and never reached the AGB phase, should not pollute with s-process elements their giant companion, which should thus never turn into a barium star. Spectra with a resolution of R ~ 86000 were obtained with the HERMES spectrograph on the 1.2-m Mercator telescope for a sample of 11 binary systems involving WD companions of various masses. We use standard 1D LTE MARCS model atmospheres coupled with the Turbospectrum radiative-transfer code to derive the atmospheric parameters using equivalent widths of FeI and Fe II lines. The abundances of s-process elements for the entire sample of 11 binary stars were derived homogeneously. The sample encompasses all levels of overabundances: from solar [s/Fe]=0 to 1.5 dex in the 2 binary systems with S-star primaries (for which dedicated MARCS model atmospheres were used). The primary components of binary systems with a WD more massive than 0.5 Msun are enriched in s-process elements. We also found a trend of increasing [s/Fe] with [C/Fe] or [(C+N)/Fe]. Our results conform to the expectation that binary systems with WD companions less massive than 0.5 Msun do not host barium stars.
IceCube collaboration has seen an unexpected population of high energy neutrinos compatible with an astrophysical origin. We consider two categories of events that can help to diagnose cosmic neutrinos: double pulse, that may allow us to clearly discriminate the cosmic component of tau neutrinos; cascades with deposited energy above 2 PeV, including events produced by electron antineutrinos at the Glashow resonance, that can be used to investigate the neutrino production mechanisms. We show that one half of the double pulse signal is due to the neutrinos spectral region already probed by IceCube. By normalizing to HESE data, we find that 10 more years are required to obtain 90% probability to observe a double pulse. The cascades above 2 PeV provide us a sensitive probe of the high energy tail of the neutrino spectrum and are potentially observable, but even in this case, the dependence on type of the source is mild. In fact we find that pp or p{\gamma} mechanisms give a difference in the number of cascades above 2 PeV of about 25 % that can be discriminated at 2{\sigma} in about 50 years of data taking.
The statistical properties of the temperature anisotropies and polarization of the of cosmic microwave background (CMB) radiation offer a powerful probe of the physics of the early universe. In recent works a statistical procedure based upon the calculation of the kurtosis and skewness of the data in patches of CMB sky-sphere has been proposed and used to investigate the large-angle deviation from Gaussianity in WMAP maps. Here we briefly address the question as to how this analysis of Gaussianity is modified if the foreground-cleaned Planck maps are considered. We show that although the foreground-cleaned Planck maps present significant deviation from Gaussianity of different degrees when a less severe mask is used, they become consistent with Gaussianity, as detected by our indicators, when masked with the union mask U73.
High angular resolution images of IRC+10216 are presented in several near infrared wavelengths spanning more than 8 years. These maps have been reconstructed from interferometric observations obtained at both Keck and the VLT, and also from stellar occultations by the rings of Saturn observed with the Cassini spacecraft. The dynamic inner regions of the circumstellar environment are monitored over eight epochs ranging between January 2000 and July 2008. The system is shown to experience substantial evolution within this period including the fading of many previously reported persistent features, some of which had been identified as the stellar photosphere. These changes are discussed in context of existing models for the nature of the underlying star and the circumstellar environment. With access to these new images, we are able to report that none of the previously identified bright spots in fact contain the star, which is buried in its own dust and not directly visible in the near infrared.
The finding of a double red clump in the luminosity function of the Milky Way bulge has been interpreted as evidence for an X-shaped structure. Recently, an alternative explanation has been suggested, where the double red clump is an effect of multiple stellar populations in a classical spheroid. In this letter we provide an observational assessment of this scenario and show that it is not consistent with the behaviour of the red clump across different lines of sight, particularly at high distances from the Galactic plane. Instead, we confirm that the shape of the red clump magnitude distribution closely follows the distance distribution expected for an X-shaped bulge at critical Galactic latitudes. We also emphasize some key observational properties of the bulge red clump that should not be neglected in the search for alternative scenarios.
Using numerical simulations we present the accretion of supersonic winds onto a rotating black hole in three dimensions. We study five representative directions of the wind with respect to the axis of rotation of the black hole and focus on the evolution and stability of the high density shock cone that is formed during the process. We explore both, the regime in which the shock cone is expected to be stable in order to confirm previous results obtained with two dimensional simulations, and the regime in which the shock cone is expected to show a flip-flop type of instability. The methods used to attempt triggering the instability were first the accumulation of numerical errors and second the explicit application of a perturbation on the velocity field after the shock-cone was formed. The result is negative, that is, we did not find the flip-flop instability within the parameter space we explored, which includes cases that are expected to be unstable.
We present the study of the structure and dynamical status of the galaxy system A523. Our analysis is based on new spectroscopic data for 132 galaxies (TNG), new photometric data (INT), X-ray data (Chandra archive), and radio data (VLA archive). We present the first measures of velocity dispersion of the galaxy population, 949 km/s, and global X-ray temperature of the hot ICM, 5.3 keV. We infer that A523 is a massive system, M200 about 7-9 10E14 Msun. Our analysis of optical data confirms the presence of two subclusters, 0.75 Mpc apart, tracing the SSW-NNE direction, finds that they are (little) separated in velocity, and identifies the two dominant galaxies (BCG1 and BCG2). We show that the X-ray surface brightness is strongly elongated towards the NNE direction, and its peak is clearly offsetted from both the BCGs, and quantify the presence of substructure. We confirm the presence of a 1.3 Mpc large central radio source, its main ESE-WNW elongation perpendicular to the optical/X-ray elongation, and the previous halo classification. We determine a large radio/X-ray peaks offset and detect evidence of polarization, being this detected in only very few radio halos. The radio/X-ray offset and polarization might be the result of having most magnetic field energy on large spatial scales, as shown by our ad hoc simulations. Most properties are consistent with scaling relations followed by other clusters hosting radio halos, but A523 is shown to be peculiar in the Pradio-Lx plane, having a higher radio power or a lower X-ray luminosity than expected. According to main optical and X-ray features, A523 can be described as a binary head--on merger after the primary collision in the SSW-NNE direction. However, both optical and radio data show some evidence in favor of a more complex cluster structure with A523 forming at the cross of two filaments, along SSW-NNE and ESE-WNW directions.
We present results obtained by applying our BAyesian HierArchical Modeling
for the Analysis of Supernova cosmology (BAHAMAS) software package to the 740
spectroscopically confirmed supernovae type Ia (SNIa) from the "Joint
Light-curve Analysis" (JLA) dataset. We simultaneously determine cosmological
parameters and standardization parameters, including host galaxy mass
corrections, residual scatter and object-by-object intrinsic magnitudes.
Combining JLA and Planck Cosmic Microwave Background data, we find significant
discrepancies in cosmological parameter constraints with respect to the
standard analysis: we find Omega_M = 0.399+/-0.027, 2.8\sigma\ higher than
previously reported and w = -0.910+/-0.045, 1.6\sigma\ higher than the standard
analysis. We determine the residual scatter to be sigma_res = 0.104+/-0.005.
We confirm (at the 95% probability level) the existence of two
sub-populations segregated by host galaxy mass, separated at
log_{10}(M/M_solar) = 10, differing in mean intrinsic magnitude by
0.055+/-0.022 mag, lower than previously reported. Cosmological parameter
constraints are however unaffected by inclusion of host galaxy mass
corrections. We find ~4\sigma\ evidence for a sharp drop in the value of the
color correction parameter, beta(z), at a redshift z_trans = 0.662+/-0.055. We
rule out some possible explanations for this behaviour, which remains
unexplained.
The most general scalar-tensor theories of gravity predict a weakening of the gravitational force inside astrophysical bodies. There is a minimum mass for hydrogen burning in stars that is set by the interplay of plasma physics and the theory of gravity. We calculate this for alternative theories of gravity, and find that it is always significantly larger than the general relativity prediction. The observation of several low mass Red Dwarf stars therefore rules out a large class of scalar-tensor gravity theories, and places strong constraints on the cosmological parameters appearing in the effective field theory of dark energy.
A relaxed primordial power spectrum (PPS) of scalar perturbations arising from inflation can impact the dark radiation constraints obtained from Cosmic Microwave Background and other cosmological measurements. If inflation produces a non-standard PPS for the initial fluctuations, a fully thermalized light sterile neutrino can be favoured by CMB observations, instead of being strongly disfavoured. In the case of a thermal axion, the constraints on the axion mass are relaxed when the PPS is different from the standard power law.
Aspects of solar flare dynamics, such as chromospheric evaporation and flare light-curves, have long been studied using one-dimensional models of plasma dynamics inside a static flare loop, subjected to some energy input. While extremely successful at explaining the observed characteristics of flares, all such models so far have specified energy input ad hoc, rather than deriving it self-consistently. There is broad consensus that flares are powered by magnetic energy released through reconnection. Recent work has generalized Petschek's basic reconnection scenario, topological change followed by field line retraction and shock heating, to permit its inclusion into a one-dimensional flare loop model. Here we compare the gas dynamics driven by retraction and shocking to those from more conventional static loop models energized by ad hoc source terms. We find significant differences during the first minute, when retraction leads to larger kinetic energies and produces higher densities at the loop top, while ad hoc heating tends to rarify the loop top. The loop-top density concentration is related to the slow magnetosonic shock, characteristic of Petschek's model, but persists beyond the retraction phase occurring in the outflow jet. This offers an explanation for observed loop-top sources of X-ray and EUV emission, with advantages over that provided by ad hoc heating scenarios. The cooling phases of the two models are, however, notably similar to one another, suggesting observations at that stage will yield little information on the nature of energy input.
Hot subdwarf stars (sdO/Bs) are the stripped cores of red giants located at the bluest extension of the horizontal branch. They constitute the dominant population of UV-bright stars in old stellar environments and are most likely formed by binary interactions. We perform the first systematic, spectroscopic analysis of a sample of those stars in the Galactic halo based on data from SDSS. In the course of this project we discovered 177 close binary candidates. A significant fraction of the sdB binaries turned out to have close substellar companions, which shows that brown dwarfs and planets can significantly influence late stellar evolution. Close hot subdwarf binaries with massive white dwarf companions on the other hand are good candidates for the progenitors of type Ia supernovae. We discovered a hypervelocity star, which not only turned out to be the fastest unbound star known in our Galaxy, but also the surviving companion of such a supernova explosion.
We examine metallicities, ages and orbital properties of halo stars in a Milky-Way like disk galaxy formed in the cosmological hydrodynamical MaGICC simulations. Halo stars were either accreted from satellites or they formed in situ in the disk or bulge of the galaxy and were then kicked up into the halo ("in situ/ kicked-up" stars). Regardless of where they formed both types show surprisingly similar orbital properties: the majority of stars of both types are on short-axis tubes with the same sense of rotation as the disk -- implying that a large fraction of satellites are accreted onto the halo with the same sense of angular momentum as the disk.
Planets less massive than Saturn tend to rapidly migrate inward in protoplanetary disks. This is the so-called type I migration. Simulations attempting to reproduce the observed properties of exoplanets show that type I migration needs to be significantly reduced over a wide region of the disk for a long time. However, the mechanism capable of suppressing type I migration over a wide region has remained elusive. The recently found turbulence-driven disk winds offer new possibilities. We investigate the effects of disk winds on the disk profile and type I migration for a range of parameters that describe the strength of disk winds. We also examine the in situ formation of close-in super-Earths in disks that evolve through disk winds. The disk profile, which is regulated by viscous diffusion and disk winds, was derived by solving the diffusion equation. We carried out a number of simulations and plot here migration maps that indicate the type I migration rate. We also performed N-body simulations of the formation of close-in super-Earths from a population of planetesimals and planetary embryos. We define a key parameter, Kw, which determines the ratio of strengths between the viscous diffusion and disk winds. For a wide range of Kw, the type I migration rate is presented in migration maps. These maps show that type I migration is suppressed over the whole close-in region when the effects of disk winds are relatively strong (Kw < 100). From the results of N-body simulations, we see that type I migration is significantly slowed down assuming Kw = 40. We also show that the results of N-body simulations match statistical orbital distributions of close-in super-Earths.
We characterize the physical properties of the cool T ~ 10^4 K circumgalactic medium surrounding z ~ 2-3 quasar host galaxies, which are predicted to evolve into present day massive ellipticals. Using a statistical sample of 14 quasar pairs with projected separation < 300 kpc and high dispersion, high S/N spectra, we find extreme kinematics with low metal ion lines typically spanning 500 km s^-1, exceeding any previously studied galactic population. The CGM is significantly enriched, even beyond the virial radius, with a median metallicity [M/H] = -0.6. The alpha/Fe abundance ratio is enhanced, suggesting that halo gas is primarily enriched by Type II supernovae. The total mass of the cool CGM is estimated to be 1.9*10^11 M_sun (R_\perp/160 kpc)^2, accounting for 1/3 of the galaxy halo baryonic budget. The ionization state of CGM gas increases with projected distance from the foreground quasars, contrary to expectation if the quasar dominates the ionizing radiation flux. However, we also found peculiarities not exhibited in the CGM of other galaxy populations. In one absorption system, we may be detecting unresolved fluorescent Ly-alpha emission, and another system shows strong NV lines. Taken together these anomalies suggest that transverse sightlines are at least in some cases possibly illuminated. We also discovered a peculiar case where detection of the CII fine structure line implies an electron density > 100 cm^-3 and subparsec scale gas clumps.
In this paper, we address a possible impact of radiative corrections from a heavy scalar field $\chi$ on the curvature perturbation $\zeta$. Integrating out $\chi$, we derive the effective action for $\zeta$, which includes the loop corrections of the heavy field $\chi$. When the mass of $\chi$ is much larger than the Hubble scale $H$, the loop corrections of $\chi$ only yield a local contribution in the effective action and hence the effective action simply gives an action for $\zeta$ in a single field model, where, as is widely known, $\zeta$ is conserved in time after the Hubble crossing time. Meanwhile, when the mass of $\chi$ is comparable to $H$, the loop corrections of $\chi$ can give a non-local contribution to the effective action. Because of the non-local contribution from $\chi$, in general, $\zeta$ may not be conserved, even if the classical background trajectory is determined only by the evolution of the inflaton. In this paper, we derive the condition that $\zeta$ is conserved in time in the presence of the radiative corrections from $\chi$. Namely, we show that when the scaling symmetry, which is a part of the diffeomorphism invariance, is preserved at the quantum level, the loop corrections of the massive field $\chi$ do not disturb the constant evolution of $\zeta$ at super Hubble scales. In this discussion, we show the Ward-Takahashi identity for the scaling symmetry, which yields a consistency relation for the correlation functions of the massive field $\chi$.
We calculate the production rate of singlet fermions from the decay of neutral or charged scalar fields in a hot plasma. We find that there are considerable thermal corrections when the temperature of the plasma exceeds the mass of the decaying scalar. We give analytic expressions for the temperature corrected production rates in the regime where the decay products are relativistic. We also study the regime of non-relativistic decay products numerically. If the scalar is a singlet, then the main correction arises from the fact that the effective mass of the decaying particle is temperature dependent. If the scalar has gauge interactions, then at least one of the decay products must also carry gauge quantum numbers and will be in thermal equilibrium in the early universe. This gives rise to Pauli blocking and a modification of the dispersion relation of these charged particles in the plasma, which affects the kinematics of the decay. Moreover, at high temperature, inelastic scatterings can contribute to the singlet fermion production rate. As a result, the production rate of singlet fermions has a non-trivial temperature dependence. This can affect the abundance and momentum distribution of sterile neutrino Dark Matter that is produced in scalar decays. Our results can be used to improve predictions for the free streaming of the Dark Matter particles, which is crucial to test the compatibility of such scenarios with cosmic structure formation.
Rotating black holes can support quasi-stationary (unstable) bound-state resonances of massive scalar fields in their exterior regions. These spatially regular scalar configurations are characterized by instability timescales which are much longer than the timescale $M$ set by the geometric size (mass) of the central black hole. It is well-known that, in the small-mass limit $\alpha\equiv M\mu\ll1$ (here $\mu$ is the mass of the scalar field), these quasi-stationary scalar resonances are characterized by the familiar hydrogenic oscillation spectrum: $\omega_{\text{R}}/\mu=1-\alpha^2/2{\bar n}^2_0$, where the integer $\bar n_0(l,n;\alpha\to0)=l+n+1$ is the principal quantum number of the bound-state resonance (here the integers $l=1,2,3,...$ and $n=0,1,2,...$ are the spheroidal harmonic index and the resonance parameter of the field mode, respectively). As it depends only on the principal resonance parameter $\bar n_0$, this small-mass ($\alpha\ll1$) hydrogenic spectrum is obviously degenerate. In this paper we go beyond the small-mass approximation and analyze the quasi-stationary bound-state resonances of massive scalar fields in rapidly-spinning Kerr black-hole spacetimes in the regime $\alpha=O(1)$. In particular, we derive the non-hydrogenic (and, in general, non-degenerate) resonance oscillation spectrum ${{\omega_{\text{R}}}/{\mu}}=\sqrt{1-(\alpha/{\bar n})^2}$, where $\bar n(l,n;\alpha)=\sqrt{(l+1/2)^2-2m\alpha+2\alpha^2}+1/2+n$ is the generalized principal quantum number of the quasi-stationary resonances. This analytically derived formula for the characteristic oscillation frequencies of the composed black-hole-massive-scalar-field system is shown to agree with direct numerical computations of the quasi-stationary bound-state resonances.
When did the universe thermalize? In this talk I review the status of this issue and its importance in establishing the expected properties of dark matter, the growth of large-scale structure, and the viability of inflation models when confronted with CMB observations. I also present a novel approach to tackling the theoretical challenges surrounding inflationary (p)reheating, which seeks to extend past work on the Effective Field Theory of Inflation to the time of reheating.
We describe how monomial chaotic inflation becomes compatible with the latest CMB data thanks to radiative corrections producing a plateau. The interactions of the inflation with other fields, required for reheating, can flatten the potential and moderate the production of primordial gravitational waves, keeping these below the current upper bound. We show that the appearance of a plateau requires that the inflaton couples to fermions and to another scalar or a gauge group. We give concrete examples of minimal particle physics models leading to plateaus for quadratic and quartic chaotic inflation. We also provide a three-parameter model-independent description of radiatively corrected inflation that is amenable to CMB analyses.
The late universe's matter distribution obeys the Copernican principle at only the coarsest of scales. The relative importance of such inhomogeneity is still not well understood. Because of the Einstein field equations' non-linear nature, some argue a non-perturbative approach is necessary to correctly model inhomogeneities and may even obviate any need for dark energy. We shall discuss an approach based on Regge calculus, a discrete approximation to general relativity: we shall discuss the Collins--Williams formulation of Regge calculus and its application to two toy universes. The first is a universe for which the continuum solution is well-established, the $\Lambda$-FLRW universe. The second is an inhomogeneous universe, the `lattice universe' wherein matter consists solely of a lattice of point masses with pure vacuum in between, a distribution more similar to that of the actual universe compared to FLRW universes. We shall discuss both regular lattices and one where one mass gets perturbed.
We obtain the critical magnetic field required for complete destruction of $S$-wave pairing in neutron matter, thereby setting limits on the pairing and superfluidity of neutrons in the crust and outer core of magnetars. We find that for fields $B \ge 10^{17}$ G the neutron fluid is non-superfluid, a result with profound consequences for the thermal, rotational, and oscillatory behavior of magnetars. Since the dineutron is not bound in vacuum, cold dilute neutron matter cannot exhibit a proper BCS-BEC crossover. Nevertheless, owing to the strongly resonant behavior of the $nn$ interaction at low densities, neutron matter shows a precursor of the BEC state, as manifested in Cooper-pair correlation lengths being comparable to the interparticle distance. We make a systematic quantitative study of this type of BCS-BEC crossover in the presence of neutron fluid spin-polarization induced by an ultra-strong magnetic field. We evaluate the Cooper pair wave-function, quasiparticle occupation numbers, and quasiparticle spectra for densities and temperatures spanning the BCS-BEC crossover region. The phase diagram of spin-polarized neutron matter is constructed and explored at different polarizations.
Neutrinos are copiously emitted from black hole accretion disks playing a fundamental role in their evolution, as well as in the production of gamma ray bursts and r-process nucleosynthesis. The black hole generates a strong gravitational field able to change the properties of the emerging neutrinos. We study the influence of the black hole spin on the structure of the neutrino surfaces, neutrino luminosities, average neutrino energies, and event counts at SuperK. We consider several disk models and provide estimates that cover different black hole efficiency scenarios. We discuss the influence of the detector's inclination with respect to the axis of the torus on neutrino properties. We find that tori around spinning black holes have larger luminosities, energies and rates compared to tori around static black holes, and that the inclination of the observer causes a reduction in the luminosities and detection rates but an increase in the average energies.
Links to: arXiv, form interface, find, astro-ph, recent, 1510, contact, help (Access key information)
Massive stars shape their surroundings with mass loss from winds during their lifetimes. Fast ejecta from supernovae, from these massive stars, shocks this circumstellar medium. Emission generated by this interaction provides a window into the final stages of stellar evolution, by probing the history of mass loss from the progenitor. Here we use Chandra and Swift x-ray observations of the type II-P/L SN 2013ej to probe the history of mass loss from its progenitor. We model the observed x-rays as emission from both heated circumstellar matter and supernova ejecta. The circumstellar density profile probed by the supernova shock reveals a history of steady mass loss during the final 400 years. The inferred mass loss rate of $2 \times 10^{-6} {\rm \; M_\odot \; yr^{-1}}$ points back to a 13 $M_\odot$ progenitor. Soon after the explosion we find significant absorption of reverse shock emission by a cooling shell. The column depth of this shell observed in absorption provides an independent and consistent measurement of the circumstellar density seen in emission. We also determine the efficiency of cosmic ray acceleration from x-rays produced by Inverse Compton scattering of optical photons by relativistic electrons. Only about 1 percent of the thermal energy is used to accelerate electrons. Our x-ray observations and modeling provides stringent tests for models of massive stellar evolution and micro-physics of shocks.
We present a study of cosmological implications of generic dark matter decays. We consider two-body and many-body decaying scenarios. In the two-body case the massive particle has a possibly relativistic kick velocity and thus possesses a dynamical equation of state. This has implications to the expansion history of the universe. We use recent observational data from the cosmic microwave background, baryon acoustic oscillations and supernovae Type Ia to obtain constraints on the lifetime of the dark matter particle. We find that for an energy splitting where more than 40% of the dark matter particle energy is transferred to massless, relativistic particles in the two-body case, or more than 50% in the many-body case, lifetimes less than the age of the universe are excluded at more than 95% confidence. When the energy splitting falls to 10% the lifetime is constrained to be more than roughly half the age.
The stellar initial mass function (IMF) is a fundamental astrophysical quantity that impacts a wide range of astrophysical problems from heavy element distribution to galactic evolution to planetary system formation. However, the origin and universality of the IMF are hotly debated both observationally and theoretically. I review recent observations of the IMF across a variety of environments. These suggest the IMF is surprisingly invariant between star-forming regions, star clusters, and spiral galaxies but that it may also vary under extreme conditions, including within the Galactic center and early type galaxies.
A large fraction of the dwarf satellite galaxies orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and seemingly kinematically coherent planar structure. Such a structure is not easily found in simulations based on the Cold Dark Matter model. Using 21 high resolution cosmological simulations based on this model we analyze in detail the kinematical structure of planes of satellites resembling the one observed around Andromeda when co-rotation is characterized by the line-of-sight velocity. At the same time, when co-rotation is inferred by the angular momenta of the satellites, the planes are in excellent agreement with the plane around the Milky Way. Furthermore, we find such planes to be common in our simulations. Investigation of the kinematics of the satellites in the plane reveals that the number of co-rotating satellites varies by 2 to 5 out of ~12 depending on the viewing angle. These variations are consistent with that obtained from a sample with random velocities. Using instead the clustering of angular momentum vectors of the satellites in the plane results in a better measure of kinematic coherence. Thus we conclude that the line-of-sight velocity as a proxy for the kinematical coherence of the plane is not a robust measure. Detailed analysis of the kinematics of our planes shows that the planes consist of ~30% chance aligned satellites. Tracking the satellites in the plane back in time reveals that the plane is a transient feature and not kinematically coherent as would appear at first sight.
We present a potentially efficient dynamical formation scenario for Low Mass X-ray Binaries (LMXBs) in the field, focusing on black-hole (BH) LMXBs. In this formation channel LMXBs are formed from wide binaries $(>1000$ AU) with a BH component and a stellar companion. The wide binary is perturbed by fly-by's of field stars and its orbit random-walks and changes over time. This diffusion process can drive the binary into a sufficiently eccentric orbit such that the binary components tidally interact at peri-center and the binary evolves to become a short period binary, which eventually evolves into an LMXB. The formation rate of LMXBs through this channel mostly depends on the number of such BH wide binaries progenitors, which in turn depends on the velocity kicks imparted to BHs (or NSs) at birth. We consider several models for the formation and survival of such wide binaries, and calculate the LMXB formation rates for each model. We find that models where BHs form through direct collapse with no/little natal kicks can give rise to high formation rates comparable with those inferred from observations. This formation scenario had several observational signatures: (1) the number density of LMXBs generally follows the background stellar density and (2) the mass function of the BH stellar companion should be comparable to the mass function of the background stellar population, likely peaking at $0.4-0.6$ M$_{\odot}$. The latter aspect, in particular, is unique to this model compared with previously suggested LMXB formation models following common envelope binary stellar evolution. We note that NS LMXBs can similarly form from wide binaries, but their formation rate through this channel is likely significantly smaller due to their much higher natal kicks.
How starburst clusters form out of molecular clouds is still an open question. In this article, I highlight some of the key constraints in this regard, that one can get from the dynamical evolutionary properties of dense stellar systems. I particularly focus on secular expansion of massive star clusters and hierarchical merging of sub-clusters, and discuss their implications vis-a-vis the observed properties of young massive clusters. The analysis suggests that residual gas expulsion is necessary for shaping these clusters as we see them today, irrespective of their monolithic or hierarchical mode of formation.
We present medium-resolution optical spectroscopy with the SOAR telescope of the O star secondary of the high-mass gamma-ray binary 1FGL J1018.6-5856 to help determine whether the primary is a neutron star or black hole. We find that the secondary has a low radial velocity semi-amplitude of 11-12 km/s, with consistent values obtained for H and He absorption lines. This low value strongly favors a neutron star primary: while a black hole cannot be excluded if the system is close to face on, such inclinations are disallowed by the observed rotation of the secondary. We also find the high-energy (X-ray and gamma-ray) flux maxima occur when the star is behind the compact object along our line of sight, inconsistent with a simple model of anisotropic inverse Compton scattering for the gamma-ray photons.
Future large-scale galaxy surveys have the potential to become leading probes for cosmology provided the influence of baryons on the total mass distribution is understood well enough. As hydrodynamical simulations strongly depend on details in the feedback implementations, no unique and robust predictions for baryonic effects currently exist. In this paper we propose a baryonic correction model that modifies the density field of dark-matter-only $N$-body simulations to mimic the effects of baryons from any underlying adopted feedback recipe. The model assumes haloes to consist of 4 components: 1- hot gas in hydrostatical equilibrium, 2- ejected gas from feedback processes, 3- central galaxy stars, and 4- adiabatically relaxed dark matter, which all modify the initial dark-matter-only density profiles. This altered mass profiles allow to define a displacement field for particles in $N$-body simulations and to modify the total density field accordingly. The main advantage of the baryonic correction model is to connect the total matter density field to the observable distribution of gas and stars in haloes, making it possible to parametrise baryonic effects on the matter power spectrum. We show that the most crucial quantities are the mass fraction of ejected gas and its corresponding ejection radius. The former controls how strongly baryons suppress the power spectrum, while the latter provides a measure of the scale where baryonic effects become important. A comparison with X-ray and Sunyaev-Zel'dovich cluster observations suggests that baryons suppress wave modes above $k\sim0.5$ h/Mpc with a maximum suppression of 10-25 percent around $k\sim 2$ h/Mpc. More detailed observations of the gas in the outskirts of groups and clusters are required to decrease the large uncertainties of these numbers.
We investigate optical, infrared, and radio active galactic nucleus (AGN) signs in the merger remnant Arp 187, which hosts luminous jets launched in the order of $10^5$ yr ago but whose present-day AGN activity is still unknown. We find AGN signs from the optical BPT diagram and infrared [OIV]25.89 $\mu$m line, originating from the narrow line regions of AGN. On the other hand, Spitzer/IRS show the host galaxy dominated spectra, suggesting that the thermal emission from the AGN torus is considerably small or already diminished. Combining the black hole mass, the upper limit of radio luminosity of the core, and the fundamental plane of the black hole enable us to estimate X-ray luminosity, which gives $<10^{40}$ erg s$^{-1}$. Those results suggest that the AGN activity of Arp 187 has already been quenched, but the narrow line region is still alive owing to the time delay of emission from the past AGN activity.
In the low redshift Universe (z<0.3), our view of galaxy evolution is primarily based on fibre optic spectroscopy surveys. Elaborate methods have been developed to address aperture effects when fixed aperture sizes only probe the inner regions for galaxies of ever decreasing redshift or increasing physical size. These aperture corrections rely on assumptions about the physical properties of galaxies. The adequacy of these aperture corrections can be tested with integral-field spectroscopic data. We use integral-field spectra drawn from 1212 galaxies observed as part of the SAMI Galaxy Survey to investigate the validity of two aperture correction methods that attempt to estimate a galaxy's total instantaneous star formation rate. We show that biases arise when assuming that instantaneous star formation is traced by broadband imaging, and when the aperture correction is built only from spectra of the nuclear region of galaxies. These biases may be significant depending on the selection criteria of a survey sample. Understanding the sensitivities of these aperture corrections is essential for correct handling of systematic errors in galaxy evolution studies.
The search for the curl component (B mode) in the cosmic microwave background (CMB) polarization induced by inflationary gravitational waves is described. The canonical single-field slow-roll model of inflation is presented, and we explain the quantum production of primordial density perturbations and gravitational waves. It is shown how these gravitational waves then give rise to polarization in the CMB. We then describe the geometric decomposition of the CMB polarization pattern into a curl-free component (E mode) and curl component (B mode) and show explicitly that gravitational waves induce B modes. We discuss the B modes induced by gravitational lensing and by Galactic foregrounds and show how both are distinguished from those induced by inflationary gravitational waves. Issues involved in the experimental pursuit of these B modes are described, and we summarize some of the strategies being pursued. We close with a brief discussion of some other avenues toward detecting/characterizing the inflationary gravitational-wave background.
In our recent investigation (Lim et al. 2015), we have shown that narrow-band photometry can be combined with low-resolution spectroscopy to effectively search for globular clusters (GCs) with supernovae (SNe) enrichments. Here we apply this technique to the metal-poor bulge GC NGC 6273, and find that the red giant branch stars in this GC are clearly divided into two distinct subpopulations having different calcium abun- dances. The Ca rich subpopulation in this GC is also enhanced in CN and CH, showing a positive correlation between them. This trend is identical to the result we found in M22, suggesting that this might be a ubiquitous nature of GCs more strongly affected by SNe in their chemical evolution. Our results suggest that NGC 6273 was massive enough to retain SNe ejecta which would place this cluster in the growing group of GCs with Galactic building block characteristics, such as {\omega} Centauri and Terzan 5.
A radial-velocity monitoring of the Carbon-Enhanced Metal-Poor (CEMP) star HE 0017+0055 over 8 years with the Nordic Optical Telescope and Mercator telescopes reveals variability with a period of 384 d and amplitude of 540$\pm27$ m s$^{-1}$, superimposed on a nearly linear long-term decline of $\sim$1 m s$^{-1}$ day$^{-1}$. High-resolution HERMES/Mercator and Keck/HIRES spectra have been used to derive elemental abundances using 1-D LTE MARCS models. A metallicity of [Fe/H] $\sim -2.4$ is found, along with s-process overabundances on the order of 2 dex (with the exception of [Y/Fe] $\sim+0.5$), and most notably overabundances of r-process elements like Sm, Eu, Dy, and Er in the range 0.9 - 2.0 dex. With [Ba/Fe] $ > 1.9$ dex and [Eu/Fe] = 2.3 dex, HE 0017+0055 is a CEMP-rs star. It appears to be a giant star below the tip of the red giant branch (RGB). The s-process pollution must therefore originate from mass transfer from a companion formerly on the AGB, now a carbon-oxygen white dwarf (WD). If the 384 d velocity variations are attributed to the WD companion, its orbit must be seen almost face-on, with $i \sim 2.3^\circ$, because the mass function is very small: $f(M_1,M_2) = (6.1\pm1.1)\times10^{-6}$ Msun. Alternatively, the WD orbital motion could be responsible for the long-term velocity variations, with a period of several decades. The 384 d variations should then be attributed either to a low-mass inner companion (perhaps a brown dwarf, depending on the orbital inclination), or to stellar pulsations. The latter possibility is made likely by the fact that similar low-amplitude velocity variations, with periods close to 1 yr, have been reported for other CEMP stars in a companion paper (Jorissen et al., 2015). A definite conclusion about the origin of the 384 d velocity variations should however await the detection of synchronous low-amplitude photometric variations.
We have obtained spectropolarimetric observations of the four images of the gravitationally lensed broad absorption line quasar H1413+117. The polarization of the microlensed image D is significantly different both in the continuum and in the broad lines from the polarization of image A which is essentially unaffected by microlensing. The observations suggest that the continuum is scattered off two regions spatially separated and producing roughly perpendicular polarizations. These results are compatible with a model in which the microlensed polarized continuum comes from a compact region located in the equatorial plane close to the accretion disk and the non-microlensed continuum from an extended region located along the polar axis.
3C318, a radio-loud quasar at z=1.574, is a subgalactic-sized radio source, and a good test-bed for the interplay between black hole and galaxy growth in the high-z Universe. Based on its IRAS, ISO, and SCUBA detections, it has long been considered as one of the most intrinsically luminous (L$_{\mathrm{IR}}$ > 10$^{13}$ L$_{\odot}$) infrared sources in the Universe. Recent far-infrared data from the Herschel Space Observatory reveal that most of the flux associated with 3C318 measured with earlier instruments in fact comes from a bright nearby source. Optical imaging and spectroscopy show that this infrared-bright source is a strongly star-forming pair of interacting galaxies at z=0.35. Adding existing Spitzer and SDSS photometry, we perform a spectral energy distribution analysis of the pair, and find that it has a combined infrared luminosity of L$_{\mathrm{IR}}$ = 1.5 $\times$ 10$^{12}$ L$_{\odot}$, comparable to other intermediate-redshift ultra-luminous infrared galaxies studied with Herschel. Isolating the emission from 3C318's host, we robustly constrain the level of star formation to a value a factor of three lower than that published earlier, which is more in line with the star formation activity found in other Herschel-detected 3CR objects at similar redshift.
We explore star-formation histories (SFHs) of galaxies based on the evolution of the star-formation rate stellar mass relation (SFR-M*). Using data from the FourStar Galaxy Evolution Survey (ZFOURGE) in combination with far-IR imaging from the Spitzer and Herschel observatories we measure the SFR-M* relation at 0.5 < z < 4. Similar to recent works we find that the average infrared SEDs of galaxies are roughly consistent with a single infrared template across a broad range of redshifts and stellar masses, with evidence for only weak deviations. We find that the SFR-M* relation is not consistent with a single power-law of the form SFR ~ M*^a at any redshift; it has a power-law slope of a~1 at low masses, and becomes shallower above a turnover mass (M_0) that ranges from 10^9.5 - 10^10.8 Msol, with evidence that M_0 increases with redshift. We compare our measurements to results from state-of-the-art cosmological simulations, and find general agreement in the slope of the SFR-M* relation albeit with systematic offsets. We use the evolving SFR-M* sequence to generate SFHs, finding that typical SFRs of individual galaxies rise at early times and decline after reaching a peak. This peak occurs earlier for more massive galaxies. We integrate these SFHs to generate mass-growth histories and compare to the implied mass-growth from the evolution of the stellar mass function. We find that these two estimates are in broad qualitative agreement, but that there is room for improvement at a more detailed level. At early times the SFHs suggest mass-growth rates that are as much as 10x higher than inferred from the stellar mass function. However, at later times the SFHs under-predict the inferred evolution, as is expected in the case of additional growth due to mergers.
In a search for periodic variation in the arrival times of pulses from 151 young, energetic pulsars, we have identified seven cases of modulation consistent with one or two harmonics of a single fundamental with time-scale 0.5-1.5 yr. We use simulations to show that these modulations are statistically significant and of high quality (sinusoidal) even when contaminated by the strong stochastic timing noise common to young pulsars. Although planetary companions could induce such modulation, the large implied masses and 2:1 mean motion resonances challenge such an explanation. Instead, the modulation is likely to be intrinsic to the pulsar, arising from quasi-periodic switching between stable magnetospheric states, and we propose that precession of the neutron star may regulate this switching.
Late on July 23, 2012, the STEREO-A spacecraft encountered a fast forward shock driven by a coronal mass ejection launched from the Sun earlier that same day. The estimated travel time of the disturbance ($\sim 20$ hrs), together with the massive magnetic field strengths measured within the ejecta ($> 100$nT), made it one of the most extreme events observed during the space era. In this study, we examine the properties of the shock wave. Because of an instrument malfunction, plasma measurements during the interval surrounding the CME were limited, and our approach has been modified to capitalize on the available measurements and suitable proxies, where possible. We were able to infer the following properties. First, the shock normal was pointing predominantly in the radial direction (${\bf n} = 0.97 {\bf e}_r -0.09 {\bf e}_t -0.23 {\bf e}_n$). Second, the angle between ${\bf n}$ and the upstream magnetic field, $\theta_{Bn}$, was estimated to be $\approx 34^{\circ}$, making the shock "quasi-parallel," and supporting the idea of an earlier "preconditioning" ICME. Third, the shock speed was estimated to be $\approx 3300$ km s$^{-1}$. Fourth, the sonic Mach number, $M_s$, for this shock was $\sim 28$. We support these results with an idealized numerical simulation of the ICME. Finally, we estimated the change in ram pressure upstream of the shock to be $\sim 5$ times larger than the pressure from the energetic particles, suggesting that this was not a standard "steady-state" cosmic-ray modified shock (CRMS). Instead it might represent an early, transient phase in the evolution of the CRMS.
Non-thermal emission has been detected in WR-stars for many years at long wavelengths spectral range, in general attributed to synchrotron emission. Two key ingredients are needed to explain such emissions, namely magnetic fields and relativistic particles. Particles can be accelerated to relativistic speeds by Fermi processes at strong shocks. Therefore, strong synchrotron emission is usually attributed to WR binarity. The magnetic field may also be amplified at shocks, however the actual picture of the magnetic field geometry, intensity, and its role on the acceleration of particles at WR binary systems is still unclear. In this work we discuss the recent developments in MHD modelling of wind-wind collision regions by means of numerical simulations, and the coupled particle acceleration processes related.
We present observational constraints on the initial-final mass relation (IFMR) using wide double white dwarfs (DWDs). We identify 65 new candidate wide DWDs within the Sloan Digital Sky Survey, bringing the number of candidate wide DWDs to 142. We then engage in a spectroscopic follow-up campaign and collect existing spectra for these objects; using these spectra, we derive masses and cooling ages for 54 hydrogen (DA) WDs in DWDs. We also identify one new DA/DB pair, four candidate DA/DC pairs, four candidate DA/DAH pairs, and one new candidate triple degenerate system. Because wide DWDs are co-eval and evolve independently, the difference in the pre-WD lifetimes should equal the difference in the WD cooling ages. We use this to develop a Bayesian hierarchical framework and construct a likelihood function to determine the probability that any particular IFMR fits a sample of wide DWDs. We then define a parametric model for the IFMR and find the best parameters indicated by our sample of DWDs. We place robust constraints on the IFMR for initial masses of 2--4 \Msun. The WD masses produced by our model for stars within this mass range differ from those predicted by semi-empirical fits to open cluster WDs. Within this mass range, where there are few constraining open cluster WDs and disagreements in the cluster ages, wide DWDs may provide more reliable constraints on the IFMR. Expanding this method to the many wide DWDs expected to be discovered by Gaia may transform our understanding of the IFMR.
OB type stars have strong ionizing radiation, and drive energetic winds. The ultraviolet (UV) radiation from ionizing stars may heat dust and ionize gas to sweep up an expanding bubble shell. This shell may be the result of feedback leading to a new generation of stars. N131 is an infrared dust bubble residing in a molecular filament. We study the formation and fragmentation of this bubble with multi-wavelength dust and gas observations. Towards the bubble N131, we analyzed archival multi-wavelength observations including 3.6, 4.5, 5.8, 8.0, 24, 70, 160, 250, 350, 500 \mu m, 1.1 mm, and 21 cm. In addition, we performed new observations of CO (2-1), CO (1-0), and ^{13}CO (1-0) with the IRAM 30-m telescope. Multi-wavelength dust and gas observations reveal a ringlike shell with compact fragments, two filamentary structures, and a secondary bubble N131-A. The bubble N131 is a rare object with a large hole at 24 \mu m and 21 cm in the direction of its center. The dust and gas clumps are compact and might have been compressed at the inner edge of the ringlike shell, while they are extended and might be pre-existing at the outer edge. The column density, excitation temperature, and velocity show a potentially hierarchical distribution from the inner to outer edge of the ringlike shell. We also detected the front and back sides of the secondary bubble N131-A in the direction of its center. The derived Lyman-continuum ionizing photon flux within N131-A is equivalent to an O9.5 star. Based on the above, we suggest that the bubble N131 might be triggered by the strong stellar winds from a group of massive stars inside the bubble. We propose a scenario in which the bubble N131 forms from the disruption of a gas filament by expansion of HII region, strong stellar winds, and fragments under self-gravity.
We estimate the age for the individual stars located at the lower part of the red giant branch from the LAMOST DR2 K giant sample. Taking into account the selection effects and the volume completeness, the age--metallicity map for the stars located between 0.3 and 1.5 kpc from the Sun is obtained. A significant substructure (denoted as the \it{narrow stripe}) located from (age, [Fe/H])$\sim$(5, 0.4) to (10 Gyr, -0.4 dex) in the age--metallicity map is clearly identified. Moreover, the \it{narrow stripe} stars are found the dominate contributors to several velocity substructures, including the well-known Hercules stream. The substantially large difference between the observed guiding-center radii and the birth radii inferred from the age--metallicity relation is evident that the \it{narrow stripe} stars have been radially migrated from about R$\sim4$ kpc to the solar neighborhood. This implies that the Hercules stream may not be owe to the resonance associated with the bar, but may be the kinematic imprint of the inner disk and later moved out due to radial migration. We estimate that the traveling speed of the radial migration are roughly 1.1$\pm0.1$ kpc Gyr$^{-1}$, equivalent with about $1.1\pm0.1$ km s$^{-1}$. This is in agreement with the median $v_R$ of $2.6^{+1.8}_{-1.9}$ km s$^{-1}$ of the \it{narrow stripe}. We also obtain that about one third stars in the solar neighborhood are radially migrated from around 4 kpc. Finally, we find that the radial migration does not lead to additional disk thickening according to the distribution of $z_{max}$.
We report on observations of the hydroxyl radical (OH) within The H{\sc I}, OH Recombination line survey (THOR) pilot region. The region is bounded approximately between Galactic coordinates l=29.2 to 31.5$^\circ$ and b=-1.0 to +1.0$^\circ$ and includes the high-mass star forming region W43. We identify 103 maser sites, including 72 with 1612\,MHz masers, 42 showing masers in either of the main line transitions at 1665 and 1667\,MHz and four showing 1720\,MHz masers. Most maser sites with either main-line or 1720\,MHz emission are associated with star formation, whereas most of the 1612\,MHz masers are associated with evolved stars. We find that nearly all of the main-line maser sites are co-spatial with an infrared source, detected by GLIMPSE. We also find diffuse OH emission, as well as OH in absorption towards selected unresolved or partially resolved sites. Extended OH absorption is found towards the well known star forming complex W43 Main.
We compute the complex 3D Fourier transform of the spatial galaxy distribution in a volume-limited sample of the Sloan Digital Sky Survey redshift survey. The direct unbinned transform yields results quite similar to those from the FFT of finely binned galaxy positions. In both cases the sampling window function is deconvolved to yield an estimate of the true 3D transform. The Fourier amplitudes resulting from this simple procedure yield power spectrum estimates consistent with those from other more complicated approaches. We display also measurements of homogeneity, isotropy and Gaussianity, based on an analysis of the complex Fourier transform that is much simpler than the multi-point methods usually employed. Our model-independent analysis avoids statistical interpretations, which have no meaning without detailed assumptions about a hypothetical process generating the initial cosmic density fluctuations.
Dome A in the Antarctic plateau is likely one of the best astronomical observing sites on Earth. The first one of three Antarctic Survey Telescope (AST3-1), a 50/68 cm Schmidt-like equatorial-mount telescope, is the first trackable telescope of China operating in Antarctica and the biggest telescope located in Antarctic inland. AST3-1 obtained huge amounts of data in 2012 and we processed the time-series parts. Here we present light curves of 29 variable stars identified from ten-day observations in 2012 with AST3-1, including 22 newly discovered variable stars. 23 of them are eclipsing binaries and the others are pulsating stars. We present the properties of the 29 variable stars, including the classifications, periods and magnitude ranges in i band. For the 17 eclipsing binaries, the phased light curves are presented with the orbital period values well determined.
A Quark-Nova (QN; the explosive transition of a Neutron star to a Quark star) occurring in the second common envelope (CE) phase of a massive binary, as described in Ouyed et al. (2015a&b), gives excellent fits to super-luminous, hydrogen-poor, Supernovae (SLSNe) with double-humped light curves including DES13S2cmm, SN 2006oz and LSQ14bdq (see {\it this http URL}). In our model, the hydrogen envelope of the less massive companion is ejected during the first CE phase while the QN explosion occurs deep inside the He-rich second CE phase after it has expanded to its equilibrium configuration at ~1200Rsun; this yields the first hump in our model. The subsequent merging of the quark star with the CO core leads to black hole formation and accretion explaining the second long-lasting hump in our model, while the collision of the QN-ejected He-rich CE with the H-rich (i.e. first) CE accounts for late emission. Here we show that our model provides an excellent fit to the recently discovered SN iPTF13ehe including the late-time lightcurve and provides an explanation for the related broad H_alpha emission. We show that LSQ14bdq and iPTF13ehe are very similar in nature and can be fit with nearly the same parameters. Because of the QN imprint (the first hump in the lightcurves), double-humped hydrogen-poor, SLSNe can potentially be used as standard candles.
(abridged) The dust content of the universe is primarily explored via its interaction with stellar photons, producing interstellar extinction. However, owing to the physical extension of the observing beam, observations may detect scattered photons, resulting in a change in the observed (or effective) extinction, depending on the spatial distribution of the dust and the resolution of the instrument. We investigate the influence of clumpy dust distributions on effective extinction toward embedded sources and those in the diffuse ISM. We use Monte Carlo radiative transfer to examine effective extinction for various geometries. By varying the number, optical depth and volume-filling factor of clumps in models of spherical shells and the diffuse ISM, we explore the evolution of extinction. Depending on the number of scattering events in the beam, the extinction curve steepens in homogeneous media and flattens in clumpy media. As a result, clumpy dust distributions can to reproduce extinction curves with arbitrary R_v, the effective ratio of total-to-selective extinction. The flattening `washes out' the 2175 \AA bump and shifts the peak to shorter wavelengths. The mean R_v of a shell correlates with the optical depth of an individual clump and the wavelength at which a clump becomes optically thick. Similar behaviour is seen for edge-on discs or tori. However, at grazing inclinations the combination of extinction and strong forward scattering results in chaotic behaviour. Caution is therefore advised when measuring extinction in, for example, AGN tori or toward SNIa or GRB afterglows. In face-on discs, the shape of the scattered continuum changes with clumpiness, however, unlike absorption features, individual features in the scattering cross-sections are preserved. Finally, we show that diffuse interstellar extinction is not modified by scattering on scales of a few kpc.
MWC 314 is one of the most luminous stars in the Milky Way. Its fundamental
parameters are similar to those of LBVs, although no large photometric
variations have been recorded.
We analyzed the radial velocity (RV) variations displayed by the absorption
lines from the star's atmosphere using high resolution optical spectra. The RV
and profile variations of some permitted and forbidden emission lines of
metallic ions and the behavior of the Balmer and HeI lines has been
investigated as well.
The RV of the absorption lines clearly shows a 60-day periodicity. A dense
coverage of the RV curve allowed us to derive accurate orbital parameters. The
RV of the FeII emission lines varies in the same way, but with a smaller
amplitude. Additionally, the intensity ratio of the blue/red peaks of these
emission lines correlates with the RV variations. The Balmer lines and {NII]
lines display a nearly constant RV and no profile variations in phase with the
orbital motion instead. The HeI-5876 \AA\ line shows a strongly variable
profile with broad and blue-shifted absorption components reaching velocities
of <-1000 km/s at some specific orbital phases. Our data strongly suggest that
the object is a binary system composed from a supergiant B[e] star and an
undetected companion. The emission lines with a non-variable RV could originate
in a circumbinary region. For the FeII emission lines we propose a simple
geometrical two-component model where a compact source of FeII emission, moving
around the center of mass, is affected by a static extra absorption that
originates from a larger area. Finally, the blue-shifted absorption in the
HeI-5876 \AA\ line could be the result of density enhancements in the primary
star wind flowing towards the companion that is best observed when projected
over the disk of the primary star.
We study the 205.9 Hz pulsations of the accreting millisecond X-ray pulsar NGC 6440 X-2 across all outbursts observed with the Rossi X-ray Timing Explorer over a period of 800 days. We find the pulsations are highly sinusoidal with a fundamental amplitude of 5%-15% rms and a second harmonic that is only occasionally detected with amplitudes of <2% rms. By connecting the orbital phase across multiple outbursts, we obtain an accurate orbital ephemeris for this source and constrain its 57 min orbital period to sub-ms precision. We do not detect an orbital period derivative to an upper limit of $ | \dot{P} | \leq 8 \times 10^{-11}$ s/s. We investigate the possibility of coherently connecting the pulse phase across all observed outbursts, but find that due to the poorly constrained systematic uncertainties introduced by a flux-dependent bias in the pulse phase, multiple statistically acceptable phase-connected timing solutions exist.
We study the role of the coronal magnetic field configuration of an active region in determining the propagation direction of a coronal mass ejection (CME). The CME occurred in the active region 11944 (S09W01) near the disk center on 2014 January 7 and was associated with an X1.2 flare. A new CME reconstruction procedure based on a polarimetric technique is adopted, which shows that the CME changed its propagation direction by around 28$^\circ$ in latitude within 2.5 R$_\odot$ and 43$^\circ$ in longitude within 6.5 R$_\odot$ with respect to the CME source region. This significant non-radial motion is consistent with the finding of M$\ddot{o}$stl et al. (2015). We use nonlinear force-free field (NLFFF) and potential field source surface (PFSS) extrapolation methods to determine the configurations of the coronal magnetic field. We also calculate the magnetic energy density distributions at different heights based on the extrapolations. Our results show that the active region coronal magnetic field has a strong influence on the CME propagation direction. This is consistent with the "channelling" by the active region coronal magnetic field itself, rather than deflection by nearby structures. These results indicate that the active region coronal magnetic field configuration has to be taken into account in order to determine CME propagation direction correctly.
In the early 1990s, contemporary interstellar dust (ISD) penetrating deep into the heliosphere was identified with the in-situ dust detector on board the Ulysses spacecraft. Between 1992 and the end of 2007 Ulysses monitored the ISD stream. The interstellar grains act as tracers of the physical conditions in the local interstellar medium surrounding our solar system. Earlier analyses of the Ulysses ISD data measured between 1992 and 1998 implied the existence of 'big' ISD grains [up to 10^-13kg]. The derived gas-to-dust-mass ratio was smaller than the one derived from astronomical observations, implying a concentration of ISD in the very local interstellar medium. We analyse the entire data set from 16 yr of Ulysses ISD measurements in interplanetary space. This paper concentrates on the overall mass distribution of ISD. An accompanying paper investigates time-variable phenomena in the Ulysses ISD data, and in a third paper we present the results from dynamical modelling of the ISD flow applied to Ulysses. We use the latest values for the interstellar hydrogen and helium densities, the interstellar helium flow speed of v_ISM,inf=23.2km/s, and the ratio of radiation pressure to gravity, beta, calculated for astronomical silicates. We find a gas-to-dust-mass ratio in the local interstellar cloud of R_g/d=193^+85_-57, and a dust density of 2.1+/-0.6x10^-24kg/m^3. For a higher inflow speed of 26km/s, the gas-to-dust-mass ratio is 20% higher, and, accordingly, the dust density is lower by the same amount. The gas-to-dust mass ratio derived from our new analysis is compatible with the value most recently determined from astronomical observations. We confirm earlier results that the very local interstellar medium contains 'big' (i.e. 1 um-sized) ISD grains. We find a dust density in the local interstellar medium that is a factor of three lower than values implied by earlier analyses.
We present the results of models that were designed to study all possible water maser transitions in the frequency range 0-1.91THz, with particular emphasis on maser transitions that may be generated in evolved-star envelopes and observed with the ALMA and SOFIA telescopes. We used tens of thousands of radiative transfer models of both spin species of H2O, spanning a considerable parameter space in number density, kinetic temperature and dust temperature. Results, in the form of maser optical depths, have been summarized in a master table, Table 6. Maser transitions identified in these models were grouped according to loci of inverted regions in the density/kinetic temperature plane, a property clearly related to the dominant mode of pumping. A more detailed study of the effect of dust temperature on maser optical depth enabled us to divide the maser transitions into three groups: those with both collisional and radiative pumping schemes (22,96,209,321,325,395,941 and 1486\,GHz), a much larger set that are predominantly radiatively pumped, and another large group with a predominantly collisional pump. The effect of accelerative and decelerative velocity shifts of up to 5km/s was found to be generally modest, with the primary effect of reducing computed maser optical depths. More subtle asymmetric effects, dependent on line overlap, include maximum gains offset from zero shift by >1km/s, but these effects were predominantly found under conditions of weak amplification. These models will allow astronomers to use multi-transition water maser observations to constrain physical conditions down to the size of individual masing clouds (size of a few astronomical units).
Blazars radiate from relativistic plasma jets with bulk Lorentz factors {\Gamma} ~ 10, closely aligned along our line of sight. In a number of blazars of the Flat Spectrum Radio Quasar type such as 3C 454.3 and 3C 279 gamma-ray flares have recently been detected with very high luminosity and little or no counterparts in the optical and soft X-ray bands. They challenge the current one-zone leptonic models of emissions from within the broad line region. The latter envisage the optical/X-ray emissions to be produced as synchrotron radiation by the same population of highly relativistic electrons in the jet that would also yield the gamma rays by inverse Compton up-scattering of surrounding soft photons. To meet the challenge we present here a model based on primary synchrotron photons emitted in the broad line region by a plasmoid moving out with the jet and scattered back toward the incoming plasmoid by an outer plasma clump acting as a mirror. We consider both a scenario based on a static mirror located outside the BLR, and an alternative provided by a moving mirror geometry. We show that mirroring phenomena can locally enhance the density and anisotropy with associated relativistic boosting of soft photons within the jet, so as to trigger bright inverse Compton gamma-ray transients from nearly steady optical/X-ray synchrotron emissions. In this picture we interpret the peculiarly asymmetric lightcurves of the recently detected gamma-ray flares from 3C 279. Our scenario provides a promising start to understand the widening class of bright and transient gamma-ray activities in blazars.
This paper is one of three companion papers presenting the results of our in-depth analysis of the interstellar neutral helium (ISN He) observations carried out using the IBEX-Lo during the first six Interstellar Boundary Explorer (IBEX) observation seasons. We derive corrections for losses due to the limited throughput of the interface buffer and determine the IBEX spin-axis pointing. We develop an uncertainty system for the data, taking into account the resulting correlations between the data points. This system includes uncertainties due to Poisson statistics, background, spin-axis determination, systematic deviation of the boresight from the prescribed position, correction for the interface buffer losses, and the expected Warm Breeze (WB) signal. Subsequently, we analyze the data from 2009 to examine the role of various components of the uncertainty system. We show that the ISN He flow parameters are in good agreement with the values obtained by the original analysis. We identify the WB as the principal contributor to the global $\chi^2$ values in previous analyses. Other uncertainties have a much milder role and their contributions are comparable to each other. The application of this uncertainty system reduced the minimum $\chi^2$ value 4-fold. The obtained $\chi^2$ value, still exceeding the expected value, suggests that either the uncertainty system may still be incomplete or the adopted physical model lacks a potentially important element, which is likely an imperfect determination of the WB parameters. The derived corrections and uncertainty system are used in the accompanying paper by Bzowski et al. in an analysis of the data from six seasons.
The Einstein Equivalence Principle is a fundamental principle of the theory of General Relativity. While this principle has been thoroughly tested with standard matter, the question of its validity in the Dark sector remains open. In this paper, we consider a general tensor-scalar theory that allows to test the equivalence principle in the Dark sector by introducing two different couplings to standard matter and to Dark matter. We constrain these couplings by considering galactic observations of strong lensing and of velocity dispersion. Our analysis shows that, in case of a violation of the Einstein Equivalence Principle, data favour violations through couplings to ordinary and Dark matters of opposite signs. At the same time, General Relativity remains perfectly compatible with observations at a 2-$\sigma$ confidence level.
Recent evidence - in particular the hard X-ray spectra obtained by NuSTAR, and the large amplitude hard X-ray variability observed when ultraluminous X-ray sources (ULXs) show soft spectra - reveals that common ULX behaviour is inconsistent with known sub-Eddington accretion modes, as would be expected for an intermediate-mass black hole (IMBH). Instead, it appears that the majority of ULXs are powered by super-Eddington accretion onto stellar-mass black holes. Here, we will review work that delves deeper into ULX spectral-timing behaviour, demonstrating it remains consistent with the expectations of super-Eddington accretion. One critical missing piece from this picture is the direct detection of the massive, radiatively-driven winds expected from ULXs as atomic emission/absorption line features in ULX spectra; we will show it is very likely these have already been detected as residuals in the soft X-ray spectra of ULXs. Finally, we will discuss ULXs that do not appear to conform to the emerging ULX behaviour patterns. In particular we discuss the implications of the identification of a good IMBH candidate as a background QSO; and the confirmation of an IMBH/ULX candidate in the galaxy NGC 2276 via the radio/X-ray fundamental plane.
This white paper examines the benefit of the upcoming James Webb Space Telescope for studies of the Solar System's four giant planets: Jupiter, Saturn, Uranus, and Neptune. JWST's superior sensitivity, combined with high spatial and spectral resolution, will enable near- and mid-infrared imaging and spectroscopy of these objects with unprecedented quality. In this paper we discuss some of the myriad scientific investigations possible with JWST regarding the giant planets. This discussion is preceded by the specifics of JWST instrumentation most relevant to giant planet observations. We conclude with identification of desired pre-launch testing and operational aspects of JWST that would greatly benefit future studies of the giant planets.
The formation and evolution of galaxies require large reservoirs of cold, neutral gas. The damped Lya systems (DLAs), seen in absorption towards distant quasars and gamma-ray bursts, are predicted to be the dominant reservoirs for this gas. Detailed properties of DLAs have been studied extensively for decades with great success. However, their size, fundamental in understanding their nature, has remained elusive, as quasar and gamma-ray-burst sightlines only probe comparatively tiny areas of the foreground DLAs. Here, we introduce a new approach to measure the full extent of DLAs in the sightlines toward extended background sources. We present the discovery of a high-column-density (log N(HI) = 21.1 +/-0.4 cm^-2) DLA at z ~ 2.4 covering 90-100% of the luminous extent of a line-of-sight background galaxy. Estimates of the size of the background galaxy range from a minimum of a few kpc^2, to ~100 kpc^2, and demonstrate that high-column density neutral gas can span continuous areas 10^8 - 10^10 times larger than previously explored in quasar or gamma-ray burst sightlines. The DLA presented here is the first from a sample of DLAs in our pilot survey that searches Lyman break and Lyman continuum galaxies at high redshift. The low luminosities, large sizes, and mass contents (>~10^6 - 10^9 M_solar) implied by this DLA and the early data suggest that DLAs contain the necessary fuel for galaxies, with many systems consistent with relatively massive, low-luminosity primeval galaxies.
The Earth's magnetic field is generated by dynamo action driven by convection in the outer core. Owing to the rapid rotation and small viscosity, the dynamical balance is believed to be between buoyancy, Coriolis and magnetic forces; inertial forces play no role. It is thus extremely important to produce explicit solutions with these features. However, from the traditional approach of solving the full governing equations at low Ekman numbers, it is not clear that the asymptotic regime has been captured. Here we adopt a complementary approach consisting of a model of rapidly rotating convection in which inertial forces are neglected from the outset. Within this framework we are able to construct a new branch of solutions in which the dynamo generates a strong magnetic field that satisfies the expected force balance. The resulting strongly magnetised convection is dramatically different to the corresponding solutions in which the magnetic field is weak.
[Abridged] Our aim is to explore the gas dynamics and the accretion process in the early phase of high-mass star formation. The inward motion of molecular gas in the massive star forming region G34.26+0.15 is investigated by using high-resolution profiles of seven transitions of ammonia at THz frequencies observed with Herschel-HIFI. The shapes and intensities of these lines are interpreted in terms of radiative transfer models of a spherical, collapsing molecular envelope. An accelerated Lambda Iteration (ALI) method is used to compute the models. The seven ammonia lines show mixed absorption and emission with inverse P-Cygni-type profiles that suggest infall onto the central source. A trend toward absorption at increasingly higher velocities for higher excitation transitions is clearly seen in the line profiles. The $J = 3\leftarrow2$ lines show only very weak emission, so these absorption profiles can be used directly to analyze the inward motion of the gas. This is the first time a multitransitional study of spectrally resolved rotational ammonia lines has been used for this purpose. Broad emission is, in addition, mixed with the absorption in the $1_0-0_0$ ortho-NH$_3$ line, possibly tracing a molecular outflow from the star forming region. The best-fitting ALI model reproduces the continuum fluxes and line profiles, but slightly underpredicts the emission and absorption depth in the ground-state ortho line $1_0-0_0$. The derived ortho-to-para ratio is approximately 0.5 throughout the infalling cloud core similar to recent findings for translucent clouds in sight lines toward W31C and W49N. We find evidence of two gas components moving inwards toward the central region with constant velocities: 2.7 and 5.3 km$\,$s$^{-1}$, relative to the source systemic velocity. The inferred mass accretion rates derived are sufficient to overcome the expected radiation pressure from G34.26+0.15.
The Fermi/LAT instrument has detected about two thousands Extragalactic High Energy (E > 100 MeV) gamma-ray sources. One of the brightest is 3FGL 1603.9-4903, associated to the radio source PMN J1603-4904. Its nature is not yet clear, it could be either a very peculiar BL Lac or a CSO (Compact Symmetric Object) radio source, considered as the early stage of a radio galaxy. The latter, if confirmed, would be the first detection in gamma-rays for this class of objects. Recently a redshift z=0.18 +/- 0.01 has been claimed on the basis of the detection of a single X-ray line at 5.44 +/- 0.05 keV interpreted as a 6.4 keV (rest frame) fluorescent line. We aim to investigate the nature of 3FGL 1603.9-4903/PMN J1603-4904 using optical to NIR spectroscopy. We observed PMN J1603-4904 with the UV-NIR VLT/X-shooter spectrograph for two hours. We extracted spectra in the VIS and NIR range that we calibrated in flux and corrected for telluric absorption and we systematically searched for absorption and emission features. The source was detected starting from ~6300 Ang down to 24000 Ang with an intensity comparable to the one of its 2MASS counterpart and a mostly featureless spectrum. The continuum lacks absorption features and thus is non-stellar in origin and likely non-thermal. On top of this spectrum we detected three emission lines that we interpret as the Halpha-[NII] complex, the [SII] 6716,6731 doublet and the [SIII] 9530 line, obtaining a redshift estimate of z= 0.2321 +/- 0.0004. The equivalent width of the Halpha-[NII] complex implies that PMN J1603-4904 does not follow the observational definition of BL Lac, the line ratios suggest that a LINER/Seyfert nucleus is powering the emission. This new redshift measurement implies that the X-ray line previously detected should be interpreted as a 6.7 keV line which is very peculiar.
The thermal Sunyaev-Zel'dovich (tSZ) effect is the inverse-Compton scattering of cosmic microwave background (CMB) photons off hot, ionized electrons, primarily located in galaxy groups and clusters. Recent years have seen immense improvement in our ability to probe cosmology and the astrophysics of the intracluster medium using the tSZ signal. Here, I describe cross-correlations of the tSZ effect measured in Planck data with gravitational lensing maps from Planck and the Canada-France-Hawaii Telescope Lensing Survey, as well as hydrodynamical simulations which show that such measurements do not probe "missing baryons," but rather the pressure of ionized gas in groups and clusters over a wide range of halo masses and redshifts. I also present recent measurements of higher-order tSZ statistics using data from the Atacama Cosmology Telescope, which yield strong constraints on the amplitude of density fluctuations. I describe stacking analyses of tSZ data from Planck, focusing on the behavior of the gas pressure in low-mass galaxy groups. I close with a prediction for the tSZ monopole, including relativistic corrections, which is the largest guaranteed spectral distortion signal for the proposed Primordial Inflation Explorer mission. The tSZ monopole will yield a direct measurement of the total thermal energy in ionized electrons in the observable universe.
A possibility is explored to account for the light curve and the low expansion velocity of the supernova SN~2011ht, a member of group of three objects showing signatures of both IIn and IIP supernovae. It is argued that the radiated energy and the expansion velocity are consistent with the low energy explosion (\approx6\times10^{49} erg) and \leq 2 M_{\odot} ejecta interacting with the circumstellar envelope of 6-8 M_{\odot} and the radius of ~2\times10^{14} cm. The test of this scenario is proposed.
This paper explores the driving of low-level hydrodynamical activity in protoplanetary-disc dead zones. A small adverse radial entropy gradient, ordinarily stabilised by rotation, excites oscillatory convection (`convective overstability') when thermal diffusion, or cooling, is neither too strong nor too weak. I revisit the linear theory of the instability, discuss its prevalence in protoplanetary discs, and show that unstable modes are exact nonlinear solutions in the local Boussinesq limit. Overstable modes cannot grow indefinitely, however, as they are subject to a secondary parametric instability that limits their amplitudes to relatively low levels. If parasites set the saturation level of the ensuing turbulence then the convective overstability is probably too weak to drive significant angular momentum transport or to generate vortices. But I also discuss an alternative, and far more vigorous, saturation route that generates radial `layers' or `zonal flows' (witnessed also in semiconvection). Numerical simulations are required to determine which outcome is favoured in realistic discs, and consequently how important the instability is for disc dynamics.
The potentially hazardous asteroid (PHA) (99942) Apophis is one of the most remarkable near-Earth asteroids (NEA) in terms of impact hazard. A good determination of its surface thermal inertia is very important in order to evaluate the Yarkovsky effect on its orbital evolution. We present thermal infrared observations obtained on January 29, 2013, with CanariCam mid-infrared camera/spectrograph attached to the Gran Telescopio CANARIAS (GTC, Roque de los Muchachos Observatory, La Palma, Spain) using the Si2-8.7, Si6-12.5, and Q1-17.65 filters with the aim of deriving Apophis' diameter ($D$), geometric albedo ($p_V$), and thermal inertia ($\Gamma$). We performed a detailed thermophysical model analysis of the GTC data combined with previously published thermal data obtained using Herschel Space Observatory PACS instrument at 70, 100, and 160 $\mu$m.The thermophysical model fit of the data favors low surface roughness solutions (within a range of roughness slope angles $rms$ between 0.1 and 0.5), and constrains the effective diameter, visible geometric albedo, and thermal inertia of Apophis to be $D_{eff} =$~380 -- 393 m, $p_V = $~0.24--0.33 (assuming absolute magnitude $H = 19.09 \pm 0.19$) and $\Gamma =$~50 -- 500 Jm$^{-2}$ s$^{-0.5}$ K$^{-1}$, respectively.
The XXIXth IAU General Assembly in Honolulu adopted IAU 2015 Resolution B2 on recommended zero points for the absolute and apparent bolometric magnitude scales. The resolution was proposed by the IAU Inter-Division A-G Working Group on Nominal Units for Stellar and Planetary Astronomy after consulting with a broad spectrum of researchers from the astronomical community. Resolution B2 resolves the long-standing absence of an internationally-adopted zero point for the absolute and apparent bolometric magnitude scales. Resolution B2 defines the zero point of the absolute bolometric magnitude scale such that a radiation source with $M_{\rm Bol}$ = 0 has luminosity L$_{\circ}$ = 3.0128e28 W. The zero point of the apparent bolometric magnitude scale ($m_{\rm Bol}$ = 0) corresponds to irradiance $f_{\circ}$ = 2.518021002e-8 W/m$^2$. The zero points were chosen so that the nominal solar luminosity (3.828e26 W) adopted by IAU 2015 Resolution B3 corresponds approximately to $M_{\rm Bol}$(Sun) = 4.74, the value most commonly adopted in recent literature. The nominal total solar irradiance (1361 W/m$^2$) adopted in IAU 2015 Resolution B3 corresponds approximately to apparent bolometric magnitude $m_{\rm bol}$(Sun) = -26.832. Implicit in the IAU 2015 Resolution B2 definition of the apparent bolometric magnitude scale is an exact definition for the parsec (648000/$\pi$ au) based on the IAU 2012 Resolution B2 definition of the astronomical unit.
X-ray emission from Young Stellar Objects (YSOs) is crucial to understand star formation. A very limited amount of X-ray results is available for the protostellar (ClassI) phase. A systematic search of transient X-ray phenomena combined with a careful evaluation of the evolutionary stage offer a widely unexplored window to our understanding of YSOs X-ray properties. Within the EXTraS project, a search for transients and variability in the whole XMM-Newton archive, we discover transient X-ray emission consistent with ISO-Oph 85, a strongly embedded YSO in the rho Ophiuchi region, not detected in previous time-averaged X-ray studies. We extract an X-ray light curve for the flare and determine its spectral parameters from XMM-Newton/EPIC (European Photon Imaging Camera) data using quantile analysis. The X-ray flare ($2500\,s$), the only one detected in the XMM-Newton archive for ISO-Oph 85, has a luminosity of $LogL_X[erg/s]=31.1$ and a spectrum consistent with a highly-absorbed one-component thermal model ($N_H=1.0^{+1.2}_{-0.5}10^{23}\,cm^{-2}$, $kT=1.15^{+2.35}_{-0.65}\,keV)$. We set an upper limit of $LogL_X[erg/s]<29.5$ to the quiescent X-ray luminosity. We build a SED with IR to mm photometry drawn from literature and mid-IR Spitzer and sub-mm Herschel photometry analysed by us, and compare it with pre-computed models. The sub-mm emission peak in the Herschel data suggests that the object is a ClassI protostar. However, the Herschel/IR position offset is larger than for other YSOs in the region, leaving some doubt on the association. This is the first X-ray flare from a YSO recognised as a candidate ClassI YSO via the analysis of its complete SED. This work shows how the analysis of the whole SED is fundamental for the classification of YSOs, and how the X-ray source detection techniques we developed can open a new era in time-resolved analysis of the X-ray emission from stars.
X-ray astronomy allows study of objects which may be associated with compact objects, i.e. neutron stars or black holes, and also may contain strong magnetic fields. Such objects are categorically non-spherical, and likely non-circular when projected on the sky. Polarization allows study of such geoemetric effects, and X-ray polarimetry is likely to become feasible for a significant number of sources in the future. A class of potential targets for future X-ray polarization observations is the high mass X-ray binaries (HMXBs), which consist of a compact object in orbit with an early type star. In this paper ws show that X-ray polarization from HMXBs has a distinct signature which depends on the source inclination and orbital phase. The presence of the X-ray source displaced from the star creates linear polarization even if the primary wind is spherically symmetric whenever the system is viewed away from conjunction. Direct X-rays dilute this polarization whenever the X-ray source is not eclipsed; at mid-eclipse the net polarization is expected to be small or zero if the wind is circularly symmetric around the line of centers. Resonance line scattering increases the scattering fraction, often by large factors, over the energy band spanned by resonance lines. Real winds are not expected to be spherically symmetric, or circularly symmetric around the line of centers, owing to the combined effects of the compact object gravity and ionization on the wind hydrodynamics. A sample calculation shows that this creates polarization fractions ranging up to tens of percent at mid-eclipse.
We present new accurate abundances for five neutron-capture (Y, La, Ce, Nd, Eu) elements in 73 classical Cepheids located across the Galactic thin disk. Individual abundances are based on high spectral resolution (R ~ 38,000) and high signal-to-noise ratio (S/N ~ 50-300) spectra collected with UVES at ESO VLT for the DIONYSOS project. Taking account for similar Cepheid abundances provided either by our group (111 stars) or available in the literature, we end up with a sample of 435 Cepheids covering a broad range in iron abundances (-1.6 < [Fe/H] < 0.6). We found, using homogeneous individual distances and abundance scales, well defined gradients for the above elements. However, the slope of the light s-process element (Y) is at least a factor of two steeper than the slopes of heavy s- (La, Ce, Nd) and r- (Eu) process elements. The s to r abundance ratio ([La/Eu]) of Cepheids shows a well defined anticorrelation with of both Eu and Fe. On the other hand, Galactic field stars attain an almost constant value and only when they approach solar iron abundance display a mild enhancement in La. The [Y/Eu] ratio shows a mild evidence of a correlation with Eu and, in particular, with iron abundance for field Galactic stars. We also investigated the s-process index - [hs/ls] - and we found a well defined anticorrelation, as expected, between [La/Y] and iron abundance. Moreover, we found a strong correlation between [La/Y] and [La/Fe] and, in particular, a clear separation between Galactic and Sagittarius red giants. Finally, the comparison between predictions for low-mass asymptotic giant branch stars and the observed [La/Y] ratio indicate a very good agreement over the entire metallicity range covered by Cepheids. However, the observed spread, at fixed iron content, is larger than predicted by current models.
We examine the evolution of the Parker instability in galactic disks using 3D numerical simulations. We consider a local Cartesian box section of a galactic disk, where gas, magnetic fields and cosmic rays are all initially in a magnetohydrostatic equilibrium. This is done for different choices of initial cosmic ray density and magnetic field. The growth rates and characteristic scales obtained from the models, as well as their dependences on the density of cosmic rays and magnetic fields, are in broad agreement with previous (linearized, ideal) analytical work. However, this non-ideal instability develops a multi-modal 3D structure, which cannot be quantitatively predicted from the earlier linearized studies. This 3D signature of the instability will be of importance in interpreting observations. As a preliminary step towards such interpretations, we calculate synthetic polarized intensity and Faraday rotation measure maps, and the associated structure functions of the latter, from our simulations; these suggest that the correlation scales inferred from rotation measure maps are a possible probe for the cosmic ray content of a given galaxy. Our calculations highlight the importance of cosmic rays in these measures, making them an essential ingredient of realistic models of the interstellar medium.
This year, a second generation of coronagraphs dedicated to high-contrast
direct imaging of exoplanets is starting operations. Among them, SPHERE,
installed at the focus of the UT3 Very Large Telescope, reaches unprecedented
contrast ratios up to $10^{-6}$ -$ 10^{-7}$, using eXtreme Adaptive Optics and
the Angular Differential Imaging (ADI) techniques.
In this paper, we present a new method called Keplerian-Stacker that improves
the detection limit of high contrast instruments like SPHERE, by up to a factor
of 10. It consists of observing a star on a long enough period to let a
hypothetical planet around that star move along its orbit. Even if in each
individual observation taken during one night, we do not detect anything, we
show that it is possible, using an optimization algorithm, to re-center the
images according to keplerian motions (ex: 10-100 images taken over a long
period of typically 1-10 years) and detect planets otherwise unreachable. This
method can be used in combination with the ADI technics (or possibly any other
high contrast data reduction method) to improve the Signal to Noise Ratio in
each individual image, and to further improve the global detection limit. It
also directly provides orbital parameters of the detected planets, as a
by-product of the optimization algorithm.
Galactic bulge microlensing surveys provide a probe of Galactic structure. We present the first field-by-field comparison between microlensing observations and the Besan\c{c}on population synthesis Galactic model. Using an updated version of the model we provide maps of optical depth, average event duration and event rate for resolved source populations and for difference imaging (DIA) events. We also compare the predicted event timescale distribution to that observed. The simulation follows the selection criteria of the MOA-II survey (Sumi et al. 2013). We modify the Besan\c{c}on model to include M dwarfs and brown dwarfs. Our best fit model requires a brown dwarf mass function slope of $-0.4$. The model provides good agreement with the observed average duration, and respectable consistency with the shape of the timescale distribution (reduced $\chi^2 \simeq 2.2$). The DIA and resolved source limiting yields bracket the observed number of events by MOA-II ($2.17\times$ and $0.83\times$ the number observed, respectively). We perform a 2-dimensional fit to the event spatial distribution to predict the optical depth and event rate across the Galactic bulge. The most serious difficulty for the model is that it provides only $\sim 50\%$ of the measured optical depth and event rate per star at low Galactic latitude around the inner bulge ($|b|<3{^\circ}$). This discrepancy most likely is associated with known under-estimated extinction and star counts in the innermost regions and therefore provides additional support for a missing inner stellar population.
Tidal forces close to massive black holes can violently disrupt stars that make a close approach. These extreme events are discovered via bright X-ray and optical/UV flares in galactic centers. Prior studies based on modeling decaying flux trends have been able to estimate broad properties, such as the mass accretion rate. Here we report the detection of flows of highly ionized X-ray gas in high-resolution X-ray spectra of a nearby tidal disruption event. Variability within the absorption-dominated spectra indicates that the gas is relatively close to the black hole. Narrow line widths indicate that the gas does not stretch over a large range of radii, giving a low volume filling factor. Modest outflow speeds of a few hundred kilometers per second are observed, significantly below the escape speed from the radius set by variability. The gas flow is consistent with a rotating wind from the inner, super-Eddington region of a nascent accretion disk, or with a filament of disrupted stellar gas near to the apocenter of an elliptical orbit. Flows of this sort are predicted by fundamental analytical theory and more recent numerical simulations.
We aim at clarifying the nature of the emission of two spatially related unidentified X-ray sources detected with XMM-Newton telescope at intermediate-low Galactic latitude. Observations reveal a point-like source aligned with elongated diffuse emission. The X-ray spectra are best-fitted by absorbed power laws with photon indices ~1.7 for the point-like and ~2.0 for the extended one. Both sources show nonthermal radio-continuum counterparts that might indicate a physical association. From the available data, we did not detect variability on the point-like source in several timescales. Two possible scenarios are analyzed: first, based on HI line absorption, assuming a Galactic origin, we infer a distance upper bound of <2 kpc, which poses a constraint on the height over the Galactic plane of <200 pc and on the linear size of the system of <2.3 pc. In this case, the X-ray luminosities are >10^32 erg/s and >7.5 x 10^32 erg/s, for the point-like and extended sources, respectively; second, an extra-Galactic nature is discussed, where the point-like source might be the core of a radio galaxy and the extended source its lobe. In this case, we compare derived fluxes, spectral indices, and spatial correlation with those typical from the radio galaxy population, showing the feasibility of this alternative astrophysical scenario. From the available observational evidence, we suggest that the most promising scenario to explain the nature of these sources is a system consisting of a one-sided radio galaxy, where the point-like source is an active galactic nucleus and the extended source corresponds to the emission from its lobe. Other possibilities include a PSR/PWN origin, where the radio/X-ray emission originates from the synchrotron cooling of relativistic particles in the PSR magnetic field or a casual alignment between two unrelated sources, such as an AGN core and a Galactic X-ray blob.
Rotation-powered "recycled" millisecond pulsars are a variety of rapidly-spinning neutron stars that typically show thermal X-ray radiation due to the heated surface of their magnetic polar caps. Detailed numerical modeling of the rotation-induced thermal X-ray pulsations observed from recycled millisecond pulsars, including all relevant relativistic and stellar atmospheric effects, has been identified as a promising approach towards an astrophysical determination of the true neutron star mass-radius relation, and by extension the state of cold matter at densities exceeding those of atomic nuclei. Herein, I review the basic model and methodology commonly used to extract information regarding neutron star structure from the pulsed X-ray radiation observed from millisecond pulsars. I also summarize the results of past X-ray observations of these objects and the prospects for precision neutron star mass-radius measurements with the upcoming Neutron Star Interior Composition Explorer (NICER) X-ray timing mission.
The Solar System contains a population of dust and small particles originating from asteroids, comets, and other bodies. These particles have been studied using a number of techniques ranging from in-situ satellite detectors to analysis of lunar microcraters to ground-based observations of zodiacal light. In this paper, we describe an approach for using the LISA Pathfinder (LPF) mission as an instrument to detect and characterize the dynamics of dust particles in the vicinity of Earth-Sun L1. Launching in late 2015, LPF is a dedicated technology demonstrator mission that will validate several key technologies for a future space-based gravitational-wave observatory. The primary science instrument aboard LPF is a precision accelerometer which we show will be capable of sensing discrete momentum impulses as small as $4\times 10^{-8}\,\textrm{N}\cdot\textrm{s}$. We then estimate the rate of such impulses resulting from impacts of micrometeoroids based on standard models of the micrometeoroid environment in the inner solar system. We find that LPF may detect dozens to hundreds of individual events corresponding to impacts of particles with masses $> 10^{-9}\,$g during LPF's roughly six-month science operations phase in a $5\times 10^5\,\textrm{km}$ by $8\times 10^5\,\textrm{km}$ Lissajous orbit around L1. In addition, we estimate the ability of LPF to characterize individual impacts by measuring quantities such as total momentum transferred, direction of impact, and location of impact on the spacecraft. Information on flux and direction provided by LPF may provide insight as to the nature and origin of the individual impact and help constrain models of the interplanetary dust complex in general. Additionally, this direct in-situ measurement of micrometeoroid impacts will be valuable to designers of future spacecraft targeting the environment around L1.
White dwarfs are the end state of most stars, including the Sun, after they exhaust their nuclear fuel. Between 1/4 and 1/2 of white dwarfs have elements heavier than helium in their atmospheres, even though these elements should rapidly settle into the stellar interiors unless they are occasionally replenished. The abundance ratios of heavy elements in white dwarf atmospheres are similar to rocky bodies in the Solar system. This and the existence of warm dusty debris disks around about 4% of white dwarfs suggest that rocky debris from white dwarf progenitors' planetary systems occasionally pollute the stars' atmospheres. The total accreted mass can be comparable to that of large asteroids in the solar system. However, the process of disrupting planetary material has not yet been observed. Here, we report observations of a white dwarf being transited by at least one and likely multiple disintegrating planetesimals with periods ranging from 4.5 hours to 4.9 hours. The strongest transit signals occur every 4.5 hours and exhibit varying depths up to 40% and asymmetric profiles, indicative of a small object with a cometary tail of dusty effluent material. The star hosts a dusty debris disk and the star's spectrum shows prominent lines from heavy elements like magnesium, aluminium, silicon, calcium, iron, and nickel. This system provides evidence that heavy element pollution of white dwarfs can originate from disrupted rocky bodies such as asteroids and minor planets.
We supplement the minimal supersymmetric standard model (MSSM) with vector-like copies of standard model particles. Such 4th generation particles can raise the Higgs boson mass to the observed value without requiring very heavy superpartners, improving naturalness and the prospects for discovering supersymmetry at the LHC. Here we show that these new particles are also motivated cosmologically: in the MSSM, pure Bino dark matter typically overcloses the Universe, but 4th generation particles open up new annihilation channels, allowing Binos to have the correct thermal relic density without resonances or co-annihilation. We show that this can be done in a sizable region of parameter space while preserving gauge coupling unification and satisfying constraints from collider, Higgs, precision electroweak, and flavor physics.
The simplest Higgs-portal dark matter model is studied in the light of dark matter self-interacting effects on the formation of large scale structures. We show the direct detection limits on the resonant and large mass regions. Finally, we also compare these limits with those at the LHC and Xenon 1T experiments.
We study the model for the magnetic field generation in a magnetar based on the magnetic field instability driven by the parity violating electroweak interaction between electrons and nucleons in the neutron star matter. Using the quantum field theory methods, we calculate the helicity flip of an electron scattering off protons in dense matter of a neutron star. The influence of the electroweak interaction between electrons and background nucleons on the helicity flip is examined. We also derive the kinetic equation for the chiral imbalance. The evolution of the magnetic field in a magnetar accounting for the correct value of the helicity flip rate is studied.
The space-based gravitational-wave detector eLISA has been selected as the ESA L3 mission, and the mission design will be finalised by the end of this decade. To prepare for mission formulation over the next few years, several outstanding and urgent questions in data analysis will be addressed using mock data challenges, informed by instrument measurements from the LISA Pathfinder satellite launching at the end of 2015. These data challenges will require accurate and computationally affordable waveform models for anticipated sources such as the extreme-mass-ratio inspirals (EMRIs) of stellar-mass compact objects into massive black holes. Previous data challenges have made use of the well-known analytic EMRI waveforms of Barack and Cutler, which are extremely quick to generate but dephase relative to more accurate waveforms within hours, due to their mismatched radial, polar and azimuthal frequencies. In this paper, we describe an augmented Barack-Cutler model that uses a frequency map to the correct Kerr frequencies, along with updated evolution equations and a simple fit to a more accurate model. The augmented waveforms stay in phase for months and may be generated with virtually no additional computational cost.
The possibility that dark matter may be in the form of a Bose-Einstein Condensate (BEC) has been extensively explored at galactic scale. In particular, good fits for the galactic rotations curves have been obtained, and upper limits for the dark matter particle mass and scattering length have been estimated. In the present paper we extend the investigation of the properties of the BEC dark matter to the galactic cluster scale, involving dark matter dominated astrophysical systems formed of thousands of galaxies each. By considering that one of the major components of a galactic cluster, the intra-cluster hot gas, is described by King's $\beta$-model, and that both intra-cluster gas and dark matter are in hydrostatic equilibrium, bound by the same total mass profile, we derive the mass and density profiles of the BEC dark matter. In our analysis we consider several theoretical models, corresponding to isothermal hot gas and zero temperature BEC dark matter, non-isothermal gas and zero temperature dark matter, and isothermal gas and finite temperature BEC, respectively. The properties of the finite temperature BEC dark matter cluster are investigated in detail numerically. We compare our theoretical results with the observational data of 106 galactic clusters. Using a least-squares fitting, as well as the observational results for the dark matter self-interaction cross section, we obtain some upper bounds for the mass and scattering length of the dark matter particle. Our results suggest that the mass of the dark matter particle is of the order of $\mu $eV, while the scattering length has values in the range of $10^{-7}$ fm.
In this work we review the role of hyperons on the properties of neutron and proto-neutron stars. In particular, we revise the so-called "hyperon puzzle", go over some of the solutions proposed to tackle it, and discuss the implications that the recent measurements of unusually high neutron star masses have on our present knowledge of hypernuclear physics. We reexamine also the role of hyperons on the cooling properties of newly born neutron stars and on the so-called r-mode instability.
The KWISP opto-mechanical force sensor has been built and calibrated in the INFN Trieste optics laboratory and is now under off-beam commissioning at CAST. It is designed to detect the pressure exerted by a flux of solar Chameleons on a thin (100 nm) Si$_3$N$_4$ micromembrane thanks to their direct coupling to matter. A thermally-limited force sensitivity of $1.5 \cdot 10^{-14}~\mbox{N}/\sqrt{\mbox{Hz}}$, corresponding to $7.5 \cdot 10^{-16}~\mbox{m}/\sqrt{\mbox{Hz}}$ in terms of displacement, has been obtained. An originally developed prototype chameleon chopper has been used in combination with the KWISP force sensor to conduct preliminary searches for solar chamaleons.
Links to: arXiv, form interface, find, astro-ph, recent, 1510, contact, help (Access key information)
In recent years a number of poorly studied chemical elements, such as phosphorus, sulphur, and strontium, have received special attention as important tracers of the Galactic chemical evolution. By exploiting the capabilities of the infrared echelle spectrograph GIANO mounted at the Telescopio Nazionale Galileo, we acquired high resolution spectra of four Galactic dwarf stars spanning the metallicity range between about one-third and twice the solar value. We performed a detailed feasibility study about the effectiveness of the P, S, and Sr line diagnostics in the Y band between 1.03 and 1.10 microm. Accurate chemical abundances have been derived using one-dimensional model atmospheres computed in local thermodynamic equilibrium (LTE). We computed the line formation assuming LTE for P, while we performed non-LTE analysis to derive S and Sr abundances. We were able to derive phosphorus abundance for three stars and an upper limit for one star, while we obtained the abundance of sulphur and strontium for all of the stars. We find [P/Fe] and [S/Fe] abundance ratios consistent with solar-scaled or slightly depleted values, while the [Sr/Fe] abundance ratios are more scattered (by +/-0.2 dex) around the solar-scaled value. This is fully consistent with previous studies using both optical and infrared spectroscopy. We verified that high-resolution, Y-band spectroscopy as provided by GIANO is a powerful tool to study the chemical evolution of P, S, and Sr in dwarf stars.
Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients, and radio emission. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger gravitational wave signals and their neutrinos and electromagnetic counterparts. In this paper, we employ the SpEC code to simulate the merger of low-mass neutron star binaries (two $1.2M_\odot$ neutron stars) for a set of three nuclear-theory based, finite temperature equations of state. We show that the frequency peaks of the post-merger gravitational wave signal are in good agreement with predictions obtained from simulations using a simpler treatment of gravity. We find, however, that only the fundamental mode of the remnant is excited for long periods of time: emission at the secondary peaks is damped on a millisecond timescale in the simulated binaries. For such low-mass systems, the remnant is a massive neutron star which, depending on the equation of state, is either permanently stable or long-lived. We observe strong excitations of l=2, m=2 modes, both in the massive neutron star and in the form of hot, shocked tidal arms in the surrounding accretion torus. We estimate the neutrino emission of the remnant using a neutrino leakage scheme and, in one case, compare these results with a gray two-moment neutrino transport scheme. We confirm the complex geometry of the neutrino emission, also observed in previous simulations with neutrino leakage, and show explicitly the presence of important differences in the neutrino luminosity, disk composition, and outflow properties between the neutrino leakage and transport schemes.
The star HII 2407 is a member of the relatively young Pleiades star cluster and was previously discovered to be a single-lined spectroscopic binary. It is newly identified here within $Kepler$/$K2$ photometric time series data as an eclipsing binary system. Mutual fitting of the radial velocity and photometric data leads to an orbital solution and constraints on fundamental stellar parameters. While the primary has arrived on the main sequence, the secondary is still pre-main-sequence and we compare our results for the $M/M_\odot$ and $R/R_\odot$ values with stellar evolutionary models. We also demonstrate that the system is likely to be tidally synchronized. Follow-up infrared spectroscopy is likely to reveal the lines of the secondary, allowing for dynamically measured masses and elevating the system to benchmark eclipsing binary status.
We conduct precise strong lensing mass modeling of four ${\it Hubble}$ Frontier Fields (HFF) clusters, Abell$~$2744, MACS$~$J0416.1$-$2403, MACS$~$J0717.5$+$3745, and MACS$~$J1149.6$+$2223, for which HFF imaging observations are completed. We construct a refined sample of more than 100 multiple images for each cluster by taking advantage of the full depth HFF images, and conduct mass modeling using the ${\small \rm GLAFIC}$ software, which assumes simply parametrized mass distributions. Our mass modeling also exploits a magnification constraint from the lensed Type Ia supernova HFF14Tom for Abell$~$2744 and positional constraints from multiple images of the lensed supernova SN Refsdal for MACS$~$J1149.6$+$2223. We find that our best-fitting mass models reproduce the observed image positions with RMS errors of $\sim 0.4$ arcsec, which are smaller than RMS errors in previous mass modeling that adopted similar numbers of multiple images. We then construct catalogs of $z\sim 6-9$ dropout galaxies behind the four clusters and estimate magnification factors for these dropout galaxies with our best-fitting mass models. The dropout sample from the four cluster fields contains $\sim 120$ galaxies at $z\gtrsim 6$, about 20 of which are predicted to be magnified by a factor of more than 10. Some of the high-redshift galaxies detected in the HFF have lensing-corrected magnitudes of $M_{\rm UV}\sim -15$ to $-14$. Our analysis demonstrates that the HFF data indeed offer an ideal opportunity to study faint high-redshift galaxies. All lensing maps produced from our mass modeling will be made available on the STScI website.
We have discovered a new H$_2$CO (formaldehyde) $1_{1,0}-1_{1,1}$ 4.82966 GHz maser in Galactic Center Cloud C, G0.38+0.04. At the time of submission, this is the eighth region containing an H$_2$CO maser detected in the Galaxy. Cloud C is one of only two sites of confirmed high-mass star formation along the Galactic Center Ridge, affirming that H$_2$CO masers are exclusively associated with high-mass star formation. This discovery led us to search for other masers, among which we found new SiO vibrationally excited masers, making this the fourth star-forming region in the Galaxy to exhibit SiO maser emission. Cloud C is also a known source of CH$_3$OH Class-II and OH maser emission. There are now two known SiO and H$_2$CO maser containing regions in the CMZ, compared to two and six respectively in the Galactic disk, while there is a relative dearth of H$_2$O and CH$_3$OH Class-II masers in the CMZ. SiO and H$_2$CO masers may be preferentially excited in the CMZ, perhaps due to higher gas-phase abundances from grain destruction and heating, or alternatively H$_2$O and CH$_3$OH maser formation may be suppressed in the CMZ. In any case, Cloud C is a new testing ground for understanding maser excitation conditions.
We present a new exploratory framework to model galaxy formation and evolution in a hierarchical universe by using machine learning (ML). Our motivations are two-fold: (1) presenting a new, promising technique to study galaxy formation, and (2) quantitatively analyzing the extent of the influence of dark matter halo properties on galaxies in the backdrop of semi-analytical models (SAMs). We use the influential Millennium Simulation and the corresponding Munich SAM to train and test various sophisticated machine learning algorithms (k-Nearest Neighbors, decision trees, random forests and extremely randomized trees). By using only essential dark matter halo physical properties for haloes of $M>10^{12} M_{\odot}$ and a partial merger tree, our model predicts the hot gas mass, cold gas mass, bulge mass, total stellar mass, black hole mass and cooling radius at z = 0 for each central galaxy in a dark matter halo for the Millennium run. Our results provide a unique and powerful phenomenological framework to explore the galaxy-halo connection that is built upon SAMs and demonstrably place ML as a promising and a computationally efficient tool to study small-scale structure formation.
There is growing observational evidence of high-redshift quasars launching energetic, fast outflows, but the effects that these have on their host galaxies is poorly understood. We employ the moving-mesh code AREPO to study the feedback from a quasar that has grown to $\sim 10^9 M_\odot$ by $z \sim 5$ and the impact that this has on its host galaxy. Our simulations use a super-Lagrangian refinement technique to increase the accuracy with which the interface of the quasar-driven wind and the surrounding gas is resolved. We find that the feedback injected in these simulations is less efficient at removing gas from the galaxy than in previous work using the same feedback strength, and that this leads to the growth of a massive, rotationally supported, star-forming disc, co-existing with a powerful quasar-driven outflow. The properties of our host galaxy, including the kinematical structure of the gaseous disc and of the outflow, are in good agreement with current observations. Upcoming ALMA and JWST observations will be an excellent test of our model and will provide further clues as to the variance in properties of high-redshift quasar hosts.
At present, the principal limitation on the resolution and contrast of astronomical imaging instruments comes from aberrations in the optical path, which may be imposed by the Earth's turbulent atmosphere or by variations in the alignment and shape of the telescope optics. These errors can be corrected physically, with active and adaptive optics, and in post-processing of the resulting image. A recently-developed adaptive optics post-processing technique, called kernel phase interferometry, uses linear combinations of phases that are self-calibrating with respect to small errors, with the goal of constructing observables that are robust against the residual optical aberrations in otherwise well-corrected imaging systems. Here we present a direct comparison between kernel phase and the more established competing techniques, aperture masking interferometry, point spread function (PSF) fitting and bispectral analysis. We resolve the alpha Ophiuchi binary system near periastron, using the Palomar 200-Inch Telescope. This is the first case in which kernel phase has been used with a full aperture to resolve a system close to the diffraction limit with ground-based extreme adaptive optics observations. Excellent agreement in astrometric quantities is found between kernel phase and masking, and kernel phase significantly outperforms PSF fitting and bispectral analysis, demonstrating its viability as an alternative to conventional non-redundant masking under appropriate conditions.
Warm jupiters are an unexpected population of extrasolar planets that are too
near to their host to have formed in situ, but distant enough to retain a
significant eccentricity in the face of tidal damping. These planets are
curiously absent around stars larger than two solar radii. We hypothesize that
the warm jupiters are migrating due to Kozai-Lidov oscillations, which leads to
transient episodes of high eccentricity and a consequent tidal decay.
As their host evolves, such planets would be rapidly dragged in or engulfed
at minimum periapse, leading to a rapid depletion of the population with
increasing stellar radius, as is observed. Using numerical simulations, we
determine the relationship between periapse distance and orbital migration rate
for planets 0.1 to 10 Jupiter masses and with orbital periods between 10 and
100 days. We find that Kozai-Lidov oscillations effectively result in planetary
removal early in the evolution of the host star, possibly accounting for the
observed deficit. While the observed eccentricity distribution is inconsistent
with the simulated distribution for an oscillating and migrating warm jupiter
population, observational biases may explain the discrepancy.
We use particle data from the Illustris simulation, combined with individual kinematic constraints on the mass of the Milky Way (MW) at specific distances from the Galactic center, to infer the radial distribution of the MW's dark matter halo mass. Our method allows us to convert any constraint on the mass of the MW within a fixed distance to a full circular velocity profile to the MW's virial radius. As primary examples, we take two recent measurements of the total mass within 50 kpc of the Galaxy -- $4.2\times 10^{11}\,M_{\odot}$ (Deason et al. 2012) and $2.9\times 10^{11}\,M_{\odot}$ (Gibbons et al. 2014) -- and find they imply very different mass profiles and stellar masses for the Galaxy. The dark-matter-only version of the Illustris simulation enables us to compute the effects of galaxy formation on such constraints on a halo-by-halo basis; on small scales, galaxy formation enhances the density relative to dark-matter-only runs, while the total mass density is approximately 20% lower at large Galactocentric distances. We are also able to quantify how current and future constraints on the mass of the MW at specific radii will be reflected in uncertainties on its virial mass: even a measurement of M(<50 kpc) with essentially perfect precision still results in a 20% uncertainty on the virial mass of the Galaxy, while a future measurement of M(<100 kpc) with 10% errors would result in the same level of uncertainty. We expect that our technique will become even more useful as (1) better kinematic constraints become available at larger distances and (2) cosmological simulations provide even more faithful representations of the observable Universe.
We use cosmological simulations from the Feedback In Realistic Environments (FIRE) project, which implement a comprehensive set of stellar feedback processes, to study ultra-violet (UV) metal line emission from the circum-galactic medium of high-redshift (z = 2-4) galaxies. Our simulations cover the halo mass range Mh~2x10^11 - 8.5x10^12 Msun at z = 2, representative of Lyman break galaxies. Of the transitions we analyze, the low-ionization C III (977 A) and Si III (1207 A) emission lines are the most luminous, with C IV (1548 A) and Si IV (1394 A) also showing interesting spatially-extended structures that should be detectable by current and upcoming integral field spectrographs such as the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope and Keck Cosmic Web Imager (KCWI). The more massive halos are on average more UV-luminous. The UV metal line emission from galactic halos in our simulations arises primarily from collisionally ionized gas and is strongly time variable, with peak-to-trough variations of up to ~2 dex. The peaks of UV metal line luminosity correspond closely to massive and energetic mass outflow events, which follow bursts of star formation and inject sufficient energy into galactic halos to power the metal line emission. The strong time variability implies that even some relatively low-mass halos may be detectable in deep observations with current generation instruments. Conversely, flux-limited samples will be biased toward halos whose central galaxy has recently experienced a strong burst of star formation.
The Next-Generation Very Large Array (ngVLA) will be critical for understanding how galaxies are built and evolve at the earliest epochs. The sensitivity and frequency coverage will allow for the detection of cold gas and dust in `normal' distant galaxies, including the low-J transitions of molecular gas tracers such as CO, HNC, and HCO+; synchrotron and free-free continuum emission; and even the exciting possibility of thermal dust emission at the highest (z~7) redshifts. In particular, by enabling the total molecular gas reservoirs to be traced to unprecedented sensitivities across a huge range of epochs simultaneously -- something no other radio or submillimeter facility will be capable of -- the detection of the crucial low-J transitions of CO in a diverse body of galaxies will be the cornerstone of ngVLA's contribution to high-redshift galaxy evolution science. The ultra-wide bandwidths will allow a complete sampling of radio SEDs, as well as the detection of emission lines necessary for spectroscopic confirmation of elusive dusty starbursts. The ngVLA will also deliver unique contributions to our understanding of cosmic magnetism and to science accessible through microwave polarimetry. Finally, the superb angular resolution will move the field beyond detection experiments and allow detailed studies of the morphology and dynamics of these systems, including dynamical modeling of disks/mergers, determining the properties of outflows, measuring black hole masses from gas disks, and resolving multiple AGN nuclei. We explore the contribution of a ngVLA to these areas and more, as well as synergies with current and upcoming facilities including ALMA, SKA, large single-dish submillimeter observatories, GMT/TMT, and JWST.
Large and magnetically complex sunspot groups are known to be associated with flares. To date, the Mount Wilson scheme has been used to classify sunspot groups based on their morphological and magnetic properties. The most flare prolific class, the delta sunspot-group, is characterised by opposite polarity umbrae within a common penumbra, separated by less than 2 degrees. In this article, we present a new system, called the Solar Monitor Active Region Tracker - Delta Finder (SMART-DF), that can be used to automatically detect and classify magnetic deltas in near-realtime. Using continuum images and magnetograms from the Helioseismic and Magnetic Imager (HMI) onboard NASA's Solar Dynamics Observatory (SDO), we first estimate distances between opposite polarity umbrae. Opposite polarity pairs having distances of less that 2 degrees are then identified, and if these pairs are found to share a common penumbra, they are identified as a magnetic delta configuration. The algorithm was compared to manual delta detections reported by the Space Weather Prediction Center (SWPC), operated by the National Oceanic and Atmospheric Administration (NOAA). SMART-DF detected 21 out of 23 active regions (ARs) that were marked as delta spots by NOAA during 2011 - 2012 (within +/- 60 degrees longitude). SMART-DF in addition detected five ARs which were not announced as delta spots by NOAA. The near-relatime operation of SMART-DF resulted in many deltas being identified in advance of NOAA's daily notification. SMART-DF will be integrated with SolarMonitor (www.solarmonitor.org) and the near-realtime information will be available to the public.
We present radial velocities, stellar parameters, and detailed abundances of 39 elements derived from high-resolution spectroscopic observations of red giant stars in the luminous, metal-poor globular cluster NGC 5824. We observe 26 stars in NGC 5824 using the Michigan/Magellan Fiber System (M2FS) and two stars using the Magellan Inamori Kyocera Echelle (MIKE) spectrograph. We derive a mean metallicity of [Fe/H]=-1.94+/-0.02 (statistical) +/-0.10 (systematic). The metallicity dispersion of this sample of stars, 0.08 dex, is in agreement with previous work and does not exceed the expected observational errors. Previous work suggested an internal metallicity spread only when fainter samples of stars were considered, so we cannot exclude the possibility of an intrinsic metallicity dispersion in NGC 5824. The M2FS spectra reveal a large internal dispersion in [Mg/Fe], 0.28 dex, which is found in a few other luminous, metal-poor clusters. [Mg/Fe] is correlated with [O/Fe] and anti-correlated with [Na/Fe] and [Al/Fe]. There is no evidence for internal dispersion among the other alpha- or Fe-group abundance ratios. Twenty-five of the 26 stars exhibit a n-capture enrichment pattern dominated by r-process nucleosynthesis ([Eu/Fe]=+0.11+/-0.12; [Ba/Eu]=-0.66+/-0.05). Only one star shows evidence of substantial s-process enhancement ([Ba/Fe]=+0.56+/-0.12; [Ba/Eu]=+0.38+/-0.14), but this star does not exhibit other characteristics associated with s-process enhancement via mass-transfer from a binary companion. The Pb and other heavy elements produced by the s-process suggest a timescale of no more than a few hundred Myr for star formation and chemical enrichment, like the complex globular clusters M2, M22, and NGC 5286.
We investigate the gamma-ray and X-ray properties of the flat spectrum radio quasar PKS 2149-306 at redshift z = 2.345. A strong gamma-ray flare from this source was detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope satellite in 2013 January, reaching on January 20 a daily peak flux of (301$\pm$36)$\times$10$^{-8}$ ph/cm$^2$/s in the 0.1-100 GeV energy range. This flux corresponds to an apparent isotropic luminosity of (1.5$\pm$0.2)$\times$10$^{50}$ erg/s, comparable to the highest values observed by a blazar so far. During the flare the increase of flux was accompanied by a significant change of the spectral properties. Moreover significant flux variations on a 6-h time-scale were observed, compatible with the light crossing time of the event horizon of the central black hole. The broad band X-ray spectra of PKS 2149-306 observed by Swift-XRT and NuSTAR are well described by a broken power-law model, with a very hard spectrum ($\Gamma$$_1$ $\sim$ 1) below the break energy, at E$_{\rm\,break}$ = 2.5-3.0 keV, and $\Gamma$$_2$ $\sim$ 1.4-1.5 above the break energy. The steepening of the spectrum below $\sim$ 3 keV may indicate that the soft X-ray emission is produced by the low-energy relativistic electrons. This is in agreement with the small variability amplitude and the lack of spectral changes in that part of the X-ray spectrum observed between the two NuSTAR and Swift joint observations. As for the other high-redshift FSRQ detected by both Fermi-LAT and Swift-BAT, the photon index of PKS 2149-306 in hard X-ray is 1.6 or lower and the average gamma-ray luminosity higher than 2$\times$10$^{48}$ erg/s.
We propose using radio line spectroscopy to detect molecular absorption lines (such as OH at 1.6-1.7 GHz) before and after the total eclipse of black widow (BW) and other short orbital period binary pulsars with low mass companions. The companion in such a binary may be ablated away by energetic particles and high energy radiation produced by the pulsar wind. The observations will probe the eclipsing wind being ablated by the pulsar and constrain the nature of the companion and its surroundings. Maser emission from the interstellar medium stimulated by a pulsar beam might also be detected from the intrabinary medium. The short temporal resolution allowed by the millisecond pulsars can probe this medium with the high angular resolution of the pulsar beam.
In light of the growing number of cosmological observations, it is important to develop versatile tools to quantify the constraining power and consistency of cosmological probes. Originally motivated from information theory, we use the relative entropy to compute the information gained by Bayesian updates in units of bits. This measure quantifies both the improvement in precision and the 'surprise', i.e. the tension arising from shifts in central values. Our starting point is a WMAP9 prior which we update with observations of the distance ladder, supernovae (SNe), baryon acoustic oscillations (BAO), and weak lensing as well as the 2015 Planck release. We consider the parameters of the flat $\Lambda$CDM concordance model and some of its extensions which include curvature and Dark Energy equation of state parameter $w$. We find that, relative to WMAP9 and within these model spaces, the probes that have provided the greatest gains are Planck (10 bits), followed by BAO surveys (5.1 bits) and SNe experiments (3.1 bits). The other cosmological probes, including weak lensing (1.7 bits) and $\rm H_0$ measures (1.7 bits), have contributed information but at a lower level. Furthermore, we do not find any significant surprise when updating the constraints of WMAP9 with any of the other experiments, meaning that they are consistent with WMAP9. However, when we choose Planck15 as the prior, we find that the weak lensing measurements of CFHTLenS produce a large surprise of 4.4 bits which is statistically significant at the 8 $\sigma$ level. We discuss how the relative entropy provides a versatile and robust framework to compare cosmological probes in the context of current and future surveys.
We use a hydro-and-radiative-transfer code to explore the behavior of a very massive star (VMS) after a giant eruption -- i.e., following a supernova impostor event. Beginning with a reasonable model for an evolved VMS, we simulate the change of state caused by a giant eruption via two methods that explicitly conserve total energy: 1. Synthetically removing outer layers of mass while reducing the energy of the inner layers. 2. Synthetically transferring energy from the core to the outer layers, an operation that automatically causes mass ejection. Our focus is on the aftermath, not the poorly-understood eruption itself. Then, using a radiation-hydrodynamic code in 1D with realistic opacities and convection, the interior disequilibrium state is followed for about 200 years. Typically the star develops a $\sim 400 ~\rm{km}~\rm{s}^{-1}$ wind with a mass loss rate that begins around $0.1 ~M_\odot~\rm{yr^{-1}}$ and gradually decreases. This outflow is driven by $\kappa$-mechanism radial pulsations. In some cases a plateau in the mass loss rate may persist about 200 years, while other cases are more like $\eta$ Car's known history. These simulations constitute a useful preliminary reconnaissance for 3D models which will be far more difficult.
This white paper discusses how a "next-generation" Very Large Array (ngVLA) operating in the frequency range 1-116 GHz could be a groundbreaking tool to study the detailed astrophysics of the "matter cycle" in the Milky Way and other galaxies. If optimized for high brightness sensitivity, the ngVLA would bring detailed microwave spectroscopy and modeling of the full radio spectral energy distribution into regular use as survey tools at resolutions of 0.1- 1 arcseconds. This wavelength range includes powerful diagnostics of density, excitation, and chemistry in the cold ISM, as well as multiple tracers of the rate of recent star formation, the magnetic field, shocks, and properties of the ionized ISM. We highlight design considerations that would make this facility revolutionary in this area, the foremost of which is a large amount of collecting area on ~km-length baselines. We also emphasize the strong case for harnessing the large proposed collecting area of the ngVLA for very long baseline applications as part of the core design. This would allow measurements of the three dimensional space motions of galaxies to beyond the Local Group and mapping of the Milky Way out to the far side of the disk. Finally, we discuss the gains from the proposed combination of very high resolution and sensitivity to thermal emission, which include observing the feeding of black holes and resolving forming protoclusters.
We report here on key science topics for the Next Generation Very Large Array in the areas of time domain, fundamental physics, and cosmology. Key science cases considered are pulsars in orbit around the Galactic Center massive black hole, Sagittarius A*, electromagnetic counterparts to gravitational waves, and astrometric cosmology. These areas all have the potential for ground-breaking and transformative discovery. Numerous other topics were discussed during the preparation of this report and some of those discussions are summarized here, as well. There is no doubt that further investigation of the science case will reveal rich and compelling opportunities.
We present multiwavelength, multi-telescope, ground-based follow-up photometry of the white dwarf WD 1145+017, that has recently been suggested to be orbited by up to six or more, short-period, low-mass, disintegrating planetesimals. We detect 9 significant dips in flux of between 10% and 30% of the stellar flux from our ground-based photometry. We observe transits deeper than 10% on average every ~3.6 hr in our photometry. This suggests that WD 1145+017 is indeed being orbited by multiple, short-period objects. Through fits to the multiple asymmetric transits that we observe, we confirm that the transit egress timescale is usually longer than the ingress timescale, and that the transit duration is longer than expected for a solid body at these short periods, all suggesting that these objects have cometary tails streaming behind them. The precise orbital periods of the planetesimals in this system are unclear from the transit-times, but at least one object, and likely more, have orbital periods of ~4.5 hours. We are otherwise unable to confirm the specific periods that have been reported, bringing into question the long-term stability of these periods. Our high precision photometry also displays low amplitude variations suggesting that dusty material is consistently passing in front of the white dwarf, either from discarded material from these disintegrating planetesimals or from the detected dusty debris disk. For the significant transits we observe, we compare the transit depths in the V- and R-bands of our multiwavelength photometry, and find no significant difference; therefore, for likely compositions the radius of single-size particles in the cometary tails streaming behind the planetesimals in this system must be ~0.15 microns or larger, or ~0.06 microns or smaller, with 2-sigma confidence.
We present some new accurate CCD photometry analysis of the white light solar corona at the time of the last 20 March 2015 total eclipse (airborne observations on a Falcon 7X and at ground-based Svalbard). We measured coronal brightness profiles taken along radial directions from 1.001 to 3 solar radii in the northern, southern and equatorial regions, after removing the F corona and the sky background. These studies allow to evaluate the density gradients, structures and temperature heterogeneity, by considering the Thomson scattering in white light of the K corona and also emissions of the EUV Fe XII 193A (1 to 2 MK) and Fe XI 171/174 (lower temperature) simultaneously observed by SDO/AIA and SWAP Proba2 space missions. Some dispersion between the regions is noticed. The limitation of the hydrostatic equilibrium assumption in the solar atmosphere is discussed as well as the contribution of the magnetic field pressure gradients as illustrated by a comparison with the model stationary magnetic corona from Predictive Sc. Inc. These results are compared with the results of the quieter 2010 total solar eclipse corona analyzed with the same method. This photometric analysis of the inner and intermediate white light corona will contribute to the preparation of the Aspiics/Proba 3 flying formation future coronagraphic mission of ESA for new investigation at time of artificial eclipses produced in Space. Note that Aspiics will also observe in the He I D3 line at 5876 A, and will record intensities of the Fe XIV line 5303A simultaneously with the analysis of the orange white- light continuum, including precise polarimetry analysis.
We summarize the design, capabilities, and some of the priority science goals of a next generation Very Large Array (ngVLA). The ngVLA is an interferometric array with 10x larger effective collecting area and 10x higher spatial resolution than the current VLA and the Atacama Large Millimeter Array (ALMA), optimized for operation in the wavelength range 0.3cm to 3cm. The ngVLA opens a new window on the Universe through ultra-sensitive imaging of thermal line and continuum emission down to milliarcecond resolution, as well as unprecedented broad band continuum polarimetric imaging of non-thermal processes. The continuum resolution will reach 9mas at 1cm, with a brightness temperature sensitivity of 6K in 1 hour. For spectral lines, the array at 1" resolution will reach 0.3K surface brightness sensitivity at 1cm and 10 km/s spectral resolution in 1 hour. These capabilities are the only means with which to answer a broad range of critical scientific questions in modern astronomy, including direct imaging of planet formation in the terrestrial-zone, studies of dust-obscured star formation and the cosmic baryon cycle down to pc-scales out to the Virgo cluster, making a cosmic census of the molecular gas which fuels star formation back to first light and cosmic reionization, and novel techniques for exploring temporal phenomena from milliseconds to years. The ngVLA is optimized for observations at wavelengths between the superb performance of ALMA at submm wavelengths, and the future SKA1 at few centimeter and longer wavelengths. This memo introduces the project. The science capabilities are outlined in a parallel series of white papers. We emphasize that this initial set of science goals are simply a starting point for the project. We invite comment on these programs, as well as new ideas, through our public forum link on the ngVLA web page https://science.nrao.edu/futures/ngvla
We use microwave temperature maps from two seasons of data from the Atacama Cosmology Telescope (ACTPol) at 146 GHz, together with the Constant Mass CMASS galaxy sample from the Baryon Oscillation Spectroscopic Survey to measure the kinematic Sunyaev-Ze\v{l}dovich (kSZ) effect over the redshift range z = 0.4 - 0.7. We use galaxy positions and the continuity equation to obtain a reconstruction of the line-of-sight velocity field. We stack the cosmic microwave background temperature at the location of each halo, weighted by the corresponding reconstructed velocity. The resulting best fit kSZ model is preferred over the no-kSZ hypothesis at 3.3sigma and 2.9sigma for two independent velocity reconstruction methods, using 25,537 galaxies over 660 square degrees. The effect of foregrounds that are uncorrelated with the galaxy velocities is expected to be well below our signal, and residual thermal Sunyaev-Ze\v{l}dovich contamination is controlled by masking the most massive clusters. Finally, we discuss the systematics involved in converting our measurement of the kSZ amplitude into the mean free electron fraction of the halos in our sample.
This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.
We present the history of investigation of the dynamical properties of pairs and groups of genetically related long-period comets (other than the Kreutz sungrazing system). Members of a comet pair or group move in nearly identical orbits and their origin as fragments of a common parent comet is unquestionable. The only variable is the time of perihelion passage, which differs from member to member considerably due primarily to an orbital-momentum increment acquired during breakup. Meter-per-second separation velocities account for gaps of years or tens of years, thanks to the orbital periods of many millennia. The physical properties of individual members may not at all be alike, as illustrated by the trio of C/1988 A1, C/1996 Q1, and C/2015 F3. We exploit orbital similarity to examine whether the celebrated and as yet unidentified object, discovered from the Lick Observatory near the Sun at sunset on 1921 August 7, happened to be a member of such a pair and to track down the long-period comet to which it could be genetically related. Our search shows that the Lick object, which could not be a Kreutz sungrazer, was most probably a companion to comet C/1847 C1 (Hind), whose perihelion distance was ~9 R_sun and true orbital period approximately 8300 years. The gap of 74.4 years between their perihelion times is consistent with a separation velocity of ~1 m/s that set the fragments apart following the parent's breakup in a general proximity of perihelion during the previous return to the Sun in the 7th millennium BCE.
We present the analysis of the entire HARPS observations of three stars that host planetary systems: HD1461, HD40307, and HD204313. The data set spans eight years and contains more than 200 nightly averaged velocity measurements for each star. This means that it is sensitive to both long-period and low-mass planets and also to the effects induced by stellar activity cycles. We modelled the data using Keplerian functions that correspond to planetary candidates and included the short- and long-term effects of magnetic activity. A Bayesian approach was taken both for the data modelling, which allowed us to include information from activity proxies such as $\log{(R'_{\rm HK})}$ in the velocity modelling, and for the model selection, which permitted determining the number of significant signals in the system. The Bayesian model comparison overcomes the limitations inherent to the traditional periodogram analysis. We report an additional super-Earth planet in the HD1461 system. Four out of the six planets previously reported for HD40307 are confirmed and characterised. We discuss the remaining two proposed signals. In particular, we show that when the systematic uncertainty associated with the techniques for estimating model probabilities are taken into account, the current data are not conclusive concerning the existence of the habitable-zone candidate HD40307 g. We also fully characterise the Neptune-mass planet that orbits HD204313 in 34.9 days.
High-resolution spectral profiles of Na I D lines from the interstellar medium towards 64 stars in the direction of the Vela supernova remnant are presented. This survey conducted mostly between 2011-12 complements an earlier survey of the same stars by Cha & Sembach done in the 1993-96 period. The interval of 15 to 18 years provides a base line to search for changes in the interstellar profiles. Dramatic disappearance of strong absorption components at low radial velocity is seen towards three stars - HD 63578, HD 68217, HD 76161 - over 15-18 years; HD 68217 and HD 76161 are associated with the Vela SNR but HD 63578 is likely associated with the wind bubble of g2 Velorum. The vanishing of these cold neutral clouds in the short time of 15 to 18 years needs some explanation. Other changes are seen in high-velocity Na D components.
Fast outflows of gas, driven by the interaction between the radio-jets and ISM of the host galaxy, are being observed in an increasing number of galaxies. One such example is the nearby radio galaxy 3C293. In this paper we present Integral Field Unit (IFU) observations taken with OASIS on the William Herschel Telescope (WHT), enabling us to map the spatial extent of the ionised gas outflows across the central regions of the galaxy. The jet-driven outflow in 3C293 is detected along the inner radio lobes with a mass outflow rate ranging from $\sim 0.05-0.17$ solar masses/yr (in ionised gas) and corresponding kinetic power of $\sim 0.5-3.5\times 10^{40}$ erg/s. Investigating the kinematics of the gas surrounding the radio jets (i.e. not directly associated with the outflow), we find line-widths broader than $300$ km/s up to 5 kpc in the radial direction from the nucleus (corresponding to 3.5 kpc in the direction perpendicular to the radio axis at maximum extent). Along the axis of the radio jet line-widths $>400$ km/s are detected out to 7 kpc from the nucleus and line-widths of $>500$ km/s at a distance of 12 kpc from the nucleus, indicating that the disturbed kinematics clearly extend well beyond the high surface brightness radio structures of the jets. This is suggestive of the cocoon structure seen in simulations of jet-ISM interaction and implies that the radio jets are capable of disturbing the gas throughout the central regions of the host galaxy in all directions.
For reliable event reconstruction of Imaging Atmospheric Cherenkov Telescopes (IACTs), calibration of the optical throughput efficiency is required. Within current facilities, this is achieved through the use of ring shaped images generated by muons. Here, a complementary approach is explored, achieving cross calibration of elements of IACT arrays through pairwise comparisons between telescopes, focussing on its applicability to the upcoming Cherenkov Telescope Array (CTA). Intercalibration of telescopes of a particular type using eventwise comparisons of shower image amplitudes has previously been demonstrated to recover the relative telescope optical responses. A method utilising the reconstructed energy as an alternative to image amplitude is presented, enabling cross calibration between telescopes of varying types within an IACT array. Monte Carlo studies for two plausible CTA layouts have shown that this calibration procedure recovers the relative telescope response efficiencies at the few percent level.
We present the results of the quest for ring-type structures on the maps observed by the Planck satellite.
We present simulations of star forming filaments incorporating - to our knowledge - the largest chemical network used to date on-the-fly in a 3D-MHD simulation. The network contains 37 chemical species and about 300 selected reaction rates. For this we use the newly developed package KROME (Grassi et al. 2014). We combine the KROME package with an algorithm which allows us to calculate the column density and attenuation of the interstellar radiation field necessary to properly model heating and ionisation rates. Our results demonstrate the feasibility of using such a complex chemical network in 3D-MHD simulations on modern supercomputers. We perform simulations with different strengths of the interstellar radiation field and the cosmic ray ionisation rate. We find that towards the centre of the filaments there is gradual conversion of hydrogen from H^+ over H to H_2 as well as of C^+ over C to CO. Moreover, we find a decrease of the dust temperature towards the centre of the filaments in agreement with recent HERSCHEL observations.
The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order $10^{-12}$erg cm$^{-2}$ s$^{-1}$ deg$^{-2}$. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-phase interstellar medium including thermal feedback by supernova explosions as well as cosmic-ray feedback. We include a time-dependent chemical network, self-shielding by gas and dust, and self-gravity. Pure thermal feedback alone is sufficient to produce the observed surface brightness, although it is very sensitive to the supernova rate. Cosmic rays suppress this sensitivity and reduce the surface brightness because they drive cooler outflows. Self-gravity has by far the largest effect because it accumulates the diffuse gas in the disk in dense clumps and filaments, so that supernovae exploding in voids can eject a large amount of hot gas into the halo. This can boost the surface brightness by several orders of magnitude. Although our simulations do not reach a steady state, all simulations produce surface brightness values of the same order of magnitude as the observations, with the exact value depending sensitively on the simulation parameters. We conclude that star formation feedback alone is sufficient to explain the origin of the hot halo gas, but measurements of the surface brightness alone do not provide useful diagnostics for the study of galactic star formation.
In this paper we investigate the opportunities provided by the James Webb Space Telescope (JWST) for significant scientific advances in the study of solar system bodies and rings using stellar occultations. The strengths and weaknesses of the stellar occultation technique are evaluated in light of JWST's unique capabilities. We identify several possible JWST occultation events by minor bodies and rings, and evaluate their potential scientific value. These predictions depend critically on accurate a priori knowledge of the orbit of JWST near the Sun-Earth Lagrange-point 2 (L2). We also explore the possibility of serendipitous stellar occultations by very small minor bodies as a by-product of other JWST observing programs. Finally, to optimize the potential scientific return of stellar occultation observations, we identify several characteristics of JWST's orbit and instrumentation that should be taken into account during JWST's development.
A gas-grain time dependent chemical code, UCL\_CHEM, has been used to investigate the possibility of using chemical tracers to differentiate between the possible formation mechanisms of brown dwarfs. In this work, we model the formation of a pre-brown dwarf core through turbulent fragmentation by following the depth-dependent chemistry in a molecular cloud through the step change in density associated with an isothermal shock and the subsequent freefall collapse once a bound core is produced. Trends in the fractional abundance of molecules commonly observed in star forming cores are then explored to find a diagnostic for identifying brown dwarf mass cores formed through turbulence. We find that the cores produced by our models would be bright in CO and NH$_3$ but not in HCO$^+$. This differentiates them from models using purely freefall collapse as such models produce cores that would have detectable transitions from all three molecules.
Aims. We present the first orbital solution for the O-type supergiant star HD 74194, which is the optical counterpart of the supergiant fast X-ray transient IGR J08408-4503. Methods. We measured the radial velocities in the optical spectrum of HD 74194, and we determined the orbital solution for the first time. We also analysed the complex H{\alpha} profile. Results. HD 74194 is a binary system composed of an O-type supergiant and a compact object in a short-period ($P=9.5436\pm0.0002$ d) and high-eccentricity ($e=0.63\pm0.03$) orbit. The equivalent width of the H{\alpha} line is not modulated entirely with the orbital period, but seems to vary in a superorbital period ($P=285\pm10$ d) nearly 30 times longer than the orbital one.
This is a follow-up study of a work by Kramers et al. (2013) on an unusual diamond-rich rock found in the SW side of the Libyan Desert Glass strewn field. This pebble, called Hypatia, is composed of almost pure carbon. Transmission Electron Microscopy and X-ray diffraction results reveal that Hypatia is made of defect-rich diamond containing lonsdaleite and deformation bands. These characteristics are compatible with an impact origin on Earth and/or in space. We analyzed concentrations and isotopic compositions of all five noble gases and nitrogen in several mg sized Hypatia samples. These data confirm that Hypatia is extra-terrestrial. The sample is rich in trapped noble gases with an isotopic composition close to the meteoritic Q component. 40Ar/36Ar ratios in individual steps are as low as 0.4. Concentrations of cosmic-ray produced 21Ne correspond to a nominal cosmic-ray exposure age of ca. 0.1 Myr if produced in a typical m-sized meteoroid. Such an atypically low nominal exposure age suggests high shielding in a considerably larger body. In addition to the Xe-Q composition, an excess of radiogenic 129Xe (from the decay of extinct 129I) is observed (129Xe/132Xe = 1.18 +/- 0.03). Two N components are present, an isotopically heavy component ({\delta}15N = +20 permil) released at low temp. and a major light component ({\delta}15N = -110 permil) at higher temp. This disequilibrium in N suggests that the diamonds in Hypatia were formed in space. Our data are broadly consistent with concentrations and isotopic compositions of noble gases in at least three different types of carbon-rich meteoritic materials. However, Hypatia does not seem to be related to any of these materials, but may have sampled a similar cosmochemical reservoir. Our study does not confirm the presence of exotic noble gases that led Kramers et al. to propose that Hypatia is a remnant of a comet that impacted the Earth.
We present the intensive spectroscopic follow up of the type Ia supernova (SN Ia) 2014J in the starburst galaxy M82. Twenty-seven optical ISIS/WHT, ACAM/WHT and IDS/INT spectra have been acquired from January 22nd to September 1st 2014. After correcting the observations for the recession velocity of M82 and for Milky Way and host galaxy extinction, we measured expansion velocities from spectral line blueshifts and pseudo-equivalent width of the strongest features in the spectra, which gives an idea on how elements are distributed with velocity within the ejecta. We position SN 2014J in the Benetti (2005), Branch et al. (2006) and Wang et al. (2009) diagrams. These diagrams are based on properties of the Si II features and provide dynamical and chemical information about the SN ejecta. The nearby low-extinguished SN 2011fe is shown as a reference for comparisons. SN 2014J is a border-line object between the Core-normal (CN) and Broad-line (BL) groups, which corresponds to an intermediate position between Low Velocity Gradient (LVG) and High Velocity Gradient (HVG) objects. SN 2014J follows the R(Si II)-Dm15 correlation, which confirms its classification as a relatively normal SN Ia. Our description of the SN Ia in terms of the evolution of the pseudo-equivalent width of various ions as well as the position in the various diagrams put this specific SN Ia into the overall sample of SN Ia.
Pluto has been observed by the New Horizons space probe to have some relatively fresh ice on the old ices covering most of the surface. Pluto was thought to consist of only a rocky core below the ice. Here I show that Pluto can have an iron core, as can also its companion Charon, which has recently been modelled to have one. The presence of an iron core means the giant impact origin calculations should be redone to include iron and thus higher temperatures. An iron core leads to the possibility of a different geology. An originally molten core becomes solid later, with contraction and a release of latent heat. The space vacated allows the upper rock layers to flow downwards at some locations at the surface of the core, and some of the ice above the rock to descend, filling the spaces left by the rock motion downwards. These phenomena can lead to the forces recently deforming the icy surface of Pluto, and in a lesser way, of Charon.
We report photon-noise limited performance of horn-coupled, aluminum lumped-element kinetic inductance detectors at millimeter wavelengths. The detectors are illuminated by a millimeter-wave source that uses an active multiplier chain to produce radiation between 140 and 160 GHz. We feed the multiplier with either amplified broadband noise or a continuous-wave tone from a microwave signal generator. We demonstrate that the detector response over a 40 dB range of source power is well-described by a simple model that considers the number of quasiparticles. The detector noise-equivalent power (NEP) is dominated by photon noise when the absorbed power is greater than approximately 1 pW, which corresponds to $\mathrm{NEP} \approx 2 \times 10^{-17} \; \mathrm{W} \; \mathrm{Hz}^{-1/2}$, referenced to absorbed power. At higher source power levels we observe the relationships between noise and power expected from the photon statistics of the source signal: $\mathrm{NEP} \propto P$ for broadband (chaotic) illumination and $\mathrm{NEP} \propto P^{1/2}$ for continuous-wave (coherent) illumination. We develop a detailed model for the device noise and demonstrate absolute calibration of the absorbed power in both source modes using the scaling of the photon noise with power.
We monitored the recent FUor 2MASS J06593158-0405277 (V960 Mon) since November 2009 at various observatories and multiple wavelengths. After the outburst by nearly 2.9 mag in $r$ around September 2014 the brightness gently fades until April 2015 by nearly 1 mag in $U$ and 0.5 mag in $z$. Thereafter the brightness at $\lambda>5000 \AA$ was constant until June 2015 while the shortest wavelengths ($U, B$) indicate a new rise, similar to that seen for the FUor V2493 Cyg (HBC722). Our near-infrared (NIR) monitoring between December 2014 and April 2015 shows a smaller outburst amplitude ($\sim$2 mag) and a smaller (0.2 $-$ 0.3 mag) post-outburst brightness decline. Optical and NIR color-magnitude diagrams indicate that the brightness decline is caused by growing extinction. The post-outburst light curves are modulated by an oscillating color-neutral pattern with a period of about 17 days and an amplitude declining from $\sim$0.08 mag in October 2014 to $\sim$0.04 mag in May 2015. The properties of the oscillating pattern lead us to suggest the presence of a close binary with eccentric orbit.
To date O2 has definitely been detected in only two sources, namely rho Oph A and Orion, reflecting the extremely low abundance of O2 in the interstellar medium. One of the sources in the HOP program is the +50 km/s Cloud in the Sgr A Complex in the centre of the Milky Way. The Herschel HIFI is used to search for the 487 and 774 GHz emission lines of O2. No O2 emission is detected towards the Sgr A +50 km/s Cloud, but a number of strong emission lines of methanol (CH3OH) and absorption lines of chloronium (H2Cl+) are observed. A 3 sigma upper limit for the fractional abundance ratio of (O2)/(H2) in the Sgr A +50 km/s Cloud is found to be X(O2) less than 5 x 10(-8). However, since we can find no other realistic molecular candidate than O2 itself, we very tentatively suggest that two weak absorption lines at 487.261 and 487.302 GHz may be caused by the 487 GHz line of O2 in two foreground spiral arm clouds. By considering that the absorption may only be apparent, the estimated upper limit to the O2 abundance of less than (10-20) x 10(-6) in these foreground clouds is very high. This abundance limit was determined also using Odin non-detection limits. If the absorption is due to a differential Herschel OFF-ON emission, the O2 fractional abundance may be of the order of (5-10) x 10(-6). With the assumption of pure absorption by foreground clouds, the unreasonably high abundance of (1.4-2.8) x 10(-4) was obtained. The rotation temperatures for CH3OH-A and CH3OH-E lines in the +50 km/s Cloud are found to be 64 and 79 K, respectively, and the fractional abundance of CH3OH is approximately 5 x 10(-7).
Most successful galaxy formation scenarios now postulate that the intense
star formation in massive, high-redshift galaxies during their major growth
period was truncated when powerful AGNs launched galaxy-wide outflows of gas
that removed large parts of the interstellar medium. The most powerful radio
galaxies at z~2 show clear signatures of such winds, but are too rare to be
good representatives of a generic phase in the evolution of all massive
galaxies at high redshift. Here we present SINFONI imaging spectroscopy of 12
radio galaxies at z~2 that are intermediate between the most powerful radio and
vigorous starburst galaxies in radio power, and common enough to represent a
generic phase in the early evolution of massive galaxies.
The kinematic properties are diverse, with regular velocity gradients with
amplitudes of Delta v=200-400 km s^-1 as in rotating disks as well as irregular
kinematics with multiple velocity jumps of a few 100 km s^-1. Line widths are
generally high, typically around FWHM=800 km s^-1, consistent with wind
velocities in hydrodynamic models. A broad H-alpha line in one target implies a
black hole mass of a few 10^9 M$_sun. The ratio of line widths, sigma, to bulk
velocity, v, is so large that even the gas in galaxies with regular velocity
fields is unlikely to be gravitationally bound. It is unclear, however, whether
the large line widths are due to turbulence or unresolved, local outflows as
are sometimes observed at low redshifts. Comparison of the kinetic energy with
the energy supply from the AGN through jet and radiation pressure suggests that
the radio source still plays a dominant role for feedback, consistent with
low-redshift radio-loud quasars.
Improving constraints on the abundance of basaltic asteroids in the main
asteroid belt is necessary for better understanding the thermal and collisional
environment in the early solar system, for more rigorously identifying the
genetic family for (4) Vesta, for determining the effectiveness of
Yarkovsky/YORP in dispersing asteroid families, and for better quantifying the
population of basaltic asteroids in the outer main belt (a greater than 2.5 AU)
that are likely unrelated to (4) Vesta.
NIR spectral observations in this work were obtained for the Vp-type
asteroids (2011) Veteraniya, (5875) Kuga, (8149) Ruff, (9147) Kourakuen, (9553)
Colas, (15237) 1988 RL6, (31414) Rotaryusa, and (32940) 1995 UW4 during August
and September 2014 utilizing the SpeX spectrograph at the NASA Infrared
Telescope Facility (IRTF), Mauna Kea, Hawaii. Spectral band parameter (band
centers, Band Area Ratios) and mineralogical analysis (pyroxene chemistry) for
each average asteroid NIR reflectance spectrum suggests a
howardite-eucrite-diogenite (HED) meteorite analog for each asteroid. (5875)
Kuga is most closely associated with the eucrite meteorites, (31414) Rotaryusa
is most closely associated with the diogenites, and the remaining other six
asteroids are most closely associated with the howardite meteorites. Along with
orbital locations in the inner main belt and in the vicinity of (4) Vesta, the
existing evidence suggests that these eight Vp-type asteroids are also likely
Vestoids.
Context. Several new techniques to determine the metallicity of M dwarfs with better precision have been developed over the last decades. However, most of these studies were based on empirical methods. In order to enable detailed abundance analysis, standard methods established for warmer solar-like stars, i.e. model-dependent methods using fitting of synthetic spectra, still need to be used. Aims. In this work we continue the reliability confirmation and development of metallicity determinations of M dwarfs using high- resolution infrared spectra. The reliability was confirmed though analysis of M dwarfs in four binary systems with FGK dwarf companions and by comparison with previous optical studies of the FGK dwarfs. Methods. The metallicity determination was based on spectra taken in the J band (1.1-1.4 {\mu}m) with the CRIRES spectrograph. In this part of the infrared, the density of stellar molecular lines is limited, reducing the amount of blends with atomic lines enabling an accurate continuum placement. Lines of several atomic species were used to determine the stellar metallicity. Results. All binaries show excellent agreement between the derived metallicity of the M dwarf and its binary companion. Our results are also in good agreement with values found in the literature. Furthermore, we propose an alternative way to determine the effective temperature of M dwarfs of spectral types later than M2 through synthetic spectral fitting of the FeH lines in our observed spectra. Conclusions. We have confirmed that a reliable metallicity determination of M dwarfs can be achieved using high-resolution infrared spectroscopy. We also note that metallicites obtained with photometric metallicity calibrations available for M dwarfs only partly agree with the results we obtain from high-resolution spectroscopy.
Galactic starburst clusters play a twin role in astrophysics, serving as laboratories for the study of stellar physics and also delineating the structure and recent star formation history of the Milky Way. In order to exploit these opportunities we have undertaken a multi-epoch spectroscopic survey of the red supergiant dominated young massive clusters thought to be present at both near and far ends of the Galactic Bar. Significant spectroscopic variability suggestive of radial pulsations was found for the yellow supergiant VdBH 222 #505. Follow-up photometric investigations revealed modulation with a period of ~23.325d; both timescale and pulsational profile are consistent with a Cepheid classification. As a consequence #505 may be recognised as one of the longest period Galactic cluster Cepheids identified to date and hence of considerable use in constraining the bright end of the period/luminosity relation at solar metallicities. In conjunction with extant photometry we infer a distance of ~6kpc for VdBH222 and an age of ~20Myr. This results in a moderate reduction in both integrated cluster mass (~2x10^4Msun) and the initial stellar masses of the evolved cluster members (~10Msun). As such, VdBH222 becomes an excellent test-bed for studying the properties of some of the lowest mass stars observed to undergo type-II supernovae. Moreover, the distance is in tension with a location of VdBH 222 at the far end of the Galactic Bar. Instead a birthsite in the near 3kpc arm is suggested; providing compelling evidence of extensive recent star formation in a region of the inner Milky Way which has hitherto been thought to be devoid of such activity.
Jupiter-family comet 15P/Finlay has been reportedly quiet in activity for over a century but has harbored two outbursts during its 2014/2015 perihelion passage. Here we present an analysis of these two outbursts using a set of cometary observations. The outbursts took place between 2014 Dec. 15.4--16.0 UT and 2015 Jan. 15.5--16.0 UT as constrained by ground-based and spacecraft observations. We find a characteristic ejection speed of $V_0=300$ to $650 \mathrm{m \cdot s^{-1}}$ for the ejecta of the first outburst and $V_0=550$ to $750 \mathrm{m \cdot s^{-1}}$ for that of the second outburst using a Monte Carlo dust model. The mass of the ejecta is calculated to be $M_\mathrm{d}=2$ to $3\times10^5 \mathrm{kg}$ for the first outburst and $M_\mathrm{d}=4$ to $5\times10^5 \mathrm{kg}$ for the second outburst, corresponds to less than $10^{-7}$ of the nucleus mass. The specific energy of the two outbursts is found to be $0.3$ to $2\times10^5 \mathrm{J \cdot kg^{-1}}$. We also revisit the long-standing puzzle of the non-detection of the hypothetical Finlayid meteor shower by performing a cued search using the 13-year data from the Canadian Meteor Orbit Radar, which does not reveal any positives. The Earth will pass the 2014/2015 outburst ejecta around 2021 Oct. 6 at 22 h UT to Oct. 7 at 1 h UT, with a chance for some significant meteor activity in the radio range, which may provide further clues to the Finlayid puzzle. A southerly radiant in the constellation of Ara will favor the observers in the southern tip of Africa.
Precision measurements of the large scale structure of the Universe require large numbers of high fidelity mock catalogs to accurately assess, and account for, the presence of systematic effects. We introduce and test a scheme for generating mock catalogs rapidly using suitably derated N-body simulations. Our aim is to reproduce the large scale structure and the gross properties of dark matter halos with high accuracy, while sacrificing the details of the internal structure of the halos. By adjusting global and local time-steps in an N-body code, we demonstrate that we recover halo masses to better than 2% and the power spectrum (both in real and redshift space, for k = 1h/Mpc) to better than 1%, while requiring a factor of 4 less CPU time. We also calibrate the redshift spacing of outputs required to generate simulated light cones. We find that outputs separated by every z = 0.05 allow us to interpolate particle positions and velocities to reproduce the real and redshift space power spectra to better than 1% (out to k = 1h/Mpc). We apply these ideas to generate a suite of simulations spanning a range of cosmologies, motivated by the Baryon Oscillation Spectroscopic Survey (BOSS) but broadly applicable to future large scale structure surveys including eBOSS and DESI. As an initial demonstration of the utility of such simulations, we calibrate the shift in the baryonic acoustic oscillation peak position as a function of galaxy bias with higher precision than has been possible so far. This paper also serves to document the simulations, which we make publicly available.
NIKA, the prototype of the NIKA2 camera, is an instrument operating at the IRAM 30m telescope that can observe the sky simultaneously at 150 and 260GHz. One of the main goals of NIKA is to measure the pressure distribution in galaxy clusters at high angular resolution using the Sunyaev-Zel'dovich (SZ) effect. Such observations have already proved to be an excellent probe of cluster pressure distributions even at high redshifts. However, an important fraction of clusters host submm and/or radio point sources that can significantly affect the reconstructed signal. Here we report <20arcsec angular resolution observations at 150 and 260GHz of the cluster MACSJ1424, which hosts both radio and submm point sources. We examine the morphological distribution of the SZ signal and compare it to other datasets. The NIKA data are combined with Herschel satellite data to study the SED of the submm point source contaminants. We then perform a joint reconstruction of the ICM electronic pressure and density by combining NIKA, Planck, XMM-Newton and Chandra data, focussing on the impact of the radio and submm sources on the reconstructed pressure profile. We find that the large-scale pressure distribution is unaffected by the point sources due to the resolved nature of the NIKA observations. The reconstructed pressure in the inner region is slightly higher when the contribution of point sources are removed. We show that it is not possible to set strong constraints on the central pressure distribution without removing accurately these contaminants. The comparison with Xray only data shows good agreement for the pressure, temperature and entropy profiles, all indicating that MACSJ1424 is a dynamically relaxed cool core system. The present observations illustrate the possibility of measuring these quantities with a relatively small integration time, even at high redshift and without Xray spectroscopy.
The Kepler-discovered Systems with Tightly-packed Inner Planets (STIPs), typically with several planets of Earth to super-Earth masses on well-aligned, sub-AU orbits may host the most common type of planets, including habitable planets, in the Galaxy. They pose a great challenge for planet formation theories, which fall into two broad classes: (1) formation further out followed by inward migration; (2) formation in situ, in the very inner regions of the protoplanetary disk. We review the pros and cons of these classes, before focusing on a new theory of sequential in situ formation from the inside-out via creation of successive gravitationally unstable rings fed from a continuous stream of small (~cm-m size) "pebbles," drifting inward via gas drag. Pebbles first collect at the pressure trap associated with the transition from a magnetorotational instability (MRI)-inactive ("dead zone") region to an inner, MRI-active zone. A pebble ring builds up that begins to dominate the local mass surface density of the disk and spawns a planet. The planet continues to grow, most likely by pebble accretion, until it becomes massive enough to isolate itself from the accretion flow via gap opening. This reduces the local gas density near the planet, leading to enhanced ionization and a retreat of the dead zone inner boundary. The process repeats with a new pebble ring gathering at the new pressure maximum associated with this boundary. We discuss the theory's predictions for planetary masses, relative mass scalings with orbital radius, and minimum orbital separations, and their comparison with observed systems. Finally, we discuss open questions, including potential causes of diversity of planetary system architectures, i.e., STIPs versus Solar System analogs.
The MEarth Project is a photometric survey systematically searching the smallest stars nearest to the Sun for transiting rocky planets. Since 2008, MEarth has taken approximately two million images of 1844 stars suspected to be mid-to-late M dwarfs. We have augmented this survey by taking nightly exposures of photometric standard stars and have utilized this data to photometrically calibrate the $MEarth$ system, identify photometric nights, and obtain an optical magnitude with $1.5\%$ precision for each M dwarf system. Each optical magnitude is an average over many years of data, and therefore should be largely immune to stellar variability and flaring. We combine this with trigonometric distance measurements, spectroscopic metallicity measurements, and 2MASS infrared magnitude measurements in order to derive a color-magnitude-metallicity relation across the mid-to-late M dwarf spectral sequence that can reproduce spectroscopic metallicity determinations to a precision of 0.1 dex. We release optical magnitudes and metallicity estimates for 1567 M dwarfs, many of which did not have an accurate determination of either prior to this work. For an additional 277 stars without a trigonometric parallax, we provide an estimate of the distance assuming solar neighborhood metallicity. We find that the median metallicity for a volume limited sample of stars within 20 parsecs of the Sun is [Fe/H] = $-0.03 \pm 0.008$, and that 29 / 565 of these stars have a metallicity of [Fe/H] = $-0.5$ or lower, similar to the low-metallicity distribution of nearby G-dwarfs. When combined with the results of ongoing and future planet surveys targeting these objects, the metallicity estimates presented here will be important in assessing the significance of any putative planet-metallicity correlation.
The absence of low energy supersymmetry in run I data at the LHC has pushed the nominal scale for supersymmetry beyond a TeV. While this is consistent with the discovery of the Higgs boson at \approx 125 GeV, simple models with scalar and gaugino mass universality are being pushed into corners of parameter space. Some possibilities within the constrained minimal supersymmetric extension of the standard model (with four parameters) are discussed along with a one parameter extension in which the Higgs soft masses are non-universal. Also discussed are 2-, 3-, and 4-parameter versions of pure gravity mediated models with a wino, Higgsino, or bino LSP respectively.
We generalize the problem of strongly interacting neutron matter by adding a periodic external modulation. This allows us to study from first principles a neutron system that is extended and inhomogeneous, with connections to the physics of both neutron-star crusts and neutron-rich nuclei. We carry out fully non-perturbative microscopic Quantum Monte Carlo calculations of the energy of neutron matter at different densities, as well as different strengths and periodicities of the external potential. In order to remove systematic errors, we examine finite-size effects and the impact of the wave function ansatz. We also make contact with energy-density functional theories of nuclei and disentangle isovector gradient contributions from bulk properties. Finally, we calculate the static density-density linear response function of neutron matter and compare it with the response of other physical systems.
We incorporate Milky Way dark matter halo profile uncertainties, as well as an accounting of diffuse gamma-ray emission uncertainties in dark matter annihilation models for the Galactic Center Extended gamma-ray excess (GCE) detected by the Fermi Gamma Ray Space Telescope. The range of particle annihilation rate and masses expand when including these unknowns. However, empirical determinations of the Milky Way halo's local density and density profile leave the signal region to be in considerable tension with dark matter annihilation searches from combined dwarf galaxy analyses. Extreme changes to the Milky Way halo, which may be possible in cases of extreme adiabatic contraction, must be adopted to escape these constraints in a dark matter annihilation model for the GCE. Dark matter annihilation models that produce the gamma-ray excess via differential mechanisms in the GCE and dwarfs may circumvent this tension.
The search for continuous gravitational waves in a wide parameter space at fixed computing cost is most efficiently done with semicoherent methods, e.g. StackSlide, due to the prohibitive computing cost of the fully coherent search strategies. Prix&Shaltev arXiv:1201.4321 have developed a semi-analytic method for finding \emph{optimal} StackSlide parameters at fixed computing cost under ideal data conditions, i.e. gap-less data and constant noise floor. In this work we consider more realistic conditions by allowing for gaps in the data and changes in noise level. We show how the sensitivity optimization can be decoupled from the data selection problem. To find optimal semicoherent search parameters we apply a numerical optimization using as example the semicoherent StackSlide search. We also describe three different data selection algorithms. Thus the outcome of the numerical optimization consists of the optimal search parameters and the selected dataset. We first test the numerical optimization procedure under ideal conditions and show that we can reproduce the results of the analytical method. Then we gradually relax the conditions on the data and find that a compact data selection algorithm yields higher sensitivity compared to a greedy data selection procedure.
Recent progress in theory, experiment and observation challenges the mean field models using the conventional Skyrme interaction, suggesting that the extension of the conventional Skyrme interaction is necessary. In this work, we construct three Skyrme interaction parameter sets, namely, eMSL07, eMSL08 and eMSL09, based on an extended Skyrme interaction which includes additional momentum and density dependent two-body forces to effectively simulate the momentum dependence of the three-body force. The three new interactions can well reproduce both the ground-state properties and isoscalar giant monopole resonance energy of finite nuclei, nicely conform to the current knowledge on the equation of state of asymmetric nuclear matter around and below saturation density $\rho_0$, eliminate the notorious unphysical instabilities of symmetric nuclear matter and pure neutron matter at densities up to about $7.5\rho_0$, and simultaneously support heavier neutron stars with mass larger than two times solar mass. The new family of the extended Skyrme interactions can thus provide a unified description for the properties of asymmetric nuclear matter from sub- to supra-saturation densities in large region of isospin values and are appropriate for the study of nuclear matter, finite nuclei and neutron stars.
Although QCD axion models are widely studied as solutions to the strong CP problem, they generically confront severe fine-tuning problems to guarantee the anomalous PQ symmetry. In this letter, we propose a simple QCD axion model without any fine-tunings. We introduce an extra dimension and a pair of extra quarks living on two branes separately, which is also charged under a bulk Abelian gauge symmetry. We assume a monopole condensation on our brane at an intermediate scale, which implies that the extra quarks develop the chiral symmetry breaking and the PQ symmetry is broken. In contrast to the original Kim's model, our model explains the origin of the PQ symmetry thanks to the extra dimension and avoids the cosmological domain wall problem because of the chiral symmetry breaking in the Abelian gauge theory.
We study the secular dynamics of extrasolar planetary systems by extending the Lagrange-Laplace theory to high order and by including the relativistic effects. We investigate the long-term evolution of the planetary eccentricities via normal form and we find an excellent agreement with direct numerical integrations. Finally we set up a simple analytic criterion that allows to evaluate the impact of the relativistic effects in the long-time evolution.
Mercury is the unique known planet that is situated in a 3:2 spin-orbit resonance nowadays. Observations and models converge to the same conclusion: the planet is presently deeply trapped in the resonance and situated at the Cassini state $1$, or very close to it. We investigate the complete non-linear stability of this equilibrium, with respect to several physical parameters, in the framework of Birkhoff normal form and Nekhoroshev stability theory. We use the same approach adopted for the 1:1 spin-orbit case with a peculiar attention to the role of Mercury's non negligible eccentricity. The selected parameters are the polar moment of inertia, the Mercury's inclination and eccentricity and the precession rates of the perihelion and node. Our study produces a bound to both the latitudinal and longitudinal librations (of 0.1 radians) for a long but finite time (greatly exceeding the age of the solar system). This is the so-called effective stability time. Our conclusion is that Mercury, placed inside the 3:2 spin-orbit resonance, occupies a very stable position in the space of these physical parameters, but not the most stable possible one.
The formalism to treat quantization and evolution of cosmological perturbations of multiple fluids is described. We first construct the Lagrangian for both the gravitational and matter parts, providing the necessary relevant variables and momenta leading to the quadratic Hamiltonian describing linear perturbations. The final Hamiltonian is obtained without assuming any equations of motions for the background variables. This general formalism is applied to the special case of two fluids, having in mind the usual radiation and matter mix which made most of our current Universe history. Quantization is achieved using an adiabatic expansion of the basis functions. This allows for an unambiguous definition of a vacuum state up to the given adiabatic order. Using this basis, we show that particle creation is well defined for a suitable choice of vacuum and canonical variables, so that the time evolution of the corresponding quantum fields is unitary. This provides constraints for setting initial conditions for an arbitrary number of fluids and background time evolution. We also show that the common choice of variables for quantization can lead to an ill-defined vacuum definition. Our formalism is not restricted to the case where the coupling between fields is small, but is only required to vary adiabatically with respect to the ultraviolet modes, thus paving the way to consistent descriptions of general models not restricted to single-field (or fluid).
We study a level crossing between the QCD axion and an axion-like particle, focusing on the recently found phenomenon, the axion roulette, where the axion-like particle runs along the potential, passing through many crests and troughs, until it gets trapped in one of the potential minima. We perform detailed numerical calculations to determine the parameter space where the axion roulette takes place, and as a result domain walls are likely formed. The domain wall network without cosmic strings is practically stable, and it is nothing but a cosmological disaster. In a certain case, one can make domain walls unstable and decay quickly by introducing an energy bias without spoiling the Peccei-Quinn solution to the strong CP problem.
We consider the construction of gauge theories of gravity, focussing in particular on the extension of local Poincar\'e invariance to include invariance under local changes of scale. We work exclusively in terms of finite transformations, which allow for a more transparent interpretation of such theories in terms of gauge fields in Minkowski spacetime. Our approach therefore differs from the usual geometrical description of locally scale-invariant Poincar\'e gauge theory (PGT) and Weyl gauge theory (WGT) in terms of Riemann--Cartan and Weyl--Cartan spacetimes, respectively. In particular, we reconsider the interpretation of the Einstein gauge and also the equations of motion of matter fields and test particles in these theories. Inspired by the observation that the PGT and WGT matter actions for the Dirac field and electromagnetic field have more general invariance properties than those imposed by construction, we go on to present a novel alternative to WGT by considering an `extended' form for the transformation law of the rotational gauge field under local dilations, which includes its `normal' transformation law in WGT as a special case. The resulting `extended' Weyl gauge theory (eWGT) has a number of interesting features that we describe in detail. In particular, we present a new scale-invariant gauge theory of gravity that accommodates ordinary matter and is defined by the most general parity-invariant eWGT Lagrangian that is at most quadratic in the eWGT field strengths, and we derive its field equations. We also consider the construction of PGTs that are invariant under local dilations assuming either the `normal' or `extended' transformation law for the rotational gauge field, but show that they are special cases of WGT and eWGT, respectively.
We calculate the properties of isospin-asymmetric nuclear matter based on chiral nucleon-nucleon (NN) and three-nucleon (3N) interactions. To this end, we develop an improved normal-ordering framework that allows to include general 3N interactions starting from a plane-wave partial-wave-decomposed form. We present results for the energy per particle for general isospin asymmetries based on a set of different Hamiltonians, study their saturation properties, the incompressibility, symmetry energy, and also provide an analytic parametrization for the energy per particle as a function of density and isospin asymmetry.
Links to: arXiv, form interface, find, astro-ph, recent, 1510, contact, help (Access key information)