Showing 1–50 of 115 results for author: Harvey, D

A deep-learning algorithm to disentangle self-interacting dark matter and AGN feedback models

Authors: David Harvey

Abstract: Different models of dark matter can alter the distribution of mass in galaxy clusters in a variety of ways. However, so can uncertain astrophysical feedback mechanisms. Here we present a Machine Learning method that ''learns'' how the impact of dark matter self-interactions differs from that of astrophysical feedback in order to break this degeneracy and make inferences on dark matter. We train a… ▽ More Different models of dark matter can alter the distribution of mass in galaxy clusters in a variety of ways. However, so can uncertain astrophysical feedback mechanisms. Here we present a Machine Learning method that ''learns'' how the impact of dark matter self-interactions differs from that of astrophysical feedback in order to break this degeneracy and make inferences on dark matter. We train a Convolutional Neural Network on images of galaxy clusters from hydro-dynamic simulations. In the idealised case our algorithm is 80% accurate at identifying if a galaxy cluster harbours collisionless dark matter, dark matter with $σ_{\rm DM}/m = 0.1$cm$^2/$g or with $σ_{DM}/m = 1$cm$^2$/g. Whilst we find adding X-ray emissivity maps does not improve the performance in differentiating collisional dark matter, it does improve the ability to disentangle different models of astrophysical feedback. We include noise to resemble data expected from Euclid and Chandra and find our model has a statistical error of < 0.01cm$^2$/g and that our algorithm is insensitive to shape measurement bias and photometric redshift errors. This method represents a new way to analyse data from upcoming telescopes that is an order of magnitude more precise and many orders faster, enabling us to explore the dark matter parameter space like never before. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: Accepted Nature Astronomy

Hydrodynamical simulations of merging galaxy clusters: giant dark matter particle colliders, powered by gravity

Authors: Ellen L. Sirks, David Harvey, Richard Massey, Kyle A. Oman, Andrew Robertson, Carlos Frenk, Spencer Everett, Ajay S. Gill, David Lagattuta, Jacqueline McCleary

Abstract: Terrestrial particle accelerators collide charged particles, then watch the trajectory of outgoing debris - but they cannot manipulate dark matter. Fortunately, dark matter is the main component of galaxy clusters, which are continuously pulled together by gravity. We show that galaxy cluster mergers can be exploited as enormous, natural dark matter colliders. We analyse hydrodynamical simulations… ▽ More Terrestrial particle accelerators collide charged particles, then watch the trajectory of outgoing debris - but they cannot manipulate dark matter. Fortunately, dark matter is the main component of galaxy clusters, which are continuously pulled together by gravity. We show that galaxy cluster mergers can be exploited as enormous, natural dark matter colliders. We analyse hydrodynamical simulations of a universe containing self-interacting dark matter (SIDM) in which all particles interact via gravity, and dark matter particles can also scatter off each other via a massive mediator. During cluster collisions, SIDM spreads out and lags behind cluster member galaxies. Individual systems can have quirky dynamics that makes them difficult to interpret. Statistically, however, we find that the mean or median of dark matter's spatial offset in many collisions can be robustly modelled, and is independent of our viewing angle and halo mass even in collisions between unequal-mass systems. If the SIDM cross-section were sigma/m = 0.1cm^2/g = 0.18 barn/GeV, the 'bulleticity' lag would be ~5 percent that of gas due to ram pressure, and could be detected at 95 percent confidence in weak lensing observations of ~100 well-chosen clusters. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 8 pages, 5 figures plus appendices. MNRAS in press

Confidence Intervals on Multivariate Normal Quantiles for Environmental Specification Development in Multi-axis Shock and Vibration Testing

Authors: Adam Watts, Thomas Thompson, Dustin Harvey

Abstract: This article describes two Monte Carlo methods for calculating confidence intervals on cumulative density function (CDF) based multivariate normal quantiles that allows for controlling the tail regions of a multivariate distribution where one is most concerned about extreme responses. The CDF based multivariate normal quantiles associated with bivariate distributions are represented as contours an… ▽ More This article describes two Monte Carlo methods for calculating confidence intervals on cumulative density function (CDF) based multivariate normal quantiles that allows for controlling the tail regions of a multivariate distribution where one is most concerned about extreme responses. The CDF based multivariate normal quantiles associated with bivariate distributions are represented as contours and for trivariate distributions represented as iso-surfaces. We first provide a novel methodology for an inverse problem, characterizing the uncertainty on the $τ^{\mathrm{th}}$ multivariate quantile probability, when using concurrent univariate quantile probabilities. The uncertainty on the $τ^{\mathrm{th}}$ multivariate quantile probability demonstrates inadequacy in univariate methods which neglect correlation between multiple variates. Limitations of traditional multivariate normal tolerance regions and simultaneous univariate tolerance methods are discussed thereby necessitating the need for confidence intervals on CDF based multivariate normal quantiles. Two Monte Carlo methods are discussed; the first calculates the CDF over a tessellated domain followed by taking a bootstrap confidence interval over the tessellated CDF. The CDF based multivariate quantiles are then estimated from the CDF confidence intervals. For the second method, only the point associated with highest probability density along the CDF based quantile is calculated, which greatly improves the computational speed compared to the first method. Monte Carlo simulation studies are used to assess the performance of the various methods. Finally, real data analysis is performed to illustrate a workflow for CDF based multivariate normal quantiles in the domain of mechanical shock and vibration to specify a minimum conservative test level for environmental specification. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 37 pages, 14 figures, 5 tables not including appendices. 45 pages in total

Water-based Quantum Dots Liquid Scintillator for Particle Physics

Authors: M. Zhao, M. Taani, J. Cole, B. Crudele, B. Zou, N. Bhuiyan, E. Chowdhury, Y. Duan, S. Fekri, D. Harvey, D. Mitra, O. Raz, A. Thompson, T. Katori, A. Rakovich

Abstract: Liquid scintillators are typically composed from organic compounds dissolved in organic solvents. However, usage of such material is often restricted due to fire safety and environmental reasons. Because of this, R\&D of water-based liquid scintillators is of extreme relevance; yet, no such scintillators have been made commercially available as yet. Here, we investigate an alternative, water-based… ▽ More Liquid scintillators are typically composed from organic compounds dissolved in organic solvents. However, usage of such material is often restricted due to fire safety and environmental reasons. Because of this, R\&D of water-based liquid scintillators is of extreme relevance; yet, no such scintillators have been made commercially available as yet. Here, we investigate an alternative, water-based quantum dots liquid scintillator. Pre-determined and controllable optical properties of the quantum dots, as well as the existence of large libraries of established protocols for their dispersion in aqueous solutions, make them an attractive option for nuclear and particle physics applications. We characterize the optical properties of water-based quantum dots liquid scintillator and find that most of its optical properties are preserved upon quantum dots' phase transfer into water, through the addition of an oleic acid hydrophilic layer. Using the developed scintillator, the time and charge responses from atmospheric muons are measured, highlighting the practical viability of water-based quantum dots liquid scintillators for nuclear and particle physics, special interest on neutrino physics. △ Less

Submitted 26 June, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comments: 23 pages, 9 figures, 3 tables, v3, accepted by JINST for publication

Weak gravitational lensing measurements of Abell 2744 using JWST and shear measurement algorithm pyRRG-JWST

Authors: David R. Harvey, Richard Massey

Abstract: We update the publicly available weak lensing shear measurement algorithm pyRRG for the James Webb Space Telescope, and apply it to UNCOVER DR1 imaging of galaxy cluster Abell 2744. At short wavelengths (< 2.5$μ$m), shear measurements are consistent between independent observations through different JWST bandpasses, and calibrated within 1.5% of those from the Hubble Space Telescope. At longer wav… ▽ More We update the publicly available weak lensing shear measurement algorithm pyRRG for the James Webb Space Telescope, and apply it to UNCOVER DR1 imaging of galaxy cluster Abell 2744. At short wavelengths (< 2.5$μ$m), shear measurements are consistent between independent observations through different JWST bandpasses, and calibrated within 1.5% of those from the Hubble Space Telescope. At longer wavelengths, shear is underestimated by ~5%, probably due to coarser pixellisation. We model the spatially varying Point Spread Function (PSF) using WebbPSF, whose moments are within 0.05 of real stars near the centre of the mosaic, where there are sufficient stars to also generate an empirical model. We measure shear from up to 162 galaxies arcminute$^2$ to derive a map of (dark plus baryonic) mass with 12 arcsecond (55 kpc) spatial resolution. All code, catalogues and maps are available from https://github.com/davidharvey1986/pyRRG. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: Accepted MNRAS

MinePlanner: A Benchmark for Long-Horizon Planning in Large Minecraft Worlds

Authors: William Hill, Ireton Liu, Anita De Mello Koch, Damion Harvey, Nishanth Kumar, George Konidaris, Steven James

Abstract: We propose a new benchmark for planning tasks based on the Minecraft game. Our benchmark contains 45 tasks overall, but also provides support for creating both propositional and numeric instances of new Minecraft tasks automatically. We benchmark numeric and propositional planning systems on these tasks, with results demonstrating that state-of-the-art planners are currently incapable of dealing w… ▽ More We propose a new benchmark for planning tasks based on the Minecraft game. Our benchmark contains 45 tasks overall, but also provides support for creating both propositional and numeric instances of new Minecraft tasks automatically. We benchmark numeric and propositional planning systems on these tasks, with results demonstrating that state-of-the-art planners are currently incapable of dealing with many of the challenges advanced by our new benchmark, such as scaling to instances with thousands of objects. Based on these results, we identify areas of improvement for future planners. Our framework is made available at https://github.com/IretonLiu/mine-pddl/. △ Less

Submitted 28 April, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: Accepted to the 6th ICAPS Workshop on the International Planning Competition (WIPC 2024)

Data downloaded via parachute from a NASA super-pressure balloon

Authors: Ellen L. Sirks, Richard Massey, Ajay S. Gill, Jason Anderson, Steven J. Benton, Anthony M. Brown, Paul Clark, Joshua English, Spencer W. Everett, Aurelien A. Fraisse, Hugo Franco, John W. Hartley, David Harvey, Bradley Holder, Andrew Hunter, Eric M. Huff, Andrew Hynous, Mathilde Jauzac, William C. Jones, Nikky Joyce, Duncan Kennedy, David Lagattuta, Jason S. -Y. Leung, Lun Li, Stephen Lishman , et al. (18 additional authors not shown)

Abstract: In April to May 2023, the superBIT telescope was lifted to the Earth's stratosphere by a helium-filled super-pressure balloon, to acquire astronomical imaging from above (99.5% of) the Earth's atmosphere. It was launched from New Zealand then, for 40 days, circumnavigated the globe five times at a latitude 40 to 50 degrees South. Attached to the telescope were four 'DRS' (Data Recovery System) cap… ▽ More In April to May 2023, the superBIT telescope was lifted to the Earth's stratosphere by a helium-filled super-pressure balloon, to acquire astronomical imaging from above (99.5% of) the Earth's atmosphere. It was launched from New Zealand then, for 40 days, circumnavigated the globe five times at a latitude 40 to 50 degrees South. Attached to the telescope were four 'DRS' (Data Recovery System) capsules containing 5 TB solid state data storage, plus a GNSS receiver, Iridium transmitter, and parachute. Data from the telescope were copied to these, and two were dropped over Argentina. They drifted 61 km horizontally while they descended 32 km, but we predicted their descent vectors within 2.4 km: in this location, the discrepancy appears irreducible below 2 km because of high speed, gusty winds and local topography. The capsules then reported their own locations to within a few metres. We recovered the capsules and successfully retrieved all of superBIT's data - despite the telescope itself being later destroyed on landing. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 12 pages

Journal ref: Aerospace 2023, 10, 960

Beyond the Ultra-deep Frontier Fields And Legacy Observations (BUFFALO): a high-resolution strong + weak-lensing view of Abell 370

Authors: Anna Niemiec, Mathilde Jauzac, Dominique Eckert, David Lagattuta, Keren Sharon, Anton M. Koekemoer, Keiichi Umetsu, Ana Acebron, Jose M. Diego, David Harvey, Eric Jullo, Vasily Kokorev, Marceau Limousin, Guillaume Mahler, Priyamvada Natarajan, Mario Nonino, Juan D. Remolina, Charles Steinhardt, Sut-Ieng Tam, Adi Zitrin

Abstract: The HST treasury program BUFFALO provides extended wide-field imaging of the six Hubble Frontier Fields galaxy clusters. Here we present the combined strong and weak-lensing analysis of Abell 370, a massive cluster at z=0.375. From the reconstructed total projected mass distribution in the 6arcmin x 6arcmin BUFFALO field-of-view, we obtain the distribution of massive substructures outside the clus… ▽ More The HST treasury program BUFFALO provides extended wide-field imaging of the six Hubble Frontier Fields galaxy clusters. Here we present the combined strong and weak-lensing analysis of Abell 370, a massive cluster at z=0.375. From the reconstructed total projected mass distribution in the 6arcmin x 6arcmin BUFFALO field-of-view, we obtain the distribution of massive substructures outside the cluster core and report the presence of a total of seven candidates, each with mass $\sim 5 \times 10^{13}M_{\odot}$. Combining the total mass distribution derived from lensing with multi-wavelength data, we evaluate the physical significance of each candidate substructure, and conclude that 5 out of the 7 substructure candidates seem reliable, and that the mass distribution in Abell 370 is extended along the North-West and South-East directions. While this finding is in general agreement with previous studies, our detailed spatial reconstruction provides new insights into the complex mass distribution at large cluster-centric radius. We explore the impact of the extended mass reconstruction on the model of the cluster core and in particular, we attempt to physically explain the presence of an important external shear component, necessary to obtain a low root-mean-square separation between the model-predicted and observed positions of the multiple images in the cluster core. The substructures can only account for up to half the amplitude of the external shear, suggesting that more effort is needed to fully replace it by more physically motivated mass components. We provide public access to all the lensing data used as well as the different lens models. △ Less

Submitted 7 July, 2023; originally announced July 2023.

Comments: 29 pages, 17 figures, 3 tables

Federated Neural Radiance Fields

Authors: Lachlan Holden, Feras Dayoub, David Harvey, Tat-Jun Chin

Abstract: The ability of neural radiance fields or NeRFs to conduct accurate 3D modelling has motivated application of the technique to scene representation. Previous approaches have mainly followed a centralised learning paradigm, which assumes that all training images are available on one compute node for training. In this paper, we consider training NeRFs in a federated manner, whereby multiple compute n… ▽ More The ability of neural radiance fields or NeRFs to conduct accurate 3D modelling has motivated application of the technique to scene representation. Previous approaches have mainly followed a centralised learning paradigm, which assumes that all training images are available on one compute node for training. In this paper, we consider training NeRFs in a federated manner, whereby multiple compute nodes, each having acquired a distinct set of observations of the overall scene, learn a common NeRF in parallel. This supports the scenario of cooperatively modelling a scene using multiple agents. Our contribution is the first federated learning algorithm for NeRF, which splits the training effort across multiple compute nodes and obviates the need to pool the images at a central node. A technique based on low-rank decomposition of NeRF layers is introduced to reduce bandwidth consumption to transmit the model parameters for aggregation. Transferring compressed models instead of the raw data also contributes to the privacy of the data collecting agents. △ Less

Submitted 1 May, 2023; originally announced May 2023.

Comments: 10 pages, 7 figures

Hydrodynamic Collapse of the Leidenfrost Vapor Layer

Authors: Dana Harvey, Justin C. Burton

Abstract: During the Leidenfrost effect, a stable vapor film can separate a hot solid from an evaporating liquid. Eventually, after formation and upon cooling, the vapor layer cannot be sustained and undergoes a violent collapse evidenced by explosive boiling. Computationally, modeling this instability involves an interplay between hydrodynamics, thermodynamics, rapid evaporation, and length-scales from… ▽ More During the Leidenfrost effect, a stable vapor film can separate a hot solid from an evaporating liquid. Eventually, after formation and upon cooling, the vapor layer cannot be sustained and undergoes a violent collapse evidenced by explosive boiling. Computationally, modeling this instability involves an interplay between hydrodynamics, thermodynamics, rapid evaporation, and length-scales from $μ$m to cm. Selective assumptions, made to reduce computational costs, have limited most previous studies to steady-state investigations. Here, we combine two-phase laminar flow, heat transfer, and evaporation in a finite-element simulation to examine the failure of Leidenfrost vapor layers during cooling. During periods of quiescence, the geometry of the vapor layer agrees well with steady-state lubrication theory. In the simulations, we report the local temperature of the solid at failure, $T_-$, which provides a lower bound for recent experimental work using the same geometric and material conditions. Surprisingly, we find that inertial forces, which are typically ignored in theoretical treatments of the vapor layer, are responsible for initiating the instability leading to failure. △ Less

Submitted 25 January, 2023; originally announced January 2023.

Comments: 16 pages, 14 figures, 3 supplemental movies

Strong gravitational lensing's 'external shear' is not shear

Authors: Amy Etherington, James W. Nightingale, Richard Massey, Sut-Ieng Tam, XiaoYue Cao, Anna Niemiec, Qiuhan He, Andrew Robertson, Ran Li, Aristeidis Amvrosiadis, Shaun Cole, Jose M. Diego, Carlos S. Frenk, Brenda L. Frye, David Harvey, Mathilde Jauzac, Anton M. Koekemoer, David J. Lagattuta, Marceau Limousin, Guillaume Mahler, Ellen Sirks, Charles L. Steinhardt

Abstract: The distribution of mass in galaxy-scale strong gravitational lenses is often modelled as an elliptical power law plus 'external shear', which notionally accounts for neighbouring galaxies and cosmic shear. We show that it does not. Except in a handful of rare systems, the best-fit values of external shear do not correlate with independent measurements of shear: from weak lensing in 45 Hubble Spac… ▽ More The distribution of mass in galaxy-scale strong gravitational lenses is often modelled as an elliptical power law plus 'external shear', which notionally accounts for neighbouring galaxies and cosmic shear. We show that it does not. Except in a handful of rare systems, the best-fit values of external shear do not correlate with independent measurements of shear: from weak lensing in 45 Hubble Space Telescope images, or in 50 mock images of lenses with complex distributions of mass. Instead, the best-fit shear is aligned with the major or minor axis of 88% of lens galaxies; and the amplitude of the external shear increases if that galaxy is disky. We conclude that 'external shear' attached to a power law model is not physically meaningful, but a fudge to compensate for lack of model complexity. Since it biases other model parameters that are interpreted as physically meaningful in several science analyses (e.g. measuring galaxy evolution, dark matter physics or cosmological parameters), we recommend that future studies of galaxy-scale strong lensing should employ more flexible mass models. △ Less

Submitted 17 November, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Comments: 13 pages, 11 figures, submitted to MNRAS and Review Response Submitted

AEROS: Oceanographic Hyperspectral Imaging and Argos-Tracking CubeSat

Authors: Sophie Prendergast, Cadence Payne, Miles Lifson, Christian Haughwout, Marcos Tieppo, Paulo Figueiredo, André Guerra, Alexander Costa, Helder Magalhães, Tiago Hormigo, Francisco Câmara, Carlos Mano, Pedro Pinheiro, Alvin D. Harvey, Bruno Macena, Luis F. Azevedo, Miguel Martin, Tiago Miranda, Eduardo Pereira, João Faria, Inês Castelão, Catarina Cecilio, Emanuel Castanho, Kerri Cahoy, Manuel Coutinho , et al. (2 additional authors not shown)

Abstract: AEROS is a 3U CubeSat pathfinder toward a future ocean-observing constellation, targeting the Portuguese Atlantic region. AEROS features a miniaturized, high-resolution Hyperspectral Imager (HSI), a 5MP RGB camera, and a Software Defined Radio (SDR). The sensor generated data will be processed and aggregated for end-users in a new web-based Data Analysis Center (DAC). The HSI has 150 spectrally co… ▽ More AEROS is a 3U CubeSat pathfinder toward a future ocean-observing constellation, targeting the Portuguese Atlantic region. AEROS features a miniaturized, high-resolution Hyperspectral Imager (HSI), a 5MP RGB camera, and a Software Defined Radio (SDR). The sensor generated data will be processed and aggregated for end-users in a new web-based Data Analysis Center (DAC). The HSI has 150 spectrally contiguous bands covering visible to near-infrared with 10 nm bandwidth. The HSI collects ocean color data to support studies of oceanographic characteristics known to influence the spatio-temporal distribution and movement behavior of marine organisms. Usage of an SDR expands AEROS's operational and communication range and allows for remote reconfiguration. The SDR receives, demodulates, and retransmits short duration messages, from sources including tagged marine organisms, autonomous vehicles, subsurface floats, and buoys. The future DAC will collect, store, process, and analyze acquired data, taking advantage of its ability to disseminate data across the stakeholders and the scientific network. Correlation of animal-borne Argos platform locations and oceanographic data will advance fisheries management, ecosystem-based management, monitoring of marine protected areas, and bio-oceanographic research in the face of a rapidly changing environment. For example, correlation of oceanographic data collected by the HSI, geolocated with supplementary images from the RGB camera and fish locations, will provide researchers with near real-time estimates of essential oceanographic variables within areas selected by species of interest. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 13 pages, 16 figures, Manuscript presented at the 73rd International Astronautical Congress, IAC 2022, Paris, France, 18 - 22 September 2022

Graphs of Linear Growth have Bounded Treewidth

Authors: Rutger Campbell, Marc Distel, J. Pascal Gollin, Daniel J. Harvey, Kevin Hendrey, Robert Hickingbotham, Bojan Mohar, David R. Wood

Abstract: A graph class $\mathcal{G}$ has linear growth if, for each graph $G \in \mathcal{G}$ and every positive integer $r$, every subgraph of $G$ with radius at most $r$ contains $O(r)$ vertices. In this paper, we show that every graph class with linear growth has bounded treewidth. A graph class $\mathcal{G}$ has linear growth if, for each graph $G \in \mathcal{G}$ and every positive integer $r$, every subgraph of $G$ with radius at most $r$ contains $O(r)$ vertices. In this paper, we show that every graph class with linear growth has bounded treewidth. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Counting points on smooth plane quartics

Authors: Edgar Costa, David Harvey, Andrew V. Sutherland

Abstract: We present efficient algorithms for counting points on a smooth plane quartic curve $X$ modulo a prime $p$. We address both the case where $X$ is defined over $\mathbb F_p$ and the case where $X$ is defined over $\mathbb Q$ and $p$ is a prime of good reduction. We consider two approaches for computing $\#X(\mathbb F_p)$, one which runs in $O(p\log p\log\log p)$ time using $O(\log p)$ space and one… ▽ More We present efficient algorithms for counting points on a smooth plane quartic curve $X$ modulo a prime $p$. We address both the case where $X$ is defined over $\mathbb F_p$ and the case where $X$ is defined over $\mathbb Q$ and $p$ is a prime of good reduction. We consider two approaches for computing $\#X(\mathbb F_p)$, one which runs in $O(p\log p\log\log p)$ time using $O(\log p)$ space and one which runs in $O(p^{1/2}\log^2\!p)$ time using $O(p^{1/2}\log p)$ space. Both approaches yield algorithms that are faster in practice than existing methods. We also present average polynomial-time algorithms for $X/\mathbb Q$ that compute $\#X(\mathbb F_p)$ for good primes $p\le N$ in $O(N\log^3\! N)$ time using $O(N)$ space. These are the first practical implementations of average polynomial-time algorithms for curves that are not cyclic covers of $\mathbb P^1$, which in combination with previous results addresses all curves of genus $g\le 3$. Our algorithms also compute Cartier-Manin/Hasse-Witt matrices that may be of independent interest. △ Less

Submitted 21 August, 2022; originally announced August 2022.

Comments: 32 pages

MSC Class: 11G40 (Primary); 14G10; 14H25 11Y16 (Secondary)

The lifetime of charged dust in the atmosphere

Authors: Josh Méndez Harper, Dana Harvey, Tianshu Huang, Jack McGrath III, David Meer, Justin C. Burton

Abstract: Windblown dust plays critical roles in numerous geophysical and biological systems, yet current models fail to explain the transport of coarse-mode particles (>5 $μ$m) to great distances from their sources. For particles larger than a few microns, electrostatic effects have been invoked to account for longer-than-predicted atmospheric residence times. Although much effort has focused on elucidatin… ▽ More Windblown dust plays critical roles in numerous geophysical and biological systems, yet current models fail to explain the transport of coarse-mode particles (>5 $μ$m) to great distances from their sources. For particles larger than a few microns, electrostatic effects have been invoked to account for longer-than-predicted atmospheric residence times. Although much effort has focused on elucidating the charging processes, comparatively little effort has been expended understanding the stability of charge on particles once electrified. Overall, electrostatic-driven transport requires that charge remain present on particles for days to weeks. Here, we present a set of experiments designed to explore the longevity of electrostatic charge on levitated airborne particles after a single charging event. Using an acoustic levitator, we measured the charge on particles of different material compositions suspended in atmospheric conditions for long periods of time. In dry environments, the total charge on particles decayed in over 1 week. The decay timescale decreased to days in humid environments. These results were independent of particle material and charge polarity. However, exposure to UV radiation could both increase or decrease the decay time depending on polarity. Our work suggests that the rate of charge decay on airborne particles is solely determined by ion capture from the air. Furthermore, using a one-dimensional sedimentation model, we predict that atmospheric dust of order 10 $μ$m will experience the largest change in residence time due to electrostatic forces. △ Less

Submitted 21 June, 2022; v1 submitted 14 June, 2022; originally announced June 2022.

Comments: Main Paper - 8 pages, 5 figures LevChargeSupplementalArxiv.pdf - is the supplemental information which is an additional 3 pages with 3 figures

Constraints on dark matter self-interaction from the internal density profiles of X-COP galaxy clusters

Authors: D. Eckert, S. Ettori, A. Robertson, R. Massey, E. Pointecouteau, D. Harvey, I. G. McCarthy

Abstract: TThe fundamental properties of the postulated dark matter (DM) affect the internal structure of gravitationally-bound structures. In the cold dark matter paradigm, DM particles interact only via gravity. Their distribution is well represented by an Einasto profile with shape parameter $α\approx0.18$, in the smallest dwarf galaxies or the most massive galaxy clusters alike. Conversely, if dark matt… ▽ More TThe fundamental properties of the postulated dark matter (DM) affect the internal structure of gravitationally-bound structures. In the cold dark matter paradigm, DM particles interact only via gravity. Their distribution is well represented by an Einasto profile with shape parameter $α\approx0.18$, in the smallest dwarf galaxies or the most massive galaxy clusters alike. Conversely, if dark matter particles self-interact via additional forces, we expect the mass density profiles of DM halos to flatten in their central regions, thereby increasing the Einasto shape parameter. We measure the structural properties of the 12 massive X-COP galaxy clusters from observations of their hot gaseous atmosphere using the X-ray observatory XMM-Newton, and of the Sunyaev-Zeldovich effect using the Planck all-sky survey. After removing morphologically disturbed systems, we measure Einasto shape parameters with mean $\langleα\rangle = 0.19 \pm 0.03$ and intrinsic scatter $σ_α=0.06$, in close agreement with the prediction of the cold dark matter paradigm. We use cosmological hydrodynamical simulations of cluster formation with self-interacting DM (BAHAMAS-SIDM) to determine how the Einasto shape parameter depends on the self-interaction cross section. We use the fitted relation to turn our measurements of $α$ into constraints on the self-interaction cross section, which imply $σ/m < 0.19$ cm$^2$/g (95% confidence level) at collision velocity $v_\mathrm{DM-DM}\sim1,000$ km/s. This is lower than the interaction cross-section required for dark matter self-interactions to solve the core-cusp problem in dwarf spheroidal galaxies, unless the cross section is a strong function of velocity. △ Less

Submitted 2 May, 2022; originally announced May 2022.

Comments: Re-submitted to A&A after revision

Journal ref: A&A 666, A41 (2022)

A deterministic algorithm for finding $r$-power divisors

Authors: David Harvey, Markus Hittmeir

Abstract: Building on work of Boneh, Durfee and Howgrave-Graham, we present a deterministic algorithm that provably finds all integers $p$ such that $p^r \mathrel| N$ in time $O(N^{1/4r+ε})$ for any $ε> 0$. For example, the algorithm can be used to test squarefreeness of $N$ in time $O(N^{1/8+ε})$; previously, the best rigorous bound for this problem was $O(N^{1/6+ε})$, achieved via the Pollard--Strassen me… ▽ More Building on work of Boneh, Durfee and Howgrave-Graham, we present a deterministic algorithm that provably finds all integers $p$ such that $p^r \mathrel| N$ in time $O(N^{1/4r+ε})$ for any $ε> 0$. For example, the algorithm can be used to test squarefreeness of $N$ in time $O(N^{1/8+ε})$; previously, the best rigorous bound for this problem was $O(N^{1/6+ε})$, achieved via the Pollard--Strassen method. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 16 pages

MSC Class: 11Y05

PyAutoLens: Open-Source Strong Gravitational Lensing

Authors: James. W. Nightingale, Richard G. Hayes, Ashley Kelly, Aristeidis Amvrosiadis, Amy Etherington, Qiuhan He, Nan Li, XiaoYue Cao, Jonathan Frawley, Shaun Cole, Andrea Enia, Carlos S. Frenk, David R. Harvey, Ran Li, Richard J. Massey, Mattia Negrello, Andrew Robertson

Abstract: Strong gravitational lensing, which can make a background source galaxy appears multiple times due to its light rays being deflected by the mass of one or more foreground lens galaxies, provides astronomers with a powerful tool to study dark matter, cosmology and the most distant Universe. PyAutoLens is an open-source Python 3.6+ package for strong gravitational lensing, with core features includi… ▽ More Strong gravitational lensing, which can make a background source galaxy appears multiple times due to its light rays being deflected by the mass of one or more foreground lens galaxies, provides astronomers with a powerful tool to study dark matter, cosmology and the most distant Universe. PyAutoLens is an open-source Python 3.6+ package for strong gravitational lensing, with core features including fully automated strong lens modeling of galaxies and galaxy clusters, support for direct imaging and interferometer datasets and comprehensive tools for simulating samples of strong lenses. The API allows users to perform ray-tracing by using analytic light and mass profiles to build strong lens systems. Accompanying PyAutoLens is the autolens workspace (see https://github.com/Jammy2211/autolens_workspace), which includes example scripts, lens datasets and the HowToLens lectures in Jupyter notebook format which introduce non experts to strong lensing using PyAutoLens. Readers can try PyAutoLens right now by going to the introduction Jupyter notebook on Binder (see https://mybinder.org/v2/gh/Jammy2211/autolens_workspace/master) or checkout the readthedocs (see https://pyautolens.readthedocs.io/en/latest/) for a complete overview of PyAutoLens's features. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: 2 pages, 1 figure

Journal ref: Journal of Open Source Software, 6(58), 2525 (2021)

A log-log speedup for exponent one-fifth deterministic integer factorisation

Authors: David Harvey, Markus Hittmeir

Abstract: Building on techniques recently introduced by the second author, and further developed by the first author, we show that a positive integer $N$ may be rigorously and deterministically factored into primes in at most \[ O\left( \frac{N^{1/5} \log^{16/5} N}{(\log\log N)^{3/5}}\right) \] bit operations. This improves on the previous best known result by a factor of $(\log \log N)^{3/5}$. Building on techniques recently introduced by the second author, and further developed by the first author, we show that a positive integer $N$ may be rigorously and deterministically factored into primes in at most \[ O\left( \frac{N^{1/5} \log^{16/5} N}{(\log\log N)^{3/5}}\right) \] bit operations. This improves on the previous best known result by a factor of $(\log \log N)^{3/5}$. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Comments: 13 pages

MSC Class: 11Y05

The impact of Self-Interacting Dark Matter on the Intrinsic Alignments of Galaxies

Authors: David Harvey, Nora Elisa Chisari, Andrew Robertson, Ian G. McCarthy

Abstract: The formation and evolution of galaxies is known to be sensitive to tidal processes leading to intrinsic correlations between their shapes and orientations. Such correlations can be measured to high significance today, suggesting that cosmological information can be extracted from them. Among the most pressing questions in particle physics and cosmology is the nature of dark matter. If dark matter… ▽ More The formation and evolution of galaxies is known to be sensitive to tidal processes leading to intrinsic correlations between their shapes and orientations. Such correlations can be measured to high significance today, suggesting that cosmological information can be extracted from them. Among the most pressing questions in particle physics and cosmology is the nature of dark matter. If dark matter is self-interacting, it can leave an imprint on galaxy shapes. In this work, we investigate whether self-interactions can produce a long-lasting imprint on intrinsic galaxy shape correlations. We investigate this observable at low redshift ($z<0.4$) using a state-of-the-art suite of cosmological hydro-dynamical simulations where the dark matter model is varied. We find that dark matter self-interactions induce a mass dependent suppression in the intrinsic alignment signal by up to 50\% out to ten's of mega-parsecs, showing that self-interactions can impact structure outside the very core of clusters. We find evidence that self-interactions have a scale-dependent impact on the intrinsic alignment signal that is sufficiently different from signatures introduced by differing baryonic physics prescriptions, suggesting that it is detectable with up-coming all-sky surveys. △ Less

Submitted 25 June, 2021; v1 submitted 5 April, 2021; originally announced April 2021.

Comments: Accepted MNRAS

Intersecting longest paths in chordal graphs

Authors: Daniel J. Harvey, Michael S. Payne

Abstract: We consider the size of the smallest set of vertices required to intersect every longest path in a chordal graph. Such sets are known as longest path transversals. We show that if $ω(G)$ is the clique number of a chordal graph $G$, then there is a transversal of order at most $4\lceil\frac{ω(G)}{5}\rceil$. We also consider the analogous question for longest cycles, and show that if $G$ is a 2-conn… ▽ More We consider the size of the smallest set of vertices required to intersect every longest path in a chordal graph. Such sets are known as longest path transversals. We show that if $ω(G)$ is the clique number of a chordal graph $G$, then there is a transversal of order at most $4\lceil\frac{ω(G)}{5}\rceil$. We also consider the analogous question for longest cycles, and show that if $G$ is a 2-connected chordal graph then there is a transversal intersecting all longest cycles of order at most $2\lceil\frac{ω(G)}{3}\rceil$. △ Less

Submitted 13 December, 2020; originally announced December 2020.

Comments: 11 pages

MSC Class: 05C38

A 4% measurement of $H_0$ using the cumulative distribution of strong-lensing time delays in doubly-imaged quasars

Authors: David Harvey

Abstract: In the advent of large scale surveys, individually modelling strong-gravitational lenses and their counterpart time-delays in order to precisely measure $H_0$ will become computationally expensive, and highly complex. A complimentary approach is to study the cumulative distribution function (CDF) of time-delays where the global population of lenses is modelled along with $H_0$. In this paper we us… ▽ More In the advent of large scale surveys, individually modelling strong-gravitational lenses and their counterpart time-delays in order to precisely measure $H_0$ will become computationally expensive, and highly complex. A complimentary approach is to study the cumulative distribution function (CDF) of time-delays where the global population of lenses is modelled along with $H_0$. In this paper we use a suite of hydro-dynamical simulations to estimate the CDF of time-delays from doubly-imaged quasars for a realistic distribution of lenses. We find that the CDFs exhibit large amounts of halo-halo variance, regulated by the density profile inner slope and the total mass within $5$kpc. With the objective of fitting to data, we compress the CDFs using Principal Component Analysis and fit a Gaussian Processes Regressor consisting of three physical features: the redshift of the lens, $z_{\rm L}$; the power law index of the halo, $α$, and the mass within $5$kpc, plus four cosmological features. Assuming a flat Universe, we fit our model to 27 doubly-imaged quasars finding $H_0=71^{+2}_{-3}$km/s/Mpc, $z_{\rm L} = 0.36_{-0.09}^{+0.2} $, $α=-1.8_{-0.1}^{+0.1}$, $\log(M(<5$kpc$)/M_\odot)=11.1_{-0.1}^{+0.1} $, $Ω_{\rm M} = 0.3_{-0.04}^{+0.04} $ and $Ω_{\rm Λ}=0.7_{-0.04}^{+0.04}$. We compare our estimates of $z_{\rm L}$ and $\log(M(<5$kpc$)/M_\odot)$ to the data and find that within the sensitivity of the data, they are not systematically biased. We generate mock CDFs and find with that the Vera Rubin Observatory (VRO) could measure $σ/H_0$ to $<3\%$, limited by the precision of the model. If we are to exploit fully VRO, we require simulations that sample a larger proportion of the lens population, with a variety of feedback models, exploring all possible systematics. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: https://www.youtube.com/watch?v=8PFgmXI8g0o ---------------------------- https://www.universiteitleiden.nl/en/news/2020/11/universe-expansion-speed-measured-in-bending-starlight

Reconciling galaxy cluster shapes, measured by theorists vs observers

Authors: David Harvey, Andrew Robertson, Sut-Ieng Tam, Mathilde Jauzac, Richard Massey, Jason Rhodes, Ian G. McCarthy

Abstract: If properly calibrated, the shapes of galaxy clusters can be used to investigate many physical processes: from feedback and quenching of star formation, to the nature of dark matter. Theorists frequently measure shapes using moments of inertia of simulated particles'. We instead create mock (optical, X-ray, strong- and weak-lensing) observations of the twenty-two most massive (… ▽ More If properly calibrated, the shapes of galaxy clusters can be used to investigate many physical processes: from feedback and quenching of star formation, to the nature of dark matter. Theorists frequently measure shapes using moments of inertia of simulated particles'. We instead create mock (optical, X-ray, strong- and weak-lensing) observations of the twenty-two most massive ($\sim10^{14.7}\,M_\odot$) relaxed clusters in the BAHAMAS simulations. We find that observable measures of shape are rounder. Even when moments of inertia are projected into 2D and evaluated at matched radius, they overestimate ellipticity by 56\% (compared to observable strong lensing) and 430\% (compared to observable weak lensing). Therefore, we propose matchable quantities and test them using observations of eight relaxed clusters from the {\emph Hubble Space Telescope} and {\emph Chandra X-Ray Observatory}. We also release our HST data reduction and lensing analysis software to the community. In real clusters, the ellipticity and orientation angle at all radii are strongly correlated. In simulated clusters, the ellipticity of inner ($<r_{\mathrm{vir}}/20$) regions becomes decoupled: for example with greater misalignment of the central cluster galaxy. This may indicate overly efficient implementation of feedback from active galactic nuclei. Future exploitation of cluster shapes as a function of radii will require better understanding of core baryonic processes. Exploitation of shapes on any scale will require calibration on simulations extended all the way to mock observations. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: Accepted to MNRAS

The Minimum Leidenfrost Temperature on Smooth Surfaces

Authors: Dana Harvey, Joshua Mendez Harper, Justin C. Burton

Abstract: During the Leidenfrost effect, a thin insulating vapor layer separates an evaporating liquid from a hot solid. Here we demonstrate that Leidenfrost vapor layers can be sustained at much lower temperatures than those required for formation. Using a high-speed electrical technique to measure the thickness of water vapor layers over smooth, metallic surfaces, we find that the explosive failure point… ▽ More During the Leidenfrost effect, a thin insulating vapor layer separates an evaporating liquid from a hot solid. Here we demonstrate that Leidenfrost vapor layers can be sustained at much lower temperatures than those required for formation. Using a high-speed electrical technique to measure the thickness of water vapor layers over smooth, metallic surfaces, we find that the explosive failure point is nearly independent of material and fluid properties, suggesting a purely hydrodynamic mechanism determines this threshold. For water vapor layers of several millimeters in size, the minimum temperature for stability is $\approx 140^{\circ}$C, corresponding to an average vapor layer thickness of 10-20$μ$m. △ Less

Submitted 4 August, 2021; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 6 pages, 4 figures, Supplemental information attached. Supplemental videos available upon request

Journal ref: Phys. Rev. Lett. 127, 104501 (2021)

An exponent one-fifth algorithm for deterministic integer factorisation

Authors: David Harvey

Abstract: Hittmeir recently presented a deterministic algorithm that provably computes the prime factorisation of a positive integer $N$ in $N^{2/9+o(1)}$ bit operations. Prior to this breakthrough, the best known complexity bound for this problem was $N^{1/4+o(1)}$, a result going back to the 1970s. In this paper we push Hittmeir's techniques further, obtaining a rigorous, deterministic factoring algorithm… ▽ More Hittmeir recently presented a deterministic algorithm that provably computes the prime factorisation of a positive integer $N$ in $N^{2/9+o(1)}$ bit operations. Prior to this breakthrough, the best known complexity bound for this problem was $N^{1/4+o(1)}$, a result going back to the 1970s. In this paper we push Hittmeir's techniques further, obtaining a rigorous, deterministic factoring algorithm with complexity $N^{1/5+o(1)}$. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 14 pages

MSC Class: 11Y05

The distribution of dark matter and gas spanning six megaparsecs around the post-merger galaxy cluster MS0451-03

Authors: Sut-Ieng Tam, Mathilde Jauzac, Richard Massey, David Harvey, Dominique Eckert, Harald Ebeling, Richard S. Ellis, Vittorio Ghirardini, Baptiste Klein, Jean-Paul Kneib, David Lagattuta, Priyamvada Natarajan, Andrew Robertson, Graham P. Smith

Abstract: Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures.… ▽ More Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions. △ Less

Submitted 28 July, 2020; v1 submitted 17 June, 2020; originally announced June 2020.

Comments: 21 pages, 13 figures, updated to match MNRAS version

From dwarf galaxies to galaxy clusters: Self-Interacting Dark Matter over 7 orders of magnitude in halo mass

Authors: Kyrylo Bondarenko, Anastasia Sokolenko, Alexey Boyarsky, Andrew Robertson, David Harvey, Yves Revaz

Abstract: In this paper we study the density profiles of self-interacting dark matter (SIDM) haloes spanning the full observable mass range, from dwarf galaxies to galaxy clusters. Using realistic simulations that model the baryonic physics relevant for galaxy formation, we compare the density profiles of haloes simulated with either SIDM or cold and collisionless dark matter (CDM) to those inferred from ob… ▽ More In this paper we study the density profiles of self-interacting dark matter (SIDM) haloes spanning the full observable mass range, from dwarf galaxies to galaxy clusters. Using realistic simulations that model the baryonic physics relevant for galaxy formation, we compare the density profiles of haloes simulated with either SIDM or cold and collisionless dark matter (CDM) to those inferred from observations of stellar velocity dispersion, gas rotation curves, weak and strong gravitational lensing, and/or X-ray maps. We make our comparison in terms of the maximal surface density of haloes, circumventing the need for semi-analytic or parametric models for dark matter density profiles. We find that the maximal surface density as a function of halo mass is well reproduced by CDM simulations that include baryons, while for SIDM with a velocity-independent cross-section of 1 cm$^2/$g, the simulated galaxy clusters have mean maximal surface densities that are below those of observed systems by an amount greater than the standard deviation of the observed maximal surface density at fixed mass. For less massive systems both CDM and SIDM agree with the observation equally well. △ Less

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

Comments: Journal version; minor changes

H0LiCOW XI. A weak lensing measurement of the external convergence in the field of the lensed quasar B1608+656 using HST and Subaru deep imaging

Authors: O. Tihhonova, F. Courbin, D. Harvey, S. Hilbert, A. Peel, C. E. Rusu, C. D. Fassnacht, V. Bonvin, P. J. Marshall, G. Meylan, D. Sluse, S. H. Suyu, T. Treu, K. C. Wong

Abstract: We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse… ▽ More We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $κ_{\rm ext}$ at the position of B1608+656 is $κ_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: Accepted MNRAS

The BUFFALO HST Survey

Authors: Charles L. Steinhardt, Mathilde Jauzac, Ana Acebron, Hakim Atek, Peter Capak, Iary Davidzon, Dominique Eckert, David Harvey, Anton M. Koekemoer, Claudia D. P. Lagos, Guillaume Mahler, Mireia Montes, Anna Niemiec, Mario Nonino, P. A. Oesch, Johan Richard, Steven A. Rodney, Matthieu Schaller, Keren Sharon, Louis-Gregory Strolger, Joseph Allingham, Adam Amara, Yannick Bah'e, Celine Boehm, Sownak Bose , et al. (70 additional authors not shown)

Abstract: The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has no… ▽ More The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters. △ Less

Submitted 13 February, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: Accepted ApJS; MAST archive will be live concurrent with publication

Distributed Possibilistic Learning in Multi-Agent Systems

Authors: Jonathan Lawry, Michael Crosscombe, David Harvey

Abstract: Possibility theory is proposed as an uncertainty representation framework for distributed learning in multi-agent systems and robot swarms. In particular, we investigate its application to the best-of-n problem where the aim is for a population of agents to identify the highest quality out of n options through local interactions between individuals and limited direct feedback from the environment.… ▽ More Possibility theory is proposed as an uncertainty representation framework for distributed learning in multi-agent systems and robot swarms. In particular, we investigate its application to the best-of-n problem where the aim is for a population of agents to identify the highest quality out of n options through local interactions between individuals and limited direct feedback from the environment. In this context we claim that possibility theory provides efficient mechanisms by which an agent can learn about the state of the world, and which can allow them to handle inconsistencies between what they and others believe by varying the level of imprecision of their own beliefs. We introduce a discrete time model of a population of agents applying possibility theory to the best-of-n problem. Simulation experiments are then used to investigate the accuracy of possibility theory in this context as well as its robustness to noise under varying amounts of direct evidence. Finally, we compare possibility theory in this context with a similar probabilistic approach. △ Less

Submitted 20 January, 2020; originally announced January 2020.

Comments: 3rd International Symposium on Swarm Behavior and Bio-Inspired Robotics (SWARM 2019)

Setting the scene for BUFFALO: A study of the matter distribution in the HFF galaxy cluster MACS J0416.1-2403 and its parallel field

Authors: Elizabeth J. Gonzalez, Martín Chalela, Mathilde Jauzac, Dominique Eckert, Matthieu Schaller, David Harvey, Anna Niemiec, Anton M. Koekemoer, David Barnes, Doug Clowe, Thomas Connor, José M. Diego, Juan D. Remolina Gonzalez, Charles L. Steinhardt

Abstract: In the context of the BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations) survey, we present a new analysis of the merging galaxy cluster MACS\,J0416.1-2403 ($z = 0.397$) and its parallel field using the data collected by the Hubble Frontier Fields (HFF) campaign. In this work, we measure the surface mass density from a weak-lensing analysis, and characterise the overall matter dis… ▽ More In the context of the BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations) survey, we present a new analysis of the merging galaxy cluster MACS\,J0416.1-2403 ($z = 0.397$) and its parallel field using the data collected by the Hubble Frontier Fields (HFF) campaign. In this work, we measure the surface mass density from a weak-lensing analysis, and characterise the overall matter distribution in both the cluster and parallel fields. The surface mass distribution derived for the parallel field shows clumpy overdensities connected by filament-like structures elongated in the direction of the cluster core. We also characterise the X-ray emission of the cluster, and compare it with the lensing mass distribution. We identify five substructures at the $>5σ$ level over the two fields, four of them being in the cluster one. Furthermore, three of them are located close to the edges of the field of view, and border issues can significantly hamper the determination of their physical parameters. Finally, we compare our results with the predicted subhalo distribution of one of the Hydrangea/C-EAGLE simulated cluster. Significant differences are obtained suggesting the simulated cluster is at a more advanced evolutionary state than MACS\,J0416.1-2403. Our results anticipate the upcoming BUFFALO observations that will link the two HFF fields, extending further the \emph{HST} coverage, and thus allowing a better characterisation of the reported substructures. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: 15 pages, 8 figures, 3 tables, submitted to MNRAS

Exploiting flux ratio anomalies to probe warm dark matter in future large scale surveys

Authors: David Harvey, Wessel Valkenburg, Amelie Tamone, Alexey Boyarsky, Frederic Courbin, Mark Lovell

Abstract: Flux ratio anomalies in strong gravitationally lensed quasars constitute a unique way to probe the abundance of non-luminous dark matter haloes, and hence the nature of dark matter. In this paper we identify double imaged quasars as a statistically efficient probe of dark matter, since they are 20 times more abundant than quadruply imaged quasars. Using N-body simulations that include realistic ba… ▽ More Flux ratio anomalies in strong gravitationally lensed quasars constitute a unique way to probe the abundance of non-luminous dark matter haloes, and hence the nature of dark matter. In this paper we identify double imaged quasars as a statistically efficient probe of dark matter, since they are 20 times more abundant than quadruply imaged quasars. Using N-body simulations that include realistic baryonic feedback, we measure the full distribution of flux ratios in doubly imaged quasars for cold (CDM) and warm dark matter (WDM) cosmologies. Through this method, we fold in two key systematics - quasar variability and line-of-sight structures. We find that WDM cosmologies predict a ~6 per cent difference in the cumulative distribution functions of flux ratios relative to CDM, with CDM predicting many more small ratios. Finally, we estimate that ~600 doubly imaged quasars will need to be observed in order to be able to unambiguously discern between CDM and the two WDM models studied here. Such sample sizes will be easily within reach of future large scale surveys such as Euclid. In preparation for this survey data we require discerning the scale of the uncertainties in modelling lens galaxies and their substructure in simulations, plus a strong understanding of the selection function of observed lensed quasars. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: Accepted MNRAS

The shape of relaxed galaxy clusters and the public release of a HST shape measurement code, pyRRG

Authors: David Harvey, Sut-Ieng Tam, Mathilde Jauzac, Richard Massey, Jason Rhodes

Abstract: We study the shape of eight dynamically relaxed galaxy clusters observed with the Hubble Space Telescope and Chandra X-Ray Observatory. Using strong and weak gravitational lensing, the shape of the Brightest Cluster Galaxy and the X-ray isophote we study the ellipticity of the cluster halo at four different radii. We find that the proxies probing the inner regions of the cluster are strongly coupl… ▽ More We study the shape of eight dynamically relaxed galaxy clusters observed with the Hubble Space Telescope and Chandra X-Ray Observatory. Using strong and weak gravitational lensing, the shape of the Brightest Cluster Galaxy and the X-ray isophote we study the ellipticity of the cluster halo at four different radii. We find that the proxies probing the inner regions of the cluster are strongly coupled with the BCG shape correlated with both the shape predicted by strong gravitational lensing and the X-ray isophote. Conversely we find no such correlation between the shape as predicted by the weak lensing and the other three probes suggesting any coupling between the inner and outer regions is weak. We also present in this paper the public release of the HST weak lensing shape measurement code pyRRG, directly available from PyPi (https://pypi.org/project/pyRRG/). This python3.7 code, based on Rhodes et al. (2000) adopts an automated star-galaxy classifier based on a Random Forest and outputs scientifically useful products such as weak lensing catalogues suitable for the mass mapping algorithm Lenstool. △ Less

Submitted 20 November, 2019; v1 submitted 14 November, 2019; originally announced November 2019.

Galaxy structure with strong gravitational lensing: decomposing the internal mass distribution of massive elliptical galaxies

Authors: James. W. Nightingale, Richard J. Massey, David R. Harvey, Andrew P. Cooper, Amy Etherington, Sut-Ieng Tam, Richard G. Hayes

Abstract: We investigate how strong gravitational lensing can test contemporary models of massive elliptical (ME) galaxy formation, by combining a traditional decomposition of their visible stellar distribution with a lensing analysis of their mass distribution. As a proof of concept, we study a sample of three ME lenses, observing that all are composed of two distinct baryonic structures, a `red' central b… ▽ More We investigate how strong gravitational lensing can test contemporary models of massive elliptical (ME) galaxy formation, by combining a traditional decomposition of their visible stellar distribution with a lensing analysis of their mass distribution. As a proof of concept, we study a sample of three ME lenses, observing that all are composed of two distinct baryonic structures, a `red' central bulge surrounded by an extended envelope of stellar material. Whilst these two components look photometrically similar, their distinct lensing effects permit a clean decomposition of their mass structure. This allows us to infer two key pieces of information about each lens galaxy: (i) the stellar mass distribution (without invoking stellar populations models) and (ii) the inner dark matter halo mass. We argue that these two measurements are crucial to testing models of ME formation, as the stellar mass profile provides a diagnostic of baryonic accretion and feedback whilst the dark matter mass places each galaxy in the context of LCDM large scale structure formation. We also detect large rotational offsets between the two stellar components and a lopsidedness in their outer mass distributions, which hold further information on the evolution of each ME. Finally, we discuss how this approach can be extended to galaxies of all Hubble types and what implication our results have for studies of strong gravitational lensing. △ Less

Submitted 12 July, 2019; v1 submitted 23 January, 2019; originally announced January 2019.

Comments: 19 pages, 4 figures. source code: https://github.com/Jammy2211/PyAutoLens

Observable tests of self-interacting dark matter in galaxy clusters: BCG wobbles in a constant density core

Authors: David Harvey, Andrew Robertson, Richard Massey, Ian G. McCarthy

Abstract: Models of Cold Dark Matter predict that the distribution of dark matter in galaxy clusters should be cuspy, centrally concentrated. Constant density cores would be strong evidence for beyond-CDM physics, such as Self-Interacting Dark Matter (SIDM). An observable consequence would be oscillations of the Brightest Cluster Galaxy (BCG) in otherwise relaxed galaxy clusters. Offset BCGs have indeed bee… ▽ More Models of Cold Dark Matter predict that the distribution of dark matter in galaxy clusters should be cuspy, centrally concentrated. Constant density cores would be strong evidence for beyond-CDM physics, such as Self-Interacting Dark Matter (SIDM). An observable consequence would be oscillations of the Brightest Cluster Galaxy (BCG) in otherwise relaxed galaxy clusters. Offset BCGs have indeed been observed - but only interpreted via a simplified, analytic model of oscillations. We compare these observations to the BAHAMAS-SIDM suite of cosmological simulations, which include SIDM and a fully hydrodynamical treatment of star formation and feedback. We predict that the median offset of BCGs increases with the SIDM cross-section, cluster mass and the amount of stellar mass within 10kpc, while CDM exhibits no trend in mass. Interpolating between the simulated cross-sections, we find that the observations (of 10 clusters) is consistent with CDM at the ~1.5$σ$ level, and prefer cross-section $σ$/m < 0.12(0.39)cm$^2$/g at 68% (95%) confidence level. This is on the verge of ruling out velocity-independent dark matter self-interactions as the solution to discrepancies between the predicted and observed behaviour of dwarf galaxies, and will be improved by larger surveys by Euclid or SuperBIT. △ Less

Submitted 13 September, 2019; v1 submitted 17 December, 2018; originally announced December 2018.

The core of the massive cluster merger MACS J0417.5-1154 as seen by VLT/MUSE

Authors: Mathilde Jauzac, Guillaume Mahler, Alastair C. Edge, Keren Sharon, Steven Gillman, Harald Ebeling, David Harvey, Johan Richard, Steven L. Hamer, Michele Fumagalli, A. Mark Swinbank, Jean-Paul Kneib, Richard Massey, Philippe Salome

Abstract: We present a multi-wavelength analysis of the core of the massive galaxy cluster MACS\,J0417.5-1154 ($z = 0.441$; MACS\;J0417). Our analysis takes advantage of VLT/MUSE observations which allow the spectroscopic confirmation of three strongly-lensed systems. One of these, nick-named \emph{The Doughnut}, consists of three complete images of a complex ring galaxy at $z = 0.8718$ and a fourth, partia… ▽ More We present a multi-wavelength analysis of the core of the massive galaxy cluster MACS\,J0417.5-1154 ($z = 0.441$; MACS\;J0417). Our analysis takes advantage of VLT/MUSE observations which allow the spectroscopic confirmation of three strongly-lensed systems. One of these, nick-named \emph{The Doughnut}, consists of three complete images of a complex ring galaxy at $z = 0.8718$ and a fourth, partial and radial image close to the Brightest Cluster Galaxy (BCG) only discernible thanks to its strong [OII] line emission. The best-fit mass model (rms of 0.38\arcsec) yields a two-dimensional enclosed mass of $M({\rm R < 200\,kpc}) = (1.77\pm0.03)\times10^{14}\,\msun$ and almost perfect alignment between the peaks of the BCG light and the dark matter of ($0.5\pm0.5$)\arcsec . Our finding that a significant misalignment results when the radial image of \emph{The Doughnut} is omitted serves as an important caveat for studies of BCG-dark matter offsets in galaxy clusters. Using \emph{Chandra} data to map the intra-cluster gas, we observe an offset between the gas and dark-matter peaks of ($1.7\pm0.5$)\arcsec, and excellent alignment of the X-ray peak with the location of optical emission line associated with the BCG. We interpret all observational evidence in the framework of on-going merger activity, noting specifically that the coincidence between the gas peak and the peak of blue light from the BCG may be evidence of dense, cold gas leading to direct star formation. We use the surface area $σ_μ$ above a given magnification factor $μ$ as a metric to estimate the lensing power of MACS\,J0417. We obtain $σ(μ> 3) = 0.22$\,arcmin$^2$ confirming MACS\,J0417 as an efficient gravitational lens. Finally, we discuss the differences between our mass model and Mahler et al. (2018). △ Less

Submitted 6 December, 2018; v1 submitted 31 October, 2018; originally announced November 2018.

Comments: 15 pages, 6 figures, 8 tables, accepted for publication in MNRAS

Observable tests of self-interacting dark matter in galaxy clusters: cosmological simulations with SIDM and baryons

Authors: Andrew Robertson, David Harvey, Richard Massey, Vincent Eke, Ian G. McCarthy, Mathilde Jauzac, Baojiu Li, Joop Schaye

Abstract: We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations.… ▽ More We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing. △ Less

Submitted 1 June, 2020; v1 submitted 12 October, 2018; originally announced October 2018.

Comments: 19 pages, 13 figures, updated to match MNRAS version

The impact of cored density profiles on the observable quantities of dwarf spheroidal galaxies

Authors: David Harvey, Yves Revaz, Andrew Robertson, Loic Hausammann

Abstract: We modify the chemo-dynamical code GEAR to simulate the impact of self-interacting dark matter on the observable quantities of 19 low mass dwarf galaxies with a variety star forming properties. We employ a relatively high, velocity independent cross-section of $σ/m = 10$cm$^2$/g and extract, in addition to integrated quantities, the total mass density profile, the luminosity profile, the line-of-s… ▽ More We modify the chemo-dynamical code GEAR to simulate the impact of self-interacting dark matter on the observable quantities of 19 low mass dwarf galaxies with a variety star forming properties. We employ a relatively high, velocity independent cross-section of $σ/m = 10$cm$^2$/g and extract, in addition to integrated quantities, the total mass density profile, the luminosity profile, the line-of-sight velocities, the chemical abundance and the star formation history. We find that despite the creation of large cores at the centre of the dark matter haloes, the impact of SIDM on the observable quantities of quenched galaxies is indiscernible, dominated mostly by the stochastic build up of the stellar matter. As such we conclude that it is impossible to make global statements on the density profile of dwarf galaxies from single or small samples. Although based mostly on quenched galaxies, this finding supports other recent work putting into question the reliability of inferred cored density profiles that are derived from observed line-of-sight velocities. △ Less

Submitted 30 August, 2018; originally announced August 2018.

Comments: MNRAS Accepted 2018 August 29. Received 2018 August 13; in original form 2018 May 29

Chromatic roots at 2 and at the Beraha number $B_{10}$

Authors: Daniel J. Harvey, Gordon F. Royle

Abstract: By the construction of suitable graphs and the determination of their chromatic polynomials, we resolve two open questions concerning real chromatic roots. First we exhibit graphs for which the Beraha number $B_{10} = (5 + \sqrt{5})/2$ is a chromatic root. As it was previously known that no other non-integer Beraha number is a chromatic root, this completes the determination of precisely which Ber… ▽ More By the construction of suitable graphs and the determination of their chromatic polynomials, we resolve two open questions concerning real chromatic roots. First we exhibit graphs for which the Beraha number $B_{10} = (5 + \sqrt{5})/2$ is a chromatic root. As it was previously known that no other non-integer Beraha number is a chromatic root, this completes the determination of precisely which Beraha numbers can be chromatic roots. Next we construct an infinite family of $3$-connected graphs such that for any $k \geqslant 1$, there is a member of the family with $q=2$ as a chromatic root of multiplicity at least $k$. The former resolves a question of Salas and Sokal [J. Statist. Pys. 104 (2001) pp. 609--699] and the latter a question of Dong and Koh [J. Graph Theory 70 (2012) pp. 262--283]. △ Less

Submitted 22 July, 2018; originally announced July 2018.

Comments: This first section of this paper was previously uploaded as arxiv:1611.01941

MSC Class: 05C31

Information estimation using nonparametric copulas

Authors: Houman Safaai, Arno Onken, Christopher D. Harvey, Stefano Panzeri

Abstract: Estimation of mutual information between random variables has become crucial in a range of fields, from physics to neuroscience to finance. Estimating information accurately over a wide range of conditions relies on the development of flexible methods to describe statistical dependencies among variables, without imposing potentially invalid assumptions on the data. Such methods are needed in cases… ▽ More Estimation of mutual information between random variables has become crucial in a range of fields, from physics to neuroscience to finance. Estimating information accurately over a wide range of conditions relies on the development of flexible methods to describe statistical dependencies among variables, without imposing potentially invalid assumptions on the data. Such methods are needed in cases that lack prior knowledge of their statistical properties and that have limited sample numbers. Here we propose a powerful and generally applicable information estimator based on non-parametric copulas. This estimator, called the non-parametric copula-based estimator (NPC), is tailored to take into account detailed stochastic relationships in the data independently of the data's marginal distributions. The NPC estimator can be used both for continuous and discrete numerical variables and thus provides a single framework for the mutual information estimation of both continuous and discrete data. By extensive validation on artificial samples drawn from various statistical distributions, we found that the NPC estimator compares well against commonly used alternatives. Unlike methods not based on copulas, it allows an estimation of information that is robust to changes of the details of the marginal distributions. Unlike parametric copula methods, it remains accurate regardless of the precise form of the interactions between the variables. In addition, the NPC estimator had accurate information estimates even at low sample numbers, in comparison to alternative estimators. The NPC estimator therefore provides a good balance between general applicability to arbitrarily shaped statistical dependencies in the data and shows accurate and robust performance when working with small sample sizes. △ Less

Submitted 25 October, 2018; v1 submitted 20 July, 2018; originally announced July 2018.

Comments: 17 pages, 13 figures

Journal ref: Phys. Rev. E 98, 053302 (2018)

Zeta functions of nondegenerate hypersurfaces in toric varieties via controlled reduction in $p$-adic cohomology

Authors: Edgar Costa, David Harvey, Kiran S. Kedlaya

Abstract: We give an interim report on some improvements and generalizations of the Abbott-Kedlaya-Roe method to compute the zeta function of a nondegenerate ample hypersurface in a projectively normal toric variety over $\mathbb{F}_p$ in linear time in $p$. These are illustrated with a number of examples including K3 surfaces, Calabi-Yau threefolds, and a cubic fourfold. The latter example is a non-special… ▽ More We give an interim report on some improvements and generalizations of the Abbott-Kedlaya-Roe method to compute the zeta function of a nondegenerate ample hypersurface in a projectively normal toric variety over $\mathbb{F}_p$ in linear time in $p$. These are illustrated with a number of examples including K3 surfaces, Calabi-Yau threefolds, and a cubic fourfold. The latter example is a non-special cubic fourfold appearing in the Ranestad-Voisin coplanar divisor on moduli space; this verifies that the coplanar divisor is not a Noether-Lefschetz divisor in the sense of Hassett. △ Less

Submitted 12 November, 2018; v1 submitted 1 June, 2018; originally announced June 2018.

Comments: 17 pages, 3 figures; v3 updated discussion regarding cubic fourfolds

MSC Class: 11G25; 14G10; 11M38; 11Y16; 14C25

Journal ref: Open Book Series 2 (2019) 221-238

Faster integer multiplication using short lattice vectors

Authors: David Harvey, Joris van der Hoeven

Abstract: We prove that $n$-bit integers may be multiplied in $O(n \log n \, 4^{\log^* n})$ bit operations. This complexity bound had been achieved previously by several authors, assuming various unproved number-theoretic hypotheses. Our proof is unconditional, and depends in an essential way on Minkowski's theorem concerning lattice vectors in symmetric convex sets. We prove that $n$-bit integers may be multiplied in $O(n \log n \, 4^{\log^* n})$ bit operations. This complexity bound had been achieved previously by several authors, assuming various unproved number-theoretic hypotheses. Our proof is unconditional, and depends in an essential way on Minkowski's theorem concerning lattice vectors in symmetric convex sets. △ Less

Submitted 22 February, 2018; originally announced February 2018.

Comments: 16 pages

MSC Class: 68W30; 68Q17; 68W40 ACM Class: G.1.0; F.2.1

Journal ref: Open Book Series 2 (2019) 293-310

Faster integer and polynomial multiplication using cyclotomic coefficient rings

Authors: David Harvey, Joris van der Hoeven

Abstract: We present an algorithm that computes the product of two n-bit integers in O(n log n (4\sqrt 2)^{log^* n}) bit operations. Previously, the best known bound was O(n log n 6^{log^* n}). We also prove that for a fixed prime p, polynomials in F_p[X] of degree n may be multiplied in O(n log n 4^{log^* n}) bit operations; the previous best bound was O(n log n 8^{log^* n}). We present an algorithm that computes the product of two n-bit integers in O(n log n (4\sqrt 2)^{log^* n}) bit operations. Previously, the best known bound was O(n log n 6^{log^* n}). We also prove that for a fixed prime p, polynomials in F_p[X] of degree n may be multiplied in O(n log n 4^{log^* n}) bit operations; the previous best bound was O(n log n 8^{log^* n}). △ Less

Submitted 11 December, 2017; originally announced December 2017.

Comments: 28 pages

MSC Class: 68W30; 68Q17; 68W40 ACM Class: G.1.0; F.2.1

Quantifying how much sensory information in a neural code is relevant for behavior

Authors: Giuseppe Pica, Eugenio Piasini, Houman Safaai, Caroline A. Runyan, Mathew E. Diamond, Tommaso Fellin, Christoph Kayser, Christopher D. Harvey, Stefano Panzeri

Abstract: Determining how much of the sensory information carried by a neural code contributes to behavioral performance is key to understand sensory function and neural information flow. However, there are as yet no analytical tools to compute this information that lies at the intersection between sensory coding and behavioral readout. Here we develop a novel measure, termed the information-theoretic inter… ▽ More Determining how much of the sensory information carried by a neural code contributes to behavioral performance is key to understand sensory function and neural information flow. However, there are as yet no analytical tools to compute this information that lies at the intersection between sensory coding and behavioral readout. Here we develop a novel measure, termed the information-theoretic intersection information $I_{II}(S;R;C)$, that quantifies how much of the sensory information carried by a neural response R is used for behavior during perceptual discrimination tasks. Building on the Partial Information Decomposition framework, we define $I_{II}(S;R;C)$ as the part of the mutual information between the stimulus S and the response R that also informs the consequent behavioral choice C. We compute $I_{II}(S;R;C)$ in the analysis of two experimental cortical datasets, to show how this measure can be used to compare quantitatively the contributions of spike timing and spike rates to task performance, and to identify brain areas or neural populations that specifically transform sensory information into choice. △ Less

Submitted 6 December, 2017; originally announced December 2017.

Journal ref: Advances in Neural Information Processing Systems 30, 3689--3699, 2017

The shape of galaxy dark matter halos in massive galaxy clusters: Insights from strong gravitational lensing

Authors: M. Jauzac, D. Harvey, R. Massey

Abstract: We assess how much unused strong lensing information is available in the deep \emph{Hubble Space Telescope} imaging and VLT/MUSE spectroscopy of the \emph{Frontier Field} clusters. As a pilot study, we analyse galaxy cluster MACS\,J0416.1-2403 ($z$$=$$0.397$, $M(R<200\,{\rm kpc})$$=$$1.6$$\times$$10^{14}\msun$), which has 141 multiple images with spectroscopic redshifts. We find that many addition… ▽ More We assess how much unused strong lensing information is available in the deep \emph{Hubble Space Telescope} imaging and VLT/MUSE spectroscopy of the \emph{Frontier Field} clusters. As a pilot study, we analyse galaxy cluster MACS\,J0416.1-2403 ($z$$=$$0.397$, $M(R<200\,{\rm kpc})$$=$$1.6$$\times$$10^{14}\msun$), which has 141 multiple images with spectroscopic redshifts. We find that many additional parameters in a cluster mass model can be constrained, and that adding even small amounts of extra freedom to a model can dramatically improve its figures of merit. We use this information to constrain the distribution of dark matter around cluster member galaxies, simultaneously with the cluster's large-scale mass distribution. We find tentative evidence that some galaxies' dark matter has surprisingly similar ellipticity to their stars (unlike in the field, where it is more spherical), but that its orientation is often misaligned. When non-coincident dark matter and baryonic halos are allowed, the model improves by 35\%. This technique may provide a new way to investigate the processes and timescales on which dark matter is stripped from galaxies as they fall into a massive cluster. Our preliminary conclusions will be made more robust by analysing the remaining five \emph{Frontier Field} clusters. △ Less

Submitted 3 April, 2018; v1 submitted 27 November, 2017; originally announced November 2017.

Comments: 6 pages, 3 figures, 2 tables - Accepted for publication in MNRAS

H0LiCOW VIII. A weak lensing measurement of the external convergence in the field of the lensed quasar HE$\,$0435$-$1223

Authors: O. Tihhonova, F. Courbin, D. Harvey, S. Hilbert, C. E. Rusu, C. D. Fassnacht, V. Bonvin, P. J. Marshall, G. Meylan, D. Sluse, S. H. Suyu, T. Treu, K. C. Wong

Abstract: We present a weak gravitational lensing measurement of the external convergence along the line of sight to the quadruply lensed quasar HE$\,$0435$-$1223. Using deep r-band images from Subaru-Suprime-Cam we observe galaxies down to a 3$σ$ limiting magnitude of $\sim 26$ mags resulting in a source galaxy density of 14 galaxies / arcmin$^2$ after redshift-based cuts. Using an inpainting technique and… ▽ More We present a weak gravitational lensing measurement of the external convergence along the line of sight to the quadruply lensed quasar HE$\,$0435$-$1223. Using deep r-band images from Subaru-Suprime-Cam we observe galaxies down to a 3$σ$ limiting magnitude of $\sim 26$ mags resulting in a source galaxy density of 14 galaxies / arcmin$^2$ after redshift-based cuts. Using an inpainting technique and Multi-Scale Entropy filtering algorithm, we find that the region in close proximity to the lens has an estimated external convergence of $κ=-0.012^{+0.020}_{-0.013}$ and is hence marginally under-dense. We also rule out the presence of any halo with a mass greater than $M_{\rm vir}=1.6\times10^{14}h^{-1}M_\odot$ (68$\%$ confidence limit). Our results, consistent with previous studies of this lens, confirm that the intervening mass along the line of sight to HE$\,$0435$-$1223 does not affect significantly the cosmological results inferred from the time delay measurements of that specific object. △ Less

Submitted 23 November, 2017; originally announced November 2017.

Comments: 14 pages, 8 figures, 2 tables

Growing a `Cosmic Beast': Observations and Simulations of MACS J0717.5+3745

Authors: M. Jauzac, D. Eckert, M. Schaller, J. Schwinn, R. Massey, Y. Bahé, C. Baugh, D. Barnes, C. Dalla Vecchia, H. Ebeling, D. Harvey, E. Jullo, S. T. Kay, J. -P. Kneib, M. Limousin, E. Medezinski, P. Natarajan, M. Nonino, A. Robertson, S. I. Tam, K. Umetsu

Abstract: We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected r… ▽ More We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters. △ Less

Submitted 27 August, 2018; v1 submitted 3 November, 2017; originally announced November 2017.

Comments: 18 pages, 6 figures, 8 tables, Accepted for publication in MNRAS

Network Controllability in the IFG Relates to Controlled Language Variability and Susceptibility to TMS

Authors: John D. Medaglia, Denise Y. Harvey, Nicole White, Danielle S. Bassett, Roy H. Hamilton

Abstract: In language production, humans are confronted with considerable word selection demands. Often, we must select a word from among similar, acceptable, and competing alternative words in order to construct a sentence that conveys an intended meaning. In recent years, the left inferior frontal gyrus (LIFG) has been identified as critical to this ability. Despite a recent emphasis on network approaches… ▽ More In language production, humans are confronted with considerable word selection demands. Often, we must select a word from among similar, acceptable, and competing alternative words in order to construct a sentence that conveys an intended meaning. In recent years, the left inferior frontal gyrus (LIFG) has been identified as critical to this ability. Despite a recent emphasis on network approaches to understanding language, how the LIFG interacts with the brain's complex networks to facilitate controlled language performance remains unknown. Here, we take a novel approach to understand word selection as a network control process in the brain. Using an anatomical brain network derived from high-resolution diffusion spectrum imaging (DSI), we computed network controllability underlying the site of transcranial magnetic stimulation in the LIFG between administrations of two word selection tasks. We find that a statistic that quantifies the LIFG's theoretically predicted control of difficult-to-reach states explains vulnerability to TMS in language tasks that vary in response (cognitive control) demands: open-response (word generation) vs. closed-response (number naming) tasks. Moreover, we find that a statistic that quantifies the LIFG's theoretically predicted control of communication across modules in the human connectome explains TMS-induced changes in open-response language task performance only. These findings establish a link between network controllability, cognitive function, and TMS effects. △ Less

Submitted 30 January, 2018; v1 submitted 14 September, 2017; originally announced September 2017.

Solitary Waves in Optical Fibers Governed by Higher Order Dispersion

Authors: Vladimir I. Kruglov, John D. Harvey

Abstract: An exact solitary wave solution is presented for the nonlinear Schrodinger equation governing the propagation of pulses in optical fibers including the effects of second, third and fourth order dispersion. The stability of this soliton-like solution with sech2 shape is proven by the sign-definiteness of the operator and an integral of the Sobolev type. The main criteria governing the existence of… ▽ More An exact solitary wave solution is presented for the nonlinear Schrodinger equation governing the propagation of pulses in optical fibers including the effects of second, third and fourth order dispersion. The stability of this soliton-like solution with sech2 shape is proven by the sign-definiteness of the operator and an integral of the Sobolev type. The main criteria governing the existence of such stable localized pulses propagating in optical fibers are also formulated. A unique feature of these soliton-like optical pulses propagating in a fiber with higher order dispersion is that their parameters satisfy efficient scaling relations. The main soliton solution term given by perturbation theory is also presented when absorption or gain is included in the nonlinear Schrodinger equation. We anticipate that this type of stable localized pulses could find practical applications in communications, slow-light devices and ultrafast lasers. △ Less

Submitted 4 September, 2017; originally announced September 2017.

Comments: 4 pages 3 Figures

Journal ref: Phys. Rev. A 98, 063811 (2018)

Dark matter dynamics in Abell 3827: new data consistent with standard Cold Dark Matter

Authors: Richard Massey, David Harvey, Jori Liesenborgs, Johan Richard, Stuart Stach, Mark Swinbank, Peter Taylor, Liliya Williams, Douglas Clowe, Frederic Courbin, Alastair Edge, Holger Israel, Mathilde Jauzac, Remy Joseph, Eric Jullo, Thomas D. Kitching, Adrienne Leonard, Julian Merten, Daisuke Nagai, James Nightingale, Andrew Robertson, Luis Javier Romualdez, Prasenjit Saha, Renske Smit, Sut Ieng Tam , et al. (1 additional authors not shown)

Abstract: We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The n… ▽ More We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications. △ Less

Submitted 14 August, 2017; originally announced August 2017.

Comments: 10 pages, 5 figures. Comments welcome