Showing 1–50 of 80 results for author: Mitchell, B

Error mitigation with stabilized noise in superconducting quantum processors

Authors: Youngseok Kim, Luke C. G. Govia, Andrew Dane, Ewout van den Berg, David M. Zajac, Bradley Mitchell, Yinyu Liu, Karthik Balakrishnan, George Keefe, Adam Stabile, Emily Pritchett, Jiri Stehlik, Abhinav Kandala

Abstract: Pre-fault tolerant quantum computers have already demonstrated the ability to estimate observable values accurately, at a scale beyond brute-force classical computation. This has been enabled by error mitigation techniques that often rely on a representative model on the device noise. However, learning and maintaining these models is complicated by fluctuations in the noise over unpredictable time… ▽ More Pre-fault tolerant quantum computers have already demonstrated the ability to estimate observable values accurately, at a scale beyond brute-force classical computation. This has been enabled by error mitigation techniques that often rely on a representative model on the device noise. However, learning and maintaining these models is complicated by fluctuations in the noise over unpredictable time scales, for instance, arising from resonant interactions between superconducting qubits and defect two-level systems (TLS). Such interactions affect the stability and uniformity of device performance as a whole, but also affect the noise model accuracy, leading to incorrect observable estimation. Here, we experimentally demonstrate that tuning of the qubit-TLS interactions helps reduce noise instabilities and consequently enables more reliable error-mitigation performance. These experiments provide a controlled platform for studying the performance of error mitigation in the presence of quasi-static noise. We anticipate that the capabilities introduced here will be crucial for the exploration of quantum applications on solid-state processors at non-trivial scales. △ Less

Submitted 5 July, 2024; v1 submitted 2 July, 2024; originally announced July 2024.

Comments: 8 pages, 4 figures (13 pages, 8 figures for supplementary material), reference numbering has been fixed at v2

Dispersive Vertex Guarding for Simple and Non-Simple Polygons

Authors: Sándor P. Fekete, Joseph S. B. Mitchell, Christian Rieck, Christian Scheffer, Christiane Schmidt

Abstract: We study the Dispersive Art Gallery Problem with vertex guards: Given a polygon $\mathcal{P}$, with pairwise geodesic Euclidean vertex distance of at least $1$, and a rational number $\ell$; decide whether there is a set of vertex guards such that $\mathcal{P}$ is guarded, and the minimum geodesic Euclidean distance between any two guards (the so-called dispersion distance) is at least $\ell$. W… ▽ More We study the Dispersive Art Gallery Problem with vertex guards: Given a polygon $\mathcal{P}$, with pairwise geodesic Euclidean vertex distance of at least $1$, and a rational number $\ell$; decide whether there is a set of vertex guards such that $\mathcal{P}$ is guarded, and the minimum geodesic Euclidean distance between any two guards (the so-called dispersion distance) is at least $\ell$. We show that it is NP-complete to decide whether a polygon with holes has a set of vertex guards with dispersion distance $2$. On the other hand, we provide an algorithm that places vertex guards in simple polygons at dispersion distance at least $2$. This result is tight, as there are simple polygons in which any vertex guard set has a dispersion distance of at most $2$. △ Less

Submitted 9 June, 2024; originally announced June 2024.

Comments: 13 pages, 14 figures; accepted at the 36th Canadian Conference on Computational Geometry (CCCG 2024)

ACM Class: F.2.2

Multirobot Watchman Routes in a Simple Polygon

Authors: Joseph S. B. Mitchell, Linh Nguyen

Abstract: The well-known \textsc{Watchman Route} problem seeks a shortest route in a polygonal domain from which every point of the domain can be seen. In this paper, we study the cooperative variant of the problem, namely the \textsc{$k$-Watchmen Routes} problem, in a simple polygon $P$. We look at both the version in which the $k$ watchmen must collectively see all of $P$, and the quota version in which t… ▽ More The well-known \textsc{Watchman Route} problem seeks a shortest route in a polygonal domain from which every point of the domain can be seen. In this paper, we study the cooperative variant of the problem, namely the \textsc{$k$-Watchmen Routes} problem, in a simple polygon $P$. We look at both the version in which the $k$ watchmen must collectively see all of $P$, and the quota version in which they must see a predetermined fraction of $P$'s area. We give an exact pseudopolynomial time algorithm for the \textsc{$k$-Watchmen Routes} problem in a simple orthogonal polygon $P$ with the constraint that watchmen must move on axis-parallel segments, and there is a given common starting point on the boundary. Further, we give a fully polynomial-time approximation scheme and a constant-factor approximation for unconstrained movement. For the quota version, we give a constant-factor approximation in a simple polygon, utilizing the solution to the (single) \textsc{Quota Watchman Route} problem. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Journal ref: 36th Canadian Conference on Computational Geometry (CCCG 2024)

Theory of quantum error mitigation for non-Clifford gates

Authors: David Layden, Bradley Mitchell, Karthik Siva

Abstract: Quantum error mitigation techniques mimic noiseless quantum circuits by running several related noisy circuits and combining their outputs in particular ways. How well such techniques work is thought to depend strongly on how noisy the underlying gates are. Weakly-entangling gates, like $R_{ZZ}(θしーた)$ for small angles $θしーた$, can be much less noisy than entangling Clifford gates, like CNOT and CZ, and t… ▽ More Quantum error mitigation techniques mimic noiseless quantum circuits by running several related noisy circuits and combining their outputs in particular ways. How well such techniques work is thought to depend strongly on how noisy the underlying gates are. Weakly-entangling gates, like $R_{ZZ}(θしーた)$ for small angles $θしーた$, can be much less noisy than entangling Clifford gates, like CNOT and CZ, and they arise naturally in circuits used to simulate quantum dynamics. However, such weakly-entangling gates are non-Clifford, and are therefore incompatible with two of the most prominent error mitigation techniques to date: probabilistic error cancellation (PEC) and the related form of zero-noise extrapolation (ZNE). This paper generalizes these techniques to non-Clifford gates, and comprises two complementary parts. The first part shows how to effectively transform any given quantum channel into (almost) any desired channel, at the cost of a sampling overhead, by adding random Pauli gates and processing the measurement outcomes. This enables us to cancel or properly amplify noise in non-Clifford gates, provided we can first characterize such gates in detail. The second part therefore introduces techniques to do so for noisy $R_{ZZ}(θしーた)$ gates. These techniques are robust to state preparation and measurement (SPAM) errors, and exhibit concentration and sensitivity--crucial features in many experiments. They are related to randomized benchmarking, and may also be of interest beyond the context of error mitigation. We find that while non-Clifford gates can be less noisy than related Cliffords, their noise is fundamentally more complex, which can lead to surprising and sometimes unwanted effects in error mitigation. Whether this trade-off can be broadly advantageous remains to be seen. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Robustly Guarding Polygons

Authors: Rathish Das, Omrit Filtser, Matthew J. Katz, Joseph S. B. Mitchell

Abstract: We propose precise notions of what it means to guard a domain "robustly", under a variety of models. While approximation algorithms for minimizing the number of (precise) point guards in a polygon is a notoriously challenging area of investigation, we show that imposing various degrees of robustness on the notion of visibility coverage leads to a more tractable (and realistic) problem for which we… ▽ More We propose precise notions of what it means to guard a domain "robustly", under a variety of models. While approximation algorithms for minimizing the number of (precise) point guards in a polygon is a notoriously challenging area of investigation, we show that imposing various degrees of robustness on the notion of visibility coverage leads to a more tractable (and realistic) problem for which we can provide approximation algorithms with constant factor guarantees. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: To appear in SoCG 2024

Optimizing Visibility-based Search in Polygonal Domains

Authors: Kien C. Huynh, Joseph S. B. Mitchell, Linh Nguyen, Valentin Polishchuk

Abstract: Given a geometric domain $P$, visibility-based search problems seek routes for one or more mobile agents ("watchmen") to move within $P$ in order to be able to see a portion (or all) of $P$, while optimizing objectives, such as the length(s) of the route(s), the size (e.g., area or volume) of the portion seen, the probability of detecting a target distributed within $P$ according to a prior distri… ▽ More Given a geometric domain $P$, visibility-based search problems seek routes for one or more mobile agents ("watchmen") to move within $P$ in order to be able to see a portion (or all) of $P$, while optimizing objectives, such as the length(s) of the route(s), the size (e.g., area or volume) of the portion seen, the probability of detecting a target distributed within $P$ according to a prior distribution, etc. The classic watchman route problem seeks a shortest route for an observer, with omnidirectional vision, to see all of $P$. In this paper we study bicriteria optimization problems for a single mobile agent within a polygonal domain $P$ in the plane, with the criteria of route length and area seen. Specifically, we address the problem of computing a minimum length route that sees at least a specified area of $P$ (minimum length, for a given area quota). We also study the problem of computing a length-constrained route that sees as much area as possible. We provide hardness results and approximation algorithms. In particular, for a simple polygon $P$ we provide the first fully polynomial-time approximation scheme for the problem of computing a shortest route seeing an area quota, as well as a (slightly more efficient) polynomial dual approximation. We also consider polygonal domains $P$ (with holes) and the special case of a planar domain consisting of a union of lines. Our results yield the first approximation algorithms for computing a time-optimal search route in $P$ to guarantee some specified probability of detection of a static target within $P$, randomly distributed in $P$ according to a given prior distribution. △ Less

Submitted 18 April, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

ACM Class: F.2.2; I.3.5

Journal ref: 19th Scandinavian Symposium on Algorithm Theory, SWAT 2024

Cloud Native Software Engineering

Authors: Brian S. Mitchell

Abstract: Cloud compute adoption has been growing since its inception in the early 2000's with estimates that the size of this market in terms of worldwide spend will increase from \$700 billion in 2021 to \$1.3 trillion in 2025. While there is a significant research activity in many areas of cloud computing technologies, we see little attention being paid to advancing software engineering practices needed… ▽ More Cloud compute adoption has been growing since its inception in the early 2000's with estimates that the size of this market in terms of worldwide spend will increase from \$700 billion in 2021 to \$1.3 trillion in 2025. While there is a significant research activity in many areas of cloud computing technologies, we see little attention being paid to advancing software engineering practices needed to support the current and next generation of cloud native applications. By cloud native, we mean software that is designed and built specifically for deployment to a modern cloud platform. This paper frames the landscape of Cloud Native Software Engineering from a practitioners standpoint, and identifies several software engineering research opportunities that should be investigated. We cover specific engineering challenges associated with software architectures commonly used in cloud applications along with incremental challenges that are expected with emerging IoT/Edge computing use cases. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Geometric Spanning Trees Minimizing the Wiener Index

Authors: A. Karim Abu-Affash, Paz Carmi, Ori Luwisch, Joseph S. B. Mitchell

Abstract: The Wiener index of a network, introduced by the chemist Harry Wiener, is the sum of distances between all pairs of nodes in the network. This index, originally used in chemical graph representations of the non-hydrogen atoms of a molecule, is considered to be a fundamental and useful network descriptor. We study the problem of constructing geometric networks on point sets in Euclidean space that… ▽ More The Wiener index of a network, introduced by the chemist Harry Wiener, is the sum of distances between all pairs of nodes in the network. This index, originally used in chemical graph representations of the non-hydrogen atoms of a molecule, is considered to be a fundamental and useful network descriptor. We study the problem of constructing geometric networks on point sets in Euclidean space that minimize the Wiener index: given a set $P$ of $n$ points in $\mathbb{R}^d$, the goal is to construct a network, spanning $P$ and satisfying certain constraints, that minimizes the Wiener index among the allowable class of spanning networks. In this work, we focus mainly on spanning networks that are trees and we focus on problems in the plane ($d=2$). We show that any spanning tree that minimizes the Wiener index has non-crossing edges in the plane. Then, we use this fact to devise an $O(n^4)$-time algorithm that constructs a spanning tree of minimum Wiener index for points in convex position. We also prove that the problem of computing a spanning tree on $P$ whose Wiener index is at most $W$, while having total (Euclidean) weight at most $B$, is NP-hard. Computing a tree that minimizes the Wiener index has been studied in the area of communication networks, where it is known as the optimum communication spanning tree problem. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: 14 pages, 7 figures

Minimum-link $C$-Oriented Paths Visiting a Sequence of Regions in the Plane

Authors: Kerem Geva, Matthew J. Katz, Joseph S. B. Mitchell, Eli Packer

Abstract: Let $E=\{e_1,\ldots,e_n\}$ be a set of $C$-oriented disjoint segments in the plane, where $C$ is a given finite set of orientations that spans the plane, and let $s$ and $t$ be two points. %(We also require that for each orientation in $C$, its opposite orientation is also in $C$.) We seek a minimum-link $C$-oriented tour of $E$, that is, a polygonal path $πぱい$ from $s$ to $t$ that visits the segmen… ▽ More Let $E=\{e_1,\ldots,e_n\}$ be a set of $C$-oriented disjoint segments in the plane, where $C$ is a given finite set of orientations that spans the plane, and let $s$ and $t$ be two points. %(We also require that for each orientation in $C$, its opposite orientation is also in $C$.) We seek a minimum-link $C$-oriented tour of $E$, that is, a polygonal path $πぱい$ from $s$ to $t$ that visits the segments of $E$ in order, such that, the orientations of its edges are in $C$ and their number is minimum. We present an algorithm for computing such a tour in $O(|C|^2 \cdot n^2)$ time. This problem already captures most of the difficulties occurring in the study of the more general problem, in which $E$ is a set of not-necessarily-disjoint $C$-oriented polygons. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: Full version of paper to appear, CIAC 2023

Integrated photoelasticity in a soft material: phase retardation, azimuthal angle and stress-optic coefficient

Authors: Yuto Yokoyama, Benjamin R. Mitchell, Ali Nassiri, Brad L. Kinsey, Yannis P. Korkolis, Yoshiyuki Tagawa

Abstract: Integrated photoelasticity is investigated for a soft material subjected to a three-dimensional stress state. In the experiment, a solid sphere is pressed against a gelatin gel (Young's modulus is about 4.2 kPa) that deforms up to 4.5 mm depending on the loading forces. The resulting photoelastic parameters (phase retardation, azimuthal angle, and stress-optic coefficient) in the gel are measured… ▽ More Integrated photoelasticity is investigated for a soft material subjected to a three-dimensional stress state. In the experiment, a solid sphere is pressed against a gelatin gel (Young's modulus is about 4.2 kPa) that deforms up to 4.5 mm depending on the loading forces. The resulting photoelastic parameters (phase retardation, azimuthal angle, and stress-optic coefficient) in the gel are measured using a polarization camera. The measured retardation and azimuth are compared with the analytical prediction based on Hertzian contact theory. Remarkably, experimental and analytical results of the photoelastic parameters show a reasonable agreement not only in the retardation but also in the azimuth that is related to the direction of principal stresses and but rarely validated in previous studies, is essential for reconstructing three-dimensional stress fields in soft materials. The stress-optic coefficient of the gelatin gel used is 3.12$\times10^{-8}$ 1/Pa. Such findings proved that integrated photoelasticity is useful for measuring the three-dimensional stress field in soft materials, which is of importance in biomedical engineering and cell printing applications. △ Less

Submitted 8 July, 2022; originally announced August 2022.

Comments: 15 pages, 13 figures

Demonstrating scalable randomized benchmarking of universal gate sets

Authors: Jordan Hines, Marie Lu, Ravi K. Naik, Akel Hashim, Jean-Loup Ville, Brad Mitchell, John Mark Kriekebaum, David I. Santiago, Stefan Seritan, Erik Nielsen, Robin Blume-Kohout, Kevin Young, Irfan Siddiqi, Birgitta Whaley, Timothy Proctor

Abstract: Randomized benchmarking (RB) protocols are the most widely used methods for assessing the performance of quantum gates. However, the existing RB methods either do not scale to many qubits or cannot benchmark a universal gate set. Here, we introduce and demonstrate a technique for scalable RB of many universal and continuously parameterized gate sets, using a class of circuits called randomized mir… ▽ More Randomized benchmarking (RB) protocols are the most widely used methods for assessing the performance of quantum gates. However, the existing RB methods either do not scale to many qubits or cannot benchmark a universal gate set. Here, we introduce and demonstrate a technique for scalable RB of many universal and continuously parameterized gate sets, using a class of circuits called randomized mirror circuits. Our technique can be applied to a gate set containing an entangling Clifford gate and the set of arbitrary single-qubit gates, as well as gate sets containing controlled rotations about the Pauli axes. We use our technique to benchmark universal gate sets on four qubits of the Advanced Quantum Testbed, including a gate set containing a controlled-S gate and its inverse, and we investigate how the observed error rate is impacted by the inclusion of non-Clifford gates. Finally, we demonstrate that our technique scales to many qubits with experiments on a 27-qubit IBM Q processor. We use our technique to quantify the impact of crosstalk on this 27-qubit device, and we find that it contributes approximately 2/3 of the total error per gate in random many-qubit circuit layers. △ Less

Submitted 10 October, 2023; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: 15 pages + appendices. v2: Corrected error in Section IIIA step 3 + minor edits. v3: Added new theory results

Journal ref: Phys. Rev. X 13, 041030 (2023)

An investigation of Hertzian contact in soft materials using photoelastic tomography

Authors: Benjamin Mitchell, Yuto Yokoyama, Ali Nassiri, Yoshiyuki Tagawa, Yannis Korkolis, Brad Kinsey

Abstract: Hertzian contact of a rigid sphere and a highly deformable soft solid is investigated using integrated photoelasticity. The experiments are performed by pressing a styrene sphere of 15 mm diameter against a 44 x 44 x 47 mm$^3$ cuboid made of 5% wt. gelatin, inside a circular polariscope, and with a range of forces. The emerging light rays are processed by considering that the retardation of each r… ▽ More Hertzian contact of a rigid sphere and a highly deformable soft solid is investigated using integrated photoelasticity. The experiments are performed by pressing a styrene sphere of 15 mm diameter against a 44 x 44 x 47 mm$^3$ cuboid made of 5% wt. gelatin, inside a circular polariscope, and with a range of forces. The emerging light rays are processed by considering that the retardation of each ray carries the cumulative effect of traversing the contact-induced axisymmetric stress field. Then, assuming Hertzian theory is valid, the retardation is analytically calculated for each ray and compared to the experimental one. Furthermore, a finite element model of the process introduces the effect of finite displacements and strains. Beyond the qualitative comparison of the retardation fields, the experimental, theoretical, and numerical results are quantitatively compared in terms of the maximum equivalent stress, surface displacement, and contact radius dimensions. A favorable agreement is found at lower force levels, where the assumptions of Hertz theory hold, whereas deviations are observed at higher force levels. A major discovery of this work is that at the maximum equivalent stress location, all three components of principal stress can be determined experimentally, and show satisfactory agreement with theoretical and numerical ones in our measurement range. This provides valuable insight into Hertzian contact problems since the maximum equivalent stress controls the initiation of plastic deformation or failure. The measured displacement and contact radii also reasonably agree with the theoretical and numerical ones. Finally, the limitations that arise due to the linearization of this problem are explored. △ Less

Submitted 8 July, 2022; originally announced July 2022.

High-Fidelity Qutrit Entangling Gates for Superconducting Circuits

Authors: Noah Goss, Alexis Morvan, Brian Marinelli, Bradley K. Mitchell, Long B. Nguyen, Ravi K. Naik, Larry Chen, Christian Jünger, John Mark Kreikebaum, David I. Santiago, Joel J. Wallman, Irfan Siddiqi

Abstract: Ternary quantum information processing in superconducting devices poses a promising alternative to its more popular binary counterpart through larger, more connected computational spaces and proposed advantages in quantum simulation and error correction. Although generally operated as qubits, transmons have readily addressable higher levels, making them natural candidates for operation as quantum… ▽ More Ternary quantum information processing in superconducting devices poses a promising alternative to its more popular binary counterpart through larger, more connected computational spaces and proposed advantages in quantum simulation and error correction. Although generally operated as qubits, transmons have readily addressable higher levels, making them natural candidates for operation as quantum three-level systems (qutrits). Recent works in transmon devices have realized high fidelity single qutrit operation. Nonetheless, effectively engineering a high-fidelity two-qutrit entanglement remains a central challenge for realizing qutrit processing in a transmon device. In this work, we apply the differential AC Stark shift to implement a flexible, microwave-activated, and dynamic cross-Kerr entanglement between two fixed-frequency transmon qutrits, expanding on work performed for the $ZZ$ interaction with transmon qubits. We then use this interaction to engineer efficient, high-fidelity qutrit CZ$^†$ and CZ gates, with estimated process fidelities of 97.3(1)% and 95.2(3)% respectively, a significant step forward for operating qutrits on a multi-transmon device. △ Less

Submitted 29 June, 2023; v1 submitted 14 June, 2022; originally announced June 2022.

Comments: 8 pages (+19 pages in supplement) 5 Figures (+7 in supplement)

On Flipping the Fréchet distance

Authors: Omrit Filtser, Mayank Goswami, Joseph S. B. Mitchell, Valentin Polishchuk

Abstract: The classical and extensively-studied Fréchet distance between two curves is defined as an inf max, where the infimum is over all traversals of the curves, and the maximum is over all concurrent positions of the two agents. In this article we investigate a "flipped" Fréchet measure defined by a sup min -- the supremum is over all traversals of the curves, and the minimum is over all concurrent pos… ▽ More The classical and extensively-studied Fréchet distance between two curves is defined as an inf max, where the infimum is over all traversals of the curves, and the maximum is over all concurrent positions of the two agents. In this article we investigate a "flipped" Fréchet measure defined by a sup min -- the supremum is over all traversals of the curves, and the minimum is over all concurrent positions of the two agents. This measure produces a notion of "social distance" between two curves (or general domains), where agents traverse curves while trying to stay as far apart as possible. We first study the flipped Fréchet measure between two polygonal curves in one and two dimensions, providing conditional lower bounds and matching algorithms. We then consider this measure on polygons, where it denotes the minimum distance that two agents can maintain while restricted to travel in or on the boundary of the same polygon. We investigate several variants of the problem in this setting, for some of which we provide linear time algorithms. Finally, we consider this measure on graphs. We draw connections between our proposed flipped Fréchet measure and existing related work in computational geometry, hoping that our new measure may spawn investigations akin to those performed for the Fréchet distance, and into further interesting problems that arise. △ Less

Submitted 21 September, 2022; v1 submitted 9 March, 2022; originally announced March 2022.

Scalable High-Performance Fluxonium Quantum Processor

Authors: Long B. Nguyen, Gerwin Koolstra, Yosep Kim, Alexis Morvan, Trevor Chistolini, Shraddha Singh, Konstantin N. Nesterov, Christian Jünger, Larry Chen, Zahra Pedramrazi, Bradley K. Mitchell, John Mark Kreikebaum, Shruti Puri, David I. Santiago, Irfan Siddiqi

Abstract: The technological development of hardware heading toward universal fault-tolerant quantum computation requires a large-scale processing unit with high performance. While fluxonium qubits are promising with high coherence and large anharmonicity, their scalability has not been systematically explored. In this work, we propose a superconducting quantum information processor based on compact high-coh… ▽ More The technological development of hardware heading toward universal fault-tolerant quantum computation requires a large-scale processing unit with high performance. While fluxonium qubits are promising with high coherence and large anharmonicity, their scalability has not been systematically explored. In this work, we propose a superconducting quantum information processor based on compact high-coherence fluxoniums with suppressed crosstalk, reduced design complexity, improved operational efficiency, high-fidelity gates, and resistance to parameter fluctuations. In this architecture, the qubits are readout dispersively using individual resonators connected to a common bus and manipulated via combined on-chip RF and DC control lines, both of which can be designed to have low crosstalk. A multi-path coupling approach enables exchange interactions between the high-coherence computational states and at the same time suppresses the spurious static ZZ rate, leading to fast and high-fidelity entangling gates. We numerically investigate the cross resonance controlled-NOT and the differential AC-Stark controlled-Z operations, revealing low gate error for qubit-qubit detuning bandwidth of up to 1 GHz. Our study on frequency crowding indicates high fabrication yield for quantum processors consisting of over thousands of qubits. In addition, we estimate low resource overhead to suppress logical error rate using the XZZX surface code. These results promise a scalable quantum architecture with high performance for the pursuit of universal quantum computation. △ Less

Submitted 5 February, 2022; v1 submitted 23 January, 2022; originally announced January 2022.

Comments: 29 pages, 15 figures

Journal ref: PRX Quantum 3, 037001 (2022)

Area-Optimal Simple Polygonalizations: The CG Challenge 2019

Authors: Erik D. Demaine, Sándor P. Fekete, Phillip Keldenich, Dominik Krupke, Joseph S. B. Mitchell

Abstract: We give an overview of theoretical and practical aspects of finding a simple polygon of minimum (Min-Area) or maximum (Max-Area) possible area for a given set of n points in the plane. Both problems are known to be NP-hard and were the subject of the 2019 Computational Geometry Challenge, which presented the quest of finding good solutions to more than 200 instances, ranging from n = 10 all the wa… ▽ More We give an overview of theoretical and practical aspects of finding a simple polygon of minimum (Min-Area) or maximum (Max-Area) possible area for a given set of n points in the plane. Both problems are known to be NP-hard and were the subject of the 2019 Computational Geometry Challenge, which presented the quest of finding good solutions to more than 200 instances, ranging from n = 10 all the way to n = 1, 000, 000. △ Less

Submitted 14 November, 2021; originally announced November 2021.

Comments: 12 pages, 9 Futures, 1 table

ACM Class: F.2.2

Homography-based Visual Servoing with Remote Center of Motion for Semi-autonomous Robotic Endoscope Manipulation

Authors: Martin Huber, John Bason Mitchell, Ross Henry, Sébastien Ourselin, Tom Vercauteren, Christos Bergeles

Abstract: The dominant visual servoing approaches in Minimally Invasive Surgery (MIS) follow single points or adapt the endoscope's field of view based on the surgical tools' distance. These methods rely on point positions with respect to the camera frame to infer a control policy. Deviating from the dominant methods, we formulate a robotic controller that allows for image-based visual servoing that require… ▽ More The dominant visual servoing approaches in Minimally Invasive Surgery (MIS) follow single points or adapt the endoscope's field of view based on the surgical tools' distance. These methods rely on point positions with respect to the camera frame to infer a control policy. Deviating from the dominant methods, we formulate a robotic controller that allows for image-based visual servoing that requires neither explicit tool and camera positions nor any explicit image depth information. The proposed method relies on homography-based image registration, which changes the automation paradigm from point-centric towards surgical-scene-centric approach. It simultaneously respects a programmable Remote Center of Motion (RCM). Our approach allows a surgeon to build a graph of desired views, from which, once built, views can be manually selected and automatically servoed to irrespective of robot-patient frame transformation changes. We evaluate our method on an abdominal phantom and provide an open source ROS Moveit integration for use with any serial manipulator. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: Accepted as conference paper at ISMR 2021

Optimal efficiency of high frequency chest wall oscillations and links with resistance and compliance in a model of the lung

Authors: Michaël Brunengo, Barrett R. Mitchell, Antonello Nicolini, Bernard Rousselet, Benjamin Mauroy

Abstract: Chest physiotherapy is a set of techniques used to help the draining of the mucus from the lung in pathological situations. The choice of the techniques, and their adjustment to the patients or to the pathologies, remains as of today largely empirical. High Frequency Chest Wall Oscillation (HFCWO) is one of these techniques, performed with a device that applies oscillating pressures on the chest.… ▽ More Chest physiotherapy is a set of techniques used to help the draining of the mucus from the lung in pathological situations. The choice of the techniques, and their adjustment to the patients or to the pathologies, remains as of today largely empirical. High Frequency Chest Wall Oscillation (HFCWO) is one of these techniques, performed with a device that applies oscillating pressures on the chest. However, there is no clear understanding of how HFCWO devices interact with the lung biomechanics. Hence, we study idealised HFCWO manipulations applied to a mathematical and numerical model of the biomechanics of the lung. The lung is represented by a fluid--structure interaction model based on an airway tree that is coupled to an homogeneous elastic medium. We show that our model is driven by two dimensionless numbers that drive the effect of the idealised HFCWO manipulation on the model of the lung. Our model allow to analyze the stress applied to an idealised mucus by the air--mucus interaction and by the airway walls deformation. This stress behaves as a buffer and has the effect of reducing the stress needed to overcome the idealised mucus yield stress. Moreover, our model predicts the existence of an optimal range of the working frequencies of HFCWO. This range is in agreement with the frequencies actually used by practitioners during HFCWO maneuvers. Finally, our model suggests that analyzing the mouth airflow during HFCWO maneuvers could allow to estimate the compliance and the hydrodynamic resistance of the lung of a patient. △ Less

Submitted 26 November, 2021; v1 submitted 5 August, 2021; originally announced August 2021.

Shortcut Hulls: Vertex-restricted Outer Simplifications of Polygons

Authors: Annika Bonerath, Jan-Henrik Haunert, Joseph S. B. Mitchell, Benjamin Niedermann

Abstract: Let $P$ be a crossing-free polygon and $\mathcal C$ a set of shortcuts, where each shortcut is a directed straight-line segment connecting two vertices of $P$. A shortcut hull of $P$ is another crossing-free polygon that encloses $P$ and whose oriented boundary is composed of elements from $\mathcal C$. Shortcut hulls find their application in geo-related problems such as the simplification of con… ▽ More Let $P$ be a crossing-free polygon and $\mathcal C$ a set of shortcuts, where each shortcut is a directed straight-line segment connecting two vertices of $P$. A shortcut hull of $P$ is another crossing-free polygon that encloses $P$ and whose oriented boundary is composed of elements from $\mathcal C$. Shortcut hulls find their application in geo-related problems such as the simplification of contour lines. We aim at a shortcut hull that linearly balances the enclosed area and perimeter. If no holes in the shortcut hull are allowed, the problem admits a straight-forward solution via shortest paths. For the more challenging case that the shortcut hull may contain holes, we present a polynomial-time algorithm that is based on computing a constrained, weighted triangulation of the input polygon's exterior. We use this problem as a starting point for investigating further variants, e.g., restricting the number of edges or bends. We demonstrate that shortcut hulls can be used for drawing the rough extent of point sets as well as for the schematization of polygons. △ Less

Submitted 25 June, 2021; originally announced June 2021.

Comments: Appears in the Proceedings of 33rd Canadian Conference on Computational Geometry, 2021

Hardware-Efficient Microwave-Activated Tunable Coupling Between Superconducting Qubits

Authors: Bradley K. Mitchell, Ravi K. Naik, Alexis Morvan, Akel Hashim, John Mark Kreikebaum, Brian Marinelli, Wim Lavrijsen, Kasra Nowrouzi, David I. Santiago, Irfan Siddiqi

Abstract: Generating high-fidelity, tunable entanglement between qubits is crucial for realizing gate-based quantum computation. In superconducting circuits, tunable interactions are often implemented using flux-tunable qubits or coupling elements, adding control complexity and noise sources. Here, we realize a tunable $ZZ$ interaction between two transmon qubits with fixed frequencies and fixed coupling, i… ▽ More Generating high-fidelity, tunable entanglement between qubits is crucial for realizing gate-based quantum computation. In superconducting circuits, tunable interactions are often implemented using flux-tunable qubits or coupling elements, adding control complexity and noise sources. Here, we realize a tunable $ZZ$ interaction between two transmon qubits with fixed frequencies and fixed coupling, induced by driving both transmons off-resonantly. We show tunable coupling over one order of magnitude larger than the static coupling, and change the sign of the interaction, enabling cancellation of the idle coupling. Further, this interaction is amenable to large quantum processors: the drive frequency can be flexibly chosen to avoid spurious transitions, and because both transmons are driven, it is resilient to microwave crosstalk. We apply this interaction to implement a controlled phase (CZ) gate with a gate fidelity of $99.43(1)\%$ as measured by cycle benchmarking, and we find the fidelity is limited by incoherent errors. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: 12 pages, 7 figures

Automatic Qubit Characterization and Gate Optimization with QubiC

Authors: Yilun Xu, Gang Huang, Jan Balewski, Ravi K. Naik, Alexis Morvan, Brad Mitchell, Kasra Nowrouzi, David I. Santiago, Irfan Siddiqi

Abstract: As the size and complexity of a quantum computer increases, quantum bit (qubit) characterization and gate optimization become complex and time-consuming tasks. Current calibration techniques require complicated and verbose measurements to tune up qubits and gates, which cannot easily expand to the large-scale quantum systems. We develop a concise and automatic calibration protocol to characterize… ▽ More As the size and complexity of a quantum computer increases, quantum bit (qubit) characterization and gate optimization become complex and time-consuming tasks. Current calibration techniques require complicated and verbose measurements to tune up qubits and gates, which cannot easily expand to the large-scale quantum systems. We develop a concise and automatic calibration protocol to characterize qubits and optimize gates using QubiC, which is an open source FPGA (field-programmable gate array) based control and measurement system for superconducting quantum information processors. We propose mutli-dimensional loss-based optimization of single-qubit gates and full XY-plane measurement method for the two-qubit CNOT gate calibration. We demonstrate the QubiC automatic calibration protocols are capable of delivering high-fidelity gates on the state-of-the-art transmon-type processor operating at the Advanced Quantum Testbed at Lawrence Berkeley National Laboratory. The single-qubit and two-qubit Clifford gate infidelities measured by randomized benchmarking are of $4.9(1.1) \times 10^{-4}$ and $1.4(3) \times 10^{-2}$, respectively. △ Less

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

Comments: Added funding agency

Leveraging Randomized Compiling for the QITE Algorithm

Authors: Jean-Loup Ville, Alexis Morvan, Akel Hashim, Ravi K. Naik, Marie Lu, Bradley Mitchell, John-Mark Kreikebaum, Kevin P. O'Brien, Joel J. Wallman, Ian Hincks, Joseph Emerson, Ethan Smith, Ed Younis, Costin Iancu, David I. Santiago, Irfan Siddiqi

Abstract: The success of the current generation of Noisy Intermediate-Scale Quantum (NISQ) hardware shows that quantum hardware may be able to tackle complex problems even without error correction. One outstanding issue is that of coherent errors arising from the increased complexity of these devices. These errors can accumulate through a circuit, making their impact on algorithms hard to predict and mitiga… ▽ More The success of the current generation of Noisy Intermediate-Scale Quantum (NISQ) hardware shows that quantum hardware may be able to tackle complex problems even without error correction. One outstanding issue is that of coherent errors arising from the increased complexity of these devices. These errors can accumulate through a circuit, making their impact on algorithms hard to predict and mitigate. Iterative algorithms like Quantum Imaginary Time Evolution are susceptible to these errors. This article presents the combination of both noise tailoring using Randomized Compiling and error mitigation with a purification. We also show that Cycle Benchmarking gives an estimate of the reliability of the purification. We apply this method to the Quantum Imaginary Time Evolution of a Transverse Field Ising Model and report an energy estimation and a ground state infidelity both below 1\%. Our methodology is general and can be used for other algorithms and platforms. We show how combining noise tailoring and error mitigation will push forward the performance of NISQ devices. △ Less

Submitted 26 October, 2021; v1 submitted 18 April, 2021; originally announced April 2021.

Comments: Data has been updated

Computing Coordinated Motion Plans for Robot Swarms: The CG:SHOP Challenge 2021

Authors: Sándor P. Fekete, Phillip Keldenich, Dominik Krupke, Joseph S. B. Mitchell

Abstract: We give an overview of the 2021 Computational Geometry Challenge, which targeted the problem of optimally coordinating a set of robots by computing a family of collision-free trajectories for a set set S of n pixel-shaped objects from a given start configuration into a desired target configuration. We give an overview of the 2021 Computational Geometry Challenge, which targeted the problem of optimally coordinating a set of robots by computing a family of collision-free trajectories for a set set S of n pixel-shaped objects from a given start configuration into a desired target configuration. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Comments: 13 pages, 8 figures, 2 tables

ACM Class: F.2.2; I.2.9

Approximating Maximum Independent Set for Rectangles in the Plane

Authors: Joseph S. B. Mitchell

Abstract: We give a polynomial-time constant-factor approximation algorithm for maximum independent set for (axis-aligned) rectangles in the plane. Using a polynomial-time algorithm, the best approximation factor previously known is $O(\log\log n)$. The results are based on a new form of recursive partitioning in the plane, in which faces that are constant-complexity and orthogonally convex are recursively… ▽ More We give a polynomial-time constant-factor approximation algorithm for maximum independent set for (axis-aligned) rectangles in the plane. Using a polynomial-time algorithm, the best approximation factor previously known is $O(\log\log n)$. The results are based on a new form of recursive partitioning in the plane, in which faces that are constant-complexity and orthogonally convex are recursively partitioned into a constant number of such faces. △ Less

Submitted 5 July, 2021; v1 submitted 1 January, 2021; originally announced January 2021.

QubiC: An open source FPGA-based control and measurement system for superconducting quantum information processors

Authors: Yilun Xu, Gang Huang, Jan Balewski, Ravi Naik, Alexis Morvan, Bradley Mitchell, Kasra Nowrouzi, David I. Santiago, Irfan Siddiqi

Abstract: As quantum information processors grow in quantum bit (qubit) count and functionality, the control and measurement system becomes a limiting factor to large scale extensibility. To tackle this challenge and keep pace with rapidly evolving classical control requirements, full control stack access is essential to system level optimization. We design a modular FPGA (field-programmable gate array) bas… ▽ More As quantum information processors grow in quantum bit (qubit) count and functionality, the control and measurement system becomes a limiting factor to large scale extensibility. To tackle this challenge and keep pace with rapidly evolving classical control requirements, full control stack access is essential to system level optimization. We design a modular FPGA (field-programmable gate array) based system called QubiC to control and measure a superconducting quantum processing unit. The system includes room temperature electronics hardware, FPGA gateware, and engineering software. A prototype hardware module is assembled from several commercial off-the-shelf evaluation boards and in-house developed circuit boards. Gateware and software are designed to implement basic qubit control and measurement protocols. System functionality and performance are demonstrated by performing qubit chip characterization, gate optimization, and randomized benchmarking sequences on a superconducting quantum processor operating at the Advanced Quantum Testbed at Lawrence Berkeley National Laboratory. The single-qubit and two-qubit process fidelities are measured to be 0.9980$\pm$0.0001 and 0.948$\pm$0.004 by randomized benchmarking. With fast circuit sequence loading capability, the QubiC performs randomized compiling experiments efficiently and improves the feasibility of executing more complex algorithms. △ Less

Submitted 27 October, 2021; v1 submitted 31 December, 2020; originally announced January 2021.

Comments: Final published version on IEEE Transactions on Quantum Engineering

Journal ref: IEEE Transactions on Quantum Engineering, vol. 2, pp. 1-11, 2021, Art no. 6003811

Randomized compiling for scalable quantum computing on a noisy superconducting quantum processor

Authors: Akel Hashim, Ravi K. Naik, Alexis Morvan, Jean-Loup Ville, Bradley Mitchell, John Mark Kreikebaum, Marc Davis, Ethan Smith, Costin Iancu, Kevin P. O'Brien, Ian Hincks, Joel J. Wallman, Joseph Emerson, Irfan Siddiqi

Abstract: The successful implementation of algorithms on quantum processors relies on the accurate control of quantum bits (qubits) to perform logic gate operations. In this era of noisy intermediate-scale quantum (NISQ) computing, systematic miscalibrations, drift, and crosstalk in the control of qubits can lead to a coherent form of error which has no classical analog. Coherent errors severely limit the p… ▽ More The successful implementation of algorithms on quantum processors relies on the accurate control of quantum bits (qubits) to perform logic gate operations. In this era of noisy intermediate-scale quantum (NISQ) computing, systematic miscalibrations, drift, and crosstalk in the control of qubits can lead to a coherent form of error which has no classical analog. Coherent errors severely limit the performance of quantum algorithms in an unpredictable manner, and mitigating their impact is necessary for realizing reliable quantum computations. Moreover, the average error rates measured by randomized benchmarking and related protocols are not sensitive to the full impact of coherent errors, and therefore do not reliably predict the global performance of quantum algorithms, leaving us unprepared to validate the accuracy of future large-scale quantum computations. Randomized compiling is a protocol designed to overcome these performance limitations by converting coherent errors into stochastic noise, dramatically reducing unpredictable errors in quantum algorithms and enabling accurate predictions of algorithmic performance from error rates measured via cycle benchmarking. In this work, we demonstrate significant performance gains under randomized compiling for the four-qubit quantum Fourier transform algorithm and for random circuits of variable depth on a superconducting quantum processor. Additionally, we accurately predict algorithm performance using experimentally-measured error rates. Our results demonstrate that randomized compiling can be utilized to leverage and predict the capabilities of modern-day noisy quantum processors, paving the way forward for scalable quantum computing. △ Less

Submitted 11 May, 2021; v1 submitted 1 October, 2020; originally announced October 2020.

Journal ref: Phys. Rev. X 11, 041039 (2021)

Qutrit randomized benchmarking

Authors: A. Morvan, V. V. Ramasesh, M. S. Blok, J. M. Kreikebaum, K. O'Brien, L. Chen, B. K. Mitchell, R. K. Naik, D. I. Santiago, I. Siddiqi

Abstract: Ternary quantum processors offer significant computational advantages over conventional qubit technologies, leveraging the encoding and processing of quantum information in qutrits (three-level systems). To evaluate and compare the performance of such emerging quantum hardware it is essential to have robust benchmarking methods suitable for a higher-dimensional Hilbert space. We demonstrate extens… ▽ More Ternary quantum processors offer significant computational advantages over conventional qubit technologies, leveraging the encoding and processing of quantum information in qutrits (three-level systems). To evaluate and compare the performance of such emerging quantum hardware it is essential to have robust benchmarking methods suitable for a higher-dimensional Hilbert space. We demonstrate extensions of industry standard Randomized Benchmarking (RB) protocols, developed and used extensively for qubits, suitable for ternary quantum logic. Using a superconducting five-qutrit processor, we find a single-qutrit gate infidelity as low as $2.38 \times 10^{-3}$. Through interleaved RB, we find that this qutrit gate error is largely limited by the native (qubit-like) gate fidelity, and employ simultaneous RB to fully characterize cross-talk errors. Finally, we apply cycle benchmarking to a two-qutrit CSUM gate and obtain a two-qutrit process fidelity of $0.82$. Our results demonstrate a RB-based tool to characterize the obtain overall performance of a qutrit processor, and a general approach to diagnose control errors in future qudit hardware. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: 6 pages (+ 2 pages supplement), 5 figures

Journal ref: Phys. Rev. Lett. 126, 210504 (2021)

Two-pathogen model with competition on clustered networks

Authors: Peter Mann, V. Anne Smith, John B. O. Mitchell, Simon Dobson

Abstract: Networks provide a mathematically rich framework to represent social contacts sufficient for the transmission of disease. Social networks are often highly clustered and fail to be locally tree-like. In this paper, we study the effects of clustering on the spread of sequential strains of a pathogen using the generating function formulation under a complete cross-immunity coupling, deriving conditio… ▽ More Networks provide a mathematically rich framework to represent social contacts sufficient for the transmission of disease. Social networks are often highly clustered and fail to be locally tree-like. In this paper, we study the effects of clustering on the spread of sequential strains of a pathogen using the generating function formulation under a complete cross-immunity coupling, deriving conditions for the threshold of coexistence of the second strain. We show that clustering reduces the coexistence threshold of the second strain and its outbreak size in Poisson networks, whilst exhibiting the opposite effects on uniform-degree models. We conclude that clustering within a population must increase the ability of the second wave of an epidemic to spread over a network. We apply our model to the study of multilayer clustered networks and observe the fracturing of the residual graph at two distinct transmissibilities. △ Less

Submitted 20 November, 2020; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: 9 pages, 5 figures

Journal ref: Phys. Rev. E 103, 062308 (2021)

Cutting Polygons into Small Pieces with Chords: Laser-Based Localization

Authors: Esther M. Arkin, Rathish Das, Jie Gao, Mayank Goswami, Joseph S. B. Mitchell, Valentin Polishchuk, Csaba D. Toth

Abstract: Motivated by indoor localization by tripwire lasers, we study the problem of cutting a polygon into small-size pieces, using the chords of the polygon. Several versions are considered, depending on the definition of the "size" of a piece. In particular, we consider the area, the diameter, and the radius of the largest inscribed circle as a measure of the size of a piece. We also consider different… ▽ More Motivated by indoor localization by tripwire lasers, we study the problem of cutting a polygon into small-size pieces, using the chords of the polygon. Several versions are considered, depending on the definition of the "size" of a piece. In particular, we consider the area, the diameter, and the radius of the largest inscribed circle as a measure of the size of a piece. We also consider different objectives, either minimizing the maximum size of a piece for a given number of chords, or minimizing the number of chords that achieve a given size threshold for the pieces. We give hardness results for polygons with holes and approximation algorithms for multiple variants of the problem. △ Less

Submitted 26 June, 2020; originally announced June 2020.

Comments: This paper will appear in ESA2020, Track A proceedings

Random graphs with arbitrary clustering and their applications

Authors: Peter Mann, V. Anne Smith, John B. O. Mitchell, Simon Dobson

Abstract: The structure of many real networks is not locally tree-like and hence, network analysis fails to characterise their bond percolation properties. In a recent paper [P. Mann, V. A. Smith, J. B. O. Mitchell, and S. Dobson, Percolation in random graphs with higher-order clustering, arXiv e-prints, p. arXiv:2006.06744, June 2020.], we developed analytical solutions to the percolation properties of ran… ▽ More The structure of many real networks is not locally tree-like and hence, network analysis fails to characterise their bond percolation properties. In a recent paper [P. Mann, V. A. Smith, J. B. O. Mitchell, and S. Dobson, Percolation in random graphs with higher-order clustering, arXiv e-prints, p. arXiv:2006.06744, June 2020.], we developed analytical solutions to the percolation properties of random networks with homogeneous clustering (clusters whose nodes are degree-equivalent). In this paper, we extend this model to investigate networks that contain clusters whose nodes are not degree-equivalent, including multilayer networks. Through numerical examples we show how this method can be used to investigate the properties of random complex networks with arbitrary clustering, extending the applicability of the configuration model and generating function formulation. △ Less

Submitted 20 November, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 11 pages, 10 figures

Journal ref: Phys. Rev. E 103, 012309 (2021)

Percolation in random graphs with higher-order clustering

Authors: Peter Mann, V. Anne Smith, John B. O. Mitchell, Simon Dobson

Abstract: Percolation theory can be used to describe the structural properties of complex networks using the generating function formulation. This mapping assumes that the network is locally tree-like and does not contain short-range loops between neighbours. In this paper we use the generating function formulation to examine clustered networks that contain simple cycles and cliques of any order. We use the… ▽ More Percolation theory can be used to describe the structural properties of complex networks using the generating function formulation. This mapping assumes that the network is locally tree-like and does not contain short-range loops between neighbours. In this paper we use the generating function formulation to examine clustered networks that contain simple cycles and cliques of any order. We use the natural generalisation to the Molloy-Reed criterion for these networks to describe their critical properties and derive an analytical description of the size of the giant component, providing solutions for Poisson and power-law networks. We find that networks comprising larger simple cycles behave increasingly more tree-like. Conversley, clustering comprised of larger cliques increasingly deviate from the tree-like solution, although the behaviour is strongly dependent on the degree-assortativity. △ Less

Submitted 20 November, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 11 pages, 6 figures

Journal ref: Phys. Rev. E 103, 012313 (2021)

Planar Bichromatic Bottleneck Spanning Trees

Authors: A. Karim Abu-Affash, Sujoy Bhore, Paz Carmi, Joseph S. B. Mitchell

Abstract: Given a set $P$ of $n$ red and blue points in the plane, a \emph{planar bichromatic spanning tree} of $P$ is a spanning tree of $P$, such that each edge connects between a red and a blue point, and no two edges intersect. In the bottleneck planar bichromatic spanning tree problem, the goal is to find a planar bichromatic spanning tree $T$, such that the length of the longest edge in $T$ is minimiz… ▽ More Given a set $P$ of $n$ red and blue points in the plane, a \emph{planar bichromatic spanning tree} of $P$ is a spanning tree of $P$, such that each edge connects between a red and a blue point, and no two edges intersect. In the bottleneck planar bichromatic spanning tree problem, the goal is to find a planar bichromatic spanning tree $T$, such that the length of the longest edge in $T$ is minimized. In this paper, we show that this problem is NP-hard for points in general position. Moreover, we present a polynomial-time $(8\sqrt{2})$-approximation algorithm, by showing that any bichromatic spanning tree of bottleneck $λらむだ$ can be converted to a planar bichromatic spanning tree of bottleneck at most $8\sqrt{2}λらむだ$. △ Less

Submitted 19 April, 2020; originally announced April 2020.

Computing Convex Partitions for Point Sets in the Plane: The CG:SHOP Challenge 2020

Authors: Erik D. Demaine, Sándor P. Fekete, Phillip Keldenich, Dominik Krupke, Joseph S. B. Mitchell

Abstract: We give an overview of the 2020 Computational Geometry Challenge, which targeted the problem of partitioning the convex hull of a given planar point set P into the smallest number of convex faces, such that no point of P is contained in the interior of a face. We give an overview of the 2020 Computational Geometry Challenge, which targeted the problem of partitioning the convex hull of a given planar point set P into the smallest number of convex faces, such that no point of P is contained in the interior of a face. △ Less

Submitted 8 April, 2020; originally announced April 2020.

Comments: 8 pages, 3 figures, 1 table

ACM Class: F.2.2

Probing a Set of Trajectories to Maximize Captured Information

Authors: Sándor P. Fekete, Alexander Hill, Dominik Krupke, Tyler Mayer, Joseph S. B. Mitchell, Ojas Parekh, Cynthia A. Phillips

Abstract: We study a trajectory analysis problem we call the Trajectory Capture Problem (TCP), in which, for a given input set ${\cal T}$ of trajectories in the plane, and an integer $k\geq 2$, we seek to compute a set of $k$ points (``portals'') to maximize the total weight of all subtrajectories of ${\cal T}$ between pairs of portals. This problem naturally arises in trajectory analysis and summarization.… ▽ More We study a trajectory analysis problem we call the Trajectory Capture Problem (TCP), in which, for a given input set ${\cal T}$ of trajectories in the plane, and an integer $k\geq 2$, we seek to compute a set of $k$ points (``portals'') to maximize the total weight of all subtrajectories of ${\cal T}$ between pairs of portals. This problem naturally arises in trajectory analysis and summarization. We show that the TCP is NP-hard (even in very special cases) and give some first approximation results. Our main focus is on attacking the TCP with practical algorithm-engineering approaches, including integer linear programming (to solve instances to provable optimality) and local search methods. We study the integrality gap arising from such approaches. We analyze our methods on different classes of data, including benchmark instances that we generate. Our goal is to understand the best performing heuristics, based on both solution time and solution quality. We demonstrate that we are able to compute provably optimal solutions for real-world instances. △ Less

Submitted 7 April, 2020; originally announced April 2020.

Excitation efficiency and limitations of the luminescence of Eu3+ ions in GaN

Authors: Dolf Timmerman, Brandon Mitchell, Shuhei Ichikawa, Jun Tatebayashi, Masaaki Ashida, Yasufumi Fujiwara

Abstract: The excitation efficiency and external luminescence quantum efficiency of trivalent Eu3+ ions doped into gallium nitride (GaN) was studied under optical and electrical excitation. For small pump fluences it was found that the excitation of Eu3+ ions is limited by an efficient carrier trap that competes in the energy transfer from the host material. For large pump fluences the limited number of hig… ▽ More The excitation efficiency and external luminescence quantum efficiency of trivalent Eu3+ ions doped into gallium nitride (GaN) was studied under optical and electrical excitation. For small pump fluences it was found that the excitation of Eu3+ ions is limited by an efficient carrier trap that competes in the energy transfer from the host material. For large pump fluences the limited number of high-efficiency Eu3+ sites, and the small excitation cross-section of the majority Eu3+ site, limit the quantum efficiency. At low temperatures under optimal excitation conditions, the external luminescence quantum efficiency reached a value of 46%. These results show the high potential for this material as an efficient light emitter, and demonstrates the importance of the excitation conditions on the light output efficiency. △ Less

Submitted 3 September, 2019; originally announced September 2019.

Journal ref: Phys. Rev. Applied 13, 014044 (2020)

Existence and hardness of conveyor belts

Authors: Molly Baird, Sara C. Billey, Erik D. Demaine, Martin L. Demaine, David Eppstein, Sándor Fekete, Graham Gordon, Sean Griffin, Joseph S. B. Mitchell, Joshua P. Swanson

Abstract: An open problem of Manuel Abellanas asks whether every set of disjoint closed unit disks in the plane can be connected by a conveyor belt, which means a tight simple closed curve that touches the boundary of each disk, possibly multiple times. We prove three main results. First, for unit disks whose centers are both $x$-monotone and $y$-monotone, or whose centers have $x$-coordinates that differ b… ▽ More An open problem of Manuel Abellanas asks whether every set of disjoint closed unit disks in the plane can be connected by a conveyor belt, which means a tight simple closed curve that touches the boundary of each disk, possibly multiple times. We prove three main results. First, for unit disks whose centers are both $x$-monotone and $y$-monotone, or whose centers have $x$-coordinates that differ by at least two units, a conveyor belt always exists and can be found efficiently. Second, it is NP-complete to determine whether disks of varying radii have a conveyor belt, and it remains NP-complete when we constrain the belt to touch disks exactly once. Third, any disjoint set of $n$ disks of arbitrary radii can be augmented by $O(n)$ "guide" disks so that the augmented system has a conveyor belt touching each disk exactly once, answering a conjecture of Demaine, Demaine, and Palop. △ Less

Submitted 20 August, 2019; originally announced August 2019.

MSC Class: 52C26 ACM Class: G.2.0

Journal ref: Electronic J. Combinatorics 27 (4), Paper 4.25, 2020

Parallel Performance of Algebraic Multigrid Domain Decomposition (AMG-DD)

Authors: Wayne B. Mitchell, Robert Strzodka, Robert D. Falgout

Abstract: Algebraic multigrid (AMG) is a widely used scalable solver and preconditioner for large-scale linear systems resulting from the discretization of a wide class of elliptic PDEs. While AMG has optimal computational complexity, the cost of communication has become a significant bottleneck that limits its scalability as processor counts continue to grow on modern machines. This paper examines the desi… ▽ More Algebraic multigrid (AMG) is a widely used scalable solver and preconditioner for large-scale linear systems resulting from the discretization of a wide class of elliptic PDEs. While AMG has optimal computational complexity, the cost of communication has become a significant bottleneck that limits its scalability as processor counts continue to grow on modern machines. This paper examines the design, implementation, and parallel performance of a novel algorithm, Algebraic Multigrid Domain Decomposition (AMG-DD), designed specifically to limit communication. The goal of AMG-DD is to provide a low-communication alternative to standard AMG V-cycles by trading some additional computational overhead for a significant reduction in communication cost. Numerical results show that AMG-DD achieves superior accuracy per communication cost compared to AMG, and speedup over AMG is demonstrated on a large GPU cluster. △ Less

Submitted 21 January, 2020; v1 submitted 25 June, 2019; originally announced June 2019.

Data Races and the Discrete Resource-time Tradeoff Problem with Resource Reuse over Paths

Authors: Rathish Das, Shih-Yu Tsai, Sharmila Duppala, Jayson Lynch, Esther M. Arkin, Rezaul Chowdhury, Joseph S. B. Mitchell, Steven Skiena

Abstract: A determinacy race occurs if two or more logically parallel instructions access the same memory location and at least one of them tries to modify its content. Races often lead to nondeterministic and incorrect program behavior. A data race is a special case of a determinacy race which can be eliminated by associating a mutual-exclusion lock or allowing atomic accesses to the memory location. Howev… ▽ More A determinacy race occurs if two or more logically parallel instructions access the same memory location and at least one of them tries to modify its content. Races often lead to nondeterministic and incorrect program behavior. A data race is a special case of a determinacy race which can be eliminated by associating a mutual-exclusion lock or allowing atomic accesses to the memory location. However, such solutions can reduce parallelism by serializing all accesses to that location. For associative and commutative updates, reducers allow parallel race-free updates at the expense of using some extra space. We ask the following question. Given a fixed budget of extra space to mitigate the cost of races in a parallel program, which memory locations should be assigned reducers and how should the space be distributed among the reducers in order to minimize the overall running time? We argue that the races can be captured by a directed acyclic graph (DAG), with nodes representing memory cells and arcs representing read-write dependencies between cells. We then formulate our optimization problem on DAGs. We concentrate on a variation of this problem where space reuse among reducers is allowed by routing extra space along a source to sink path of the DAG and using it in the construction of reducers along the path. We consider two reducers and the corresponding duration functions (i.e., reduction time as a function of space budget). We generalize our race-avoiding space-time tradeoff problem to a discrete resource-time tradeoff problem with general non-increasing duration functions and resource reuse over paths. For general DAGs, the offline problem is strongly NP-hard under all three duration functions, and we give approximation algorithms. We also prove hardness of approximation for the general resource-time tradeoff problem and give a pseudo-polynomial time algorithm for series-parallel DAGs. △ Less

Submitted 19 April, 2019; originally announced April 2019.

The Balanced Connected Subgraph Problem

Authors: Sujoy Bhore, Sourav Chakraborty, Satyabrata Jana, Joseph S. B. Mitchell, Supantha Pandit, Sasanka Roy

Abstract: The problem of computing induced subgraphs that satisfy some specified restrictions arises in various applications of graph algorithms and has been well studied. In this paper, we consider the following Balanced Connected Subgraph (shortly, BCS) problem. The input is a graph $G=(V,E)$, with each vertex in the set $V$ having an assigned color, "red" or "blue". We seek a maximum-cardinality subset… ▽ More The problem of computing induced subgraphs that satisfy some specified restrictions arises in various applications of graph algorithms and has been well studied. In this paper, we consider the following Balanced Connected Subgraph (shortly, BCS) problem. The input is a graph $G=(V,E)$, with each vertex in the set $V$ having an assigned color, "red" or "blue". We seek a maximum-cardinality subset $V'\subseteq V$ of vertices that is color-balanced (having exactly $|V'|/2$ red nodes and $|V'|/2$ blue nodes), such that the subgraph induced by the vertex set $V'$ in $G$ is connected. We show that the BCS problem is NP-hard, even for bipartite graphs $G$ (with red/blue color assignment not necessarily being a proper 2-coloring). Further, we consider this problem for various classes of the input graph $G$, including, e.g., planar graphs, chordal graphs, trees, split graphs, bipartite graphs with a proper red/blue $2$-coloring, and graphs with diameter $2$. For each of these classes either we prove NP-hardness or design a polynomial time algorithm. △ Less

Submitted 24 September, 2018; originally announced September 2018.

Comments: 15 pages, 3 figures

Don't Rock the Boat: Algorithms for Balanced Dynamic Loading and Unloading

Authors: Sándor P. Fekete, Sven von Höveling, Joseph S. B. Mitchell, Christian Rieck, Christian Scheffer, Arne Schmidt, James R. Zuber

Abstract: We consider dynamic loading and unloading problems for heavy geometric objects. The challenge is to maintain balanced configurations at all times: minimize the maximal motion of the overall center of gravity. While this problem has been studied from an algorithmic point of view, previous work only focuses on balancing the final center of gravity; we give a variety of results for computing balanced… ▽ More We consider dynamic loading and unloading problems for heavy geometric objects. The challenge is to maintain balanced configurations at all times: minimize the maximal motion of the overall center of gravity. While this problem has been studied from an algorithmic point of view, previous work only focuses on balancing the final center of gravity; we give a variety of results for computing balanced loading and unloading schemes that minimize the maximal motion of the center of gravity during the entire process. In particular, we consider the one-dimensional case and distinguish between loading and unloading. In the unloading variant, the positions of the intervals are given, and we search for an optimal unloading order of the intervals. We prove that the unloading variant is NP-complete and give a 2.7-approximation algorithm. In the loading variant, we have to compute both the positions of the intervals and their loading order. We give optimal approaches for several variants that model different loading scenarios that may arise, e.g., in the loading of a transport ship with containers. △ Less

Submitted 17 January, 2018; v1 submitted 18 December, 2017; originally announced December 2017.

Comments: 23 pages, 3 figures, full version of conference paper to appear in LATIN 2018

ACM Class: F.2.2

Hi, how can I help you?: Automating enterprise IT support help desks

Authors: Senthil Mani, Neelamadhav Gantayat, Rahul Aralikatte, Monika Gupta, Sampath Dechu, Anush Sankaran, Shreya Khare, Barry Mitchell, Hemamalini Subramanian, Hema Venkatarangan

Abstract: Question answering is one of the primary challenges of natural language understanding. In realizing such a system, providing complex long answers to questions is a challenging task as opposed to factoid answering as the former needs context disambiguation. The different methods explored in the literature can be broadly classified into three categories namely: 1) classification based, 2) knowledge… ▽ More Question answering is one of the primary challenges of natural language understanding. In realizing such a system, providing complex long answers to questions is a challenging task as opposed to factoid answering as the former needs context disambiguation. The different methods explored in the literature can be broadly classified into three categories namely: 1) classification based, 2) knowledge graph based and 3) retrieval based. Individually, none of them address the need of an enterprise wide assistance system for an IT support and maintenance domain. In this domain the variance of answers is large ranging from factoid to structured operating procedures; the knowledge is present across heterogeneous data sources like application specific documentation, ticket management systems and any single technique for a general purpose assistance is unable to scale for such a landscape. To address this, we have built a cognitive platform with capabilities adopted for this domain. Further, we have built a general purpose question answering system leveraging the platform that can be instantiated for multiple products, technologies in the support domain. The system uses a novel hybrid answering model that orchestrates across a deep learning classifier, a knowledge graph based context disambiguation module and a sophisticated bag-of-words search system. This orchestration performs context switching for a provided question and also does a smooth hand-off of the question to a human expert if none of the automated techniques can provide a confident answer. This system has been deployed across 675 internal enterprise IT support and maintenance projects. △ Less

Submitted 2 November, 2017; originally announced November 2017.

Comments: To appear in IAAI 2018

Bayesian Inference over the Stiefel Manifold via the Givens Representation

Authors: Arya A Pourzanjani, Richard M Jiang, Brian Mitchell, Paul J Atzberger, Linda R Petzold

Abstract: We introduce an approach based on the Givens representation for posterior inference in statistical models with orthogonal matrix parameters, such as factor models and probabilistic principal component analysis (PPCA). We show how the Givens representation can be used to develop practical methods for transforming densities over the Stiefel manifold into densities over subsets of Euclidean space. We… ▽ More We introduce an approach based on the Givens representation for posterior inference in statistical models with orthogonal matrix parameters, such as factor models and probabilistic principal component analysis (PPCA). We show how the Givens representation can be used to develop practical methods for transforming densities over the Stiefel manifold into densities over subsets of Euclidean space. We show how to deal with issues arising from the topology of the Stiefel manifold and how to inexpensively compute the change-of-measure terms. We introduce an auxiliary parameter approach that limits the impact of topological issues. We provide both analysis of our methods and numerical examples demonstrating the effectiveness of the approach. We also discuss how our Givens representation can be used to define general classes of distributions over the space of orthogonal matrices. We then give demonstrations on several examples showing how the Givens approach performs in practice in comparison with other methods. △ Less

Submitted 2 November, 2019; v1 submitted 25 October, 2017; originally announced October 2017.

Network Optimization on Partitioned Pairs of Points

Authors: Esther M. Arkin, Aritra Banik, Paz Carmi, Gui Citovsky, Su Jia, Matthet J. Katz, Tyler Mayer, Joseph S. B. Mitchell

Abstract: Given $n$ pairs of points, $\mathcal{S} = \{\{p_1, q_1\}, \{p_2, q_2\}, \dots, \{p_n, q_n\}\}$, in some metric space, we study the problem of two-coloring the points within each pair, red and blue, to optimize the cost of a pair of node-disjoint networks, one over the red points and one over the blue points. In this paper we consider our network structures to be spanning trees, traveling salesman… ▽ More Given $n$ pairs of points, $\mathcal{S} = \{\{p_1, q_1\}, \{p_2, q_2\}, \dots, \{p_n, q_n\}\}$, in some metric space, we study the problem of two-coloring the points within each pair, red and blue, to optimize the cost of a pair of node-disjoint networks, one over the red points and one over the blue points. In this paper we consider our network structures to be spanning trees, traveling salesman tours or matchings. We consider several different weight functions computed over the network structures induced, as well as several different objective functions. We show that some of these problems are NP-hard, and provide constant factor approximation algorithms in all cases. △ Less

Submitted 2 October, 2017; originally announced October 2017.

TSP With Locational Uncertainty: The Adversarial Model

Authors: Gui Citovsky, Tyler Mayer, Joseph S. B. Mitchell

Abstract: In this paper we study a natural special case of the Traveling Salesman Problem (TSP) with point-locational-uncertainty which we will call the {\em adversarial TSP} problem (ATSP). Given a metric space $(X, d)$ and a set of subsets $R = \{R_1, R_2, ... , R_n\} : R_i \subseteq X$, the goal is to devise an ordering of the regions, $σしぐま_R$, that the tour will visit such that when a single point is chos… ▽ More In this paper we study a natural special case of the Traveling Salesman Problem (TSP) with point-locational-uncertainty which we will call the {\em adversarial TSP} problem (ATSP). Given a metric space $(X, d)$ and a set of subsets $R = \{R_1, R_2, ... , R_n\} : R_i \subseteq X$, the goal is to devise an ordering of the regions, $σしぐま_R$, that the tour will visit such that when a single point is chosen from each region, the induced tour over those points in the ordering prescribed by $σしぐま_R$ is as short as possible. Unlike the classical locational-uncertainty-TSP problem, which focuses on minimizing the expected length of such a tour when the point within each region is chosen according to some probability distribution, here, we focus on the {\em adversarial model} in which once the choice of $σしぐま_R$ is announced, an adversary selects a point from each region in order to make the resulting tour as long as possible. In other words, we consider an offline problem in which the goal is to determine an ordering of the regions $R$ that is optimal with respect to the "worst" point possible within each region being chosen by an adversary, who knows the chosen ordering. We give a $3$-approximation when $R$ is a set of arbitrary regions/sets of points in a metric space. We show how geometry leads to improved constant factor approximations when regions are parallel line segments of the same lengths, and a polynomial-time approximation scheme (PTAS) for the important special case in which $R$ is a set of disjoint unit disks in the plane. △ Less

Submitted 17 May, 2017; originally announced May 2017.

Comments: To appear, International Symposium on Computational Geometry (SoCG 2017)

ACM Class: F.2.2; G.2.2

Symmetric Assembly Puzzles are Hard, Beyond a Few Pieces

Authors: Erik D. Demaine, Matias Korman, Jason S. Ku, Joseph S. B. Mitchell, Yota Otachi, André van Renssen, Marcel Roeloffzen, Ryuhei Uehara, Yushi Uno

Abstract: We study the complexity of symmetric assembly puzzles: given a collection of simple polygons, can we translate, rotate, and possibly flip them so that their interior-disjoint union is line symmetric? On the negative side, we show that the problem is strongly NP-complete even if the pieces are all polyominos. On the positive side, we show that the problem can be solved in polynomial time if the num… ▽ More We study the complexity of symmetric assembly puzzles: given a collection of simple polygons, can we translate, rotate, and possibly flip them so that their interior-disjoint union is line symmetric? On the negative side, we show that the problem is strongly NP-complete even if the pieces are all polyominos. On the positive side, we show that the problem can be solved in polynomial time if the number of pieces is a fixed constant. △ Less

Submitted 7 April, 2019; v1 submitted 7 March, 2017; originally announced March 2017.

A PTAS for TSP with Neighborhoods Among Fat Regions in the Plane

Authors: Joseph S. B. Mitchell

Abstract: The Euclidean TSP with neighborhoods (TSPN) problem seeks a shortest tour that visits a given collection of $n$ regions ({\em neighborhoods}). We present the first polynomial-time approximation scheme for TSPN for a set of regions given by arbitrary disjoint fat regions in the plane. This improves substantially upon the known approximation algorithms, and is the first PTAS for TSPN on regions of n… ▽ More The Euclidean TSP with neighborhoods (TSPN) problem seeks a shortest tour that visits a given collection of $n$ regions ({\em neighborhoods}). We present the first polynomial-time approximation scheme for TSPN for a set of regions given by arbitrary disjoint fat regions in the plane. This improves substantially upon the known approximation algorithms, and is the first PTAS for TSPN on regions of non-comparable sizes. Our result is based on a novel extension of the $m$-guillotine method. The result applies to regions that are "fat" in a very weak sense: each region $P_i$ has area $Ωおめが([diam(P_i)]^2)$, but is otherwise arbitrary. △ Less

Submitted 5 March, 2017; originally announced March 2017.

Comments: Manuscript that was published in SODA 2007, with updates made after the proceedings version/talk

Approximation algorithms for TSP with neighborhoods in the plane

Authors: Adrian Dumitrescu, Joseph S. B. Mitchell

Abstract: In the Euclidean TSP with neighborhoods (TSPN), we are given a collection of n regions (neighborhoods) and we seek a shortest tour that visits each region. As a generalization of the classical Euclidean TSP, TSPN is also NP-hard. In this paper, we present new approximation results for the TSPN, including (1) a constant-factor approximation algorithm for the case of arbitrary connected neighborhood… ▽ More In the Euclidean TSP with neighborhoods (TSPN), we are given a collection of n regions (neighborhoods) and we seek a shortest tour that visits each region. As a generalization of the classical Euclidean TSP, TSPN is also NP-hard. In this paper, we present new approximation results for the TSPN, including (1) a constant-factor approximation algorithm for the case of arbitrary connected neighborhoods having comparable diameters; and (2) a PTAS for the important special case of disjoint unit disk neighborhoods (or nearly disjoint, nearly-unit disks). Our methods also yield improved approximation ratios for various special classes of neighborhoods, which have previously been studied. Further, we give a linear-time O(1)-approximation algorithm for the case of neighborhoods that are (infinite) straight lines. △ Less

Submitted 5 March, 2017; originally announced March 2017.

Comments: Manuscript that was published in J. Algorithms (2003), with brief clarifying remarks added (2014)

Journal ref: Journal of Algorithms, Volume 48, Issue 1, August 2003, Pages 135-159. (special issue for SODA 2001)

Universal Guard Problems

Authors: Sándor P. Fekete, Qian Li, Joseph S. B. Mitchell, Christian Scheffer

Abstract: We provide a spectrum of results for the Universal Guard Problem, in which one is to obtain a small set of points ("guards") that are "universal" in their ability to guard any of a set of possible polygonal domains in the plane. We give upper and lower bounds on the number of universal guards that are always sufficient to guard all polygons having a given set of n vertices, or to guard all polygon… ▽ More We provide a spectrum of results for the Universal Guard Problem, in which one is to obtain a small set of points ("guards") that are "universal" in their ability to guard any of a set of possible polygonal domains in the plane. We give upper and lower bounds on the number of universal guards that are always sufficient to guard all polygons having a given set of n vertices, or to guard all polygons in a given set of k polygons on an n-point vertex set. Our upper bound proofs include algorithms to construct universal guard sets of the respective cardinalities. △ Less

Submitted 27 March, 2017; v1 submitted 24 November, 2016; originally announced November 2016.

Comments: 28 pages, 19 figures, full version of extended abstract that appeared in the 27th International Symposium on Algorithms and Computation (ISAAC 2016), 32:1-32:13

ACM Class: F.2.2

DUSTiNGS III: Distribution of Intermediate-Age and Old Stellar Populations in Disks and Outer Extremities of Dwarf Galaxies

Authors: Kristen B. W. McQuinn, Martha L. Boyer, Mallory B. Mitchell, Evan D. Skillman, R. D. Gehrz, Martin A. T. Groenewegen, Iain McDonald, G. C. Sloan, Jacco Th. van Loon, Patricia A. Whitelock, Albert A. Zijlstra

Abstract: We have traced the spatial distributions of intermediate-age and old stars in nine dwarf galaxies in the distant parts of the Local Group, using multi-epoch 3.6 and 4.5 micron data from the DUST in Nearby Galaxies with Spitzer (DUSTiNGS) survey. Using complementary optical imaging from the Hubble Space Telescope, we identify the tip of the red giant branch (TRGB) in the 3.6 micron photometry, sepa… ▽ More We have traced the spatial distributions of intermediate-age and old stars in nine dwarf galaxies in the distant parts of the Local Group, using multi-epoch 3.6 and 4.5 micron data from the DUST in Nearby Galaxies with Spitzer (DUSTiNGS) survey. Using complementary optical imaging from the Hubble Space Telescope, we identify the tip of the red giant branch (TRGB) in the 3.6 micron photometry, separating thermally-pulsating asymptotic giant branch (TP-AGB) stars from the larger red giant branch (RGB) populations. Unlike the constant TRGB in the I-band, at 3.6 micron the TRGB magnitude varies by ~0.7 mag, making it unreliable as a distance indicator. The intermediate-age and old stars are well mixed in two-thirds of the sample with no evidence of a gradient in the ratio of the intermediate-age to old stellar populations outside the central ~1-2'. Variable AGB stars are detected in the outer extremities of the galaxies, indicating that chemical enrichment from these dust-producing stars may occur in the outer regions of galaxies with some frequency. Theories of structure formation in dwarf galaxies must account for the lack of radial gradients in intermediate-age populations and the presence of these stars in the outer extremities of dwarfs. Finally, we identify unique features in individual galaxies, such as extended tidal features in Sex A and Sag DIG and a central concentration of AGB stars in the inner regions of NGC 185 and NGC 147. △ Less

Submitted 15 November, 2016; originally announced November 2016.

Comments: 27 pages, 21 figures, 6 tables

Uncertainties in the attribution of greenhouse gas warming and implications for climate prediction

Authors: Gareth S. Jones, Peter A. Stott, John F. B. Mitchell

Abstract: Using optimal detection techniques with climate model simulations, most of the observed increase of near surface temperatures over the second half of the twentieth century is attributed to anthropogenic influences. However, the partitioning of the anthropogenic influence to individual factors, such as greenhouse gases and aerosols, is much less robust. Differences in how forcing factors are applie… ▽ More Using optimal detection techniques with climate model simulations, most of the observed increase of near surface temperatures over the second half of the twentieth century is attributed to anthropogenic influences. However, the partitioning of the anthropogenic influence to individual factors, such as greenhouse gases and aerosols, is much less robust. Differences in how forcing factors are applied, in their radiative influence and in models' climate sensitivities, substantially influence the response patterns. We find standard optimal detection methodologies cannot fully reconcile this response diversity. By selecting a set of experiments to enable the diagnosing of greenhouse gases and the combined influence of other anthropogenic and natural factors, we find robust detections of well mixed greenhouse gases across a large ensemble of models. Of the observed warming over the 20th century of 0.65K/century we find, using a multi model mean not incorporating pattern uncertainty, a well mixed greenhouse gas warming of 0.87 to 1.22K/century. This is partially offset by cooling from other anthropogenic and natural influences of -0.54 to -0.22K/century. Although better constrained than recent studies, the attributable trends across climate models are still wide, with implications for observational constrained estimates of transient climate response. Some of the uncertainties could be reduced in future by having more model data to better quantify the simulated estimates of the signals and natural variability, by designing model experiments more effectively and better quantification of the climate model radiative influences. Most importantly, how model pattern uncertainties are incorporated into the optimal detection methodology should be improved. △ Less

Submitted 16 June, 2016; originally announced June 2016.

Comments: 33 pages, 15 figures, and supporting information. Accepted for publication in Journal of Geophysical Research - Atmospheres, 2016