1902.03252
High performance computing for gravitational lens modeling: single vs double precision on GPUs and CPUs
Rexroth, et al
Strong gravitational lensing is a powerful probe of cosmology and the dark matter distribution. Efficient lensing software is already a necessity to fully use its potential and the performance demands will only increase with the upcoming generation of telescopes. In this paper, we study the possible impact of High Performance Computing techniques on a performance-critical part of the widely used lens modeling software LENSTOOL. We implement the algorithm once as a highly optimized CPU version and once with graphics card acceleration for a simple parametric lens model. In addition, we study the impact of finite machine precision on the lensing algorithm. While double precision is the default choice for scientific applications, we find that single precision can be sufficiently accurate for our purposes and lead to a big speedup. Therefore we develop and present a mixed precision algorithm which only uses double precision when necessary. We measure the performance of the different implementations and find that the use of High Performance Computing Techniques dramatically improves the code performance both on CPUs and GPUs. Compared to the current LENSTOOL implementation on 12 CPU cores, we obtain speedup factors of up to 170. We achieve this optimal performance by using our mixed precision algorithm on a high-end GPU which is common in modern supercomputers. We also show that these techniques reduce the energy consumption by up to 98%. Furthermore, we demonstrate that a highly competitive speedup can be reached with consumer GPUs. While they are an order of magnitude cheaper than the high-end graphics cards, they are rarely used for scientific computations due to their low double precision performance. Our mixed precision algorithm unlocks their full potential. The consumer GPU delivers a speedup which is only a factor of four lower than the best speedup achieved by a high-end GPU.
1902.03261
Type Ia Supernovae are Excellent Standard Candles in the Near-Infrared
Avelino, et al
1902.03523
Propelling interplanetary spacecraft utilizing water-steam
Martinez, Thangavelautham
Water has been identified as a critical resource both to sustain human-life but also for use in propulsion, attitude-control, power, thermal and radiation pro-tection systems. Water may be obtained off-world through In-Situ Resource Utilization (ISRU) in the course of human or robotic space exploration that replace materials that would otherwise be shipped from Earth. Water has been highlighted by many in the space community as a credible solution for affordable/sustainable exploration. Water can be extracted from the Moon, C-class Near Earth Objects (NEOs), surface of Mars and Martian Moons Phobos and Deimos and from the surface of icy, rugged terrains of Ocean Worlds. However, use of water for propulsion faces some important technological barriers. A technique to use water as a propellant is to electrolyze it into hydrogen and oxygen that is then pulse-detonated. High-efficiency electrolysis requires use of platinum-catalyst based fuel cells. Even trace elements of sulfur and carbon monoxide found on planetary bodies can poison these cells making them unusable. In this work, we develop steam-based propulsion that avoids the technological barriers of electrolyzing impure water as propellant. Using a solar concentrator, heat is used to extract the water which is then condensed as a liquid and stored. Steam is then formed using the solar thermal reflectors to concentrate the light into a nanoparticle-water mix. This solar thermal heating (STH) process converts 80 to 99% of the in-coming light into heat.
1902.03624
QSOs acting as gravitational lenses: halo mass and projected mass density profile at $z\sim0.7$
Bonavera, et al
Magnification bias is a gravitational lensing effect that is normally overlooked because it is considered sub-optimal in comparison with the lensing shear. Thanks to the demonstrated optimal characteristics of the sub-millimetre galaxies (SMGs) for lensing analysis, in this work we were able to measure the magnification bias produced by a sample of QSOs acting as lenses, $0.2<z<1.0$, on the SMGs observed by Herschel at $1.2<z<4.0$. Two different methodologies were successfully applied: the traditional cross-correlation function approach and the stacking technique. The second one was found to be more robust for analysing the strong lensing regime ($<20-30$ arcsec in our case) and provides, in addition, a density map of the lensing effect. From the halo modelling of the cross-correlation function, the QSOs host halo mass was estimated to be greater than $\log_{10}{(M_{min}/M_\odot)} > 13.6_{-0.4}^{+0.9}$, also confirmed by the mass density profile analysis ($M_{NFW}=1.6_{-0.5}^{+2.1}\times10^{14} M_\odot$). These mass values indicate that we are observing the lensing effect of a cluster size halo signposted by the QSOs. Moreover, we were able to estimate the lensing convergence, $\kappa(\theta)$, for our magnification bias measurements down to a few kpcs. The derived mass density profile shows a core in its centre and it is compatible with a Navarro-Frank-White (NFW) profile. The core radius, $r_c=30_{-10}^{+14}$ kpc, was estimated using a cored Singular Isothermal Sphere profile that fits even better the data. In addition, we were not able to detect any baryonic/stellar component at very small scales. Both results could indicate the effect of a strong baryonic feedback produced by the active galactic nuclei in the centre of the QSOs responsible of the removal of the expected cusp.
1902.03928
We report on the relation between the mass of supermassive black holes (SMBHs, M_BH) and that of hosting dark matter haloes (M_h) for 49 z ~ 6 QSOs with [CII]158um velocity-width measurements. Here, we estimate M_h assuming that the rotation velocity from FWHM_CII is equal to the circular velocity of the halo; we have tested this procedure using z ~ 3 QSOs which also have clustering-based M_h estimates. We find that a vast majority of the z ~ 6 SMBHs are more massive than expected from the local M_BH - M_h relation, with one third of the sample by factors >~ 10^2. The median mass ratio of the sample, M_BH/M_h = 6 x 10^{-4}, means that 0.4% of the baryons in haloes are locked up in SMBHs. The mass growth rates of our SMBHs amount to ~ 10% of the SFRs, or ~ 1% of the mean baryon accretion rates, of the hosting galaxies. A large fraction of the hosting galaxies are consistent with average galaxies in terms of SFR and perhaps of stellar mass and size. Our study indicates that the growth of SMBHs (M_BH ~ 10^{8-10} Msun) in luminous z ~ 6 QSOs greatly precedes that of hosting haloes owing to efficient gas accretion even under normal star formation activities, although we cannot rule out the possibility that undetected SMBHs have local M_BH/M_h ratios. This preceding growth is in contrast to much milder evolution of the stellar-to-halo mass ratio.
During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|\beta_{He}-\beta_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before.
The Local Group contains a great number of dwarf irregulars and spheroidals, for which the spectroscopy of individual stars can be obtained. Thus, the chemical evolution of these galaxies can be traced, with the only need of finding populations spanning a large age range and such that we can accurately derive the composition. Planetary nebulae (PNe) are old- and intermediate-age star remnants and their chemical abundances can be obtained up to 3-4 Mpc. H II regions, which are brighter and much easily detected, represent galaxies young content. PNe and H II regions share similar spectroscopic features and are analysed in the same way. Both are among the best tracers of the chemical evolution allowing to draw the chemical time line of nearby galaxies. The focus in this review are the PN and H II region populations as constraints to the chemical evolution models and the mass-metallicity relation of the local universe.
We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We investigate the impact of different magneto-genesis scenarios, both, primordial and astrophysical, on the propagation of cosmic rays. Our study shows that different scenarios of magneto-genesis do not have a large impact on the anisotropy measurements. The distribution of nearby sources causes anisotropy at very high energies, independent of the magnetic field model. We compare our results to the dipole signal measured by the Pierre Auger Observatory. All our models could reproduce the observed dipole amplitude with a pure iron injection composition. This is due to clustering of secondary nuclei in direction of nearby sources of heavy nuclei. A light injection composition is disfavoured by the non-observation of anisotropy at energies of 4 - 8 EeV.
One of the most universal trends in science and technology today is the growth of large teams in all areas, as solitary researchers and small teams diminish in prevalence1,2,3. Increases in team size have been attributed to the specialization of scientific activities3, improvements in communication technology4,5, or the complexity of modern problems that require interdisciplinary solutions6,7,8. This shift in team size raises the question of whether and how the character of the science and technology produced by large teams differs from that of small teams. Here we analyse more than 65 million papers, patents and software products that span the period 1954–2014, and demonstrate that across this period smaller teams have tended to disrupt science and technology with new ideas and opportunities, whereas larger teams have tended to develop existing ones. Work from larger teams builds on more-recent and popular developments, and attention to their work comes immediately. By contrast, contributions by smaller teams search more deeply into the past, are viewed as disruptive to science and technology and succeed further into the future—if at all. Observed differences between small and large teams are magnified for higher-impact work, with small teams known for disruptive work and large teams for developing work. Differences in topic and research design account for a small part of the relationship between team size and disruption; most of the effect occurs at the level of the individual, as people move between smaller and larger teams. These results demonstrate that both small and large teams are essential to a flourishing ecology of science and technology, and suggest that, to achieve this, science policies should aim to support a diversity of team sizes.
In the laboratory, the photo-and thermochemical evolution of ices, made of simple molecules of astrophysical relevance, always leads to the formation of semi-refractory water-soluble organic residues. Targeted searches for specific molecules do reveal the notable presence of two families of important molecular ''bricks of life'': amino acids, key molecules in metabolism, and sugars, including ribose, the backbone of RNA molecules which support the genetic information in all living entities. Most of these molecules are indeed found in primitive carbonaceous meteorites and their implication in prebiotic chemistry at the surface of the early Earth must be seriously considered. These molecules are, almost all, chiral. In meteorites, some amino acids do show significant enantiomeric excesses, practically exclusively of the L-form. In our experiments, we investigate the role of circularly polarized light obtained from the DESIRS beamline of the synchrotron SOLEIL, a light commonly observed in regions of star formation, in order to generate an initial symmetry breaking in chiral amino acids produced and then indeed detected in our samples. We present first a brief global description of the chemical evolution of the Galaxy. Then, using our laboratory simulations, we suggest the importance of cosmic ices in the build-up of complex organic matter, including enantioenrichment at the surface of telluric planets like the Earth, thus establishing a link between astrochemistry and astrobiology.
Our statistical understanding of galaxy evolution is fundamentally driven by objects that lie above the surface-brightness limits of current wide-area surveys (mu ~ 23 mag arcsec^-2). While both theory and small, deep surveys have hinted at a rich population of low-surface-brightness galaxies (LSBGs) fainter than these limits, their formation remains poorly understood. We use Horizon-AGN, a cosmological hydrodynamical simulation to study how LSBGs, and in particular the population of ultra-diffuse galaxies (UDGs; mu > 24.5 mag arcsec^-2), form and evolve over time. For M* > 10^8 MSun, LSBGs contribute 47, 7 and 6 per cent of the local number, mass and luminosity densities respectively (~85/11/10 per cent for M* > 10^7 MSun). Today's LSBGs have similar dark-matter fractions and angular momenta to high-surface-brightness galaxies (HSBGs; mu < 23 mag arcsec^-2), but larger effective radii (x2.5 for UDGs) and lower fractions of dense, star-forming gas (more than x6 less in UDGs than HSBGs). LSBGs originate from the same progenitors as HSBGs at z > 2. However, LSBG progenitors form stars more rapidly at early epochs. The higher resultant rate of supernova-energy injection flattens their gas-density profiles, which, in turn, creates shallower stellar profiles that are more susceptible to tidal processes. After z ~ 1, tidal perturbations broaden LSBG stellar distributions and heat their cold gas, creating the diffuse, largely gas-poor LSBGs seen today. In clusters, ram-pressure stripping provides an additional mechanism that assists in gas removal in LSBG progenitors. Our results offer insights into the formation of a galaxy population that is central to a complete understanding of galaxy evolution, and which will be a key topic of research using new and forthcoming deep-wide surveys.
1902.04585
Cosmic voids uncovered -- first-order statistics of depressions in the biased density field
Ronconi, et al
Cosmic voids are the major volume component in the matter distribution of the Universe. They posses great potential for constraining dark energy as well as for testing theories of gravity. Nevertheless, in spite of their growing popularity as cosmological probes, a gap of knowledge between cosmic void observations and theory still persists. In particular, the void size function models proposed in literature have been proven unsuccessful in reproducing the results obtained from cosmological simulations in which cosmic voids are detected from biased tracers of the density field, undermining the possibility of using them as cosmological probes. The goal of this work is to cover this gap. In particular, we make use of the findings of a previous work in which we have improved the void selection procedure, presenting an algorithm that redefines the void ridges and, consequently, their radius. By applying this algorithm, we validate the volume conserving model of the void size function on a set of unbiased simulated density field tracers. We highlight the difference in the internal structure between voids selected in this way and those identified by the popular VIDE void finder. We also extend the validation of the model to the case of biased tracers. We find that a relation exists between the tracer used to sample the underlying dark matter density field and its unbiased counterpart. Moreover, we demonstrate that, as long as this relation is accounted for, the size function is a viable approach for studying cosmology with cosmic voids.
1902.03663
Weak lensing cosmology with convolutional neural networks on noisy data
Ribli, et al
Weak gravitational lensing is one of the most promising cosmological probes of the late universe. Several large ongoing (DES, KiDS, HSC) and planned (LSST, EUCLID, WFIRST) astronomical surveys attempt to collect even deeper and larger scale data on weak lensing. Due to gravitational collapse, the distribution of dark matter is non-Gaussian on small scales. However, observations are typically evaluated through the two-point correlation function of galaxy shear, which does not capture non-Gaussian features of the lensing maps. Previous studies attempted to extract non-Gaussian information from weak lensing observations through several higher-order statistics such as the three-point correlation function, peak counts or Minkowski-functionals. Deep convolutional neural networks (CNN) emerged in the field of computer vision with tremendous success, and they offer a new and very promising framework to extract information from 2 or 3-dimensional astronomical data sets, confirmed by recent studies on weak lensing. We show that a CNN is able to yield significantly stricter constraints of ($\sigma_8, \Omega_m$) cosmological parameters than the power spectrum using convergence maps generated by full N-body simulations and ray-tracing, at angular scales and shape noise levels relevant for future observations. In a scenario mimicking LSST or Euclid, the CNN yields 2.4-2.8 times smaller credible contours than the power spectrum, and 3.5-4.2 times smaller at noise levels corresponding to a deep space survey such as WFIRST. We also show that at shape noise levels achievable in future space surveys the CNN yields 1.4-2.1 times smaller contours than peak counts, a higher-order statistic capable of extracting non-Gaussian information from weak lensing maps.
Weak lensing cosmology with convolutional neural networks on noisy data
Ribli, et al
Weak gravitational lensing is one of the most promising cosmological probes of the late universe. Several large ongoing (DES, KiDS, HSC) and planned (LSST, EUCLID, WFIRST) astronomical surveys attempt to collect even deeper and larger scale data on weak lensing. Due to gravitational collapse, the distribution of dark matter is non-Gaussian on small scales. However, observations are typically evaluated through the two-point correlation function of galaxy shear, which does not capture non-Gaussian features of the lensing maps. Previous studies attempted to extract non-Gaussian information from weak lensing observations through several higher-order statistics such as the three-point correlation function, peak counts or Minkowski-functionals. Deep convolutional neural networks (CNN) emerged in the field of computer vision with tremendous success, and they offer a new and very promising framework to extract information from 2 or 3-dimensional astronomical data sets, confirmed by recent studies on weak lensing. We show that a CNN is able to yield significantly stricter constraints of ($\sigma_8, \Omega_m$) cosmological parameters than the power spectrum using convergence maps generated by full N-body simulations and ray-tracing, at angular scales and shape noise levels relevant for future observations. In a scenario mimicking LSST or Euclid, the CNN yields 2.4-2.8 times smaller credible contours than the power spectrum, and 3.5-4.2 times smaller at noise levels corresponding to a deep space survey such as WFIRST. We also show that at shape noise levels achievable in future space surveys the CNN yields 1.4-2.1 times smaller contours than peak counts, a higher-order statistic capable of extracting non-Gaussian information from weak lensing maps.
1902.03786
A universal constant for dark matter-baryon interplay
Chan
Recent studies point out that there exists some rough scaling relations for dark matter and some tight connections between dark matter and baryons. However, most of the relations and tight connections can only be found in galaxies, but not in galaxy clusters. In this article, we consider a new expression that can characterize the properties of dark matter-baryon interplay for both galactic and galaxy cluster scales. By using the archival observational data of galaxies and galaxy clusters, we show that the value $K=\bar{n}_D \bar{n}_Br_oV/v^4$ is almost a constant and scale independent within the optical radius $r_o$, where $\bar{n}_D$ is the average dark matter number density, $\bar{n}_B$ is the average baryon number density, $v$ is the characteristic velocity and $V$ is the interacting volume. This would be the first universal relation between dark matter and baryons on both galactic and galaxy cluster scales. We anticipate this result to be a starting point to explain the small-scale problem and the scaling relations for dark matter in galaxies. The constant $K$ discovered may reveal some underlying global interaction between dark matter and baryons.
1902.03928
The degree of fine-tuning in our Universe -- and Others
Adams
Both fundamental constants that describe the laws of physics and cosmological parameters that determine the cosmic properties must fall within a range of values in order for the universe to develop astrophysical structures and ultimately support life. This paper reviews current constraints on these quantities. The standard model of particle physics contains both coupling constants and particle masses, and the allowed ranges of these parameters are discussed first. We then consider cosmological parameters, including the total energy density, the vacuum energy density, the baryon-to-photon ratio, the dark matter contribution, and the amplitude of primordial density fluctuations. These quantities are constrained by the requirements that the universe lives for a long time, emerges from the BBN epoch with an acceptable chemical composition, and can successfully produce galaxies. On smaller scales, stars and planets must be able to form and function. The stars must have sufficiently long lifetimes and hot surface temperatures. The planets must be massive enough to maintain an atmosphere, small enough to remain non-degenerate, and contain enough particles to support a complex biosphere. These requirements place constraints on the gravitational constant, the fine structure constant, and composite parameters that specify nuclear reaction rates. We consider specific instances of possible fine-tuning in stars, including the triple alpha reaction that produces carbon, as well as the effects of unstable deuterium and stable diprotons. For all of these issues, viable universes exist over a range of parameter space, which is delineated herein. Finally, for universes with significantly different parameters, new types of astrophysical processes can generate energy and support habitability.
1902.04029
Quantifying tension: interpreting the DES evidence ratio
Handley, Lemos
We provide a new interpretation for the Bayes factor combination used in the DES first year analysis to quantify the tension between the Dark Energy Survey (DES) and Planck datasets. The ratio quantifies a Bayesian confidence in our ability to combine the datasets. This interpretation is prior-dependent, with wider prior widths boosting the confidence. We therefore propose that if there are any reasonable priors which reduce the confidence to below unity, then we cannot assert that the datasets are compatible. Computing the evidence ratios for the DES first year analysis and Planck, given that narrower priors drop the confidence to below unity, we conclude that DES and Planck are, in a Bayesian sense, incompatible under $\Lambda$CDM. Additionally we compute ratios which confirm the consensus that measurements of the acoustic scale by the Baryon Oscillation Spectroscopic Survey (BOSS) are compatible with Planck, whilst direct measurements of the acceleration rate of the Universe by the SH0ES collaboration are not. We propose a modification to the Bayes ratio which removes the prior dependency using Kullback-Leibler divergences, and using this statistical test find Planck in strong tension with SH0ES in moderate tension with DES, and in no tension with BOSS. We propose this statistic as the optimal way to compare datasets, ahead of the next DES data releases, as well as future surveys. Finally, as an element of these calculations, we introduce in a cosmological setting the Bayesian model dimensionality, which is a parameterisation-independent measure of the number of parameters that a given dataset constrains.
Quantifying tension: interpreting the DES evidence ratio
Handley, Lemos
We provide a new interpretation for the Bayes factor combination used in the DES first year analysis to quantify the tension between the Dark Energy Survey (DES) and Planck datasets. The ratio quantifies a Bayesian confidence in our ability to combine the datasets. This interpretation is prior-dependent, with wider prior widths boosting the confidence. We therefore propose that if there are any reasonable priors which reduce the confidence to below unity, then we cannot assert that the datasets are compatible. Computing the evidence ratios for the DES first year analysis and Planck, given that narrower priors drop the confidence to below unity, we conclude that DES and Planck are, in a Bayesian sense, incompatible under $\Lambda$CDM. Additionally we compute ratios which confirm the consensus that measurements of the acoustic scale by the Baryon Oscillation Spectroscopic Survey (BOSS) are compatible with Planck, whilst direct measurements of the acceleration rate of the Universe by the SH0ES collaboration are not. We propose a modification to the Bayes ratio which removes the prior dependency using Kullback-Leibler divergences, and using this statistical test find Planck in strong tension with SH0ES in moderate tension with DES, and in no tension with BOSS. We propose this statistic as the optimal way to compare datasets, ahead of the next DES data releases, as well as future surveys. Finally, as an element of these calculations, we introduce in a cosmological setting the Bayesian model dimensionality, which is a parameterisation-independent measure of the number of parameters that a given dataset constrains.
1902.0
1902.04165
Black vs. Dark: rapid growth of supermassive black holes in dark matter haloes at z~6
Shimasaku, Izumi
1902.04193
Test of Einstein equivalence principle near the Galactic center supermassive black hole
GRAVITY Collaboration, et al
1902.04291
What do planetary nebulae and H II regions reveal about the chemical evolution of nearby dwarf galaxies?
Conçalves
1902.04408
Propagation of UHECRs in the local Universe and origin of cosmic magnetic fields
Hackstein, et al
Nature
Large teams develop and small teams disrupt science and technology
Wu, Wang, Evans
1902.04575
Interstellar ices: a possible scenario for symmetry reading of extraterrestrial chiral organic molecules of prebiotic interest
D’hendecourt, et al
In the laboratory, the photo-and thermochemical evolution of ices, made of simple molecules of astrophysical relevance, always leads to the formation of semi-refractory water-soluble organic residues. Targeted searches for specific molecules do reveal the notable presence of two families of important molecular ''bricks of life'': amino acids, key molecules in metabolism, and sugars, including ribose, the backbone of RNA molecules which support the genetic information in all living entities. Most of these molecules are indeed found in primitive carbonaceous meteorites and their implication in prebiotic chemistry at the surface of the early Earth must be seriously considered. These molecules are, almost all, chiral. In meteorites, some amino acids do show significant enantiomeric excesses, practically exclusively of the L-form. In our experiments, we investigate the role of circularly polarized light obtained from the DESIRS beamline of the synchrotron SOLEIL, a light commonly observed in regions of star formation, in order to generate an initial symmetry breaking in chiral amino acids produced and then indeed detected in our samples. We present first a brief global description of the chemical evolution of the Galaxy. Then, using our laboratory simulations, we suggest the importance of cosmic ices in the build-up of complex organic matter, including enantioenrichment at the surface of telluric planets like the Earth, thus establishing a link between astrochemistry and astrobiology.
1902.04580
The formation and evolution of low-surface-brightness galaxies
Martin, et al
1902.04585
Cosmic voids uncovered -- first-order statistics of depressions in the biased density field
Ronconi, et al
Cosmic voids are the major volume component in the matter distribution of the Universe. They posses great potential for constraining dark energy as well as for testing theories of gravity. Nevertheless, in spite of their growing popularity as cosmological probes, a gap of knowledge between cosmic void observations and theory still persists. In particular, the void size function models proposed in literature have been proven unsuccessful in reproducing the results obtained from cosmological simulations in which cosmic voids are detected from biased tracers of the density field, undermining the possibility of using them as cosmological probes. The goal of this work is to cover this gap. In particular, we make use of the findings of a previous work in which we have improved the void selection procedure, presenting an algorithm that redefines the void ridges and, consequently, their radius. By applying this algorithm, we validate the volume conserving model of the void size function on a set of unbiased simulated density field tracers. We highlight the difference in the internal structure between voids selected in this way and those identified by the popular VIDE void finder. We also extend the validation of the model to the case of biased tracers. We find that a relation exists between the tracer used to sample the underlying dark matter density field and its unbiased counterpart. Moreover, we demonstrate that, as long as this relation is accounted for, the size function is a viable approach for studying cosmology with cosmic voids.
1902.04651
Spin evolution and feedback of supermassive black holes in cosmological simulations
Bustamante, Springel
It is well established that the properties of supermassive black holes and their host galaxies are correlated through scaling relations. While hydrodynamical cosmological simulations have begun to account for the co-evolution of BHs and galaxies, they typically have neglected the BH spin, even though it may play an important role in modulating the growth and feedback of BHs. Here we introduce a new sub-grid model for the BH spin evolution in the moving-mesh code {\small AREPO} in order to improve the physical faithfulness of the BH modelling in galaxy formation simulations. We account for several different channels of spin evolution, in particular gas accretion through a Shakura-Sunyaev $\alpha$-disc, chaotic accretion, and BH mergers. For BH feedback, we extend the IllustrisTNG model, which considers two different BH feedback modes, a thermal quasar mode for high accretion states and a kinetic mode for low Eddington ratios, with a self-consistent accounting of spin-dependent radiative efficiencies and thus feedback strength. We find that BHs with mass $M_{\rm{bh}}\lesssim 10^{8}\, {\rm M}_\odot$ reach high spin values as they typically evolve in the coherent gas accretion regime. On the other hand, BHs with mass $M_{\rm{bh}}\gtrsim 10^{8}\, {\rm M}_\odot$ have lower spins as BH mergers become more frequent, and their accretion discs fragment due to self-gravity, inducing chaotic accretion. We also explore the hypothesis that the transition between the quasar and kinetic feedback modes is mediated by the accretion mode of the BH disc itself, i.e.~the kinetic feedback mode is activated when the disc enters the self-gravity regime. We find excellent agreement between the galaxy and BH populations for this approach and the fiducial TNG model with no spin evolution. Furthermore, our new approach alleviates a tension in the galaxy morphology -- colour relation of the original TNG model.
1902.004720
A limited habitable zone for complex life
Schweiterman, et al
The habitable zone (HZ) is defined as the range of distances from a host star within which liquid water, a key requirement for life, may exist at a planet's surface. Substantially more CO2 than present in Earth's modern atmosphere is required to maintain clement temperatures for most of the HZ, with concentrations of several bars required at the outer edge. However, most complex aerobic life on Earth is precluded by CO2 concentrations of just a small fraction of a bar. At the same time, most of the HZ volume resides in proximity to K and M dwarfs, which are more numerous than Sun-like G dwarfs but are predicted to promote greater abundances of CO in the atmospheres of orbiting planets, a highly toxic gas for complex aerobic organisms with circulatory systems. Here we show that the HZ for complex aerobic life is significantly limited relative to that for simple microbial life. We use 1-D radiative-convective climate and photochemical models to circumscribe the Habitable Zone for Complex Life (HZCL) based on known toxicity limits for a range of complex organisms. We find that for CO2 tolerances of 0.005-0.05 bar, the HZCL is only ~20-28% as wide as the traditional HZ for a Sun-like star and that CO concentrations may limit complex life throughout the entire HZ of the coolest M dwarfs. These results cast new light on the likely distribution of complex life in the universe and have important practical ramifications for the search for exoplanet biosignatures and technosignatures.
1902.04872
What determines the shape of the local (z<0.1) infrared galaxy luminosity function?
Symeonidis, Page
We investigate what shapes the infrared luminosity function of local galaxies by comparing it to the local infrared AGN luminosity function. The former corresponds to emission from dust heated by stars and AGN, whereas the latter includes emission from AGN-heated dust only. Our results show that infrared emission from AGN starts mixing into the galaxy luminosity function in the luminous infrared galaxy (LIRG) regime and becomes significant in the ultraluminous infrared galaxy (ULIRG) regime, with the luminosity above which local ULIRGs become AGN-dominated being in the log(LIR/Lsun)~12.2-12.7 range. We propose that as a result of the AGN contribution, the infrared galaxy luminosity function has a flatter high luminosity slope than UV/optical galaxy luminosity functions. Furthermore, we note that the increased AGN contribution as a function of LIR is reflected in the average dust temperature (Tdust) of local galaxies, and may be responsible for the local LIR-Tdust relation. However, although our results show that AGN play a central role in defining the properties of local ULIRGs, we find that the dominant power source in the local ULIRG population is star-formation.
What determines the shape of the local (z<0.1) infrared galaxy luminosity function?
Symeonidis, Page
We investigate what shapes the infrared luminosity function of local galaxies by comparing it to the local infrared AGN luminosity function. The former corresponds to emission from dust heated by stars and AGN, whereas the latter includes emission from AGN-heated dust only. Our results show that infrared emission from AGN starts mixing into the galaxy luminosity function in the luminous infrared galaxy (LIRG) regime and becomes significant in the ultraluminous infrared galaxy (ULIRG) regime, with the luminosity above which local ULIRGs become AGN-dominated being in the log(LIR/Lsun)~12.2-12.7 range. We propose that as a result of the AGN contribution, the infrared galaxy luminosity function has a flatter high luminosity slope than UV/optical galaxy luminosity functions. Furthermore, we note that the increased AGN contribution as a function of LIR is reflected in the average dust temperature (Tdust) of local galaxies, and may be responsible for the local LIR-Tdust relation. However, although our results show that AGN play a central role in defining the properties of local ULIRGs, we find that the dominant power source in the local ULIRG population is star-formation.
1902.04877
Estimating the integrated bispectrum from weak lensing maps
Munshi, et al
1902.05089
Hunting for the dark matter wake induced by the Large Magellanic Cloud
Garavito-Camargo, et al
Satellite galaxies are predicted to generate gravitational density wakes as they orbit within the dark matter (DM) halos of their hosts, causing their orbits to decay over time. The recent infall of the Milky Way's (MW) most massive satellite galaxy, the Large Magellanic Cloud (LMC), affords us the unique opportunity to study this process in action. In this work, we present high-resolution ($m_{dm} = 4 \times 10^4 M_{\odot}$ ) N-body simulations of the MW-LMC interaction over the past 2 Gyr. We quantify the impact of the LMC's passage on the density and kinematics of the MW's DM halo and the observability of these structures in the MW's stellar halo. The LMC is found to generate pronounced Local and Global wakes in both the DM and stellar halos, leads to both local overdensities and distinct kinematic patterns that should be observable with ongoing and future surveys. Specifically, the Global Wake will result in redshifted radial velocities of stars in the North and blueshifts in the South, at distances larger than 45 kpc. The Local Wake traces the orbital path of the LMC through the halo (50-200 kpc), resulting in a stellar overdensity with a distinct, tangential kinematic pattern that persists to the present day. The detection of the MW's halo response will constrain: the infall mass of the LMC and its orbital trajectory, the mass of the MW, and it may inform us about the nature of the dark matter particle itself.
Hunting for the dark matter wake induced by the Large Magellanic Cloud
Garavito-Camargo, et al
Satellite galaxies are predicted to generate gravitational density wakes as they orbit within the dark matter (DM) halos of their hosts, causing their orbits to decay over time. The recent infall of the Milky Way's (MW) most massive satellite galaxy, the Large Magellanic Cloud (LMC), affords us the unique opportunity to study this process in action. In this work, we present high-resolution ($m_{dm} = 4 \times 10^4 M_{\odot}$ ) N-body simulations of the MW-LMC interaction over the past 2 Gyr. We quantify the impact of the LMC's passage on the density and kinematics of the MW's DM halo and the observability of these structures in the MW's stellar halo. The LMC is found to generate pronounced Local and Global wakes in both the DM and stellar halos, leads to both local overdensities and distinct kinematic patterns that should be observable with ongoing and future surveys. Specifically, the Global Wake will result in redshifted radial velocities of stars in the North and blueshifts in the South, at distances larger than 45 kpc. The Local Wake traces the orbital path of the LMC through the halo (50-200 kpc), resulting in a stellar overdensity with a distinct, tangential kinematic pattern that persists to the present day. The detection of the MW's halo response will constrain: the infall mass of the LMC and its orbital trajectory, the mass of the MW, and it may inform us about the nature of the dark matter particle itself.
1902.05140
Discovery of strongly-lensed gravitational waves - implications for the LSST observing strategy
Smith, et al
LSST's wide-field of view and sensitivity will revolutionize studies of the transient sky by finding extraordinary numbers of new transients every night. The recent discovery of a kilonova counterpart to LIGO/Virgo's first detection of gravitational waves (GWs) from a double neutron star (NS-NS) merger also creates an exciting opportunity for LSST to offer a Target of Opportunity (ToO) mode of observing. We have been exploring the possibility of detecting strongly lensed GWs, that would enable new tests of GR, extend multi-messenger astronomy out to $z\gtrsim1$, and deliver a new class of sub-millisecond precision time-delay constraints on lens mass distributions. We forecast that the rate of detection of lensed NS-NS mergers in the 2020s will be $\sim0.1$ per Earth year, that the typical source will be at $z\simeq2$, and that the multiply-imaged kilonova counterpart will have a magnitude of ${\rm AB}\simeq25.4$ in $g/r/i$-band filters - i.e. fainter than the sensitivity of a single LSST WFD visit. We therefore advocate (1) creating a flexible and efficient Target of Opportunity programme within the LSST observing strategy that is capable of discovering sources fainter than single-visit depth, and (2) surveying the entire observable extragalactic sky as rapidly as possible in the WFD survey. The latter will enable a very broad range of early science that relies on wide survey area for detection of large samples of objects and/or maximizing the fraction of sky over which reference imaging is available. For example, it will enable prompt discovery of a uniform and all-sky sample of galaxy/group/cluster-scale lenses that will underpin LSST strong-lensing science. This white paper complements submissions from DESC, SLSC, and TVSSC, that discuss kilonova, GW, and strong lensing.
1902.05141
Strong lensing considerations for the LSST observing strategy
Verma, et al
Strong gravitational lensing enables a wide range of science: probing cosmography; testing dark matter models; understanding galaxy evolution; and magnifying the faint, small and distant Universe. However to date exploiting strong lensing as a tool for these numerous cosmological and astrophysical applications has been severely hampered by limited sample size. LSST will drive studies of strongly lensed galaxies, galaxy groups and galaxy clusters into the statistical age. Time variable lensing events, e.g. measuring cosmological time delays from strongly lensed supernovae and quasars, place the strongest constraints on LSST's observing strategy and have been considered in the DESC observing strategy white papers. Here we focus on aspects of `static' lens discovery that will be affected by the observing strategy. In summary, we advocate (1) ensuring comparable (sub-arcsecond) seeing in the g-band as in r and i to facilitate discovery of gravitational lenses, and (2) initially surveying the entire observable extragalactic sky as rapidly as possible to enable early science spanning a broad range of static and transient interests.
1902.05165
Non-thermal emission from the interaction of magnetized exoplanets with the wind of their host star
Wang, Loeb
We study the non-thermal emission from the interaction between magnetized Jupiter-like exoplanets and the wind from their host star. The supersonic motion of planets through the wind forms a bow shock that accelerates electrons which produces non-thermal radiation across a broad wavelength range. We discuss three wind mass loss rates: $\dot{M}_{\rm w}\sim10^{-14}$, $10^{-9}$, $10^{-6}\,M_{\odot}\,\rm yr^{-1}$ corresponding to solar-type, T Tauri and massive O/B type stars, respectively. We find that the expected radio synchrotron emission from a Jupiter-like planet is detectable by the Jansky Very Large Array and the Square Kilometer Array at ~1-10 GHz out to a distance ~ 100 pc, whereas the infrared emission is detectable by the James Webb Space Telescope out to a similar distance. Inverse Compton scattering of the stellar radiation results in X-ray emission detectable by Chandra X-ray Observatory out to ~ 150 pc. Finally, we apply our model to the upper limit constraints on V380 Tau, the first star-hot Jupiter system observed in radio wavelength. Our bow shock model provides constraints on the magnetic field, the interplanetary medium and the non-thermal emission efficiency in V380 Tau.
1902.05294
Solar cycle activity: an early prediction for cycle #25
Sello
Solar activity forecasting is an important topic for numerous scientific and technological areas, such as space mission operations, electric power transmission lines, power transformation stations and earth geophysical and climatic impact. Nevertheless, the well-known difficulty is how to accurately predict, on the basis of various recorded solar activity indices, the complete evolution of future solar cycles, due to highly complex dynamical and stochastic processes involved, mainly related to interaction of different components of internal magnetic fields. There are two main distinct classes of solar cycle prediction methods: the precursor-like ones and the mathematical-numerical ones. The main characteristic of precursor techniques, both purely solar and geomagnetic, is their physical basis. Conversely, the non-precursor methods use different mathematical and/or numerical properties of the known temporal evolution of solar activity indices to extract useful information for predicting future activity. For current solar cycle #24 we obtained fairly good statistical performances from both precursor and purely numerical methods, such as the so-called solar precursor and nonlinear ones. To further check the performances of these prediction techniques, we compared the early predictions for the next solar cycle #25. Preliminary results support some coherence of the prediction methods considered and confirm the current trend of a relatively low solar activity.
1902.05440
Kernel formalism applied to Fourier based wave front sensing in presence of residual phases
Fauvarque, et al
In this paper, we show how it is possible to understand Fourier-based Wave Front Sensors (WFS) as Linear Integral Operators. Indeed optical laws allow to identify quantities called "Kernels" which are the continuous version of well-known interaction matrices used in Adaptive Optics. The main purpose of this article is to understand the dependence of these Kernels regarding to the optical parameters of the WFSs: the entrance pupil geometry, the filtering mask and the tip/tilt modulation. Kernels are then linked to the classical performance criterion which is the sensitivity. This approach also allows to take into account time-varying phase residuals. As it turns out, they have physically the same status as the modulation: they are both "pupil plane light shaping". The only significant difference is their time-dependent behaviour. Modulation is highly regular and set by the experimenter while phase residuals are erratic and inevitable. Fortunately, the latter may be described thanks to their statistics, namely their "structure function" which naturally arises in the Kernel framework. These results are essential to tackle the Pyramid Wave Front Sensor's optical gain problem which is at the present time the major factor which constrains the Pyramid WFS performance and operability. A second part focuses on the special case of "convolutional Kernels". They drastically simplify the mathematical formulation since the WFS's input and output are linked by a convolution product. Such a formalism firstly provides a fast diagnostic tool to study Fourier-based WFSs but also suggests natural algorithms to reconstruct the phase, e.g. Wigner deconvolution. We pay special attention to the assumptions required to be in such a regime (which does not apply in general). Fortunately, we show that a wise choice of the modulation allows to improve its validity.
No comments:
Post a Comment