The gamma-ray burst (GRB) jet powers the afterglow emission by shocking the surrounding medium, and radio afterglow can now be routinely observed to almost a year after the explosion. Long-duration GRBs are accompanied by supernovae (SNe) that typically contain much more energy than the GRB jet. Here we consider the fact that the SN blast wave will also produce its own afterglow, which will peak at much later time (since it is non-relativistic), when the SN blast wave transitions from a coasting phase to a decelerating Sedov-Taylor phase. We predict that this component will peak generally a few tens of years after the explosion and it will outshine the GRB powered afterglow well-before its peak emission. In the case of GRB 030329, where the external density is constrained by the $\sim 10$-year coverage of the radio GRB afterglow, the radio emission is predicted to start rising over the next decade and to continue to increase for the following decades up to a level of $\sim 0.5$ mJy. Detection of the SN-powered radio emission will greatly advance our knowledge of particle acceleration in $ \sim 0.1$c shocks.
An accurate knowledge of the dark matter distribution in the Milky Way is of crucial importance for galaxy formation studies and current searches for particle dark matter. In this paper we set new dynamical constraints on the Galactic dark matter profile by comparing the observed rotation curve, updated with a comprehensive compilation of kinematic tracers, with that inferred from a wide range of observation-based morphologies of the bulge, disc and gas. The generalised Navarro-Frenk-White (NFW) and Einasto dark matter profiles are fitted to the data in order to determine the favoured ranges of local density, slope and scale radius. For a representative baryonic model, we find a local dark matter density 0.420+0.021-0.018 (2 sigma) +- 0.025 GeV/cm^3 (0.420+0.019-0.021 (2 sigma) +- 0.026 GeV/cm^3) for NFW (Einasto), where the second error is an estimate of the systematic due to baryonic modelling. The main sources of uncertainty inside and outside the solar circle are baryonic modelling and rotation curve measurements, respectively. Astronomical observations over the coming years are expected to reduce uncertainties on both fronts.
Theoretical simulations and observations at different angular resolutions have shown that magnetic fields have a central role in massive star formation. Like in low-mass star formation, the magnetic field in massive young stellar objects can either be oriented along the outflow axis or randomly. Measuring the magnetic field at milliarcsecond resolution (10-100 au) around a substantial number of massive young stellar objects permits determining with a high statistical significance whether the direction of the magnetic field is correlated with the orientation of the outflow axis or not. In late 2012, we started a large VLBI campaign with the European VLBI Network to measure the linearly and circularly polarized emission of 6.7 GHz methanol masers around a sample of massive star-forming regions. This paper focuses on the first seven observed sources, G24.78+0.08, G25.65+1.05, G29.86-0.04, G35.03+0.35, G37.43+1.51, G174.20-0.08, and G213.70-12.6. For all these sources, molecular outflows have been detected in the past. We detected a total of 176 methanol masing cloudlets toward the seven massive star-forming regions, 19% of which show linearly polarized emission. The methanol masers around the massive young stellar object MM1 in G174.20-0.08 show neither linearly nor circularly polarized emission. The linear polarization vectors are well ordered in all the other massive young stellar objects. We measured significant Zeeman splitting toward both A1 and A2 in G24.78+0.08, and toward G29.86-0.04 and G213.70-12.6. By considering all the 19 massive young stellar objects reported in the literature for which both the orientation of the magnetic field at milliarcsecond resolution and the orientation of outflow axes are known, we find evidence that the magnetic field (on scales 10-100 au) is preferentially oriented along the outflow axes.
We study the interaction between the atmospheres of Venus-like, non-magnetized exoplanets orbiting an M-dwarf star, and the stellar wind using a multi-species Magnetohydrodynaic (MHD) model. We focus our investigation on the effect of enhanced stellar wind and enhanced EUV flux as the planetary distance from the star decreases. Our simulations reveal different topologies of the planetary space environment for sub- and super-Alfvenic stellar wind conditions, which could lead to dynamic energy deposition in to the atmosphere during the transition along the planetary orbit. We find that the stellar wind penetration for non-magnetized planets is very deep, up to a few hundreds of kilometers. We estimate a lower limit for the atmospheric mass-loss rate and find that it is insignificant over the lifetime of the planet. However, we predict that when accounting for atmospheric ion acceleration, a significant amount of the planetary atmosphere could be eroded over the course of a billion years.
We describe a new method of combining optical and infrared photometry to select Luminous Red Galaxies (LRGs) at redshifts $z > 0.6$. We explore this technique using a combination of optical photometry from CFHTLS and HST, infrared photometry from the WISE satellite, and spectroscopic or photometric redshifts from the DEEP2 Galaxy Redshift Survey or COSMOS. We present a variety of methods for testing the success of our selection, and present methods for optimization given a set of rest-frame color and redshift requirements. We have tested this selection in two different regions of the sky, the COSMOS and Extended Groth Strip (EGS) fields, to reduce the effect of cosmic/sample variance. We have used these methods to assemble large samples of LRGs for two different ancillary programs as a part of the SDSS-III/ BOSS spectroscopic survey. This technique is now being used to select $\sim$600,000 LRG targets for SDSS-IV/eBOSS, which began observations in Fall 2014, and will be adapted for the proposed DESI survey. We have found these methods can select high-redshift LRGs efficiently with minimal stellar contamination; this is extremely difficult to achieve with selections that rely on optical photometry alone.
This paper presents an overview of results obtained during the CAWSES II period on the short term variability of the Sun and how it affects the near Earth space environment. CAWSES II was planned to examine the behavior of the solar terrestrial system as the solar activity climbed to its maximum phase in solar cycle 24. After a deep minimum following cycle 23, the Sun climbed to a very weak maximum in terms of the sunspot number in cycle 24 (MiniMax24), so many of the results presented here refer to this weak activity in comparison with cycle 23. The short term variability that has immediate consequence to Earth and geospace manifests as solar eruptions from closed field regions and high speed streams from coronal holes. Both electromagnetic (flares) and mass emissions (coronal mass ejections, CMEs) are involved in solar eruptions, while coronal holes result in high speed streams that collide with slow wind forming the so called corotating interaction regions (CIRs). Fast CMEs affect Earth via leading shocks accelerating energetic particles and creating large geomagnetic storms. CIRs and their trailing high speed streams (HSSs), on the other hand, are responsible for recurrent small geomagnetic storms and extended (days) of auroral zone activity, respectively. The latter lead to the acceleration of relativistic magnetospheric killer electrons. One of the major consequences of the weak solar activity is the altered physical state of the heliosphere that has serious implications for the shock-driving and storm causing properties of CMEs. Finally, a discussion is presented on extreme space weather events prompted by the 2012 July 23 super storm event that occurred on the backside of the Sun. Many of these studies were enabled by the simultaneous availability of remote-sensing and in situ observations from multiple vantage points with respect to the Sun Earth line.
The cosmic microwave background (CMB) radiation provides a remarkable window onto the early universe, revealing its composition and structure. In these lectures we review and discuss the physics underlying the main features of the CMB.
PSR B1259$-$63/SS 2883 is a $\gamma$-ray binary system composed of a radio pulsar in a long (1236.7 days) and elliptical ($e\sim0.87$) orbit around a Be star. In its 2010 periastron passage, multiwavelength emission from radio to TeV was observed, and an unexpected GeV flare was detected by the Fermi Large Area Telescope (LAT). Here we present the results of the LAT monitoring of PSR B1259$-$63 during its most recent 2014 periastron passage. We confirm that the GeV flare is recurrent within the orbit. The comparison of the 2014 and 2010 periastron passages shows overall similarities of flare durations, average flux levels, and spectra. In contrast, the detailed time evolutions of the two flares present interesting differences. Indeed, the light curves of the two flares show both a different structure and peak energy flux ($9.6 \pm1.8 \times 10^{-10}$ erg cm$^{-2}$ s$^{-1}$ and $7.1 \pm1.3 \times 10^{-10}$ erg cm$^{-2}$ s$^{-1}$, respectively in 2010 and 2014). While the tail of the 2010 flare the flux decayed exponentially, in 2014 it persisted at a high level. The interpretation of these differences as well as of the flare themselves is subject of debate.
We specify the range to which perturbations penetrate a planetesimal system. Such perturbations can originate from massive planets or from encounters with other stars. The latter can have an origin in the star cluster in which the planetary system was born, or from random encounters once the planetary system has escaped its parental cluster. The probability of a random encounter, either in a star cluster or in the Galactic field depends on the local stellar density, the velocity dispersion and the time spend in that environment. By adopting order of magnitude estimates we argue that the majority of planetary systems born in open clusters will have a {\em Parking zone}, in which planetesimals are affected by encounters in their parental star cluster but remain unperturbed after the star has left the cluster. Objects found in this range of semi-major axis and eccentricity preserve the memory of the encounter that last affected their orbits, and they can therefore be used to reconstruct this encounter. Planetary systems born in a denser environment, such as in a globular cluster are unlikely to have a Parking zone. We further argue that some planetary systems may have a {\em Frozen zone}, in which orbits are not affected either by the more inner massive planets or by external influences. Objects discovered in this zone will have preserved information about their formation in their orbital parameters.
We discuss the detectability of gravitationally bounded pairs of gas-giant planets (which we call "binary planets") in extrasolar planetary systems that are formed through orbital instability followed by planet-planet dynamical tides during their close encounters, based on the results of N-body simulations by Ochiai, Nagasawa and Ida (Paper I). Paper I showed that the formation probability of a binary is as much as $\sim 10\%$ for three giant planet systems that undergo orbital instability, and after post-capture long-term tidal evolution, the typical binary separation is 3--5 times the sum of physical radii of the planets. The binary planets are stable during main sequence lifetime of solar-type stars, if the stellarcentric semimajor axis of the binary is larger than 0.3 AU. We show that detecting modulations of transit light curves is the most promising observational method to detect binary planets. Since the likely binary separations are comparable to the stellar diameter, the shape of the transit light curve is different from transit to transit, depending on the phase of the binary's orbit. The transit durations and depth for binary planet transits are generally longer and deeper than those for the single planet case. We point out that binary planets could exist among the known inflated gas giant planets or objects classified as false positive detections at orbital radii > 0.3 AU, propose a binary planet explanation for the CoRoT candidate SRc01 E2 1066, and show that binary planets are likely to be present in, and could be detected using Kepler-quality data.
Fermi satellite discovered that cosmological gamma-ray bursts (GRBs) are accompanied by long GeV flashes. In two GRBs, an optical counterpart of the GeV flash has been detected. Recent work suggests that the GeV+optical flash is emitted by the external blast wave from the explosion in a medium loaded with copious $e^\pm$ pairs. The full light curve of the flash is predicted by a first-principle radiative transfer simulation and can be tested against observations. Here we examine a sample of 7 bursts with best GeV+optical data and test the model. We find that the observed light curves are in agreement with the theoretical predictions and allow us to measure three parameters for each burst: the Lorentz factor of the explosion, its isotropic kinetic energy, and the external density. With one possible exception of GRB 090510 (which is the only short burst in the sample) the ambient medium is consistent with a wind from a Wolf-Rayet progenitor. The wind density parameter $A=\rho r^2$ varies in the sample around $10^{11}$g/cm. The initial Lorentz factor of the blast wave varies from 200 to 540 and correlates with the burst luminosity. Radiative efficiency of the prompt emission in the sample is between 0.1 and 0.8. For the two bursts with detected optical flash, GRB 120711A and GRB 130427A, we also estimate the magnetization of the external blast wave. Remarkably, the model reproduces the entire optical light curve of GRB 120711A (with its sharp peak, fast decay, plateau, and break) as well as the GeV data. The spectrum of GeV flashes is predicted to extend above 0.1 TeV, where they can be detected by ground-based Cherenkov telescopes.
We present a three epoch survey for transient and variables in the extended Chandra Deep Field South at 5.5 GHz with the Australia Telescope Compact Array. A region covering $\sim$0.3 deg$^{2}$ was observed on timescales of 2.5 months and 2.5 years and typical sensitivities 12.1 $-$ 17.1 $\mu$Jy beam$^{-1}$ (1$\sigma$) were achieved. This survey represents the deepest search for transient and variable radio sources at 5.5 GHz. In total 124 sources were detected above the 5.5$\sigma$ level. One highly variable radio source was found with $\Delta S > 50%$ implying a surface density of $\sim$3 deg$^{-2}$. A further three radio sources were found with lower levels of variability equating to a surface density of $\sim$13 deg$^{-2}$ above a detection threshold of 82.3 $\mu$Jy. All of the variable sources have inverted radio spectra (between 1.4 and 5.5 GHz) and are associated with active galactic nuclei. We conclude that these variables are young gigahertz peaked-spectrum sources with active and self-absorbed radio jets. We explore the variability completeness of this sample and conclude that the fairly low levels of variability would only be detectable in 3$-$25% of all sources within the field. No radio transients were detected in this survey and we place an upper limit on the surface density of transient events $ < 7.5$ deg$^{-2}$ above a detection threshold of 68.8~$\mu$Jy.
KaVA (KVN and VERA Array) is a new combined VLBI array with KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). First, we briefly review the imaging capabilities of KaVA array which actually achieves more than three times better dynamic range than that achieved by VERA alone. The KaVA images clearly show detailed structures of extended radio jets in AGNs. Next, we represent the key science program to be led by KaVA AGN sub working group. We will conduct the monitoring observations of Sgr A* and M87 because of the largeness of their central super-massive black hole angular sizes. The main science goals of the program are (i) testing magnetically-driven-jet paradigm by mapping velocity fields of the M87 jet, and (ii) obtaining tight constraints on physical properties of radio emitting region in Sgr A*.
We present the second part of our detailed analysis of the hot sdO and spectroscopic standard star BD+28 4211, in which we focus on the optical spectrum. In the first part of our study, we determined the abundances of some 11 metals detected in the atmosphere of BD+28 4211 using UV spectra of the star and corroborated the fundamental parameters estimated in past studies (Teff $\sim$ 82,000 K, log g $\sim$ 6.2, and solar N(He)/N(H)). In this work, we aim at rederiving these secured parameters on the sole basis of high-quality optical spectra. A first grid of non-LTE line-blanketed model atmospheres, including metals with the abundances derived from the UV spectrum, does not give satisfactory results when we apply a standard simultaneous fitting procedure to the observed H and He lines of our optical spectra. The line profiles are not finely reproduced and the resulting effective temperatures, in particular, are too low by $\sim$10,000 K. We next investigate the probable cause of this failure, that is, the importance of missing opacity sources on the atmospheric stratification. We compare line profiles computed from models with artificially boosted metallicities, from solar abundances to 15$\times$ these values. We find that the structural effects saturate for a metallicity of $\sim$10x solar, and use this to compute a second full grid of models and synthetic spectra. This metal-enriched grid allows us to achieve significantly improved spectral fits with models having the expected parameters. Our test case thus reveals that there is still a need for models with enhanced metallicity for better estimating the atmospheric parameters of objects such as hot subdwarfs and hot white dwarfs if only optical spectra are available.
One of the most important task of the Gamma-Ray Burst field is the classification of the bursts. Many researches have proven the existence of the third kind (intermediate duration) of GRBs in the BATSE data. Recent works have analyzed BeppoSax and Swift observations and can also identify three types of GRBs in the data sets. However, the class memberships are probabilistic we have enough observed redshifts to calculate the redshift and spatial distribution of the intermediate GRBs. They are significantly farther than the short bursts and seems to be closer than the long ones.
This article represents a short review of the variability characteristics of young stellar objects. Variability is a key property of young stars. Two major origins may be distinguished: a scaled-up version of the magnetic activity seen on main-sequence stars and various processes related to circumstellar disks, accretion and outflows.
Methyl acetate (CH_3COOCH_3) has been recently observed by IRAM 30 m radio telescope in Orion though the presence of its deuterated isotopomers is yet to be confirmed. We therefore study the properties of various forms of methyl acetate, namely, CH_3COOCH_3, CH_2DCOOCH_3 and CH_3COOCH_2D. Our simulation reveals that these species could be produced efficiently both in gas as well as in ice phases. Production of methyl acetate could follow radical-radical reaction between acetyl (CH_3CO) and methoxy (CH_3O) radicals. To predict abundances of CH_3COOCH_3 along with its two singly deuterated isotopomers and its two isomers (ethyl formate and hydroxyacetone), we prepare a gas-grain chemical network to study chemical evolution of these molecules. Since gas phase rate coefficients for methyl acetate and its related species were unknown, either we consider similar rate coefficients for similar types of reactions (by following existing data bases) or we carry out quantum chemical calculations to estimate the unknown rate coefficients. For the surface reactions, we use adsorption energies of reactants from some earlier studies. Moreover, we perform quantum chemical calculations to obtain spectral properties of methyl acetate in infrared and sub-millimeter regions. We prepare two catalog files for the rotational transitions of CH_2DCOOCH_3 and CH_3COOCH_2D in JPL format, which could be useful for their detection in regions of interstellar media where CH_3COOCH_3 has already been observed.
It is long debated if pre-biotic molecules are indeed present in the interstellar medium. Despite substantial works pointing to their existence, pre-biotic molecules are yet to be discovered with a complete confidence. In this paper, our main aim is to study the chemical evolution of interstellar adenine under various circumstances. We prepare a large gas-grain chemical network by considering various pathways for the formation of adenine. Majumdar et al. (2012) proposed that in the absence of adenine detection, one could try to trace two precursors of adenine, namely, HCCN and NH_2CN. Recently Merz et al. (2014), proposed another route for the formation of adenine in interstellar condition. They proposed two more precursor molecules. But it was not verified by any accurate gas-grain chemical model. Neither was it known if the production rate would be high or low. Our paper fills this important gap. We include this new pathways to find that the contribution through this pathways for the formation of Adenine is the most dominant one in the context of interstellar medium. We propose that observers may look for the two precursors (C_3NH and HNCNH) in the interstellar media which are equally important for predicting abundances of adenine. We perform quantum chemical calculations to find out spectral properties of adenine and its two new precursor molecules in infrared, ultraviolet and sub-millimeter region. Our present study would be useful for predicting abundance of adenine.
$H_2$ is the most abundant interstellar species. Its deuterated forms ($HD$ and $D_2$) are also significantly abundant. Huge abundances of these molecules could be explained by considering the chemistry occurring on the interstellar dust. Because of its simplicity, Rate equation method is widely used to study the formation of grain-surface species. However, since recombination efficiency of formation of any surface species are heavily dependent on various physical and chemical parameters, Monte Carlo method would be best method suited to take care of randomness of the processes. We perform Monte Carlo simulation to study the formation of $H_2$, $HD$ and $D_2$ on interstellar ices. Adsorption energies of surface species are the key inputs for the formation of any species on interstellar dusts but binding energies of deuterated species are yet to known with certainty. A zero point energy correction exists between hydrogenated and deuterated species which should be considered while modeling the chemistry on the interstellar dusts. Following some earlier studies, we consider various sets of adsorption energies to study the formation of these species in diverse physical circumstances. As expected, noticeable difference in these two approaches (Rate equation method and Monte Carlo method) is observed for production of these simple molecules on interstellar ices. We introduce two factors, namely, $S_f$ and $\beta$ to explain these discrepancies: $S_f$ is a scaling factor, which could be used to correlate discrepancies between Rate equation and Monte Carlo methods. $\beta$ factor indicates the formation efficiency under various circumstances. Higher values of $\beta$ indicates a lower production efficiency. We found that $\beta$ increases with a decrease in rate of accretion from gas phase to grain phase.
We present a novel method for interpreting observations of high-mass X-ray binaries (HMXBs) based on a combination of spectroscopic data and numerical results from a radiation hydrodynamic model of stellar winds. We calculate synthetic Doppler tomograms of predicted emission in low/hard and high/soft X-ray states and compare them with Doppler tomograms produced using spectra of Cygnus X-1, a prototype of HMXBs. Emission from HMXBs is determined by local conditions within the circumstellar medium, namely density, temperature, and ionization state. These quantities depend strongly on the X-ray state of the systems. By increasing intensity of an X-ray emission produced by the compact companion in the HMXB model, we achieved a complete redistribution of the circumstellar medium in the vicinity of the modelled system. These changes (which simulate the transitions between two major spectral states) are also apparent in the synthetic Doppler tomograms which are in good agreement with the observations.
We apply our two-dimensional (2D), radially self-similar steady-state accretion flow model to the analysis of hydrodynamic simulation results of supercritical accretion flows. Self-similarity is checked and the input parameters for the model calculation, such as advective factor and heat capacity ratio, are obtained from time-averaged simulation data. Solutions of the model are then calculated and compared with the simulation results. We find that in the converged region of the simulation, excluding the part too close to the black hole, the radial distribution of azimuthal velocity $v_\phi$, density $\rho$ and pressure $p$ basically follows the self-similar assumptions, i.e. they are roughly proportional to $r^{-0.5}$, $r^{-n}$, and $r^{-(n+1)}$, respectively, where $n\sim0.85$ for the mass injection rate of $1000L_\mathrm{E}/c^2$, and $n\sim0.74$ for $3000L_\mathrm{E}/c^2$. The distribution of $v_r$ and $v_\theta$ agrees less with self-similarity, possibly due to convective motions in the $r\theta$ plane. The distribution of velocity, density and pressure in $\theta$ direction obtained by the steady model agrees well with the simulation results {within the calculation boundary of the steady model}. Outward mass flux in the simulations is overall directed toward polar angle of 0.8382 rad ($\sim 48.0^\circ$) for $1000L_\mathrm{E}/c^2$, and 0.7852 rad ($\sim 43.4^\circ$) for $3000L_\mathrm{E}/c^2$, and $\sim$94\% of the mass inflow are driven away as outflow, while outward momentum and energy fluxes are focused around the polar axis. Part of these fluxes lie in the region that are not calculated by the steady model, and special attention should be paid when the model is applied.
The locally observed cosmic ray spectrum has several puzzling features, such as the excess of positrons and antiprotons above $\sim 20$ GeV and the discrepancy in the slopes of the spectra of cosmic ray protons and heavier nuclei in the TeV-PeV energy range. We show that these features are consistently explained by a nearby source which was active $\sim 2$ Myr ago and has injected $(1-2)\times 10^{50}$ erg in cosmic rays. The transient nature of the source and its overall energy budget point to the supernova origin of this local cosmic ray source. The age of the supernova suggests that the local cosmic ray injection was produced by the same supernova that has deposited $^{60}$Fe isotopes in the deep ocean crust.
Aims. Forbush decrease (FD) is a transient decrease followed by a gradual recovery in the observed galactic cosmic ray intensity. We seek to understand the relationship between the FDs and near-Earth interplanetary magnetic field (IMF) enhancements associated with solar coronal mass ejections (CMEs). Methods. We use muon data at cutoff rigidities ranging from 14 to 24 GV from the GRAPES-3 tracking muon telescope to identify FD events. We select those FD events that have a reasonably clean profile, and magnitude > 0.25%. We use IMF data from ACE/WIND spacecrafts. We look for correlations between the FD profile and that of the one hour averaged IMF. We ask if the diffusion of high energy protons into the large scale magnetic field is the cause of the lag observed between the FD and the IMF. Results. The enhancement of the IMF associated with FDs occurs mainly in the shock-sheath region, and the turbulence level in the magnetic field is also enhanced in this region. The observed FD profiles look remarkably similar to the IMF enhancement profiles. The FDs typically lag the IMF enhancement by a few hours. The lag corresponds to the time taken by high energy protons to diffuse into the magnetic field enhancement via cross-field diffusion. Conclusions. Our findings show that high rigidity FDs associated with CMEs are caused primarily by the cumulative diffusion of protons across the magnetic field enhancement in the turbulent sheath region between the shock and the CME.
The expansion of Galactic HII regions can trigger the formation of a new generation of stars. However, little is know about the physical conditions that prevail in these regions. We study the physical conditions that prevail in specific zones towards expanding HII regions that trace representative media such as the photodissociation region, the ionized region, and condensations with and without ongoing star formation. We use the SPIRE Fourier Transform Spectrometer (FTS) on board $Herschel$ to observe the HII region RCW 120. Continuum and lines are observed in the $190-670\,\mu$m range. Line intensities and line ratios are obtained and used as physical diagnostics of the gas. We used the Meudon PDR code and the RADEX code to derive the gas density and the radiation field at nine distinct positions including the PDR surface and regions with and without star-formation activity. For the different regions we detect the atomic lines [NII] at $205\,\mu$m and [CI] at $370$ and $609\,\mu$m, the $^{12}{\rm CO}$ ladder between the $J=4$ and $J=13$ levels and the $^{13}{\rm CO}$ ladder between the $J=5$ and $J=14$ levels, as well as CH$ ^{+} $ in absorption. We find gas temperatures in the range $45-250\,$K for densities of $10^4-10^6\,{\rm cm}^{-3}$, and a high column density on the order of $N_{{\rm H}}\sim10^{22}\,{\rm cm}^{-2}$ that is in agreement with dust analysis. The ubiquitousness of the atomic and CH$ ^{+} $ emission suggests the presence of a low-density PDR throughout RCW 120. High-excitation lines of CO indicate the presence of irradiated dense structures or small dense clumps containing young stellar objects, while we also find a less dense medium ($N_{{\rm H}}\sim10^{20}\,{\rm cm}^{-2}$) with high temperatures ($80-200\,$K).
Aims: We investigate the electron acceleration in convective electric fields of cascading magnetic reconnection in a flaring solar corona and show the resulting hard X-ray (HXR) radiation spectra caused by Bremsstrahlung for the coronal source. Methods: We perform test particle calculation of electron motions in the framework of a guiding center approximation. The electromagnetic fields and their derivatives along electron trajectories are obtained by linearly interpolating the results of high-resolution adaptive mesh refinement (AMR) MHD simulations of cascading magnetic reconnection. Hard X-ray (HXR) spectra are calculated using an optically thin Bremsstrahlung model. Results: Magnetic gradients and curvatures in cascading reconnection current sheet accelerate electrons: trapped in magnetic islands, precipitating to the chromosphere and ejected into the interplanetary space. The final location of an electron is determined by its initial position, pitch angle and velocity. These initial conditions also influence electron acceleration efficiency. Most of electrons have enhanced perpendicular energy. Trapped electrons are considered to cause the observed bright spots along coronal mass ejection CME-trailing current sheets as well as the flare loop-top HXR emissions.
I present empirical measurements of the rate of relaxation in N-body simulations of stable spherical systems and distinguish two separate causes of relaxation: two-body effects, which cause energy exchange between particles of differing masses, and collective oscillations driven by shot noise that lead to an enhanced energy diffusion rate. I use four different methods to compute the gravitational field, and a 100-fold range in the numbers of particles in each case. I find the rate at which energy is exchanged between particles of differing masses does not depend at all on the force determination method, but I do find the energy diffusion rate, which is substantially enhanced by collective modes, is marginally lower when a field method is used. The relaxation rate in 3D is virtually independent of the method used because it is dominated by distant encounters; any method to estimate the gravitational field that correctly captures the contributions from distant particles must also capture their statistical fluctuations and the collective modes they drive.
We use the HB morphology of 48 Galactic GCs to study the radial distributions of the different stellar populations known to exist in globular clusters. Assuming that the (extremely) blue HB stars correspond to stars enriched in Helium and light elements, we compare the radial distributions of stars selected according to colour on the HB to trace the distribution of the secondary stellar populations in globular clusters. Unlike other cases, our data show that the populations are well mixed in 80% of the cases studied. This provides some constraints on the mechanisms proposed to pollute the interstellar medium in young globular clusters.
The statistical study of voids in the matter distribution promises to be an important tool for precision cosmology, but there are known discrepancies between theoretical models of voids and the voids actually found in large simulations or galaxy surveys. The empirical properties of observed voids are also not well understood. In this paper we study voids in an N-body simulation, using the ZOBOV watershed algorithm. As in other studies, we use sets of subsampled dark matter particles as tracers to identify voids, but we use the full-resolution simulation output to measure dark matter densities at the identified locations. Voids span a wide range of sizes and densities, but there is a clear trend towards larger voids containing deeper density minima, a trend which is expected for all watershed void finders. We also find that the tracer density at void locations is smaller than the true density, and that this relationship depends on the sampling density of tracers. We show that fitting functions given in the literature fail to match the density profiles of voids either quantitatively or qualitatively. The average enclosed density contrast within watershed voids varies widely with both the size of the void and the minimum density within it, but is always far from the shell-crossing threshold expected from theoretical models. Voids with deeper density minima also show much broader density profiles. We discuss the implications of these results for the excursion set approach to modelling such voids.
We report upper limits to the radio and X-ray emission from the newly discovered ultracool dwarf binary WISE J104915.57$-$531906.1 (Luhman 16AB). As the nearest ultracool dwarf binary (2 pc), its proximity offers a hefty advantage to studying plasma processes in ultracool dwarfs which are more similar in gross properties (radius, mass, temperature) to the solar system giant planets than stars. The radio and X-ray emission upper limits from the Australia Telescope Compact Array (ATCA) and Chandra observations, each spanning multiple rotation periods, provide the deepest fractional radio and X-ray luminosities to date on an ultracool dwarf, with $\log{(L_{\rm r,\nu}/L_{\rm bol}) [Hz^{-1}]} < -18.1$ (5.5 GHz), $\log{(L_{\rm r,\nu}/L_{\rm bol}) [Hz^{-1}]} < -17.9$ (9 GHz), and $\log{(L_{\rm x}/L_{\rm bol})} < -5.7$. While the radio upper limits alone do not allow for a constraint on the magnetic field strength, we limit the size of any coherently emitting region in our line of sight to less than 0.2\% of the radius of one of the brown dwarfs. Any source of incoherent emission must span less than about 20\% of the brown dwarf radius, assuming magnetic field strengths of a few tens to a few hundred Gauss. The fast rotation and large amplitude photometric variability exhibited by the T dwarf in the Luhman 16AB system are not accompanied by enhanced nonthermal radio emission, nor enhanced heating to coronal temperatures, as observed on some higher mass ultracool dwarfs, confirming the expected decoupling of matter and magnetic field in cool neutral atmospheres.
We study the contribution of thermal and non-thermal processes to the inverse Compton emission of the radio galaxy M 87 by modelling its broad-band emission. Through this we aim to derive insight into where within the AGN the X-ray, gamma-ray, and VHE emission is produced. We have analysed all available INTEGRAL IBIS/ISGRI data on M 87, spanning almost 10 years, to set an upper limit to the average hard X-ray flux of $f(20 - 60 \rm \, keV) < 3\times 10^{-12}$ $\rm \, erg \, cm^{-2} \, s^{-1}$, using several techniques beyond the standard analysis which are also presented here. We also analysed hard X-ray data from Suzaku/PIN taken late November 2006, and we report the first hard X-ray detection of M 87 with a flux of $f(20 - 60 \rm \, keV) = 10^{-11}\rm \, erg \, cm^{-2} \, s^{-1} $. In addition we analyse data from Fermi/LAT, INTEGRAL/JEM-X, and Suzaku/XIS. We collected historical radio/IR/optical and VHE data and combined them with the X-ray and gamma-ray data, to create broad-band spectral energy distributions for the average low-flux state and the flaring state. The resulting spectral energy distributions are modelled by applying a single-zone SSC model with a jet angle of theta = 15 degrees. We also show that modelling the core emission of M 87 using a single-zone synchrotron self-Compton model does represent the SED, suggesting that the core emission is dominated by a BL Lac type AGN core. Using SED modelling we also show that the hard X-ray emission detected in 2006 is likely due to a flare of the jet knot HST-1, rather than being related to the core.
The intensity mapping of the [CII] 157.7 $\rm \mu$m fine-structure emission line represents an ideal experiment to probe star formation activity in galaxies, especially in those that are too faint to be individually detected. Here, we investigate the feasibility of such an experiment for $z > 5$ galaxies. We construct the $L_{\rm CII} - M_{\rm h}$ relation from observations and simulations, then generate mock [CII] intensity maps by applying this relation to halo catalogs built from large scale N-body simulations. Maps of the extragalactic far-infrared (FIR) continuum, referred to as "foreground", and CO rotational transition lines and [CI] fine-structure lines referred to as "contamination", are produced as well. We find that, at 316 GHz (corresponding to $z_{\rm CII} = 5$), the mean intensities of the extragalactic FIR continuum, [CII] signal, all CO lines from $J=1$ to 13 and two [CI] lines are $\sim 3\times10^5$ Jy sr$^{-1}$, $\sim 1200$ Jy sr$^{-1}$, $\sim 800$ Jy sr$^{-1}$ and $\sim 100$ Jy sr$^{-1}$, respectively. We discuss a method that allows us to subtract the FIR continuum foreground by removing a spectrally smooth component from each line of sight, and to suppress the CO/[CI] contamination by discarding pixels that are bright in contamination emission. The $z > 5$ [CII] signal comes mainly from halos in the mass range $10^{11-12} \,M_\odot$; as this mass range is narrow, intensity mapping is an ideal experiment to investigate these early galaxies. In principle such signal is accessible to a ground-based telescope with a 6 m aperture, 150 K system temperature, a $128\times128$ pixels FIR camera in 5000 hr total integration time, however it is difficult to perform such an experiment by using currently available telescopes.
The young nearby solar-type star HD 91962 is a rare quadruple system where three companions revolve around the main component with periods of 170.3 days, 8.85 years, and 205 years. The two outer orbits are nearly co-planar, and all orbits have small eccentricities. We refine the visual orbit of the outer pair, determine the combined spectro-interferometric orbit of the middle 8.8-yr pair and the spectroscopic orbit of the inner binary. The middle and inner orbits are likely locked in a 1:19 resonance, the ratio of the outer and middle periods is ~23. The masses of all components are estimated (inside-out: 1.14, 0.32, 0.64, 0.64 solar mass), the dynamical parallax is 27.4+-0.6 mas. We speculate that this multiple system originated from collapse of an isolated core and that the companions migrated in a dissipative disk. Other multiple systems with similar features (coplanarity, small eccentricity, and period ratio around 20) are known.
We take advantage of the wealth of rotation measures data contained in the NVSS catalogue to derive new, statistically robust, upper limits on the strength of extragalactic magnetic fields. We simulate the extragalactic contribution to the rotation measures for a given field strength and correlation length, by assuming that the electron density follows the distribution of Lyman-$\alpha$ clouds. Based on the observation that rotation measures from low-luminosity distant radio sources do not exhibit any trend with redshift, while the extragalactic contribution instead grows with distance, we constrain fields with Mpc coherence length to be below 1.2 nG at the $2\sigma$ level, and fields coherent across the entire observable Universe below 0.5 nG. These limits do not depend on the particular origin of these cosmological fields.
In this work, early-time photometry of supernova (SN) 2014J is presented, extending the AAVSO CCD database to prediscovery dates. The applicability of NASA's small robotic MicroObservatory Network telescopes for photometric measurements is evaluated. Prediscovery and postdiscovery photometry of SN 2014J is measured from images taken by two different telescopes of the network, and is compared to measurements from the Katzman Automatic Imaging Telescope and the Itagaki Observatory. In the early light-curve phase (which exhibits stable spectral behavior with constant color indices), these data agree with reasonably high accuracy (better than 0.05 mag around maximum brightness, and 0.15 mag at earlier times). Owing to the changing spectral energy distribution of the SN and the different spectral characteristics of the systems used, differences increase after maximum light. We augment light curves of SN 2014J downloaded from the American Association of Variable Star Observers (AAVSO) online database with these data, and consider the complete brightness evolution of this important Type Ia SN. Furthermore, the first detection presented here (Jan. 15.427, 2014) appears to be one of the earliest observations of SN 2014J yet published, taken less than a day after the SN exploded.
We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons/second, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to: 1) construct kilometer scale phase screens necessary for wide-field telescopes, 2) reproduce atmospheric point-spread-function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, 3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and 4) recover system effective area predicted from analytic photometry integrals. This new code, the photon simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope (LSST) design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.
Self-interactions of dark matter particles can potentially lead to an observable separation between the dark matter halo and the stars of a galaxy moving through a region of large dark matter density. Such a separation has recently been observed in a galaxy falling into the core of the galaxy cluster Abell 3827. We estimated the DM self-interaction cross section needed to reproduce the observed effects and find that the sensitivity of Abell 3827 has been significantly overestimated in a previous study. Our corrected estimate is $\tilde{\sigma}/m_\text{DM} \sim 3\:\text{cm}^2\:\text{g}^{-1}$ when self-interactions result in an effective drag force and $\sigma/m_\text{DM} \sim 1.5\:\text{cm}^2\:\text{g}^{-1}$ for the case of contact interactions, in some tension with previous upper bounds.
The Atacama Large mm-submm Array (ALMA) is, by far, the most sensitive mm/submm telescope in the World. The ALMA Phasing Project (APP) will allow us to phase-up all the ALMA antennas and use them as one single VLBI station. This will be a key component of the Event Horizon Telescope (EHT), a Global VLBI array at millimeter wavelengths. A problem in the APP is the calibration and conversion of the polarization channels. Most VLBI stations record their signals in a circular basis, but the ALMA receivers record in a linear basis. The strategy that will be followed in the phased-ALMA VLBI observations will be to correlate in "mixed" basis (i.e., linear versus circular) and convert the visibilities to a pure circular basis after the correlation. We have developed an algorithm to perform such a polarization conversion of the VLBI visibilities. In these proceedings, we present the basics of the PolConvert algorithm and discuss on the polarization conversion in the general case were single dishes (besides phased arrays) record with linear receivers in VLBI observations. We show some results of PolConvert applied to realistic simulations, as well as a test with real VLBI observations at 86\,GHz between the Onsala radiotelescope (recording in linear basis) and the Effelsberg radiotelescope (recording in circular basis).
A recent article by Wojtak {\it et al} (arXiv:1504.00178) pointed out that the local gravitational redshift, despite its smallness ($\sim 10^{-5}$), can have a noticeable ($\sim 1\%$) systematic effect on our cosmological parameter measurements. The authors studied a few extended cosmological models (non-flat $\Lambda$CDM, $w$CDM, and $w_0$-$w_a$CDM) with a mock supernova dataset. We repeat this calculation and find that the $\sim 1\%$ biases are due to strong degeneracy between cosmological parameters. When Cosmic Microwave Background (CMB) data are added to break the degeneracy, the biases due to local gravitational redshift are negligible ($\lesssim 0.1 \sigma$).
Gravitational waves are perturbations in the spacetime that propagate at the speed of light. The study of such phenomenon is interesting because many cosmological processes and astrophysical objects, such as binary systems, are potential sources of gravitational radiation and can have their emissions detected in the near future by the next generation of interferometric detectors. Concerning the astrophysical objects, an interesting case is when there are several sources emitting in such a way that there is a superposition of signals, resulting in a smooth spectrum which spans a wide range of frequencies, the so-called stochastic background. In this paper, we are concerned with the stochastic backgrounds generated by compact binaries (i.e. binary systems formed by neutron stars and black holes) in the coalescing phase. In particular, we obtain such backgrounds by employing a new method developed in our previous studies.
Ionization of the Earth's atmosphere by sunlight forms a complex, multi-layered plasma environment within the Earth's magnetosphere, the innermost layers being the ionosphere and plasmasphere. The plasmasphere is believed to be embedded with cylindrical density structures (ducts) aligned along the Earth's magnetic field, but direct evidence for these remains scarce. Here we report the first direct wide-angle observation of an extensive array of field-aligned ducts bridging the upper ionosphere and inner plasmasphere, using a novel ground-based imaging technique. We establish their heights and motions by feature-tracking and parallax analysis. The structures are strikingly organized, appearing as regularly-spaced, alternating tubes of overdensities and underdensities strongly aligned with the Earth's magnetic field. These findings represent the first direct visual evidence for the existence of such structures.
We embed general $f(R)$ inflationary models in minimal supergravity plus matter, a single chiral superfield $\Phi$, with or without another superfield $S$, via a Jordan frame Einstein+scalar description. In particular, inflationary models like a generalized Starobinsky one are analyzed and constraints on them are found. We also embed the related models of conformal inflation, also described as Jordan frame Einstein+scalar models, in particular the conformal inflation from the Higgs model, and analyze the inflationary constraints on them.
We compare the general effective theory of one-body dark matter nucleon interactions to current direct detection experiments in a global multidimensional statistical analysis. We derive exclusion limits on the 28 isoscalar and isovector coupling constants of the theory, and show that current data place interesting constraints on dark matter-nucleon interaction operators usually neglected in this context. We characterize the interference patterns that can arise in dark matter direct detection from pairs of dark matter-nucleon interaction operators, or from isoscalar and isovector components of the same operator. We find that commonly neglected destructive interference effects weaken standard direct detection exclusion limits by up to one order of magnitude in the coupling constants.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
In current cosmological models, galaxies form from the gravitational collapse of small perturbations in the matter distribution. This process involves both a hierarchy of merging structures and smooth accretion, so that early galaxies are predicted to be morphologically irregular, clumpy, and compact. This is supported by recent observational data on samples of galaxies at redshift $z=8$ and beyond. The volumes accessible to these studies, both computational and observational are however thousands of times smaller than those that will be probed by upcoming telescopes, such as WFIRST. As a result, studies so far have never been able to reach the realm of massive galaxies. Whether among the myriad tiny proto-galaxies there exists a population with similarities to present day galaxies is an open question. Here we show, using BlueTides, the first hydrodynamic simulation large enough to resolve the relevant scales, that the first massive galaxies to form are in fact predicted to have extensive rotationally-supported disks and resemble in some ways Milky-way types seen at much lower redshifts. From a kinematic analysis of a statistical sample of 216 galaxies at redshift $z=8-10$ we have found that disk galaxies make up 70% of the population of galaxies with stellar mass $10^{10} M_\odot$ or greater. Cold Dark Matter cosmology therefore makes specific predictions for the population of large galaxies 500 million years after the Big Bang. We argue that wide-field satellite telescopes will in the near future discover these first massive disk galaxies. The simplicity of their structure and formation history should make possible new tests of cosmology.
We introduce the BlueTides simulation and report initial results for the luminosity functions of the first galaxies and AGN, and their contribution to reionization. BlueTides was run on the BlueWaters cluster at NCSA from $z=99$ to $z=8.0$ and includes 2$\times$7040$^3$ particles in a $400$Mpc/h per side box, making it the largest hydrodynamic simulation ever performed at high redshift. BlueTides includes a pressure-entropy formulation of smoothed particle hydrodynamics, gas cooling, star formation (including molecular hydrogen), black hole growth and models for stellar and AGN feedback processes. The star formation rate density in the simulation is a good match to current observational data at $z\sim 8-10$. We find good agreement between observations and the predicted galaxy luminosity function in the currently observable range $-18\le M_{\mathrm UV} \le -22.5$ with some dust extinction required to match the abundance of brighter objects. BlueTides implements a patchy reionization model that produces a fluctuating UV background. BlueTides predicts number counts for galaxies fainter than current observational limits which are consistent with extrapolating the faint end slope of the luminosity function with a power law index $\alpha\sim -1.8$ at $z\sim 8$ and redshift dependence of $\alpha\sim (1+z)^{-0.4}$. The AGN population has a luminosity function well fit by a power law with a slope $\alpha\sim -2.4$ that compares favourably with the deepest CANDELS-Goods fields. We investigate how these luminosity functions affect the progress of reionization, and find that a high Lyman-$\alpha$ escape fraction ($f_\mathrm{esc} \sim 0.5$) is required if galaxies dominate the ionising photon budget during reionization. Smaller galaxy escape fractions imply a large contribution from faint AGN (down to $M_\mathrm{UV}=-12$) which results in a rapid reionization, disfavoured by current observations.
We infer the UV luminosities of Local Group galaxies at early cosmic times ($z \sim 2$ and $z \sim 7$) by combining stellar population synthesis modeling with star formation histories derived from deep color-magnitude diagrams constructed from Hubble Space Telescope (HST) observations. Our analysis provides a basis for understanding high-$z$ galaxies - including those that may be unobservable even with the James Webb Space Telescope (JWST) - in the context of familiar, well-studied objects in the very low-$z$ Universe. We find that, at the epoch of reionization, all Local Group dwarfs were less luminous than the faintest galaxies detectable in deep HST observations of blank fields. We predict that JWST will observe $z \sim 7$ progenitors of galaxies similar to the Large Magellanic Cloud today; however, the HST Frontier Fields initiative may already be observing such galaxies, highlighting the power of gravitational lensing. Consensus reionization models require an extrapolation of the observed blank-field luminosity function at $z \approx 7$ by at least two orders of magnitude in order to maintain reionization. This scenario requires the progenitors of the Fornax and Sagittarius dwarf spheroidal galaxies to be contributors to the ionizing background at $z \sim 7$. Combined with numerical simulations, our results argue for a break in the UV luminosity function from a faint-end slope of $\alpha \sim -2$ at $M_{\rm UV} < -13$ to $\alpha \sim -1.2$ at lower luminosities. Applied to photometric samples at lower redshifts, our analysis suggests that HST observations in lensing fields at $z \sim 2$ are capable of probing galaxies with luminosities comparable to the expected progenitor of Fornax.
Pulsars are famous for their rotational stability. Most of them steadily spin down and display a highly repetitive pulse shape. But some pulsars experience timing irregularities such as nulling, intermittency, mode changing and timing noise. As changes in the pulse shape are often correlated with timing irregularities, precession is a possible cause of these phenomena. Whereas pulsar magnetospheres are filled with plasma, most pulsar precession studies were carried out within the vacuum approximation and neglected the effects of magnetospheric currents and charges. Recent numerical simulations of plasma-filled pulsar magnetospheres provide us with a detailed quantitative description of magnetospheric torques exerted on the pulsar surface. In this paper, we present the study of neutron star evolution using these new torque expressions. We show that they lead to (1) much slower long-term evolution of pulsar parameters and (2) much less extreme solutions for these parameters than the vacuum magnetosphere models. To facilitate the interpretation of observed pulsar timing residuals, we derive an analytic model that (1) describes the time evolution of non-spherical pulsars and (2) translates the observed pulsar timing residuals into the geometrical parameters of the pulsar. We apply this model to two pulsars with very different temporal behaviours. For the pulsar B1828-11, we demonstrate that the timing residual curves allow two pulsar geometries: one with stellar deformation pointing along the magnetic axis and one along the rotational axis. For the Crab pulsar, we use the model show that the recent observation of its magnetic and rotational axes moving away from each other can be explained by precession.
Photodissociation regions (PDRs) contain a large fraction of all of the interstellar matter in galaxies. Classical examples include the boundaries between ionized regions and molecular clouds in regions of massive star formation, marking the point where all of the photons energetic enough to ionize hydrogen have been absorbed. In this paper we determine the physical properties of the PDRs associated with the star forming regions IRAS 23133+6050 and S 106 and present them in the context of other Galactic PDRs associated with massive star forming regions. We employ Herschel PACS and SPIRE spectroscopic observations to construct a full 55-650 {\mu}m spectrum of each object from which we measure the PDR cooling lines, other fine- structure lines, CO lines and the total far-infrared flux. These measurements are then compared to standard PDR models. Subsequently detailed numerical PDR models are compared to these predictions, yielding additional insights into the dominant thermal processes in the PDRs and their structures. We find that the PDRs of each object are very similar, and can be characterized by a two-phase PDR model with a very dense, highly UV irradiated phase (n $\sim$ 10^6 cm^(-3), G$_0$ $\sim$ 10^5) interspersed within a lower density, weaker radiation field phase (n $\sim$ 10^4 cm^(-3), G$_0$ $\sim$ 10^4). We employed two different numerical models to investigate the data, firstly we used RADEX models to fit the peak of the $^{12}$CO ladder, which in conjunction with the properties derived yielded a temperature of around 300 K. Subsequent numerical modeling with a full PDR model revealed that the dense phase has a filling factor of around 0.6 in both objects. The shape of the $^{12}$CO ladder was consistent with these components with heating dominated by grain photoelectric heating. An extra excitation component for the highest J lines (J > 20) is required for S 106.
To measure the properties of both components of the RS CVn binary sigma Geminorum (sigma Gem), we directly detect the faint companion, measure the orbit, obtain model-independent masses and evolutionary histories, detect ellipsoidal variations of the primary caused by the gravity of the companion, and measure gravity darkening. We detect the companion with interferometric observations obtained with the Michigan InfraRed Combiner (MIRC) at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array with a primary-to-secondary H-band flux ratio of 270+/-70. A radial velocity curve of the companion was obtained with spectra from the Tillinghast Reflector Echelle Spectrograph (TRES) on the 1.5-m Tillinghast Reflector at Fred Lawrence Whipple Observatory (FLWO). We additionally use new observations from the Tennessee State University Automated Spectroscopic and Photometric Telescopes (AST and APT, respectively). From our orbit, we determine model-independent masses of the components (M_1 = 1.28 +/- 0.07 M_Sun, M_2 = 0.73 +/- 0.03 M_Sun), and estimate a system age of 5 -/+ 1 Gyr. An average of the 27-year APT light curve of sigma Gem folded over the orbital period (P = 19.6027 +/- 0.0005 days) reveals a quasi-sinusoidal signature, which has previously been attributed to active longitudes 180 deg apart on the surface of sigma Gem. With the component masses, diameters, and orbit, we find that the predicted light curve for ellipsoidal variations due to the primary star partially filling its Roche lobe potential matches well with the observed average light curve, offering a compelling alternative explanation to the active longitudes hypothesis. Measuring gravity darkening from the light curve gives beta < 0.1, a value slightly lower than that expected from recent theory.
We report the discovery of two super-Earth mass planets orbiting the nearby K0.5 dwarf HD 7924 which was previously known to host one small planet. The new companions have masses of 7.9 and 6.4 M$_\oplus$, and orbital periods of 15.3 and 24.5 days. We perform a joint analysis of high-precision radial velocity data from Keck/HIRES and the new Automated Planet Finder Telescope (APF) to robustly detect three total planets in the system. We refine the ephemeris of the previously known planet using five years of new Keck data and high-cadence observations over the last 1.3 years with the APF. With this new ephemeris, we show that a previous transit search for the inner-most planet would have covered 70% of the predicted ingress or egress times. Photometric data collected over the last eight years using the Automated Photometric Telescope shows no evidence for transits of any of the planets, which would be detectable if the planets transit and their compositions are hydrogen-dominated. We detect a long-period signal that we interpret as the stellar magnetic activity cycle since it is strongly correlated with the Ca II H and K activity index. We also detect two additional short-period signals that we attribute to rotationally-modulated starspots and a one month alias. The high-cadence APF data help to distinguish between the true orbital periods and aliases caused by the window function of the Keck data. The planets orbiting HD 7924 are a local example of the compact, multi-planet systems that the Kepler Mission found in great abundance.
We report the results of a deep SCUBA-2 850- and 450-$\mu$m survey for dust-obscured ultra-luminous infrared galaxies (U/LIRGs) in the field of the z=1.46 cluster XCS J2215.9-1738. We detect a striking overdensity of sub-millimeter sources coincident with the core of this cluster: $\sim 3-4 \times$ higher than expected in a blank field. We use the likely radio and mid-infrared counterparts to show that the bulk of these sub-millimeter sources have spectroscopic or photometric redshifts which place them in the cluster and that their multi-wavelength properties are consistent with this association. The average far-infrared luminosities of these galaxies are $(1.0\pm0.1) \times 10^{12} L_{\odot}$, placing them on the U/LIRG boundary. Using the total star formation occurring in the obscured U/LIRG population within the cluster we show that the resulting mass-normalized star-formation rate for this system supports previous claims of a rapid increase in star-formation activity in cluster cores out to $z\sim1.5$, which must be associated with the on-going formation of the early-type galaxies which reside in massive clusters today.
We study the effect of new emerging solar active regions on the large-scale magnetic environment of existing regions. We first present a theoretical approach to quantify the "interaction energy" between new and pre-existing regions as the difference between (i) the summed magnetic energies of their individual potential fields and (ii) the energy of their superposed potential fields. We expect that this interaction energy can, depending upon the relative arrangements of newly emerged and pre-existing magnetic flux, indicate the existence of "topological" free magnetic energy in the global coronal field that is independent of any "internal" free magnetic energy due to coronal electric currents flowing within the newly emerged and pre-existing flux systems. We then examine the interaction energy in two well-studied cases of flux emergence, but find that the predicted energetic perturbation is relatively small compared to energies released in large solar flares. Next, we present an observational study on the influence of the emergence of new active regions on flare statistics in pre-existing active regions, using NOAA's Solar Region Summary and GOES flare databases. As part of an effort to precisely determine the emergence time of active regions in a large event sample, we find that emergence in about half of these regions exhibits a two-stage behavior, with an initial gradual phase followed by a more rapid phase. Regarding flaring, we find that newly emerging regions produce a significant increase in the occurrence rate of X- and M-class flares in pre-existing regions. Given the relative weakness of the interaction energy, this effect suggests that perturbations in the large-scale magnetic field, such as topology changes invoked the "breakout" model of coronal mass ejections, might play a significant role in the occurrence of some flares.
When measuring star formation rates using ultraviolet light, correcting for dust extinction is a critical step. However, with the variety of dust extinction curves to choose from, the extinction correction is quite uncertain. Here, we use Swift/UVOT to measure the extinction curve for star-forming regions in the SMC and M33. We find that both the slope of the curve and the strength of the 2175 Angstrom bump vary across both galaxies. In addition, as part of our modeling, we derive a detailed recent star formation history for each galaxy.
Using observations obtained with the LOw Fequency ARray (LOFAR), the
Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA)
data, we have traced the radio emission to large scales in the complex source
4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial
resolution (~4"), the source was known to have two inner radio lobes spanning
31 kpc and a diffuse, low-brightness extension running parallel to them, offset
by about 11 kpc (in projection).
At 62 MHz, we detect the radio emission of this structure extending out to
210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure
appears to have a helical morphology.
We have derived the characteristics of the radio spectral index across the
source. We show that the source morphology is most likely the result of at
least two episodes of AGN activity separated by a dormant period of around 35
Myr.
The AGN is hosted by one of the galaxies located in the cluster core of Abell
407. We propose that it is intermittently active as it moves in the dense
environment in the cluster core. Using LOFAR, we can trace the relic plasma
from that episode of activity out to greater distances from the core than ever
before.
Using the the WSRT, we detect HI in absorption against the center of the
radio source. The absorption profile is relatively broad (FWHM of 288 km/s),
similar to what is found in other clusters.
Understanding the duty cycle of the radio emission as well as the triggering
mechanism for starting (or restarting) the radio-loud activity can provide
important constraints to quantify the impact of AGN feedback on galaxy
evolution. The study of these mechanisms at low frequencies using morphological
and spectral information promises to bring new important insights in this
field.
In the theory of star formation, the first hydrostatic core (FHSC) phase is a
critical step in which a condensed object emerges from a prestellar core. This
step lasts about one thousand years, a very short time compared with the
lifetime of prestellar cores, and therefore is hard to detect unambiguously.
We present IRAM Plateau de Bure observations of the Barnard 1b dense
molecular core, combining detections of H2CO and CH3OH spectral lines and dust
continuum at 2.3" resolution (~ 500 AU). The two compact cores B1b-N and B1b-S
are detected in the dust continuum at 2mm, with fluxes that agree with their
spectral energy distribution. Molecular outflows associated with both cores are
detected. They are inclined relative to the direction of the magnetic field, in
agreement with predictions of collapse in turbulent and magnetized gas with a
ratio of mass to magnetic flux somewhat higher than the critical value, \mu ~ 2
- 7. The outflow associated with B1b-S presents sharp spatial structures, with
ejection velocities of up to ~ 7 kms from the mean velocity. Its dynamical age
is estimated to be ~2000 yrs. The B1b-N outflow is smaller and slower, with a
short dynamical age of ~1000 yrs. The B1b-N outflow mass, mass-loss rate, and
mechanical luminosity agree well with theoretical predictions of FHSC. These
observations confirm the early evolutionary stage of B1b-N and the slightly
more evolved stage of B1b-S.
In recent years, the Lyman-$\alpha$ absorption observed in the spectra of high-redshift quasars has been used as a tracer of large-scale structure by means of the three-dimensional Lyman-$\alpha$ forest auto-correlation function at redshift $z\simeq 2.3$, but the need to fit the quasar continuum in every absorption spectrum introduces a broadband distortion that is difficult to correct and causes a systematic error for measuring any broadband properties. We describe a $k$-space model for this broadband distortion based on a multiplicative correction to the power spectrum of the transmitted flux fraction that suppresses power on scales corresponding to the typical length of a Lyman-$\alpha$ forest spectrum. Implementing the distortion model in fits for the baryon acoustic oscillation (BAO) peak position in the Lyman-$\alpha$ forest auto-correlation, we find that the fitting method recovers the input values of the linear bias parameter $b_{F}$ and the redshift-space distortion parameter $\beta_{F}$ for mock data sets with a systematic error of less than 0.5\%. Applied to the auto-correlation measured for BOSS Data Release 11, our method improves on the previous treatment of broadband distortions in BAO fitting by providing a better fit to the data using fewer parameters and reducing the statistical errors on $\beta_{F}$ and the combination $b_{F}(1+\beta_{F})$ by more than a factor of seven. The measured values at redshift $z=2.3$ are $\beta_{F}=1.39^{+0.11\ +0.24\ +0.38}_{-0.10\ -0.19\ -0.28}$ and $b_{F}(1+\beta_{F})=-0.374^{+0.007\ +0.013\ +0.020}_{-0.007\ -0.014\ -0.022}$ (1$\sigma$, 2$\sigma$ and 3$\sigma$ statistical errors). Our fitting software and the input files needed to reproduce our main results are publicly available.
We report the late-time evolution of Type IIb Supernova (SN IIb) 2013df. SN 2013df showed a dramatic change in its spectral features at ~1 year after the explosion. Early on it showed typical characteristics shared by SNe IIb/Ib/Ic dominated by metal emission lines, while later on it was dominated by broad and flat-topped Halpha and He I emissions. The late-time spectra are strikingly similar to SN IIb 1993J, which is the only previous example clearly showing the same transition. This late-time evolution is fully explained by a change in the energy input from the 56Co decay to the interaction between the SN ejecta and dense circumstellar matter (CSM). The mass loss rate is derived to be (~5.4 +- 3.2) x 10^{-5} Msun/yr (for the wind velocity of ~20 km/s, similar to SN 1993J but larger than SN IIb 2011dh by an order of magnitude. The striking similarity between SNe 2013df and 1993J in the (candidate) progenitors and the CSM environments, and the contrast in these natures to SN 2011dh, infer that there is a link between the natures of the progenitor and the mass loss: SNe IIb with a more extended progenitor have experienced a much stronger mass loss in the final centuries toward the explosion. It might indicate that SNe IIb from a more extended progenitor are the explosions during a strong binary interaction phase, while those from a less extended progenitor have a delay between the strong binary interaction and the explosion.
Interpreting the spectra of brown dwarfs is key to determining the fundamental physical and chemical processes occurring in their atmospheres. Powerful Bayesian atmospheric retrieval tools have recently been applied to both exoplanet and brown dwarf spectra to tease out the thermal structures and molecular abundances to understand those processes. In this manuscript we develop a significantly upgraded retrieval method and apply it to the SpeX spectral library data of two benchmark late T-dwarfs, Gl570D and HD3651B, to establish the validity of our upgraded forward model parameterization and Bayesian estimator. Our retrieved metallicities, gravities, and effective temperature are consistent with the metallicity and presumed ages of the systems. We add the carbon-to-oxygen ratio as a new dimension to benchmark systems and find good agreement between carbon-to-oxygens ratio derived in the brown dwarfs and the host stars. Furthermore, we have for the first time unambiguously determined the presence of ammonia in the low-resolution spectra of these two late T-dwarfs. We also show that the retrieved results are not significantly impacted by the possible presence of clouds, though some quantities are significantly impacted by uncertainties in photometry. This investigation represents a watershed study in establishing the utility of atmospheric retrieval approaches on brown dwarf spectra.
We present the stellar and gas kinematics of DDO 46 and DDO 168 from the LITTLE THINGS survey and determine their respective Vmax/sigma_z,0 values. We used the KPNO's 4-meter telescope with the Echelle spectrograph as a long-slit spectrograph. We acquired spectra of DDO 168 along four position angles by placing the slit over the morphological major and minor axes and two intermediate position angles. However, due to poor weather conditions during our observing run for DDO 46, we were able to extract only one useful data point from the morphological major axis. We determined a central stellar velocity dispersion perpendicular to the disk, sigma_z,0, of 13.5+/-8 km/s for DDO 46 and <sigma_z,0> of 10.7+/-2.9 km/s for DDO 168. We then derived the maximum rotation speed in both galaxies using the LITTLE THINGS HI data. We separated bulk motions from non-circular motions using a double Gaussian decomposition technique and applied a tilted-ring model to the bulk velocity field. We corrected the observed HI rotation speeds for asymmetric drift and found a maximum velocity, Vmax, of 77.4 +/- 3.7 and 67.4 +/- 4.0 km/s for DDO 46 and DDO 168, respectively. Thus, we derived a kinematic measure, Vmax/sigma_z,0, of 5.7 +/- 0.6 for DDO 46 and 6.3 +/- 0.3 for DDO 168. Comparing these values to ones determined for spiral galaxies, we find that DDO 46 and DDO 168 have Vmax/sigma_z,0 values indicative of thin disks, which is in contrast to minor-to-major axis ratio studies.
We use the integrated polarized radio emission at 1.4 GHz ($\Pi_{\rm 1.4\,GHz}$) from a large sample of AGN (796 sources at redshifts $z<0.7$) to study the large-scale magnetic field properties of radio galaxies in relation to the host galaxy accretion state. We find a fundamental difference in $\Pi_{\rm 1.4\,GHz}$ between radiative-mode AGN (i.e. high-excitation radio galaxies, HERGs, and radio-loud QSOs) and jet-mode AGN (i.e. low-excitation radio galaxies, LERGs). While LERGs can achieve a wide range of $\Pi_{\rm 1.4\,GHz}$ (up to $\sim$$30\%$), the HERGs and radio-loud QSOs are limited to $\Pi_{\rm 1.4\,GHz} \lesssim 15\%$. A difference in $\Pi_{\rm 1.4\,GHz}$ is also seen when the sample is divided at 0.5% of the total Eddington-scaled accretion rate, where the weakly accreting sources can attain higher values of $\Pi_{\rm 1.4\,GHz}$. We do not find any clear evidence that this is driven by intrinsic magnetic field differences of the different radio morphological classes. Instead, we attribute the differences in $\Pi_{\rm 1.4\,GHz}$ to the local environments of the radio sources, in terms of both the ambient gas density and the magnetoionic properties of this gas. Thus, not only are different large-scale gaseous environments potentially responsible for the different accretion states of HERGs and LERGs, we argue that the large-scale magnetised environments may also be important for the formation of powerful AGN jets. Upcoming high angular resolution and broadband radio polarization surveys will provide the high precision Faraday rotation measure and depolarization data required to robustly test this claim.
KCDC, the KASCADE Cosmic-ray Data Centre, is a web portal, where data of astroparticle physics experiments will be made available for the interested public. The KASCADE experiment, financed by public money, was a large-area detector for the measurement of high-energy cosmic rays via the detection of air showers. KASCADE and its extension KASCADE-Grande stopped finally the active data acquisition of all its components including the radio EAS experiment LOPES end of 2012 after more than 20 years of data taking. In a first release, with KCDC we provide to the public the measured and reconstructed parameters of more than 160 million air showers. In addition, KCDC provides the conceptional design, how the data can be treated and processed so that they are also usable outside the community of experts in the research field. Detailed educational examples make a use also possible for high-school students and early stage researchers.
We study the correlation between gamma-ray and radio band radiation for 123 blazars, using the Fermi-LAT first source catalog (1FGL) and the RATAN-600 data obtained at the same period of time (within a few months). We found an apparent positive correlation for BL Lac and flat-spectrum radio quasar (FSRQ) sources from our sample through testing the value of the Pearson product-moment correlation coefficient. The BL Lac objects show higher values of the correlation coefficient than FSRQs at all frequencies, except 21.7 GHz, and at all bands, except $10-100$ GeV, typically at high confidence level (> 99%). At higher gamma-ray energies the correlation weakens and even becomes negative for BL Lacs and FSRQs. For BL Lac blazars, the correlation of the fluxes appeared to be more sensitive to the considered gamma-ray energy band, than to the frequency, while for FSRQ sources the correlation changed notably both with the considered radio frequency and gamma-ray energy band. We used a data randomization method to quantify the significance of the computed correlation coefficients. We find that the statistical significance of the correlations we obtained between the flux densities at all frequencies and the photon flux in all gamma-ray bands below 3 GeV is high for BL Lacs (chance probability $\sim 10^{-3} - 10^{-7}$). The correlation coefficient is high and significant for the $0.1-0.3$ GeV band and low and insignificant for the $10-100$ GeV band for both types of blazars for all considered frequencies.
Visual observations of meteor showers have been carried out for many decades by astronomers. Modern developments in imaging systems have advanced our knowledge about the shower dynamics and their origin. We have made an attempt to make a portable video recording system to observe meteor showers. The system consists of a camera capable of recording the rapid motion of meteors entering the Earths atmosphere. The camera is interfaced with a laptop using an open source software like VirtualDub. The initial testing was carried out while observing Geminid meteor shower in December, 2010 from the base of Mahuli fort (Lat: $19.47{\deg}$ N, Long: $73.26{\deg}$ E) near Asangaon, India. Here, we present few of the meteors recorded during this event, followed by a preliminary analysis of the shower. The portable video recording system enables to capture meteors at a remote location. This system will strengthen traditional visual observation methods used by astronomers in India. Future coordinated studies with a multi-station approach using such systems will assist in deriving the parameters associated with meteor shower activity and its impact on the Earth's atmosphere. Therefore, we propose long term and simultaneous multi-station video observations in the Indian subcontinent for continuous monitoring and better understanding of meteor dynamics.
Results from long-term multicolor optical photometric observations of the pre-main sequence stars FHO 26, FHO 27, FHO 28, FHO 29 and V1929 Cyg collected during the period from June 1997 to December 2014 are presented. The objects are located in the dense molecular cloud L935, named "Gulf of Mexico", in the field between the North America and Pelican nebulae. All stars from our study exhibit strong photometric variability in all optical passbands. Using our BVRI observations and data published by other authors, we tried to define the reasons for the observed brightness variations. The presented paper is a part of our long-term photometric study of the young stellar objects in the region of "Gulf of Mexico".
We present an application of a particular machine-learning method (Boosted Decision Trees, BDTs using AdaBoost) to separate stars and galaxies in photometric images using their catalog characteristics. BDTs are a well established machine learning technique used for classification purposes. They have been widely used specially in the field of particle and astroparticle physics, and we use them here in an optical astronomy application. This algorithm is able to improve from simple thresholding cuts on standard separation variables that may be affected by local effects such as blending, badly calculated background levels or which do not include information in other bands. The improvements are shown using the Sloan Digital Sky Survey Data Release 9, with respect to the type photometric classifier. We obtain an improvement in the impurity of the galaxy sample of a factor 2-4 for this particular dataset, adjusting for the same efficiency of the selection. Another main goal of this study is to verify the effects that different input vectors and training sets have on the classification performance, the results being of wider use to other machine learning techniques.
In this paper for studying the anisotropic strange quark stars, we assume that the radial pressure inside the anisotropic star is a superposition of pressure in an isotropic case plus a Gaussian perturbation term. Considering a proportionality between electric charge density and the density of matter, we solve the TOV equation for different cases numerically. Our results indicate that anisotropy increases the maximum mass $M_{max}$ and also its corresponding radius $R$ for a typical strange quark star. According to our calculations, an anisotropy amplitude of $A=3\times10^{33}Nm^{-2}$ with a standard deviation of $\sigma=3\times10^{3}m$ leads to a neutron star of 1.97$M_{\odot}$. Furthermore, electric charge not only increases the maximum mass and its corresponding radius, but also raises up the anisotropy factor. We can see that the tangential pressure $p_{t}$ and anisotropy factor $\Delta$ unlike the radial pressure $p_{r}$ have a maximum on the surface and this maximum increases by adding electric charge effect. However, we show that anisotropy can be more effective than electric charge in rasing maximum mass of strange quark stars.
Over the past six years, the Fermi Large Area Telescope has detected more than 150 gamma-ray pulsars, discovering a variety of light curve trends and classes. Such diversity hints at distinct underlying magnetospheric and/or emission geometries. We implemented an offset-dipole magnetic field, with an offset characterised by parameters epsilon and magnetic azimuthal angle, in an existing geometric pulsar modelling code which already includes static and retarded vacuum dipole fields. We use these different magnetic field solutions in conjunction with standard emission geometries, namely the two-pole caustic and outer gap models (the latter only for non-offset dipoles), and construct intensity maps and light curves for several pulsar parameters. We compare our model light curves to the Vela data from the second pulsar catalogue of Fermi. We use a refined chi-square grid search method for finding best-fit light curves for each of the different models. Our best fit is for the retarded vacuum dipole field and the outer gap model.
In this Letter, we investigate how galaxy mass assembly mode depends on stellar mass $M_{\ast}$, using a large sample of $\sim$10, 000 low redshift galaxies. Our galaxy sample is selected to have SDSS $R_{90}>5\arcsec.0$, which allows the measures of both the integrated and the central NUV$-r$ color indices. We find that: in the $M_{\ast}-($ NUV$-r$) green valley, the $M_{\ast}<10^{10}~M_{\sun}$ galaxies mostly have positive or flat color gradients, while most of the $M_{\ast}>10^{10.5}~M_{\sun}$ galaxies have negative color gradients. When their central $D_{n}4000$ index values exceed 1.6, the $M_{\ast}<10^{10.0}~M_{\sun}$ galaxies have moved to the UV red sequence, whereas a large fraction of the $M_{\ast}>10^{10.5}~M_{\sun}$ galaxies still lie on the UV blue cloud or the green valley region. We conclude that the main galaxy assembly mode is transiting from "the outside-in" mode to "the inside-out" mode at $M_{\ast}< 10^{10}~M_{\sun}$ and at $M_{\ast}> 10^{10.5}~M_{\sun}$. We argue that the physical origin of this is the compromise between the internal and the external process that driving the star formation quenching in galaxies. These results can be checked with the upcoming large data produced by the on-going IFS survey projects, such as CALIFA, MaNGA and SAMI in the near future.
We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (Ultra) Luminous Infrared Galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited (v2=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness >5e13 Lsun/kpc2. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, v2=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses and physical conditions of (U)LIRG nuclei when H2 column densities exceed 1e24 cm-2. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, the ground vibrational state (v=0), J=3-2 and 4-3 rotational lines of HCN and HCO+ fail to probe the highly enshrouded, compact nuclear regions due to strong self- and continuum absorption. The HCN and HCO+ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early type spiral LIRGs and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.
We present an analytical model for light echoes (LEs) coming from circumstellar material (CSM) around Type Ia Supernovae (SNe Ia). Using this model we find two spectral signatures at 4100 {\AA} and 6200 {\AA} that are useful to identify LEs during the Lira law phase (between 35 and 80 days after maximum light) coming from nearby CSM at distances of 0.01-0.25 pc. We analyze a sample of 89 SNe Ia divided in two groups according to their B-V decline rate during the Lira law phase, and search for LEs from CSM interaction in the group of SNe with steeper slopes by comparing their spectra with our LE model. We find that a model with LEs + pure extinction from interstellar material (ISM) fits better the observed spectra than a pure ISM extinction model that is constant in time, but we find that a decreasing extinction alone explains better the observations without the need of LEs, possibly implying dust sublimation due to the radiation from the SN.
We describe a "spatio-spectral" deconvolution algorithm for wide-band imaging in radio interferometry. In contrast with the existing multi-frequency reconstruction algorithms, the proposed method does not rely on a model of the sky-brightness spectral distribution. This non-parametric approach can be of particular interest for the new generation of low frequency radiotelescopes. The proposed solution formalizes the reconstruction problem as a convex optimization problem with spatial and spectral regularizations. The efficiency of this approach has been already proven for narrow-band image reconstruction and the present contribution can be considered as its extension to the multi-frequency case. Because the number of frequency bands multiplies the size of the inverse problem, particular attention is devoted to the derivation of an iterative large scale optimization algorithm. It is shown that the main computational bottleneck of the approach, which lies in the resolution of a linear system, can be efficiently overcome by a fully parallel implementation w.r.t. the frequencies, where each processor reconstructs a narrow-band image. All the other optimization steps are extremely fast. A parallel implementation of the algorithm in Julia is publicly available at https://github.com/andferrari. Preliminary simulations illustrate the performances of the method and its ability to reconstruct complex spatio-spectral structures.
We consider the generation of primordial magnetic fields in a class of bouncing models when the electromagnetic action is coupled non-minimally to a scalar field that, say, drives the background evolution. For scale factors that have the power law form at very early times and non-minimal couplings which are simple powers of the scale factor, one can easily show that scale invariant spectra for the magnetic fields can arise {\it before the bounce} for certain values of the indices involved. It will be interesting to examine if these power spectra retain their shape {\it after the bounce}. However, analytical solutions for the Fourier modes of the electromagnetic vector potential across the bounce are difficult to obtain. In this work, with the help of a new time variable that we introduce, which we refer to as the ${\rm e}$-${\cal N}$-fold, we investigate these scenarios numerically. Imposing the initial conditions on the modes in the contracting phase, we numerically evolve the modes across the bounce and evaluate the spectra of the electric and magnetic fields at suitably late times. As one could have intuitively expected, though the complete spectra depend on the details of the bounce, we find that, under the original conditions, scale invariant spectra of the magnetic fields do arise for wavenumbers much smaller than the scale associated with the bounce. We also show that magnetic fields which correspond to observed strengths today can be generated for specific values of the parameters. But, we find that, at the bounce, the backreaction due to the electromagnetic modes that have been generated can be significantly large calling into question the viability of the model. We briefly discuss the implications of our results.
We present Keck/MOSFIRE observations of UV metal lines in four bright gravitationally-lensed z~6-8 galaxies behind the cluster Abell 1703. The spectrum of A1703-zd6, a highly-magnified star forming galaxy with a Lyman-alpha redshift of z=7.045, reveals a confident detection of the nebular CIV emission line (unresolved with FWHM < 125 km/s). UV metal lines are not detected in the three other galaxies. At z~2-3, nebular CIV emission is observed in just 1% of UV-selected galaxies. The presence of strong CIV emission in one of the small sample of galaxies targeted in this paper may indicate hard ionizing spectra are more common at z~7. The total estimated equivalent width of the CIV doublet (38 A) and CIV/Lyman-alpha flux ratio (0.3) are comparable to measurements of narrow-lined AGNs. Photoionization models show that the nebular CIV line can also be reproduced by a young stellar population, with very hot metal poor stars dominating the photon flux responsible for triply ionizing carbon. Regardless of the origin of the CIV, we show that the ionizing spectrum of A1703-zd6 is different from that of typical galaxies at z~2, producing more H ionizing photons per unit 1500A luminosity and a larger flux density at 30-50 eV. If such extreme radiation fields are typical in UV-selected systems at z>7, it would indicate that reionization-era galaxies are more efficient ionizing agents than previously thought. Alternatively, we suggest that the small sample of Lyman-alpha emitters at z>7 may trace a rare population with intense radiation fields capable of ionizing their surrounding hydrogen distribution. Additional constraints on high ionization emission lines in galaxies with and without Lyman-alpha detections will help clarify whether hard ionizing spectra are common in the reionization era.
We present measurements of both scale- and time-dependent deviations from the standard gravitational field equations. These late-time modifications are introduced separately for relativistic and non-relativistic particles, by way of the parameters $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ using two bins in both scale and time, with transition wavenumber $0.01$ Mpc$^{-1}$ and redshift 1. We emphasize the use of two dynamical probes to constrain this set of parameters, galaxy power spectrum multipoles and the direct peculiar velocity power spectrum, which probe fluctuations on different scales. The multipole measurements are derived from the WiggleZ and BOSS Data Release 11 CMASS galaxy redshift surveys and the velocity power spectrum is measured from the velocity sub-sample of the 6-degree Field Galaxy Survey. We combine with additional cosmological probes including baryon acoustic oscillations, Type Ia SNe, the cosmic microwave background (CMB), lensing of the CMB, and the temperature--galaxy cross-correlation. Using a Markov Chain Monte Carlo likelihood analysis, we find the inferred best-fit parameter values of $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ to be consistent with the standard model at the $95\%$ confidence level. Furthermore, accounting for the Alcock-Paczynski effect, we perform joint fits for the expansion history and growth index gamma; we measure $\gamma = 0.665 \pm 0.0669$ ($68\%$ C.L) for a fixed expansion history, and $\gamma = 0.73^{+0.08}_{-0.10}$ ($68\%$ C.L) when the expansion history is allowed to deviate from $\Lambda$CDM. With a fixed expansion history the inferred value is consistent with GR at the $95\%$ C.L; alternatively, a $2\sigma$ tension is observed when the expansion history is not fixed, this tension is worsened by the combination of growth and SNe data.
The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.
The precise spectra of Cosmic Ray (CR) electrons and positrons have been published by the measurement of AMS-02. It is reasonable to regard the difference between the electrons and positrons spectra ( $\triangle \Phi= \Phi_{e^-}-\Phi_{e^+}$ ) as being dominated by primary electrons. Noticing that the resulting electron spectrum shows no sign of spectral softening above 20 GeV, which is in contrast with the prediction of standard model. In this work, we generalize the analytic one dimensional two-halo model of diffusion to a three dimensional realistic calculation by implementing a spatial variant diffusion coefficients in DRAGON package. As a result, we can reproduce the spectral hardening of protons observed by several experiments, and predict an excess of high energy primary electrons which agrees with the measurement reasonably well. Unlike the break spectrum obtained for protons, the model calculation predicts a smooth electron excess and thus slightly over predicts the flux from tens of GeV to 100GeV. To understand this issue, further experimental and theoretical studies are necessary.
The power spectrum of the cosmic microwave background from both the Planck and WMAP data exhibits a slight dip in for multipoles in the range of l=10-30. We show that such a dip could be the result of resonant creation of a massive particle that couples to the inflaton field. For our best-fit models, epochs of resonant particle creation reenters the horizon at wave numbers of k* ~ 0.00011 (h/Mpc). The amplitude and location of these features correspond to the creation of a number of degenerate fermion species of mass ~ 15 times the planck mass during inflation with a coupling constant between the inflaton field and the created fermion species of near unity. Although the evidence is marginal, if this interpretation is correct, this could be one of the first observational hints of new physics at the Planck scale.
We systematically show that in potential driven generalized G-inflation models, quantum corrections coming from new physics at the strong coupling scale can be avoided, while producing observable tensor modes. The effective action can be approximated by the tree level action, and as a result, these models are internally consistent, despite the fact that we introduced new mass scales below the energy scale of inflation. Although observable tensor modes are produced with sub-strong coupling scale field excursions, this is not an evasion of the Lyth bound, since the models include higher-derivative non-canonical kinetic terms, and effective rescaling of the field would result in super-Planckian field excursions. We argue that the enhanced kinetic term of the inflaton screens the interactions with other fields, keeping the system weakly coupled during inflation.
We explore the origin of the observed Lya and other UV lines from galaxies at z<3.7 by detailed modelling of the spectra. The objects are chosen among those showing a) UV-optical-near-IR lines, b) only UV lines and c) those showing Lya in the UV and a few optical lines. We also present UV line predictions for a sample of galaxies in the 0.0686<z<0.8829 range. The sample of galaxies including Lya observations in their spectra does not show peculiar physical conditions of the emitting gas, nor abnormal element abundances. However, the high velocity (Vs>1000 km/s) component of the emitting gas is accompanied by relatively low preshock densities (n0 ~100-400 cm^-3) leading in some cases to broad forbidden lines. Some spectra are best reproduced by shock dominated models in which the photoionizing source is hidden or absent. Within more than 50 galaxies modelled in this work, only a few spectra from galaxies at z~2.5 correspond to a starburst temperature Ts>10^5 K, similar to that found in galaxies showing some activity.
Galaxy distances and derived radial peculiar velocity catalogs constitute valuable datasets to study the dynamics of the Local Universe. However, such catalogs suffer from biases whose effects increase with the distance. Malmquist biases and lognormal error distribution affect the catalogs. Velocity fields of the Local Universe reconstructed with these catalogs present a spurious overall infall onto the Local Volume if they are not corrected for biases. Such an infall is observed in the reconstructed velocity field obtained when applying the BayesianWiener-Filter technique to the raw second radial peculiar velocity catalog of the Cosmicflows project. In this paper, an iterative method to reduce spurious non-Gaussianities in the radial peculiar velocity distribution, to retroactively derive overall better distance estimates resulting in a minimization of the effects of biases, is presented. This method is tested with mock catalogs. To control the cosmic variance, mocks are built out of different cosmological constrained simulations which resemble the Local Universe. To realistically reproduce the effects of biases, the mocks are constructed to be look-alikes of the second data release of the Cosmicflows project, with respect to the size, distribution of data and distribution of errors. Using a suite of mock catalogs, the outcome of the correction is verified to be affected neither by the added error realization, nor by the datapoint selection, nor by the constrained simulation. Results are similar for the different tested mocks. After correction, the general infall is satisfactorily suppressed. The method allows us to obtained catalogs which together with the Wiener-Filter technique give reconstructions approximating non biased velocity fields at 100-150 km/s (2-3 Mpc/h in terms of linear displacement), the linear theory threshold.
We present properties of moderately massive clusters of galaxies detected by the newly developed Hyper Suprime-Cam on the Subaru telescope using weak gravitational lensing. Eight peaks exceeding a S/N ratio of 4.5 are identified on the convergence S/N map of a 2.3 square degree field observed during the early commissioning phase of the camera. Multi-color photometric data is used to generate optically selected clusters using the CAMIRA algorithm. The optical cluster positions were correlated with the peak positions from the convergence map. All eight significant peaks have optical counterparts. The velocity dispersion of clusters are evaluated by adopting the Singular Isothemal Sphere (SIS) fit to the tangential shear profiles, yielding virial mass estimates, M500c, of the clusters which range from 2.7x10^13 to 4.4x10^14 solar mass. The number of peaks is considerably larger than the average number expected from LambdaCDM cosmology but this is not extremely unlikely if one takes the large sample variance in the small field into account. We could, however, safely argue that the peak count strongly favours the recent Planck result suggesting high sigma8$value of 0.83. The ratio of stellar mass to the dark matter halo mass shows a clear decline as the halo mass increases. If the gas mass fraction, fg, in halos is universal, as has been suggested in the literature, the observed baryon mass in stars and gas shows a possible deficit compared with the total baryon density estimated from the baryon oscillation peaks in anisotropy of the cosmic microwave background.
Context: It appears that most (if not all) massive stars are born in multiple systems. At the same time, the most massive binaries are hard to find due to their low numbers throughout the Galaxy and the implied large distances and extinctions. AIMS: We want to study: [a] LS III +46 11, identified in this paper as a very massive binary; [b] another nearby massive system, LS III +46 12; and [c] the surrounding stellar cluster, Berkeley 90. Methods: Most of the data used in this paper are multi-epoch high-S/N optical spectra though we also use Lucky Imaging and archival photometry. The spectra are reduced with devoted pipelines and processed with our own software, such as a spectroscopic-orbit code, CHORIZOS, and MGB. Results: LS III +46 11 is identified as a new very-early-O-type spectroscopic binary [O3.5 If* + O3.5 If*] and LS III +46 12 as another early O-type system [O4.5 V((f))]. We measure a 97.2-day period for LS III +46 12 and derive minimum masses of 38.80$\pm$0.83 M_Sol and 35.60$\pm$0.77 M_Sol for its two stars. We measure the extinction to both stars, estimate the distance, search for optical companions, and study the surrounding cluster. In doing so, a variable extinction is found as well as discrepant results for the distance. We discuss possible explanations and suggest that LS III +46 12 may be a hidden binary system, where the companion is currently undetected.
Compact elliptical galaxies form a rare class of stellar system (~30 presently known) characterized by high stellar densities and small sizes and often harboring metal-rich stars. They were thought to form through tidal stripping of massive progenitors, until two isolated objects were discovered where massive galaxies performing the stripping could not be identified. By mining astronomical survey data, we have now found 195 compact elliptical galaxies in all types of environment. They all share similar dynamical and stellar population properties. Dynamical analysis for nonisolated galaxies demonstrates the feasibility of their ejection from host clusters and groups by three-body encounters, which is in agreement with numerical simulations. Hence, isolated compact elliptical and isolated quiescent dwarf galaxies are tidally stripped systems that ran away from their hosts.
This work investigates the main mechanism(s) that regulate the specific star formation rate (SSFR) in nearby galaxies, cross-correlating two proxies of this quantity -- the equivalent width of the \Ha\ line and the $(u-r)$ colour -- with other physical properties (mass, metallicity, environment, morphology, and the presence of close companions) in a sample of $\sim82500$ galaxies extracted from the Sloan Digital Sky Survey (SDSS). The existence of a relatively tight `ageing sequence' in the colour-equivalent width plane favours a scenario where the secular conversion of gas into stars (i.e. `nature') is the main physical driver of the instantaneous SSFR and the gradual transition from a `chemically primitive' (metal-poor and intensely star-forming) state to a `chemically evolved' (metal-rich and passively evolving) system. Nevertheless, environmental factors (i.e. `nurture') are also important. In the field, galaxies may be temporarily affected by discrete `quenching' and `rejuvenation' episodes, but such events show little statistical significance in a probabilistic sense, and we find no evidence that galaxy interactions are, on average, a dominant driver of star formation. Although visually classified mergers tend to display systematically higher EW(H$\alpha$) and bluer $(u-r)$ colours for a given luminosity, most galaxies with high SSFR have uncertain morphologies, which could be due to either internal or external processes. Field galaxies of early and late morphological types are consistent with the gradual `ageing' scenario, with no obvious signatures of a sudden decrease in their SSFR. In contrast, star formation is significantly reduced and sometimes completely quenched on a short time scale in dense environments, where many objects are found on a `quenched sequence' in the colour-equivalent width plane.
It has been shown that muon flux intensities calculated in terms of the EPOS LHC and EPOS 1.99 models at the energy of 10^4 GeV exceed the data of the classical experiments L3+Cosmic, MACRO and LVD on the spectra of atmospheric muons by a factor of 1.9 and below these data at the same energy by a factor of 1.8 in case of the QGSJET II-03 model. It has been concluded that these tested models overestimate (underestimate in case of QGSJET II-03 model) the production of secondary particles with the highest energies in interactions of hadrons by a factor of ~1.5. The LHCf and TOTEM accelerator experiments show also this type of disagreements with these model predictions at highest energies of secondary particles.
Gamma-ray bursts (GRBs) are widely proposed as an effective probe to trace the Hubble diagram of the Universe in high redshift range. However, the calibration of GRBs is not as easy as that of type-Ia supernovae (SNe Ia). Most calibrating methods at present take use one or some of the empirical luminosity corrections, e.g., Amati relation. One of the underlying assumptions of these calibrating methods is that the empirical correlation is universal over all redshifts. In this paper, we check to what extent this assumption holds. Assuming that SNe Ia exactly trace the Hubble diagram of the Universe, we re-investigate the Amati relation for low redshift ($z<1.4$) and high redshift ($z>1.4$) GRBs, respectively. It is found that the Amati relation of low-$z$ GRBs differs from that of high-$z$ GRBs at more than $3\sigma$ confidence level. This result is insensitive to cosmological models.
The Fermi Bubbles have been imaged in sub-TeV gamma rays at Fermi-LAT, and, if their origin is hadronic, they might have been seen with low statistics in $\sim 0.1- 1$ PeV neutrinos at IceCube. We discuss the detectability of these objects at the new High Altitude Water Cherenkov (HAWC) gamma ray detector. HAWC will view the North Bubble for $\sim 2-3$ hours a day, and will map its spectrum at 0.1-100 TeV. For the hard primary proton spectrum required to explain five events at IceCube, a high significance detection at HAWC will be achieved in less than 30 days. The combination of results at HAWC and IceCube will substantiate the hadronic model, or constrain its spectral parameters.
We study the linear and nonlinear evolution of the tearing instability on thin current sheets by means of two-dimensional numerical simulations, within the framework of compressible, resistive magnetohydrodynamics. In particular we analyze the behavior of current sheets whose inverse aspect ratio scales with the Lundquist number $S$ as $S^{-1/3}$. This scaling has been recently recognized to yield the threshold separating fast, ideal reconnection, with an evolution and growth which are independent of $S$ provided this is high enough, as it should be natural having the ideal case as a limit for $S\to\infty$. Our simulations confirm that the tearing instability growth rate can be as fast as $\gamma\approx 0.6\,{\tau_A}^{-1}$, where $\tau_A$ is the ideal Alfv\'enic time set by the macroscopic scales, for our least diffusive case with $S=10^7$. The expected instability dispersion relation and eigenmodes are also retrieved in the linear regime, for the values of $S$ explored here. Moreover, in the nonlinear stage of the simulations we observe secondary events obeying the same critical scaling with $S$, here calculated on the \emph{local}, much smaller lengths, leading to increasingly faster reconnection. These findings strongly support the idea that in a fully dynamic regime, as soon as current sheets develop, thin and reach this critical threshold in their aspect ratio, the tearing mode is able to trigger plasmoid formation and reconnection on the local (ideal) Alfv\'enic timescales, as required to explain the explosive flaring activity often observed in solar and astrophysical plasmas.
Thanks to the large collecting area (3 x ~1500 cm$^2$ at 1.5 keV) and wide field of view (30' across in full field mode) of the X-ray cameras on board the European Space Agency X-ray observatory XMM-Newton, each individual pointing can result in the detection of hundreds of X-ray sources, most of which are newly discovered. Recently, many improvements in the XMM-Newton data reduction algorithms have been made. These include enhanced source characterisation and reduced spurious source detections, refined astrometric precision of sources, greater net sensitivity for source detection and the extraction of spectra and time series for fainter sources, with better signal-to-noise. Further, almost 50% more observations are in the public domain compared to 2XMMi-DR3, allowing the XMM-Newton Survey Science Centre (XMM-SSC) to produce a much larger and better quality X-ray source catalogue. The XMM-SSC has developed a pipeline to reduce the XMM-Newton data automatically and using improved calibration a new catalogue version has been made from XMM-Newton data made public by 2013 Dec. 31 (13 years of data). Manual screening ensures the highest data quality. This catalogue is known as 3XMM. In the latest release, 3XMM-DR5, there are 565962 X-ray detections comprising 396910 unique X-ray sources. For the 133000 brightest sources, spectra and lightcurves are provided. For all detections, the positions on the sky, a measure of the detection quality, and an evaluation of variability is provided, along with the fluxes and count rates in 7 X-ray energy bands, the total 0.2-12 keV band counts, and four hardness ratios. To identify the detections, a cross correlation with 228 catalogues is also provided for each X-ray detection. 3XMM-DR5 is the largest X-ray source catalogue ever produced. Thanks to the large array of data products, it is an excellent resource in which to find new and extreme objects.
We present polarimetric maser observations with the Australia Telescope Compact Array (ATCA) of excited-state hydroxyl (OH) masers. We observed 30 fields of OH masers in full Stokes polarization with the Compact Array Broadband Backend (CABB) at both the 6030 and 6035 MHz excited-state OH transitions, and the 6668-MHz methanol maser transition, detecting 70 sites of maser emission. Amongst the OH we found 112 Zeeman pairs, of which 18 exhibited candidate {\pi} components. This is the largest single full polarimetric study of multiple sites of star formation for these frequencies, and the rate of 16% {\pi} components clearly indicates the {\pi} component exists, and is comparable to the percentage recently found for ground-state transitions. This significant percentage of {\pi} components, with consistent proportions at both ground- and excited-state transitions, argues against Faraday rotation suppressing the {\pi} component emission. Our simultaneous observations of methanol found the expected low level of polarisation, with no circular detected, and linear only found at the less than or equal to 10% level for the brightest sources.
Our discovery of 1SWASP J093010.78+533859.5 as a probable doubly eclipsing quadruple system containing a contact binary with P~0.23 d and a detached binary with P~1.31 d was announced in 2013. Subsequently Koo et al. confirmed the detached binary spectroscopically and identified a fifth set of static spectral lines at its location, corresponding to a further non-eclipsing component of the system. Here we present new spectroscopic and photometric observations, allowing confirmation of the contact binary and improved modelling of all four eclipsing components. The detached binary is found to contain components of masses 0.837(8) and 0.674(7) M_sol, with radii of 0.832(18) and 0.669(18) R_sol and effective temperatures of 5185(-20,+25) and 4325(-15,+20) K respectively, the contact system has masses 0.86(2) and 0.341(11) M_sol, radii of 0.79(4) and 0.52(5) R_sol respectively, and a common T_eff of 4700(50) K. The fifth star is of similar temperature and spectral type to the primaries in the two binaries. Long-term photometric observations indicate the presence of a spot on one component of the detached binary, moving at an apparent rate of approximately one rotation every two years. Both binaries have consistent system velocities around -11 to -12 km/s, which match the average radial velocity of the fifth star, consistent distance estimates for both subsystems of d=78(3) and d=73(4) pc are also found, and (with some further assumptions) of d=83(9) pc for the fifth star. These findings strongly support the claim that both binaries (and very probably all five stars) are gravitationally bound in a single system. The consistent angles of inclination found for the two binaries (88.2(3) and 86(4) degrees) may also indicate that they originally formed by fragmentation (~9-10 Gyr ago) from a single protostellar disk and subsequently remained in the same orbital plane.
We investigate the mid-infrared properties of the largest (42 objects) sample of radio-loud narrow-line Seyfert 1 (RL NLS1) collected to date, using data from the Wide-field Infrared Survey Explorer (WISE). We analyse the mid-IR colours of these objects and compare them to what is expected from different combinations of AGN and galaxy templates. We find that, in general, the host-galaxy emission gives an importan contribution to the observed mid-IR flux in particular at the longest wavelengths (W3, at 12micron, and W4, at 22micron). In about half of the sources (22 objects) we observe a very red mid-IR colour (W4-W3>2.5) that can be explained only using a starburst galaxy template (M82). Using the 22micron luminosities, corrected for the AGN contribution, we have then estimated the star-formation rate for 20 of these "red" RL NLS1, finding values ranging from 10 to 500 Msun/y. For the RL NLS1 showing bluer colours, instead, we cannot exclude the presence of a star-forming host galaxy although, on average, we expect a lower star-formation rate. Studying the radio (1.4GHz) to mid-IR (22micron) flux ratios of the RL NLS1 in the sample we found that in ~10 objects the star-forming activity could represent the most important component also at radio frequencies, in addition (or in alternative) to the relativistic jet. We conclude that both the mid-IR and the radio emission of RL NLS1 are a mixture of different components, including the relativistic jet, the dusty torus and an intense star-forming activity.
The Super Wide Angle Search for Planets (SuperWASP) is a whole-sky high-cadence optical survey which has searched for exoplanetary transit signatures since 2004. Its archive contains long-term light curves for ~30 million 8-15 V magnitude stars, making it a valuable serendipitous resource for variable star research. We have concentrated on the evidence it provides for eclipsing binaries, in particular those exhibiting orbital period variations, and have developed custom tools to measure periods precisely and detect period changes reliably. Amongst our results are: a collection of 143 candidate contact or semi-detached eclipsing binaries near the short-period limit in the main sequence binary period distribution; a probable hierarchical triple exhibiting dramatic sinusoidal period variations; a new doubly-eclipsing quintuple system; and new evidence for period change or stability in 12 post-common-envelope eclipsing binaries, which may support the existence of circumbinary planets in such systems. A large-scale search for period changes in ~14000 SuperWASP eclipsing binary candidates also yields numerous examples of sinusoidal period change, suggestive of tertiary companions, and may allow us to constrain the frequency of triple systems amongst low-mass stars.
We investigate the cosmological production of gravitational waves for a nonsingular flat cosmology driven by a decaying vacuum energy density evolving as $\rho_{\text{vac}}(H) = \rho_b + H^{3}/H_I$, where $\rho_b$ is the bare vacuum energy density, $H$ is the Hubble parameter and $H_I$ is the primordial inflationary scale. This model can be interpreted as a particular case of the class recently discussed by Perico et al. (Phys. Rev. D 88, 063531, 2013) which is termed complete in the sense that the cosmic evolution occurs between two extreme de Sitter stages (early and late time de Sitter phases). The gravitational wave equation is derived and its time-dependent part numerically integrated since the primordial de Sitter stage. The transition from the early de Sitter to the radiation phase is smooth (no exit problem) and the generated spectrum of gravitons is compared with the standard calculations where an abrupt transition is assumed. It is found that the stochastic background of gravitons is very similar to the one predicted by the cosmic concordance model plus inflation except in the limit of higher frequencies ($\nu \gtrsim 100$ kHz). This remarkable signature of a decaying vacuum cosmology combined with the proposed high frequency gravitational wave detectors of improved sensitivity may provide in the future a crucial test for inflationary mechanisms.
Jellyfish galaxies are galaxies that exhibit tentacles of debris material suggestive of gas stripping. We have conducted the first systematic search for jellyfish galaxies at low-z (z=0.04-0.07) in different environments. We have visually inspected B and V-band images and identified 241+153 candidates in 41+31 galaxy clusters of the OMEGAWINGS+WINGS sample and 99 candidates in groups and lower mass structures in the PM2GC sample. This large sample is well suited for follow-up studies of the gas and for a detailed analysis of the environments where such episodes of gas stripping occur. We present here the atlas of jellyfish candidates, a first analysis of their environment and their basic properties, such as morphologies, star formation rates and galaxy stellar masses. Jellyfish candidates are found in all clusters and at all clustercentric radii, and their number does not correlate with the cluster velocity dispersion or X-ray luminosity. Interestingly, convincing cases of jellyfish candidates are also found in groups and lower mass haloes (10^{11}-10^{14} M_sun). All the candidates are disky, have stellar masses ranging from log M/M_sun < 9 to > 11.5 and the majority of them form stars, at a rate that is on average a factor of 2 higher compared to non-stripped galaxies of similar mass. The few post-starburst and passive candidates have weak tentacles. We conclude that the jellyfish phenomenon is ubiquitous in clusters and can be present even in groups and low mass haloes. Further studies will reveal the physics of the gas stripping and clarify the mechanisms at work.
We employ the Delta-variance analysis and study the turbulent gas dynamics of simulated molecular clouds (MCs). Our models account for a simplified treatment of time-dependent chemistry and the non-isothermal nature of the gas. We investigate simulations using three different initial mean number densities of n_0 = 30, 100 and 300 cm^{-3} that span the range of values typical for MCs in the solar neighbourhood. Furthermore, we model the CO line emission in a post-processing step using a radiative transfer code. We evaluate Delta-variance spectra for centroid velocity (CV) maps as well as for integrated intensity and column density maps for various chemical components: the total, H2 and 12CO number density and the integrated intensity of both the 12CO and 13CO (J = 1 -> 0) lines. The spectral slopes of the Delta-variance computed on the CV maps for the total and H2 number density are significantly steeper compared to the different CO tracers. We find slopes for the linewidth-size relation ranging from 0.4 to 0.7 for the total and H2 density models, while the slopes for the various CO tracers range from 0.2 to 0.4 and underestimate the values for the total and H2 density by a factor of 1.5-3.0. We demonstrate that optical depth effects can significantly alter the Delta-variance spectra. Furthermore, we report a critical density threshold of ~100 cm^{-3} at which the Delta-variance slopes of the various CO tracers change sign. We thus conclude that carbon monoxide traces the total cloud structure well only if the average cloud density lies above this limit.
We use N-body simulations to study the matter distribution in disformal gravity. The disformal model studied here is a conformally coupled symmetron field with an additional exponential disformal term. We conduct cosmological simulations with the aim to find the impact of the new disformal terms in the matter power spectrum, halo mass function and radial profile of the scalar field. This is done by calculating the disformal geodesic equation and the equation of motion for the scalar field, then implementing them into the N-body code ISIS, which is a modified gravity version of the code RAMSES. The presence of a conformal symmetron field increases both the power spectrum and mass function compared to standard gravity on small scales. Our main result is that the newly added disformal terms tend to counteract this effects and can make the evolution slightly closer to standard gravity. We finally show that the disformal terms give rise to oscillations of the scalar field in the centre of the dark matter haloes.
By proposing a pure leptonic radiation model, we study the potential gamma-ray emissions from jets of the low-mass X-ray binaries. In this model, the relativistic electrons that are accelerated in the jets are responsible for radiative outputs. Nevertheless, dynamics of jets are dominated by the magnetic and proton-matter kinetic energies. The model involves all kinds of related radiative processes and considers the evolution of relativistic electrons along the jet by numerically solving the kinetic equation. Numerical results show that the spectral energy distributions can extend up to TeV bands, in which synchrotron radiation and synchrotron self-Compton scattering are dominant components. As an example, we apply the model to the low-mass X-ray binary GX 339-4. The results can not only reproduce the currently available observations from GX 339-4, but also predict detectable radiation at GeV and TeV bands by Fermi and CTA telescopes. The future observations with Fermi and CTA can be used to test our model, which could be employed to distinguish the origin of X-ray emissions.
We present new integral field spectroscopy of the gravitationally lensed broad absorption line (BAL) quasar H1413+117, covering the ultraviolet to visible rest-frame spectral range. We observe strong microlensing signatures in lensed image D, and we use this microlensing to simultaneously constrain both the broad emission and broad absorption line gas. By modeling the lens system over the range of probable lensing galaxy redshifts and using on a new argument based on the wavelength-independence of the broad line lensing magnifications, we determine that there is no significant broad line emission from smaller than ~20 light days. We also perform spectral decomposition to derive the intrinsic broad emission line (BEL) and continuum spectrum, subject to BAL absorption. We also reconstruct the intrinsic BAL absorption profile, whose features allow us to constrain outflow kinematics in the context of a disk-wind model. We find a very sharp, blueshifted onset of absorption of 1,500 km/s in both C IV and N V that may correspond to an inner edge of a disk-wind's radial outflow. The lower ionization Si IV and Al III have higher-velocity absorption onsets, consistent with a decreasing ionization parameter with radius in an accelerating outflow. There is evidence of strong absorption in the BEL component which indicates a high covering factor for absorption over two orders of magnitude in outflow radius.
We have investigated the evaporation of close-in exoplanets irradiated by ionizing photons. We find that the properties of the flow are controlled by the ratio of the recombination time to the flow time-scale. When the recombination time-scale is short compared to the flow time-scale the the flow is in approximate local ionization equilibrium with a thin ionization front, where the photon mean free path is short compared to flow scale. In this "recombination limited" flow the mass-loss scales roughly with the square root of the incident flux. When the recombination time is long compared to the flow time-scale the ionization front becomes thick and encompasses the entire flow, with the mass-loss rate scaling linearly with flux. If the planet's potential is deep the flow is approximately "energy-limited"; however, if the planet's potential is shallow we identify a new limiting mass-loss regime, which we term "photon-limited". In this scenario the mass-loss rate is purely limited by the incoming flux of ionizing photons. We have developed a new numerical approach that takes into account the frequency dependence of the incoming ionizing spectrum and performed a large suite of 1D simulations to characterise UV driven mass-loss around low mass planets. We find the flow is "recombination-limited" at high fluxes but becomes "energy-limited" at low fluxes; however, the transition is broad occurring over several order of magnitude in flux. Finally, we point out the transitions between the different flow types does not occur at a single flux value, but depends on the planet's properties, with higher mass planets becoming "energy-limited" at lower fluxes.
Though megamasers are used to probe extragalactic phenomena, questions remain regarding their production and connection to galactic processes. The observation that water and hydroxyl megamasers rarely coexist in the same galaxy has given rise to the hypothesis that the two megamaser species appear in different phases of nuclear activity or that somehow their presence distinguishes different physical conditions in the medium. However, simultaneous hydroxyl and water megamaser emission has recently been detected in IC 694. Studies of this object are underway but, because many megamasers have not been surveyed for emission in the other molecule, it remains unclear whether IC 694 occupies a narrow phase of galaxy evolution or whether the relationship between megamaser species and galactic processes is more complicated than previously believed. In this paper, we present preliminary results of a systematic search for 22 GHz water maser emission among OH megamaser hosts to identify additional objects hosting both megamaser species to constrain the frequency of the phenomenon. Our work roughly doubles the number of galaxies searched for emission in both molecules which host at least one confirmed maser. We report a definitive ($> 8 \sigma$) detection of water emission toward II Zw 96, firmly establishing it as the second object to co-host both water and hydroxyl megamasers after IC 694. We find high luminosity, narrow features in the water feature in II Zw 96, evidence that at least a part of the maser emission is from an accreting AGN. Additionally, we report evidence for a water megamaser in IRAS 15179+3956. All dual megamaser candidates appear in merging galaxy systems suggestive that megamaser coexistance may signal a brief phase along the merger sequence. Our observations provide possible evidence for an exclusion of H$_2$O kilomasers among OH megamaser hosts.
We discuss the renormalization of the initial value problem in Nonequilibrium Quantum Field Theory within a simple, yet instructive, example and show how to obtain a renormalized time evolution for the two-point functions of a scalar field and its conjugate momentum at all times. The scheme we propose is applicable to systems that are initially far from equilibrium and compatible with non-secular approximation schemes which capture thermalization. It is based on Kadanoff-Baym equations for non-Gaussian initial states, complemented by usual vacuum counterterms. We explicitly demonstrate how various cutoff-dependent effects peculiar to nonequilibrium systems, including time-dependent divergences or initial-time singularities, are avoided by taking an initial non-Gaussian three-point vacuum correlation into account.
Due to Earth's revolution around the Sun, the expected scattering rate in direct dark matter searches is annually modulated. This modulation is expected to differ between experiments when given as a function of recoil energy $E_\text{R}$, e.g. due to the gravitational focusing effect of the Sun; a better variable to compare results among experiments employing different targets is the minimum speed $v_\text{min}$ a dark matter particle must have to impart a recoil energy $E_\text{R}$ to a target nucleus. It is widely believed that the modulation expressed as a function of $v_\text{min}$ is common to all experiments, irrespective of the dark matter distribution. We point out that the annual modulation as a function of $v_\text{min}$, and in particular the times at which the rate is maximum and minimum, could be very different depending on the detector material. This would be an indication of a scattering cross section with non-factorizable velocity and target material dependence. Observing an annual modulation with at least two different target elements would be necessary to identify this type of cross section.
It is thought that the spacetime geometry around black hole candidates is described by the Kerr solution, but an observational confirmation is still missing. Today, the continuum-fitting method and the analysis of the iron K$\alpha$ line cannot unambiguously test the Kerr paradigm because of the degeneracy among the parameters of the system, in the sense that it is impossible with current X-ray data to distinguish a Kerr black hole from a non-Kerr object with different spin and observed from a different viewing angle. In this paper, we study the possibility of testing the Kerr nature of black hole candidates with X-ray spectropolarimetric measurements. As a preliminary study, we employ a model with some simplifications, but our conclusions can unlikely change with a more sophisticated description. We find that -- even in the case of high quality data -- it is impossible to test the Kerr metric and the problem is still the strong correlation between the spin and possible deviations from the Kerr geometry. The possibility of combining spectropolarimetric measurements with those of the continuum and of the iron line does not seem very helpful, because the correlation between the estimates of the spin and of the deformation parameter in the three techniques is very similar. In the end, only good measurements of the iron line may test the Kerr metric, with the other two techniques capable of providing a consistency check.
We investigate which Jordan frame $F(R)$ gravity can describe a Type IV singular bouncing cosmological evolution, with special emphasis given near the point at which the Type IV singularity occurs. The cosmological bounce is chosen in such a way so that the bouncing point coincides exactly with Type IV singularity point. The stability of the resulting $F(R)$ gravity is examined and in addition, we study the Einstein frame scalar-tensor theory counterpart of the resulting Jordan frame $F(R)$ gravity. Also, by assuming that the Jordan frame metric is chosen in such a way so that, when conformally transformed in the Einstein frame, it yields a quasi de Sitter or de Sitter Friedmann-Robertson-Walker metric, we study the observational indexes which turn out to be consistent with Planck 2015 data in the case of the Einstein frame scalar theory. Finally, we study the behavior of the effective equation of state corresponding to the Type IV singular bounce and after we compare the resulting picture with other bouncing cosmologies, we critically discuss the implications of our analysis.
We undertake a hydrodynamical study of a magnetised cosmic fluid between the end of the leptonic era and the beginning of the radiation-dominated epoch. We assume this fluid to be the source of a Bianchi I model and to be a mixture of tightly coupled primordial radiation, neutrinos, baryons, electrons and positrons, together with a gas of already decoupled dark matter WIMPS and an already existing magnetic field. The interaction of this field with the tightly coupled gas mixture is described by suitable equations of state that are appropriate for the particle species of the mixture. Comparison of our results with those of previous studies based on an FLRW framework reveals that the effects of the anisotropy of the magnetic field on the evolution of the main thermodynamical variables are negligible, thus validating these studies, though subtle differences are found in the evolution of the magnetic field itself. For larger field intensities we find quantitative and qualitative differences from the FLRW based analysis. Our approach and our results may provide interesting guidelines in potential situations in which non-perturbative methods are required to study the interaction between magnetic fields and the cosmic fluid.
We present Higgs-Portal dark matter (DM) models to explain the reported Galactic Center GeV gamma-ray excess. Naive effective theories are inconsistent with direct detection constraint for the relevant parameter range. Simple extended models with dark gauge symmetries can easily accommodate the gamma-ray excess through the Higgs-Portal coupling while satisfying various constraints.
Objects that are on the threshold of forming the horizon but never collapse are called quasi-black holes (QBHs). We discuss the properties of the general spherically symmetric QBH metric without addressing its material source, including its limiting cases as the corresponding small parameter tends to zero. We then show that QBHs can exist among self-gravitating configurations of electromagnetic and dilatonic scalar fields without matter. These general results are illustrated by explicit examples of exact solutions.
Meta learning uses information from base learners (e.g. classifiers or estimators) as well as information about the learning problem to improve upon the performance of a single base learner. For example, the Bayes error rate of a given feature space, if known, can be used to aid in choosing a classifier, as well as in feature selection and model selection for the base classifiers and the meta classifier. Recent work in the field of f-divergence functional estimation has led to the development of simple and rapidly converging estimators that can be used to estimate various bounds on the Bayes error. We estimate multiple bounds on the Bayes error using an estimator that applies meta learning to slowly converging plug-in estimators to obtain the parametric convergence rate. We compare the estimated bounds empirically on simulated data and then estimate the tighter bounds on features extracted from an image patch analysis of sunspot continuum and magnetogram images.
We study cosmological evolution after inflation in models with non-minimal derivative coupling to gravity. The background dynamics is solved and particle production associated with rapidly oscillating Hubble parameter is studied in detail. In addition, production of gravitons through the non-minimal derivative coupling with the inflaton is studied. We also find that the sound speed squared of the scalar perturbation oscillates between positive and negative values when the non-minimal derivative coupling dominates over the minimal kinetic term. This may lead to an instability of this model. We point out that the particle production rates are the same as those in the Einstein gravity with the minimal kinetic term, if we require the sound speed squared is positive definite.
We consider dark matter in a minimal extension of the Standard Model (SM) which breaks electroweak symmetry dynamically and leads to a complete unification of the SM and technicolor coupling constants. The unification scale is determined to be $M_{\rm U} \approx 2.2 \times 10^{15}$ GeV and the unified coupling $\alpha_{\rm U} \approx 0.0304$. Moreover, unification strongly suggest that the technicolor sector of the model must become strong at the scale of ${\cal O}$(TeV). The model also contains a tightly constrained sector of mixing neutral fields stabilized by a discrete symmetry. We find the lightest of these states can be DM with a mass in the range $m_{\rm DM} \approx 30-800$ GeV. We find a large set of parameters that satisfy all available constraints from colliders and from dark matter search experiments. However, most of the available parameter space is within the reach of the next generation of DM search experiments. The model is also sensitive to a modest improvement in the measurement of the precision electroweak parameters.
We study a model of accidental inflation in type IIB string theory where inflation occurs near the inflection point of a small K\"ahler modulus. A racetrack structure helps to alleviate the known concern that string-loop corrections may spoil K\"ahler Moduli Inflation unless having a significant suppression via the string coupling or a special brane setup. Also, the hierarchy of gauge group ranks required for the separation between moduli stabilization and inflationary dynamics is relaxed. The relaxation becomes more significant when we use the recently proposed D-term generated racetrack model.
The AMS-02 collaboration has just released the cosmic antiproton to proton ratio $\bar{p}/p$ with a high precision up to $\sim 450$ GeV. In this work, we calculate the secondary antiprotons generated by cosmic ray interactions with the interstellar medium taking into account the uncertainties from the cosmic ray propagation. The $\bar{p}/p$ ratio predicted by these processes shows some tension with the AMS-02 data in some regions of propagation parameters, but the excess is not significant. We then try to derive upper bounds on the dark matter annihilation cross section from the $\bar{p}/p$ data or signal regions favored by the data. It is shown that the constraint derived by the AMS-02 data is similar to that from Fermi-LAT observations of dwarf galaxies. The signal region for dark matter is usually required $m_\chi \sim O(10)$ TeV and $\left<\sigma v\right>\sim\mathcal{O}(10^{-23})~\cm^3~\sec^{-1}$.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
Blazars are radio-loud active galactic nuclei well known for their non thermal emission spanning a wide range of frequencies. The Roma-BZCAT is, to date, the most comprehensive list of these sources. We performed the cross-match of several catalogs obtained from recent surveys at different frequencies to search for new blazars. We cross-matched the 1$^{st}$ Swift-XRT Point Source catalog with the spectroscopic sample of the 9$^{th}$ Data Release of the Sloan Digital Sky Survey. Then, we performed further cross-matches with the catalogs corresponding to the Faint Images of the Radio Sky at Twenty cm survey and to the AllWISE Data release, focusing on sources with infrared colors similar to those of confirmed $\gamma$-ray blazars included in the Second Fermi-LAT catalog. As a result, we obtained a preliminary list of objects with all the elements needed for a proper blazar classification according to the prescriptions of the Roma-BZCAT. We carefully investigated additional properties such as their morphology and the slope of their spectral energy distribution in the radio domain, the features shown in their optical spectrum, and the luminosity in the soft X rays to exclude generic active galactic nuclei and focus on authentic blazar-like sources. At the end of our screening we obtained a list of 15 objects with firmly established blazar properties.
We report results from deep observations (750 ks) of Tycho's supernova remnant (SNR) with NuSTAR. Using these data, we produce narrow-band images over several energy bands to identify the regions producing the hardest X-rays and to search for radioactive decay line emission from 44Ti. We find that the hardest (>10 keV) X-rays are concentrated in the southwest of Tycho, where recent Chandra observations have revealed high emissivity "stripes" associated with particles accelerated to the knee of the cosmic-ray spectrum. We do not find evidence of 44Ti, and we set tight limits on its presence which exclude the reported Swift/BAT and INTEGRAL detections and correspond to an upper-limit 44Ti mass of M44 < 8.4e-5 Msun for a distance of 2.3 kpc. We perform spatially resolved spectroscopic analysis of sixty-six regions across Tycho. We map the best-fit rolloff frequency of the hard X-ray spectra, and we compare these results to measurements of the shock expansion and ambient density. We find that the highest energy electrons are accelerated at the lowest densities and in the fastest shocks, with a steep dependence of the roll-off frequency with shock velocity. Such a dependence is predicted by models where the maximum energy of accelerated electrons is limited by the age of the SNR rather than by synchrotron losses, but this scenario requires far lower magnetic field strengths than those derived from observations in Tycho. One way to reconcile these discrepant findings is through shock obliquity effects, and future observational work is necessary to explore the role of obliquity in the particle acceleration process.
We present new near-infrared spectroscopic observations of the outer edges of the young stellar cluster around the supermassive black hole at the Galactic center. The observations show a break in the surface-density profile of young stars at approximately 13 arcsec (0.52 pc). These observations spectroscopically confirm previous suggestions of a break based on photometry. Using Gemini North's Near-Infrared Integral Field Spectrometer (NIFS) we are able to detect and separate early- and late-type stars with a 75% completeness at Ks = 15.5. We sample a region with radii between 7" to 23" (0.28 pc to 0.92 pc) from Sgr A*, and present new spectral classifications of 144 stars brighter than Ks = 15.5, where 140 stars are late-type (> 1 Gyr) and only four stars are early-type (young, 4-6 Myr). A broken power-law fit of the early-type surface-density matches well with our data and previously published values. The projected surface-density of late-type stars is also measured and found to be consistent with previous results. We find that the observed early-type surface-density profile is inconsistent with the theory of the young stars originating from a tightly bound infalling cluster, as no significant trail of young stars is found at radii above 13". We also note that either a simple disk instability criterion or a cloud-cloud collision could explain the location of the outer edge, though we lack information to make conclusive remarks on either alternative. If this break in surface-density represents an edge to the young stellar cluster it would set an important scale for the most recent episode of star formation at the Galactic center.
We present 0.95-1.80 $\mu$m spectroscopy of the $\sim$12-27 $M_{\rm Jup}$ companion orbiting the faint ($R$$\sim$13.6), young ($\sim$120 Myr) M-dwarf 2MASS J01225093--2439505 ("2M0122--2439 B") at 1.5 arcsecond separation (50 AU). Our coronagraphic long-slit spectroscopy was obtained with the new high contrast imaging platform VLT-SPHERE during Science Verification. The unique long-slit capability of SPHERE enables spectral resolution an order of magnitude higher than other extreme AO exoplanet imaging instruments. With a low mass, cool temperature, and very red colors, 2M0122-2439 B occupies a particularly important region of the substellar color-magnitude diagram by bridging the warm directly imaged hot planets with late-M/early-L spectral types (e.g. $\beta$ Pic b and ROXs 42Bb) and the cooler, dusty objects near the L/T transition (e.g. HR 8799bcde and 2MASS 1207b). We fit BT-Settl atmospheric models to our $R$$\approx$350 spectrum and find $T_{\rm eff}$=1600$\pm$100 K and $\log(g)$=4.5$\pm$0.5 dex. Visual analysis of our 2M0122-2439 B spectrum suggests a spectral type L3-L4, and we resolve shallow $J$-band alkali lines, confirming its low gravity and youth. Specifically, we use the Allers & Liu (2013) spectral indices to quantitatively measure the strength of the FeH, VO, KI, spectral features, as well as the overall $H$-band shape. Using these indices, along with the visual spectral type analysis, we classify 2M0122-2439 B as an intermediate gravity (INT-G) object with spectral type L3.7$\pm$1.0.
We report on simultaneous observations of the magnetar SGR J1745-2900 at frequencies $\nu = 2.57$ to $225\,\rm{GHz}$ using the Nancay 94-m equivalent, Effelsberg 100-m, and IRAM 30-m radio telescopes. We detect SGR J1745-2900 up to 225 GHz, the highest radio frequency detection of pulsed emission from a neutron star to date. Strong single pulses are also observed from 4.85 up to 154 GHz. At the millimetre band we see significant flux density and spectral index variabilities on time scales of tens of minutes, plus variability between days at all frequencies. Additionally, SGR J1745-2900 was observed at a different epoch at frequencies 296 to 472 GHz using the APEX 12-m radio telescope, with no detections. Over the period MJD 56859.83-56862.93 the fitted spectrum yields a spectral index of $\left<\alpha\right> = -0.4 \pm 0.1$ for a reference flux density $\left< S_{154} \right> = 1.1 \pm 0.2\rm{\,mJy}$ (with $S_{\nu} \propto {\nu}^{\alpha})$, a flat spectrum alike those of the other radio loud magnetars. These results show that strongly magnetized neutron stars can be effective radio emitters at frequencies notably higher to what was previously known and that pulsar searches in the Galactic Centre are possible in the millimetre band.
We present the first numerical simulations in coupled dark energy cosmologies with high enough resolution to investigate the effects of the coupling on galactic and sub-galactic scales. We choose two constant couplings and a time-varying coupling function and we run simulations of three Milky-Way-size halos ($\sim$10$^{12}$M$_{\odot}$), a lower mass halo (6$\times$10$^{11}$M$_{\odot}$) and a dwarf galaxy halo (5$\times$10$^{9}$M$_{\odot}$). We resolve each halo with several millions dark matter particles. On all scales the coupling causes lower halo concentrations and a reduced number of substructures with respect to LCDM. We show that the reduced concentrations are not due to different formation times, but they are related to the extra terms that appear in the equations describing the gravitational dynamics. On the scale of the Milky Way satellites, we show that the lower concentrations can help in reconciling observed and simulated rotation curves, but the coupling values necessary to have a significant difference from LCDM are outside the current observational constraints. On the other hand, if other modifications to the standard model allowing a higher coupling (e.g. massive neutrinos) are considered, coupled dark energy can become an interesting scenario to alleviate the small-scale issues of the LCDM model.
Bayesian inference is often used in cosmology and astrophysics to derive constraints on model parameters from observations. This approach relies on the ability to compute the likelihood of the data given a choice of model parameters. In many practical situations, the likelihood function may however be unavailable or intractable due to non-gaussian errors, non-linear measurements processes, or complex data formats such as catalogs and maps. In these cases, the simulation of mock data sets can often be made through forward modeling. We discuss how Approximate Bayesian Computation (ABC) can be used in these cases to derive an approximation to the posterior constraints using simulated data sets. This technique relies on the sampling of the parameter set, a distance metric to quantify the difference between the observation and the simulations and summary statistics to compress the information in the data. We first review the principles of ABC and discuss its implementation using a Population Monte-Carlo (PMC) algorithm. We test the performance of the implementation using a Gaussian toy model. We then apply the ABC technique to the practical case of the calibration of image simulations for wide field cosmological surveys. We find that the ABC analysis is able to provide reliable parameter constraints for this problem and is therefore a promising technique for other applications in cosmology and astrophysics. Our implementation of the ABC PMC method is made available via a public code release.
We present the results of the decade-long M31 observation from the Wendelstein Calar Alto Pixellensing Project (WeCAPP). WeCAPP has monitored M31 from 1997 till 2008 in both R- and I-filters, thus provides the longest baseline of all M31 microlensing surveys. The data are analyzed with the difference imaging analysis, which is most suitable to study variability in crowded stellar fields. We extracted light curves based on each pixel, and devised selection criteria that are optimized to identify microlensing events. This leads to 10 new events, and sums up to a total of 12 microlensing events from WeCAPP, for which we derive their timescales, flux excesses, and colors from their light curves. The color of the lensed stars fall between (R-I) = 0.56 to 1.36, with a median of 1.0 mag, in agreement with our expectation that the sources are most likely bright, red stars at post main-sequence stage. The event FWHM timescales range from 0.5 to 14 days, with a median of 3 days, in good agreement with predictions based on the model of Riffeser et al. (2006).
We have newly identified a substantial number of type 1 active galactic nuclei (AGN) featuring weak broad-line regions (BLRs) at z < 0.2 from detailed analysis of galaxy spectra in the Sloan Digital Sky Survey Data Release 7. These objects predominantly show a stellar continuum but also a broad H-alpha emission line, indicating the presence of a low-luminosity AGN oriented so that we are viewing the central engine directly without significant obscuration. These accreting black holes have previously eluded detection due to their weak nature. The new BLR AGNs we found increased the number of known type 1 AGNs by 49%. Some of these new BLR AGNs were detected at the Chandra X-ray Observatory, and their X-ray properties confirm that they are indeed type 1 AGN. Based on our new and more complete catalogue of type 1 AGNs, we derived the type 1 fraction of AGNs as a function of [OIII] 5007 emission luminosity and explored the possible dilution effect on the obscured AGN due to star-formation. The new type 1 AGN fraction shows much more complex behavior with respect to black hole mass and bolometric luminosity than suggested by the existing receding torus model. The type 1 AGN fraction is sensitive to both of these factors, and there seems to be a sweet spot (ridge) in the diagram of black hole mass and bolometric luminosity. Furthermore, we present a hint that the Eddington ratio plays a role in determining the opening angles.
We reason that, without physical fine-tuning, neither the supermassive black holes (SMBHs) nor the stellar bulges can self-regulate or inter-regulate by driving away already fallen cold gas to produce the observed correlation between them. We suggest an alternative scenario where the observed mass ratios of the SMBHs to bulges reflect the angular momentum distribution of infallen gas such that the mass reaching the stable accretion disc is a small fraction of that reaching the bulge region, averaged over the cosmological time scales. We test this scenario using high resolution, large-scale cosmological hydrodynamic simulations (without AGN feedback), assuming the angular momentum distribution of gas landing in the bulge region to yield a Mestel disc that is supported by independent simulations resolving the Bondi radii of SMBHs. A mass ratio of $0.1-0.3\%$ between the very low angular momentum gas that free-falls to the sub-parsec region to accrete to the SMBH and the overall star formation rate is found. This ratio is found to increase with increasing redshift to within a factor of $\sim 2$, suggesting that the SMBH to bulge ratio is nearly redshift independent, with a modest increase with redshift, a testable prediction. Furthermore, the duty cycle of active galactic nuclei (AGN) with high Eddington ratios is expected to increase significantly with redshift. Finally, while SMBHs and bulges are found to coevolve on $\sim 30-150$Myr time scales or longer, there is indication that, on shorer time scales, the SMBH accretion rate and star formation may be less correlated.
Supergiant Fast X-ray Transients (SFXTs) are HMXBs with OB supergiant companions. I review the results of the Swift SFXT Project, which since 2007 has been exploiting Swift's capabilities in a systematic study of SFXTs and supergiant X-ray binaries (SGXBs) by combining follow-ups of outbursts, when detailed broad-band spectroscopy is possible, with long-term monitoring campaigns, when the out-of-outburst fainter states can be observed. This strategy has led us to measure their duty cycles as a function of luminosity, to extract their differential luminosity distributions in the soft X-ray domain, and to compare, with unprecedented detail, the X-ray variability in these different classes of sources. I also discuss the "seventh year crisis", the challenges that the recent Swift observations are making to the prevailing models attempting to explain the SFXT behaviour.
We propose a method to generate `genetically-modified' (GM) initial conditions for high-resolution simulations of galaxy formation in a cosmological context. Building on the Hoffman-Ribak algorithm, we start from a reference simulation with fully random initial conditions, then make controlled changes to specific properties of a single halo (such as its mass and merger history). The algorithm demonstrably makes minimal changes to other properties of the halo and its environment, allowing us to isolate the impact of a given modification. As a significant improvement over previous work, we are able to calculate the abundance of the resulting objects relative to the $\Lambda$CDM reference cosmology. Our approach can be applied to a wide range of cosmic structures and epochs; here we study two problems as a proof-of-concept. First, we investigate the change in density profile and concentration as the collapse time of three individual halos are varied at fixed final mass, showing good agreement with previous statistical studies using large simulation suites. Second, we modify the $z=0$ mass of halos to show that our theoretical abundance calculations correctly recover the halo mass function. The results demonstrate that the technique is robust, opening the way to controlled experiments in galaxy formation using hydrodynamic zoom simulations.
We have studied in detail the 0.15-15 GHz radio spectrum of the gamma-ray binary LS 5039 to look for a possible turnover and absorption mechanisms at low frequencies, and to constrain the physical properties of its emission. We have analysed two archival VLA monitorings, all the available archival GMRT data and a coordinated quasi-simultaneous observational campaign conducted in 2013 with GMRT and WSRT. The data show that the radio emission of LS 5039 is persistent on day, week and year timescales, with a variability $\lesssim 25~\%$ at all frequencies, and no signature of orbital modulation. The obtained spectra reveal a power-law shape with a curvature below 5 GHz and a turnover at $\sim0.5$ GHz, which can be reproduced by a one-zone model with synchrotron self-absorption plus Razin effect. We obtain a coherent picture for a size of the emitting region of $\sim0.85~\mathrm{mas}$, setting a magnetic field of $B\sim20~\mathrm{mG}$, an electron density of $n_{\rm e}\sim4\times10^5~{\rm cm^{-3}}$ and a mass-loss rate of $\dot M\sim5\times10^{-8}~{\rm M_{\odot} yr^{-1}}$. These values imply a significant mixing of the stellar wind with the relativistic plasma outflow from the compact companion. At particular epochs the Razin effect is negligible, implying changes in the injection and the electron density or magnetic field. The Razin effect is reported for first time in a gamma-ray binary, giving further support to the young non-accreting pulsar scenario.
We present the results of a search for cold gas at high redshift along QSO lines-of-sight carried out without any a priori assumption on the neutral atomic-hydrogen (HI) content of the absorbers. To do this, we systematically looked for neutral-carbon (CI) 1560,1656 transition lines in low-resolution QSO spectra from the SDSS database. We built up a sample of 66 CI absorbers with redshifts 1.5<z<3.1 and equivalent widths 0.1<W_r(1560)<1.7 A. The completeness limit of our survey is W_r,lim(1560)~0.4 A. CI systems stronger than that are more than one hundred-times rarer than DLAs at z_abs=2.5. The number of CI systems per unit redshift increases significantly below z=2. We suggest that the CI absorbers are closely related to the process of star formation and the production of dust in galaxies. We derive the HI content of the CI systems and find that a majority of them are sub-DLAs with N(HI)~10^20 atoms cm^-2. The dust content of these absorbers is yet significant as seen from the redder optical colours of the background QSOs and their reddened SEDs. The overall N(HI) distribution of CI systems is relatively flat however. As a consequence, among the CI systems classifying as DLAs there is a probable excess of strong DLAs with log N(HI)>21 compared to systematic DLA surveys. We study empirical relations between W_r(CI), N(HI), E(B-V) and the strength of the 2175 A extinction feature, the latter being detected in about 30% of the CI absorbers. We show that the 2175 A feature is weak compared to Galactic lines-of-sight exhibiting the same amount of reddening. This is probably the consequence of current or past star formation in the vicinity of the CI systems. We also find that the strongest CI systems tend to have the largest amounts of dust and that the metallicity of the gas and its molecular fraction is likely to be high in a large number of cases.
We propose a new method to estimate the photometric redshift of galaxies by using the full galaxy image in each measured band. This method draws from the latest techniques and advances in machine learning, in particular Deep Neural Networks. We pass the entire multi-band galaxy image into the machine learning architecture to obtain a redshift estimate that is competitive with the best existing standard machine learning techniques. The standard techniques estimate redshifts using post-processed features, such as magnitudes and colours, which are extracted from the galaxy images and are deemed to be salient by the user. This new method removes the user from the photometric redshift estimation pipeline. However we do note that Deep Neural Networks require many orders of magnitude more computing resources than standard machine learning architectures.
We present stellar metallicities in Leo I, Leo II, IC 1613, and Phoenix dwarf
galaxies derived from medium (F390M) and broad (F555W, F814W) band photometry
using the Wide Field Camera 3 (WFC3) instrument aboard the Hubble Space
Telescope. We measured metallicity distribution functions (MDFs) in two ways,
1) matching stars to isochrones in color-color diagrams, and 2) solving for the
best linear combination of synthetic populations to match the observed
color-color diagram. The synthetic technique reduces the effect of photometric
scatter, and produces MDFs 30-50 % narrower than the MDFs produced from
individually matched stars. We fit the synthetic and individual MDFs to
analytical chemical evolution models (CEM) to quantify the enrichment and the
effect of gas flows within the galaxies. Additionally, we measure stellar
metallicity gradients in Leo I and II. For IC 1613 and Phoenix our data do not
have the radial extent to confirm a metallicity gradient for either galaxy.
We find the MDF of Leo I (dwarf spheroidal) to be very peaked with a steep
metal rich cutoff and an extended metal poor tail, while Leo II (dwarf
spheroidal), Phoenix (dwarf transition) and IC 1613 (dwarf irregular) have
wider, less peaked MDFs than Leo I. A simple CEM is not the best fit for any of
our galaxies, therefore we also fit the `Best Accretion Model' of Lynden-Bell
1975. For Leo II, IC 1613 and Phoenix we find similar accretion parameters for
the CEM, even though they all have different effective yields, masses, star
formation histories and morphologies. We suggest that the dynamical history of
a galaxy is reflected in the MDF, where broad MDFs are seen in galaxies that
have chemically evolved in relative isolation and narrowly peaked MDFs are seen
in galaxies that have experienced more complicated dynamical interactions
concurrent with their chemical evolution.
We present a new photometric catalog of 326 candidate globular clusters (GCs) in the nearby spiral galaxy M101, selected from B, V, and I Hubble Space Telescope Advanced Camera for Surveys images. The luminosity function (LF) of these clusters has an unusually large number of faint sources compared with GCLFs in many other spiral galaxies. Accordingly, we separate and compare the properties of "bright" (M_V < -6.5) versus "faint" (M_V > -6.5; one magnitude fainter than the expected GC peak) clusters within our sample. The LF of the bright clusters is well fit by a peaked distribution similar to those observed in the Milky Way (MW) and other galaxies. These bright clusters also have similar size (r_{eff}) and spatial distributions as MW GCs. The LF of the faint clusters, on the other hand, is well described by a power law, dN(L_V)/dL_V proportional to L_V^alpha with alpha = -2.6 plus or minus 0.3, similar to those observed for young and intermediate-age cluster systems in star forming galaxies. We find that the faint clusters have larger typical r_{eff} than the bright clusters, and have a flatter surface density profile, being more evenly distributed, as we would expect for clusters associated with the disk. We use the shape of the LF and predictions for mass-loss driven by two-body relaxation to constrain the ages of the faint clusters. Our results are consistent with two populations of old star clusters in M101: a bright population of halo clusters and a fainter, possibly younger, population of old disk clusters.
We provide predictions on small-scale cosmological density power spectrum from supernova lensing dispersion. Parameterizing the primordial power spectrum with running $\alpha$ and running of running $\beta$ of the spectral index, we exclude large positive $\alpha$ and $\beta$ parameters which induce too large lensing dispersions over current observational upper bound. We ran cosmological N-body simulations of collisionless dark matter particles to investigate non-linear evolution of the primordial power spectrum with positive running parameters. The initial small-scale enhancement of the power spectrum is largely erased when entering into the non-linear regime. For example, even if the linear power spectrum at $k>10h {\rm Mpc}^{-1}$ is enhanced by $1-2$ orders of magnitude, the enhancement much decreases to a factor of $2-3$ at late time ($z \leq 1.5$). Therefore, the lensing dispersion induced by the dark matter fluctuations weakly constrains the running parameters. When including baryon-cooling effects (which strongly enhance the small-scale clustering), the constraint is comparable or tighter than the PLANCK constraint, depending on the UV cut-off. Further investigations of the non-linear matter spectrum with baryonic processes is needed to reach a firm constraint.
We present a comprehensive study of interstellar X-ray extinction using the extensive Chandra supernova remnant archive and use our results to refine the empirical relation between the hydrogen column density and optical extinction. In our analysis, we make use of the large, uniform data sample to assess various systematic uncertainties in the measurement of the interstellar X-ray absorption. Specifically, we address systematic uncertainties that originate from (i) the emission models used to fit supernova remnant spectra, (ii) the spatial variations within individual remnants, (iii) the physical conditions of the remnant such as composition, temperature, and non-equilibrium regions, and (iv) the model used for the absorption of X-rays in the interstellar medium. Using a Bayesian framework to quantify these systematic uncertainties, and combining the resulting hydrogen column density measurements with the measurements of optical extinction toward the same remnants, we find the empirical relation NH = (2.87+/-0.12) x 10^21 AV cm^(-2), which is significantly higher than the previous measurements.
A hot, dusty torus located around the outer edge of the broad-line region of AGNs is a fundamental ingredient in unified AGN models. While the existence of circumnuclear dust around AGNs at pc-scale radii is now widely accepted, questions about the origin, evolution and long-term stability of these dust tori remain unsettled.\\ We used reverberation mapping of the hot circumnuclear dust in the Seyfert 1 galaxy NGC 4151, to monitor its temperature and reverberation lag as a function of the varying accretion disk brightness. We carried out multiband, multiepoch photometric observations of the nucleus of NGC 4151 in the z,Y,J,H, and K bands for 29 epochs from 2010 January to 2014 June, supported by new near-infrared and optical spectroscopic observations, and archived WISE data.\\ We see no signatures of dust destruction due to sublimation in our data, since they show no increase in the hot dust reverberation delay directly correlated with substantial accretion disk flux increases in the observed period. Instead, we find that the hot dust in NGC 4151 appears to merely heat up, and the hot dust temperature closely tracks the accretion disk luminosity variations. We find indications of a decreased reverberation delay within the observed period from t = 42.5 +/- 4.0 days in 2010 to t = 29.6 +/- 1.7 days in 2013-2014. Such a varying reverberation radius on longer timescales would explain the intrinsic scatter observed in the radius-luminosity relation of dust around AGNs.\\ Our observations rule out that a second, larger dust component within a 100-light-day radius from the source contributes significantly to the observed near-infrared flux in this galaxy.
We verified the off-axis jet model of X-ray flashes (XRFs) and examined a discovery of off-axis orphan gamma-ray burst (GRBs) afterglows. The XRF sample was selected on the basis of the following three factors: (1) a constraint on the lower peak energy of the prompt spectrum $E^{src}_{obs}$, (2) redshift measurements, and (3) multi-color observations of an earlier (or brightening) phase. XRF020903 was the only sample selected basis of these criteria. A complete optical multi-color afterglow light curve of XRF020903 obtained from archived data and photometric results in literature showed an achromatic brightening around 0.7 days. An off-axis jet model with a large observing angle (0.21 rad, which is twice the jet opening half-angle, $\theta_{jet}$) can naturally describe the achromatic brightening and the prompt X-ray spectral properties. This result indicates the existence of off-axis orphan GRB afterglow light curves. Events with a larger viewing angle ($>\sim2\theta_{jet}$) could be discovered using an 8-m class telescope with wide field imagers such as Subaru Hyper-Suprime-Cam and the Large Synoptic Survey Telescope.
Aristotelian electrodynamics (AE) describes the regime of a plasma with a very strong electric field that is not shorted out, with charge current determined completely by pair production and the balance of Lorentz 4-force against curvature radiation reaction. Here it is shown how the principal null directions and associated eigenvalues of the field tensor govern AE, and how force-free electrodynamics arises smoothly from AE when the eigenvalues (and therefore the electric field in some frame) vanish. A criterion for validity of AE and force-free electrodynamics is proposed in terms of a pair of "field curvature scalars" formed from the first derivative of the principal null directions.
We perform a detailed comparison of the phase-space density traced by the particle distribution in Gadget simulations to the result obtained with a spherical Vlasov solver using the splitting algorithm. The systems considered are apodized H\'enon spheres with two values of the virial ratio, R ~ 0.1 and 0.5. After checking that spherical symmetry is well preserved by the N-body simulations, visual and quantitative comparisons are performed. In particular we introduce new statistics, correlators and entropic estimators, based on the likelihood of whether N-body simulations actually trace randomly the Vlasov phase-space density. When taking into account the limits of both the N-body and the Vlasov codes, namely collective effects due to the particle shot noise in the first case and diffusion and possible nonlinear instabilities due to finite resolution of the phase-space grid in the second case, we find a spectacular agreement between both methods, even in regions of phase-space where nontrivial physical instabilities develop. However, in the colder case, R=0.1, it was not possible to prove actual numerical convergence of the N-body results after a number of dynamical times, even with N=10$^8$ particles.
We present a novel method for particle splitting in smoothed particle hydrodynamics simulations. Our method utilizes the Voronoi diagram for a given particle set to determine the position of fine daughter particles. We perform several test simulations to compare our method with a conventional splitting method in which the daughter particles are placed isotropically over the local smoothing length. We show that, with our method, the density deviation after splitting is reduced by a factor of about two compared with the conventional method. Splitting would smooth out the anisotropic density structure if the daughters are distributed isotropically, but our scheme allows the daughter particles to trace the original density distribution with length scales of the mean separation of their parent. We apply the particle splitting to simulations of the primordial gas cloud collapse. The thermal evolution is accurately followed to the hydrogen number density of 10^12 /cc. With the effective mass resolution of ~10^-4 Msun after the multi-step particle splitting, the protostellar disk structure is well resolved. We conclude that the method offers an efficient way to simulate the evolution of an interstellar gas and the formation of stars.
The Atacama Large Millimeter/submillimeter Array (ALMA), an international partnership of Europe, North America and East Asia in cooperation with the Republic of Chile, is the largest astronomical project in existence. While ALMA's capabilities are ramping up, Early Science observations have started. The ALMA Archive is at the center of the operations of the telescope array and is designed to manage the 200 TB of data that will be taken each year, once the observatory is in full operations. We briefly describe design principles. The second part of this paper focuses on how astronomy is likely to evolve as the amount and complexity of data taken grows. We argue that in the future observatories will compete for astronomers to work with their data, that observatories will have to reorient themselves to from providing good data only to providing an excellent end-to-end user-experience with all its implications, that science-grade data-reduction pipelines will become an integral part of the design of a new observatory or instrument and that all this evolution will have a deep impact on how astronomers will do science. We show how ALMA's design principles are in line with this paradigm.
Galaxy-scale outflows, which are thought to provide the link connecting the central black hole to its host galaxy, are now starting to be observed. However, the physical origin of the mechanism driving the observed outflows, whether due to energy-driving or radiation-driving, is still debated; and in some cases, it is not clear whether the central source is an active galactic nucleus (AGN) or a nuclear starburst. Here we study the role of radiation pressure on dust in driving galactic-scale AGN outflows, and analyse the dynamics of the outflowing shell as a function of the underlying physical parameters. We show that high-velocity outflows ($\gtrsim$1000 km/s) with large momentum flux ($\gtrsim 10 L/c$) can be obtained, by taking into account the effects of radiation trapping. In particular, the high observed values of the momentum boosts can be reproduced, provided that the shell is initially optically thick to the reprocessed infrared radiation. Alternatively, the inferred measurements of the momentum flux may be significantly biased by AGN variability. In this context, the observations of powerful outflows on kiloparsec scales, with no or weak signs of ongoing nuclear activity at the present time, could be re-interpreted as relics of past AGN episodes.
Strange stars are one of the possible compact stellar objects that can be formed after a supernova collapse. We consider a model of strange star having an inner core in the color-flavor locked phase surmounted by a crystalline color superconducting layer. These two phases constitute the {\it quarksphere}, which we assume to be the largest and heaviest part of the strange star. The next layer consists of standard nuclear matter forming a ionic crust, hovering on the top of the quarksphere and prevented from falling by a strong dipolar electric field. The dipolar electric field arises because quark matter is confined in the quarksphere by the strong interaction, but electrons can leak outside forming a few hundreds Fermi thick electron layer separating the ionic crust from the underlying quark matter. The ionic matter and the crystalline color superconducting matter constitute two electromagnetically coupled crust layers. We study the torsional oscillations of these two layers. Remarkably, we find that if a fraction larger than $10^{-4}$ of the energy of a Vela-like glitch is conveyed to a torsional oscillation, the ionic crust will likely break. The reason is that the very rigid and heavy crystalline color superconducting crust layer will absorb only a small fraction of the glitch energy, leading to a large amplitude torsional oscillation of the ionic crust.
We investigate the impact of high optical depth on the HI surface density (Sigma_HI) saturation observed in the Perseus molecular cloud. We use Arecibo HI emission and absorption measurements obtained toward 26 radio continuum sources to derive the spin temperature and optical depth of individual HI components along each line of sight. The derived properties are used to estimate the correction for high optical depth. We examine two different methods for the correction, Gaussian decomposition and isothermal approximation methods, and find that they are consistent as having the maximum correction factor of ~1.2 likely due to the relatively low optical depth and insignificant contribution from the diffuse radio continuum emission for Perseus. We apply the correction to the HI column density image derived in the optically thin approximation on a pixel-by-pixel basis, and find that the total HI mass increases by only ~10%. Using the corrected HI column density image and far-infrared data from the Improved Reprocessing of the IRAS Survey, we then derive the H2 column density across the cloud on ~0.4 pc scales. For five dark and star-forming regions in Perseus, the HI surface density is uniform with Sigma_HI ~ 7-9 solar mass/pc2, in agreement with the minimum HI surface density required for shielding H2 against photodissociation. As a result, Sigma_H2/Sigma_HI and Sigma_HI+Sigma_H2 show a remarkably tight relation. Our results are consistent with predictions for H2 formation in steady state and chemical equilibrium, and suggest that H2 formation is mainly responsible for the HI saturation in Perseus. We also compare the optically thick HI with the observed "CO-dark" gas, and find that the optically thick HI only accounts for ~20% of the "CO-dark" gas in Perseus.
We present a fast and user friendly exoplanet transit light curve modelling package PyTransit, implementing optimised versions of the Gimen\'ez and the Mandel & Agol transit models. The package offers an object-oriented Python interface to access the two models implemented natively in Fortran with OpenMP parallelisation. A partial OpenCL version of the quadratic Mandel-Agol model is also included for GPU-accelerated computations. The aim of PyTransit is to facilitate the analysis of photometric time series of exoplanet transits consisting of hundreds of thousands of datapoints, and of multi-passband transit light curves from spectrophotometric observations, as a part of a researcher's programming toolkit for building complex, problem-specific, analyses.
We have undertaken a Spitzer campaign to measure the IR structures and spectra of low-redshift 3CRR radio galaxies. The results show that the 3.6 - 160 micron infrared properties vary systematically with integrated source power, and so demonstrate that contemporary core activity is characteristic of the behaviour of sources over their lifetimes. IR synchrotron emission is seen from jets and hotspots in some cases. Thermal emission is found from a jet/gas interaction in NGC 7385. Most of the near-IR integrated colours of the low-redshift 3CRR radio galaxies are similar to those of passive galaxies, so that IR colours are poor indicators of radio activity.
We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the GRB originated from PopIII stars, while the detection would be easier if the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first phase of the Square Kilometer Array (SKA1-low). On the other hand, more standard GRBs are too dim to be detected even with ten times the sensitivity of SKA1-low, and only in the most optimistic case can a S/N larger than a few be reached at z>9.
We present a new numerical tool for solving the special relativistic ideal MHD equations that is based on the combination of the following three key features: (i) a one-step ADER discontinuous Galerkin (DG) scheme that allows for an arbitrary order of accuracy in both space and time, (ii) an a posteriori subcell finite volume limiter that is activated to avoid spurious oscillations at discontinuities without destroying the natural subcell resolution capabilities of the DG finite element framework and finally (iii) a space-time adaptive mesh refinement (AMR) framework with time-accurate local time-stepping. The divergence-free character of the magnetic field is instead taken into account through the so-called 'divergence-cleaning' approach. The convergence of the new scheme is verified up to 5th order in space and time and the results for a sample of significant numerical tests including shock tube problems, the RMHD rotor problem and the Orszag-Tang vortex system are shown. We also consider a simple case of the relativistic Kelvin-Helmholtz instability with a magnetic field, emphasizing the potential of the new method for studying turbulent RMHD flows. We discuss the advantages of our new approach when the equations of relativistic MHD need to be solved with high accuracy within various astrophysical systems.
We present the results of precision gamma-ray timing measurements of the binary millisecond pulsar PSR J2339$-$0533, an irradiating system of "redback" type, using data from the Fermi Large Area Telescope. We describe an optimized analysis method to determine a long-term phase-coherent timing solution spanning more than six years, including a measured eccentricity of the binary orbit and constraints on the proper motion of the system. A major result of this timing analysis is the discovery of an extreme variation of the nominal 4.6-hour orbital period $P_{\rm orb}$ over time, showing alternating epochs of decrease and increase. We inferred a cyclic modulation of $P_{\rm orb}$ with an approximate cycle duration of 4.2 years and a modulation amplitude of $\Delta P_{\rm orb}/ P_{\rm orb} = 2.3 \times 10^{-7}$. Considering different possible physical causes, the observed orbital-period modulation most likely results from a variable gravitational quadrupole moment of the companion star due to cyclic magnetic activity in its convective zone.
We consider the local dynamics of a realistic neutron star core, including composition gradients, superfluidity and thermal effects. The main focus is on the gravity g-modes, which are supported by composition stratification and thermal gradients. We derive the equations that govern this problem in full detail, paying particular attention to the input that needs to be provided through the equation of state and distinguishing between normal and superfluid regions. The analysis highlights a number of key issues that should be kept in mind whenever equation of state data is compiled from nuclear physics for use in neutron star calculations. We provide explicit results for a particular stellar model and a specific nucleonic equation of state, making use of cooling simulations to show how the local wave spectrum evolves as the star ages. Our results show that the composition gradient is effectively dominated by the muons whenever they are present. When the star cools below the superfluid transition, the support for g-modes at lower densities (where there are no muons) is entirely thermal. We confirm the recent suggestion that the g-modes in this region may be unstable, but our results indicate that this instability will be weak and would only be present for a brief period of the star's life. Our analysis accounts for the presence of thermal excitations encoded in entrainment between the entropy and the superfluid component. Finally, we discuss the complete spectrum, including the normal sound waves and, in superfluid regions, the second sound.
We determined accurate positions for 3000 of the "faint blue stars" in the PHL (Palomar-Haro-Luyten) and Ton/TonS catalogues. These were published from 1957 to 1962, and, aimed at finding new white dwarfs, provide approximate positions for about 10750 blue stellar objects. Some of these "stars" had become known as quasars, a type of objects unheard-of before 1963. We derived subarcsec positions from a comparison of published finding charts with images from the first-epoch Digitized Sky Survey. Numerous objects are now well known, but unfortunately neither their PHL or Ton numbers, nor their discoverers, are recognized in current databases. A comparison with modern radio, IR, UV and X-ray surveys leads us to suggest that the fraction of extragalactic objects in the PHL and Ton catalogues is at least 15 per cent. However, because we failed to locate the original PHL plates or finding charts, it may be impossible to correctly identify the remaining 7726 PHL objects.
We introduce a novel Earth-like planet surface temperature model (ESTM) for habitability studies based on the spatial-temporal distribution of planetary surface temperatures. The ESTM adopts a surface Energy Balance Model complemented by: radiative-convective atmospheric column calculations, a set of physically-based parameterizations of meridional transport, and descriptions of surface and cloud properties more refined than in standard EBMs. The parameterization is valid for rotating terrestrial planets with shallow atmospheres and moderate values of axis obliquity (epsilon >= 45^o). Comparison with a 3D model of atmospheric dynamics from the literature shows that the equator-to-pole temperature differences predicted by the two models agree within ~5K when the rotation rate, insolation, surface pressure and planet radius are varied in the intervals 0.5 <= Omega/Omega_o <= 2, 0.75 <= S/S_o <= 1.25, 0.3 <= p/(1 bar) <= 10, and 0.5 <= R/R_o <= 2, respectively. The ESTM has an extremely low computational cost and can be used when the planetary parameters are scarcely known (as for most exoplanets) and/or whenever many runs for different parameter configurations are needed. Model simulations of a test-case exoplanet (Kepler-62e) indicate that an uncertainty in surface pressure within the range expected for terrestrial planets may impact the mean temperature by ~60 K. Within the limits of validity of the ESTM, the impact of surface pressure is larger than that predicted by uncertainties in rotation rate, axis obliquity, and ocean fractions. We discuss the possibility of performing a statistical ranking of planetary habitability taking advantage of the flexibility of the ESTM.
Standing slow mode waves in hot flaring loops are exclusively observed in spectrometers and are used to diagnose the magnetic field strength and temperature of the loop structure. Due to the lack of spatial information, the longitudinal mode cannot be effectively identified. In this study, we simulate standing slow mode waves in flaring loops and compare the synthesized line emission properties with SUMER spectrographic and SDO/AIA imaging observations. We find that the emission intensity and line width oscillations are a quarter period out of phase with Doppler shift velocity both in time and spatial domain, which can be used to identify a standing slow mode wave from spectroscopic observations. However, the longitudinal overtones could be only measured with the assistance of imagers. We find emission intensity asymmetry in the positive and negative modulations, this is because the contribution function pertaining to the atomic emission process responds differently to positive and negative temperature variations. One may detect \textbf{half} periodicity close to the loop apex, where emission intensity modulation is relatively small. The line-of-sight projection affects the observation of Doppler shift significantly. A more accurate estimate of the amplitude of velocity perturbation is obtained by de-projecting the Doppler shift by a factor of $1-2\theta/\pi$ rather than the traditionally used $\cos\theta$. \textbf{If a loop is heated to the hotter wing, the intensity modulation could be overwhelmed by background emission, while the Doppler shift velocity could still be detected to a certain extent.
We present results of a search for giant radio galaxies (GRGs) with a projected largest linear size in excess of 1 Mpc. We designed a computational algorithm to identify contiguous emission regions, large and elongated enough to serve as GRG candidates, and applied it to the entire 1.4-GHz NRAO VLA Sky survey (NVSS). In a subsequent visual inspection of 1000 such regions we discovered 15 new GRGs, as well as many other candidate GRGs, some of them previously reported, for which no redshift was known. Our follow-up spectroscopy of 25 of the brighter hosts using two 2.1-m telescopes in Mexico, and four fainter hosts with the 10.4-m Gran Telescopio Canarias (GTC), yielded another 24 GRGs. We also obtained higher-resolution radio images with the Karl G. Jansky Very Large Array for GRG candidates with inconclusive radio structures in NVSS.
We investigate the influence of stellar migration caused by minor mergers (mass ratio from 1:70 to 1:8) on the radial distribution of chemical abundances in the disks of Milky Way-like galaxies during the last four Gyr. A GPU-based pure N-body tree-code model without hydrodynamics and star formation was used. We computed a large set of mergers with different initial satellite masses, positions, and orbital velocities. We find that there is no significant metallicity change at any radius of the primary galaxy in the case of accretion of a low-mass satellite of 10$^9$ M$_{\odot}$ (mass ratio 1:70) except for the special case of prograde satellite motion in the disk plane of the host galaxy. The accretion of a satellite of a mass $\gtrsim3\times10^9$ M$_{\odot}$ (mass ratio 1:23) results in an appreciable increase of the chemical abundances at galactocentric distances larger than $\sim10$ kpc. The radial abundance gradient flattens in the range of galactocentric distances from 5 to 15 kpc in the case of a merger with a satellite with a mass $\gtrsim3\times10^9$ M$_{\odot}$. There is no significant change in the abundance gradient slope in the outer disk (from $\sim15$ kpc up to 25 kpc) in any merger while the scatter in metallicities at a given radius significantly increases for most of the satellite's initial masses/positions compared to the case of an isolated galaxy. This argues against attributing the break (flattening) of the abundance gradient near the optical radius observed in the extended disks of Milky Way-like galaxies only to merger-induced stellar migration.
A new fitting methodology is presented which is equally well suited for the estimation of low-, medium-, and high-degree mode parameters from $m$-averaged solar oscillation power spectra of widely differing spectral resolution. This method, which we call the "Windowed, MuLTiple-Peak, averaged spectrum", or WMLTP Method, constructs a theoretical profile by convolving the weighted sum of the profiles of the modes appearing in the fitting box with the power spectrum of the window function of the observing run using weights from a leakage matrix that takes into account both observational and physical effects, such as the distortion of modes by solar latitudinal differential rotation. We demonstrate that the WMLTP Method makes substantial improvements in the inferences of the properties of the solar oscillations in comparison with a previous method that employed a single profile to represent each spectral peak. We also present an inversion for the internal solar structure which is based upon 6,366 modes that we have computed using the WMLTP method on the 66-day long 2010 SOHO/MDI Dynamics Run. To improve both the numerical stability and reliability of the inversion we developed a new procedure for the identification and correction of outliers in a frequency data set. We present evidence for a pronounced departure of the sound speed in the outer half of the solar convection zone and in the subsurface shear layer from the radial sound speed profile contained in Model~S of Christensen-Dalsgaard and his collaborators that existed in the rising phase of Solar Cycle~24 during mid-2010.
Context: The magnetic field is a crucial ingredient of neutron stars. It governs the physics of accretion and of the resulting high-energy emission in accreting pulsars. Studies of the cyclotron resonant scattering features (CRSFs) seen as absorption lines in the X-ray spectra of the pulsars permit direct measuremets of the field strength. Aims: From an analysis of a number of pointed observations with different instruments, the energy of CRSF, Ecyc, has recently been found to decay in Her X-1, which is one of the best-studied accreting pulsars. We present our analysis of a homogeneous and almost uninterrupted monitoring of the line energy with Swift/BAT. Methods: We analyzed the archival Swift/BAT observations of Her X-1 from 2005 to 2014. The data were used to measure the CRSF energy averaged over several months. Results: The analysis confirms the long-term decay of the line energy. The downward trend is highly significant and consistent with the trend measured with the pointed observations: dEcyc/dt ~-0.3 keV per year. Conclusions: The decay of Ecyc either indicates a local evolution of the magnetic field structure in the polar regions of the neutron star or a geometrical displacement of the line-forming region due to long-term changes in the structure of the X-ray emitting region. The shortness of the observed timescale of the decay, -Ecyc/(dEcyc/dt) ~ 100 yr, suggests that trend reversals and/or jumps of the line energy might be observed in the future.
Extra-galactic sources of photons have been used to constrain space-time quantum fluctuations in the Universe. In these proposals, the fundamental "fuzziness" of distance caused by space-time quantum fluctuations has been directly identified with fluctuations in optical paths. Phase-front corrugations deduced from these optical-path fluctuations are then applied to light from extra-galactic point sources, and used to constrain various models of quantum gravity. However, when a photon propagates in three spatial dimensions, it does not follow a specific ray, but rather samples a finite, three-dimensional region around that ray --- thereby averaging over space-time quantum fluctuations all through that region. We use a simple, random-walk type model to demonstrate that, once the appropriate wave optics is applied, the averaging of neighboring space-time fluctuations will cause much less distortion to the phase front. In our model, the extra suppression factor due to diffraction is the wave length in units of the Planck length, which is at least $10^{29}$ for astronomical observations.
The NIST database lists several Mn II lines that were observed in the laboratory but not classified. They cannot be used in spectrum synthesis because their atomic line data are unknown. These lines are concentrated in the 2380-2700 A interval. We aimed to assign energy levels and log gf values to these lines. Semi-empirical line data for Mn II computed by Kurucz were used to synthesize the ultraviolet spectrum of the slow-rotating, HgMn star HD 175640. The spectrum was compared with the high-resolution spectrum observed with the HST-STIS equipment. A UVES spectrum covering the 3050-10000 A region was also examined. We determined a total of 73 new energy levels, 58 from the STIS spectrum of HD 175640 and another 15 from the UVES spectrum. The new energy levels give rise to numerous new computed lines. We have identified more than 50% of the unclassified lines listed in the NIST database and have changed the assignement of another 24 lines. An abundance analysis of the star HD 175640, based on the comparison of observed and computed ultraviolet spectra in the 1250-3040 A interval, is the by-product of this study on Mn II.
Intensity mapping experiments survey the spectrum of diffuse line radiation rather than detect individual objects at high signal-to-noise. Spectral maps of unresolved atomic and molecular line radiation contain three-dimensional information about the density and environments of emitting gas, and efficiently probe cosmological volumes out to high redshift. Intensity mapping survey volumes also contain all other sources of radiation at the frequencies of interest. Continuum foregrounds are typically ~10^2-10^3 times brighter than the cosmological signal. The instrumental response to bright foregrounds will produce new spectral degrees of freedom that are not known in advance, nor necessarily spectrally smooth. The intrinsic spectra of foregrounds may also not be well-known in advance. We describe a general class of quadratic estimators to analyze data from single-dish intensity mapping experiments, and determine contaminated spectral modes from the data itself. The key attribute of foregrounds is not that they are spectrally smooth, but instead that they have fewer bright spectral degrees of freedom than the cosmological signal. Spurious correlations between the signal and foregrounds produce additional bias. Compensation for signal attenuation must estimate and correct this bias. A successful intensity mapping experiment will control instrumental systematics that spread variance into new modes, and it must observe a large enough volume that contaminant modes can be determined independently from the signal on scales of interest.
In this paper we describe a simple numerical approach which allows to study the structure of steady-state axisymmetric relativistic jets using one-dimensional time-dependent simulations. It is based on the fact that for narrow jets with v~c the steady-state equations of relativistic magnetohydrodynamics can be accurately approximated by the one-dimensional time-dependent equations after the substitution z=ct. Since only the time-dependent codes are now publicly available this is a valuable and efficient alternative to the development of a high-specialized code for the time-independent equations. The approach is also much cheaper and more robust compared to the relaxation method. We tested this technique against numerical and analytical solutions found in literature as well as solutions we obtained using the relaxation method and found it sufficiently accurate. In the process, we discovered the reason for the failure of the self-similar analytical model of the jet reconfinement in relatively flat atmospheres and elucidated the nature of radial oscillations of steady-state jets.
We investigated the response of the solar atmosphere to non-thermal electron beam heating using the radiative transfer and hydrodynamics modelling code RADYN. The temporal evolution of the parameters that describe the non-thermal electron energy distribution were derived from hard X-ray observations of a particular flare, and we compared the modelled and observed parameters. The evolution of the non-thermal electron beam parameters during the X1.5 solar flare on 2011 March 9 were obtained from analysis of RHESSI X-ray spectra. The RADYN flare model was allowed to evolve for 110 seconds, after which the electron beam heating was ended, and was then allowed to continue evolving for a further 300s. The modelled flare parameters were compared to the observed parameters determined from extreme-ultraviolet spectroscopy. The model produced a hotter and denser flare loop than that observed and also cooled more rapidly, suggesting that additional energy input in the decay phase of the flare is required. In the explosive evaporation phase a region of high-density cool material propagated upward through the corona. This material underwent a rapid increase in temperature as it was unable to radiate away all of the energy deposited across it by the non-thermal electron beam and via thermal conduction. A narrow and high-density ($n_{e} \le 10^{15}$ cm$^{-3}$) region at the base of the flare transition region was the source of optical line emission in the model atmosphere. The collision-stopping depth of electrons was calculated throughout the evolution of the flare, and it was found that the compression of the lower atmosphere may permit electrons to penetrate farther into a flaring atmosphere compared to a quiet Sun atmosphere.
The Kepler Mission has discovered thousands of planets with radii $<4 R_\oplus$, paving the way for the first statistical studies of the dynamics, formation, and evolution of these sub-Neptunes and super-Earths. Planetary masses are an important physical property for these studies, and yet the vast majority of Kepler planet candidates do not have theirs measured. A key concern for these studies is therefore how to map the measured radii to mass estimates for this Earth-to-Neptune size range where there are no Solar System analogs. Previous works have derived deterministic, one-to-one relationships between radius and mass. However, if these planets span a range of compositions as expected, then an intrinsic scatter about this relationship must exist in the population. Here we present the first probabilistic mass-radius relationship (M-R relation) evaluated within a Bayesian framework, which both quantifies this intrinsic dispersion and the uncertainties on the M-R relation parameters. We analyze how the details depend on the radius range of the sample, and on the method used to provide the mass measurements. Assuming that the M-R relation can be described as a power law with a dispersion that is constant and normally distributed, we find that $M/M_\oplus=2.7(R/R_\oplus)^{1.3}$ and a scatter in mass of $1.9M_\oplus$ is the "best-fit" probabilistic M-R relation for the sample of RV-measured transiting sub-Neptunes ($R_{pl}<4 R_\oplus$).
The variability of young stellar objects (YSO) changes their brightness and color preventing a proper classification in traditional color-color and color magnitude diagrams. We have explored the feasibility of the flux variation gradient (FVG) method for YSOs, using $H$ and $K$ band monitoring data of the star forming region RCW\,38 obtained at the University Observatory Bochum in Chile. Simultaneous multi-epoch flux measurements follow a linear relation $F_{H}=\alpha + \beta \cdot F_{K}$ for almost all YSOs with large variability amplitude. The slope $\beta$ gives the mean $HK$ color temperature $T_{var}$ of the varying component. Because $T_{var}$ is hotter than the dust sublimation temperature, we have tentatively assigned it to stellar variations. If the gradient does not meet the origin of the flux-flux diagram, an additional non- or less-varying component may be required. If the variability amplitude is larger at the shorter wavelength, e.g. $\alpha < 0$, this component is cooler than the star (e.g. a circumstellar disk); vice versa, if $\alpha > 0$, the component is hotter like a scattering halo or even a companion star. We here present examples of two YSOs, where the $HK$ FVG implies the presence of a circumstellar disk; this finding is consistent with additional data at $J$ and $L$. One YSO shows a clear $K$-band excess in the $JHK$ color-color diagram, while the significance of a $K$-excess in the other YSO depends on the measurement epoch. Disentangling the contributions of star and disk it turns out that the two YSOs have huge variability amplitudes ($\sim 3-5$\,mag). The $HK$ FVG analysis is a powerful complementary tool to analyze the varying components of YSOs and worth further exploration of monitoring data at other wavelengths.
I present a catalog of point source objects towards NGC 1333, resolving a wide variety of confusion about source names (and occasionally positions) in the literature. I incorporate data from optical to radio wavelengths, but focus most of the effort on being complete and accurate from J (1.25 um) to 24 um. The catalog encompasses 52 deg<RA<52.5 deg and 31 deg<Dec<31.6 deg. Cross-identifications include those from more than 25 papers and catalogs from 1994-2014, primarily those in wide use as origins of nomenclature. Gaps in our knowledge are identified, with the most important being a lack of spectroscopy for spectral types or even confirmation of youth and/or cluster membership. I fit a slope to the spectral energy distribution (SED) between 2 and 24 um for the members (and candidate members) to obtain an SED classification, and compare the resulting classes to those for the same sources in the literature, and for an SED fit between 2 and 8 um. While there are certainly differences, for the majority of the sources, there is good agreement.
The current detectors of gamma-ray emission have too poor resolution to determine whether this emission is produced in the jet or in the core, specially of low luminous, non-blazar AGNs (as radio galaxies). In recent works it has been found that the power released by events of turbulent fast magnetic reconnection in the core region of these sources is more than sufficient to reproduce the observed gamma-ray luminosities. Besides, 3D MHD simulations with test particles have demonstrated that a first-order Fermi process within reconnection sites with embedded turbulence results very efficient particle acceleration rates. Employing this acceleration mechanism and the model above, and considering the relevant leptonic and hadronic loss processes in the core region, we computed the spectral energy distribution (SED) of radio galaxies for which very high energy (VHE) emission has been detected (namely, M87, Cen A, Per A, and IC 310). We found that these match very well specially with the VHE observations, therefore strengthening the conclusions above in favour of a core emission origin for the VHE emission of these sources. The model also naturally explains the observed very fast variability of the VHE emission.
Mg is often traced in late-type stars using lines of neutral magnesium, which is expected to be subject to departures from LTE. The astrophysical importance of Mg as well as its relative simplicity from an atomic physics point of view, makes it a prime target and test bed for detailed ab initio non-LTE modelling in stellar atmospheres. For the low-lying states of Mg i, electron collision data were calculated using the R-matrix method. Calculations for collisional broadening by neutral hydrogen were also performed where data were missing. These calculations, together with data from the literature, were used to build a model atom. First, the modelling was tested by comparisons with observed spectra of benchmark stars with well-known parameters. Second, the spectral line behaviour and uncertainties were explored by extensive experiments in which sets of collisional data were changed or removed. The modelled spectra agree well with observed spectra. The line-to-line scatter in the derived abundances shows improvements compared to LTE. The observed Mg emission features at 7 and 12 microns in the spectra of the Sun and Arcturus are reasonably well reproduced. The modelling predicts non-LTE abundance corrections in dwarfs, both solar metallicity and metal-poor, to be very small (of order 0.01 dex), even smaller than found in previous studies. In giants, corrections vary greatly between lines, but can be as large as 0.4 dex. Our results emphasise the need for accurate data of Mg collisions with both electrons and H atoms for precise non-LTE predictions of stellar spectra, but demonstrate that such data can be calculated and that ab initio non-LTE modelling without resort to free parameters is possible. Grids of departure coefficients and abundance corrections for a range of stellar parameters are planned for a forthcoming paper.
HARPS spectra with S/N > 600 for 21 solar twin stars are used to determine very precise (sigma ~ 0.01 dex) differential abundances of C, O, Na, Mg, Al, Si, S, Ca, Ti, Cr, Fe, Ni, Zn, and Y in order to see how well [X/Fe] is correlated with elemental condensation temperature, Tc. In addition, precise (sigma < 0.8 Gyr) stellar ages are obtained by interpolating between Yonsei-Yale isochrones in the logg - Teff diagram. It is confirmed that the ratio between refractory and volatile elements is lower in the Sun than in most of the solar twins, but for many stars, the relation between [X/Fe] and Tc is not well defined. For several elements there is, instead, an astonishingly tight correlation between [X/Fe] and stellar age with amplitudes up to 0.2 dex over an age interval of 8 Gyr in contrast to the lack of correlation between [Fe/H] and age. While [Mg/Fe] increases with age, the s-process element yttrium shows the opposite behavior so that [Y/Mg] can be used as a sensitive chronometer for Galactic evolution. [Na/Fe] and [Ni/Fe] are not well correlated with stellar age, but define a tight Ni-Na relation similar to that previously found for more metal-poor stars. These results provide new constraints on supernovae yields and Galactic evolution. Furthermore, it is found that the C/O ratio evolves very little with time, which is of interest for discussions of the composition of exoplanets.
Neutron star (binary neutron star and neutron star - black hole) mergers are believed to produce short-duration gamma-ray bursts. They are also believed to be the dominant source of gravitational waves to be detected by the advanced LIGO and the dominant source of the heavy r-process elements in the universe. Whether or not these mergers produce short-duration GRBs depends sensitively on the fate of the core of the remnant (whether, and how quickly, it forms a black hole). In this paper, we combine the results of merger calculations and equation of state studies to determine the fate of the cores of neutron star mergers. Using population studies, we can determine the distribution of these fates to compare to observations. We find that black hole cores form quickly only for equations of state that predict maximum non-rotating neutron star masses below 2.3-2.4 solar masses. If quick black hole formation is essential in producing gamma-ray bursts, LIGO observed rates compared to GRB rates could be used to constrain the equation of state for dense nuclear matter.
Thermoelastic distortion resulting from optical absorption by transmissive and reflective optics can cause unacceptable changes in optical systems that employ high power beams. In advanced-generation laser-interferometric gravitational wave detectors for example, optical absorption is expected to result in wavefront distortions that would compromise the sensitivity of the detector; thus necessitating the use of adaptive thermal compensation. Unfortunately, these systems have long thermal time constants and so predictive feed-forward control systems could be required - but the finite-element analysis is computationally expensive. We describe here the use of the Betti-Maxwell elastodynamic reciprocity theorem to calculate the response of linear elastic bodies (optics) to heating that has arbitrary spatial distribution. We demonstrate using a simple example, that it can yield accurate results in computational times that are significantly less than those required for finite-element analyses.
We report new experimental results obtained on three different laser facilities that show directed laser-driven relativistic electron-positron jets with up to 30 times larger yields than previously obtained and a quadratic (~ E^2) dependence of the positron yield on the laser energy. This favorable scaling stems from a combination of higher energy electrons due to increased laser intensity and the recirculation of MeV electrons in the mm-thick target. Based on this scaling, first principles simulations predict the possibility of using such electron-positron jets, produced at upcoming high-energy laser facilities, to probe the physics of relativistic collisionless shocks in the laboratory.
A sinusoidally time-varying pattern for the values of the Newton's constant of gravitation $G$ measured in Earth-based laboratories over the latest decades has been recently reported in the literature. Its amplitude and period amount to $A_G=1.619\times 10^{-14} \textrm{kg}^{-1} \textrm{m}^3 \textrm{s}^{-2}, P_G=5.899 \textrm{yr}$, respectively. Given the fundamental role played by $G$ in the currently accepted theory of gravitation and the attempts to merge it with quantum mechanics, it is important to put to the test the hypothesis that the aforementioned harmonic variation may pertain $G$ itself in a direct and independent way. The bounds on $\dot G/G$ existing in the literature may not be extended straightforwardly to the present case since they were inferred by considering just secular variations. Thus, we numerically integrated the ad-hoc modified equations of motion of the major bodies of the Solar System by finding that the orbits of the planets would be altered by an unacceptably larger amount in view of the present-day high accuracy astrometric measurements. In the case of Saturn, its geocentric right ascension $\alpha$, declination $\delta$ and range $\rho$ would be affected up to $10^4-10^5$ milliarcseconds and $10^5$ km, respectively; the present-day residuals of such observables are as little as about $4$ milliarcseconds and $10^{-1}$ km, respectively.
We propose and study a new class of of superconducting detectors which are sensitive to O(meV) electron recoils from dark matter-electron scattering. Such devices could detect dark matter as light as the warm dark matter limit, mX > keV. We compute the rate of dark matter scattering off free electrons in a (superconducting) metal, including the relevant Pauli blocking factors. We demonstrate that classes of dark matter consistent with all astrophysical and terrestrial constraints could be detected by such detectors with a moderate size exposure.
Antideuterons are a potential messenger for dark matter annihilation or decay in our own galaxy, with very low backgrounds expected from astrophysical processes. The standard coalescence model of antideuteron formation, while simple to implement, is shown to be under considerable strain by recent data from the LHC. We suggest a new empirically based model, with only one free parameter, which is better able to cope with these data, and we explore the consequences of the model for dark matter searches.
In this work, the Friedman equations for hadronic matter in the Robertson-Walker metric in the early Universe are obtained. We consider the hadronic phase, formed after the hadronization of the quark-gluon plasma, that means times from 10^{-6}s to 1s. The set of equations is derived and the behavior of the system is studied considering one approximate analytical solution.
We study the nature of constraints and count the number of degrees of freedom in the non-projectable version of the U(1) extension of Ho\v{r}ava-Lifshitz gravity, using the standard method of Hamiltonian analysis in the classical field theory. This makes it possible for us to investigate the condition under which the scalar graviton is absent in a fully nonlinear level. We show that the scalar graviton does not exist at the classical level if and only if two specific coupling constants are exactly zero. The operators corresponding to these two coupling constants are marginal for any values of the dynamical critical exponent of the Lifshitz scaling and thus should be generated by quantum corrections even if they are eliminated from the bare action. We thus conclude that the theory in general contains the scalar graviton.
An interaction between dark matter and dark energy is usually introduced by a phenomenological modification of the matter conservation equations, while the Einstein equations are left unchanged. Starting from some general and fundamental considerations, in this work it is shown that a coupling in the dark sector is likely to introduce new terms also in the gravitational dynamics. Specifically in the cosmological background equations a bulk dissipative pressure, characterizing viscous effects and able to suppress structure formation at small scales, should appear from the dark coupling. At the level of the perturbations the analysis presented in this work reveals instead the difficulties in properly defining the dark sector interaction from a phenomenological perspective.
We investigate the cosmological behavior in a universe governed by time asymmetric extensions of general relativity, which is a novel modified gravity based on the addition of new, time-asymmetric, terms on the Hamiltonian framework, in a way that the algebra of constraints and local physics remain unchanged. Nevertheless, at cosmological scales these new terms can have significant effects that can alter the universe evolution, both at early and late times, and the freedom in the choice of the involved modification function makes the scenario able to produce a huge class of cosmological behaviors. For basic ansatzes of modification, we perform a detailed dynamical analysis, extracting the stable late time solutions. Amongst others, we find that the universe can result in dark-energy dominated, accelerating solutions, even in the absence of an explicit cosmological constant, in which the dark energy can be quintessence-like, phantom-like, or behave as an effective cosmological constant. Moreover, it can result to matter-domination, or to a Big Rip, or experience the sequence from matter to dark energy domination. Finally, these scenarios can easily satisfy the observational and phenomenological requirements. Hence, time asymmetric cosmology can be a good candidate for the description of the universe.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
We present a strong lensing modeling technique based on versatile basis sets for the lens and source planes. Our method uses high performance Monte Carlo algorithms, allows for an adaptive build up of complexity and bridges the gap between parametric and pixel based reconstruction methods. We apply our method to a HST image of the strong lens system RXJ1131-1231 and show that our method finds a reliable solution and is able to detect substructure in the lens and source planes simultaneously. Using mock data we show that our method is sensitive to sub-clumps with masses four orders of magnitude smaller than the main lens, which corresponds to about $10^8 M_{\odot}$, without prior knowledge on the position and mass of the sub-clump. The modelling approach is flexible and maximises automation to facilitate the analysis of the large number of strong lensing systems expected in upcoming wide field surveys. The resulting search for dark sub-clumps in these systems, without mass-to-light priors, offers promise for probing physics beyond the standard model in the dark matter sector.
Interactions with neutrons and protons play a crucial role for the neutrino opacity of matter in the supernova core. Their current implementation in many simulation codes, however, is rather schematic and ignores not only modifications for the correlated nuclear medium of the nascent neutron star, but also free-space corrections from nucleon recoil, weak magnetism or strange quarks, which can easily add up to changes of several 10% for neutrino energies in the spectral peak. In the Garching supernova simulations with the Prometheus-Vertex code, such sophistications have been included for a long time except for the strange-quark contributions to the nucleon spin, which affect neutral-current neutrino scattering. We demonstrate on the basis of a 20 Msun progenitor star that a moderate strangeness-dependent contribution of g_a^s = -0.2 to the axial-vector coupling constant g_a = 1.26 can turn an unsuccessful three-dimensional (3D) model into a successful explosion. Such a modification is well compatible with current experimental limits and reduces the neutral-current scattering opacity of neutrons, which dominate in the medium around and above the neutrinosphere. This leads to increased luminosities and mean energies of all neutrino species and strengthens the neutrino-energy deposition in the heating layer. Higher nonradial kinetic energy in the gain layer signals enhanced buoyancy activity that enables the onset of the explosion at ~300 ms after bounce, in contrast to the model with vanishing strangeness contributions to neutrino-nucleon scattering. Our results demonstrate the close proximity to explosion of the previously published, unsuccessful 3D models of the Garching group.
The fact that the clustering and concentration of dark matter halos depend not only on their mass, but also the formation epoch, is a prominent, albeit subtle, feature of the cold dark matter structure formation theory, and is known as assembly bias. At low mass scales ($\sim 10^{12}\,h^{-1}M_\odot$), early-forming halos are predicted to be more strongly clustered than the late-forming ones. In this study we aim to robustly detect the signature of assembly bias observationally, making use of formation time indicators of central galaxies in low mass halos as a proxy for the halo formation history. Weak gravitational lensing is employed to ensure our early- and late-forming halo samples have similar masses, and are free of contamination of satellites from more massive halos. For the two formation time indicators used (resolved star formation history and current specific star formation rate), we do not find convincing evidence of assembly bias. For a pair of early- and late-forming galaxy samples with mean mass $M_{200c} \approx 9\times 10^{11}\,h^{-1}M_\odot$, the relative bias is $1.00\pm 0.12$. We attribute the lack of detection to the possibilities that either the current measurements of these indicators are too noisy, or they do not correlate well with the halo formation history. Alternative proxies for the halo formation history that should perform better are suggested for future studies.
Chandra X-ray observations of the nearby brightest cluster galaxy M87 resolve the hot gas structure across the Bondi accretion radius of the central supermassive black hole, a measurement possible in only a handful of systems but complicated by the bright nucleus and jet emission. By stacking only short frame-time observations to limit pileup, and after subtracting the nuclear PSF, we analysed the X-ray gas properties within the Bondi radius at 0.12-0.22 kpc (1.5-2.8 arcsec), depending on the black hole mass. Within 2 kpc radius, we detect two significant temperature components, which are consistent with constant values of 2 keV and 0.9 keV down to 0.15 kpc radius. No evidence was found for the expected temperature increase within ~0.25 kpc due to the influence of the SMBH. Within the Bondi radius, the density profile is consistent with $\rho\propto r^{-1}$. The lack of a temperature increase inside the Bondi radius suggests that the hot gas structure is not dictated by the SMBH's potential and, together with the shallow density profile, shows that the classical Bondi rate may not reflect the accretion rate onto the SMBH. If this density profile extends in towards the SMBH, the mass accretion rate onto the SMBH could be at least two orders of magnitude less than the Bondi rate, which agrees with Faraday rotation measurements for M87. We discuss the evidence for outflow from the hot gas and the cold gas disk and for cold feedback, where gas cooling rapidly from the hot atmosphere could feed the cirumnuclear disk and fuel the SMBH. At 0.2 kpc radius, the cooler X-ray temperature component represents ~20% of the total X-ray gas mass and, by losing angular momentum to the hot gas component, could provide a fuel source of cold clouds within the Bondi radius.
We present a simulation of the formation of the earliest Population II stars, starting from cosmological initial conditions and ending when metals created in the first supernovae are incorporated into a collapsing gas-cloud. This occurs after a supernova blast-wave collides with a nearby mini-halo, inducing further turbulence that efficiently mixes metals into the dense gas in the center of the halo. The gas that first collapses has been enriched to a metallicity of Z ~ 2e-5 Zsun. Due to the extremely low metallicity, collapse proceeds similarly to metal-free gas until dust cooling becomes efficient at high densities, causing the cloud to fragment into a large number of low mass objects. This external enrichment mechanism provides a plausible origin for the most metal-poor stars observed, such as SMSS J031300.36-670839.3, that appear to have formed out of gas enriched by a single supernova. This mechanism operates on shorter timescales than the time for low-mass mini-halos (M < 5e5 Msun) to recover their gas after experiencing a supernova. As such, metal-enriched stars will likely form first via this channel if the conditions are right for it to occur. We identify a number of other externally enriched halos that may form stars in this manner. These halos have metallicities as high as 0.01 Zsun, suggesting that some members of the first generation of metal-enriched stars may be hiding in plain sight in current stellar surveys.
We combine previously published interferometric and single-dish data of relatively nearby massive dense cores that are actively forming stars to test whether their `fragmentation level' is controlled by turbulent or thermal support. We find no clear correlation between the fragmentation level and velocity dispersion, nor between the observed number of fragments and the number of fragments expected when the gravitationally unstable mass is calculated including various prescriptions for `turbulent support'. On the other hand, the best correlations are found for the case of pure thermal Jeans fragmentation, for which we infer a core formation efficiency in the range 12-34%, consistent with previous works. We conclude that the dominant factor determining the fragmentation level of star-forming massive dense cores seems to be thermal Jeans fragmentation.
As is usual in dwarf spheroidal galaxies, today the Local Group galaxy Ursa Minor is depleted of its gas content. How this galaxy lost its gas is still a matter of debate. To study the history of gas loss in Ursa Minor, we conducted the first three-dimensional hydrodynamical simulations of this object, assuming that the gas loss was driven by galactic winds powered only by type II supernovae (SNe II). The initial gas setup and supernova (SN) rates used in our simulations are mainly constrained by the inferred star formation history and the observed velocity dispersion of Ursa Minor. After 3 Gyr of evolution, we found that the gas removal efficiency is higher when the SN rate is increased, and also when the initial mean gas density is lowered. The derived mass-loss rates are systematically higher in the central regions (<300 pc), even though such a relationship has not been strictly linear in time and in terms of the galactic radius. The filamentary structures induced by Rayleigh-Taylor instabilities and the concentric shells related to the acoustic waves driven by SNe can account for the inferred mass losses from the simulations. Our results suggest that SNe II are able to transfer most of the gas from the central region outward to the galactic halo. However, other physical mechanisms must be considered in order to completely remove the gas at larger radii.
The elemental compositions of hot Jupiters are informative relics of planet formation that can help us answer long-standing questions regarding the origin and formation of giant planets. Here, I present the main conclusions from a comprehensive atmospheric retrieval survey of eight hot Jupiters with detectable molecular absorption in their near-infrared transmission spectra. I analyze the eight transmission spectra using the newly-developed, self-consistent atmospheric retrieval framework, SCARLET. Unlike previous methods, SCARLET combines the physical and chemical consistency of complex atmospheric models with the statistical treatment of observational uncertainties known from atmospheric retrieval techniques. I find that all eight hot Jupiters consistently require carbon-to-oxygen ratios (C/O) below 0.9. The finding of C/O<0.9 is highly robust for HD209458b, WASP-12b, WASP-19b, HAT-P-1b, and XO-1b. For HD189733b, WASP-17b, and WASP-43b, I find that the published WFC3 transmission spectra favor C/O<0.9 at greater than 95% confidence. I further show that the water abundances on all eight hot Jupiters are consistent with solar composition. The relatively small depth of the detected water absorption features is due to the presence of clouds, not due to a low water abundance as previously suggested for HD209458b. The presence of a thick cloud deck is inferred for HD209458b and WASP-12b. HD189733b may host a similar cloud deck, rather than the previously suggested Rayleigh hazes, if star spots affect the observed spectrum. The approach taken in SCARLET can be regarded as a new pathway to interpreting spectral observations of planetary atmospheres. In this work, including our prior knowledge of H-C-N-O chemistry enables me to constrain the C/O ratio without detecting a single carbon-bearing molecule.
Collisional excitation rate coefficients play an important role in the dynamics of energy transfer in the interstellar medium. In particular, accurate rotational excitation rates are needed to interpret microwave and infrared observations of the interstellar gas for nonlocal thermodynamic equilibrium line formation. Theoretical cross sections and rate coefficients for collisional deexcitation of rotationally excited HF in the vibrational ground state are reported. The quantum-mechanical close-coupling approach implemented in the nonreactive scattering code MOLSCAT was applied in the cross section and rate coefficient calculations on an accurate 2D HF-He potential energy surface. Estimates of rate coefficients for H and H$_2$ colliders were obtained from the HF-He collisional data with a reduced-potential scaling approach. The calculation of state-to-state rotational quenching cross sections for HF due to He with initial rotational levels up to $j=20$ were performed for kinetic energies from 10$^{-5}$ to 15000 cm$^{-1}$. State-to-state rate coefficients for temperatures between 0.1 and 3000 K are also presented. The comparison of the present results with previous work for lowly-excited rotational levels reveals significant differences. In estimating HF-H$_2$ rate coefficients, the reduced-potential method is found to be more reliable than the standard reduced-mass approach. The current state-to-state rate coefficient calculations are the most comprehensive to date for HF-He collisions. We attribute the differences between previously reported and our results to differences in the adopted interaction potential energy surfaces. The new He rate coefficients can be used in a variety of applications. The estimated H$_2$ and H collision rates can also augment the smaller datasets previously developed for H$_2$ and electrons.
We report results of a search for CIII] $\lambda \lambda$1907,1909 {\AA} emission using Keck's MOSFIRE spectrograph in a sample of 7 $z_{phot}\sim7-8$ candidates ($H\sim27$) lensed by the Hubble Frontier Field cluster Abell 2744. Earlier work has suggested the promise of using the CIII] doublet for redshift confirmation of galaxies in the reionization era given $Ly\alpha$ ($\lambda$1216 {\AA}) is likely attenuated by the neutral intergalactic medium. The primary challenge of this approach is the feasibility of locating CIII] emission without advanced knowledge of the spectroscopic redshift. With an integration time of 5 hours in the H-band, we reach a $5\sigma$ median flux limit (in between the skylines) of $1.5\times10^{-18}$ ergs cm$^{-2}$ sec$^{-1}$ but no convincing CIII] emission was found. We also incorporate preliminary measurements from two other CLASH/HFF clusters in which, similarly, no line was detected, but these were observed to lesser depth. Using the known distribution of OH emission and the photometric redshift likelihood distribution of each lensed candidate, we present statistical upper limits on the mean total CIII] rest-frame equivalent width for our $z\simeq7-8$ sample. For a signal/noise ratio of 5, we estimate the typical CIII] doublet rest-frame equivalent width is, with 95\% confidence, $<26\pm5$ {\AA}. Although consistent with the strength of earlier detections in brighter objects at $z\simeq6-7$, our study illustrates the necessity of studying more luminous or strongly-lensed examples prior to the launch of the James Webb Space Telescope.
The recent cosmological observations are in good agreement with the scalar spectral index $n_s$ with $n_s-1\sim -2/N$, where $N$ is the number of e-foldings. Quadratic chaotic model, Starobinsky model and Higgs inflation or $\alpha$-attractors connecting them are typical examples predicting such a relation. We consider the problem in the opposite: given $n_s$ as a function of $N$, what is the inflaton potential $V(\phi)$. We find that for $n_s-1=-2/N$, $V(\phi)$ is either $\tanh^2(\gamma\phi/2)$ ("T-model") or $\phi^2$ (chaotic inflation) to the leading order in the slow-roll approximation. $\gamma$ is the ratio of $1/V$ at $N\rightarrow \infty$ to the slope of $1/V$ at a finite $N$ and is related to "$\alpha$" in the $\alpha$-attractors by $\gamma^2=2/3\alpha$. The tensor-to-scalar ratio $r$ is $r=8/N(\gamma^2 N +1) $. The implications for the reheating temperature are also discussed. We also derive formulas for $n_s-1=-p/N$. Although $r$ depends on a parameter, the running of the spectral index is independent of it, which can be used as a consistency check of the assumed relation of $n_s(N)$.
We present the identification of potential members of nearby Galactic globular clusters using radial velocities from the RAdial Velocity Experiment Data Release 4 (RAVE-DR4) survey database. Our identifications are based on three globular clusters -- NGC 3201, NGC 5139 ($\omega$ Cen) and NGC 362 -- all of which are shown to have |RV|>100 km/s. The identification of globular cluster stars in RAVE DR4 data offers a unique opportunity to test the precision and accuracy of the stellar parameters determined with the currently available Stellar Parameter Pipelines (SPPs) used in the survey, as globular clusters are ideal testbeds for the validation of stellar atmospheric parameters, abundances, distances and ages. For both NGC 3201 and $\omega$ Cen, there is compelling evidence for numerous members (> 10) in the RAVE database; in the case of NGC 362 the evidence is more ambiguous, and there may be significant foreground and/or background contamination in our kinematically-selected sample. A comparison of the RAVE-derived stellar parameters and abundances with published values for each cluster and with BASTI isochrones for ages and metallicities from the literature reveals overall good agreement, with the exception of the apparent underestimation of surface gravities for giants, in particular for the most metal-poor stars. Moreover, if the selected members are part of the main body of each cluster our results would also suggest that the distances from Binney et al. 2013, where only isochrones more metal-rich than -0.9 dex were used, are typically underestimated by ~ 40% with respect to the published distances for the clusters, while the distances from Zwitter et al. 2010 show stars ranging from 1 to ~ 6.5 kpc -- with indications of a trend toward higher distances at lower metallicities -- for the three clusters analysed in this study.
The population of Low Surface Brightness (LSB) galaxies is crucial for understanding the extremes of galaxy formation and evolution of the universe. As LSB galaxies are mostly rich in gas (HI), the alpha.40-SDSS DR7 sample is absolutely one of the best survey combinations to select a sample of them in the local Universe. Since the sky backgrounds are systematically overestimated for galaxy images by the SDSS photometric pipeline, particularly for luminous galaxies or galaxies with extended low surface brightness outskirts, in this paper, we above all estimated the sky backgrounds of SDSS images in the alpha.40-SDSS DR7 sample, using a precise method of sky subtraction. Once subtracting the sky background, we did surface photometry with the Kron elliptical aperture and fitted geometric parameters with an exponential profile model for each galaxy image. Basing on the photometric and geometric results, we further calculated the B-band central surface brightness, mu_{0}(B), for each galaxy and ultimately defined a sample of LSB galaxies consisting of 1129 galaxies with mu_{0}(B) > 22.5 mag arcsec^{-2} and the minor-to-major axis ratio b/a > 0.3. This HI-selected LSB galaxy sample from the alpha.40-SDSS DR7 is a relatively unbiased sample of gas-rich and disk-dominated LSB galaxies, which is complete both in HI observation and the optical magnitude within the limit of SDSS DR7 photometric survey. We made analysis on optical and radio 21cm HI properties of this LSB galaxy sample. Additionally, we statistically investigated the environment of our LSB galaxies, and found that up to 92% of the total LSB galaxies have less than 8 neighbouring galaxies, which strongly evidenced that LSB galaxies prefer to reside in the low-density environment.
The Five-hundred-meter Aperture Spherical Radio Telescope (FAST) is supported by a cable-net structure, whose change in form leads to a stress range of approximately 500MPa. This stress range is more than twice the standard authorized value. The cable-net structure is thus the most critical and fragile part of the FAST reflector system. In this study, we first search for a more appropriate deformation strategy that reduces the stress amplitude generated by the form-changing operation. Second, we roughly estimate the tracking trajectory of the telescope during its service life, and conduct an extensive numerical investigation to assess the fatigue resistance requirements. Finally, we develop a new type of steel cable system that meets that cable requirements for FAST construction.
This paper introduces the data cubes from GHIGLS, our deep Green Bank Telescope surveys of the 21-cm line emission of HI in targeted fields at intermediate Galactic latitude. The GHIGLS fields together cover over 800 square degrees at 9.55' spatial resolution. The HI spectra have an effective velocity resolution about 1.0 km/s and cover at least -450 < v < +250 km/s. GHIGLS highlights that even at intermediate Galactic latitude the interstellar medium is very complex. Spatial structure of the HI is quantified through power spectra of maps of the column density, NHI. For our featured representative field, centered on the North Ecliptic Pole, the scaling exponents in power-law representations of the power spectra of NHI maps for low, intermediate, and high velocity gas components (LVC, IVC, and HVC) are -2.90 +/- 0.03, -2.55 +/- 0.04, and -2.66 +/- 0.06, respectively. After Gaussian decomposition of the line profiles, NHI maps were also made corresponding to the broad and narrow line components in the LVC range; for these maps the exponents are -2.8 +/- 0.1 and -2.2 +/- 0.2, respectively, the latter revealing more small scale structure in the cold neutral medium (CNM). There is evidence that filamentary structure in the HI CNM is oriented parallel to the Galactic magnetic field. The power spectrum analysis also offers insight into the various contributions to uncertainty in the data. The effect of 21-cm line opacity on the GHIGLS NHI maps is estimated.
The computed $\log(\tau/{\rm y})$-$\log(m/m_\odot)$ relation for the stellar initial mass range, 0.6-120.0, and the stellar initial metallicity range, 0.0004-0.0500, tabulated in an earlier attempt (Portinari et al. 1998) is fitted to a good extent by a four-parameter curve, expressed by a double power-law, for assigned stellar initial metallicity, which can be reduced to a three-parameter curve, expressed by a single power-law, for the whole set of stellar initial metallicities. The relative errors do not exceed about 2% and 4%, respectively. The extent to which the interpolation curve, expressed by a single power-law, can be extrapolated towards both high-mass and low-mass stars, is also investigated. High-mass star lifetimes are understimated by a factor less than 2 up to $m/m_\odot= 1000$ and by a fiducial factor less than 4 up to infinite. Low-mass star lifetimes are overstimated by a factor of about 3 down to $m/m_\odot=0.25$ and by an unacceptably large factor down to $m/m_\odot=0.08$. The interpolation curve, expressed by a single power-law, is related (in differential form) to a generalization of the equation of the exponential decay, which could be the starting point for a theoretical interpretation. As a simple application, the star mass fraction of a single star generation with stellar initial mass function defined by a power-law, is plotted vs. the logarithmic stellar lifetime. The star mass fraction declines in time at a decreasing rate for mild stellar initial mass function and at an increasing rate for steep stellar initial mass function, where a linear trend is exhibited for a value of the exponent close to the Salpeter's value, equal to -2.35.
The {\it Kepler} space mission and its {\it K2} extension provide photometric
time series data with unprecedented accuracy. These data challenge our current
understanding of the metallic-lined A stars (Am stars) for what concerns the
onset of pulsations in their atmospheres. It turns out that the predictions of
current diffusion models do not agree with observations. To understand this
discrepancy, it is of crucial importance to obtain ground-based spectroscopic
observations of Am stars in the {\it Kepler} and {\it K2} fields in order to
determine the best estimates of the stellar parameters.
In this paper, we present a detailed analysis of high-resolution
spectroscopic data for seven stars previously classified as Am stars. We
determine the effective temperatures, surface gravities, projected rotational
velocities, microturbulent velocities and chemical abundances of these stars
using spectral synthesis. These spectra were obtained with {\it CAOS}, a new
instrument recently installed at the observing station of the Catania
Astrophysical Observatory on Mt. Etna. Three stars have already been observed
during quarters Q0-Q17, namely: HD\,180347, HD\,181206, and HD\,185658, while
HD\,43509 was already observed during {\it K2} C0 campaign.
We confirm that HD\,43509 and HD\,180347 are Am stars, while HD 52403,
HD\,50766, HD\,58246, HD\,181206 and HD\,185658 are marginal Am stars. By means
of non-LTE analysis, we derived oxygen abundances from O{\sc
I}$\lambda$7771--5{\AA} triplet and we also discussed the results obtained with
both non-LTE and LTE approaches.
We present observations of M31LRN 2015 (MASTER OT J004207.99+405501.1), discovered in M31 in January 2015, and identified as a rare and enigmatic luminous red nova (LRN). Spectroscopic and photometric observations obtained by the Liverpool Telescope showed the LRN becoming extremely red as it faded from its M(V) = -9.4 +/- 0.2 peak. Early spectra showed strong Halpha emission that weakened over time as a number of absorption features appeared, including Na I D and Ba II. At later times strong TiO absorption bands were also seen. A search of archival Hubble Space Telescope data revealed a luminous red source to be the likely progenitor system, with pre-outburst Halpha emission also detected in ground-based data. The outburst of M31LRN 2015 shows many similarities, both spectroscopically and photometrically, with that of V838 Mon, the best studied LRN. We finally discuss the possible progenitor scenarios.
AA Dor is a close, totally eclipsing, post common-envelope binary with an
sdOB-type primary and an extremely low-mass secondary, located close to the
mass limit of stable central hydrogen burning. Within error limits, it may
either be a brown dwarf or a late M-type dwarf.
We aim to extract the secondary's contribution to the phase-dependent
composite spectra. The spectrum and identified lines of the secondary decide on
its nature.
In January 2014, we measured the phase-dependent spectrum of AA Dor with
XSHOOTER over one complete orbital period. Since the secondary's rotation is
presumable synchronized with the orbital period, its surface strictly divides
into a day and night side. Therefore, we may obtain the spectrum of its cool
side during its transit and of its hot, irradiated side close to its
occultation. We developed the Virtual Observatory (VO) tool TLISA to search for
weak lines of a faint companion in a binary system.
We identified 53 spectral lines of the secondary in the ultraviolet-blue,
visual, and near-infrared XSHOOTER spectra that are strongest close to its
occultation. We identified 57 (20 additional) lines in available UVES
(Ultraviolet and Visual Echelle Spectrograph) spectra from 2001. The lines are
mostly from C II-III and O II, typical for a low-mass star that is irradiated
and heated by the primary. We verified the orbital period of P = 22597.033201
+/- 0.00007 s and determined the orbital velocity Ksec = 232.9 (+16.6 / -6.5)
km/s of the secondary. The mass of the secondary is Msec = 0.081 (+0.018 /
-0.010) Msun and, hence, it is not possible to reliably determine a brown dwarf
or an M-type dwarf nature.
Although we identified many emission lines of the secondary's irradiated
surface, the resolution and signal-to-noise ratio of our UVES and XSHOOTER
spectra are not good enough to extract a good spectrum of the secondary's
nonirradiated hemisphere.
We argue that the intracluster medium (ICM) of the cooling flow galaxy group NGC 5813 is more likely to be heated by mixing of post-shock jets' gas residing in hot bubbles with the ICM, than by shocks or turbulent-heating. Efficient shocks-heating of the inner cooling flow region would have caused overheating of the rest of the ICM, leading to its ejection from the group. Heating by mixing, that was found to be much more efficient than turbulent-heating and shocks-heating, hence, rescues the outer ICM of NGC 5813 from its predestined fate according to cooling flow feedback scenarios that are based on heating by shocks.
Multi-object astronomical adaptive-optics (MOAO) is now a mature wide-field
observation mode to enlarge the adaptive-optics-corrected field in a few
specific locations over tens of arc-minutes.
The work-scope provided by open-loop tomography and pupil conjugation is
amenable to a spatio-angular Linear-Quadratic Gaussian (SA-LQG) formulation
aiming to provide enhanced correction across the field with improved
performance over static reconstruction methods and less stringent computational
complexity scaling laws.
Starting from our previous work [1], we use stochastic time-progression
models coupled to approximate sparse measurement operators to outline a
suitable SA-LQG formulation capable of delivering near optimal correction.
Under the spatio-angular framework the wave-fronts are never explicitly
estimated in the volume,providing considerable computational savings on
10m-class telescopes and beyond.
We find that for Raven, a 10m-class MOAO system with two science channels,
the SA-LQG improves the limiting magnitude by two stellar magnitudes when both
Strehl-ratio and Ensquared-energy are used as figures of merit. The
sky-coverage is therefore improved by a factor of 5.
The coalescence of compact binaries containing neutron stars or black holes is one of the most promising signals for advanced ground-based laser interferometer gravitational-wave detectors, with the first direct detections expected over the next few years. The rate of binary coalescences and the distribution of component masses is highly uncertain, and population synthesis models predict a wide range of plausible values. Poorly constrained parameters in population synthesis models correspond to poorly understood astrophysics at various stages in the evolution of massive binary stars, the progenitors of binary neutron star and binary black hole systems. These include effects such as supernova kick velocities, parameters governing the energetics of common envelope evolution and the strength of stellar winds. Observing multiple binary black hole systems through gravitational waves will allow us to infer details of the astrophysical mechanisms that lead to their formation. We simulate gravitational-wave observations from a series of population synthesis models including the effects of known selection biases, measurement errors and cosmology. We compare the predictions arising from different models and show that, with observations (or the lack of them) from the early runs of the advanced LIGO and Virgo detectors, we will begin to be able to pin down details of the astrophysical processes that lead to black hole binary formation.
Generally, the oblateness of a planet or moon is what causes rings to settle into its equatorial plane. However, the recent suggestion that a ring system might exist (or have existed) about Rhea, a moon whose shape includes a strong prolate component pointed toward Saturn, raises the question of whether rings around a triaxial primary can be stable. We study the role of prolateness in the behavior of rings around Rhea and extend our results to similar problems such as possible rings around exoplanets. Using a Hamiltonian approach, we point out that the dynamical behavior of ring particles is governed by three different time scales: the orbital period of the particles, the rotation period of the primary, and the precession period of the particles' orbital plane. In the case of Rhea, two of these are well separated from the third, allowing us to average the Hamiltonian twice. To study the case of slow rotation of the primary, we also carry out numerical simulations of a thin disk of particles undergoing secular effects and damping. For Rhea, the averaging reduces the Hamiltonian to an oblate potential, under which rings would be stable only in the equatorial plane. This is not the case for Iapetus; rather, it is the lack of a prolate component to its shape that allows Iapetus to host rings. Plausible exoplanets should mostly be in the same regime as Rhea, though other outcomes are possible. The numerical simulations indicate that, even when the double averaging is irrelevant, rings settle in the equatorial plane on an approximately constant time scale.
We analyse intensity variations, as measured by the Atmospheric Imaging Assembly (AIA) in the 171 {\AA} passband, in two coronal loops embedded within a single coronal magnetic arcade. We detect oscillations in the fundamental mode with periods of roughly 2 minutes and decay times of 5 minutes. The oscillations were initiated by interaction of the arcade with a large wavefront issuing from a flare site. Further, the power spectra of the oscillations evince signatures consistent with oblique propagation to the field lines and for the existence of a 2-D waveguide instead of a 1-D one.
Observations of several nearby galaxies of regular magnetic fields reveal magnetic arms situated between the material arms. The nature of these magnetic arms is a topic of active debate. Previously we found a hint that taking into account the effects of injections of small-scale magnetic fields generated, e.g., by turbulent dynamo action, into the large-scale galactic dynamo can result in magnetic arm formation. We now investigate the joint roles of an arm/interarm turbulent diffusivity contrast and injections of small-scale magnetic field on the formation of large-scale magnetic field ("magnetic arms") in the interarm region. We use the relatively simple "no-$z$" model for the galactic dynamo. This involves projection on to the galactic equatorial plane of the azimuthal and radial magnetic field components; the field component orthogonal to the galactic plane is estimated from the solenoidality condition. We find that addition of diffusivity gradients to the effect of magnetic field injections makes the magnetic arms much more pronounced. In particular, the regular magnetic field component becomes larger in the interarm space compared to that within the material arms.The joint action of the turbulent diffusivity contrast and small-scale magnetic field injections (with the possible participation of other effects previously suggested) appears to be a plausible explanation for the phenomenon of magnetic arms.
We present a set of convective dynamo simulations in rotating spherical fluid shells based on an anelastic approximation of compressible fluids. The simulations extend into a "buoyancy-dominated" regime where the buoyancy forcing is dominant while the Coriolis force is no longer balanced by pressure gradients and strong anti-solar differential rotation develops as a result. Dynamos in this regime are strongly dominated by dipole components but at the same time their magnetic energies are relatively small compared to the corresponding kinetic energies of the flow. Despite being relatively weak the self-sustained magnetic fields are able to reverse the direction of differential rotation to solar-like. We find that the convection in the buoyancy-dominated regime is significantly stronger near the pole than in the equatorial region, leading to non-oscillatory dipolar dynamo solutions. The results are obtained with a new simulation code for modelling of convection and MHD dynamo generation in rotating spherical shells under the anelastic approximation. The model equations and the new code are also presented here along with code validation results based on benchmark dynamo models.
In this work, we establish an ICME list from 1996 to 2014 based on the in-situ observations from the WIND and ACE satellites. Based on this ICME list, we extend the statistical analysis of the ICMEs to the solar maximum phase of solar cycle 24th. The analysis of the annual variations of the properties of ICMEs show that the number of ICMEs, the number of shocks, the percentage of ICMEs drove shocks, the magnetic field and plasma properties of ICMEs are well correlated with the solar cycle variation. The number of MCs do not show any correlation with sunspot number. But, the percentage of the MCs in ICMEs show good anti-correlation with the sunspot number. By comparison the parameters of MCs with None-MC ICMEs, we found that the MCs are stronger than the None-MC ICMEs . In addition, we compare the parameters of ICMEs with and without shocks. It is found that the ICMEs with shocks are much stronger than the ICME without shocks. Meanwhile, we discuss the distribution of the magnetic field and solar wind plasmas parameters of the sheath regions of ICMEs at first time. We find that the magnetic field and solar wind velocity in the sheath region are higher than them in the ejecta of ICMEs from statistical point of view.
The origin of large magnetic fields in the Universe remains currently unknown. We investigate here a mechanism before recombination based on known physics. The source of the vorticity is due to the changes in the photon distribution function caused by the fluctuations in the background photons. We show that the magnetic field generated in the MHD limit, due to the Coulomb scattering, is of the order $10^{-49}$ G. We explicitly show that the magnetic fields generated from this process are sustainable and are not erased by resistive diffusion. We compare the results with current observations and discuss the implications.
Based on the Sloan Digital Sky Survey DR 7, we investigate the environment, morphology and stellar population of bulgeless low surface brightness (LSB) galaxies in a volume-limited sample with redshift ranging from 0.024 to 0.04 and $M_r$ $\leq$ $-18.8$. The local density parameter $\Sigma_5$ is used to trace their environments. We find that, for bulgeless galaxies, the surface brightness does not depend on the environment. The stellar populations are compared for bulgeless LSB galaxies in different environments and for bulgeless LSB galaxies with different morphologies. The stellar populations of LSB galaxies in low density regions are similar to those of LSB galaxies in high density regions. Irregular LSB galaxies have more young stars and are more metal-poor than regular LSB galaxies. These results suggest that the evolution of LSB galaxies may be driven by their dynamics including mergers rather than by their large scale environment.
In this paper we combine WFC3/UVIS F275W, F336W, and F438W data from the "UV Legacy Survey of Galactic Globular Clusters: Shedding Light on Their Populations and Formation" (GO-13297) HST Treasury program with F606W, F625W, F658N, and F814W ACS archive data for a multi-wavelength study of the globular cluster NGC 6352. In the color-magnitude and two-color diagrams obtained with appropriate combination of the photometry in the different bands we separate two distinct stellar populations and trace them from the main sequence to the subgiant, red giant, horizontal and asymptotic giant branches. We infer that the two populations differ in He by Delta Y=0.029+/-0.006. With a new method, we also estimate the age difference between the two sequences. Assuming no difference in [Fe/H] and [alpha/Fe], and the uncertainties on Delta Y, we found a difference in age between the two populations of 10+/-120 Myr. If we assume [Fe/H] and [alpha/Fe] differences of 0.02 dex (well within the uncertainties of spectroscopic measurements), the total uncertainty in the relative age rises to ~300 Myr.
We present high-speed, three-colour photometry of the eclipsing dwarf nova PHL 1445, which, with an orbital period of 76.3 min, lies just below the period minimum of ~82 min for cataclysmic variable stars. Averaging four eclipses reveals resolved eclipses of the white dwarf and bright spot. We determined the system parameters by fitting a parameterised eclipse model to the averaged lightcurve. We obtain a mass ratio of q = 0.087 +- 0.006 and inclination i = 85.2 +- 0.9 degrees. The primary and donor masses were found to be Mw = 0.73 +- 0.03 Msun and Md = 0.064 +- 0.005 Msun, respectively. Through multicolour photometry a temperature of the white dwarf of Tw = 13200 +- 700 K and a distance of 220 +- 50 pc were determined. The evolutionary state of PHL 1445 is uncertain. We are able to rule out a significantly evolved donor, but not one that is slightly evolved. Formation with a brown dwarf donor is plausible; though the brown dwarf would need to be no older than 600 Myrs at the start of mass transfer, requiring an extremely low mass ratio (q = 0.025) progenitor system. PHL 1445 joins SDSS 1433 as a sub-period minimum CV with a substellar donor. These existence of two such systems raises an alternative possibility; that current estimates for the intrinsic scatter and/or position of the period minimum may be in error.
Photometry alone is not sufficient to unambiguously distinguish between ultra-faint star clusters and dwarf galaxies because of their overlap in morphological properties. Here we report on VLT/GIRAFFE spectra of candidate member stars in two recently discovered ultra-faint satellites Reticulum 2 and Horologium 1, obtained as part of the ongoing Gaia-ESO Survey. We identify 18 members in Reticulum 2 and 5 in Horologium 1. We find Reticulum 2 to have a velocity dispersion of ~3.22 km/s, implying a M/L ratio of ~ 500. We have inferred stellar parameters for all candidates and we find Reticulum 2 to have a mean metallicity of [Fe/H] = -2.46+/-0.1, with an intrinsic dispersion of ~ 0.29, and is alpha-enhanced to the level of [alpha/Fe]~0.4. We conclude that Reticulum 2 is a dwarf galaxy. We also report on the serendipitous discovery of four stars in a previously unknown stellar substructure near Reticulum 2 with [Fe/H] ~ -2 and V_hel ~ 220 km/s, far from the systemic velocity of Reticulum 2. For Horologium 1 we infer a velocity dispersion of 4.9^{+2.8}_{-0.9} km/s and a consequent M/L ratio of ~ 600, leading us to conclude that Horologium 1 is also a dwarf galaxy. Horologium 1 is slightly more metal-poor than Reticulum 2 [Fe/H] = -2.76 +/- 0.1 and is similarly alpha-enhanced: [alpha/Fe] ~ 0.3. Despite a large error-bar, we also measure a significant spread of metallicities of 0.17 dex which strengthen the evidence that Horologium 1 is indeed a dwarf galaxy. The line-of-sight velocity of Reticulum 2 is offset by some 100 km/s from the prediction of the orbital velocity of the LMC, thus making its association with the Cloud uncertain. However, at the location of Horologium 1, both the backward integrated LMC's orbit and the LMC's halo are predicted to have radial velocities similar to that of the dwarf. Therefore, it is likely that Horologium 1 is or once was a member of the Magellanic Family.
We present X-ray timing and spectral analyses of simultaneous 150 ks Nuclear Spectroscopic Telescope Array (NuSTAR) and Suzaku X-ray observations of the Seyfert 1.5 galaxy NGC 4151. We disentangle the continuum emission, absorption, and reflection properties of the active galactic nucleus (AGN) by applying inner accretion disk reflection and absorption-dominated models. With a time-averaged spectral analysis, we find strong evidence for relativistic reflection from the inner accretion disk. We find that relativistic emission arises from a highly ionized inner accretion disk with a steep emissivity profile, which suggests an intense, compact illuminating source. We find a preliminary, near-maximal black hole spin a>0.9 accounting for statistical and systematic modeling errors. We find a relatively moderate reflection fraction with respect to predictions for the lamp post geometry, in which the illuminating corona is modeled as a point source. Through a time-resolved spectral analysis, we find that modest coronal and inner disk reflection flux variation drives the spectral variability during the observations. We discuss various physical scenarios for the inner disk reflection model, and we find that a compact corona is consistent with the observed features.
In this work, we show that a large class of models with a composite dark sector undergo a strong first order phase transition in the early universe, which could lead to a detectable gravitational wave signal. We summarise the basic conditions for a strong first order phase transition for SU(N) dark sectors with n_f flavours, calculate the gravitational wave spectrum and show that, depending on the dark confinement scale, it can be detected at eLISA or in pulsar timing array experiments. The gravitational wave signal provides a unique test of the gravitational interactions of a dark sector, and we discuss the complementarity with conventional searches for new dark sectors. The discussion includes Twin Higgs and SIMP models as well as symmetric and asymmetric composite dark matter scenarios.
In this note, we amend the incorrect discussion in Nucl. Phys. B 886 (2014) 569 [1] concerning the numerical examples considered there. In particular, we discuss the viability of minimal radiative models of Resonant Leptogenesis and prove that no asymmetry can be generated at O(h^4) in these scenarios. We present a minimal modification of the model considered in [1], where electroweak-scale right-handed Majorana neutrinos can easily accommodate both successful leptogenesis and observable signatures at Lepton Number and Flavour Violation experiments. The importance of the fully flavour-covariant rate equations, as developed in [1], for describing accurately the generation of the asymmetry is reconfirmed.
Most extreme-mass-ratio-inspirals of small compact objects into supermassive black holes end with a fast plunge from an eccentric last stable orbit. For rapidly rotating black holes such fast plunges may be studied in the context of the Kerr/CFT correspondence because they occur in the near-horizon region where dynamics are governed by the infinite dimensional conformal symmetry. In this paper we use conformal transformations to analytically solve for the radiation emitted from fast plunges into near-extreme Kerr black holes. We find perfect agreement between the gravity and CFT computations.
We consider inflation in a universe with a positive cosmological constant and a nonminimally coupled scalar field, in which the field couples both quadratically and quartically to the Ricci scalar. When considered in the Einstein frame and when the nonminimal couplings are negative, the field starts in slow roll and inflation ends with an asymptotic value of the principal slow roll parameter, $\epsilon_E=4/3$. Graceful exit can be achieved by suitably (tightly) coupling the scalar field to matter, such that at late time the total energy density reaches the scaling of matter, $\epsilon_E=\epsilon_m$. Quite generically the model produces a red spectrum of scalar cosmological perturbations and a small amount of gravitational radiation. With a suitable choice of the nonminimal couplings, the spectral slope can be as large as $n_s\simeq 0.955$, which is about one standard deviation away from the central value measured by the Planck satellite. The model can be ruled out by future measurements if any of the following is observed: (a) the spectral index of scalar perturbations is $n_s>0.960$; (b) the amplitude of tensor perturbations is above about $r\sim 10^{-2}$; (c) the running of the spectral index of scalar perturbations is positive.
We propose a dark matter explanation to simultaneously account for the excess of antiproton-to-proton and positron power spectra observed in the AMS-02 experiment while having the right dark matter relic abundance and satisfying the current direct search bounds. We extend the Higgs triplet model with a hidden gauge symmetry of $SU(2)_X$ that is broken to $Z_3$ by a quadruplet scalar field, rendering the associated gauge bosons stable weakly-interacting massive particle dark matter candidates. By coupling the complex Higgs triplet and the $SU(2)_X$ quadruplet, the dark matter candidates can annihilate into triplet Higgs bosons each of which in turn decays into lepton or gauge boson final states. Such a mechanism gives rise to correct excess of positrons and antiprotons with an appropriate choice of the triplet vacuum expectation value. Besides, the model provides a link between neutrino mass and dark matter phenomenology.
Classical chaos is often characterized as exponential divergence of nearby trajectories. In many interesting cases these trajectories can be identified with geodesic curves. We define here the entropy by $S = \ln \chi (x)$ with $\chi(x)$ being the distance between two nearby geodesics. We derive an equation for the entropy which by transformation to a Ricatti-type equation becomes similar to the Jacobi equation. We further show that the geodesic equation for a null geodesic in a double warped space time leads to the same entropy equation. By applying a Robertson-Walker metric for a flat three-dimensional Euclidian space expanding as a function of time, we again reach the entropy equation stressing the connection between the chosen entropy measure and time. We finally turn to the Raychaudhuri equation for expansion, which also is a Ricatti equation similar to the transformed entropy equation. Those Ricatti-type equations have solutions of the same form as the Jacobi equation. The Raychaudhuri equation can be transformed to a harmonic oscillator equation, and it has been shown that the geodesic deviation equation of Jacobi is essentially equivalent to that of a harmonic oscillator. The Raychaudhuri equations are strong geometrical tools in the study of General Relativity and Cosmology. We suggest a refined entropy measure applicable in Cosmology and defined by the average deviation of the geodesics in a congruence.
Astroinformatics is a new impact area in the world of astronomy, occasionally called the final frontier, where several astrophysicists, statisticians and computer scientists work together to tackle various data intensive astronomical problems. Exponential growth in the data volume and increased complexity of the data augments difficult questions to the existing challenges. Classical problems in Astronomy are compounded by accumulation of astronomical volume of complex data, rendering the task of classification and interpretation incredibly laborious. The presence of noise in the data makes analysis and interpretation even more arduous. Machine learning algorithms and data analytic techniques provide the right platform for the challenges posed by these problems. A diverse range of open problem like star-galaxy separation, detection and classification of exoplanets, classification of supernovae is discussed. The focus of the paper is the applicability and efficacy of various machine learning algorithms like K Nearest Neighbor (KNN), random forest (RF), decision tree (DT), Support Vector Machine (SVM), Na\"ive Bayes and Linear Discriminant Analysis (LDA) in analysis and inference of the decision theoretic problems in Astronomy. The machine learning algorithms, integrated into ASTROMLSKIT, a toolkit developed in the course of the work, have been used to analyze HabCat data and supernovae data. Accuracy has been found to be appreciably good.
Experiments searching for weak interacting massive particles with noble gases such as liquid argon require very low detection thresholds for nuclear recoils. A determination of the scintillation efficiency is crucial to quantify the response of the detector at low energy. We report the results obtained with a small liquid argon cell using a monoenergetic neutron beam produced by a deuterium-deuterium fusion source. The light yield relative to electrons was measured for six argon recoil energies between 11 and 120 keV at zero electric drift field.
We propose a new scenario of baryogenesis, in which annihilation of axion domain walls generates a sizable baryon asymmetry. Successful baryogenesis is possible for a wide range of the axion mass and decay constant, $m \simeq 10^8 -10^{13}$ GeV and $f \simeq 10^{13} - 10^{16}$ GeV. Baryonic isocurvature perturbations are significantly suppressed in our model, in contrast to various spontaneous baryogenesis scenarios in the slow-roll regime. In particular, the axion domain wall baryogenesis is consistent with high-scale inflation which generates a large tensor-to-scalar ratio within the reach of future CMB B-mode experiments. We also discuss the gravitational waves produced by the domain wall annihilation and its implications for the future gravitational wave experiments.
We consider the possibility to produce a bouncing universe in the framework of scalar-tensor gravity models in which the scalar field potential may be negative, and even unbounded from below. We find a set of viable solutions with nonzero measure in the space of initial conditions passing a bounce, even in the presence of a radiation component, and approaching a constant gravitational coupling afterwards. Hence we have a model with a minimal modification of gravity in order to produce a bounce in the early universe with gravity tending dynamically to GR after the bounce.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)
We present radio and X-ray observations of the nearby Type IIb Supernova 2013df in NGC4414 from 10 to 250 days after the explosion. The radio emission showed a peculiar soft-to-hard spectral evolution. We present a model in which inverse Compton cooling of synchrotron emitting electrons can account for the observed spectral and light curve evolution. A significant mass loss rate, $\dot{M} \approx 8 \times 10^{-5}\,\rm M_{\odot}/yr$ for a wind velocity of 10 km/s, is estimated from the detailed modeling of radio and X-ray emission, which are primarily due to synchrotron and bremsstrahlung, respectively. We show that SN 2013df is similar to SN 1993J in various ways. The shock wave speed of SN 2013df was found to be average among the radio supernovae; $v_{sh}/c \sim 0.07$. We did not find any significant deviation from smooth decline in the light curve of SN 2013df. One of the main results of our self-consistent multiband modeling is the significant deviation from energy equipartition between magnetic fields and relativistic electrons behind the shock. We estimate $\epsilon_{e} = 200 \epsilon_{B}$. In general for Type IIb SNe, we find that the presence of bright optical cooling envelope emission is linked with free-free radio absorption and bright thermal X-ray emission. This finding suggests that more extended progenitors, similar to that of SN 2013df, suffer from substantial mass loss in the years before the supernova.
We have considered the star formation properties of 1616 isolated galaxies
from the 2MASS XSC selected sample (2MIG) with the FUV GALEX magnitudes. This
sample was then compared with corresponding properties of isolated galaxies
from the Local Orphan Galaxies catalogue (LOG) and paired galaxies.
We found that different selection algorithms define different populations of
isolated galaxies. The population of the LOG catalogue, selected from
non-clustered galaxies in the Local Supercluster volume, mostly consists of
low-mass spiral and late type galaxies. The SSFR upper limit in isolated and
paired galaxies does not exceed the value of ~dex(-9.4). This is probably
common for galaxies of differing activity and environment (at least at z<0.06).
The fractions of quenched galaxies are nearly twice as high in the paired
galaxy sample as in the 2MIG isolated galaxy sample. From the behaviour of
(S)SFR vs. M_* relations we deduced that the characteristic value influencing
evolutionary processes is the galaxy mass. However the environmental influence
is notable: paired massive galaxies with logM_*>11.5 have higher (S)SFR than
isolated galaxies. Our results suggest that the environment helps to trigger
the star formation in the highest mass galaxies. We found that the fraction of
AGN in the paired sample is only a little higher than in our isolated galaxy
sample. We assume that AGN phenomenon is probably defined by secular galaxy
evolution.
We present a study of 33 planet-candidate host stars for which asteroseismic observations have sufficiently high signal-to-noise ratio to allow extraction of individual pulsation frequencies. We implement a new Bayesian scheme that is flexible in its input to process individual oscillation frequencies, combinations of them, and average asteroseismic parameters, and derive robust fundamental properties for these targets. Applying this scheme to grids of evolutionary models yields stellar properties with median statistical uncertainties of 1.1% (radius), 1.7% (density), 3.3% (mass), 4.4% (distance), and 14% (age), making this the exoplanet host-star sample with the most precise and uniformly determined fundamental parameters to date. We assess the systematics from changes in the solar abundances and mixing-length parameter, showing that they are smaller than the statistical errors. We also determine the stellar properties with three other fitting algorithms and explore the systematics arising from using different evolution and pulsation codes, resulting in 1% in density and radius, and 2% and 7% in mass and age, respectively. We identify the initial helium abundance as a source of systematics comparable to our statistical uncertainties, and discuss future prospects for constraining this parameter by combining asteroseismology and data from space missions. Finally we compare our derived properties with those obtained using the global average asteroseismic observables along with effective temperature and metallicity, finding an excellent level of agreement. Owing to selection effects, our results show that the majority of the high signal-to-noise ratio asteroseismic Kepler host stars are older than the Sun.
We report constraints on the three-dimensional orbital architecture for all four planets known to orbit the nearby M dwarf Gliese 876 (=GJ 876) based solely on Doppler measurements and demanding long-term orbital stability. Our dataset incorporates publicly available radial velocities taken with the ELODIE and CORALIE spectrographs, HARPS, and Keck HIRES as well as previously unpublished HIRES RVs. We first quantitatively assess the validity of the planets thought to orbit GJ 876 by computing the Bayes factors for a variety of different coplanar models using an importance sampling algorithm. We confirm that a four-planet model is indeed preferred over a three-planet model. Next, we apply a Newtonian MCMC algorithm (RUNDMC, Nelson et al. 2014) to perform a Bayesian analysis of the planet masses and orbits using an n-body model that allows each planet to take on its own orbit in three-dimensional space. Based on the radial velocities alone, the mutual inclinations for the outer three resonant planets are constrained to $\Phi_{cb}=2.8\pm^{1.7}_{1.3}$ degrees for the "c" and "b" pair and $\Phi_{be}=10.3\pm^{6.3}_{5.1}$ degrees for the "b" and "e" pair. We integrate the equations of motion of a sample of initial conditions drawn from our posterior for $10^7$ years. We identify dynamically unstable models and find that the GJ 876 planets must be roughly coplanar ($\Phi_{cb}=1.41\pm^{0.62}_{0.57}$ degrees and $\Phi_{be} = 3.9\pm^{2.0}_{1.9}$ degrees), indicating the amount of planet-planet scattering in the system has been low. We investigate the distribution of the respective resonant arguments of each planet pair and find that at least one resonant argument for each planet pair and the Laplace argument librate. The libration amplitudes in our three-dimensional orbital model supports the idea of the outer-three planets having undergone significant past disk migration.
We investigate the effect of dimensionality on the transition to explosion in neutrino-driven core-collapse supernovae. Using parameterized hydrodynamic simulations of the stalled supernova shock in one-, two- (2D), and three spatial dimensions (3D), we systematically probe the extent to which hydrodynamic instabilities alone can tip the balance in favor of explosion. In particular, we focus on systems that are well into the regimes where the Standing Accretion Shock Instability (SASI) or neutrino-driven convection dominate the dynamics, and characterize the difference between them. We find that SASI-dominated models can explode with up to ~20% lower neutrino luminosity in 3D than in 2D, with the magnitude of this difference decreasing with increasing resolution. This improvement in explosion conditions originates in the ability of spiral modes to generate more non-radial kinetic energy than a single sloshing mode, increasing the size of the average shock radius, and hence generating better conditions for the formation of large-scale, high-entropy bubbles. In contrast, convection-dominated explosions show a smaller difference in their critical heating rate between 2D and 3D (<15%), in agreement with previous studies. The ability of our numerical implementation to maintain arbitrary symmetries is quantified with a set of SASI-based tests. We discuss implications for the diversity of explosion paths in a realistic supernova environment.
We present initial results from the low-latitude Galactic plane region of the High Time Resolution Universe pulsar survey conducted at the Parkes 64-m radio telescope. We discuss the computational challenges arising from the processing of the terabyte-sized survey data. Two new radio interference mitigation techniques are introduced, as well as a partially-coherent segmented acceleration search algorithm which aims to increase our chances of discovering highly-relativistic short-orbit binary systems, covering a parameter space including potential pulsar-black hole binaries. We show that under a constant acceleration approximation, a ratio of data length over orbital period of ~0.1 results in the highest effectiveness for this search algorithm. From the 50 per cent of data processed thus far, we have re-detected 435 previously known pulsars and discovered a further 60 pulsars, two of which are fast-spinning pulsars with periods less than 30ms. PSR J1101-6424 is a millisecond pulsar whose heavy white dwarf (WD) companion and short spin period of 5.1ms indicate a rare example of full-recycling via Case A Roche lobe overflow. PSR J1757-27 appears to be an isolated recycled pulsar with a relatively long spin period of 17ms. In addition, PSR J1244-6359 is a mildly-recycled binary system with a heavy WD companion, PSR J1755-25 has a significant orbital eccentricity of 0.09, and PSR J1759-24 is likely to be a long-orbit eclipsing binary with orbital period of the order of tens of years. Comparison of our newly-discovered pulsar sample to the known population suggests that they belong to an older population. Furthermore, we demonstrate that our current pulsar detection yield is as expected from population synthesis.
We present the BaLROG (Bars in Low Redshift Optical Galaxies) sample of 16 morphologically distinct barred spirals to characterise observationally the influence of bars on nearby galaxies. Each galaxy is a mosaic of several pointings observed with the IFU spectrograph SAURON leading to a tenfold sharper spatial resolution (~100 pc) compared to ongoing IFU surveys. In this paper we focus on the kinematic properties. We calculate the bar strength Qb from classical torque analysis using 3.6 {\mu}m Spitzer (S4G) images, but also develop a new method based solely on the kinematics. A correlation between the two measurements is found and backed up by N-body simulations, verifying the measurement of Qb . We find that bar strengths from ionised gas kinematics are ~2.5 larger than those measured from stellar kinematics and that stronger bars have enhanced influence on inner kinematic features. We detect that stellar angular momentum "dips" at 0.2$\pm$0.1 bar lengths and half of our sample exhibits an anti-correlation of h3 - stellar velocity (v/{\sigma}) in these central parts. An increased flattening of the stellar {\sigma} gradient with increasing bar strength supports the notion of bar-induced orbit mixing. These measurements set important constraints on the spatial scales, namely an increasing influence in the central regions (0.1-0.5 bar lengths), revealed by kinematic signatures due to bar-driven secular evolution in present day galaxies.
We report on both high-precision photometry from the MOST space telescope and ground-based spectroscopy of the triple system delta Ori A consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7d, and a more distant tertiary (O9 IV P > 400 yrs). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for 3 weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the ?first time in non-phased data. From the spectroscopy we have a well constrained radial velocity curve of Aa1. While we are unable to recover radial velocity variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries which show evidence for tidally induced pulsations.
We present fluxes in both neutral carbon [CI] lines at the centers of 76 galaxies with FIR luminosities between 10^{9} and 10^{12} L(o) obtained with Herschel-SPIRE and with ground-based facilities, along with the J=7-6, J=4-3, J=2-1 12CO and J=2-1 13CO line fluxes. We investigate whether these lines can be used to characterize the molecular ISM of the parent galaxies in simple ways and how the molecular gas properties define the model results. In most starburst galaxies, the [CI]/13CO flux ratio is much higher than in Galactic star-forming regions, and it is correlated to the total FIR luminosity. The [CI](1-0)/CO(4-3), the [CI](2-1) (2-1)/CO(7-6), and the [CI] (2-1)/(1-0) flux ratios are also correlated, and trace the excitation of the molecular gas. In the most luminous infrared galaxies (LIRGs), the ISM is fully dominated by dense and moderately warm gas clouds that appear to have low [C]/[CO] and [13CO]/[12CO] abundances. In less luminous galaxies, emission from gas clouds at lower densities becomes progressively more important, and a multiple-phase analysis is required to determine consistent physical characteristics. Neither the CO nor the [CI] velocity-integrated line fluxes are good predictors of H2 column densities in individual galaxies, and X(CI) conversion factors are not superior to X(CO) factors. The methods and diagnostic diagrams outlined in this paper also provide a new and relatively straightforward means of deriving the physical characteristics of molecular gas in high-redshift galaxies up to z=5, which are otherwise hard to determine.
The detection of high energy neutrinos ($10^{15}-10^{20}$ eV or $1-10^{5}$ PeV) is an important step toward understanding the most energetic cosmic accelerators and would enable tests of fundamental physics at energy scales that cannot easily be achieved on Earth. In this energy range, there are two expected populations of neutrinos: the astrophysical flux observed with IceCube at lower energies ($\sim1$ PeV) and the predicted cosmogenic flux at higher energies ($\sim10^{18}$ eV). Radio detector arrays such as RICE, ANITA, ARA, and ARIANNA exploit the Askaryan effect and the radio transparency of glacial ice, which together enable enormous volumes of ice to be monitored with sparse instrumentation. We describe here the design for a phased radio array that would lower the energy threshold of radio techniques to the PeV scale, allowing measurement of the astrophysical flux observed with IceCube over an extended energy range. Meaningful energy overlap with optical Cherenkov telescopes could be used for energy calibration. The phased radio array design would also provide more efficient coverage of the large effective volume required to discover cosmogenic neutrinos.
The Kleinmann-Low nebula in Orion (Orion-KL) is the nearest example of a
high-mass star-forming environment. For the first time, we complemented 1.3 mm
Submillimeter Array (SMA) interferometric line survey with IRAM 30 m
single-dish observations of the Orion-KL region. Covering a 4 GHz bandwidth in
total, this survey contains over 160 emission lines from 20 species (25
isotopologues), including 10 complex organic molecules (COMs).
At a spatial resolution of 1200 AU, the continuum substructures are resolved.
Extracting the spectra from individual substructures and providing the
intensity-integrated distribution map for each species, we studied the
small-scale chemical variations in this region. Our main results are: (1) We
identify lines from the low-abundance COMs CH3COCH3 and CH3CH2OH, as well as
tentatively detect CH3CHO and long carbon-chains C6H and HC7N. (2) We find that
while most COMs are segregated by type, peaking either towards the hot core
(e.g., N-bearing species) or the compact ridge (e.g., O-bearing species like
HCOOCH3, CH3OCH3, the distributions of others do not follow this segregated
structure (e.g., CH3CH2OH, CH3OH, CH3COCH3). (3) We find a second velocity
component of HNCO, 34SO2, and SO lines, which may be associated with a strong
shock event in the low-velocity outflow. (4) Temperatures and molecular
abundances show large gradients between central condensations and the outflow
regions, illustrating a transition between hot molecular core and
shock-chemistry dominated regimes.
Our observations of spatially resolved chemical variations in Orion-KL
provide the nearest reference source for hot molecular core and outflow
chemistry, which will be an important example for interpreting the chemistry of
more distant HMSFRs.
The Swift Gamma Ray Burst Explorer has proven to be an incredible platform for studying the multiwavelength properties of supernova explosions. In its first ten years, Swift has observed over three hundred supernovae. The ultraviolet observations reveal a complex diversity of behavior across supernova types and classes. Even amongst the standard candle type Ia supernovae, ultraviolet observations reveal distinct groups. When the UVOT data is combined with higher redshift optical data, the relative populations of these groups appear to change with redshift. Among core-collapse supernovae, Swift discovered the shock breakout of two supernovae and the Swift data show a diversity in the cooling phase of the shock breakout of supernovae discovered from the ground and promptly followed up with Swift. Swift observations have resulted in an incredible dataset of UV and X-ray data for comparison with high-redshift supernova observations and theoretical models. Swift's supernova program has the potential to dramatically improve our understanding of stellar life and death as well as the history of our universe.
We present models that predict spectra of old- and intermediate-aged stellar populations at 2.51\AA\ (FWHM) with varying [\alpha/Fe] abundance. The models are based on the MILES library and on corrections from theoretical stellar spectra. The models employ recent [Mg/Fe] determinations for the MILES stars and BaSTI scaled-solar and \alpha-enhanced isochrones. We compute models for a suite of IMF shapes and slopes, covering a wide age/metallicity range. Using BaSTI, we also compute "base models" matching The Galactic abundance pattern. We confirm that the \alpha-enhanced models show a flux excess with respect to the scaled-solar models blue-ward $\sim$4500\AA, which increases with age and metallicity. We also confirm that both [MgFe] and [MgFe]' indices are [\alpha/Fe]-insensitive. We show that the sensitivity of the higher order Balmer lines to [\alpha/Fe] resides in their pseudo-continua, with narrower index definitions yielding lower sensitivity. We confirm that the \alpha-enhanced models yield bluer (redder) colours in the blue (red) spectral range. To match optical colours of massive galaxies we require both \alpha-enhancement and a bottom-heavy IMF. The comparison of Globular Cluster line-strengths with our predictions match the [Mg/Fe] determinations from their individual stars. We obtain good fits to both full spectra and indices of galaxies with varying [\alpha/Fe]. Using thousands of SDSS galaxy spectra we obtain a linear relation between a proxy for the abundance, [Z$_{\rm Mg}$/Z$_{\rm Fe}$]$_{\rm SS(BaSTI)}$, using solely scaled-solar models and the [Mg/Fe] derived with models with varying abundance ([Mg/Fe]=0.59[Z$_{\rm Mg}$/Z$_{\rm Fe}$]$_{\rm SS(BaSTI)}$). Finally we provide a user-friendly, web-based facility, which allows composite populations with varying IMF and [\alpha/Fe] ( this http URL ).
The Canada-France Ecliptic Plane Survey discovered four trans-Neptunian objects with semi-major axes near the 5:1 resonance, revealing a large and previously undetected intrinsic population. Three of these objects are currently resonant with Neptune, and the fourth is consistent with being an object that escaped the resonance at some point in the past. The non-resonant object may be representative of a detached population that is stable at slightly lower semi-major axes than the 5:1 resonance. We generated clones of these objects by resampling the astrometric uncertainty and examined their behavior over a 4.5 Gyr numerical simulation. The majority of the clones of the three resonant objects (>90%) spend a total of 10^7 years in resonance during their 4.5 Gyr integrations; most clones experience multiple periods of resonance capture. Our dynamical integrations reveal an exchange between the 5:1 resonance, the scattering objects, and other large semi-major axis resonances, especially the 4:1, 6:1, and 7:1. The multiple capture events and relatively short resonance lifetimes after capture suggest that these objects are captured scattering objects that stick in the 5:1 resonance. These 5:1 resonators may be representative of a temporary population, requiring regular contributions from a source population. We examined the dynamical characteristics (inclination, eccentricity, resonant island, libration amplitude) of the detected objects and their clones in order to provide an empirical model of the orbit structure of the 5:1 resonance. This resonance is dynamically hot and includes primarily symmetric librators. Given our orbit model, the intrinsic population necessary for the detection of these three objects in the 5:1 resonance is 1900(+3300 -1400, 95% confidence) with H_g< 8 and e > 0.5.
Radio core dominance, the rest-frame ratio of core to lobe luminosity, has been widely used as a measure of Doppler boosting of a quasar's radio jets and hence of the inclination of the central engine's spin axis to the line of sight. However, the use of the radio lobe luminosity in the denominator (essentially to try and factor out the intrinsic power of the central engine) has been criticized and other proxies for the intrinsic engine power have been proposed. These include the optical continuum luminosity, and the luminosity of the narrow-line region. Each is plausible, but so far none has been shown to be clearly better than the others. In this paper we evaluate four different measures of core dominance using a new sample of 126 radio loud quasars, carefully selected to be as free as possible of orientation bias, together with high quality VLA images and optical spectra from the SDSS. We find that normalizing the radio core luminosity by the optical continuum luminosity yields a demonstrably superior orientation indicator. In addition, by comparing the equivalent widths of broad emission lines in our orientation-unbiased sample to those of sources in the MOJAVE program, we show that the beamed optical synchrotron emission from the jets is not a significant component of the optical continuum for the sources in our sample. We also discuss future applications of these results.
We report results from a continuation of our searches for high field magnetic white dwarfs paired in a detached binary with non degenerate companions. We made use of the Sloan Digital Sky Survey DR7 catalog of Kleinman et al. (2013) with 19,712 spectroscopically-identified white dwarfs. These include 1,735 white dwarf plus M dwarf detached pairs (almost 10\% of the Kleinman at al.'s list). No new pairs were found, although we did recover the polar (AM~Herculis system) ST\,LMi in a low state of accretion. With the larger sample the original situation reported ten years ago remains intact now at a much higher level of statistical significance: in the selected SDSS sample, high field magnetic white dwarfs are not found in an apparently-detached pairing with an M dwarf, unless they are a magnetic CV in a low state of accretion. This finding strengthens the case that the fields in the isolated high field magnetic white dwarfs are generated by stellar mergers but also raises questions on the nature of the progenitors of the magnetic CVs.
By performing a global magneto-hydrodynamical simulation for the Milky Way with an axisymmetric gravitational potential, we propose that spatially dependent amplification of magnetic fields possibly explains the observed noncircular motion of the gas in the Galactic center region. The radial distribution of the rotation frequency in the bulge region is not monotonic in general. The amplification of the magnetic field is enhanced in regions with stronger differential rotation, because magnetorotational instability and field-line stretching are more effective. The strength of the amplified magnetic field reaches >~ 0.5 mG, and radial flows of the gas are excited by the inhomogeneous transport of angular momentum through turbulent magnetic field that is amplified in a spatially dependent manner. In addition, the magnetic pressure-gradient force also drives radial flows in a similar manner. As a result, the simulated position-velocity diagram exhibits a time-dependent asymmetric parallelogram-shape owing to the intermittency of the magnetic turbulence; the present model provides a viable alternative to the bar-potential-driven model for the parallelogram-shape of the central molecular zone. This is a natural extension into the central few 100 pc of the magnetic activity, which is observed as molecular loops at radii from a few 100 pc to 1 kpc. Furthermore, the time-averaged net gas flow is directed outward, whereas the flows are highly time-dependent, which we discuss from a viewpoint of the outflow from the bulge.
Context. Solar-like oscillations exhibit a regular pattern of frequencies.
This pattern is dominated by the small and large frequency separations between
modes. The accurate determination of these parameters is of great interest,
because they give information about e.g. the evolutionary state and the mass of
a star.
Aims. We want to develop a robust method to determine the large and small
frequency separations for time series with low signal-tonoise ratio. For this
purpose, we analyse a time series of the Sun from the GOLF instrument aboard
SOHO and a time series of the star KIC 5184732 from the NASA Kepler satellite
by employing a combination of Fourier and Hilbert transform.
Methods. We use the analytic signal of filtered stellar oscillation time
series to compute the signal envelope. Spectral analysis of the signal envelope
then reveals frequency differences of dominant modes in the periodogram of the
stellar time series.
Results. With the described method the large frequency separation $\Delta\nu$
can be extracted from the envelope spectrum even for data of poor
signal-to-noise ratio. A modification of the method allows for an overview of
the regularities in the periodogram of the time series.
Meridional flow is thought to play a very important role in the dynamics of
the solar convection zone; however, because of its relatively small amplitude,
precisely measuring it poses a significant challenge. Here we present a
complete time-distance helioseismic analysis of about two years of ground-based
GONG Doppler data to retrieve the meridional circulation profile for modest
latitudes, in an attempt to corroborate results from other studies. We use an
empirical correction to the travel times due to an unknown center-to-limb
systematic effect. The helioseismic inversion procedure is first tested and
reasonably validated on artificial data from a large-scale numerical
simulation, followed by a test to broadly recover the solar differential
rotation found from global seismology.
From GONG data, we measure poleward photospheric flows at all latitudes with
properties that are comparable with earlier studies, and a shallow equatorward
flow about $65$\,Mm beneath the surface, in agreement with recent findings from
HMI data. No strong evidence of multiple circulation cells in depth nor
latitude is found, yet the whole phase space has not yet been explored. Tests
of mass flux conservation are then carried out on the inferred GONG and HMI
flows and compared to a fiducial numerical baseline from models, and we find
that the continuity equation is poorly satisfied.
While the two disparate data sets do give similar results for about the outer
$15\%$ of the interior radius, the total inverted circulation pattern appears
to be unphysical in terms of mass conservation when interpreted over modest
time scales. We can likely attribute this to both the influence of realization
noise and subtle effects in the data and measurement procedure.
In this paper we review the current status of research on the observational
and theoretical characteristics of isolated and binary magnetic white dwarfs
(MWDs).
Magnetic fields of isolated MWDs are observed to lie in the range 10^3-10^9G.
While the upper limit cutoff appears to be real, the lower limit is more
difficult to investigate. The incidence of magnetism below a few 10^3G still
needs to be established by sensitive spectropolarimetric surveys conducted on
8m class telescopes.
Highly magnetic WDs tend to exhibit a complex and non-dipolar field structure
with some objects showing the presence of higher order multipoles. There is no
evidence that fields of highly magnetic WDs decay over time, which is
consistent with the estimated Ohmic decay times scales of ~10^11 yrs. MWDs, as
a class, also appear to be more massive than their weakly or non-magnetic
counterparts.
MWDs are also found in binary systems where they accrete matter from a
low-mass donor star. These binaries, called magnetic Cataclysmic Variables
(MCVs) and comprise about 20-25\% of all known CVs. Zeeman and cyclotron
spectroscopy of MCVs have revealed the presence of fields in the range $\sim
7-230$\,MG. Complex field geometries have been inferred in the high field MCVs
(the polars) whilst magnetic field strength and structure in the lower field
group (intermediate polars, IPs) are much harder to establish.
The origin of fields in MWDs is still being debated. While the fossil field
hypothesis remains an attractive possibility, field generation within the
common envelope of a binary system has been gaining momentum, since it would
explain the absence of MWDs paired with non-degenerate companions and also the
lack of relatively wide pre-MCVs.
Enormous progress has been made on observing stellar magnetism in stars from
the main sequence through to compact objects. Recent data have thrown into
sharper relief the vexed question of the origin of stellar magnetic fields,
which remains one of the main unanswered questions in astrophysics. In this
chapter we review recent work in this area of research. In particular, we look
at the fossil field hypothesis which links magnetism in compact stars to
magnetism in main sequence and pre-main sequence stars and we consider why its
feasibility has now been questioned particularly in the context of highly
magnetic white dwarfs. We also review the fossil versus dynamo debate in the
context of neutron stars and the roles played by key physical processes such as
buoyancy, helicity, and superfluid turbulence,in the generation and stability
of neutron star fields.
Independent information on the internal magnetic field of neutron stars will
come from future gravitational wave detections. Thus we maybe at the dawn of a
new era of exciting discoveries in compact star magnetism driven by the opening
of a new, non-electromagnetic observational window.
We also review recent advances in the theory and computation of
magnetohydrodynamic turbulence as it applies to stellar magnetism and dynamo
theory. These advances offer insight into the action of stellar dynamos as well
as processes whichcontrol the diffusive magnetic flux transport in stars.
We present the study of the effects of high energy cosmic rays (CRs) over the astrophysical ices, observed toward the embedded class I protostar Elias 29, by using computational modeling and laboratory data. Its spectrum was observed with {\it Infrared Space Observatory - ISO}, covering 2.3 - 190 $\mu$m. The modeling employed the three-dimensional Monte Carlo radiative transfer code RADMC-3D (Dullemond et al. 2012) and laboratory data of bombarded ice grains by CRs analogs, and unprocessed ices (not bombarded). We are assuming that Elias 29 has a self-irradiated disk with inclination $i =$ 60$^{\circ}$, surrounded by an envelope with bipolar cavity. The results show that absorption features toward Elias 29, are better reproduced by assuming a combination between unprocessed astrophysical ices at low temperature (H$_2$O, CO, CO$_2$) and bombarded ices (H$_2$O:CO$_2$) by high energy CRs. Evidences of the ice processing around Elias 29 can be observed by the good fitting around 5.5-8.0 $\mu$m, by polar and apolar ice segregation in 15.15-15.25 $\mu$m, and by presence of the CH$_4$ and HCOOH ices. Given that non-nitrogen compounds were employed in this work, we assume that absorption around 5.5-8.0 $\mu$m should not be associated with NH$_4^+$ ion (Shutte & Khanna2003), but more probably with aliphatic ethers (e.g. R1-OCH$_2$-R2), CH$_3$CHO and related species. The results obtained in this paper are important, because they show that the environment around protostars is better modeled considering processed samples and, consequently, demonstrates the chemical evolution of the astrophysical ices.
We propose an in-situ formation model for inverse-polarity solar prominence and demonstrate it using self-consistent 2.5-dimensional magnetohydrodynamics simulations, including thermal conduction along magnetic fields and optically thin radiative cooling. The model enables us to form cool dense plasma clouds inside a flux rope by radiative condensation, which is regarded as an inverse-polarity prominence. Radiative condensation is triggered by changes in the magnetic topology, i.e., formation of the flux rope from the sheared arcade field, and by thermal imbalance due to the dense plasma trapped inside the flux rope. The flux rope is created by imposing converging and shearing motion on the arcade field. Either when the footpoint motion is in the anti-shearing direction or when heating is proportional to local density, the thermal state inside the flux rope becomes cooling-dominant, leading to radiative condensation. By controlling the temperature of condensation, we investigate the relationship between the temperature and density of prominence and derive a scaling formula for this relationship. This formula suggests that the proposed model reproduce the observed density of prominence, which is 10-100 times larger than the coronal density. Moreover, the time evolution of the extreme ultraviolet emission synthesized by combining our simulation results with the response function of the Solar Dynamics Observatory Atmospheric Imaging Assembly filters agrees with the observed temporal and spatial intensity shift among multi-wavelength during in-situ condensation.
Type-I X-ray bursts on the surface of a neutron star are a unique probe to the accretion in X-ray binary systems. However, we know little about the feedback of the burst emission upon accretion. Hard X-ray shortages and enhancements of the persistent emission at soft X-rays have been observed. To put these findings in context with the aim of understanding the possible mechanism underneath, we investigated 68 bursts seen by RXTE from the clocked burster GS 1826--238. We diagnosed jointly the burst influence at both soft and hard X-rays, and found that the observations can be described as the CompTT model with variable normalization, electron temperature and optical depth. Putting these results in a scenario of coronal Compton cooling via the burst emission would lead to a shortage of the cooling power, which may suggest that additional consideration like the influence of the burst on the corona formation should be accounted for as well.
In this paper we discuss the influence of gravitational instabilities in massive protostellar discs on the dynamics of dust grains. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation, we have computed the evolution of the dust in a quasi-static gas density structure typical of self-gravitating disc. For different grain size distributions we have investigated the capability of spiral arms to trap particles. We have run 3D radiative transfer simulations in order to construct maps of the expected emission at (sub-)millimetre and near-infrared wavelengths. Finally, we have simulated realistic observations of our disc models at (sub-)millimetre and near-infrared wavelengths as they may appear with the Atacama Large Millimetre/sub-millimetre Array (ALMA) and the High-Contrast Coronographic Imager for Adaptive Optics (HiCIAO) in order to investigate whether there are observational signatures of the spiral structure. We find that the pressure inhomogeites induced by gravitational instabilities produce a non-negligible dynamical effect on centimetre sized particles leading to significant overdensities in spiral arms. We also find that the spiral structure is readily detectable by ALMA over a wide range of (sub-)millimetre wavelengths and by HiCIAO in near-infrared scattered light for non-face-on discs located in the Ophiucus star-forming region. In addition, we find clear spatial spectral index variations across the disc, revealing that the dust trapping produces a migration of large grains that can be potentially investigated through multi-wavelenghts observations in the (sub-)millimetric. Therefore, the spiral arms observed to date in protoplanetary disc might be interpreted as density waves induced by the development of gravitational instabilities.
Chinese Spectral Radio Heliograph (CSRH) is an advanced aperture synthesis solar radio heliograph, developed by National Astronomical Observatories, Chinese Academy of Sciences independently. It consists of 100 reflector antennas, which are grouped into two antenna arrays (CSRH-I and CSRH-II) for low and high frequency bands respectively. The frequency band of CSRH-I is 0.4-2GHz and for CSRH-II, the frequency band is 2-15GHz. In the antenna and feed system, CSRH uses an Eleven feed to receive signals coming from the Sun, the radiation pattern with lower side lobe and back lobe of the feed is well radiated. The characteristics of gain G and antenna noise temperature T effect the quality of solar radio imaging. For CSRH, measured G is larger than 60 dBi and $ T $ is less than 120K, after CSRH-I was established, we have successfully captured a solar radio burst between 1.2-1.6GHz on November 12, 2010 through this instrument and this event was confirmed through the observation of Solar Broadband Radio Spectrometer (SBRS) at 2.84GHz and Geostationary Operational Environmental Satellite (GOES). In addition, an image obtained from CSRH-I clearly reveals the profile of the solar radio burst. The other observational work is the imaging of Fengyun-2E geosynchronous satellite which is assumed to be a point source. This data processing method indicates that, the method of deleting errors about dirty image could be used for processing other surface sources.
Quasi-constant heating at the footpoints of loops leads to evaporation and condensation cycles of the plasma: thermal non-equilibrium (TNE). This phenomenon is believed to play a role in the formation of prominences and coronal rain. However, it is often discarded to be involved in the heating of warm loops as the models do not reproduce observations. Recent simulations have shown that these inconsistencies with observations may be due to oversimplifications of the geometries of the models. In addition, our recent observations reveal that long-period intensity pulsations (several hours) are common in solar coronal loops. These periods are consistent with those expected from TNE. The aim of this paper is to derive characteristic physical properties of the plasma for some of these events to test the potential role of TNE in loop heating. We analyzed three events in detail using the six EUV coronal channels of SDO/AIA. We performed both a Differential Emission Measure (DEM) and a time-lag analysis, including a new method to isolate the relevant signal from the foreground and background emission. For the three events, the DEM undergoes long-period pulsations, which is a signature of periodic heating even though the loops are captured in their cooling phase, as is the bulk of the active regions. We link long-period intensity pulsations to new signatures of loop heating with strong evidence for evaporation and condensation cycles. We thus witness simultaneously widespread cooling and TNE. Finally, we discuss the implications of our new observations for both static and impulsive heating models.
We consider generic models of quintessence and we investigate the influence of massive neutrino matter with field-dependent masses on the matter power spectrum. In case of minimally coupled neutrino matter, we examine the effect in tracker models with inverse power-law and double exponential potentials. We present detailed investigations for the scaling field with a steep exponential potential, non-minimally coupled to massive neutrino matter, and we derive constraints on field-dependent neutrino masses from the observational data.
Plasma velocity and magnetic field measurements from the Voyager 2 mission are used to study solar wind turbulence in the slow solar wind at two different heliocentric distances, 5 and 29 astronomical units, sufficiently far apart to provide information on the radial evolution of this turbulence. The magnetic helicity and the cross-helicity, which express the correlation between the plasma velocity and the magnetic field, are used to characterize the turbulence. Wave number spectra are computed by means of the Taylor hypothesis applied to time resolved single point Voyager 2 measurements. The overall picture we get is complex and difficult to interpret. A substantial decrease of the cross-helicity at smaller scales (over 1-3 hours of observation) with increasing heliocentric distance is observed. At 5 AU the only peak in the probability density of the normalized residual energy is negative, near -0.5. At 29 AU the probability density becomes doubly peaked, with a negative peak at -0.5 and a smaller peak at a positive values of about 0.7. A decrease of the cross-helicity for increasing heliocentric distance is observed, together with a reduction of the unbalance toward the magnetic energy of the energy of the fluctuations. For the smaller scales, we found that at 29 AU the normalized polarization is small and positive on average (about 0.1), it is instead zero at 5 AU. For the larger scales, the polarization is low and positive at 5 AU (average around 0.1) while it is negative (around - 0.15) at 29 AU.
We examine the fate of a dead radio source in which jet injection from central engine has stopped at early stage of its evolution ($t = t_j \lesssim 10^5$ yr). To this aim, we theoretically evaluate the evolution of the emission from both lobe and shell which are composed of shocked jet matter and shocked ambient medium, respectively. Based on a simple dynamical model of expanding lobe and shell, we clarify how the broadband spectrum of each component evolves before and after the cessation of the jet activity. It is shown that the spectrum is strongly dominated by the lobe emission while the jet is active ($t \leq t_j$). On the other hand, once the jet activity has ceased ($t > t_j$), the lobe emission fades out rapidly, since fresh electrons are no longer supplied from the jet. Meanwhile, shell emission only shows gradual decrease, since accelerated electrons are continuously supplied from the bow shock that is propagating into the ambient medium. As a result, overall emission from the shell overwhelms that from the lobe at wide range of frequencies from radio up to gamma-ray soon after the jet activity has ceased. Our result predicts a new class of dead radio sources that are dominated by shell emission. We suggest that the emission from the shell can be probed in particular at a radio wavelengths with the Square Kilometer Array (SKA) phase 1.
The early history of the Vogt-Russell theorem is retraced following its route starting at the realization of a correlation between mass and luminosity of binary and pulsating stars, through the embossing of this observation into a theorem, and finally to the emerging first signs of its failure to serve as a theorem in the strict mathematical sense of the word.
In this article, I have studied the cosmoparticle constraints on a generic class of large field ($|\Delta\phi|>M_{p}$) and small field ($|\Delta\phi|<M_{p}$) models of brane inflationary magnetogenesis aka "braneflamagnetogenesis" from: (1) tensor-to-scalar ratio ($r$), (2) reheating, (3) leptogenesis and (4) baryogenesis in case of Randall-Sundrum single braneworld gravity (RSII) framework. I also establish a direct connection between the magnetic field at the present epoch ($B_{0}$) and primordial gravity waves ($r$), which give a precise estimate of non-vanishing CP asymmetry ($\epsilon_{CP}$) in leptogenesis and baryon asymmetry ($\eta_{B}$) in baryogenesis scenario respectively. Further assuming the conformal invariance to be restored after inflation in the framework of RSII, I have explicitly shown that the requirement of the sub-dominant feature of large scale coherent magnetic field after inflation gives two fold non-trivial characteristic constraints- on equation of state parameter ($w$) and the corresponding energy scale during reheating ($\rho^{1/4}_{rh}$) epoch. Hence giving the proposal for avoiding the contribution of back-reaction from the magnetic field I have established a bound on the generic reheating characteristic parameter ($R_{rh}$) and its rescaled version ($R_{sc}$), to achieve large scale magnetic field within the prescribed setup and further apply the CMB constraints as obtained from recently observed Planck 2015 data and Planck+BICEP2+Keck Array joint constraints. Using all these derived results I have shown that it is possible to put further stringent constraints on various classes of large and small field inflationary models to break the degeneracy between various cosmophenomenological parameters within the framework of RSII. Finally, I have studied the consequences from two specific models of brane inflation- monomial and hilltop.
Deep R-band CCD linear polarimetry collected for fields with lines-of-sight toward the Lupus I molecular cloud is used to investigate the properties of the magnetic field within this molecular cloud. The observed sample contains about 7000 stars, almost 2000 of them with polarization signal-to-noise ratio larger than 5. These data cover almost the entire main molecular cloud and also sample two diffuse infrared patches in the neighborhood of Lupus I. The large scale pattern of the plane-of-sky projection of the magnetic field is perpendicular to the main axis of Lupus I, but parallel to the two diffuse infrared patches. A detailed analysis of our polarization data combined with the Herschel/SPIRE 350 um dust emission map shows that the principal filament of Lupus I is constituted by three main clumps acted by magnetic fields having different large-scale structure properties. These differences may be the reason for the observed distribution of pre- and protostellar objects along the molecular cloud and its apparent evolutive stage. On the other hand, assuming that the magnetic field is composed by a large-scale and a turbulent components, we find that the latter is rather similar in all three clumps. The estimated plane-of-sky component of the large-scale magnetic field ranges from about 70 uG to 200 uG in these clumps. The intensity increases towards the Galactic plane. The mass-to-magnetic flux ratio is much smaller than unity, implying that Lupus I is magnetically supported on large scales.
Classical novae eruptions are possible sources of lithium formation and gamma-rays emission. The remnant systems of novae eruptions can also become Type Ia supernovae. The contribution of novae to these phenomena depends on nova rates that are not well established for the Galaxy. Here, we directly measure the Galactic bulge nova rate of $13.9 \pm 2.6$ yr$^{-1}$. This measurement is much more accurate than any previous measurement of this kind thanks to multiple years of bulge monitoring by the OGLE survey. Our sample consists of 39 novae eruptions, $\sim$1/3 of which are OGLE-based discoveries. The long-term monitoring allows us to not only measure the nova rate but also to study in detail the light curves of 39 eruptions and over 80 post-nova candidates. We measured orbital periods for 9 post-novae and 9 novae, in 14 cases we procured the first estimates. The OGLE survey is very sensitive to the frequently erupting recurrent novae. We did not found any object similar to M31 2008-12a, which erupts once a year. The lack of detection indicates a very small number of them in the Galactic bulge.
We present an analytic model to estimate the capabilities of space missions dedicated to the search for biosignatures in the atmosphere of rocky planets located in the habitable zone of nearby stars. Relations between performance and mission parameters such as mirror diameter, distance to targets, and radius of planets, are obtained. Two types of instruments are considered: coronagraphs observing in the visible, and nulling interferometers in the thermal infrared. Missions considered are: single-pupil coronagraphs with a 2.4 m primary mirror, and formation flying interferometers with 4 x 0.75 m collecting mirrors. The numbers of accessible planets are calculated as a function of {\eta}earth. When Kepler gives its final estimation for {\eta}earth, the model will permit a precise assessment of the potential of each instrument. Based on current estimations, {\eta}earth = 10% around FGK stars and 50% around M stars, the coronagraph could study in spectroscopy only ~1.5 relevant planets, and the interferometer ~14.0. These numbers are obtained under the major hypothesis that the exozodiacal light around the target stars is low enough for each instrument. In both cases, a prior detection of planets is assumed and a target list established. For the long-term future, building both types of spectroscopic instruments, and using them on the same targets, will be the optimal solution because they provide complementary information. But as a first affordable space mission, the interferometer looks the more promising in term of biosignature harvest.
We study the impact of possible spiral-arm distributions of Galactic
cosmic-ray sources on the flux of various cosmic-ray nuclei throughout our
Galaxy. We investigate model cosmic-ray spectra at the nominal position of the
sun and at different positions within the Galaxy. The modelling is performed
using the recently introduced numerical cosmic ray propagation code
\textsc{Picard}. Assuming non-axisymmetric cosmic ray source distributions
yields new insights on the behaviour of primary versus secondary nuclei.
We find that primary cosmic rays are more strongly confined to the vicinity
of the sources, while the distribution of secondary cosmic rays is much more
homogeneous compared to the primaries. This leads to stronger spatial variation
in secondary to primary ratios when compared to axisymmetric source
distribution models. A good fit to the cosmic-ray data at Earth can be
accomplished in different spiral-arm models, although leading to decisively
different spatial distributions of the cosmic-ray flux. This results in very
different cosmic ray anisotropies, where even a good fit to the data becomes
possible. Consequently, we advocate directions to seek best fit propagation
parameters that take into account the higher complexity introduced by the
spiral-arm structure on the cosmic-ray distribution. We specifically
investigate whether the flux at Earth is representative for a large fraction of
the Galaxy. The variance among possible spiral-arm models allows us to quantify
the spatial variation of the cosmic-ray flux within the Galaxy in presence of
non-axisymmetric source distributions.
Fast brightness variations are a unique tool to probe the innermost regions of active galactic nuclei (AGN). These variations are called microvariability or intra-night variability, and this phenomenon has been monitored in samples of blazars and unobscured AGNs. Detecting optical microvariations in targets hidden by the obscuring torus is a challenging task because the region responsible for the variations is hidden from our sight. However, there have been reports of fast variations in obscured Seyfert galaxies in X-rays, which rises the question whether microvariations can also be detected in obscured AGNs in the optical regime. Because the expected variations are very small and can easily be lost within the noise, the analysis requires a statistical approach. We report the use of a one-way analysis of variance, ANOVA, with which we searched for microvariability. ANOVA was successfully employed in previous studies of unobscured AGNs. As a result, we found microvariable events during three observing blocks: in two we observed the same object (Mrk 477), and in another, J0759+5050. The results on Mrk 477 confirm previous findings. However, since Mrk 477 is quite a peculiar target with hidden broad-line regions, we cannot rule out the possibility that we have serendipitously chosen a target prone to variations.
Using data from the DEEP2 galaxy redshift survey and the All Wavelength Extended Groth Strip International Survey we obtain stacked X-ray maps of galaxies at 0.7 < z < 1.0 as a function of stellar mass. We compute the total X-ray counts of these galaxies and show that in the soft band (0.5--2,kev) there exists a significant correlation between galaxy X-ray counts and stellar mass at these redshifts. The best-fit relation between X-ray counts and stellar mass can be characterized by a power law with a slope of 0.58 +/- 0.1. We do not find any correlation between stellar mass and X-ray luminosities in the hard (2--7,kev) and ultra-hard (4--7,kev) bands. The derived hardness ratios of our galaxies suggest that the X-ray emission is degenerate between two spectral models, namely point-like power-law emission and extended plasma emission in the interstellar medium. This is similar to what has been observed in low redshift galaxies. Using a simple spectral model where half of the emission comes from power-law sources and the other half from the extended hot halo we derive the X-ray luminosities of our galaxies. The soft X-ray luminosities of our galaxies lie in the range 10^39-8x10^40, ergs/s. Dividing our galaxy sample by the criteria U-B > 1, we find no evidence that our results for X-ray scaling relations depend on optical color.
Using a flexible galactic model with variable stellar velocity anisotropy, I apply the classical Jeans mass modeling approach to the five dwarf spheroidal galaxies with the largest homogeneous datasets of stellar line-of-sight velocities (between 330 and 2500 stars per galaxy) -- Carina, Fornax, Leo~I, Sculptor, and Sextans. I carry out an exhaustive model parameter search, assigning absolute probabilities to each parameter combination. My main finding is that there is a well defined radius (unique for each galaxy) where the Jeans analysis constraints on the enclosed mass are tightest, and are much better than the constraints at previously suggested radii (e.g. 300~pc). For Carina, Fornax, Leo~I, Sculptor, and Sextans the enclosed DM mass is $0.94 \pm 0.20$ (at 410~pc), $7.1 \pm 0.9$ (at 925~pc), $1.75 \pm 0.20$ (at 390~pc), $2.59 \pm 0.42$ (at 435~pc), and $2.3 \pm 0.9$ (at 1035~pc), respectively (two-sigma uncertainties; in $10^7$~M$_\odot$ units). Local DM density has the tightest constraints at smaller (and also unique for each galaxy) radii. The largest central DM density constraint is for Sculptor: $\rho_0\gtrsim 0.09$~$M_\odot$~pc$^{-3}$ (at two-sigma level). I show that the DM density logarithmic slope is totally unconstrained by the Jeans analysis at all the radii probed by the data (and not just at the center, as was demonstrated before). Stellar velocity anisotropy has only very weak constraints. In particular, pure central tangential anisotropy is ruled out at better than two sigma level for three dwarfs, and the data are consistent with the global stellar velocity isotropy for all the five galaxies.
Molecular hydrogen is the most abundant molecule in the Universe. It is thought that a large portion of H2 forms by association of hydrogen atoms to polycyclic aromatic hydrocarbons (PAHs). We model the influence of PAHs on total H2 formation rates in photodissociation regions (PDRs) and assess the effect of these formation rates on the total cloud structure. We set up a chemical kinetic model at steady state in a PDR environment and included adiative transfer to calculate the chemistry at different depths in the PDR. This model includes known dust grain chemistry for the formation of H2 and a H2 formation mechanism on PAHs. Since H2 formation on PAHs is impeded by thermal barriers, this pathway is only efficient at higher temperatures (T > 200 K). At these temperatures the conventional route of H2 formation via H atoms physisorbed on dust grains is no longer feasible, so the PAH mechanism enlarges the region where H2 formation is possible. We find that PAHs have a significant influence on the structure of PDRs. The extinction at which the transition from atomic to molecular hydrogen occurs strongly depends on the presence of PAHs, especially for PDRs with a strong external radiation field. A sharp spatial transition between fully dehydrogenated PAHs on the outside of the cloud and normally hydrogenated PAHs on the inside is found. As a proof of concept, we use coronene to show that H2 forms very efficiently on PAHs, and that this process can reproduce the high H2 formation rates derived in several PDRs.
We report the identification and study of an unusual soft state of the black hole low-mass X-ray binary GRO J1655-40, observed during its 2005 outburst by the Rossi X-ray Timing Explorer. Chandra X-ray grating observations have revealed a high mass-outflow accretion disc wind in this state, and we show that the broadband X-ray spectrum is remarkably similar to that observed in the so-called `hypersoft' state of the high mass X-ray binary Cyg X-3, which possesses a strong stellar wind from a Wolf-Rayet secondary. The power-spectral density (PSD) of GRO J1655-40 shows a bending power-law shape, similar to that of canonical soft states albeit with larger fractional rms. However, the characteristic bend-frequency of the PSD is strongly correlated with the X-ray flux, such that the bend-frequency increases by two decades for less than a factor 2 increase in flux. The strong evolution of PSD bend-frequency for very little change in flux or X-ray spectral shape seems to rule out the suppression of high-frequency variability by scattering in the wind as the origin of the PSD bend. Instead, we suggest that the PSD shape is intrinsic to the variability process and may be linked to the evolution of the scale-height in a slim disc. An alternative possibility is that variability is introduced by variable absorption and scattering in the wind. We further argue that the hypersoft state in GRO J1655-40 and Cyg X-3 is associated with accretion close to or above the Eddington limit.
We have developed a modified form of the equations of smoothed particle magnetohydrodynamics which are stable in the presence of very steep density gradients. Using this formalism, we have performed simulations of the collapse of magnetised molecular cloud cores to form protostars and drive outflows. Our stable formalism allows for smaller sink particles (< 5 AU) than used previously and the investigation of the effect of varying the angle, {\theta}, between the initial field axis and the rotation axis. The nature of the outflows depends strongly on this angle: jet-like outflows are not produced at all when {\theta} > 30{\deg}, and a collimated outflow is not sustained when {\theta} > 10{\deg}. No substantial outflows of any kind are produced when {\theta} > 60{\deg}. This may place constraints on the geometry of the magnetic field in molecular clouds where bipolar outflows are seen.
We have discovered a luminous light echo around the normal Type II-Plateau Supernova (SN) 2012aw in Messier 95 (M95; NGC 3351), detected in images obtained approximately two years after explosion with the Wide Field Channel 3 on-board the Hubble Space Telescope (HST) by the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). The multi-band observations span from the near-ultraviolet through the optical (F275W, F336W, F438W, F555W, and F814W). The apparent brightness of the echo at the time was ~21--22 mag in all of these bands. The echo appears circular, although less obviously as a ring, with an inhomogeneous surface brightness, in particular, a prominent enhanced brightness to the southeast. The SN itself was still detectable, particularly in the redder bands. We are able to model the light echo as the time-integrated SN light scattered off of diffuse interstellar dust in the SN environment. We have assumed that this dust is analogous to that in the Milky Way with R_V=3.1. The SN light curves that we consider also include models of the unobserved early burst of light from the SN shock breakout. Our analysis of the echo suggests that the distance from the SN to the scattering dust elements along the echo is ~45 pc. The implied visual extinction for the echo-producing dust is consistent with estimates made previously from the SN itself. Finally, our estimate of the SN brightness in F814W is fainter than that measured for the red supergiant star at the precise SN location in pre-SN images, possibly indicating that the star has vanished and confirming it as the likely SN progenitor.
We review our state of knowledge of coronal element abundance anomalies in
the Sun and stars. We concentrate on the first ionization potential (FIP)
effect observed in the solar corona and slow-speed wind, and in the coronae of
solar-like dwarf stars, and the "inverse FIP" effect seen in the corona of
stars of later spectral type; specifically M dwarfs. These effects relate to
the enhancement or depletion, respectively, in coronal abundance with respect
to photospheric values of elements with FIP below about 10~eV. They are
interpreted in terms of the ponderomotive force due to the propagation and/or
reflection of magnetohydrodynamic waves in the chromosphere. This acts on
chromospheric ions, but not neutrals, and so can lead to ion-neutral
fractionation.
A detailed description of the model applied to closed magnetic loops, and to
open field regions is given, accounting for the observed difference in solar
FIP fractionation between the slow and fast wind. It is shown that such a model
can also account for the observed depletion of helium in the solar wind. The
helium depletion is sensitive to the chromospheric altitude where ion-neutral
separation occurs, and the behavior of the helium abundance in the closed
magnetic loop strongly suggests that the waves have a coronal origin. This, and
other similar inferences may be expected to have a strong bearing on theories
of solar coronal heating.
Chromospheric waves originating from below as acoustic waves mode convert,
mainly to fast mode waves, can also give rise to ion-neutral separation.
Depending on the geometry of the magnetic field, this can result in FIP or
Inverse FIP effects. We argue that such configurations are more likely to occur
in later-type stars (known to have stronger field in any case), and that this
explains the occurrence of the Inverse FIP effect in M dwarfs.
We consider the possibility that the universe is made of a single dark fluid described by a logotropic equation of state $P=A\ln(\rho/\rho_*)$, where $\rho$ is the rest-mass density, $\rho_*$ is a reference density, and $A$ is the logotropic temperature. The energy density $\epsilon$ is the sum of two terms: a rest-mass energy term $\rho c^2$ that mimics dark matter and an internal energy term $u(\rho)=-P(\rho)-A$ that mimics dark energy. This decomposition leads to a natural, and physical, unification of dark matter and dark energy, and elucidates their mysterious nature. The logotropic model depends on a single parameter $B=A/\rho_{\Lambda}c^2$ where $\rho_{\Lambda}$ is the cosmological density. For $B=0$, we recover the $\Lambda$CDM model. Using cosmological constraints, we find that $0\le B\le 0.09425$. We consider the possibility that dark matter halos are described by the same logotropic equation of state. When $B>0$, pressure gradients prevent gravitational collapse and provide halo density cores instead of cuspy density profiles, in agreement with the observations. The universal rotation curve of logotropic dark matter halos is consistent with the observational Burkert profile up to the halo radius. Interestingly, if we assume that all the dark matter halos have the same logotropic temperature $B$, we find that their surface density $\Sigma=\rho_0 r_h$ is constant. This result is in agreement with the observations where it is found that $\Sigma_0=141\, M_{\odot}/{\rm pc}^2$ for dark matter halos differing by several orders of magnitude in size. Using this observational result, we obtain $B=3.53\times 10^{-3}$. Assuming that $\rho_*=\rho_P$, where $\rho_P$ is the Planck density, we predict $B=3.53\times 10^{-3}$, in perfect agreement with the value obtained from the observations.
The metric outside an isolated object made up of ordinary matter is bound to be the classical Schwarzschild vacuum solution of General Relativity. Nevertheless, some solutions are known (e.g. Morris-Thorne wormholes) that do not match Schwarzschild asymptotically. On a phenomenological point of view, gravitational lensing in metrics falling as $1/r^q$ has recently attracted great interest. In this work, we explore the conditions on the source matter for constructing static spherically symmetric metrics exhibiting an arbitrary power-law as Newtonian limit. For such space-times we also derive the expressions of gravitational redshift and force on probe masses, which, together with light deflection, can be used in astrophysical searches of non-Schwarzschild objects made up of exotic matter. Interestingly, we prove that even a minimally coupled scalar field with a power-law potential can support non-Schwarzschild metrics with arbitrary asymptotic behaviour.
Inflationary models including vector fields have attracted a great deal of attention over the past decade. Such an interest owes to the fact that they might contribute to, or even be fully responsible for, the curvature perturbation imprinted in the CMB. However, the necessary breaking of the vector field's conformal invariance during inflation is not without problems. In recent years it has been realized that a number of instabilities endangering the consistency of the theory arise when the conformal invariance is broken by means of a non-minimal coupling to gravity. In this paper we consider a massive vector field non-minimally coupled to gravity through the Gauss-Bonnet invariant, and investigate whether the vector can obtain a nearly scale-invariant perturbation spectrum while evading the emergence of perturbative instabilities. We find that the strength of the coupling must be extremely small if the vector field is to have a chance to contribute to the total curvature perturbation.
We discuss the inviscid 2D instability recently uncovered by Ilin & Morgulis (2013) in the context of irrotational Taylor-Couette flow with a radial flow imposed. By finding a simplier rectilinear example of the instability - the sheared half plane, the minimal ingredients for the instability are identified and the destabilizing/stabilizing effect of the inflow/outflow boundaries clarified. The instability - christened `boundary inflow instability' here - is of critical layer type where this layer is either at the inflow wall and the growth rate is $O(\eta^{1/2})$ (as found by Ilin & Morgulis 2013), or in the interior of the flow and the growth rate is $O(\eta \log(1/\eta) )$ where $\eta$ measures the (small) inflow-to-tangential-flow ratio. The instability is robust to changes in the rotation profile even to those which are very Rayleigh-stable and the addition of further physics such as viscosity, 3-dimensionality and compressibility but is sensitive to the boundary condition imposed on the tangential velocity field at the inflow boundary. Both the primary bifurcation to 2D states and secondary bifurcations to 3D states are found to be supercritical. Assuming an accretion flow driven by molecular viscosity only so $\eta=O(Re^{-1})$, the instability is not immediately relevant for accretion disks since the critical threshold is $O(Re^{-2/3})$ and the inflow boundary conditions are more likely to be stress-free than non-slip. However, the analysis presented here does highlight the potential for mass entering a disk to disrupt the orbiting flow if this mass flux possesses vorticity.
It is shown whether Higgs inflation can be saved with high-scale supersymmetry critically depends on the magnitude of non-minimal coupling constant $\xi$. For small $\xi \leq 500$, the threshold correction at scale $M_{P}/\xi$ is constrained in high precision.Its magnitude is in the narrow range of $(-0.03, -0.02)$ and $(-0.05, -0.04)$ for the wino and higgsino/singlino dark matter, respectively. While in the large $\xi$-region with $\xi \geq 10^{4}$, such high-scale supersymmetry is excluded by too large threshold correction as required by Higgs inflation.
In this paper we seek for relevant information on the asymptotic cosmological dynamics of the Brans--Dicke theory of gravity for several self-interaction potentials. By means of the simplest tools of the dynamical systems theory, it is shown that the general relativity de Sitter solution is an attractor of the Jordan frame (dilatonic) Brans--Dicke theory only for the exponential potential $U(\vphi)\propto\exp\vphi$, which corresponds to the quadratic potential $V(\phi)\propto\phi^2$ in terms of the original Brans--Dicke field $\phi=\exp\vphi$, or for potentials which asymptote to $\exp\vphi$. At the stable de Sitter critical point, as well as at the stiff-matter equilibrium configurations, the dilaton is necessarily massless. We find bounds on the Brans--Dicke coupling constant $\omega_\textsc{bd}$, which are consistent with well-known results.
Links to: arXiv, form interface, find, astro-ph, recent, 1504, contact, help (Access key information)