Showing 151–200 of 596 results for author: Lewis, G

The GALAH survey: A census of lithium-rich giant stars

Authors: Sarah Martell, Jeffrey Simpson, Adithya Balasubramaniam, Sven Buder, Sanjib Sharma, Marc Hon, Dennis Stello, Yuan-Sen Ting, Martin Asplund, Joss Bland-Hawthorn, Gayandhi De Silva, Ken Freeman, Michael Hayden, Janez Kos, Geraint Lewis, Karin Lind, Daniel Zucker, Tomaz Zwitter, Simon Campbell, Klemen Cotar, Jonathan Horner, Benjamin Montet, Rob Wittenmyer

Abstract: We investigate the properties of 1262 red giant stars with high photospheric abundances of lithium observed by the GALAH and \Ktwo-HERMES surveys, and discuss them in the context of proposed mechanisms for lithium enrichment and re-depletion in giant stars. We confirm that Li-rich giants are rare, making up only 1.2 per cent of our giant star sample. We use stellar parameters from the third public… ▽ More We investigate the properties of 1262 red giant stars with high photospheric abundances of lithium observed by the GALAH and \Ktwo-HERMES surveys, and discuss them in the context of proposed mechanisms for lithium enrichment and re-depletion in giant stars. We confirm that Li-rich giants are rare, making up only 1.2 per cent of our giant star sample. We use stellar parameters from the third public data release from the GALAH survey and a Bayesian isochrone analysis to divide the sample into first-ascent red giant branch and red clump stars, and confirm these classifications using asteroseismic data from \Ktwo. We find that red clump stars are 2.5 times as likely to be lithium-rich as red giant branch stars, in agreement with other recent work. The probability for a star to be lithium-rich is affected by a number of factors, though the causality in those correlations is not entirely clear. We show for the first time that primary and secondary red clump stars have distinctly different lithium enrichment patterns. The data set discussed here is large and heterogeneous in terms of evolutionary phase, metallicity, rotation rate and mass. We expect that if the various mechanisms that have been proposed for lithium enrichment in evolved stars are in fact active, they should all contribute to this sample of lithium-rich giants at some level. △ Less

Submitted 10 May, 2021; v1 submitted 3 June, 2020; originally announced June 2020.

Comments: 17 pages, 12 figures, MNRAS accepted

OzDES multi-object fibre spectroscopy for the Dark Energy Survey: Results and second data release

Authors: C. Lidman, B. E. Tucker, T. M. Davis, S. A. Uddin, J. Asorey, K. Bolejko, D. Brout, J. Calcino, D. Carollo, A. Carr, M. Childress, J. K. Hoormann, R. J. Foley, L. Galbany, K. Glazebrook, S. R. Hinton, R. Kessler, A. G. Kim, A. King, A. Kremin, K. Kuehn, D. Lagattuta, G. F. Lewis, E. Macaulay, U. Malik , et al. (79 additional authors not shown)

Abstract: We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the… ▽ More We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30,000 sources, some as faint as $r_{\mathrm AB}=24$ mag, and 375,000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to $z\sim1.2$ and the masses of hundreds of black holes out to $z\sim4$. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs. △ Less

Submitted 31 May, 2020; originally announced June 2020.

Comments: 20 pages, 16 figures. Accepted for publication by MNRAS. Data release available at https://datacentral.org.au

Report number: Report-no: FERMILAB-PUB-20-029-AE; DES 2019-0517

K2-HERMES II. Planet-candidate properties from K2 Campaigns 1-13

Authors: Robert A. Wittenmyer, Jake T. Clark, Sanjib Sharma, Dennis Stello, Jonathan Horner, Stephen R. Kane, Catherine P. Stevens, Duncan J. Wright, Lorenzo Spina, Klemen Cotar, Martin Asplund, Joss Bland-Hawthorn, Sven Buder, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken Freeman, Janez Kos, Geraint Lewis, Jane Lin, Karin Lind, Sarah L. Martell, Jeffrey D. Simpson, Daniel B. Zucker, Tomaz Zwitter

Abstract: Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present complete results for planet-candidate hosts from the K2-HERMES survey, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28… ▽ More Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present complete results for planet-candidate hosts from the K2-HERMES survey, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra for more than 30,000 K2 stars. We present complete host-star parameters and planet-candidate radii for 224 K2 candidate planets from C1-C13. Our results cast severe doubt on 30 K2 candidates, as we derive unphysically large radii, larger than 2 $R_{Jup}$. This work highlights the importance of obtaining accurate, precise, and self-consistent stellar parameters for ongoing large planet search programs - something that will only become more important in the coming years, as TESS begins to deliver its own harvest of exoplanets. △ Less

Submitted 21 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS. Full versions of online-only Tables 2 and 3 available on request. Same for Table 4, in exchange for the trick to make MNRAS style break long tables across pages

Studying Type II supernovae as cosmological standard candles using the Dark Energy Survey

Authors: T. de Jaeger, L. Galbany, S. González-Gaitán, R. Kessler, A. V. Filippenko, F. Förster, M. Hamuy, P. J. Brown, T. M. Davis, C. P. Gutiérrez, C. Inserra, G F. Lewis, A. Möller, D. Scolnic, M. Smith, D. Brout, D. Carollo, R. J. Foley, K. Glazebrook, S. R. Hinton, E. Macaulay, B. Nichol, M. Sako, N. E. Sommer, B. E. Tucker , et al. (63 additional authors not shown)

Abstract: Despite vast improvements in the measurement of the cosmological parameters, the nature of dark energy and an accurate value of the Hubble constant (H$_0$) in the Hubble-Lemaître law remain unknown. To break the current impasse, it is necessary to develop as many independent techniques as possible, such as the use of Type II supernovae (SNe II). The goal of this paper is to demonstrate the utility… ▽ More Despite vast improvements in the measurement of the cosmological parameters, the nature of dark energy and an accurate value of the Hubble constant (H$_0$) in the Hubble-Lemaître law remain unknown. To break the current impasse, it is necessary to develop as many independent techniques as possible, such as the use of Type II supernovae (SNe II). The goal of this paper is to demonstrate the utility of SNe II for deriving accurate extragalactic distances, which will be an asset for the next generation of telescopes where more-distant SNe II will be discovered. More specifically, we present a sample from the Dark Energy Survey Supernova Program (DES-SN) consisting of 15 SNe II with photometric and spectroscopic information spanning a redshift range up to 0.35. Combining our DES SNe with publicly available samples, and using the standard candle method (SCM), we construct the largest available Hubble diagram with SNe II in the Hubble flow (70 SNe II) and find an observed dispersion of 0.27 mag. We demonstrate that adding a colour term to the SN II standardisation does not reduce the scatter in the Hubble diagram. Although SNe II are viable as distance indicators, this work points out important issues for improving their utility as independent extragalactic beacons: find new correlations, define a more standard subclass of SNe II, construct new SN II templates, and dedicate more observing time to high-redshift SNe II. Finally, for the first time, we perform simulations to estimate the redshift-dependent distance-modulus bias due to selection effects. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: 39 pages, 22 figures, 10 tables, Accepted for publication in MNRAS

The Host Galaxies of Rapidly Evolving Transients in the Dark Energy Survey

Authors: P. Wiseman, M. Pursiainen, M. Childress, E. Swann, M. Smith, L. Galbany, C. Lidman, T. M. Davis, C. P. Gutiérrez, A. Möller, B. P. Thomas, C. Frohmaier, R. J. Foley, S. R. Hinton, L. Kelsey, R. Kessler, G. F. Lewis, M. Sako, D. Scolnic, M. Sullivan, M. Vincenzi, T. M. C. Abbott, M. Aguena, S. Allam, J. Annis , et al. (50 additional authors not shown)

Abstract: Rapidly evolving transients (RETs), also termed fast blue optical transients (FBOTs), are a distinct class of astrophysical event. They are characterised by lightcurves that decline much faster than common classes supernovae (SNe), span vast ranges in peak luminosity and can be seen to redshifts greater than 1. Their evolution on fast timescales has hindered high quality follow-up observations, an… ▽ More Rapidly evolving transients (RETs), also termed fast blue optical transients (FBOTs), are a distinct class of astrophysical event. They are characterised by lightcurves that decline much faster than common classes supernovae (SNe), span vast ranges in peak luminosity and can be seen to redshifts greater than 1. Their evolution on fast timescales has hindered high quality follow-up observations, and thus their origin and explosion/emission mechanism remains unexplained. In this paper we define the largest sample of RETs to date, comprising 106 objects from the Dark Energy Survey, and perform the most comprehensive analysis of RET host galaxies. Using deep-stacked photometry and emission-lines from OzDES spectroscopy, we derive stellar masses and star-formation rates (SFRs) for 49 host galaxies, and metallicities for 37. We find that RETs explode exclusively in star-forming galaxies and are thus likely associated with massive stars. Comparing RET hosts to samples of host galaxies of other explosive transients as well as field galaxies, we find that RETs prefer galaxies with high specific SFRs, indicating a link to young stellar populations, similar to stripped-envelope SNe. RET hosts appear to show a lack of chemical enrichment, their metallicities akin to long duration gamma-ray bursts and superluminous SN host galaxies. There are no clear relationships between properties of the host galaxies and the peak magnitudes or decline rates of the transients themselves. △ Less

Submitted 18 May, 2020; originally announced May 2020.

Comments: 16 pages, 15 figures, submitted to MNRAS

Report number: FERMILAB-PUB-20-206-AE; DES-2020-0534

The SAMI Galaxy Survey: Gas velocity dispersions in low-$z$ star-forming galaxies and the drivers of turbulence

Authors: Mathew R. Varidel, Scott M. Croom, Geraint F. Lewis, Deanne B. Fisher, Karl Glazebrook, Barbara Catinella, Luca Cortese, Mark R. Krumholz, Joss Bland-Hawthorn, Julia J. Bryant, Brent Groves, Sarah Brough, Christoph Federrath, Jon S. Lawrence, Nuria P. Lorente, Matt S. Owers, Samuel N. Richards, Ángel R. López-Sánchez, Sarah M. Sweet, Jesse van de Sande, Sam P. Vaughan

Abstract: We infer the intrinsic ionised gas kinematics for 383 star-forming galaxies across a range of integrated star-formation rates (SFR $\in [10^{-3}, 10^2]$ M$_\odot$ yr$^{-1}$) at $z \lesssim 0.1$ using a consistent 3D forward-modelling technique. The total sample is a combination of galaxies from the SAMI Galaxy Survey and DYNAMO survey. For typical low-$z$ galaxies taken from the SAMI Galaxy Survey… ▽ More We infer the intrinsic ionised gas kinematics for 383 star-forming galaxies across a range of integrated star-formation rates (SFR $\in [10^{-3}, 10^2]$ M$_\odot$ yr$^{-1}$) at $z \lesssim 0.1$ using a consistent 3D forward-modelling technique. The total sample is a combination of galaxies from the SAMI Galaxy Survey and DYNAMO survey. For typical low-$z$ galaxies taken from the SAMI Galaxy Survey, we find the vertical velocity dispersion ($σ_{v, z}$) to be positively correlated with measures of star-formation rate, stellar mass, HI gas mass, and rotational velocity. The greatest correlation is with star-formation rate surface density ($Σ_\text{SFR}$). Using the total sample, we find $σ_{v, z}$ increases slowly as a function of integrated star-formation rate in the range SFR $\in$ [$10^{-3}$, 1] M$_\odot$ yr$^{-1}$ from $17\pm3$ km s$^{-1}$ to $24\pm5$ km s$^{-1}$ followed by a steeper increase up to $σ_{v, z}$ $\sim 80$ km s$^{-1}$ for SFR $\gtrsim 1$ M$_\odot$ yr$^{-1}$. This is consistent with recent theoretical models that suggest a $σ_{v, z}$ floor driven by star-formation feedback processes with an upturn in $σ_{v, z}$ at higher SFR driven by gravitational transport of gas through the disc. △ Less

Submitted 17 May, 2020; v1 submitted 11 May, 2020; originally announced May 2020.

Comments: 27 pages, 10 figures, accepted for MNRAS

The GALAH survey: Multiple stars and our Galaxy. I. A comprehensive method for deriving properties of FGK binary stars

Authors: G. Traven, S. Feltzing, T. Merle, M. Van der Swaelmen, K. Čotar, R. Church, T. Zwitter, Y. -S. Ting, C. Sahlholdt, M. Asplund, J. Bland-Hawthorn, G. De Silva, K. Freeman, S. Martell, S. Sharma, D. Zucker, S. Buder, A. Casey, V. D'Orazi, J. Kos, G. Lewis, J. Lin, K. Lind, J. Simpson, D. Stello , et al. (2 additional authors not shown)

Abstract: Binary stellar systems form a large fraction of the Galaxy's stars. They are useful as laboratories for studying the physical processes taking place within stars, and must be correctly taken into account when observations of stars are used to study the structure and evolution of the Galaxy. We present a sample of 12760 well-characterised double-lined spectroscopic binaries that are appropriate for… ▽ More Binary stellar systems form a large fraction of the Galaxy's stars. They are useful as laboratories for studying the physical processes taking place within stars, and must be correctly taken into account when observations of stars are used to study the structure and evolution of the Galaxy. We present a sample of 12760 well-characterised double-lined spectroscopic binaries that are appropriate for statistical studies of the binary populations. They were detected as SB2s using a t-distributed stochastic neighbour embedding (t-SNE) classification and a cross-correlation analysis of GALAH spectra. This sample consists mostly of dwarfs, with a significant fraction of evolved stars and several dozen members of the giant branch. To compute parameters of the primary and secondary star ($T_{\rm eff[1,2]}$, $\log g_{[1,2]}$, [Fe/H], $V_{r[1,2]}$, $v_{\rm mic[1,2]}$, $v_{\rm broad[1,2]}$, $R_{[1,2]}$, and $E(B-V)$), we used a Bayesian approach that includes a parallax prior from Gaia DR2, spectra from GALAH, and apparent magnitudes from APASS, Gaia DR2, 2MASS, and WISE. The derived stellar properties and their distributions show trends that are expected for a population of close binaries (a $<$ 10 AU) with mass ratios $0.5 \leq q \leq 1$. The derived metallicity of these binary stars is statistically lower than that of single dwarf stars from the same magnitude-limited sample. △ Less

Submitted 30 April, 2020; originally announced May 2020.

Comments: Accepted for publication in A&A

MSC Class: 85-XX

Journal ref: A&A 638, A145 (2020)

The first Hubble diagram and cosmological constraints using superluminous supernova

Authors: C. Inserra, M. Sullivan, C. R. Angus, E. Macaulay, R. C. Nichol, M. Smith, C. Frohmaier, C. P. Gutiérrez, M. Vicenzi, A. Möller, D. Brout, P. J. Brown, T. M. Davis, C. B. D'Andrea, L. Galbany, R. Kessler, A. G. Kim, Y. -C. Pan, M. Pursiainen, D. Scolnic, B. P. Thomas, P. Wiseman, T. M. C. Abbott, J. Annis, S. Avila , et al. (66 additional authors not shown)

Abstract: We present the first Hubble diagram of superluminous supernovae (SLSNe) out to a redshift of two, together with constraints on the matter density, $Ω_{\rm M}$, and the dark energy equation-of-state parameter, $w(\equiv p/ρ)$. We build a sample of 20 cosmologically useful SLSNe~I based on light curve and spectroscopy quality cuts. We confirm the robustness of the peak decline SLSN~I standardization… ▽ More We present the first Hubble diagram of superluminous supernovae (SLSNe) out to a redshift of two, together with constraints on the matter density, $Ω_{\rm M}$, and the dark energy equation-of-state parameter, $w(\equiv p/ρ)$. We build a sample of 20 cosmologically useful SLSNe~I based on light curve and spectroscopy quality cuts. We confirm the robustness of the peak decline SLSN~I standardization relation with a larger dataset and improved fitting techniques than previous works. We then solve the SLSN model based on the above standardisation via minimisation of the $χ^2$ computed from a covariance matrix which includes statistical and systematic uncertainties. For a spatially flat $Λ$CDM cosmological model, we find $Ω_{\rm M}=0.38^{+0.24}_{-0.19}$, with a rms of 0.27 mag for the residuals of the distance moduli. For a $w_0w_a$CDM cosmological model, the addition of SLSNe~I to a `baseline' measurement consisting of Planck temperature together with type Ia supernovae, results in a small improvement in the constraints of $w_0$ and $w_a$ of 4\%. We present simulations of future surveys with 868 and 492 SLSNe I (depending on the configuration used) and show that such a sample can deliver cosmological constraints in a flat $Λ$CDM model with the same precision (considering only statistical uncertainties) as current surveys that use type Ia supernovae, while providing a factor 2-3 improvement in the precision of the constraints on the time variation of dark energy, $w_0$ and $w_a$. This paper represents the proof-of-concept for superluminous supernova cosmology, and demonstrates they can provide an independent test of cosmology in the high-redshift ($z>1$) universe. △ Less

Submitted 8 April, 2021; v1 submitted 25 April, 2020; originally announced April 2020.

Comments: 16 pages, 9 figures. MNRAS accepted

Fundamental relations for the velocity dispersion of stars in the Milky Way

Authors: Sanjib Sharma, Michael R. Hayden, Joss Bland-Hawthorn, Dennis Stello, Sven Buder, Joel C. Zinn, Thomas Kallinger, Martin Asplund, Gayandhi M. De Silva, Valentina Dorazi, Ken Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah Martell, Jeffrey D. Simpson, Rob A. Wittenmyer, Daniel B. Zucker, Tomaz Zwitter, Boquan Chen, Klemen Cotar, James Esdaile, Marc Hon, Jonathan Horner , et al. (12 additional authors not shown)

Abstract: We explore the fundamental relations governing the radial and vertical velocity dispersions of stars in the Milky Way, from combined studies of complementary surveys including GALAH, LAMOST, APOGEE, the NASA $Kepler$ and K2 missions, and $Gaia$ DR2. We find that different stellar samples, even though they target different tracer populations and employ a variety of age estimation techniques, follow… ▽ More We explore the fundamental relations governing the radial and vertical velocity dispersions of stars in the Milky Way, from combined studies of complementary surveys including GALAH, LAMOST, APOGEE, the NASA $Kepler$ and K2 missions, and $Gaia$ DR2. We find that different stellar samples, even though they target different tracer populations and employ a variety of age estimation techniques, follow the same set of fundamental relations. We provide the clearest evidence to date that, in addition to the well-known dependence on stellar age, the velocity dispersions of stars depend on orbital angular momentum $L_z$, metallicity and height above the plane $|z|$, and are well described by a multiplicatively separable functional form. The dispersions have a power-law dependence on age with exponents of 0.441$\pm 0.007$ and 0.251$\pm 0.006$ for $σ_z$ and $σ_R$ respectively, and the power law is valid even for the oldest stars. For the solar neighborhood stars, the apparent break in the power law for older stars, as seen in previous studies, is due to the anti-correlation of $L_z$ with age. The dispersions decrease with increasing $L_z$ until we reach the Sun's orbital angular momentum, after which $σ_z$ increases (implying flaring in the outer disc) while $σ_R$ flattens. The dispersions increase with decreasing metallicity, suggesting that the dispersions increase with birth radius. The dispersions also increase linearly with $|z|$. The same set of relations that work in the solar neighborhood also work for stars between $3<R/{\rm kpc}<20$. Finally, the high-[$α$/Fe] stars follow the same relations as the low-[$α$/Fe] stars. △ Less

Submitted 14 April, 2020; originally announced April 2020.

Comments: 19 pages, 15 figures, submitted to MNRAS

Predicting the vulnerability of spacecraft components: modelling debris impact effects through vulnerable-zones

Authors: Mirko Trisolini, Hugh G. Lewis, Camilla Colombo

Abstract: The space environment around the Earth is populated by more than 130 million objects of 1 mm in size and larger, and future predictions shows that this amount is destined to increase, even if mitigation measures are implemented at a far better rate than today. These objects can hit and damage a spacecraft or its components. It is thus necessary to assess the risk level for a satellite during its m… ▽ More The space environment around the Earth is populated by more than 130 million objects of 1 mm in size and larger, and future predictions shows that this amount is destined to increase, even if mitigation measures are implemented at a far better rate than today. These objects can hit and damage a spacecraft or its components. It is thus necessary to assess the risk level for a satellite during its mission lifetime. Few software packages perform this analysis, and most of them employ time-consuming ray-tracing methodology, where particles are randomly sampled from relevant distributions. In addition, they tend not to consider the risk associated with the secondary debris clouds. The paper presents the development of a vulnerability assessment model, which relies on a fully statistical procedure: the debris fluxes are directly used combining them with the concept of the vulnerable zone, avoiding the random sampling the debris fluxes. A novel methodology is presented to predict damage to internal components. It models the interaction between the components and the secondary debris cloud through basic geometric operations, considering mutual shielding and shadowing between internal components. The methodologies are tested against state-of-the-art software for relevant test cases, comparing results on external structures and internal components. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: Article accepted for pubblication in Advances in Space Research

Journal ref: Adv.Space.Res. 65. 11 (2020). 2692-2710

A unified framework for 21cm tomography sample generation and parameter inference with Progressively Growing GANs

Authors: Florian List, Geraint F. Lewis

Abstract: Creating a database of 21cm brightness temperature signals from the Epoch of Reionisation (EoR) for an array of reionisation histories is a complex and computationally expensive task, given the range of astrophysical processes involved and the possibly high-dimensional parameter space that is to be probed. We utilise a specific type of neural network, a Progressively Growing Generative Adversarial… ▽ More Creating a database of 21cm brightness temperature signals from the Epoch of Reionisation (EoR) for an array of reionisation histories is a complex and computationally expensive task, given the range of astrophysical processes involved and the possibly high-dimensional parameter space that is to be probed. We utilise a specific type of neural network, a Progressively Growing Generative Adversarial Network (PGGAN), to produce realistic tomography images of the 21cm brightness temperature during the EoR, covering a continuous three-dimensional parameter space that models varying X-ray emissivity, Lyman band emissivity, and ratio between hard and soft X-rays. The GPU-trained network generates new samples at a resolution of $\sim 3'$ in a second (on a laptop CPU), and the resulting global 21cm signal, power spectrum, and pixel distribution function agree well with those of the training data, taken from the 21SSD catalogue \citep{Semelin2017}. Finally, we showcase how a trained PGGAN can be leveraged for the converse task of inferring parameters from 21cm tomography samples via Approximate Bayesian Computation. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Comments: 15 pages, 8+1 figures, accepted by MNRAS

Double/Debiased Machine Learning for Dynamic Treatment Effects via g-Estimation

Authors: Greg Lewis, Vasilis Syrgkanis

Abstract: We consider the estimation of treatment effects in settings when multiple treatments are assigned over time and treatments can have a causal effect on future outcomes or the state of the treated unit. We propose an extension of the double/debiased machine learning framework to estimate the dynamic effects of treatments, which can be viewed as a Neyman orthogonal (locally robust) cross-fitted versi… ▽ More We consider the estimation of treatment effects in settings when multiple treatments are assigned over time and treatments can have a causal effect on future outcomes or the state of the treated unit. We propose an extension of the double/debiased machine learning framework to estimate the dynamic effects of treatments, which can be viewed as a Neyman orthogonal (locally robust) cross-fitted version of $g$-estimation in the dynamic treatment regime. Our method applies to a general class of non-linear dynamic treatment models known as Structural Nested Mean Models and allows the use of machine learning methods to control for potentially high dimensional state variables, subject to a mean square error guarantee, while still allowing parametric estimation and construction of confidence intervals for the structural parameters of interest. These structural parameters can be used for off-policy evaluation of any target dynamic policy at parametric rates, subject to semi-parametric restrictions on the data generating process. Our work is based on a recursive peeling process, typical in $g$-estimation, and formulates a strongly convex objective at each stage, which allows us to extend the $g$-estimation framework in multiple directions: i) to provide finite sample guarantees, ii) to estimate non-linear effect heterogeneity with respect to fixed unit characteristics, within arbitrary function spaces, enabling a dynamic analogue of the RLearner algorithm for heterogeneous effects, iii) to allow for high-dimensional sparse parameterizations of the target structural functions, enabling automated model selection via a recursive lasso algorithm. We also provide guarantees for data stemming from a single treated unit over a long horizon and under stationarity conditions. △ Less

Submitted 16 June, 2021; v1 submitted 17 February, 2020; originally announced February 2020.

Supernova Siblings: Assessing the Consistency of Properties of Type Ia Supernovae that Share the Same Parent Galaxies

Authors: D. Scolnic, M. Smith, A. Massiah, P. Wiseman, D. Brout, R. Kessler, T. M. Davis, R. J. Foley, L. Galbany, S. R. Hinton, R. Hounsell, L. Kelsey, C. Lidman, R. Morgan, R. C. Nichol, A. Möller, B. Popovic, M. Sako, M. Sullivan, B. P. Thomas, T. M. C. Abbott, M. Aguena, S. Allam, J. Annis, S. Avila , et al. (60 additional authors not shown)

Abstract: While many studies have shown a correlation between properties of the light curves of Type Ia SN (SNe Ia) and properties of their host galaxies, it remains unclear what is driving these correlations. We introduce a new direct method to study these correlations by analyzing `parent' galaxies that host multiple SNe Ia 'siblings'. Here, we search the Dark Energy Survey SN sample, one of the largest s… ▽ More While many studies have shown a correlation between properties of the light curves of Type Ia SN (SNe Ia) and properties of their host galaxies, it remains unclear what is driving these correlations. We introduce a new direct method to study these correlations by analyzing `parent' galaxies that host multiple SNe Ia 'siblings'. Here, we search the Dark Energy Survey SN sample, one of the largest samples of discovered SNe, and find 8 galaxies that hosted two likely Type Ia SNe. Comparing the light-curve properties of these SNe and recovered distances from the light curves, we find no better agreement between properties of SNe in the same galaxy as any random pair of galaxies, with the exception of the SN light-curve stretch. We show at $2.8σ$ significance that at least 1/2 of the intrinsic scatter of SNe Ia distance modulus residuals is not from common host properties. We also discuss the robustness with which we could make this evaluation with LSST, which will find $100\times$ more pairs of galaxies, and pave a new line of study on the consistency of Type Ia supernovae in the same parent galaxies. Finally, we argue that it is unlikely some of these SNe are actually single, lensed SN with multiple images. △ Less

Submitted 5 February, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: Submitted to ApJL. Comments welcome. Fermilab id included

DES16C3cje: A low-luminosity, long-lived supernova

Authors: C. P. Gutiérrez, M. Sullivan, L. Martinez, M. C. Bersten, C. Inserra, M. Smith, J. P. Anderson, Y. -C. Pan, A. Pastorello, L. Galbany, P. Nugent, C. R. Angus, C. Barbarino, T. -W. Chen, T. M. Davis, M. Della Valle, R. J. Foley, M. Fraser, C. Frohmaier, S. González-Gaitán, G. F. Lewis, M. Gromadzki, E. Kankare, R. Kokotanekova, J. Kollmeier , et al. (67 additional authors not shown)

Abstract: We present DES16C3cje, a low-luminosity, long-lived type II supernova (SN II) at redshift 0.0618, detected by the Dark Energy Survey (DES). DES16C3cje is a unique SN. The spectra are characterized by extremely narrow photospheric lines corresponding to very low expansion velocities of $\lesssim1500$ km s$^{-1}$, and the light curve shows an initial peak that fades after 50 days before slowly rebri… ▽ More We present DES16C3cje, a low-luminosity, long-lived type II supernova (SN II) at redshift 0.0618, detected by the Dark Energy Survey (DES). DES16C3cje is a unique SN. The spectra are characterized by extremely narrow photospheric lines corresponding to very low expansion velocities of $\lesssim1500$ km s$^{-1}$, and the light curve shows an initial peak that fades after 50 days before slowly rebrightening over a further 100 days to reach an absolute brightness of M$_r\sim -15.5$ mag. The decline rate of the late-time light curve is then slower than that expected from the powering by radioactive decay of $^{56}$Co but is comparable to that expected from accretion power. Comparing the bolometric light curve with hydrodynamical models, we find that DES16C3cje can be explained by either i) a low explosion energy (0.11 foe) and relatively large $^{56}$Ni production of 0.075 M$_{\odot}$ from a $\sim15$ M$_{\odot}$ red supergiant progenitor typical of other SNe II, or ii) a relatively compact $\sim40$ M$_{\odot}$ star, explosion energy of 1 foe, and 0.08 M$_{\odot}$ of $^{56}$Ni. Both scenarios require additional energy input to explain the late-time light curve, which is consistent with fallback accretion at a rate of $\sim0.5\times{10^{-8}}$ M$_{\odot}$ s$^{-1}$. △ Less

Submitted 29 May, 2020; v1 submitted 30 January, 2020; originally announced January 2020.

Comments: Accepted in MNRAS. 17 pages, 6 figures

First Cosmology Results using Type Ia Supernovae from the Dark Energy Survey: The Effect of Host Galaxy Properties on Supernova Luminosity

Authors: M. Smith, M. Sullivan, P. Wiseman, R. Kessler, D. Scolnic, D. Brout, C. B. D'Andrea, T. M. Davis, R. J. Foley, C. Frohmaier, L. Galbany, R. R. Gupta, C. P. Gutiérrez, S. R. Hinton, L. Kelsey, C. Lidman, E. Macaulay, A. Möller, R. C. Nichol, P. Nugent, A. Palmese, M. Pursiainen, M. Sako, R. C. Thomas, B. E. Tucker , et al. (66 additional authors not shown)

Abstract: We present improved photometric measurements for the host galaxies of 206 spectroscopically confirmed type Ia supernovae discovered by the Dark Energy Survey Supernova Program (DES-SN) and used in the first DES-SN cosmological analysis. Fitting spectral energy distributions to the $griz$ photometric measurements of the DES-SN host galaxies, we derive stellar masses and star-formation rates. For th… ▽ More We present improved photometric measurements for the host galaxies of 206 spectroscopically confirmed type Ia supernovae discovered by the Dark Energy Survey Supernova Program (DES-SN) and used in the first DES-SN cosmological analysis. Fitting spectral energy distributions to the $griz$ photometric measurements of the DES-SN host galaxies, we derive stellar masses and star-formation rates. For the DES-SN sample, when considering a 5D ($z$, $x_1$, $c$, $α$, $β$) bias correction, we find evidence of a Hubble residual `mass step', where SNe Ia in high mass galaxies ($>10^{10} \textrm{M}_{\odot}$) are intrinsically more luminous (after correction) than their low mass counterparts by $γ=0.040\pm0.019$mag. This value is larger by $0.031$mag than the value found in the first DES-SN cosmological analysis. This difference is due to a combination of updated photometric measurements and improved star formation histories and is not from host-galaxy misidentification. When using a 1D (redshift-only) bias correction the inferred mass step is larger, with $γ=0.066\pm0.020$mag. The 1D-5D $γ$ difference for DES-SN is $0.026\pm0.009$mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the $x_1$ component of the 5D distance-bias correction. To better understand this effect, we include an intrinsic correlation between light-curve width and stellar mass in simulated SN Ia samples. We show that a 5D fit recovers $γ$ with $-9$mmag bias compared to a $+2$mmag bias for a 1D fit. This difference can explain part of the discrepancy seen in the data. Improvements in modeling correlations between galaxy properties and SN is necessary to determine the implications for $γ$ and ensure unbiased precision estimates of the dark energy equation-of-state as we enter the era of LSST. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: 27 pages, 13 figures; Submitted to MNRAS

Report number: FERMILAB-PUB-20-003-AE; DES-2018-0402

Abundances in the Milky Way across five nucleosynthetic channels from 4 million LAMOST stars

Authors: Adam Wheeler, Melissa Ness, Sven Buder, Joss Bland-Hawthorn, Gayandhi De Silva, Michael Hayden, Janez Kos, Geraint F. Lewis, Sarah Martell, Sanjib Sharma, Jeffrey D. Simpson, D. B. Zucker, Thomaž Zwitter

Abstract: Large stellar surveys are revealing the chemodynamical structure of the Galaxy across a vast spatial extent. However, the many millions of low-resolution spectra observed to date are yet to be fully exploited. We employ The Cannon, a data-driven approach to estimating abundances, to obtain detailed abundances from low-resolution (R = 1800) LAMOST spectra, using the GALAH survey as our reference. W… ▽ More Large stellar surveys are revealing the chemodynamical structure of the Galaxy across a vast spatial extent. However, the many millions of low-resolution spectra observed to date are yet to be fully exploited. We employ The Cannon, a data-driven approach to estimating abundances, to obtain detailed abundances from low-resolution (R = 1800) LAMOST spectra, using the GALAH survey as our reference. We deliver five (for dwarfs) or six (for giants) estimated abundances representing five different nucleosynthetic channels, for 3.9 million stars, to a precision of 0.05 - 0.23 dex. Using wide binary pairs, we demonstrate that our abundance estimates provide chemical discriminating power beyond metallicity alone. We show the coverage of our catalogue with radial, azimuthal and dynamical abundance maps, and examine the neutron capture abundances across the disk and halo, which indicate different origins for the in-situ and accreted halo populations. LAMOST has near-complete Gaia coverage and provides an unprecedented perspective on chemistry across the Milky Way. △ Less

Submitted 11 June, 2020; v1 submitted 22 January, 2020; originally announced January 2020.

Comments: 17 pages + refs and appendices, to appear in ApJ. Data: https://doi.org/10.7910/DVN/5VWKMC

Supernova Host Galaxies in the Dark Energy Survey: I. Deep Coadds, Photometry, and Stellar Masses

Authors: P. Wiseman, M. Smith, M. Childress, L. Kelsey, A. Möller, R. R. Gupta, E. Swann, C. R. Angus, D. Brout, T. M. Davis, R. J. Foley, C. Frohmaier, L. Galbany, C. P. Gutiérrez, R. Kessler, G. F. Lewis, C. Lidman, E. Macaulay, R. C. Nichol, M. Pursiainen, M. Sako, D. Scolnic, N. E. Sommer, M. Sullivan, B. E. Tucker , et al. (53 additional authors not shown)

Abstract: The five-year Dark Energy Survey supernova programme (DES-SN) is one of the largest and deepest transient surveys to date in terms of volume and number of supernovae. Identifying and characterising the host galaxies of transients plays a key role in their classification, the study of their formation mechanisms, and the cosmological analyses. To derive accurate host galaxy properties, we create dep… ▽ More The five-year Dark Energy Survey supernova programme (DES-SN) is one of the largest and deepest transient surveys to date in terms of volume and number of supernovae. Identifying and characterising the host galaxies of transients plays a key role in their classification, the study of their formation mechanisms, and the cosmological analyses. To derive accurate host galaxy properties, we create depth-optimised coadds using single-epoch DES-SN images that are selected based on sky and atmospheric conditions. For each of the five DES-SN seasons, a separate coadd is made from the other 4 seasons such that each SN has a corresponding deep coadd with no contaminating SN emission. The coadds reach limiting magnitudes of order $\sim 27$ in $g$-band, and have a much smaller magnitude uncertainty than the previous DES-SN host templates, particularly for faint objects. We present the resulting multi-band photometry of host galaxies for samples of spectroscopically confirmed type Ia (SNe Ia), core-collapse (CCSNe), and superluminous (SLSNe) as well as rapidly evolving transients (RETs) discovered by DES-SN. We derive host galaxy stellar masses and probabilistically compare stellar-mass distributions to samples from other surveys. We find that the DES spectroscopically confirmed sample of SNe Ia selects preferentially fewer high mass hosts at high redshift compared to other surveys, while at low redshift the distributions are consistent. DES CCSNe and SLSNe hosts are similar to other samples, while RET hosts are unlike the hosts of any other transients, although these differences have not been disentangled from selection effects. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: 18 pages, 17 figures. Submitted to MNRAS

Report number: Fermilab PUB-20-002-AE; DES-2019-0511

Cosmological Signatures of Dark Sector Physics: The Evolution of Haloes and Spin Alignment

Authors: Absem W. Jibrail, Pascal J. Elahi, Geraint F. Lewis

Abstract: The standard cosmological paradigm currently lacks a detailed account of physics in the dark sector, the dark matter and energy that dominate cosmic evolution. In this paper, we consider the distinguishing factors between three alternative models - warm dark matter, quintessence and coupled dark matter-energy - and \lcdm{} through numerical simulations of cosmological structure formation. Key halo… ▽ More The standard cosmological paradigm currently lacks a detailed account of physics in the dark sector, the dark matter and energy that dominate cosmic evolution. In this paper, we consider the distinguishing factors between three alternative models - warm dark matter, quintessence and coupled dark matter-energy - and \lcdm{} through numerical simulations of cosmological structure formation. Key halo statistics - halo spin/velocity alignment between large-scale structure and neighboring haloes, halo formation time and migration - were compared across cosmologies within the redshift range 0$\leq$$z$$\leq$2.98. We found the alignment of halo motion and spin to large-scale structures and neighbouring haloes to be similar in all cosmologies for a range of redshifts. The search was extended to low density regions, avoiding non-linear disturbances of halo spins, yet very similar alignment trends were found between cosmologies which are difficult to characterize and use as a probe of cosmology. We found haloes in quintessence cosmologies form earlier than their \lcdm{} counterparts. Relating this to the fact that such haloes originate in high density regions, such findings could hold clues to distinguishing factors for the quintessence cosmology from the standard model. Although in general, halo statistics are not an accurate probe of the dark sector physics. △ Less

Submitted 22 December, 2019; originally announced December 2019.

Comments: 15 pages, 13 figures

Voluntary Disclosure and Personalized Pricing

Authors: S. Nageeb Ali, Greg Lewis, Shoshana Vasserman

Abstract: Central to privacy concerns is that firms may use consumer data to price discriminate. A common policy response is that consumers should be given control over which firms access their data and how. Since firms learn about a consumer's preferences based on the data seen and the consumer's disclosure choices, the equilibrium implications of consumer control are unclear. We study whether such measure… ▽ More Central to privacy concerns is that firms may use consumer data to price discriminate. A common policy response is that consumers should be given control over which firms access their data and how. Since firms learn about a consumer's preferences based on the data seen and the consumer's disclosure choices, the equilibrium implications of consumer control are unclear. We study whether such measures improve consumer welfare in monopolistic and competitive markets. We find that consumer control can improve consumer welfare relative to both perfect price discrimination and no personalized pricing. First, consumers can use disclosure to amplify competitive forces. Second, consumers can disclose information to induce even a monopolist to lower prices. Whether consumer control improves welfare depends on the disclosure technology and market competitiveness. Simple disclosure technologies suffice in competitive markets. When facing a monopolist, a consumer needs partial disclosure possibilities to obtain any welfare gains. △ Less

Submitted 17 August, 2020; v1 submitted 10 December, 2019; originally announced December 2019.

A SkyMapper view of the Large Magellanic Cloud: The dynamics of stellar populations

Authors: Zhen Wan, Magda Guglielmo, Geraint F. Lewis, Dougal Mackey, Rodrigo A. Ibata

Abstract: We present the first SkyMapper stellar population analysis of the Large Magellanic Cloud (hereafter LMC),including the identification of 3578 candidate Carbon Stars through their extremely red $g-r$ colours. Coupled with Gaia astrometry, we analyse the distribution and kinematics of this Carbon Star population, finding the LMC to be centred at… ▽ More We present the first SkyMapper stellar population analysis of the Large Magellanic Cloud (hereafter LMC),including the identification of 3578 candidate Carbon Stars through their extremely red $g-r$ colours. Coupled with Gaia astrometry, we analyse the distribution and kinematics of this Carbon Star population, finding the LMC to be centred at $(R.A., Dec.) = (80.90^{\circ}\pm{0.29}, -68.74^{\circ}\pm{0.12})$, with a bulk proper motion of $(μ_α,μ_δ) = (1.878\pm0.007,0.293\pm0.018) \mathrm{mas\ yr^{-1}}$ and a disk inclination of $i = 25.6^{\circ}\pm1.1$ at position angle $θ= 135.6^{\circ}\pm 3.3^{\circ}$. We complement this study with the identification and analysis of additional stellar populations, finding that the dynamical centre for Red Giant Branch (RGB) stars is similar to that seen for the Carbon Stars, whereas for young stars the dynamical centre is significantly offset from the older populations. This potentially indicates that the young stars were formed as a consequence of a strong tidal interaction, probably with the Small Magellanic Cloud (SMC). In terms of internal dynamics, the tangential velocity profile increases linearly within $\sim3\ \mathrm{kpc}$, after which it maintains an approximately constant value of $V_{rot} = 83.6\pm 1.7 \mathrm{km\ s^{-1}}$ until $\sim7 \mathrm{kpc}$. With an asymmetric drift correction, we estimate the mass within $7 \mathrm{kpc}$ to be $M_{\rm LMC}(<7\mathrm{kpc}) = (2.5\pm0.1)\times10^{10}\ \mathrm{M}_{\odot}$ and within the tidal radius ($\sim 30\ \mathrm{kpc}$) to be $M_{\rm LMC}(<30\mathrm{kpc}) = (1.06 \pm 0.32)\times10^{11}\ \mathrm{M}_{\odot}$, consistent with other recent measurements. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: 14 pages, 12 figures. Accepted by MNRAS

A Geometric Probe of Cosmology: I. Gravitational Lensing Time Delays and Quasar Reverberation Mapping

Authors: Angela L. H. Ng, Geraint F. Lewis

Abstract: We present a novel, purely geometric probe of cosmology based on measurements of differential time delays between images of strongly lensed quasars due to finite source effects. Our approach is solely dependent on cosmology via a ratio of angular diameter distances, the image separation, and the source size. It thereby entirely avoids the challenges of lens modelling that conventionally limit time… ▽ More We present a novel, purely geometric probe of cosmology based on measurements of differential time delays between images of strongly lensed quasars due to finite source effects. Our approach is solely dependent on cosmology via a ratio of angular diameter distances, the image separation, and the source size. It thereby entirely avoids the challenges of lens modelling that conventionally limit time delay cosmography, and instead entails the lensed reverberation mapping of the quasar Broad Line Region. We demonstrate that differential time delays are measurable with short cadence spectroscopic monitoring of lensed quasars, through the timing of kinematically identified features within the broad emission lines. This provides a geometric determination of an angular diameter distance ratio complementary to standard probes, and as a result is a potentially powerful new method of constraining cosmology. △ Less

Submitted 6 December, 2019; originally announced December 2019.

Comments: 8 pages, 6 figures, Accepted MNRAS

On the origin of the asymmetric dwarf galaxy distribution around Andromeda

Authors: Zhen Wan, William H. Oliver, Geraint F. Lewis, Justin I. Read, Michelle L. M. Collins

Abstract: The dwarf galaxy distribution surrounding M31 is significantly anisotropic in nature. Of the 30 dwarf galaxies in this distribution, 15 form a disc-like structure and 23 are contained within the hemisphere facing the Milky Way. Using a realistic local potential, we analyse the conditions required to produce and maintain these asymmetries. We find that some dwarf galaxies are required to have highl… ▽ More The dwarf galaxy distribution surrounding M31 is significantly anisotropic in nature. Of the 30 dwarf galaxies in this distribution, 15 form a disc-like structure and 23 are contained within the hemisphere facing the Milky Way. Using a realistic local potential, we analyse the conditions required to produce and maintain these asymmetries. We find that some dwarf galaxies are required to have highly eccentric orbits in order to preserve the presence of the hemispherical asymmetry with an appropriately large radial dispersion. Under the assumption that the dwarf galaxies originate from a single association or accretion event, we find that the initial size and specific energy of that association must both be relatively large in order to produce the observed hemispherical asymmetry. However if the association was large in physical size, the very high-energy required would enable several dwarf galaxies to escape from the M31 and be captured by the Milky Way. Furthermore, we find that associations that result in this structure have total specific energies concentrated around $E = V_{esc}^{2} - V_{init}^{2} \sim 200^2$ -- $300^2\ \rm{km^2\ s^{-2}}$, implying that the initial velocity and initial position needed to produce the structure are strongly correlated. The overlap of initial conditions required to produce the radial dispersion, angular dispersion, and the planar structure is small and suggests that either they did not originate from a single accretion event, or that these asymmetric structures are short-lived. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: 12 pages, 10 figures, accepted by MNRAS

Microlensing and Photon Bunching: The impact of decoherence

Authors: Geraint F. Lewis, Peter Tuthill

Abstract: Gravitational microlensing within the Galaxy offers the prospect of probing the details of distant stellar sources, as well as revealing the distribution of compact (and potentially non-luminous) masses along the line-of-sight. Recently, it has been suggested that additional constraints on the lensing properties can be determined through the measurement of the time delay between images through the… ▽ More Gravitational microlensing within the Galaxy offers the prospect of probing the details of distant stellar sources, as well as revealing the distribution of compact (and potentially non-luminous) masses along the line-of-sight. Recently, it has been suggested that additional constraints on the lensing properties can be determined through the measurement of the time delay between images through the correlation of the bunching of photon arrival times; an application of the Hanbury-Brown Twiss effect. In this paper, we revisit this analysis, examining the impact of decoherence of the radiation from a spatially extended source along the multiple paths to an observer. The result is that, for physically reasonable situations, such decoherence completely erases any correlation that could otherwise be used to measure the gravitational lensing time delay. Indeed, the divergent light paths traverse extremely long effective baselines at the lens plane, corresponding to extremes of angular resolving power well beyond those attainable with any terrestrial technologies; the drawback being that few conceivable celestial objects would be sufficiently compact with high enough surface brightness to yield usable signals. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: 4 pages, 1 figure, accepted for publication in MNRAS

Report number: GFL-001

Quantum Computation with Machine-Learning-Controlled Quantum Stuff

Authors: Lucien Hardy, Adam G. M. Lewis

Abstract: We describe how one may go about performing quantum computation with arbitrary "quantum stuff", as long as it has some basic physical properties. Imagine a long strip of stuff, equipped with regularly spaced wires to provide input settings and to read off outcomes. After showing how the corresponding map from settings to outcomes can be construed as a quantum circuit, we provide a machine learning… ▽ More We describe how one may go about performing quantum computation with arbitrary "quantum stuff", as long as it has some basic physical properties. Imagine a long strip of stuff, equipped with regularly spaced wires to provide input settings and to read off outcomes. After showing how the corresponding map from settings to outcomes can be construed as a quantum circuit, we provide a machine learning algorithm to tomographically "learn" which settings implement the members of a universal gate set. At optimum, arbitrary quantum gates, and thus arbitrary quantum programs, can be implemented using the stuff. △ Less

Submitted 29 November, 2019; originally announced November 2019.

Comments: 13 pages, 3 figures, 1 table, 3 algorithms

Classical Simulations of Quantum Field Theory in Curved Spacetime I: Fermionic Hawking-Hartle Vacua from a Staggered Lattice Scheme

Authors: Adam G. M. Lewis, Guifré Vidal

Abstract: We numerically compute renormalized expectation values of quadratic operators in a quantum field theory (QFT) of free Dirac fermions in curved two-dimensional (Lorentzian) spacetime. First, we use a staggered-fermion discretization to generate a sequence of lattice theories yielding the desired QFT in the continuum limit. Numerically-computed lattice correlators are then used to approximate, throu… ▽ More We numerically compute renormalized expectation values of quadratic operators in a quantum field theory (QFT) of free Dirac fermions in curved two-dimensional (Lorentzian) spacetime. First, we use a staggered-fermion discretization to generate a sequence of lattice theories yielding the desired QFT in the continuum limit. Numerically-computed lattice correlators are then used to approximate, through extrapolation, those in the continuum. Finally, we use so-called point-splitting regularization and Hadamard renormalization to remove divergences, and thus obtain finite, renormalized expectation values of quadratic operators in the continuum. As illustrative applications, we show how to recover the Unruh effect in flat spacetime and how to compute renormalized expectation values in the Hawking-Hartle vacuum of a Schwarzschild black hole and in the Bunch-Davies vacuum of an expanding universe described by de Sitter spacetime. Although here we address a non-interacting QFT using free fermion techniques, the framework described in this paper lays the groundwork for a series of subsequent studies involving simulation of interacting QFTs in curved spacetime by tensor network techniques. △ Less

Submitted 25 August, 2020; v1 submitted 29 November, 2019; originally announced November 2019.

Comments: 29 pages, 9 figures

Journal ref: Quantum 4, 351 (2020)

The GALAH Survey: Temporal Chemical Enrichment of the Galactic Disk

Authors: Jane Lin, Martin Asplund, Yuan-Sen Ting, Luca Casagrande, Sven Buder, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Janez Kos, K Lind, Sarah L. Martell, Sanjib Sharma, Jeffrey D. Simpson, Tomaž Zwitter, Daniel B. Zucker, Ivan Minchev, Klemen Čotar, Michael Hayden, Jonti Horner, Geraint F. Lewis, Thomas Nordlander, Rosemary F. G. Wyse, Maruša Žerjal

Abstract: We present isochrone ages and initial bulk metallicities ($\rm [Fe/H]_{bulk}$, by accounting for diffusion) of 163,722 stars from the GALAH Data Release 2, mainly composed of main sequence turn-off stars and subgiants ($\rm 7000 K>T_{eff}>4000 K$ and $\rm log g>3$ dex). The local age-metallicity relationship (AMR) is nearly flat but with significant scatter at all ages; the scatter is even higher… ▽ More We present isochrone ages and initial bulk metallicities ($\rm [Fe/H]_{bulk}$, by accounting for diffusion) of 163,722 stars from the GALAH Data Release 2, mainly composed of main sequence turn-off stars and subgiants ($\rm 7000 K>T_{eff}>4000 K$ and $\rm log g>3$ dex). The local age-metallicity relationship (AMR) is nearly flat but with significant scatter at all ages; the scatter is even higher when considering the observed surface abundances. After correcting for selection effects, the AMR appear to have intrinsic structures indicative of two star formation events, which we speculate are connected to the thin and thick disks in the solar neighborhood. We also present abundance ratio trends for 16 elements as a function of age, across different $\rm [Fe/H]_{bulk}$ bins. In general, we find the trends in terms of [X/Fe] vs age from our far larger sample to be compatible with studies based on small ($\sim$ 100 stars) samples of solar twins but we now extend it to both sub- and super-solar metallicities. The $α$-elements show differing behaviour: the hydrostatic $α$-elements O and Mg show a steady decline with time for all metallicities while the explosive $α$-elements Si, Ca and Ti are nearly constant during the thin disk epoch (ages $\lessapprox $ 12 Gyr). The s-process elements Y and Ba show increasing [X/Fe] with time while the r-process element Eu have the opposite trend, thus favouring a primary production from sources with a short time-delay such as core-collapse supernovae over long-delay events such as neutron star mergers. △ Less

Submitted 12 November, 2019; originally announced November 2019.

The Canada-France Imaging Survey: Reconstructing the Milky Way Star Formation History from its White Dwarf Population

Authors: Nicholas J. Fantin, Patrick Côté, Alan W. McConnachie, Pierre Bergeron, Jean-Charles Cuillandre, Stephen D. J. Gwyn, Rodrigo A. Ibata, Guillaume F. Thomas, Raymond G. Carlberg, Sébastien Fabbro, Misha Haywood, Ariane Lançon, Geraint F. Lewis, Khyati Malhan, Nicolas F. Martin, Julio F. Navarro, Douglas Scott, Else Starkenburg

Abstract: As the remnants of stars with initial masses $\lesssim$ 8 M$_{\odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over… ▽ More As the remnants of stars with initial masses $\lesssim$ 8 M$_{\odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over $\sim$4500 deg$^2$ using a reduced proper motion diagram. We develop a new white dwarf population synthesis code that returns mock observations of the Galactic field white dwarf population for a given star formation history, while simultaneously taking into account the geometry of the Milky Way, survey parameters, and selection effects. We use this model to derive the star formation histories of the thin disk, thick disk, and stellar halo. Our results show that the Milky Way disk began forming stars (11.3 $\pm$ 0.5) Gyr ago, with a peak rate of (8.8 $\pm$ 1.4) M$_{\odot}$yr$^{-1}$ at (9.8 $\pm$ 0.4) Gyr, before a slow decline to a constant rate until the present day --- consistent with recent results suggesting a merging event with a satellite galaxy. Studying the residuals between the data and best-fit model shows evidence for a slight increase in star formation over the past 3 Gyr. We fit the local fraction of helium-atmosphere white dwarfs to be (21 $\pm$ 3) %. Incorporating this methodology with data from future wide-field surveys such as LSST, Euclid, CASTOR, and WFIRST should provide an unprecedented view into the formation of the Milky Way at its earliest epoch through its white dwarfs. △ Less

Submitted 6 November, 2019; originally announced November 2019.

The GALAH Survey: Chemically tagging the Fimbulthul stream to the globular cluster $ω$ Centauri

Authors: Jeffrey D. Simpson, Sarah L. Martell, Gary Da Costa, Jonathan Horner, Rosemary F. G. Wyse, Yuan-Sen Ting, Martin Asplund, Joss Bland-Hawthorn, Sven Buder, Gayandhi M. De Silva, Ken C. Freeman, Janez Kos, Geraint F. Lewis, Karin Lind, Sanjib Sharma, Daniel B. Zucker, Tomaž Zwitter, Klemen Čotar, Peter L. Cottrell, Thomas Nordlander

Abstract: Using kinematics from Gaia and the large elemental abundance space of the second data release of the GALAH survey, we identify two new members of the Fimbulthul stellar stream, and chemically tag them to massive, multi-metallic globular cluster $ω$ Centauri. Recent analysis of the second data release of Gaia had revealed the Fimbulthul stellar stream in the halo of the Milky Way. It had been propo… ▽ More Using kinematics from Gaia and the large elemental abundance space of the second data release of the GALAH survey, we identify two new members of the Fimbulthul stellar stream, and chemically tag them to massive, multi-metallic globular cluster $ω$ Centauri. Recent analysis of the second data release of Gaia had revealed the Fimbulthul stellar stream in the halo of the Milky Way. It had been proposed that the stream is associated with the $ω$ Centauri, but this proposition relied exclusively upon the kinematics and metallicities of the stars to make the association. In this work, we find our two new members of the stream to be metal-poor stars that are enhanced in sodium and aluminium, typical of second population globular cluster stars, but not otherwise seen in field stars. Furthermore, the stars share the s-process abundance pattern seen in $ω$ Centauri, which is rare in field stars. Apart from one star within 1.5 deg of $ω$ Centauri, we find no other stars observed by GALAH spatially near $ω$ Centauri or the Fimbulthul stream that could be kinematically and chemically linked to the cluster. Chemically tagging stars in the Fimbulthul stream to $ω$ Centauri confirms the earlier work, and further links this tidal feature in the Milky Way halo to $ω$ Centauri. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: 11 pages, 6 figures, accepted to MNRAS

A data-driven model of nucleosynthesis with chemical tagging in a lower-dimensional latent space

Authors: Andrew R. Casey, John C. Lattanzio, Aldeida Aleti, David L. Dowe, Joss Bland-Hawthorn, Sven Buder, Geraint F. Lewis, Sarah L. Martell, Thomas Nordlander, Jeffrey D. Simpson, Sanjib Sharma, Daniel B. Zucker

Abstract: Chemical tagging seeks to identify unique star formation sites from present-day stellar abundances. Previous techniques have treated each abundance dimension as being statistically independent, despite theoretical expectations that many elements can be produced by more than one nucleosynthetic process. In this work we introduce a data-driven model of nucleosynthesis where a set of latent factors (… ▽ More Chemical tagging seeks to identify unique star formation sites from present-day stellar abundances. Previous techniques have treated each abundance dimension as being statistically independent, despite theoretical expectations that many elements can be produced by more than one nucleosynthetic process. In this work we introduce a data-driven model of nucleosynthesis where a set of latent factors (e.g., nucleosynthetic yields) contribute to all stars with different scores, and clustering (e.g., chemical tagging) is modelled by a mixture of multivariate Gaussians in a lower-dimensional latent space. We use an exact method to simultaneously estimate the factor scores for each star, the partial assignment of each star to each cluster, and the latent factors common to all stars, even in the presence of missing data entries. We use an information-theoretic Bayesian principle to estimate the number of latent factors and clusters. Using the second Galah data release we find that six latent factors are preferred to explain N = 2,566 stars with 17 chemical abundances. We identify the rapid- and slow-neutron capture processes, as well as latent factors consistent with Fe-peak and α-element production, and another where K and Zn dominate. When we consider N ~ 160,000 stars with missing abundances we find another 7 factors, as well as 16 components in latent space. Despite these components showing separation in chemistry that is explained through different yield contributions, none show significant structure in their positions or motions. We argue that more data, and joint priors on cluster membership that are constrained by dynamical models, are necessary to realise chemical tagging at a galactic-scale. We release software that allows for model parameters to be optimised in seconds given a fixed number of latent factors, components, and $10^7$ abundance measurements. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: Accepted to ApJ

Demisability and survivability sensitivity to design-for-demise techniques

Authors: Mirko Trisolini, Hugh G. Lewis, Camilla Colombo

Abstract: The paper is concerned with examining the effects that design-for-demise solutions can have not only on the demisability of components, but also on their survivability that is their capability to withstand impacts from space debris. First two models are introduced. A demisability model to predict the behaviour of spacecraft components during the atmospheric re-entry and a survivability model to as… ▽ More The paper is concerned with examining the effects that design-for-demise solutions can have not only on the demisability of components, but also on their survivability that is their capability to withstand impacts from space debris. First two models are introduced. A demisability model to predict the behaviour of spacecraft components during the atmospheric re-entry and a survivability model to assess the vulnerability of spacecraft structures against space debris impacts. Two indices that evaluate the level of demisability and survivability are also proposed. The two models are then used to study the sensitivity of the demisability and of the survivability indices as a function of typical design-for-demise options. The demisability and the survivability can in fact be influenced by the same design parameters in a competing fashion that is while the demisability is improved, the survivability is worsened and vice versa. The analysis shows how the design-for-demise solutions influence the demisability and the survivability independently. In addition, the effect that a solution has simultaneously on the two criteria is assessed. Results shows which, among the design-for-demise parameters mostly influence the demisability and the survivability. For such design parameters maps are presented, describing their influence on the demisability and survivability indices. These maps represent a useful tool to quickly assess the level of demisability and survivability that can be expected from a component, when specific design parameters are changed. △ Less

Submitted 10 October, 2019; originally announced October 2019.

Comments: Accepted for publication in Acta Astronautica

The Pristine Inner Galaxy Survey (PIGS) I: Tracing the kinematics of metal-poor stars in the Galactic bulge

Authors: Anke Arentsen, Else Starkenburg, Nicolas F. Martin, Vanessa Hill, Rodrigo Ibata, Andrea Kunder, Mathias Schultheis, Kim A. Venn, Daniel B. Zucker, David Aguado, Ray Carlberg, Jonay I. González Hernández, Carmela Lardo, Nicolas Longeard, Khyati Malhan, Julio F. Navarro, Ruben Sánchez-Janssen, Federico Sestito, Guillaume Thomas, Kris Youakim, Geraint F. Lewis, Jeffrey D. Simpson, Zhen Wan

Abstract: Our Galaxy is known to contain a central boxy/peanut-shaped bulge, yet the importance of a classical, pressure-supported component within the central part of the Milky Way is still being debated. It should be most visible at low metallicity, a regime that has not yet been studied in detail. Using metallicity-sensitive narrow-band photometry, the Pristine Inner Galaxy Survey (PIGS) has collected a… ▽ More Our Galaxy is known to contain a central boxy/peanut-shaped bulge, yet the importance of a classical, pressure-supported component within the central part of the Milky Way is still being debated. It should be most visible at low metallicity, a regime that has not yet been studied in detail. Using metallicity-sensitive narrow-band photometry, the Pristine Inner Galaxy Survey (PIGS) has collected a large sample of metal-poor ([Fe/H] < -1.0) stars in the inner Galaxy to address this open question. We use PIGS to trace the metal-poor inner Galaxy kinematics as function of metallicity for the first time. We find that the rotational signal decreases with decreasing [Fe/H], until it becomes negligible for the most metal-poor stars. Additionally, the velocity dispersion increases with decreasing metallicity for -3.0 < [Fe/H] < -0.5, with a gradient of -44 $\pm$ 4 km$\,$s$^{-1}\,$dex$^{-1}$. These observations may signal a transition between Galactic components of different metallicities and kinematics, a different mapping onto the boxy/peanut-shaped bulge for former disk stars of different metallicities and/or the secular dynamical and gravitational influence of the bar on the pressure-supported component. Our results provide strong constraints on models that attempt to explain the properties of the inner Galaxy. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: 5 pages + appendices, accepted to MNRAS Letters

Component Mismatches Are a Critical Bottleneck to Fielding AI-Enabled Systems in the Public Sector

Authors: Grace A. Lewis, Stephany Bellomo, April Galyardt

Abstract: The use of machine learning or artificial intelligence (ML/AI) holds substantial potential toward improving many functions and needs of the public sector. In practice however, integrating ML/AI components into public sector applications is severely limited not only by the fragility of these components and their algorithms, but also because of mismatches between components of ML-enabled systems. Fo… ▽ More The use of machine learning or artificial intelligence (ML/AI) holds substantial potential toward improving many functions and needs of the public sector. In practice however, integrating ML/AI components into public sector applications is severely limited not only by the fragility of these components and their algorithms, but also because of mismatches between components of ML-enabled systems. For example, if an ML model is trained on data that is different from data in the operational environment, field performance of the ML component will be dramatically reduced. Separate from software engineering considerations, the expertise needed to field an ML/AI component within a system frequently comes from outside software engineering. As a result, assumptions and even descriptive language used by practitioners from these different disciplines can exacerbate other challenges to integrating ML/AI components into larger systems. We are investigating classes of mismatches in ML/AI systems integration, to identify the implicit assumptions made by practitioners in different fields (data scientists, software engineers, operations staff) and find ways to communicate the appropriate information explicitly. We will discuss a few categories of mismatch, and provide examples from each class. To enable ML/AI components to be fielded in a meaningful way, we will need to understand the mismatches that exist and develop practices to mitigate the impacts of these mismatches. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: Presented at AAAI FSS-19: Artificial Intelligence in Government and Public Sector, Arlington, Virginia, USA

Hadamard Renormalization of a 2-Dimensional Dirac Field

Authors: Adam G. M. Lewis

Abstract: The Hadamard renormalization procedure is applied to a free, massive Dirac field $ψ$ on a 2 dimensional Lorentzian spacetime. This yields the state-independent divergent terms in the Hadamard bispinor $G^{(1)}(x, x') = \frac{1}{2} \left\langle \left[ \barψ(x'), ψ(x) \right] \right\rangle$ as $x$ and $x'$ are brought together along the unique geodesic connecting them. Subtracting these divergent te… ▽ More The Hadamard renormalization procedure is applied to a free, massive Dirac field $ψ$ on a 2 dimensional Lorentzian spacetime. This yields the state-independent divergent terms in the Hadamard bispinor $G^{(1)}(x, x') = \frac{1}{2} \left\langle \left[ \barψ(x'), ψ(x) \right] \right\rangle$ as $x$ and $x'$ are brought together along the unique geodesic connecting them. Subtracting these divergent terms within the limit assigns $G^{(1)}(x, x')$, and thus any operator expressed in terms of it, a finite value at the coincident point $x' = x$. In this limit, one obtains a quadratic operator instead of a bispinor. The procedure is thus used to assign finite values to various quadratic operators, including the stress-energy tensor. Results are presented covariantly, in a conformally-flat coordinate chart at purely spatial separations, and in the Minkowski metric. These terms can be directly subtracted from combinations of $G^{(1)}(x, x')$ - themselves obtained, for example, from a numerical simulation - to obtain finite expectation values defined in the continuum. △ Less

Submitted 7 March, 2020; v1 submitted 12 October, 2019; originally announced October 2019.

Comments: 8 pages. A few corrections made, and a new operator is now renormalized

Journal ref: Phys. Rev. D 101, 125019 (2020)

Spacecraft design optimisation for demise and survivability

Authors: Mirko Trisolini, Hugh G. Lewis, Camilla Colombo

Abstract: Among the mitigation measures introduced to cope with the space debris issue there is the de-orbiting of decommissioned satellites. Guidelines for re-entering objects call for a ground casualty risk no higher than 0.0001. To comply with this requirement, satellites can be designed through a design-for-demise philosophy. Still, a spacecraft designed to demise has to survive the debris-populated spa… ▽ More Among the mitigation measures introduced to cope with the space debris issue there is the de-orbiting of decommissioned satellites. Guidelines for re-entering objects call for a ground casualty risk no higher than 0.0001. To comply with this requirement, satellites can be designed through a design-for-demise philosophy. Still, a spacecraft designed to demise has to survive the debris-populated space environment for many years. The demisability and the survivability of a satellite can both be influenced by a set of common design choices such as the material selection, the geometry definition, and the position of the components. Within this context, two models have been developed to analyse the demise and the survivability of satellites. Given the competing nature of the demisability and the survivability, a multi-objective optimisation framework was developed, with the aim to identify trade-off solutions for the preliminary design of satellites. As the problem is nonlinear and involves the combination of continuous and discrete variables, classical derivative based approaches are unsuited and a genetic algorithm was selected instead. The genetic algorithm uses the developed demisability and survivability criteria as the fitness functions of the multi-objective algorithm. The paper presents a test case, which considers the preliminary optimisation of tanks in terms of material, geometry, location, and number of tanks for a representative Earth observation mission. The configuration of the external structure of the spacecraft is fixed. Tanks were selected because they are sensitive to both design requirements: they represent critical components in the demise process and impact damage can cause the loss of the mission because of leaking and ruptures. The results present the possible trade off solutions, constituting the Pareto front obtained from the multi-objective optimisation. △ Less

Submitted 11 October, 2019; originally announced October 2019.

Comments: Paper accepted for publication in Aerospace Science and Technology

Dwarfs or giants? Stellar metallicities and distances in the Canada-France-Imaging-Survey from $ugrizG$ multi-band photometry

Authors: Guillaume F. Thomas, Nicholaas Annau, Alan McConnachie, Sebastien Fabbro, Hossen Teimoorinia, Patrick Côté, Jean-Charles Cuillandre, Stephen Gwyn, Rodrigo A. Ibata, Else Starkenburg, Raymond Carlberg, Benoit Famaey, Nicholas Fantin, Laura Ferrarese, Vincent Hénault-Brunnet, Jaclyn Jensen, Ariane Lancon, Geraint F. Lewis, Nicolas F. Martin, Julio F. Navarro, Céline Reylé, Rubén Sánchez-Janssen

Abstract: We present a new fully data-driven algorithm that uses photometric data from the Canada-France-Imaging-Survey (CFIS; $u$), Pan-STARRS 1 (PS1; $griz$), and Gaia ($G$) to discriminate between dwarf and giant stars and to estimate their distances and metallicities. The algorithm is trained and tested using the SDSS/SEGUE spectroscopic dataset and Gaia photometric/astrometric dataset. At [Fe/H]… ▽ More We present a new fully data-driven algorithm that uses photometric data from the Canada-France-Imaging-Survey (CFIS; $u$), Pan-STARRS 1 (PS1; $griz$), and Gaia ($G$) to discriminate between dwarf and giant stars and to estimate their distances and metallicities. The algorithm is trained and tested using the SDSS/SEGUE spectroscopic dataset and Gaia photometric/astrometric dataset. At [Fe/H]$<-1.2$, the algorithm succeeds in identifying more than 70% of the giants in the training/test set, with a dwarf contamination fraction below 30% (with respect to the SDSS/SEGUE dataset). The photometric metallicity estimates have uncertainties better than 0.2 dex when compared with the spectroscopic measurements. The distances estimated by the algorithm are valid out to a distance of at least $\sim 80$ kpc without requiring any prior on the stellar distribution, and have fully independent uncertainities that take into account both random and systematic errors. These advances allow us to estimate these stellar parameters for approximately 12 million stars in the photometric dataset. This will enable studies involving the chemical mapping of the distant outer disc and the stellar halo, including their kinematics using the Gaia proper motions. This type of algorithm can be applied in the Southern hemisphere to the first release of LSST data, thus providing an almost complete view of the external components of our Galaxy out to at least $\sim 80$ kpc. Critical to the success of these efforts will be ensuring well-defined spectroscopic training sets that sample a broad range of stellar parameters with minimal biases. A catalogue containing the training/test set and all relevant parameters within the public footprint of CFIS is available online. △ Less

Submitted 9 October, 2019; v1 submitted 7 October, 2019; originally announced October 2019.

Comments: 27 pages. 21 Figures. Accepted for publication in AAS. Training/testing sample available at: http://www.canfar.net/storage/list/gthomas/dwarfOrGiants.fits

The Globular Cluster Population of NGC 1052-DF2: Evidence for Rotation

Authors: Geraint F. Lewis, Brendon J. Brewer, Zhen Wan

Abstract: Based upon the kinematics of ten globular clusters, it has recently been claimed that the ultra-diffuse galaxy, NCD 1052-DF2, lacks a significant quantity of dark matter. Dynamical analyses have generally assumed that this galaxy is pressure supported, with the relatively small velocity dispersion of the globular cluster population indicating the deficit of dark matter. However, the presence of a… ▽ More Based upon the kinematics of ten globular clusters, it has recently been claimed that the ultra-diffuse galaxy, NCD 1052-DF2, lacks a significant quantity of dark matter. Dynamical analyses have generally assumed that this galaxy is pressure supported, with the relatively small velocity dispersion of the globular cluster population indicating the deficit of dark matter. However, the presence of a significant rotation of the globular cluster population could substantially modify this conclusion. Here we present the discovery of such a signature of rotation in the kinematics of NGC 1052-DF2's globular clusters, with a velocity amplitude of $\sim12.44^{+4.40}_{-5.16}$ km/s, which, through Bayesian model comparison, represents a marginally better fit to the available kinematic data; note that this rotation is distinct from, and approximately perpendicular to, the recently identified rotation of the stellar component of NGC 1052-DF2. Assuming this truly represents an underlying rotation, it is shown that the determined mass depends upon the inclination of the rotational component and, with a moderate inclination, the resultant mass to light ratio can exceed $M/L\sim10$. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS Letters

Report number: GFL-001

Two major accretion epochs in M31 from two distinct populations of globular clusters

Authors: Dougal Mackey, Geraint F. Lewis, Brendon J. Brewer, Annette M. N. Ferguson, Jovan Veljanoski, Avon P. Huxor, Michelle L. M. Collins, Patrick Côté, Rodrigo A. Ibata, Mike J. Irwin, Nicolas Martin, Alan W. McConnachie, Jorge Peñarrubia, Nial Tanvir, Zhen Wan

Abstract: Large galaxies grow through the accumulation of dwarf galaxies. In principle it is possible to trace this growth history using the properties of a galaxy's stellar halo. Previous investigations of the galaxy M31 (Andromeda) have shown that outside a radius of 25 kpc the population of halo globular clusters is rotating in alignment with the stellar disk, as are more centrally located clusters. The… ▽ More Large galaxies grow through the accumulation of dwarf galaxies. In principle it is possible to trace this growth history using the properties of a galaxy's stellar halo. Previous investigations of the galaxy M31 (Andromeda) have shown that outside a radius of 25 kpc the population of halo globular clusters is rotating in alignment with the stellar disk, as are more centrally located clusters. The M31 halo also contains coherent stellar substructures, along with a smoothly distributed stellar component. Many of the globular clusters outside 25 kpc are associated with the most prominent substructures, while others are part of the smooth halo. Here we report a new analysis of the kinematics of these globular clusters. We find that the two distinct populations are rotating with perpendicular orientations. The rotation axis for the population associated with the smooth halo is aligned with the rotation axis for the plane of dwarf galaxies that encircles M31. We interpret these separate cluster populations as arising from two major accretion epochs, likely separated by billions of years. Stellar substructures from the first epoch are gone, but those from the more recent second epoch still remain. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: Authors' version of a Letter published in Nature on October 3rd, 2019

A black box for dark sector physics: Predicting dark matter annihilation feedback with conditional GANs

Authors: Florian List, Ishaan Bhat, Geraint F. Lewis

Abstract: Traditionally, incorporating additional physics into existing cosmological simulations requires re-running the cosmological simulation code, which can be computationally expensive. We show that conditional Generative Adversarial Networks (cGANs) can be harnessed to predict how changing the underlying physics alters the simulation results. To illustrate this, we train a cGAN to learn the impact of… ▽ More Traditionally, incorporating additional physics into existing cosmological simulations requires re-running the cosmological simulation code, which can be computationally expensive. We show that conditional Generative Adversarial Networks (cGANs) can be harnessed to predict how changing the underlying physics alters the simulation results. To illustrate this, we train a cGAN to learn the impact of dark matter annihilation feedback (DMAF) on the gas density distribution. The predicted gas density slices are visually difficult to distinguish from their real brethren and the peak counts differ by less than 10 per cent for all test samples (the average deviation is < 3 per cent). Finally, we invert the problem and show that cGANs are capable of endowing smooth density distributions with realistic substructure. The cGAN does however have difficulty generating new knots as well as creating/eliminating bubble-like structures. We conclude that trained cGANs can be an effective approach to provide mock samples of cosmological simulations incorporating DMAF physics from existing samples of standard cosmological simulations of the evolution of cosmic structure. △ Less

Submitted 30 October, 2019; v1 submitted 1 October, 2019; originally announced October 2019.

Comments: 11 pages, 9 figures

Journal ref: MNRAS 490(3), 3134 (2019)

Comparing the Quenching Times of Faint M31 and Milky Way Satellite Galaxies

Authors: Daniel R. Weisz, Nicolas F. Martin, Andrew E. Dolphin, Saundra M. Albers, Michelle L. M. Collins, Annette M. N. Ferguson, Geraint F. Lewis, 8 Dougal Mackey, Alan McConnachie, R. Michael Rich, Evan D. Skillman

Abstract: We present the star formation histories (SFHs) of 20 faint M31 satellites ($-12 \lesssim M_V \lesssim -6$) that were measured by modeling sub-horizontal branch (HB) depth color-magnitude diagrams constructed from Hubble Space Telescope (HST) imaging. Reinforcing previous results, we find that virtually all galaxies quenched between 3 and 9 Gyr ago, independent of luminosity, with a notable concent… ▽ More We present the star formation histories (SFHs) of 20 faint M31 satellites ($-12 \lesssim M_V \lesssim -6$) that were measured by modeling sub-horizontal branch (HB) depth color-magnitude diagrams constructed from Hubble Space Telescope (HST) imaging. Reinforcing previous results, we find that virtually all galaxies quenched between 3 and 9 Gyr ago, independent of luminosity, with a notable concentration $3-6$ Gyr ago. This is in contrast to the Milky Way (MW) satellites, which are generally either faint with ancient quenching times or luminous with recent ($<3$ Gyr) quenching times. We suggest that systematic differences in the quenching times of M31 and MW satellites may be a reflection of the varying accretion histories of M31 and the MW. This result implies that the formation histories of low-mass satellites may not be broadly representative of low-mass galaxies in general. Among the M31 satellite population we identify two distinct groups based on their SFHs: one with exponentially declining SFHs ($τ\sim 2$ Gyr) and one with rising SFHs with abrupt quenching. We speculate how these two groups could be related to scenarios for a recent major merger involving M31. The Cycle 27 HST Treasury survey of M31 satellites will provide well-constrained ancient SFHs to go along with the quenching times we measure here. The discovery and characterization of M31 satellites with $M_V \gtrsim -6$ would help quantify the relative contributions of reionization and environment to quenching of the lowest-mass satellites. △ Less

Submitted 24 September, 2019; originally announced September 2019.

Comments: 7 pages, 3 figures, 1 table. Submitted to ApJL and posted after addressing the referee's report. Comments welcome

A Dwarf Disrupting -- Andromeda XXVII and the North West Stream

Authors: Janet Preston, Michelle L. M. Collins, Rodrigo A. Ibata, Erik J. Tollerud, R. Michael Rich, Ana Bonaca, Alan W. McConnachie, Dougal Mackey, Geraint F. Lewis, Nicolas F. Martin, Jorge Peñarrubia, Scott C. Chapman, Maxime Delorme

Abstract: We present a kinematic and spectroscopic analysis of 38 red giant branch stars, in 7 fields, spanning the dwarf spheroidal galaxy Andromeda XXVII and the upper segment of the North West Stream. Both features are located in the outer halo of the Andromeda galaxy at a projected radius of 50-80 kpc, with the stream extending for $\sim$3$^{\circ}$ on the sky. Our data is obtained as part of the PAndAS… ▽ More We present a kinematic and spectroscopic analysis of 38 red giant branch stars, in 7 fields, spanning the dwarf spheroidal galaxy Andromeda XXVII and the upper segment of the North West Stream. Both features are located in the outer halo of the Andromeda galaxy at a projected radius of 50-80 kpc, with the stream extending for $\sim$3$^{\circ}$ on the sky. Our data is obtained as part of the PAndAS survey and enables us to confirm that Andromeda XXVII's heliocentric distance is 827 $\pm$ 47 kpc and spectroscopic metallicity is -2.1$^{+0.4}_{-0.5}$. We also re-derive Andromeda XXVII's kinematic properties, measuring a systemic velocity = -526.1$^{+10.0}_{-11.0}${\kms} and a velocity dispersion that we find to be non-Gaussian but for which we derive a formal value of 27.0$^{+2.2}_{-3.9}${\kms}. In the upper segment of the North West Stream we measure mean values for the metallicity = -1.8$\pm$0.4, systemic velocity = -519.4 $\pm$4.0{\kms} and velocity dispersion = 10.0$\pm$4.0{\kms}. We also detect a velocity gradient of 1.7$\pm$0.3 {\kms} kpc$^{-1}$ on an infall trajectory towards M31. With a similar gradient, acting in the same direction, in the lower segment we suggest that the North West Stream is not a single structure. As the properties of the upper segment of the North West Stream and Andromeda XXVII are consistent within ~90\% confidence limits, it is likely that the two are related and plausible that Andromeda XXVII is the progenitor of this stream. △ Less

Submitted 20 September, 2019; originally announced September 2019.

Comments: 15 pages, 10 figures

A rogues gallery of Andromeda's dwarf galaxies II. Precise Distances to 17 Faint Satellites

Authors: Daniel R. Weisz, Andrew E. Dolphin, Nicolas F. Martin, Saundra M. Albers, Michelle L. M. Collins, Annette M. N. Ferguson, Geraint F. Lewis, A. Dougal Mackey, Alan McConnachie, R. Michael Rich, Evan D. Skillman

Abstract: We present new horizontal branch (HB) distance measurements to 17 of the faintest known M31 satellites ($-6 \lesssim M_{V} \lesssim -13$) based on deep Hubble Space Telescope (HST) imaging. The color-magnitude diagrams extend $\sim$1-2 magnitudes below the HB, which provides for well-defined HBs, even for faint galaxies in which the tip of the red giant branch (TRGB) is sparsely populated. We dete… ▽ More We present new horizontal branch (HB) distance measurements to 17 of the faintest known M31 satellites ($-6 \lesssim M_{V} \lesssim -13$) based on deep Hubble Space Telescope (HST) imaging. The color-magnitude diagrams extend $\sim$1-2 magnitudes below the HB, which provides for well-defined HBs, even for faint galaxies in which the tip of the red giant branch (TRGB) is sparsely populated. We determine distances across the sample to an average precision of 4% ($\sim 30$~kpc at $800$~kpc). We find that the majority of these galaxies are in good agreement, though slightly farther (0.1-0.2 mag) when compared to recent ground-based TRGB distances. Two galaxies (And~IX and And~XVII) have discrepant HST and ground-based distances by $\sim 0.3$ mag ($\sim 150$~kpc), which may be due to contamination from Milky Way foreground stars and/or M31 halo stars in sparsely populated TRGB regions. We use the new distances to update the luminosities and structural parameters for these 17 M31 satellites. The new distances do not substantially change the spatial configuration of the M31 satellite system. We comment on future prospects for precise and accurate HB distances for faint galaxies in the Local Group and beyond. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: 8 pages, 3 figures. Accepted for publication in MNRAS

Is there a Cosmological Basis for E = mc^2?

Authors: Geraint F. Lewis

Abstract: It has recently been claimed that relativity's most famous equation, E = mc^2, has a cosmological basis, representing the gravitational binding energy for a particle to escape from the origin to a gravitational horizon of the universe. In this paper, I examine these claims in detail, concluding that they result from a misinterpretation of motion of particles in the cosmological space-time, and an… ▽ More It has recently been claimed that relativity's most famous equation, E = mc^2, has a cosmological basis, representing the gravitational binding energy for a particle to escape from the origin to a gravitational horizon of the universe. In this paper, I examine these claims in detail, concluding that they result from a misinterpretation of motion of particles in the cosmological space-time, and an incorrect application of 4-vectors. Finally, I demonstrate that the origin of E = mc^2 comes from its usual relativistic interpretation, namely that it is the energy of a particle as seen in its own rest-frame. △ Less

Submitted 25 August, 2019; originally announced August 2019.

Comments: 6 pages, accepted for publication in General Relativity and Gravitation

Report number: GFL-001

A Novel Scheme for Dark Matter Annihilation Feedback in Cosmological Simulations

Authors: Florian List, Nikolas Iwanus, Pascal J. Elahi, Geraint F. Lewis

Abstract: We present a new self-consistent method for incorporating dark matter annihilation feedback (DMAF) in cosmological N-body simulations. The power generated by DMAF is evaluated at each dark matter (DM) particle which allows for flexible energy injection into the surrounding gas based on the specific DM annihilation model under consideration. Adaptive, individual time steps for gas and DM particles… ▽ More We present a new self-consistent method for incorporating dark matter annihilation feedback (DMAF) in cosmological N-body simulations. The power generated by DMAF is evaluated at each dark matter (DM) particle which allows for flexible energy injection into the surrounding gas based on the specific DM annihilation model under consideration. Adaptive, individual time steps for gas and DM particles are supported and a new time-step limiter, derived from the propagation of a Sedov--Taylor blast wave, is introduced. We compare this donor-based approach with a receiver-based approach used in recent studies and illustrate the differences by means of a toy example. Furthermore, we consider an isolated halo and a cosmological simulation and show that for these realistic cases, both methods agree well with each other. The extension of our implementation to scenarios such as non-local energy injection, velocity-dependent annihilation cross-sections, and DM decay is straightforward. △ Less

Submitted 19 September, 2019; v1 submitted 15 August, 2019; originally announced August 2019.

Comments: 18 pages, 10 (+2) figures

Journal ref: MNRAS 489(3), 4217 (2019)

Probing the independence within the dark sector in the fluid approximation

Authors: Lawrence Dam, Krzysztof Bolejko, Geraint F. Lewis

Abstract: The standard model of cosmology is based on two unknown dark components that are uncoupled from each other. In this paper we investigate whether there is evidence for an interaction between these components of cold dark matter (CDM) and dark energy (DE). In particular, we focus on a minimal extension and reconstruct the interaction history at low-redshifts non-parametrically using a variation of t… ▽ More The standard model of cosmology is based on two unknown dark components that are uncoupled from each other. In this paper we investigate whether there is evidence for an interaction between these components of cold dark matter (CDM) and dark energy (DE). In particular, we focus on a minimal extension and reconstruct the interaction history at low-redshifts non-parametrically using a variation of the commonly used principal component analysis. Although we focus on the interaction in the dark sector, any significant deviation from the standard model that changes the expansion history of the Universe, should leave imprints detectable by our analysis. Thus, detecting signatures of interaction could also be indicative of other non-standard phenomena even if they are not the results of the interaction. It is thus interesting to note that the results presented in this paper do not provide support for the interaction in the dark sector, although the uncertainty is still quite large. In so far as interaction is present but undetectable using current data, we show from a Fisher forecast that forthcoming LSST and DESI surveys will be able to constrain a DM-DE coupling at $20\%$ precision --- enough to falsify the non-interacting scenario, assuming the presence of a modest amount of interaction. △ Less

Submitted 25 November, 2019; v1 submitted 6 August, 2019; originally announced August 2019.

Comments: 33 pages, 11 figures, 1 table; v2: minor changes and additions, conclusions unchanged, accepted for publication in JCAP

Journal ref: JCAP12(2019)030

The Great Escape: Discovery of a nearby 1700 km/s star ejected from the Milky Way by Sgr A*

Authors: Sergey E. Koposov, Douglas Boubert, Ting S. Li, Denis Erkal, Gary S. Da Costa, Daniel B. Zucker, Alexander P. Ji, Kyler Kuehn, Geraint F. Lewis, Dougal Mackey, Jeffrey D. Simpson, Nora Shipp, Zhen Wan, Vasily Belokurov, Joss Bland-Hawthorn, Sarah L. Martell, Thomas Nordlander, Andrew B. Pace, Gayandhi M. De Silva, Mei-Yu Wang

Abstract: We present the serendipitous discovery of the fastest Main Sequence hyper-velocity star (HVS) by the Southern Stellar Stream Spectroscopic Survey (S5). The star S5-HVS1 is a $\sim 2.35$ M$_\odot$ A-type star located at a distance of $\sim 9$ kpc from the Sun and has a heliocentric radial velocity of $1017\pm 2.7$ km/s without any signature of velocity variability. The current 3-D velocity of the s… ▽ More We present the serendipitous discovery of the fastest Main Sequence hyper-velocity star (HVS) by the Southern Stellar Stream Spectroscopic Survey (S5). The star S5-HVS1 is a $\sim 2.35$ M$_\odot$ A-type star located at a distance of $\sim 9$ kpc from the Sun and has a heliocentric radial velocity of $1017\pm 2.7$ km/s without any signature of velocity variability. The current 3-D velocity of the star in the Galactic frame is $1755\pm50$ km/s. When integrated backwards in time, the orbit of the star points unambiguously to the Galactic Centre, implying that S5-HVS1 was kicked away from Sgr A* with a velocity of $\sim 1800$ km/s and travelled for $4.8$ Myr to its current location. This is so far the only HVS confidently associated with the Galactic Centre. S5-HVS1 is also the first hyper-velocity star to provide constraints on the geometry and kinematics of the Galaxy, such as the Solar motion $V_{y,\odot}= 246.1\pm 5.3$ km/s or position $R_0=8.12\pm 0.23$ kpc. The ejection trajectory and transit time of S5-HVS1 coincide with the orbital plane and age of the annular disk of young stars at the Galactic centre, and thus may be linked to its formation. With the S5-HVS1 ejection velocity being almost twice the velocity of other hyper-velocity stars previously associated with the Galactic Centre, we question whether they have been generated by the same mechanism or whether the ejection velocity distribution has been constant over time. △ Less

Submitted 11 November, 2019; v1 submitted 26 July, 2019; originally announced July 2019.

Comments: Accepted to MNRAS; 16 pages

Proper Motions of Stellar Streams Discovered in the Dark Energy Survey

Authors: N. Shipp, T. S. Li, A. B. Pace, D. Erkal, A. Drlica-Wagner, B. Yanny, V. Belokurov, W. Wester, S. E. Koposov, G. F. Lewis, J. D. Simpson, Z. Wan, D. B. Zucker, S. L. Martell, M. Y. Wang

Abstract: We cross-match high-precision astrometric data from Gaia DR2 with accurate multi-band photometry from the Dark Energy Survey (DES) DR1 to confidently measure proper motions for nine stellar streams in the DES footprint: Aliqa Uma, ATLAS, Chenab, Elqui, Indus, Jhelum, Phoenix, Tucana III, and Turranburra. We determine low-confidence proper motion measurements for four additional stellar streams: Ra… ▽ More We cross-match high-precision astrometric data from Gaia DR2 with accurate multi-band photometry from the Dark Energy Survey (DES) DR1 to confidently measure proper motions for nine stellar streams in the DES footprint: Aliqa Uma, ATLAS, Chenab, Elqui, Indus, Jhelum, Phoenix, Tucana III, and Turranburra. We determine low-confidence proper motion measurements for four additional stellar streams: Ravi, Wambelong, Willka Yaku, and Turbio. We find evidence for a misalignment between stream tracks and the systemic proper motion of streams that may suggest a systematic gravitational influence from the Large Magellanic Cloud. These proper motions, when combined with radial velocity measurements, will allow for detailed orbit modeling which can be used to constrain properties of the LMC and its affect on nearby streams, as well as global properties of the Milky Way's gravitational potential. △ Less

Submitted 22 July, 2019; originally announced July 2019.

Comments: 22 pages, 9 figures, 6 tables; submitted to ApJ

The Southern Stellar Stream Spectroscopic Survey (${S}^5$): Overview, Target Selection, Data Reduction, Validation, and Early Science

Authors: T. S. Li, S. E. Koposov, D. B. Zucker, G. F. Lewis, K. Kuehn, J. D. Simpson, A. P. Ji, N. Shipp, Y. -Y. Mao, M. Geha, A. B. Pace, A. D. Mackey, S. Allam, D. L. Tucker, G. S. Da Costa, D. Erkal, J. D. Simon, J. R. Mould, S. L. Martell, Z. Wan, G. M. De Silva, K. Bechtol, E. Balbinot, V. Belokurov, J. Bland-Hawthorn , et al. (7 additional authors not shown)

Abstract: We introduce the Southern Stellar Stream Spectroscopy Survey (${S}^5$), an on-going program to map the kinematics and chemistry of stellar streams in the Southern Hemisphere. The initial focus of ${S}^5$ has been spectroscopic observations of recently identified streams within the footprint of the Dark Energy Survey (DES), with the eventual goal of surveying streams across the entire southern sky.… ▽ More We introduce the Southern Stellar Stream Spectroscopy Survey (${S}^5$), an on-going program to map the kinematics and chemistry of stellar streams in the Southern Hemisphere. The initial focus of ${S}^5$ has been spectroscopic observations of recently identified streams within the footprint of the Dark Energy Survey (DES), with the eventual goal of surveying streams across the entire southern sky. Stellar streams are composed of material that has been tidally striped from dwarf galaxies and globular clusters and hence are excellent dynamical probes of the gravitational potential of the Milky Way, as well as providing a detailed snapshot of its accretion history. Observing with the 3.9-m Anglo-Australian Telescope's 2-degree-Field fibre positioner and AAOmega spectrograph, and combining the precise photometry of DES DR1 with the superb proper motions from $Gaia$ DR2, allows us to conduct an efficient spectroscopic survey to map these stellar streams. So far ${S}^5$ has mapped 9 DES streams and 3 streams outside of DES; the former are the first spectroscopic observations of these recently discovered streams. In addition to the stream survey, we use spare fibres to undertake a Milky Way halo survey and a low-redshift galaxy survey. This paper presents an overview of the ${S}^5$ program, describing the scientific motivation for the survey, target selection, observation strategy, data reduction and survey validation. Finally, we describe early science results on stellar streams and Milky Way halo stars drawn from the survey. Updates on ${S}^5$, including future public data release, can be found at \url{http://s5collab.github.io}. △ Less

Submitted 11 November, 2019; v1 submitted 22 July, 2019; originally announced July 2019.

Comments: 25 pages, 14 figures (1 in appendix), 3 tables (1 in appendix). Published on MNRAS. See also paper from Shipp et al. 2019, which measures the proper motion of the DES streams

Direct detection of the cosmic expansion: the redshift drift and the flux drift

Authors: Krzysztof Bolejko, Chengyi Wang, Geraint F. Lewis

Abstract: The redshift drift, the change of cosmological redshift with time, is a direct consequence of the expansion of the Universe. Thus the measurement of the cosmological redshift drift will offer a direct test of our models of cosmology. The magnitude of the effect is very small, i.e. the spectral shift is of order of $10^{-10} - 10^{-9}$ over the period of a decade, but the next generation facilities… ▽ More The redshift drift, the change of cosmological redshift with time, is a direct consequence of the expansion of the Universe. Thus the measurement of the cosmological redshift drift will offer a direct test of our models of cosmology. The magnitude of the effect is very small, i.e. the spectral shift is of order of $10^{-10} - 10^{-9}$ over the period of a decade, but the next generation facilities such as ELT and SKA will be able to directly detect the expansion of our Universe by the year 2040. In this paper we focus on detectebility of this effect, including strategies of overcoming the kinematic contamination of the cosmological signal. We also show that the redshift drift directly impacts the change of flux. Thus apart from the redshift drift, measurements of the flux drift will provide an additional tool of detecting the expansion of the universe, including its acceleration. We discuss the strategies of detecting the flux drift and show that by including the flux drift signal to the redshift drift signal we boost the chances of a direct detection of the expansion of the Universe. We show that if only the stability of flux is at the level of $ΔF/F \sim 10^{-6}$ then the SKA1-mid Array should be able to detect these effects, before the ETL and the full SKA. Thus, by including the flux drift into the SKA1-mid Array's analysis pipeline, we could be able to provide by mid-2030s a direct evidence of the expansion of the universe including its accelerating phase. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 14 pages, 7 figures

Machine Learning Estimation of Heterogeneous Treatment Effects with Instruments

Authors: Vasilis Syrgkanis, Victor Lei, Miruna Oprescu, Maggie Hei, Keith Battocchi, Greg Lewis

Abstract: We consider the estimation of heterogeneous treatment effects with arbitrary machine learning methods in the presence of unobserved confounders with the aid of a valid instrument. Such settings arise in A/B tests with an intent-to-treat structure, where the experimenter randomizes over which user will receive a recommendation to take an action, and we are interested in the effect of the downstream… ▽ More We consider the estimation of heterogeneous treatment effects with arbitrary machine learning methods in the presence of unobserved confounders with the aid of a valid instrument. Such settings arise in A/B tests with an intent-to-treat structure, where the experimenter randomizes over which user will receive a recommendation to take an action, and we are interested in the effect of the downstream action. We develop a statistical learning approach to the estimation of heterogeneous effects, reducing the problem to the minimization of an appropriate loss function that depends on a set of auxiliary models (each corresponding to a separate prediction task). The reduction enables the use of all recent algorithmic advances (e.g. neural nets, forests). We show that the estimated effect model is robust to estimation errors in the auxiliary models, by showing that the loss satisfies a Neyman orthogonality criterion. Our approach can be used to estimate projections of the true effect model on simpler hypothesis spaces. When these spaces are parametric, then the parameter estimates are asymptotically normal, which enables construction of confidence sets. We applied our method to estimate the effect of membership on downstream webpage engagement on TripAdvisor, using as an instrument an intent-to-treat A/B test among 4 million TripAdvisor users, where some users received an easier membership sign-up process. We also validate our method on synthetic data and on public datasets for the effects of schooling on income. △ Less

Submitted 5 June, 2019; v1 submitted 24 May, 2019; originally announced May 2019.

Semi-Parametric Efficient Policy Learning with Continuous Actions

Authors: Mert Demirer, Vasilis Syrgkanis, Greg Lewis, Victor Chernozhukov

Abstract: We consider off-policy evaluation and optimization with continuous action spaces. We focus on observational data where the data collection policy is unknown and needs to be estimated. We take a semi-parametric approach where the value function takes a known parametric form in the treatment, but we are agnostic on how it depends on the observed contexts. We propose a doubly robust off-policy estima… ▽ More We consider off-policy evaluation and optimization with continuous action spaces. We focus on observational data where the data collection policy is unknown and needs to be estimated. We take a semi-parametric approach where the value function takes a known parametric form in the treatment, but we are agnostic on how it depends on the observed contexts. We propose a doubly robust off-policy estimate for this setting and show that off-policy optimization based on this estimate is robust to estimation errors of the policy function or the regression model. Our results also apply if the model does not satisfy our semi-parametric form, but rather we measure regret in terms of the best projection of the true value function to this functional space. Our work extends prior approaches of policy optimization from observational data that only considered discrete actions. We provide an experimental evaluation of our method in a synthetic data example motivated by optimal personalized pricing and costly resource allocation. △ Less

Submitted 20 July, 2019; v1 submitted 24 May, 2019; originally announced May 2019.