Showing 1–50 of 73 results for author: Townsend, J

Entropy Coding of Unordered Data Structures

Authors: Julius Kunze, Daniel Severo, Giulio Zani, Jan-Willem van de Meent, James Townsend

Abstract: We present shuffle coding, a general method for optimal compression of sequences of unordered objects using bits-back coding. Data structures that can be compressed using shuffle coding include multisets, graphs, hypergraphs, and others. We release an implementation that can easily be adapted to different data types and statistical models, and demonstrate that our implementation achieves state-of-… ▽ More We present shuffle coding, a general method for optimal compression of sequences of unordered objects using bits-back coding. Data structures that can be compressed using shuffle coding include multisets, graphs, hypergraphs, and others. We release an implementation that can easily be adapted to different data types and statistical models, and demonstrate that our implementation achieves state-of-the-art compression rates on a range of graph datasets including molecular data. △ Less

Submitted 16 August, 2024; originally announced August 2024.

Comments: Published at ICLR 2024

Ferroelectric phases and phase transitions in CsGeBr$_3$ induced by mechanical load

Authors: Joshua Townsend, Ravi Kashikar, S. Lisenkov, I. Ponomareva

Abstract: First-principles-based atomistic simulations are used to reveal ferroelectric phases and phase transitions induced in a semiconductor ferroelectric, CsGeBr$_3$, by external loads: hydrostatic pressure, uniaxial and biaxial stresses, and misfit strain. Hydrostatic pressure was found to suppress the Curie point at the rate -0.45$T_C(0)$ K/GPa, where $T_C(0)$ is the zero pressure Curie temperature. S… ▽ More First-principles-based atomistic simulations are used to reveal ferroelectric phases and phase transitions induced in a semiconductor ferroelectric, CsGeBr$_3$, by external loads: hydrostatic pressure, uniaxial and biaxial stresses, and misfit strain. Hydrostatic pressure was found to suppress the Curie point at the rate -0.45$T_C(0)$ K/GPa, where $T_C(0)$ is the zero pressure Curie temperature. Stresses and misfit strains were found to induce additional ferroelectric phase transitions and phases not available under normal conditions. We find that tensile load significantly enhances both the Curie temperature and spontaneous polarization, while compressive load has the opposite effect but with the difference that the Curie temperature is only slightly suppressed. The isothermal dependencies of polarization on pressure and stresses are highly nonlinear, which could result in large nonlinear piezoelectric responses. The phase diagrams reveal the diversity of the phases accessible through mechanical load, which include tetragonal, orthorhombic and monoclinic symmetries in addition to the rhombohedral and cubic ones realizable under normal conditions. We believe that this work reveals the potential of Ge-based halide perovskites for applications in energy converting devices, which is especially significant in the current pursuit of environmental friendly lead-free technologies. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Discovery of a Rare Eclipsing Be/X-ray Binary System, Swift J010902.6-723710 = SXP 182

Authors: Thomas M. Gaudin, Jamie A. Kennea, Malcolm J. Coe, Itumeleng M. Monageng, Andrzej Udalski, Lee J. Townsend, David A. H. Buckley, Phil A. Evans

Abstract: We report on the discovery of Swift J010902.6-723710, a rare eclipsing Be/X-ray Binary system by the Swift SMC Survey (S-CUBED). Swift J010902.6-723710 was discovered via weekly S-CUBED monitoring observations when it was observed to enter a state of X-ray outburst on 10 October 2023. X-ray emission was found to be modulated by a 182s period. Optical spectroscopy is used to confirm the presence of… ▽ More We report on the discovery of Swift J010902.6-723710, a rare eclipsing Be/X-ray Binary system by the Swift SMC Survey (S-CUBED). Swift J010902.6-723710 was discovered via weekly S-CUBED monitoring observations when it was observed to enter a state of X-ray outburst on 10 October 2023. X-ray emission was found to be modulated by a 182s period. Optical spectroscopy is used to confirm the presence of a highly-inclined circumstellar disk surrounding a B0-0.5Ve optical companion. Historical UV and IR photometry are then used to identify strong eclipse-like features re-occurring in both light curves with a 60.623 day period, which is adopted as the orbital period of the system. Eclipsing behavior is found to be the result of a large accretion disk surrounding the neutron star. Eclipses are produced when the disk passes in front of the OBe companion, blocking light from both the stellar surface and circumstellar disk. This is only the third Be/X-ray Binary to have confirmed eclipses. We note that this rare behavior provides an important opportunity to constrain the physical parameters of a Be/X-ray Binary with greater accuracy than is possible in non-eclipsing systems. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: 10 pages, 12 figures. To be published in the Astrophysical Journal Letters

The 2022 super-Eddington outburst of the source SMC X-2

Authors: Malcolm J. Coe, Jamie A. Kennea, Itumeleng M. Monageng, Lee J. Townsend, David A. H. Buckley, Maia Williams, Andrzej Udalski, Phil A. Evans

Abstract: SMC X-2 exhibits X-ray outburst behaviour that makes it one of the most luminous X-ray sources in the Small Magellanic Cloud. In the last decade it has undergone two such massive outbursts - in 2015 and 2022. The first outburst is well reported in the literature, but the 2022 event has yet to be fully described and discussed. That is the goal of this paper. In particular, the post-peak characteris… ▽ More SMC X-2 exhibits X-ray outburst behaviour that makes it one of the most luminous X-ray sources in the Small Magellanic Cloud. In the last decade it has undergone two such massive outbursts - in 2015 and 2022. The first outburst is well reported in the literature, but the 2022 event has yet to be fully described and discussed. That is the goal of this paper. In particular, the post-peak characteristics of the two events are compared. This reveals clear similarities in decay profiles, believed to be related to different accretion mechanisms occurring at different times as the outbursts evolve. The Hα emission line indicates that the Be disc undergoes complex structural variability, with evidence of warping as a result of its interaction with the neutron star. The detailed observations reported here will be important for modelling such interactions in this kind of binary systems. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: Accepted for publication in MNRAS

4U 1210-64: a new member of the rare intermediate-mass X-ray binary subclass

Authors: Itumeleng M. Monageng, Vanessa A. McBride, Julia Alfonso-Garzon, Lee J. Townsend, Joel B. Coley, Benjamin Montesinos, Robin H. D. Corbet, Katja Pottschmidt

Abstract: 4U 1210-64 is a peculiar X-ray binary that was first discovered in 1978 by the Uhuru satellite. The analysis of the X-ray data revealed a 6.7-day orbital period and an additional long-term modulation that is manifested as low and high flux states. Based on the previous classification of the donor star from the analysis of its optical spectra, the system has been suggested to be a high-mass X-ray b… ▽ More 4U 1210-64 is a peculiar X-ray binary that was first discovered in 1978 by the Uhuru satellite. The analysis of the X-ray data revealed a 6.7-day orbital period and an additional long-term modulation that is manifested as low and high flux states. Based on the previous classification of the donor star from the analysis of its optical spectra, the system has been suggested to be a high-mass X-ray binary. We re-visit the optical classification where, based on the spectra from the Southern African Large Telescope (SALT), we conclude that the donor star is of spectral class A8 III-IV, making it a member of the rare intermediate-mass X-ray binaries. We perform radial velocity analysis using the SALT spectra where we consider circular and eccentric orbits. From the mass function derived and the mass constraints of the donor star, we demonstrate that a neutron star is favoured as the compact object in the binary system. We show, for the first time, the folded optical lightcurves, whose shape is interpreted to be due to a combination of ellipsoidal variations, irradiation of the donor star, and mutual eclipses of the star and accretion disk. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: Accepted for publication in MNRAS

Random Edge Coding: One-Shot Bits-Back Coding of Large Labeled Graphs

Authors: Daniel Severo, James Townsend, Ashish Khisti, Alireza Makhzani

Abstract: We present a one-shot method for compressing large labeled graphs called Random Edge Coding. When paired with a parameter-free model based on Pólya's Urn, the worst-case computational and memory complexities scale quasi-linearly and linearly with the number of observed edges, making it efficient on sparse graphs, and requires only integer arithmetic. Key to our method is bits-back coding, which is… ▽ More We present a one-shot method for compressing large labeled graphs called Random Edge Coding. When paired with a parameter-free model based on Pólya's Urn, the worst-case computational and memory complexities scale quasi-linearly and linearly with the number of observed edges, making it efficient on sparse graphs, and requires only integer arithmetic. Key to our method is bits-back coding, which is used to sample edges and vertices without replacement from the edge-list in a way that preserves the structure of the graph. Optimality is proven under a class of random graph models that are invariant to permutations of the edges and of vertices within an edge. Experiments indicate Random Edge Coding can achieve competitive compression performance on real-world network datasets and scales to graphs with millions of nodes and edges. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: Published at ICML 2023

Verified Reversible Programming for Verified Lossless Compression

Authors: James Townsend, Jan-Willem van de Meent

Abstract: Lossless compression implementations typically contain two programs, an encoder and a decoder, which are required to be inverse to one another. We observe that a significant class of compression methods, based on asymmetric numeral systems (ANS), have shared structure between the encoder and decoder -- the decoder program is the 'reverse' of the encoder program -- allowing both to be simultaneousl… ▽ More Lossless compression implementations typically contain two programs, an encoder and a decoder, which are required to be inverse to one another. We observe that a significant class of compression methods, based on asymmetric numeral systems (ANS), have shared structure between the encoder and decoder -- the decoder program is the 'reverse' of the encoder program -- allowing both to be simultaneously specified by a single, reversible function. To exploit this, we have implemented a small reversible language, embedded in Agda, which we call 'Flipper' (available at https://github.com/j-towns/flipper). Agda supports formal verification of program properties, and the compiler for our reversible language (which is implemented as an Agda macro), produces not just an encoder/decoder pair of functions but also a proof that they are inverse to one another. Thus users of the language get formal verification 'for free'. We give a small example use-case of Flipper in this paper, and plan to publish a full compression implementation soon. △ Less

Submitted 23 November, 2022; v1 submitted 2 November, 2022; originally announced November 2022.

Gamma-ray Eclipses and Orbital Modulation Transitions in the Candidate Redback 4FGL J1702.7-5655

Authors: R. H. D. Corbet, L. Chomiuk, J. B. Coley, G. Dubus, P. G. Edwards, N. Islam, V. A. McBride, J. Stevens, J. Strader, S. J. Swihart, L. J. Townsend

Abstract: Observations with the Fermi Large Area Telescope (LAT) of the gamma-ray source 4FGL J1702.7-5655, previously classified as a candidate millisecond pulsar, show highly-significant modulation at a period of 0.2438033 days (~ 5.85 hours). Further examination of the folded light curve indicates the presence of narrow eclipses, suggesting this is a redback binary system. An examination of the long-term… ▽ More Observations with the Fermi Large Area Telescope (LAT) of the gamma-ray source 4FGL J1702.7-5655, previously classified as a candidate millisecond pulsar, show highly-significant modulation at a period of 0.2438033 days (~ 5.85 hours). Further examination of the folded light curve indicates the presence of narrow eclipses, suggesting this is a redback binary system. An examination of the long-term properties of the modulation over 13 years of LAT observations indicates that the orbital modulation of the gamma-rays changed from a simple eclipse before early 2013, to a broader, more easily detected, quasi-sinusoidal modulation. In addition, the time of the eclipse shifts to ~0.05 later in phase. This change in the orbital modulation properties is, however, not accompanied by a significant overall change in gamma-ray flux or spectrum. The quasi-sinusoidal component peaks ~0.5 out of phase with the eclipse, which would indicate inferior conjunction of the compact object in the system. Swift X-ray Telescope observations reveal a possible X-ray counterpart within the LAT error ellipse. However, radio observations obtained with the Australia Telescope Compact Array do not detect a source in the region. 4FGL J1702.7-5655 appears to have changed its state in 2013, perhaps related to changes in the intrabinary shock in the system. We discuss how the properties of 4FGL J1702.7-5655 compare to other binary millisecond pulsars that have exhibited orbital modulation in gamma rays. △ Less

Submitted 23 May, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Comments: Accepted for publication in the Astrophysical Journal

Disentangling the neighbouring pulsars SXP 15.3 and SXP 305

Authors: Itumeleng M. Monageng, Malcolm J. Coe, Lee J. Townsend, Silas G. T. Laycock, Jamie A. Kennea, Ankur Roy, Andrzej Udalski, Sayantan Bhattacharya, Dimitris M. Christodoulou, David A. H. Buckley, Phil A. Evans

Abstract: SXP 15.3 and SXP 305 are two Be X-ray binaries in the Small Magellanic Cloud that are spatially separated by ~7 arcsec. The small separation between these sources has, in the past, resulted in confusion about the origin of the emission from the combined region. We present long-term optical and X-ray monitoring results of both sources, where we study the historic and recent behaviour. In particular… ▽ More SXP 15.3 and SXP 305 are two Be X-ray binaries in the Small Magellanic Cloud that are spatially separated by ~7 arcsec. The small separation between these sources has, in the past, resulted in confusion about the origin of the emission from the combined region. We present long-term optical and X-ray monitoring results of both sources, where we study the historic and recent behaviour. In particular, from data collected as part of the S-CUBED project we see repeating X-ray outbursts from the combined region of the two sources in the recent lightcurve from the Neil Gehrels Swift Observatory, and we investigate the origin of this emission. Using the H-alpha emission line from the Southern African Large Telescope (SALT) and photometric flux from the Optical Gravitational Lensing Experiment (OGLE) to study the changes in the size and structure of the Be disc, we demonstrate that the X-ray emission likely originates from SXP 15.3. Timing analysis reveals unusual behaviour, where the optical outburst profile shows modulation at twice the frequency of the X-ray outbursts. We consider either of these periodicities being the true orbital period in SXP 15.3 and propose models based on the geometric orientations of the Be disc and neutron star to explain the physical origin of the outbursts. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: Accepted for publication in MNRAS

Parallel Neural Local Lossless Compression

Authors: Mingtian Zhang, James Townsend, Ning Kang, David Barber

Abstract: The recently proposed Neural Local Lossless Compression (NeLLoC), which is based on a local autoregressive model, has achieved state-of-the-art (SOTA) out-of-distribution (OOD) generalization performance in the image compression task. In addition to the encouragement of OOD generalization, the local model also allows parallel inference in the decoding stage. In this paper, we propose two paralleli… ▽ More The recently proposed Neural Local Lossless Compression (NeLLoC), which is based on a local autoregressive model, has achieved state-of-the-art (SOTA) out-of-distribution (OOD) generalization performance in the image compression task. In addition to the encouragement of OOD generalization, the local model also allows parallel inference in the decoding stage. In this paper, we propose two parallelization schemes for local autoregressive models. We discuss the practicalities of implementing the schemes and provide experimental evidence of significant gains in compression runtime compared to the previous, non-parallel implementation. △ Less

Submitted 26 June, 2022; v1 submitted 13 January, 2022; originally announced January 2022.

Adaptive Optimization with Examplewise Gradients

Authors: Julius Kunze, James Townsend, David Barber

Abstract: We propose a new, more general approach to the design of stochastic gradient-based optimization methods for machine learning. In this new framework, optimizers assume access to a batch of gradient estimates per iteration, rather than a single estimate. This better reflects the information that is actually available in typical machine learning setups. To demonstrate the usefulness of this generaliz… ▽ More We propose a new, more general approach to the design of stochastic gradient-based optimization methods for machine learning. In this new framework, optimizers assume access to a batch of gradient estimates per iteration, rather than a single estimate. This better reflects the information that is actually available in typical machine learning setups. To demonstrate the usefulness of this generalized approach, we develop Eve, an adaptation of the Adam optimizer which uses examplewise gradients to obtain more accurate second-moment estimates. We provide preliminary experiments, without hyperparameter tuning, which show that the new optimizer slightly outperforms Adam on a small scale benchmark and performs the same or worse on larger scale benchmarks. Further work is needed to refine the algorithm and tune hyperparameters. △ Less

Submitted 30 November, 2021; originally announced December 2021.

Comments: 9 pages, 1 figure, 3 tables

Identification of an X-ray Pulsar in the BeXRB system IGR J18219$-$1347

Authors: B. O'Connor, E. Gogus, D. Huppenkothen, C. Kouveliotou, N. Gorgone, L. J. Townsend, A. Calamida, A. Fruchter, D. A. H. Buckley, M. G. Baring, J. A. Kennea, G. Younes, Z. Arzoumanian, E. Bellm, S. B. Cenko, K. Gendreau, J. Granot, C. Hailey, F. Harrison, D. Hartmann, L. Kaper, A. Kutyrev, P. O. Slane, D. Stern, E. Troja , et al. (3 additional authors not shown)

Abstract: We report on observations of the candidate Be/X-ray binary IGR J18219$-$1347 with \textit{Swift}/XRT, \textit{NuSTAR}, and \textit{NICER} during Type-I outbursts in March and June 2020. Our timing analysis revealed the spin period of a neutron star with $P_\textrm{spin}=52.46$ s. This periodicity, combined with the known orbital period of $72.4$ d, indicates that the system is a BeXRB. Furthermore… ▽ More We report on observations of the candidate Be/X-ray binary IGR J18219$-$1347 with \textit{Swift}/XRT, \textit{NuSTAR}, and \textit{NICER} during Type-I outbursts in March and June 2020. Our timing analysis revealed the spin period of a neutron star with $P_\textrm{spin}=52.46$ s. This periodicity, combined with the known orbital period of $72.4$ d, indicates that the system is a BeXRB. Furthermore, by comparing the infrared counterpart's spectral energy distribution to known BeXRBs, we confirm this classification and set a distance of approximately $10-15$ kpc for the source. The source's broadband X-ray spectrum ($1.5-50$ keV) is described by an absorbed power-law with photon index $Γ$\,$\sim$\,$0.5$ and cutoff energy at $\sim$\,$13$ keV. △ Less

Submitted 1 February, 2022; v1 submitted 24 November, 2021; originally announced November 2021.

Comments: Revised manuscript. Accepted to ApJ

The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview

Authors: Devin Crichton, Moumita Aich, Adam Amara, Kevin Bandura, Bruce A. Bassett, Carlos Bengaly, Pascale Berner, Shruti Bhatporia, Martin Bucher, Tzu-Ching Chang, H. Cynthia Chiang, Jean-Francois Cliche, Carolyn Crichton, Romeel Dave, Dirk I. L. de Villiers, Matt A. Dobbs, Aaron M. Ewall-Wice, Scott Eyono, Christopher Finlay, Sindhu Gaddam, Ken Ganga, Kevin G. Gayley, Kit Gerodias, Tim Gibbon, Austin Gumba , et al. (75 additional authors not shown)

Abstract: The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 40… ▽ More The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array. △ Less

Submitted 17 January, 2022; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: 24 pages, 6 figures. Updated to match published version

Journal ref: J. of Astronomical Telescopes, Instruments, and Systems, 8(1), 011019 (2022)

Swift J011511.0-725611: Discovery of a rare Be Star / White Dwarf binary system in the SMC

Authors: J. A. Kennea, M. J. Coe, P. A. Evans, L. J. Townsend, Z. A. Campbell, A. Udalski

Abstract: We report on the discovery of Swift J011511.0-725611, a rare Be X-ray binary system (BeXRB) with a White Dwarf (WD) compact object, in the Small Magellanic Cloud (SMC) by S-CUBED, a weekly X-ray/UV survey of the SMC by the Neil Gehrels Swift Observatory. Observations show an approximately 3 month outburst from Swift J011511.0-725611, the first detected by S-CUBED since it began in 2016 June. Swift… ▽ More We report on the discovery of Swift J011511.0-725611, a rare Be X-ray binary system (BeXRB) with a White Dwarf (WD) compact object, in the Small Magellanic Cloud (SMC) by S-CUBED, a weekly X-ray/UV survey of the SMC by the Neil Gehrels Swift Observatory. Observations show an approximately 3 month outburst from Swift J011511.0-725611, the first detected by S-CUBED since it began in 2016 June. Swift J011511.0-725611 shows super-soft X-ray emission, indicative of a White Dwarf compact object, which is further strengthened by the presence of an 0.871 keV edge, commonly attributed to O viii K-edge in the WD atmosphere. Spectroscopy by SALT confirms the Be nature of the companion star, and long term light-curve by OGLE finds both the signature of a circumstellar disk in the system at outburst time, and the presence of a 17.4 day periodicity, likely the orbital period of the system. Swift J011511.0-725611 is suggested to be undergoing a Type-II outburst, similar to the previously reported SMC Be White Dwarf binary (BeWD), Swift J004427.3-734801. It is likely that the rarity of known BeWD is in part due to the difficulty in detecting such outbursts due to both their rarity, and their relative faintness compared to outbursts in Neutron Star BeXRBs. △ Less

Submitted 11 September, 2021; originally announced September 2021.

Comments: 8 pages, 10 figures, accepted for publication in MNRAS

Compressing Multisets with Large Alphabets using Bits-Back Coding

Authors: Daniel Severo, James Townsend, Ashish Khisti, Alireza Makhzani, Karen Ullrich

Abstract: Current methods which compress multisets at an optimal rate have computational complexity that scales linearly with alphabet size, making them too slow to be practical in many real-world settings. We show how to convert a compression algorithm for sequences into one for multisets, in exchange for an additional complexity term that is quasi-linear in sequence length. This allows us to compress mult… ▽ More Current methods which compress multisets at an optimal rate have computational complexity that scales linearly with alphabet size, making them too slow to be practical in many real-world settings. We show how to convert a compression algorithm for sequences into one for multisets, in exchange for an additional complexity term that is quasi-linear in sequence length. This allows us to compress multisets of exchangeable symbols at an optimal rate, with computational complexity decoupled from the alphabet size. The key insight is to avoid encoding the multiset directly, and instead compress a proxy sequence, using a technique called `bits-back coding'. We demonstrate the method experimentally on tasks which are intractable with previous optimal-rate methods: compression of multisets of images and JavaScript Object Notation (JSON) files. Code for our experiments is available at https://github.com/facebookresearch/multiset-compression. △ Less

Submitted 27 February, 2023; v1 submitted 15 July, 2021; originally announced July 2021.

Journal ref: IEEE Journal on Selected Areas in Information Theory, 2023

SRGA J124404.1-632232/SRGU J124403.8-632231: a new X-ray pulsar discovered in the all-sky survey by SRG

Authors: V. Doroshenko, R. Staubert, C. Maitra, A. Rau, F. Haberl, A. Santangelo, A. Schwope, J. Wilms, D. A. H. Buckley, A. Semena, I. Mereminskiy, A. Lutovinov, M. Gromadzki, L. J. Townsend, I. M. Monageng

Abstract: Ongoing all-sky surveys by the the eROSITA and the Mikhail Pavlinsky ART-XC telescopes on-board the Spectrum Roentgen Gamma (SRG) mission have already revealed over a million of X-ray sources. One of them, SRGA J124404.1-632232/SRGU J124403.8-632231, was detected as a new source in the third (of the planned eight) consecutive X-ray surveys by ART-XC. Based on the properties of the identified optic… ▽ More Ongoing all-sky surveys by the the eROSITA and the Mikhail Pavlinsky ART-XC telescopes on-board the Spectrum Roentgen Gamma (SRG) mission have already revealed over a million of X-ray sources. One of them, SRGA J124404.1-632232/SRGU J124403.8-632231, was detected as a new source in the third (of the planned eight) consecutive X-ray surveys by ART-XC. Based on the properties of the identified optical counterpart it was classified as a candidate X-ray binary (XRB). We report on the follow-up observations of this source with Nuclear Spectroscopic Telescope Array (NuSTAR), Neil Gehrels Swift Observatory (Swift), and the Southern African Large Telescope (SALT), which allowed us to unambiguously confirm the initial identification and establish SRGU J124403.8-632231 as a new X-ray pulsar with a spin period of ~538 s and a Be-star companion, making it one of the first Galactic X-ray pulsars discovered by SRG. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

Journal ref: A&A 661, A21 (2022)

Identification of SRGt 062340.2-265715 as a bright, strongly variable, novalike cataclysmic variable

Authors: Axel Schwope, David A. H. Buckley, Adela Kawka, Ole König, Alexander Lutovinov, Chandreyee Maitra, Ilya Mereminskiy, James Miller-Jones, Manuel Pichardo Marcano, Arne Rau, Andrei Semena, Lee J. Townsend, Jörn Wilms

Abstract: We report the identification and follow-up of the transient SRG 062340.2-265715 detected with both instruments on board the Spektrum-Roentgen-Gamma mission. Optical spectroscopy of the G=12.5 counterpart firmly classifies the object as a novalike cataclysmic variable (CV) at a distance of 495 pc. A highly significant TESS period of 3.941 hours, tentatively identified with the orbital period of the… ▽ More We report the identification and follow-up of the transient SRG 062340.2-265715 detected with both instruments on board the Spektrum-Roentgen-Gamma mission. Optical spectroscopy of the G=12.5 counterpart firmly classifies the object as a novalike cataclysmic variable (CV) at a distance of 495 pc. A highly significant TESS period of 3.941 hours, tentatively identified with the orbital period of the binary, could not be found when the object was reobserved with TESS two years later. The newer high-cadence TESS data revealed quasi-periodic oscillations around 25 min, while ground-based photometry indicated periodic variability at 32 min. Located in very sparsely populated regions of color-magnitude diagrams involving X-ray and optical magnitudes and colors, the new object could be an X-ray underluminous magnetic CV, an intermediate polar, or an overluminous nonmagnetic CV. The lack of uniquely identified spin and orbital periods prevents a final classification. The site of X-ray production in the system, L(X, bol) = 4.8 x 10^{32} erg/s, remains to be understood given its high variability on long and short timescales. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: 9 pages, 16 figures, to appear on A&A, Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

Lossless Compression with Latent Variable Models

Authors: James Townsend

Abstract: We develop a simple and elegant method for lossless compression using latent variable models, which we call 'bits back with asymmetric numeral systems' (BB-ANS). The method involves interleaving encode and decode steps, and achieves an optimal rate when compressing batches of data. We demonstrate it firstly on the MNIST test set, showing that state-of-the-art lossless compression is possible using… ▽ More We develop a simple and elegant method for lossless compression using latent variable models, which we call 'bits back with asymmetric numeral systems' (BB-ANS). The method involves interleaving encode and decode steps, and achieves an optimal rate when compressing batches of data. We demonstrate it firstly on the MNIST test set, showing that state-of-the-art lossless compression is possible using a small variational autoencoder (VAE) model. We then make use of a novel empirical insight, that fully convolutional generative models, trained on small images, are able to generalize to images of arbitrary size, and extend BB-ANS to hierarchical latent variable models, enabling state-of-the-art lossless compression of full-size colour images from the ImageNet dataset. We describe 'Craystack', a modular software framework which we have developed for rapid prototyping of compression using deep generative models. △ Less

Submitted 22 April, 2021; v1 submitted 21 April, 2021; originally announced April 2021.

The Be/neutron star system Swift J004929.5-733107 in the Small Magellanic Cloud -- X-ray characteristics and optical counterpart candidates

Authors: M. J. Coe, J. A. Kennea, P. A. Evans, L. J. Townsend, A. Udalski, I. M. Monageng, D. A. H. Buckley

Abstract: Swift J004929.5-733107 is an X-ray source in the Small Magellanic Cloud (SMC) that has been reported several times, but the optical counterpart has been unclear due to source confusion in a crowded region of the SMC. Previous works proposed [MA93] 302 as the counterpart, however we show here, using data obtained from the S-CUBED project, that the X-ray positio is inconsistent with that object. Ins… ▽ More Swift J004929.5-733107 is an X-ray source in the Small Magellanic Cloud (SMC) that has been reported several times, but the optical counterpart has been unclear due to source confusion in a crowded region of the SMC. Previous works proposed [MA93] 302 as the counterpart, however we show here, using data obtained from the S-CUBED project, that the X-ray positio is inconsistent with that object. Instead we propose a previously unclassified object which has all the indications of being a newly identified Be star exhibiting strong HU emission. Evidence for the presence of significant I-band variability strongly suggest that this is, in fact, a Be type star with a large circumstellar disk. Over 18 years worth of optical monitoring by the OGLE project reveal a periodic modulation at a period of 413d, probably the binary period of the system. A SALT optical spectrum shows strong Balmer emission and supports a proposed spectral classification of B1-3 III-IVe. The X-ray data obtained from the S-CUBED project reveal a time-averaged spectrum well fitted by a photon index = 0.93 pm 0.16. Assuming the known distance to the SMC the flux corresponds to a luminosity 10E35 erg/s. All of these observational facts suggest that this is confirmed as a Be star-neutron star X-ray binary (BeXRB) in the SMC, albeit one with an unusually long binary period at the limits of the Corbet Diagram. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 9 pages 12 figures

Lossless compression with state space models using bits back coding

Authors: James Townsend, Iain Murray

Abstract: We generalize the 'bits back with ANS' method to time-series models with a latent Markov structure. This family of models includes hidden Markov models (HMMs), linear Gaussian state space models (LGSSMs) and many more. We provide experimental evidence that our method is effective for small scale models, and discuss its applicability to larger scale settings such as video compression. We generalize the 'bits back with ANS' method to time-series models with a latent Markov structure. This family of models includes hidden Markov models (HMMs), linear Gaussian state space models (LGSSMs) and many more. We provide experimental evidence that our method is effective for small scale models, and discuss its applicability to larger scale settings such as video compression. △ Less

Submitted 4 May, 2021; v1 submitted 18 March, 2021; originally announced March 2021.

Improving Lossless Compression Rates via Monte Carlo Bits-Back Coding

Authors: Yangjun Ruan, Karen Ullrich, Daniel Severo, James Townsend, Ashish Khisti, Arnaud Doucet, Alireza Makhzani, Chris J. Maddison

Abstract: Latent variable models have been successfully applied in lossless compression with the bits-back coding algorithm. However, bits-back suffers from an increase in the bitrate equal to the KL divergence between the approximate posterior and the true posterior. In this paper, we show how to remove this gap asymptotically by deriving bits-back coding algorithms from tighter variational bounds. The key… ▽ More Latent variable models have been successfully applied in lossless compression with the bits-back coding algorithm. However, bits-back suffers from an increase in the bitrate equal to the KL divergence between the approximate posterior and the true posterior. In this paper, we show how to remove this gap asymptotically by deriving bits-back coding algorithms from tighter variational bounds. The key idea is to exploit extended space representations of Monte Carlo estimators of the marginal likelihood. Naively applied, our schemes would require more initial bits than the standard bits-back coder, but we show how to drastically reduce this additional cost with couplings in the latent space. When parallel architectures can be exploited, our coders can achieve better rates than bits-back with little additional cost. We demonstrate improved lossless compression rates in a variety of settings, especially in out-of-distribution or sequential data compression. △ Less

Submitted 14 June, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

ERIC: Extracting Relations Inferred from Convolutions

Authors: Joe Townsend, Theodoros Kasioumis, Hiroya Inakoshi

Abstract: Our main contribution is to show that the behaviour of kernels across multiple layers of a convolutional neural network can be approximated using a logic program. The extracted logic programs yield accuracies that correlate with those of the original model, though with some information loss in particular as approximations of multiple layers are chained together or as lower layers are quantised. We… ▽ More Our main contribution is to show that the behaviour of kernels across multiple layers of a convolutional neural network can be approximated using a logic program. The extracted logic programs yield accuracies that correlate with those of the original model, though with some information loss in particular as approximations of multiple layers are chained together or as lower layers are quantised. We also show that an extracted program can be used as a framework for further understanding the behaviour of CNNs. Specifically, it can be used to identify key kernels worthy of deeper inspection and also identify relationships with other kernels in the form of the logical rules. Finally, we make a preliminary, qualitative assessment of rules we extract from the last convolutional layer and show that kernels identified are symbolic in that they react strongly to sets of similar images that effectively divide output classes into sub-classes with distinct characteristics. △ Less

Submitted 19 October, 2020; originally announced October 2020.

Comments: Accepted for poster presentation at ACCV (Asian Conference on Computer Vision) 2020

HI intensity mapping with the MIGHTEE survey: power spectrum estimates

Authors: Sourabh Paul, Mario G. Santos, Junaid Townsend, Matt J. Jarvis, Natasha Maddox, Jordan D. Collier, Bradley S. Frank, Russ Taylor

Abstract: Intensity mapping (IM) with neutral hydrogen is a promising avenue to probe the large scale structure of the Universe. In this paper, we demonstrate that using the 64-dish MeerKAT radio telescope as a connected interferometer, it is possible to make a statistical detection of HI in the post-reionization Universe. With the MIGHTEE (MeerKAT International GHz Tiered Extragalactic Exploration) survey… ▽ More Intensity mapping (IM) with neutral hydrogen is a promising avenue to probe the large scale structure of the Universe. In this paper, we demonstrate that using the 64-dish MeerKAT radio telescope as a connected interferometer, it is possible to make a statistical detection of HI in the post-reionization Universe. With the MIGHTEE (MeerKAT International GHz Tiered Extragalactic Exploration) survey project observing in the L-band ($856 < ν< 1712$ MHz, $z < 0.66$), we can achieve the required sensitivity to measure the HI IM power spectrum on quasi-linear scales, which will provide an important complementarity to the single-dish IM MeerKAT observations. We present a purpose-built simulation pipeline that emulates the MIGHTEE observations and forecast the constraints that can be achieved on the HI power spectrum at $z = 0.27$ for $k > 0.3$ $\rm{Mpc}^{-1}$ using the foreground avoidance method. We present the power spectrum estimates with the current simulation on the COSMOS field that includes contributions from HI, noise and point source models constructed from the observed MIGHTEE data. The results from our \textit{visibility} based pipeline are in qualitative agreement to the already available MIGHTEE data. This paper demonstrates that MeerKAT can achieve very high sensitivity to detect HI with the full MIGHTEE survey on quasi-linear scales (signal-to-noise ratio $> 7$ at $k=0.49$ $\rm{Mpc}^{-1}$) which are instrumental in probing cosmological quantities such as the spectral index of fluctuation, constraints on warm dark matter, the quasi-linear redshift space distortions and the measurement of the HI content of the Universe up to $z\sim 0.5$. △ Less

Submitted 15 April, 2021; v1 submitted 28 September, 2020; originally announced September 2020.

Comments: 13 pages, 8 figures, accepted for publication in MNRAS

The newly discovered Be/X-ray binary Swift J004516.6-734703 in the SMC: witnessing the emergence of a circumstellar disc

Authors: J. A. Kennea, M. J. Coe, P. A. Evans, I. M. Monageng, L. J. Townsend, M. H. Siegel, A. Udalski, D. A. H. Buckley

Abstract: We report on the discovery of Swift J004516.6-734703, a Be/X-ray binary system by the Swift SMC Survey, S-CUBED. Swift J004516.6-734703, or SXP 146.6, was found to be exhibiting a bright (~10^37 erg/s) X-ray outburst in 2020 June 18. The historical UV and IR light-curves from OGLE and Swift/UVOT showed that after a long period of steady brightness, it experienced a significant brightening beginnin… ▽ More We report on the discovery of Swift J004516.6-734703, a Be/X-ray binary system by the Swift SMC Survey, S-CUBED. Swift J004516.6-734703, or SXP 146.6, was found to be exhibiting a bright (~10^37 erg/s) X-ray outburst in 2020 June 18. The historical UV and IR light-curves from OGLE and Swift/UVOT showed that after a long period of steady brightness, it experienced a significant brightening beginning around 2019 March. This IR/UV rise is likely the signature of the formation of a circumstellar disc, confirmed by the presence of strong a Hα line in SALT spectroscopy, that was not previously present. Periodicity analysis of the OGLE data reveals a plausible 426 day binary period, and in the X-ray a pulsation period of 146.6s period is detected. The onset of X-ray emission from Swift J004516.6-734703 is likely the signature of a Type-I outburst from the first periastron passage of the neutron star companion through the newly formed circumstellar disc. We note that the formation of the circumstellar disc began at the predicted time of the previous periastron passage, suggesting its formation was spurred by tidal interaction with the neutron star. △ Less

Submitted 25 August, 2020; v1 submitted 22 June, 2020; originally announced June 2020.

Comments: 6 Pages, 7 Figures, 1 Table. Accepted by MNRAS

The effect of extended closure of red-light areas on COVID-19 transmission in India

Authors: Abhishek Pandey, Sudhakar V. Nuti, Pratha Sah, Chad R. Wells, Alison P. Galvani, Jeffrey P. Townsend

Abstract: The novel coronavirus disease (COVID-19) pandemic has resulted in over 200,000 cases in India. Thus far, India has implemented lockdown measures to curb disease transmission. However, commercial sex work in red-light areas (RLAs) has potential to lead to COVID-19 resurgence after lockdown. We developed a model of COVID-19 transmission in RLAs, evaluating the impact of extended RLA closure compared… ▽ More The novel coronavirus disease (COVID-19) pandemic has resulted in over 200,000 cases in India. Thus far, India has implemented lockdown measures to curb disease transmission. However, commercial sex work in red-light areas (RLAs) has potential to lead to COVID-19 resurgence after lockdown. We developed a model of COVID-19 transmission in RLAs, evaluating the impact of extended RLA closure compared with RLA reopening on cases, hospitalizations, and mortality rates within the RLAs of five major Indian cities, within the cities, and across India. Closure lowered transmission at all scales. More than 90% of cumulative cases and deaths among RLA residents of Kolkata, Pune, and Nagpur could be averted by the time the epidemic would peak under a re-opening scenario. Across India, extended closure of RLAs would benefit the population at large, delaying the peak of COVID-19 cases by 8 to 23 days, and avert 32% to 60.2% of cumulative cases and 43% to 67.6% of cumulative deaths at the peak of the epidemic. Extended closure of RLAs until better prevention and treatment strategies are developed would benefit public health in India. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Simultaneous multi-telescope observations of FRB 121102

Authors: M. Caleb, B. W. Stappers, T. D. Abbott, E. D. Barr, M. C. Bezuidenhout, S. J. Buchner, M. Burgay, W. Chen, I. Cognard, L. N. Driessen, R. Fender, G. H. Hilmarsson, J. Hoang, D. M. Horn, F. Jankowski, M. Kramer, D. R. Lorimer, M. Malenta, V. Morello, M. Pilia, E. Platts, A. Possenti, K. M. Rajwade, A. Ridolfi, L. Rhodes , et al. (7 additional authors not shown)

Abstract: We present 11 detections of FRB 121102 in ~3 hours of observations during its 'active' period on the 10th of September 2019. The detections were made using the newly deployed MeerTRAP system and single pulse detection pipeline at the MeerKAT radio telescope in South Africa. Fortuitously, the Nancay radio telescope observations on this day overlapped with the last hour of MeerKAT observations and r… ▽ More We present 11 detections of FRB 121102 in ~3 hours of observations during its 'active' period on the 10th of September 2019. The detections were made using the newly deployed MeerTRAP system and single pulse detection pipeline at the MeerKAT radio telescope in South Africa. Fortuitously, the Nancay radio telescope observations on this day overlapped with the last hour of MeerKAT observations and resulted in 4 simultaneous detections. The observations with MeerKAT's wide band receiver, which extends down to relatively low frequencies (900-1670 MHz usable L-band range), have allowed us to get a detailed look at the complex frequency structure, intensity variations and frequency-dependent sub-pulse drifting. The drift rates we measure for the full-band and sub-banded data are consistent with those published between 600-6500 MHz with a slope of -0.147 +/- 0.014 ms^-1. Two of the detected bursts exhibit fainter 'precursors' separated from the brighter main pulse by ~28 ms and ~34 ms. A follow-up multi-telescope campaign on the 6th and 8th October 2019 to better understand these frequency drifts and structures over a wide and continuous band was undertaken. No detections resulted, indicating that the source was 'inactive' over a broad frequency range during this time. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS

Orbital and super-orbital periods in ULX pulsars, disc-fed HMXBs, Be/X-ray binaries and double-periodic variables

Authors: L. J. Townsend, P. A. Charles

Abstract: We present evidence for a simple linear relationship between the orbital period and super-orbital period in ultra-luminous X-ray (ULX) pulsars, akin to what is seen in the population of disc-fed neutron star super-giant X-ray binary and Be/X-ray binary systems. We argue that the most likely cause of this relationship is the modulation of precessing hot spots or density waves in an accretion or cir… ▽ More We present evidence for a simple linear relationship between the orbital period and super-orbital period in ultra-luminous X-ray (ULX) pulsars, akin to what is seen in the population of disc-fed neutron star super-giant X-ray binary and Be/X-ray binary systems. We argue that the most likely cause of this relationship is the modulation of precessing hot spots or density waves in an accretion or circumstellar disc by the binary motion of the system, implying a physical link between ULX pulsars and high-mass X-ray binary (HMXB) pulsars. This hypothesis is supported by recent studies of Galactic and Magellanic Cloud HMXBs accreting at super-Eddington rates, and the position of ULX pulsars on the spin period--orbital period diagram of HMXBs. An interesting secondary relationship discovered in this work is the apparent connection between disc-fed HMXBs, ULXs and a seemingly unrelated group of early-type binaries showing so-called "double-periodic" variability. We suggest that these systems are good candidates to be the direct progenitors of Be/X-ray binaries. △ Less

Submitted 28 April, 2020; originally announced April 2020.

Comments: 5 pages, accepted for publication in MNRAS Letters

SOFIA-EXES Observations of Betelgeuse during the Great Dimming of 2019/2020

Authors: Graham M. Harper, Curtis N. DeWitt, Matthew J. Richter, Edward F. Guinan, Richard Wasatonic, Nils Ryde, Edward, J. Montiel, Amanda J. Townsend

Abstract: In 2019 October Betelgeuse began a decline in V-band brightness that went beyond the minimum expected from its quasi-periodic ~420 day cycle, becoming the faintest in recorded photometric history. Observations obtained in 2019 December with VLT/SPHERE (Montarges 2020) have shown that the southern half of the star has become markedly fainter than in 2019 January indicating that a major change has o… ▽ More In 2019 October Betelgeuse began a decline in V-band brightness that went beyond the minimum expected from its quasi-periodic ~420 day cycle, becoming the faintest in recorded photometric history. Observations obtained in 2019 December with VLT/SPHERE (Montarges 2020) have shown that the southern half of the star has become markedly fainter than in 2019 January indicating that a major change has occurred in, or near, the photosphere. We present SOFIA-EXES high spectral-resolution observations of [Fe II] 25.99 mic and [S I] 25.25 mic emission lines from Betelgeuse obtained during the unprecedented 2020 February V-band brightness minimum to investigate potential changes in the circumstellar flow. These spectra are compared to observations obtained in 2015 and 2017 when the V magnitude was typical of brighter phases. We find only very small changes in the gas velocities reflected by either of the line profiles, no significant changes in the flux to continuum ratios, and hence no significant changes in the [Fe II]/[S I] flux ratios. There is evidence that absorption features have appeared in the 2020 continuum. The Alfven wave-crossing time from the upper-photosphere is sufficiently long that one would not expect a change in the large scale magnetic field to reach the circumstellar [Fe II] and [S I] line forming regions, 3 < R(R*) < 20. However, the light-crossing time is of order a few hours and a reduction in luminosity may reduce the dust-gas heating rate and [O I] 63.19 mic emission which has contributions from R > 20R*, where significant circumstellar oxygen-rich dust is observed. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 20 pages, 5 figures, ApJL Accepted

A tutorial on the range variant of asymmetric numeral systems

Authors: James Townsend

Abstract: This paper is intended to be a brief and accessible introduction to the range variant of asymmetric numeral systems (ANS), a system for lossless compression of sequences which can be used as a drop in replacement for arithmetic coding (AC). Because of the relative simplicity of ANS, we are able to provide enough mathematical detail to rigorously prove that ANS attains a compression rate close to t… ▽ More This paper is intended to be a brief and accessible introduction to the range variant of asymmetric numeral systems (ANS), a system for lossless compression of sequences which can be used as a drop in replacement for arithmetic coding (AC). Because of the relative simplicity of ANS, we are able to provide enough mathematical detail to rigorously prove that ANS attains a compression rate close to the Shannon limit. Pseudo-code, intuitive interpretation and diagrams are given alongside the mathematical derivations. A working Python demo which accompanies this tutorial is available at https://raw.githubusercontent.com/j-towns/ans-notes/master/rans.py. △ Less

Submitted 7 October, 2020; v1 submitted 24 January, 2020; originally announced January 2020.

HiLLoC: Lossless Image Compression with Hierarchical Latent Variable Models

Authors: James Townsend, Thomas Bird, Julius Kunze, David Barber

Abstract: We make the following striking observation: fully convolutional VAE models trained on 32x32 ImageNet can generalize well, not just to 64x64 but also to far larger photographs, with no changes to the model. We use this property, applying fully convolutional models to lossless compression, demonstrating a method to scale the VAE-based 'Bits-Back with ANS' algorithm for lossless compression to large… ▽ More We make the following striking observation: fully convolutional VAE models trained on 32x32 ImageNet can generalize well, not just to 64x64 but also to far larger photographs, with no changes to the model. We use this property, applying fully convolutional models to lossless compression, demonstrating a method to scale the VAE-based 'Bits-Back with ANS' algorithm for lossless compression to large color photographs, and achieving state of the art for compression of full size ImageNet images. We release Craystack, an open source library for convenient prototyping of lossless compression using probabilistic models, along with full implementations of all of our compression results. △ Less

Submitted 20 December, 2019; originally announced December 2019.

The Shock Physics of Giant Impacts: Key Requirements for the Equations of State

Authors: Sarah T. Stewart, Erik J. Davies, Megan S. Duncan, Simon J. Lock, Seth Root, Joshua P. Townsend, Richard G. Kraus, Razvan Caracas, Stein B. Jacobsen

Abstract: We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an upda… ▽ More We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg$_2$SiO$_4$), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model. △ Less

Submitted 16 January, 2020; v1 submitted 10 October, 2019; originally announced October 2019.

Comments: M-ANEOS source code public release (https://github.com/isale-code/M-ANEOS); forsterite EOS public release (https://github.com/ststewart/aneos-forsterite-2019); manuscript accepted (no changes): Stewart, S., et al. (accepted). In J. Lane, T. Germann, and M. Armstrong (Eds.), 21st Biennial APS Conference on Shock Compression of Condensed Matter (SCCM19). AIP Publishing

Discovery of the Galactic High-Mass Gamma-ray Binary 4FGL J1405.1-6119

Authors: R. H. D. Corbet, L. Chomiuk, M. J. Coe, J. B. Coley, G. Dubus, P. G. Edwards, P. Martin, V. A. McBride, J. Stevens, J. Strader, L. J. Townsend

Abstract: We report the identification from multi-wavelength observations of the Fermi Large Area Telescope (LAT) source 4FGL J1405.1-6119 (= 3FGL J1405.4-6119) as a high-mass gamma-ray binary. Observations with the LAT show that gamma-ray emission from the system is modulated at a period of 13.7135 +/- 0.0019 days, with the presence of two maxima per orbit with different spectral properties. X-ray observat… ▽ More We report the identification from multi-wavelength observations of the Fermi Large Area Telescope (LAT) source 4FGL J1405.1-6119 (= 3FGL J1405.4-6119) as a high-mass gamma-ray binary. Observations with the LAT show that gamma-ray emission from the system is modulated at a period of 13.7135 +/- 0.0019 days, with the presence of two maxima per orbit with different spectral properties. X-ray observations using the Neil Gehrels Swift Observatory X-ray Telescope (XRT) show that X-ray emission is also modulated at this period, but with a single maximum that is closer to the secondary lower-energy gamma-ray maximum. A radio source, coincident with the X-ray source, is also found from Australia Telescope Compact Array (ATCA) observations, and the radio emission is modulated on the gamma-ray period with similar phasing to the X-ray emission. A large degree of interstellar obscuration severely hampers optical observations, but a near-infrared counterpart is found. Near-infrared spectroscopy indicates an O6 III spectral classification. This is the third gamma-ray binary to be discovered with the Fermi LAT from periodic modulation of the gamma-ray emission, the other two sources also have early O star, rather than Be star, counterparts. We consider at what distances we can detect such modulated gamma-ray emission with the LAT, and examine constraints on the gamma-ray binary population of the Milky Way. △ Less

Submitted 28 August, 2019; originally announced August 2019.

Comments: Accepted for publication in the Astrophysical Journal

Practical Lossless Compression with Latent Variables using Bits Back Coding

Authors: James Townsend, Tom Bird, David Barber

Abstract: Deep latent variable models have seen recent success in many data domains. Lossless compression is an application of these models which, despite having the potential to be highly useful, has yet to be implemented in a practical manner. We present `Bits Back with ANS' (BB-ANS), a scheme to perform lossless compression with latent variable models at a near optimal rate. We demonstrate this scheme by… ▽ More Deep latent variable models have seen recent success in many data domains. Lossless compression is an application of these models which, despite having the potential to be highly useful, has yet to be implemented in a practical manner. We present `Bits Back with ANS' (BB-ANS), a scheme to perform lossless compression with latent variable models at a near optimal rate. We demonstrate this scheme by using it to compress the MNIST dataset with a variational auto-encoder model (VAE), achieving compression rates superior to standard methods with only a simple VAE. Given that the scheme is highly amenable to parallelization, we conclude that with a sufficiently high quality generative model this scheme could be used to achieve substantial improvements in compression rate with acceptable running time. We make our implementation available open source at https://github.com/bits-back/bits-back . △ Less

Submitted 15 January, 2019; originally announced January 2019.

Towards a systematically improvable many-body description of antiferromagnetic iron oxide

Authors: Joshua P Townsend, Raymond C Clay III, Thomas R Mattsson, Eric Neuscamman, Luning Zhao, Ken Esler, Ronald E Cohen, Luke Shulenburger

Abstract: We report variational and fixed-node diffusion quantum Monte Carlo (QMC) calculations of anti-ferromagnetic iron oxide (FeO) in the ground state B1 crystal structure. The goal of this study was a systematic investigation of the sensitivity of several ground state properties to a variety of QMC wave function generation techniques including advanced wave functions such as multi-determinant expansion… ▽ More We report variational and fixed-node diffusion quantum Monte Carlo (QMC) calculations of anti-ferromagnetic iron oxide (FeO) in the ground state B1 crystal structure. The goal of this study was a systematic investigation of the sensitivity of several ground state properties to a variety of QMC wave function generation techniques including advanced wave functions such as multi-determinant expansions and backflow transformations. We found that the predicted lattice distortion was largely controlled by the choice of single particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. However, the absolute magnetic moment was remarkably insensitive to the method of wave function construction. QMC estimates of total spin density indicate that in addition to strong electronic correlation of the Fe $3d$ states, charge transfer may be an important but challenging piece of physics to accurately capture within existing QMC methods. Finally, we highlight the need for advanced and systematically improvable many-body wave functions suitable for accurately describing challenging real systems. △ Less

Submitted 6 June, 2018; v1 submitted 4 June, 2018; originally announced June 2018.

Comments: Further calculations are in progress. This version does not accurately represent the results of our investigation

A Theoretical Study of Process Dependence for Critical Statistics in Standard Serial Models and Standard Parallel Models

Authors: Ru Zhang, Yanjun Liu, James T Townsend

Abstract: Critical parts of the definitions of standard serial and standard parallel modes refer to stochastic independence. Standard serial models are defined by stochastic independence and identical distributions of their processing times. Processing times in the serial models are identical to the intercompletion time statistics. Similarly, standard parallel models assume stochastically independent and id… ▽ More Critical parts of the definitions of standard serial and standard parallel modes refer to stochastic independence. Standard serial models are defined by stochastic independence and identical distributions of their processing times. Processing times in the serial models are identical to the intercompletion time statistics. Similarly, standard parallel models assume stochastically independent and identical processing times. Their processing times are equivalent to the statistic known as total completion times. Little is known about what standard serial models can predict for the total completion time or what standard parallel models can predict for the intercompletion times. Here we demonstrate that standard serial models possess a tendency to predict a positive dependence for the total completion times with that always being true in the case of a single processing order. However, with mixtures of processing orders, standard serial models may predict negative dependence of the total completion times. Comparably, standard parallel models typically predict neither independence of the intercompletion times nor identical distributions. In fact, standard parallel models predict increasing intercompletion times as the individual channels continue to finish. Nevertheless, dramatically increasing hazard functions of the channels can defeat that tendency. And, standard parallel models can predict intercompletion time independence but only when individual channel distributions are exponential. Finally, we use these and ancillary mathematical results to conclude that standard serial and standard parallel models can never perfectly mimic one another. Therefore, our findings set the stage for explicit model testing between these classes. △ Less

Submitted 11 February, 2019; v1 submitted 19 March, 2018; originally announced March 2018.

Comments: arXiv admin note: substantial text overlap with arXiv:1712.00528

Sub 20 meV Schottky barriers in metal/MoTe2 junctions

Authors: Nicola J. Townsend, Iddo Amit, Monica F. Craciun, Saverio Russo

Abstract: The newly emerging class of atomically-thin materials has shown a high potential for the realisation of novel electronic and optoelectronic components. Amongst this family, semiconducting transition metal dichalcogenides (TMDCs) are of particular interest. While their band gaps are compatible with those of conventional solid state devices, they present a wide range of exciting new properties that… ▽ More The newly emerging class of atomically-thin materials has shown a high potential for the realisation of novel electronic and optoelectronic components. Amongst this family, semiconducting transition metal dichalcogenides (TMDCs) are of particular interest. While their band gaps are compatible with those of conventional solid state devices, they present a wide range of exciting new properties that is bound to become a crucial ingredient in the future of electronics. To utilise these properties for the prospect of electronics in general, and long-wavelength-based photodetectors in particular, the Schottky barriers formed upon contact with a metal and the contact resistance that arises at these interfaces have to be measured and controlled. We present experimental evidence for the formation of Schottky barriers as low as 10 meV between MoTe2 and metal electrodes. By varying the electrode work functions, we demonstrate that Fermi level pinning due to metal induced gap states at the interfaces occurs at 0.14 eV above the valence band maximum. In this configuration, thermionic emission is observed for the first time at temperatures between 40 K and 75 K. Finally, we discuss the ability to tune the barrier height using a gate electrode. △ Less

Submitted 12 March, 2018; originally announced March 2018.

Comments: 26 pages, 13 figures, 2D Materials (2018)

QMCPACK : An open source ab initio Quantum Monte Carlo package for the electronic structure of atoms, molecules, and solids

Authors: Jeongnim Kim, Andrew Baczewski, Todd D. Beaudet, Anouar Benali, M. Chandler Bennett, Mark A. Berrill, Nick S. Blunt, Edgar Josue Landinez Borda, Michele Casula, David M. Ceperley, Simone Chiesa, Bryan K. Clark, Raymond C. Clay III, Kris T. Delaney, Mark Dewing, Kenneth P. Esler, Hongxia Hao, Olle Heinonen, Paul R. C. Kent, Jaron T. Krogel, Ilkka Kylanpaa, Ying Wai Li, M. Graham Lopez, Ye Luo, Fionn D. Malone , et al. (23 additional authors not shown)

Abstract: QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a s… ▽ More QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit (CPU) and graphical processing unit (GPU) systems. We detail the program's capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://www.qmcpack.org . △ Less

Submitted 4 April, 2018; v1 submitted 19 February, 2018; originally announced February 2018.

Journal ref: J. Phys.: Condens. Matter 30 195901 (2018)

A Theoretical Study of Process Dependence for Standard Two-Process Serial Models and Standard Two-Process Parallel Models

Authors: Ru Zhang, Yanjun Liu, James T. Townsend

Abstract: In this article we differentiate and characterize the standard two-process serial models and the standard two process parallel models by investigating the behavior of (conditional) distributions of the total completion times and survivals of intercompletion times without assuming any particular forms for the distributions of processing times. We address our argument through mathematical proofs and… ▽ More In this article we differentiate and characterize the standard two-process serial models and the standard two process parallel models by investigating the behavior of (conditional) distributions of the total completion times and survivals of intercompletion times without assuming any particular forms for the distributions of processing times. We address our argument through mathematical proofs and computational methods. It is found that for the standard two-process serial models, positive dependence between the total completion times does not hold if no specific distributional forms are imposed to the processing times. By contrast, for the standard two-process parallel models the total completion times are independent. According to different nature of process dependence, one can distinguish a standard two process serial model from a standard two-process parallel model. We also find that in standard two-process parallel models the monotonicity of survival function of the intercompletion time of stage 2 conditional on the completion of stage 1 depends on the monotonicity of the hazard function of processing time. We also find that the survival of intercompletion time(s) from stage 1 to stage 2 is increasing when the ratio of hazard function meets certain criterion. Then the empirical finding that the intercompletion time is grown with the growth of the number of recalled words can be accounted by standard parallel models. We also find that if the cumulative hazard function is concave or linear, the survival from stage 1 to stage 2 is increasing. △ Less

Submitted 1 December, 2017; originally announced December 2017.

Long-term optical and X-ray variability of the Be/X-ray binary H 1145-619: discovery of an on-going retrograde density wave

Authors: J. Alfonso-Garzón, J. Fabregat, P. Reig, J. J. E. Kajava, C. Sánchez-Fernández, L. J. Townsend, J. M. Mas-Hesse, S. M. Crawford, P. Kretschmar, M. J. Coe

Abstract: Multiwavelength monitoring of Be/X-ray binaries is crucial to understand the mechanisms producing their outbursts. H 1145-619 is one of these systems, which has recently displayed X-ray activity. We investigate the correlation between the optical emission and the X-ray activity to predict the occurrence of new X-ray outbursts from the inferred state of the circumstellar disc. We have performed a m… ▽ More Multiwavelength monitoring of Be/X-ray binaries is crucial to understand the mechanisms producing their outbursts. H 1145-619 is one of these systems, which has recently displayed X-ray activity. We investigate the correlation between the optical emission and the X-ray activity to predict the occurrence of new X-ray outbursts from the inferred state of the circumstellar disc. We have performed a multiwavelength study of H 1145-619 from 1973 to 2017 and present here a global analysis of its variability over the last 40 years. We have used optical spectra from the SAAO, SMARTS and SALT telescopes and optical photometry from INTEGRAL/OMC and ASAS. We also used X-ray observations from INTEGRAL/JEM-X, and IBIS to generate the light curves and combined them with Swift/XRT to extract the X-ray spectra. In addition, we have compiled archival observations and measurements from the literature to complement these data. Comparing the evolution of the optical continuum emission with the Halpha line variability, we have identified three different patterns of optical variability: Global increases and decreases of the optical brightness (observed from 1982 to 1994 and from 2009 to 2017) that can be explained by the dissipation and replenishment of the circumstellar disc, superorbital variations with a period of Psuperorb~590 days (observed in 2002-2009) which seems to be related with the circumstellar disc, and optical outbursts (observed in 1998-1999 and 2002-2005) that we interpret as mass ejections. We have discovered the presence of a retrograde one-armed density wave, which appeared in 2016 and is still present in the circumstellar disc. We have carried out the most complete long-term optical study of the Be/X-ray binary H 1145-619 in correlation with its X-ray activity. We found for the fist time the presence of a retrograde density perturbation in the circumstellar disc of a Be/X-ray binary. △ Less

Submitted 20 September, 2017; v1 submitted 4 September, 2017; originally announced September 2017.

Comments: Accepted on 2 September 2017, 13 pages, 12 figures

Journal ref: A&A 607, A52 (2017)

Role of Charge Traps in the Performance of Atomically-Thin Transistors

Authors: Iddo Amit, Tobias J. Octon, Nicola J. Townsend, Francesco Reale, C. David Wright, Cecilia Mattevi, Monica F. Craciun, Saverio Russo

Abstract: Transient currents in atomically thin MoTe$_2$ field-effect transistor are measured during cycles of pulses through the gate electrode. The transients are analyzed in light of a newly proposed model for charge trapping dynamics that renders a time-dependent change in threshold voltage the dominant effect on the channel hysteretic behavior over emission currents from the charge traps. The proposed… ▽ More Transient currents in atomically thin MoTe$_2$ field-effect transistor are measured during cycles of pulses through the gate electrode. The transients are analyzed in light of a newly proposed model for charge trapping dynamics that renders a time-dependent change in threshold voltage the dominant effect on the channel hysteretic behavior over emission currents from the charge traps. The proposed model is expected to be instrumental in understanding the fundamental physics that governs the performance of atomically thin FETs and is applicable to the entire class of atomically thin-based devices. Hence, the model is vital to the intelligent design of fast and highly efficient opto-electronic devices. △ Less

Submitted 16 March, 2017; originally announced March 2017.

Comments: 25 pages, 7 figures, including the supporting information

Journal ref: Advanced Materials Early View (2017), 1605598

The 2016 super-Eddington outburst of SMC X-3: X-ray and optical properties and system parameters

Authors: L. J. Townsend, J. A. Kennea, M. J. Coe, V. A. McBride, D. A. H. Buckley, P. A. Evans, A. Udalski

Abstract: On 2016 July 30 (MJD 57599), observations of the Small Magellanic Cloud by Swift/XRT found an increase in X-ray counts coming from a position consistent with the Be/X-ray binary pulsar SMC X-3. Follow-up observations on 2016 August 3 (MJD 57603) and 2016 August 10 (MJD 57610) revealed a rapidly increasing count rate and confirmed the onset of a new X-ray outburst from the system. Further monitorin… ▽ More On 2016 July 30 (MJD 57599), observations of the Small Magellanic Cloud by Swift/XRT found an increase in X-ray counts coming from a position consistent with the Be/X-ray binary pulsar SMC X-3. Follow-up observations on 2016 August 3 (MJD 57603) and 2016 August 10 (MJD 57610) revealed a rapidly increasing count rate and confirmed the onset of a new X-ray outburst from the system. Further monitoring by Swift began to uncover the enormity of the outburst, which peaked at 1.2 x 10^39 erg/s on 2016 August 25 (MJD 57625). The system then began a gradual decline in flux that was still continuing over 5 months after the initial detection. We explore the X-ray and optical behaviour of SMC X-3 between 2016 July 30 and 2016 December 18 during this super-Eddington outburst. We apply a binary model to the spin-period evolution that takes into account the complex accretion changes over the outburst, to solve for the orbital parameters. Our results show SMC X-3 to be a system with a moderately low eccentricity amongst the Be/X-ray binary systems and to have a dynamically determined orbital period statistically consistent with the prominent period measured in the OGLE optical light curve. Our optical and X-ray derived ephemerides show that the peak in optical flux occurs roughly 6 days after periastron. The measured increase in I-band flux from the counterpart during the outburst is reflected in the measured equivalent width of the H-alpha line emission, though the H-alpha emission itself seems variable on sub-day time-scales, possibly due to the NS interacting with an inhomogeneous disc. △ Less

Submitted 21 July, 2017; v1 submitted 9 January, 2017; originally announced January 2017.

Comments: 10 pages, 8 figures, accepted for publication in MNRAS. This article replaces a previously submitted version

A new foreperiod effect on single-trial phase coherence. Part I: existence and relevance

Authors: Joaquin Rapela, Marissa Westerfield, Jeanne Townsend, Scott Makeig

Abstract: Expecting events in time leads to more efficient behavior. A remarkable early finding in the study of temporal expectancy is the foreperiod effect on reaction times; i.e., the influence or reaction time of the time period between a warning signal and an imperative stimulus to which subjects are instructed to respond as quickly as possible. Recently it has been shown that the phase of oscillatory a… ▽ More Expecting events in time leads to more efficient behavior. A remarkable early finding in the study of temporal expectancy is the foreperiod effect on reaction times; i.e., the influence or reaction time of the time period between a warning signal and an imperative stimulus to which subjects are instructed to respond as quickly as possible. Recently it has been shown that the phase of oscillatory activity preceding stimulus presentation is related to behavior. Here we connect both of these findings by reporting a novel foreperiod effect on the inter-trial phase coherence of the electroencephalogram (EEG) triggered by stimuli to which subjects are instructed not to respond. Inter-trial phase coherence has been used to describe regularities in phases of groups of trials time locked to an event of interest. We propose a single-trial measure of inter-trial phase coherence and prove its soundness. Equipped with this measure, and using a multivariate decoding method, we demonstrate that the foreperiod duration in and audiovisual attention-shifting task modulates single-trial phase coherence. In principle, this modulation could be an artifact of the decoding method used to detect it. We show that this is not the case, since the modulation can also be observed using a simple averaging method. We show that the strength of this modulation correlates with subject behavior (both error rates and mean-reaction times). We anticipate that the new foreperiod effect on inter-trial phase coherence, and the decoding method used here to detect it, will be important tools to understand cognition at the single-trial level. In Part II of this manuscript, we support this claim, by showing that changes in attention modulate the strength of the new foreperiod effect on a trial-by-trial basis. △ Less

Submitted 26 August, 2016; originally announced November 2016.

First multicolour polarimetry of TeV gamma-ray binary HESS J0632+057 close to periastron passage

Authors: R. V. Yudin, S. B. Potter, L. J. Townsend

Abstract: We present the results of UBVRI polarimetry of the TeV gamma-ray binary HESS J0632+057 obtained on 2015 March 24 (JD 2457106) and 2015 December 12 (JD 2457369). The detected polarisation values of HESS J0632+057, just after periastron passage (March 24), are higher than all previously published values (p(V) ~ 4.2%) and the position angle (theta ~ 171-172 degrees) is also different by ~6-10 degrees… ▽ More We present the results of UBVRI polarimetry of the TeV gamma-ray binary HESS J0632+057 obtained on 2015 March 24 (JD 2457106) and 2015 December 12 (JD 2457369). The detected polarisation values of HESS J0632+057, just after periastron passage (March 24), are higher than all previously published values (p(V) ~ 4.2%) and the position angle (theta ~ 171-172 degrees) is also different by ~6-10 degrees from previously published values. The data obtained just before the subsequent periastron passage (December 12) show statistically lower polarisation in all photometric bands (p(V) ~ 3.9%) and a different position angle theta ~ 167-168 degrees. From observations of a nearby field star, the interstellar component of the measured polarisation was estimated as p_is(V) ~ 0.65% and theta_is ~ 153 degrees. This estimate was used with the previous "V"-band estimation by "field-stars method" (p(V) ~ 2% and theta_is ~ 165 degrees) of Yudin (2014) to identify the wavelength dependence of the intrinsic polarisation in HESS J0632+057. It was found that after subtraction of the interstellar component (for both p_is estimates), the wavelength dependence of the intrinsic polarisation in HESS J0632+057 is essentially flat. We propose that the formation of an additional source of polarisation or some perturbation of circumstellar material at this orbital phase can explain the changes in the level of polarisation in HESS J0632+057 close to periastron passage. △ Less

Submitted 15 October, 2016; originally announced October 2016.

Comments: 7 pages, 4 figures, Accepted for publication in MNRAS on 2016 October 14

Orbital and Superorbital Monitoring of the Be/X-ray binary A0538-66: constraints on the system parameters

Authors: A. F. Rajoelimanana, P. A. Charles, P. J. Meintjes, L. J. Townsend, M. P. E. Schurch, A. Udalski

Abstract: We combine the decade long photometry of the Be/X-ray binary system A0538-66 provided by the MACHO and OGLE IV projects with high resolution SALT spectroscopy to provide detailed constraints on the orbital parameters and system properties. The ~420d superorbital modulation is present throughout, but has reduced in amplitude in recent years. The well-defined 16.6409d orbital outbursts, which were a… ▽ More We combine the decade long photometry of the Be/X-ray binary system A0538-66 provided by the MACHO and OGLE IV projects with high resolution SALT spectroscopy to provide detailed constraints on the orbital parameters and system properties. The ~420d superorbital modulation is present throughout, but has reduced in amplitude in recent years. The well-defined 16.6409d orbital outbursts, which were a strong function of superorbital phase in the MACHO data (not occurring at all at superorbital maximum), are present throughout the OGLE IV coverage. However, their amplitude reduces during superorbital maximum. We have refined the orbital period and ephemeris of the optical outburst based on ~25 yrs light curves to HJD = 2455674.48 +/- 0.03 + n*16.6409 +/- 0.0003d. Our SALT spectra reveal a B1 III star with vsini of 285 km/s from which we have derived an orbital radial velocity curve which confirms the high eccentricity of e = 0.72 +/- 0.14. Furthermore, the mass function indicates that, unless the neutron star far exceeds the canonical 1.44 Msun, the donor must be significantly undermassive for its spectral type. We discuss the implications of the geometry and our derived orbital solution on the observed behaviour of the system. △ Less

Submitted 30 September, 2016; originally announced September 2016.

Comments: 13 pages, 15 figures, 2 tables, accepted for publication in MNRAS

An integrated and multi-purpose microscope for the characterization of atomically thin optoelectronic devices

Authors: Adolfo De Sanctis, Gareth F. Jones, Nicola J. Townsend, Monica F. Craciun, Saverio Russo

Abstract: Optoelectronic devices based on graphene and other two-dimensional (2D) materials, such as transition metal dichalcogenides (TMDs) are the focus of wide research interest. The characterization these emerging atomically thin materials and devices strongly relies on a set of measurements involving both optical and electronic instrumentation ranging from scanning photocurrent mapping to Raman and pho… ▽ More Optoelectronic devices based on graphene and other two-dimensional (2D) materials, such as transition metal dichalcogenides (TMDs) are the focus of wide research interest. The characterization these emerging atomically thin materials and devices strongly relies on a set of measurements involving both optical and electronic instrumentation ranging from scanning photocurrent mapping to Raman and photoluminescence (PL) spectroscopy. Furthermore, proof-of-concept devices are usually fabricated from micro-meter size flakes, requiring microscopy techniques to characterize them. Current state-of-the-art commercial instruments offer the ability to characterize individual properties of these materials with no option for the in situ characterization of a wide enough range of complementary optical and electrical properties. Presently, the requirement to switch atomically-thin materials from one system to another often radically affects the properties of these uniquely sensitive materials through atmospheric contamination. Here, we present an integrated, multi-purpose instrument dedicated to the optical and electrical characterization of devices based on 2D materials which is able to perform low frequency electrical measurements, scanning photocurrent mapping, Raman, absorption and PL spectroscopy in one single set-up. We characterize this apparatus by performing multiple measurements on graphene, transition metal dichalcogenides (TMDs) and Si. The performance and resolution of each individual measurement technique is found to be equivalent to that of commercially available instruments. Contrary to nowadays commercial systems, a significant advantage of the developed instrument is that for the first time the integration of a wide range of complementary opto-electronic and spectroscopy characterization techniques is demonstrated in a single compact unit. △ Less

Submitted 26 April, 2017; v1 submitted 21 September, 2016; originally announced September 2016.

Comments: To appear in Review of Scientific Instruments (2017). 29 pages, 9 figures

Journal ref: Review of Scientific Instruments 88, 055102 (2017)

A Luminous Gamma-ray Binary in the Large Magellanic Cloud

Authors: R. H. D. Corbet, L. Chomiuk, M. J. Coe, J. B. Coley, G. Dubus, P. G. Edwards, P. Martin, V. A. McBride, J. Stevens, J. Strader, L. J. Townsend, A. Udalski

Abstract: Gamma-ray binaries consist of a neutron star or a black hole interacting with a normal star to produce gamma-ray emission that dominates the radiative output of the system. Only a handful of such systems have been previously discovered, all within our Galaxy. Here we report the discovery with the Fermi Large Area Telescope (LAT) of a luminous gamma-ray binary in the Large Magellanic Cloud from a s… ▽ More Gamma-ray binaries consist of a neutron star or a black hole interacting with a normal star to produce gamma-ray emission that dominates the radiative output of the system. Only a handful of such systems have been previously discovered, all within our Galaxy. Here we report the discovery with the Fermi Large Area Telescope (LAT) of a luminous gamma-ray binary in the Large Magellanic Cloud from a search for periodic modulation in all sources in the third Fermi LAT catalog. This is the first such system to be found outside the Milky Way. The system has an orbital period of 10.3 days and is associated with a massive O5III star located in the supernova remnant DEM L241, previously identified as the candidate high-mass X-ray binary (HMXB) CXOU J053600.0-673507. X-ray and radio emission are also modulated on the 10.3 day period, but are in anti-phase with the gamma-ray modulation. Optical radial velocity measurements suggest that the system contains a neutron star. The source is significantly more luminous than similar sources in the Milky Way at radio, optical, X-ray and gamma-ray wavelengths. The detection of this extra-galactic system, but no new Galactic systems raises the possibility that the predicted number of gamma-ray binaries in our Galaxy has been overestimated, and that HMXBs may be born containing relatively slowly rotating neutron stars. △ Less

Submitted 23 August, 2016; originally announced August 2016.

Comments: 11 pages, 6 figures. Accepted for publication in The Astrophysical Journal

2S1553-542: a Be/X-ray binary pulsar on the far side of the Galaxy

Authors: Alexander A. Lutovinov, David A. H. Buckley, Lee J. Townsend, Sergey S. Tsygankov, Jamie Kennea

Abstract: We report the results of a comprehensive analysis of X-ray (Chandra and Swift observatories), optical (Southern African Large Telescope, SALT) and near-infrared (the VVV survey) observations of the Be/X-ray binary pulsar 2S1553-542. Accurate coordinates for the X-ray source are determined and are used to identify the faint optical/infrared counterpart for the first time. Using VVV and SALTICAM pho… ▽ More We report the results of a comprehensive analysis of X-ray (Chandra and Swift observatories), optical (Southern African Large Telescope, SALT) and near-infrared (the VVV survey) observations of the Be/X-ray binary pulsar 2S1553-542. Accurate coordinates for the X-ray source are determined and are used to identify the faint optical/infrared counterpart for the first time. Using VVV and SALTICAM photometry, we have constructed the spectral energy distribution (SED) for this star and found a moderate NIR excess that is typical for Be stars and arises due to the presence of circumstellar material (disk). A comparison of the SED with those of known Be/X-ray binaries has allowed us to estimate the spectral type of the companion star as B1-2V and the distance to the system as $>15$ kpc. This distance estimation is supported by the X-ray data and makes 2S1553-542 one of the most distant X-ray binaries within the Milky Way, residing on the far side in the Scutum-Centaurus arm or even further. △ Less

Submitted 27 July, 2016; originally announced July 2016.

Comments: 8 pages, 4 figures, accepted for publication in MNRAS

Pymanopt: A Python Toolbox for Optimization on Manifolds using Automatic Differentiation

Authors: James Townsend, Niklas Koep, Sebastian Weichwald

Abstract: Optimization on manifolds is a class of methods for optimization of an objective function, subject to constraints which are smooth, in the sense that the set of points which satisfy the constraints admits the structure of a differentiable manifold. While many optimization problems are of the described form, technicalities of differential geometry and the laborious calculation of derivatives pose a… ▽ More Optimization on manifolds is a class of methods for optimization of an objective function, subject to constraints which are smooth, in the sense that the set of points which satisfy the constraints admits the structure of a differentiable manifold. While many optimization problems are of the described form, technicalities of differential geometry and the laborious calculation of derivatives pose a significant barrier for experimenting with these methods. We introduce Pymanopt (available at https://pymanopt.github.io), a toolbox for optimization on manifolds, implemented in Python, that---similarly to the Manopt Matlab toolbox---implements several manifold geometries and optimization algorithms. Moreover, we lower the barriers to users further by using automated differentiation for calculating derivative information, saving users time and saving them from potential calculation and implementation errors. △ Less

Submitted 8 September, 2016; v1 submitted 10 March, 2016; originally announced March 2016.

Journal ref: Journal of Machine Learning Research, 17(137):1-5, 2016 ( https://jmlr.org/papers/v17/16-177.html )

The VC-Dimension of Similarity Hypotheses Spaces

Authors: Mark Herbster, Paul Rubenstein, James Townsend

Abstract: Given a set $X$ and a function $h:X\longrightarrow\{0,1\}$ which labels each element of $X$ with either $0$ or $1$, we may define a function $h^{(s)}$ to measure the similarity of pairs of points in $X$ according to $h$. Specifically, for $h\in \{0,1\}^X$ we define $h^{(s)}\in \{0,1\}^{X\times X}$ by $h^{(s)}(w,x):= \mathbb{1}[h(w) = h(x)]$. This idea can be extended to a set of functions, or hypo… ▽ More Given a set $X$ and a function $h:X\longrightarrow\{0,1\}$ which labels each element of $X$ with either $0$ or $1$, we may define a function $h^{(s)}$ to measure the similarity of pairs of points in $X$ according to $h$. Specifically, for $h\in \{0,1\}^X$ we define $h^{(s)}\in \{0,1\}^{X\times X}$ by $h^{(s)}(w,x):= \mathbb{1}[h(w) = h(x)]$. This idea can be extended to a set of functions, or hypothesis space $\mathcal{H} \subseteq \{0,1\}^X$ by defining a similarity hypothesis space $\mathcal{H}^{(s)}:=\{h^{(s)}:h\in\mathcal{H}\}$. We show that ${vc-dimension}(\mathcal{H}^{(s)}) \in Θ({vc-dimension}(\mathcal{H}))$. △ Less

Submitted 25 February, 2015; originally announced February 2015.

Comments: 6 pages

SXP 5.05 = IGR J00569-7226 : using X-rays to explore the structure of a Be stars circumstellar disk

Authors: M. J. Coe, E. S. Bartlett, A. J. Bird, F. Haberl, J. A. Kennea, V. A. McBride, L. J. Townsend, A. Udalski

Abstract: On MJD 56590-1 (2013 Oct 25-26) observations of the Magellanic Clouds by the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) observatory discovered a previously-unreported bright, flaring X-ray source. This source was initially given the identification IGR J00569-7226. Subsequent multi-wavelength observations identified the system as new Be/X-ray binary system in the Small Magellanic Cl… ▽ More On MJD 56590-1 (2013 Oct 25-26) observations of the Magellanic Clouds by the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) observatory discovered a previously-unreported bright, flaring X-ray source. This source was initially given the identification IGR J00569-7226. Subsequent multi-wavelength observations identified the system as new Be/X-ray binary system in the Small Magellanic Cloud. Follow-up X-ray observations by Swift and XMM-Newton revealed an X-ray pulse period of 5.05s and that the system underwent regular occulation/eclipse behaviour every 17d. This is the first reported eclipsing Be/X-ray binary system in the SMC, and only the second such system known to date. Furthermore, the nature of the occultation makes it possible to use the neutron star to X-ray the circumstellar disk, thereby, for the first time, revealing direct observational evidence for its size and clumpy structure. Swift timing measurements allowed for the binary solution to be calculated from the Doppler shifted X-ray pulsations. This solution suggests this is a low eccentricity binary relative to others measured in the SMC. Finally it is interesting to note that the mass determined from this dynamical method for the Be star (approx 13 solar masses) is significantly different from that inferred from the spectroscopic classification of B0.2Ve (approx 16 solar masses) - an effect that has been noted for some other high mass X-ray binary (HMXB) systems. △ Less

Submitted 10 December, 2014; originally announced December 2014.

Comments: 19 pages, accepted by MNRAS