Showing 1–50 of 233 results for author: Ross, N

Exact Synthesis of Multiqutrit Clifford-Cyclotomic Circuits

Authors: Andrew N. Glaudell, Neil J. Ross, John van de Wetering, Lia Yeh

Abstract: It is known that the matrices that can be exactly represented by a multiqubit circuit over the Toffoli+Hadamard, Clifford+$T$, or, more generally, Clifford-cyclotomic gate set are precisely the unitary matrices with entries in the ring $\mathbb{Z}[1/2, ζ_k]$, where $k$ is a positive integer that depends on the gate set and $ζ_k$ is a primitive $2^k$-th root of unity. In the present paper, we estab… ▽ More It is known that the matrices that can be exactly represented by a multiqubit circuit over the Toffoli+Hadamard, Clifford+$T$, or, more generally, Clifford-cyclotomic gate set are precisely the unitary matrices with entries in the ring $\mathbb{Z}[1/2, ζ_k]$, where $k$ is a positive integer that depends on the gate set and $ζ_k$ is a primitive $2^k$-th root of unity. In the present paper, we establish an analogous correspondence for qutrits. We define the multiqutrit Clifford-cyclotomic gate set of degree $3^k$ by extending the classical qutrit gates $X$, $CX$, and $CCX$ with the Hadamard gate $H$ and the $T_k$ gate $T_k=\mathrm{diag}(1,ω_k, ω_k^2)$, where $ω_k$ is a primitive $3^k$-th root of unity. This gate set is equivalent to the qutrit Toffoli+Hadamard gate set when $k=1$, and to the qutrit Clifford+$T_k$ gate set when $k>1$. We then prove that a $3^n\times 3^n$ unitary matrix $U$ can be represented by an $n$-qutrit circuit over the Clifford-cyclotomic gate set of degree $3^k$ if and only if the entries of $U$ lie in the ring $\mathbb{Z}[1/3,ω_k]$. △ Less

Submitted 18 June, 2024; v1 submitted 13 May, 2024; originally announced May 2024.

Comments: v1: To appear in the proceedings of the 21st International Conference on Quantum Physics and Logic (QPL 2024). Only change to v2: Appendix figures now display. Only change to v3: Fix typo in metadata. Changes to v4: Minor changes to incorporate reviewer comments, edits to bibliography

Catalysing Completeness and Universality

Authors: Aleks Kissinger, Neil J. Ross, John van de Wetering

Abstract: A catalysis state is a quantum state that is used to make some desired operation possible or more efficient, while not being consumed in the process. Recent years have seen catalysis used in state-of-the-art protocols for implementing magic state distillation or small angle phase rotations. In this paper we will see that we can also use catalysis to prove that certain gate sets are computationally… ▽ More A catalysis state is a quantum state that is used to make some desired operation possible or more efficient, while not being consumed in the process. Recent years have seen catalysis used in state-of-the-art protocols for implementing magic state distillation or small angle phase rotations. In this paper we will see that we can also use catalysis to prove that certain gate sets are computationally universal, and to extend completeness results of graphical languages to larger fragments. In particular, we give a simple proof of the computational universality of the CS+Hadamard gate set using the catalysis of a $T$ gate using a CS gate, which sidesteps the more complicated analytic arguments of the original proof by Kitaev. This then also gives us a simple self-contained proof of the computational universality of Toffoli+Hadamard. Additionally, we show that the phase-free ZH-calculus can be extended to a larger complete fragment, just by using a single catalysis rule (and one scalar rule). △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 12+5 pages

Exact Synthesis of Multiqubit Clifford-Cyclotomic Circuits

Authors: Matthew Amy, Andrew N. Glaudell, Shaun Kelso, William Maxwell, Samuel S. Mendelson, Neil J. Ross

Abstract: Let $n\geq 8$ be divisible by 4. The Clifford-cyclotomic gate set $\mathcal{G}_n$ is the universal gate set obtained by extending the Clifford gates with the $z$-rotation $T_n = \mathrm{diag}(1,ζ_n)$, where $ζ_n$ is a primitive $n$-th root of unity. In this note, we show that, when $n$ is a power of 2, a multiqubit unitary matrix $U$ can be exactly represented by a circuit over $\mathcal{G}_n$ if… ▽ More Let $n\geq 8$ be divisible by 4. The Clifford-cyclotomic gate set $\mathcal{G}_n$ is the universal gate set obtained by extending the Clifford gates with the $z$-rotation $T_n = \mathrm{diag}(1,ζ_n)$, where $ζ_n$ is a primitive $n$-th root of unity. In this note, we show that, when $n$ is a power of 2, a multiqubit unitary matrix $U$ can be exactly represented by a circuit over $\mathcal{G}_n$ if and only if the entries of $U$ belong to the ring $\mathbb{Z}[1/2,ζ_n]$. We moreover show that $\log(n)-2$ ancillas are always sufficient to construct a circuit for $U$. Our results generalize prior work to an infinite family of gate sets and show that the limitations that apply to single-qubit unitaries, for which the correspondence between Clifford-cyclotomic operators and matrices over $\mathbb{Z}[1/2,ζ_n]$ fails for all but finitely many values of $n$, can be overcome through the use of ancillas. △ Less

Submitted 12 April, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Gaussian random field approximation via Stein's method with applications to wide random neural networks

Authors: Krishnakumar Balasubramanian, Larry Goldstein, Nathan Ross, Adil Salim

Abstract: We derive upper bounds on the Wasserstein distance ($W_1$), with respect to $\sup$-norm, between any continuous $\mathbb{R}^d$ valued random field indexed by the $n$-sphere and the Gaussian, based on Stein's method. We develop a novel Gaussian smoothing technique that allows us to transfer a bound in a smoother metric to the $W_1$ distance. The smoothing is based on covariance functions constructe… ▽ More We derive upper bounds on the Wasserstein distance ($W_1$), with respect to $\sup$-norm, between any continuous $\mathbb{R}^d$ valued random field indexed by the $n$-sphere and the Gaussian, based on Stein's method. We develop a novel Gaussian smoothing technique that allows us to transfer a bound in a smoother metric to the $W_1$ distance. The smoothing is based on covariance functions constructed using powers of Laplacian operators, designed so that the associated Gaussian process has a tractable Cameron-Martin or Reproducing Kernel Hilbert Space. This feature enables us to move beyond one dimensional interval-based index sets that were previously considered in the literature. Specializing our general result, we obtain the first bounds on the Gaussian random field approximation of wide random neural networks of any depth and Lipschitz activation functions at the random field level. Our bounds are explicitly expressed in terms of the widths of the network and moments of the random weights. We also obtain tighter bounds when the activation function has three bounded derivatives. △ Less

Submitted 30 April, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: To appear in Applied and Computational Harmonic Analysis

Improved Synthesis of Toffoli-Hadamard Circuits

Authors: Matthew Amy, Andrew N. Glaudell, Sarah Meng Li, Neil J. Ross

Abstract: The matrices that can be exactly represented by a circuit over the Toffoli-Hadamard gate set are the orthogonal matrices of the form $M/ \sqrt{2}{}^k$, where $M$ is an integer matrix and $k$ is a nonnegative integer. The exact synthesis problem for this gate set is the problem of constructing a circuit for a given such matrix. Existing methods produce circuits consisting of $O(2^n \log(n)k)$ gates… ▽ More The matrices that can be exactly represented by a circuit over the Toffoli-Hadamard gate set are the orthogonal matrices of the form $M/ \sqrt{2}{}^k$, where $M$ is an integer matrix and $k$ is a nonnegative integer. The exact synthesis problem for this gate set is the problem of constructing a circuit for a given such matrix. Existing methods produce circuits consisting of $O(2^n \log(n)k)$ gates, where $n$ is the dimension of the matrix. In this paper, we provide two improved synthesis methods. First, we show that a technique introduced by Kliuchnikov in 2013 for Clifford+$T$ circuits can be straightforwardly adapted to Toffoli-Hadamard circuits, reducing the complexity of the synthesized circuit from $O(2^n \log(n)k)$ to $O(n^2 \log(n)k)$. Then, we present an alternative synthesis method of similarly improved cost, but whose application is restricted to circuits on no more than three qubits. Our results also apply to orthogonal matrices over the dyadic fractions, which correspond to circuits using the 2-qubit gate $H\otimes H$, rather than the usual single-qubit Hadamard gate $H$. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Catalytic Embeddings of Quantum Circuits

Authors: M. Amy, M. Crawford, A. N. Glaudell, M. L. Macasieb, S. S. Mendelson, N. J. Ross

Abstract: If a set $\mathbb{G}$ of quantum gates is countable, then the operators that can be exactly represented by a circuit over $\mathbb{G}$ form a strict subset of the collection of all unitary operators. When $\mathbb{G}$ is universal, one circumvents this limitation by resorting to repeated gate approximations: every occurrence of a gate which cannot be exactly represented over $\mathbb{G}$ is replac… ▽ More If a set $\mathbb{G}$ of quantum gates is countable, then the operators that can be exactly represented by a circuit over $\mathbb{G}$ form a strict subset of the collection of all unitary operators. When $\mathbb{G}$ is universal, one circumvents this limitation by resorting to repeated gate approximations: every occurrence of a gate which cannot be exactly represented over $\mathbb{G}$ is replaced by an approximating circuit. Here, we introduce catalytic embeddings, which provide an alternative to repeated gate approximations. With catalytic embeddings, approximations are relegated to the preparation of a fixed number of reusable resource states called catalysts. Because the catalysts only need to be prepared once, when catalytic embeddings exist they always produce shorter circuits, in the limit of increasing gate count and target precision. In the present paper, we lay the foundations of the theory of catalytic embeddings and we establish several of their structural properties. In addition, we provide methods to design catalytic embeddings, showing that their construction can be reduced to that of a single fixed matrix when the gates involved have entries in well-behaved rings of algebraic numbers. Finally, we showcase some concrete examples and applications. Notably, we show that catalytic embeddings generalize a technique previously used to implement the Quantum Fourier Transform over the Clifford+$T$ gate set with $O(n)$ gate approximations. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Stochastic Recurrent Neural Networks for Modelling Astronomical Time Series: Advantages and Limitations

Authors: Xinyue Sheng, Matt Nicholl, Nicholas Ross

Abstract: This paper reviews the Stochastic Recurrent Neural Network (SRNN) as applied to the light curves of Active Galactic Nuclei by Sheng et al. (2022). Astronomical data have inherent limitations arising from telescope capabilities, cadence strategies, inevitable observing weather conditions, and current understanding of celestial objects. When applying machine learning methods, it is vital to understa… ▽ More This paper reviews the Stochastic Recurrent Neural Network (SRNN) as applied to the light curves of Active Galactic Nuclei by Sheng et al. (2022). Astronomical data have inherent limitations arising from telescope capabilities, cadence strategies, inevitable observing weather conditions, and current understanding of celestial objects. When applying machine learning methods, it is vital to understand the effects of data limitations on our analysis and ability to make inferences. We take Sheng et al. (2022) as a case study, and illustrate the problems and limitations encountered in implementing the SRNN for simulating AGN variability as seen by the Rubin Observatory. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: To appear in proceedings of EAS 2022 S11 session on machine learning in astronomy

A Model for Understanding and Reducing Developer Burnout

Authors: Bianca Trinkenreich, Klaas-Jan Stol, Igor Steinmacher, Marco Gerosa, Anita Sarma, Marcelo Lara, Michael Feathers, Nicholas Ross, Kevin Bishop

Abstract: Job burnout is a type of work-related stress associated with a state of physical or emotional exhaustion that also involves a sense of reduced accomplishment and loss of personal identity. Burnt out can affect one's physical and mental health and has become a leading industry concern and can result in high workforce turnover. Through an empirical study at Globant, a large multi-national company, w… ▽ More Job burnout is a type of work-related stress associated with a state of physical or emotional exhaustion that also involves a sense of reduced accomplishment and loss of personal identity. Burnt out can affect one's physical and mental health and has become a leading industry concern and can result in high workforce turnover. Through an empirical study at Globant, a large multi-national company, we created a theoretical model to evaluate the complex interplay among organizational culture, work satisfaction, and team climate, and how they impact developer burnout. We conducted a survey of developers in software delivery teams (n=3281) to test our model and analyzed the data using structural equation modeling, moderation, and multi-group analysis. Our results show that Organizational Culture, Climate for Learning, Sense of Belonging, and Inclusiveness are positively associated with Work Satisfaction, which in turn is associated with Reduced Burnout. Our model generated through a large-scale survey can guide organizations in how to reduce workforce burnout by creating a climate for learning, inclusiveness in teams, and a generative organizational culture where new ideas are welcome, information is actively sought and bad news can be shared without fear. △ Less

Submitted 24 January, 2023; v1 submitted 22 January, 2023; originally announced January 2023.

Journal ref: 45th IEEE/ACM International Conference on Software Engineering, Software Engineering in Practice (ICSE SEIP 2023)

Scalable fabrication of hemispherical solid immersion lenses in silicon carbide through grayscale hard-mask lithography

Authors: Christiaan Bekker, Muhammad Junaid Arshad, Pasquale Cilibrizzi, Charalampos Nikolatos, Peter Lomax, Graham S. Wood, Rebecca Cheung, Wolfgang Knolle, Neil Ross, Brian Gerardot, Cristian Bonato

Abstract: Grayscale lithography allows the creation of micrometer-scale features with spatially-controlled height in a process that is fully compatible with standard lithography. Here, solid immersion lenses are demonstrated in silicon carbide using a novel fabrication protocol combining grayscale lithography and hard-mask techniques to allow nearly hemispherical lenses of 5 $μ$m radius to be etched into th… ▽ More Grayscale lithography allows the creation of micrometer-scale features with spatially-controlled height in a process that is fully compatible with standard lithography. Here, solid immersion lenses are demonstrated in silicon carbide using a novel fabrication protocol combining grayscale lithography and hard-mask techniques to allow nearly hemispherical lenses of 5 $μ$m radius to be etched into the substrate. The technique is highly scalable and compatible with CMOS technology, and device aspect ratios can be tuned after resist patterning by controlling the chemistry of the subsequent dry etch. These results provide a low-cost, high-throughput and industrially-relevant alternative to focused ion beam milling for the creation of high-aspect-ratio, rounded microstructures. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 9 pages, 4 + 2 supplementary figures

Journal ref: Appl. Phys. Lett. 122, 173507 (2023)

On the Lambek embedding and the category of product-preserving presheaves

Authors: Peng Fu, Kohei Kishida, Neil J. Ross, Peter Selinger

Abstract: It is well-known that the category of presheaf functors is complete and cocomplete, and that the Yoneda embedding into the presheaf category preserves products. However, the Yoneda embedding does not preserve coproducts. It is perhaps less well-known that if we restrict the codomain of the Yoneda embedding to the full subcategory of limit-preserving functors, then this embedding preserves colimits… ▽ More It is well-known that the category of presheaf functors is complete and cocomplete, and that the Yoneda embedding into the presheaf category preserves products. However, the Yoneda embedding does not preserve coproducts. It is perhaps less well-known that if we restrict the codomain of the Yoneda embedding to the full subcategory of limit-preserving functors, then this embedding preserves colimits, while still enjoying most of the other useful properties of the Yoneda embedding. We call this modified embedding the Lambek embedding. The category of limit-preserving functors is known to be a reflective subcategory of the category of all functors, i.e., there is a left adjoint for the inclusion functor. In the literature, the existence of this left adjoint is often proved non-constructively, e.g., by an application of Freyd's adjoint functor theorem. In this paper, we provide an alternative, more constructive proof of this fact. We first explain the Lambek embedding and why it preserves coproducts. Then we review some concepts from multi-sorted algebras and observe that there is a one-to-one correspondence between product-preserving presheaves and certain multi-sorted term algebras. We provide a construction that freely turns any presheaf functor into a product-preserving one, hence giving an explicit definition of the left adjoint functor of the inclusion. Finally, we sketch how to extend our method to prove that the subcategory of limit-preserving functors is also reflective. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Symbolic Synthesis of Clifford Circuits and Beyond

Authors: Matthew Amy, Owen Bennett-Gibbs, Neil J. Ross

Abstract: Path sums are a convenient symbolic formalism for quantum operations with applications to the simulation, optimization, and verification of quantum protocols. Unlike quantum circuits, path sums are not limited to unitary operations, but can express arbitrary linear ones. Two problems, therefore, naturally arise in the study of path sums: the unitarity problem and the extraction problem. The former… ▽ More Path sums are a convenient symbolic formalism for quantum operations with applications to the simulation, optimization, and verification of quantum protocols. Unlike quantum circuits, path sums are not limited to unitary operations, but can express arbitrary linear ones. Two problems, therefore, naturally arise in the study of path sums: the unitarity problem and the extraction problem. The former is the problem of deciding whether a given path sum represents a unitary operator. The latter is the problem of constructing a quantum circuit, given a path sum promised to represent a unitary operator. In this paper, we show that the unitarity problem is co-NP-hard in general, but that it is in P when restricted to Clifford path sums. We then provide an algorithm to synthesize a Clifford circuit from a unitary Clifford path sum. The circuits produced by our extraction algorithm are of the form C1-H-C2, where C1 and C2 are Hadamard-free circuits and H is a layer of Hadamard gates. We also provide a heuristic generalization of our extraction algorithm to arbitrary path sums. While this algorithm is not guaranteed to succeed, it often succeeds and typically produces natural looking circuits. Alongside applications to the optimization and decompilation of quantum circuits, we demonstrate the capability of our algorithm by synthesizing the standard quantum Fourier transform directly from a path sum. △ Less

Submitted 15 November, 2023; v1 submitted 29 April, 2022; originally announced April 2022.

Comments: In Proceedings QPL 2022, arXiv:2311.08375

Journal ref: EPTCS 394, 2023, pp. 343-362

Proto-Quipper with dynamic lifting

Authors: Peng Fu, Kohei Kishida, Neil J. Ross, Peter Selinger

Abstract: Quipper is a functional programming language for quantum computing. Proto-Quipper is a family of languages aiming to provide a formal foundation for Quipper. In this paper, we extend Proto-Quipper-M with a construct called dynamic lifting, which is present in Quipper. By virtue of being a circuit description language, Proto-Quipper has two separate runtimes: circuit generation time and circuit exe… ▽ More Quipper is a functional programming language for quantum computing. Proto-Quipper is a family of languages aiming to provide a formal foundation for Quipper. In this paper, we extend Proto-Quipper-M with a construct called dynamic lifting, which is present in Quipper. By virtue of being a circuit description language, Proto-Quipper has two separate runtimes: circuit generation time and circuit execution time. Values that are known at circuit generation time are called parameters, and values that are known at circuit execution time are called states. Dynamic lifting is an operation that enables a state, such as the result of a measurement, to be lifted to a parameter, where it can influence the generation of the next portion of the circuit. As a result, dynamic lifting enables Proto-Quipper programs to interleave classical and quantum computation. We describe the syntax of a language we call Proto-Quipper-Dyn. Its type system uses a system of modalities to keep track of the use of dynamic lifting. We also provide an operational semantics, as well as an abstract categorical semantics for dynamic lifting based on enriched category theory. We prove that both the type system and the operational semantics are sound with respect to our categorical semantics. Finally, we give some examples of Proto-Quipper-Dyn programs that make essential use of dynamic lifting. △ Less

Submitted 8 November, 2022; v1 submitted 27 April, 2022; originally announced April 2022.

A Biset-Enriched Categorical Model for Proto-Quipper with Dynamic Lifting

Authors: Peng Fu, Kohei Kishida, Neil J. Ross, Peter Selinger

Abstract: Quipper and Proto-Quipper are a family of quantum programming languages that, by their nature as circuit description languages, involve two runtimes: one at which the program generates a circuit and one at which the circuit is executed, normally with probabilistic results due to measurements. Accordingly, the language distinguishes two kinds of data: parameters, which are known at circuit generati… ▽ More Quipper and Proto-Quipper are a family of quantum programming languages that, by their nature as circuit description languages, involve two runtimes: one at which the program generates a circuit and one at which the circuit is executed, normally with probabilistic results due to measurements. Accordingly, the language distinguishes two kinds of data: parameters, which are known at circuit generation time, and states, which are known at circuit execution time. Sometimes, it is desirable for the results of measurements to control the generation of the next part of the circuit. Therefore, the language needs to turn states, such as measurement outcomes, into parameters, an operation we call dynamic lifting. The goal of this paper is to model this interaction between the runtimes by providing a general categorical structure enriched in what we call "bisets". We demonstrate that the biset-enriched structure achieves a proper semantics of the two runtimes and their interaction, by showing that it models a variant of Proto-Quipper with dynamic lifting. The present paper deals with the concrete categorical semantics of this language, whereas a companion paper deals with the syntax, type system, operational semantics, and abstract categorical semantics. △ Less

Submitted 15 November, 2023; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: In Proceedings QPL 2022, arXiv:2311.08375

Journal ref: EPTCS 394, 2023, pp. 302-342

LSST Cadence Strategy Evaluations for AGN Time-series Data in Wide-Fast-Deep Field

Authors: Xinyue Sheng, Nicholas Ross, Matt Nicholl

Abstract: Machine learning is a promising tool to reconstruct time-series phenomena, such as variability of active galactic nuclei (AGN), from sparsely-sampled data. Here we use three Continuous Auto-Regressive Moving Average (CARMA) representations of AGN variability -- the Damped Random Walk (DRW) and (over/under-)Damped Harmonic Oscillator (DHO) -- to simulate 10-year AGN light curves as they would appea… ▽ More Machine learning is a promising tool to reconstruct time-series phenomena, such as variability of active galactic nuclei (AGN), from sparsely-sampled data. Here we use three Continuous Auto-Regressive Moving Average (CARMA) representations of AGN variability -- the Damped Random Walk (DRW) and (over/under-)Damped Harmonic Oscillator (DHO) -- to simulate 10-year AGN light curves as they would appear in the upcoming Vera Rubin Observatory Legacy Survey of Space and Time (LSST), and provide a public tool to generate these for any survey cadence. We investigate the impact on AGN science of five proposed cadence strategies for LSST's primary Wide-Fast-Deep (WFD) survey. We apply for the first time in astronomy a novel Stochastic Recurrent Neural Network (SRNN) algorithm to reconstruct input light curves from the simulated LSST data, and provide a metric to evaluate how well SRNN can help recover the underlying CARMA parameters. We find that the light curve reconstruction is most sensitive to the duration of gaps between observing season, and that of the proposed cadences, those that change the balance between filters, or avoid having long gaps in the {g}-band perform better. Overall, SRNN is a promising means to reconstruct densely sampled AGN light curves and recover the long-term Structure Function of the DRW process (SF$_\infty$) reasonably well. However, we find that for all cadences, CARMA/SRNN models struggle to recover the decorrelation timescale ($τ$) due to the long gaps in survey observations. This may indicate a major limitation in using LSST WFD data for AGN variability science. △ Less

Submitted 19 April, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: accepted by MNRAS

Qutrit metaplectic gates are a subset of Clifford+T

Authors: Andrew Glaudell, Neil J. Ross, John van de Wetering, Lia Yeh

Abstract: A popular universal gate set for quantum computing with qubits is Clifford+T, as this can be readily implemented on many fault-tolerant architectures. For qutrits, there is an equivalent T gate, that, like its qubit analogue, makes Clifford+T approximately universal, is injectable by a magic state, and supports magic state distillation. However, it was claimed that a better gate set for qutrits mi… ▽ More A popular universal gate set for quantum computing with qubits is Clifford+T, as this can be readily implemented on many fault-tolerant architectures. For qutrits, there is an equivalent T gate, that, like its qubit analogue, makes Clifford+T approximately universal, is injectable by a magic state, and supports magic state distillation. However, it was claimed that a better gate set for qutrits might be Clifford+R, where R=diag(1,1,-1) is the metaplectic gate, as certain protocols and gates could more easily be implemented using the R gate than the T gate. In this paper we show that when we have at least two qutrits, the qutrit Clifford+R unitaries form a strict subset of the Clifford+T unitaries, by finding a direct decomposition of $R \otimes \mathbb{I}$ as a Clifford+T circuit and proving that the T gate cannot be exactly synthesized in Clifford+R. This shows that in fact the T gate is at least as powerful as the R gate, up to a constant factor. Moreover, we additionally show that it is impossible to find a single-qutrit Clifford+T decomposition of the R gate, making our result tight. △ Less

Submitted 21 February, 2022; v1 submitted 18 February, 2022; originally announced February 2022.

Comments: 10 pages + 4 page appendix

Journal ref: Proceedings TQC 2022. pp. 12:1--12:15

Stein's method, Gaussian processes and Palm measures, with applications to queueing

Authors: A. D. Barbour, Nathan Ross, Guangqu Zheng

Abstract: We develop a general approach to Stein's method for approximating a random process in the path space $D([0,T]\to R^d)$ by a real continuous Gaussian process. We then use the approach in the context of processes that have a representation as integrals with respect to anunderlying point process, deriving a general quantitative Gaussian approximation. The error bound is expressed in terms of coupling… ▽ More We develop a general approach to Stein's method for approximating a random process in the path space $D([0,T]\to R^d)$ by a real continuous Gaussian process. We then use the approach in the context of processes that have a representation as integrals with respect to anunderlying point process, deriving a general quantitative Gaussian approximation. The error bound is expressed in terms of couplings of the original process to processes generated from the reduced Palm measures associated with the point process. As applications, we study certain $\text{GI}/\text{GI}/\infty$ queues in the "heavy traffic" regime. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: 40 pages

MSC Class: Primary 60G15; 60G55; secondary 60K25; 60F25

Journal ref: Ann. Appl. Probab. 33(5): 3835-3871 (2023)

Generators and Relations for Real Stabilizer Operators

Authors: Justin Makary, Neil J. Ross, Peter Selinger

Abstract: Real stabilizer operators, which are also known as real Clifford operators, are generated, through composition and tensor product, by the Hadamard gate, the Pauli Z gate, and the controlled-Z gate. We introduce a normal form for real stabilizer circuits and show that every real stabilizer operator admits a unique normal form. Moreover, we give a finite set of relations that suffice to rewrite any… ▽ More Real stabilizer operators, which are also known as real Clifford operators, are generated, through composition and tensor product, by the Hadamard gate, the Pauli Z gate, and the controlled-Z gate. We introduce a normal form for real stabilizer circuits and show that every real stabilizer operator admits a unique normal form. Moreover, we give a finite set of relations that suffice to rewrite any real stabilizer circuit to its normal form. △ Less

Submitted 12 September, 2021; originally announced September 2021.

Comments: In Proceedings QPL 2021, arXiv:2109.04886

Journal ref: EPTCS 343, 2021, pp. 14-36

Understanding the roles of electronic effect in CO on Pt-Sn alloy surface via band structure measurements

Authors: Jongkeun Jung, Sungwoo Kang Laurent Nicolai, Jisook Hong, Jan Minár, Inkyung Song, Wonshik Kyung, Soohyun Cho, Beomseo Kim, Jonathan D. Denlinger, Francisco J. C. S. Aires, Eric Ehret, Philip N. Ross, Jihoon Shim, Slavomir Nemšák, Doyoung Noh, Seungwu Han, Changyoung Kim, Bongjin S. Mun

Abstract: Using angle-resolved photoemission spectroscopy, we show the direct evidence of charge transfer between adsorbed molecules and metal substrate, i.e. chemisorption of CO on Pt(111) and Pt-Sn/Pt(111) 2x2 surfaces. The observed band structure shows a unique signature of charge transfer as CO atoms are adsorbed,revealing the roles of specific orbital characters participating in the chemisorption proce… ▽ More Using angle-resolved photoemission spectroscopy, we show the direct evidence of charge transfer between adsorbed molecules and metal substrate, i.e. chemisorption of CO on Pt(111) and Pt-Sn/Pt(111) 2x2 surfaces. The observed band structure shows a unique signature of charge transfer as CO atoms are adsorbed,revealing the roles of specific orbital characters participating in the chemisorption process. As the coverage of CO increases, the degree of charge transfer between CO and Pt shows clear difference to that of Pt-Sn. With comparison to DFT calculation results, the observed distinct features in the band structure are interpreted as backdonation bonding states of Pt molecular orbital to the 2π orbital of CO. Furthermore, the change in the surface charge concentration, measured from the Fermi surface area, shows Pt surface has a larger charge concentration change than Pt-Sn surface upon CO adsorption. The difference in the charge concentration change between Pt and Pt-Sn surfaces reflects the degree of electronic effects during CO adsorption on Pt-Sn. △ Less

Submitted 9 August, 2021; originally announced August 2021.

Comments: 15 pages, 9 figures, 2 tables

Stein's method, smoothing and functional approximation

Authors: A. D. Barbour, Nathan Ross, Guangqu Zheng

Abstract: Stein's method for Gaussian process approximation can be used to bound the differences between the expectations of smooth functionals $h$ of a càdlàg random process $X$ of interest and the expectations of the same functionals of a well understood target random process $Z$ with continuous paths. Unfortunately, the class of smooth functionals for which this is easily possible is very restricted. Her… ▽ More Stein's method for Gaussian process approximation can be used to bound the differences between the expectations of smooth functionals $h$ of a càdlàg random process $X$ of interest and the expectations of the same functionals of a well understood target random process $Z$ with continuous paths. Unfortunately, the class of smooth functionals for which this is easily possible is very restricted. Here, we prove an infinite dimensional Gaussian smoothing inequality, which enables the class of functionals to be greatly expanded -- examples are Lipschitz functionals with respect to the uniform metric, and indicators of arbitrary events -- in exchange for a loss of precision in the bounds. Our inequalities are expressed in terms of the smooth test function bound, an expectation of a functional of $X$ that is closely related to classical tightness criteria, a similar expectation for $Z$, and, for the indicator of a set $K$, the probability $\mathbb{P}(Z \in K^θ\setminus K^{-θ})$ that the target process is close to the boundary of $K$. △ Less

Submitted 4 September, 2023; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: Ver 4: 33 pages, added Example 1.10, improved some details and discussion, and streamlined some notation; Ver3: 28 pages, corrected a mistake in Lemma 1.10 and added discussion and details, only superficial changes to the main result; Ver2: 21 pages, additional discussion and details; Ver1: 17 pages

Journal ref: Electronic Journal of Probability 2024, Vol. 29, paper no. 20, 1-29

Generators and Relations for the Group On(Z[1/2])

Authors: Sarah Meng Li, Neil J. Ross, Peter Selinger

Abstract: We give a finite presentation by generators and relations for the group O_n(Z[1/2]) of n-dimensional orthogonal matrices with entries in Z[1/2]. We then obtain a similar presentation for the group of n-dimensional orthogonal matrices of the form M/sqrt(2)^k, where k is a nonnegative integer and M is an integer matrix. Both groups arise in the study of quantum circuits. In particular, when the dime… ▽ More We give a finite presentation by generators and relations for the group O_n(Z[1/2]) of n-dimensional orthogonal matrices with entries in Z[1/2]. We then obtain a similar presentation for the group of n-dimensional orthogonal matrices of the form M/sqrt(2)^k, where k is a nonnegative integer and M is an integer matrix. Both groups arise in the study of quantum circuits. In particular, when the dimension is a power of 2, the elements of the latter group are precisely the unitary matrices that can be represented by a quantum circuit over the universal gate set consisting of the Toffoli gate, the Hadamard gate, and the computational ancilla. △ Less

Submitted 12 September, 2021; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: In Proceedings QPL 2021, arXiv:2109.04886

Journal ref: EPTCS 343, 2021, pp. 210-264

The phase/state duality in reversible circuit design

Authors: Matthew Amy, Neil J. Ross

Abstract: The reversible implementation of classical functions accounts for the bulk of most known quantum algorithms. As a result, a number of reversible circuit constructions over the Clifford+$T$ gate set have been developed in recent years which use both the state and phase spaces, or $X$ and $Z$ bases, to reduce circuit costs beyond what is possible at the strictly classical level. We study and general… ▽ More The reversible implementation of classical functions accounts for the bulk of most known quantum algorithms. As a result, a number of reversible circuit constructions over the Clifford+$T$ gate set have been developed in recent years which use both the state and phase spaces, or $X$ and $Z$ bases, to reduce circuit costs beyond what is possible at the strictly classical level. We study and generalize two particular classes of these constructions: relative phase circuits, including Giles and Selinger's multiply-controlled $iX$ gates and Maslov's $4$ qubit Toffoli gate, and measurement-assisted circuits, including Jones' Toffoli gate and Gidney's temporary logical-AND. In doing so, we introduce general methods for implementing classical functions up to phase and for measurement-assisted termination of temporary values. We then apply these techniques to find novel $T$-count efficient constructions of some classical functions in space-constrained regimes, notably multiply-controlled Toffoli gates and temporary products. △ Less

Submitted 5 November, 2021; v1 submitted 27 May, 2021; originally announced May 2021.

Comments: 16 pages, comments welcome

Journal ref: Phys. Rev. A 104, 052602 (2021)

A Systematic Survey for z < 0.04 Changing-Look AGNs

Authors: Madhooshi R. Senarath, Michael J. I. Brown, Michelle E. Cluver, Thomas H. Jarrett, Christian Wolf, Nicholas P. Ross, John R. Lucey, Vaishali Parkash, Wei J. Hon

Abstract: We have conducted a systematic survey for z $<$ 0.04 active Galactic nuclei (AGNs) that may have changed spectral class over the past decade. We use SkyMapper, Pan-STARRS and the Véron-Cetty & Véron (2010) catalogue to search the entire sky for these ``changing-look'' AGNs using a variety of selection methods, where Pan-STARRS has a coverage of 3$π$ steradians (sky north of Declination… ▽ More We have conducted a systematic survey for z $<$ 0.04 active Galactic nuclei (AGNs) that may have changed spectral class over the past decade. We use SkyMapper, Pan-STARRS and the Véron-Cetty & Véron (2010) catalogue to search the entire sky for these ``changing-look'' AGNs using a variety of selection methods, where Pan-STARRS has a coverage of 3$π$ steradians (sky north of Declination $-30^\circ$) and SkyMapper has coverage of $\sim$ 21,000$~\rm{deg^2}$ (sky south of Declination $0^\circ$). We use small aperture photometry to measure how colour and flux have changed over time, where a change may indicate a change in spectral type. Optical colour and flux are used as a proxy for changing H$α$ equivalent width, while WISE 3.4 $μ$m flux is used to look for changes in the hot dust component. We have identified four AGNs with varying spectra selected using our optical colour selection method. Three AGNs were confirmed from recent observations with WiFeS on the 2.3 m telescope at Siding Spring and the other was identified from archival spectra alone. From this, we identify two new changing look AGNs; NGC 1346 and 2MASX J20075129-1108346. We also recover Mrk 915 and Mrk 609, which are known to have varying spectra in the literature, but they do not meet our specific criteria for changing look AGNs. △ Less

Submitted 15 February, 2021; originally announced February 2021.

Comments: Accepted for publication in the Monthly Notices of the Royal Astronomical Society; 15 pages, 9 figures and 4 tables

Characterising the shape, size and orientation of cloud-feeding coherent boundary layer structures

Authors: Leif Denby, Steven J. Böing, Douglas J. Parker, Andrew N. Ross, Steven M. Tobias

Abstract: This paper presents two techniques for characterisation of cloud-feeding coherent boundary layer structures through analysis of large-eddy simulations of shallow cumulus clouds, contrasting conditions with and without ambient shear. The first technique is a generalisation of the two-point correlation function where the correlation length-scale as well as orientation can be extracted. The second te… ▽ More This paper presents two techniques for characterisation of cloud-feeding coherent boundary layer structures through analysis of large-eddy simulations of shallow cumulus clouds, contrasting conditions with and without ambient shear. The first technique is a generalisation of the two-point correlation function where the correlation length-scale as well as orientation can be extracted. The second technique decomposes the vertical transport by coherent structures by the shape, size and orientation of these structures. It is found that the structures dominating the vertical flux are plume-like in character (extending from the surface into cloud), show small width/thickness asymmetry and rise near-vertically in the absence of ambient wind. The planar stretching and tilting of boundary layer structures caused by the introduction of ambient shear is also quantified, demonstrating the general applicability of the techniques for future study of other boundary layer patterns. △ Less

Submitted 17 May, 2021; v1 submitted 17 August, 2020; originally announced August 2020.

Comments: 26 pages, 14 figures, submitted to QJRMS

Quasars at intermediate redshift are not special; but they are often satellites

Authors: Shadab Alam, Nicholas P. Ross, Sarah Eftekharzadeh, John A. Peacock, Johan Comparat, Adam D. Myers, Ashley J. Ross

Abstract: Understanding the links between the activity of supermassive black holes (SMBH) at the centres of galaxies and their host dark matter haloes is a key question in modern astrophysics. The final data release of the SDSS-IV eBOSS provides the largest contemporary spectroscopic sample of galaxies and QSOs. Using this sample and covering the redshift interval $z=0.7-1.1$, we have measured the clusterin… ▽ More Understanding the links between the activity of supermassive black holes (SMBH) at the centres of galaxies and their host dark matter haloes is a key question in modern astrophysics. The final data release of the SDSS-IV eBOSS provides the largest contemporary spectroscopic sample of galaxies and QSOs. Using this sample and covering the redshift interval $z=0.7-1.1$, we have measured the clustering properties of the eBOSS QSOs, Emission Line Galaxies (ELGs) and Luminous Red Galaxies (LRGs). We have also measured the fraction of QSOs as a function of the overdensity defined by the galaxy population. Using these measurements, we investigate how QSOs populate and sample the galaxy population, and how the host dark-matter haloes of QSOs sample the underlying halo distribution. We find that the probability of a galaxy hosting a QSO is independent of the host dark matter halo mass of the galaxy. We also find that about 60\% of eBOSS QSOs are hosted by LRGs and about 20-40\% of QSOs are hosted by satellite galaxies. We find a slight preference for QSOs to populate satellite galaxies over central galaxies. This is connected to the host halo mass distribution of different types of galaxies. Based on our analysis, QSOs should be hosted by a very broad distribution of haloes, and their occurrence should be modulated only by the efficiency of galaxy formation processes. △ Less

Submitted 24 September, 2021; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: 15 pages, 6 Figures, Data and code from this analysis is available at https://www.roe.ac.uk/~salam/GQOD/ , Final published version with minor changes in text

Candidate Electromagnetic Counterpart to the Binary Black Hole Merger Gravitational Wave Event S190521g

Authors: M. J. Graham, K. E. S. Ford, B. McKernan, N. P. Ross, D. Stern, K. Burdge, M. Coughlin, S. G. Djorgovski, A. J. Drake, D. Duev, M. Kasliwal, A. A. Mahabal, S. van Velzen, J. Belicki, E. C. Bellm, R. Burruss, S. B. Cenko, V. Cunningham, G. Helou, S. R. Kulkarni, F. J. Masci, T. Prince, D. Reiley, H. Rodriguez, B. Rusholme , et al. (2 additional authors not shown)

Abstract: We report the first plausible optical electromagnetic (EM) counterpart to a (candidate) binary black hole (BBH) merger. Detected by the Zwicky Transient Facility (ZTF), the EM flare is consistent with expectations for a kicked BBH merger in the accretion disk of an active galactic nucleus (AGN), and is unlikely ($<O(0.01\%$)) due to intrinsic variability of this source. The lack of color evolution… ▽ More We report the first plausible optical electromagnetic (EM) counterpart to a (candidate) binary black hole (BBH) merger. Detected by the Zwicky Transient Facility (ZTF), the EM flare is consistent with expectations for a kicked BBH merger in the accretion disk of an active galactic nucleus (AGN), and is unlikely ($<O(0.01\%$)) due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernovae and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing or a tidal disruption event, are ruled out or constrained to be $<O(0.1\%$). If the flare is associated with S190521g, we find plausible values of: total mass $ M_{\rm BBH} \sim 100 M_{\odot}$, kick velocity $v_k \sim 200\, {\rm km}\, {\rm s}^{-1}$ at $θ\sim 60^{\circ}$ in a disk with aspect ratio $H/a \sim 0.01$ (i.e., disk height $H$ at radius $a$) and gas density $ρ\sim 10^{-10}\, {\rm g}\, {\rm cm}^{-3}$. The merger could have occurred at a disk migration trap ($a \sim 700\, r_{g}$; $r_g \equiv G M_{\rm SMBH} / c^2$, where $M_{\rm SMBH}$ is the mass of the AGN supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a re-encountering with the disk in $\sim 1.6\, {\rm yr}\, (M_{\rm SMBH}/10^{8}M_{\odot})\, (a/10^{3}r_{g})^{3/2}$. △ Less

Submitted 24 June, 2020; originally announced June 2020.

Comments: 9 pages, 4 figures, accepted for publication in Physical Review Letters (June 25, 2020)

Journal ref: Phys. Rev. Lett. 124, 251102

An outflow powers the optical rise of the nearby, fast-evolving tidal disruption event AT2019qiz

Authors: M. Nicholl, T. Wevers, S. R. Oates, K. D. Alexander, G. Leloudas, F. Onori, A. Jerkstrand, S. Gomez, S. Campana, I. Arcavi, P. Charalampopoulos, M. Gromadzki, N. Ihanec, P. G. Jonker, A. Lawrence, I. Mandel, S. Schulze, P. Short, J. Burke, C. McCully, D. Hiramatsu, D. A. Howell, C. Pellegrino, H. Abbot, J. P. Anderson , et al. (20 additional authors not shown)

Abstract: At 66 Mpc, AT2019qiz is the closest optical tidal disruption event (TDE) to date, with a luminosity intermediate between the bulk of the population and iPTF16fnl. Its proximity allowed a very early detection and triggering of multiwavelength and spectroscopic follow-up well before maximum light. The velocity dispersion of the host galaxy and fits to the TDE light curve indicate a black hole mass… ▽ More At 66 Mpc, AT2019qiz is the closest optical tidal disruption event (TDE) to date, with a luminosity intermediate between the bulk of the population and iPTF16fnl. Its proximity allowed a very early detection and triggering of multiwavelength and spectroscopic follow-up well before maximum light. The velocity dispersion of the host galaxy and fits to the TDE light curve indicate a black hole mass $\approx 10^6$ M$_\odot$, disrupting a star of $\approx 1$ M$_\odot$. Comprehensive UV, optical and X-ray data shows that the early optical emission is dominated by an outflow, with a luminosity evolution $L \propto t^2$, consistent with a photosphere expanding at constant velocity ($\gtrsim 2000$ km s$^{-1}$), and a line-forming region producing initially blueshifted H and He II profiles with $v=3000-10000$ km s$^{-1}$. The fastest optical ejecta approach the velocity inferred from radio detections (modelled in a forthcoming companion paper from K.~D.~Alexander et al.), thus the same outflow may be responsible for both the fast optical rise and the radio emission -- the first time this connection has been observed in a TDE. The light curve rise begins $29 \pm 2$ days before maximum light, peaking when the photosphere reaches the radius where optical photons can escape. The photosphere then undergoes a sudden transition, first cooling at constant radius then contracting at constant temperature. At the same time, the blueshifts disappear from the spectrum and Bowen fluorescence lines (N III) become prominent, implying a source of far-UV photons, while the X-ray light curve peaks at $\approx 10^{41}$ erg s$^{-1}$. Assuming that these X-rays are from prompt accretion, the size and mass of the outflow are consistent with the reprocessing layer needed to explain the large optical to X-ray ratio in this and other optical TDEs, possibly favouring accretion-powered over collision-powered outflow models. △ Less

Submitted 14 September, 2020; v1 submitted 3 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS

A tutorial introduction to quantum circuit programming in dependently typed Proto-Quipper

Authors: Peng Fu, Kohei Kishida, Neil J. Ross, Peter Selinger

Abstract: We introduce dependently typed Proto-Quipper, or Proto-Quipper-D for short, an experimental quantum circuit programming language with linear dependent types. We give several examples to illustrate how linear dependent types can help in the construction of correct quantum circuits. Specifically, we show how dependent types enable programming families of circuits, and how dependent types solve the p… ▽ More We introduce dependently typed Proto-Quipper, or Proto-Quipper-D for short, an experimental quantum circuit programming language with linear dependent types. We give several examples to illustrate how linear dependent types can help in the construction of correct quantum circuits. Specifically, we show how dependent types enable programming families of circuits, and how dependent types solve the problem of type-safe uncomputation of garbage qubits. We also discuss other language features along the way. △ Less

Submitted 12 December, 2020; v1 submitted 17 May, 2020; originally announced May 2020.

Comments: Added a section on related work and a paragraph explaining qubit initialization and termination

Journal ref: LNCS 12227:153-168 (2020)

The Rockerverse: Packages and Applications for Containerization with R

Authors: Daniel Nüst, Dirk Eddelbuettel, Dom Bennett, Robrecht Cannoodt, Dav Clark, Gergely Daroczi, Mark Edmondson, Colin Fay, Ellis Hughes, Lars Kjeldgaard, Sean Lopp, Ben Marwick, Heather Nolis, Jacqueline Nolis, Hong Ooi, Karthik Ram, Noam Ross, Lori Shepherd, Péter Sólymos, Tyson Lee Swetnam, Nitesh Turaga, Charlotte Van Petegem, Jason Williams, Craig Willis, Nan Xiao

Abstract: The Rocker Project provides widely used Docker images for R across different application scenarios. This article surveys downstream projects that build upon the Rocker Project images and presents the current state of R packages for managing Docker images and controlling containers. These use cases cover diverse topics such as package development, reproducible research, collaborative work, cloud-ba… ▽ More The Rocker Project provides widely used Docker images for R across different application scenarios. This article surveys downstream projects that build upon the Rocker Project images and presents the current state of R packages for managing Docker images and controlling containers. These use cases cover diverse topics such as package development, reproducible research, collaborative work, cloud-based data processing, and production deployment of services. The variety of applications demonstrates the power of the Rocker Project specifically and containerisation in general. Across the diverse ways to use containers, we identified common themes: reproducible environments, scalability and efficiency, and portability across clouds. We conclude that the current growth and diversification of use cases is likely to continue its positive impact, but see the need for consolidating the Rockerverse ecosystem of packages, developing common practices for applications, and exploring alternative containerisation software. △ Less

Submitted 17 August, 2020; v1 submitted 28 January, 2020; originally announced January 2020.

Comments: Source code for article available at https://github.com/nuest/rockerverse-paper/ Updated version includes some new paragraphs and corrections throughout the text; full diff available at https://github.com/nuest/rockerverse-paper/compare/preprint.v2...preprint.v3

MSC Class: 68N01 ACM Class: D.2.6; D.2.7; K.6.3

Journal ref: The R Journal (2020), 12:1, pages 437-461

Optimal Two-Qubit Circuits for Universal Fault-Tolerant Quantum Computation

Authors: Andrew N. Glaudell, Neil J. Ross, Jacob M. Taylor

Abstract: We study two-qubit circuits over the Clifford+CS gate set, which consists of the Clifford gates together with the controlled-phase gate CS=diag(1,1,1,i). The Clifford+CS gate set is universal for quantum computation and its elements can be implemented fault-tolerantly in most error-correcting schemes through magic state distillation. Since non-Clifford gates are typically more expensive to perform… ▽ More We study two-qubit circuits over the Clifford+CS gate set, which consists of the Clifford gates together with the controlled-phase gate CS=diag(1,1,1,i). The Clifford+CS gate set is universal for quantum computation and its elements can be implemented fault-tolerantly in most error-correcting schemes through magic state distillation. Since non-Clifford gates are typically more expensive to perform in a fault-tolerant manner, it is often desirable to construct circuits that use few CS gates. In the present paper, we introduce an efficient and optimal synthesis algorithm for two-qubit Clifford+CS operators. Our algorithm inputs a Clifford+CS operator U and outputs a Clifford+CS circuit for U, which uses the least possible number of CS gates. Because the algorithm is deterministic, the circuit it associates to a Clifford+CS operator can be viewed as a normal form for that operator. We give an explicit description of these normal forms and use this description to derive a worst-case lower bound of 5log(1/epsilon)+O(1) on the number of CS gates required to epsilon-approximate elements of SU(4). Our work leverages a wide variety of mathematical tools that may find further applications in the study of fault-tolerant quantum circuits. △ Less

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

Comments: 23 pages. Significant revisions to Sections 5 and 6. Updated introduction and conclusion. Link to open access implementation provided

The Bolometric Quasar Luminosity Function at z = 0-7

Authors: Xuejian Shen, Philip F. Hopkins, Claude-André Faucher-Giguère, D. M. Alexander, Gordon T. Richards, Nicholas P. Ross, R. C. Hickox

Abstract: In this paper, we provide updated constraints on the bolometric quasar luminosity function (QLF) from $z=0$ to $z=7$. The constraints are based on an observational compilation that includes observations in the rest-frame IR, B band, UV, soft and hard X-ray in past decades. Our method follows Hopkins et al. 2007 with an updated quasar SED model and bolometric and extinction corrections. The new bes… ▽ More In this paper, we provide updated constraints on the bolometric quasar luminosity function (QLF) from $z=0$ to $z=7$. The constraints are based on an observational compilation that includes observations in the rest-frame IR, B band, UV, soft and hard X-ray in past decades. Our method follows Hopkins et al. 2007 with an updated quasar SED model and bolometric and extinction corrections. The new best-fit bolometric quasar luminosity function behaves qualitatively different from the Hopkins et al. 2007 model at high redshift. Compared with the old model, the number density normalization decreases towards higher redshift and the bright-end slope is steeper at $z\gtrsim 2$. Due to the paucity of measurements at the faint end, the faint end slope at $z\gtrsim 5$ is quite uncertain. We present two models, one featuring a progressively steeper faint-end slope at higher redshift and the other featuring a shallow faint-end slope at $z\gtrsim 5$. Further multi-band observations of the faint-end QLF are needed to distinguish between these models. The evolutionary pattern of the bolometric QLF can be interpreted as an early phase likely dominated by the hierarchical assembly of structures and a late phase likely dominated by the quenching of galaxies. We explore the implications of this model on the ionizing photon production by quasars, the CXB spectrum, the SMBH mass density and mass functions. The predicted hydrogen photoionization rate contributed by quasars is subdominant during the epoch of reionization and only becomes important at $z\lesssim 3$. The predicted CXB spectrum, cosmic SMBH mass density and SMBH mass function are generally consistent with existing observations. △ Less

Submitted 15 May, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: Accepted by MNRAS. The code and data are available at https://bitbucket.org/ShenXuejian/quasarlf/src/master/

Through Andromeda and Beyond: AGN Optical Transient 'Sharov21' Revisited

Authors: A. Bruce, A. Lawrence, D. McLeod, Nicholas P. Ross

Abstract: We revisit a notable AGN known as `Sharov21', seen to undergo a dramatic outburst in 1992, brightening by a factor of thirty over a period of approximately one year. A simple microlensing model fit to the event lightcurve provides a constraint on the distance of the lensing object which is consistent with the distance to M31, strongly suggesting that this is the correct explanation. Archival XMM/H… ▽ More We revisit a notable AGN known as `Sharov21', seen to undergo a dramatic outburst in 1992, brightening by a factor of thirty over a period of approximately one year. A simple microlensing model fit to the event lightcurve provides a constraint on the distance of the lensing object which is consistent with the distance to M31, strongly suggesting that this is the correct explanation. Archival XMM/Hubble/Spitzer data show that this AGN can be considered an otherwise unremarkable type-I AGN. Our analysis of the expected rate of background AGN being microlensed by a factor of two or more due to stellar-mass objects in M31 shows that events of this nature should only occur on average every half century. It is thus perhaps surprising that we have uncovered evidence for two more events that are qualitatively similar. A systematic search for new and archival events, with follow-up spectroscopy, is thus warranted. △ Less

Submitted 12 December, 2019; originally announced December 2019.

Comments: 14 pages, 13 figures

The first high-redshift changing-look quasars

Authors: Nicholas P. Ross, Matthew J. Graham, Giorgio Calderone, K. E. Saavik Ford, Barry McKernan, Daniel Stern

Abstract: We report on three redshift $z>2$ quasars with dramatic changes in their C IV emission lines, the first sample of changing-look quasars (CLQs) at high redshift. This is also the first time the changing-look behaviour has been seen in a high-ionisation emission line. SDSS J1205+3422, J1638+2827, and J2228+2201 show interesting behaviour in their observed optical light curves, and subsequent spectro… ▽ More We report on three redshift $z>2$ quasars with dramatic changes in their C IV emission lines, the first sample of changing-look quasars (CLQs) at high redshift. This is also the first time the changing-look behaviour has been seen in a high-ionisation emission line. SDSS J1205+3422, J1638+2827, and J2228+2201 show interesting behaviour in their observed optical light curves, and subsequent spectroscopy shows significant changes in the C IV broad emission line, with both line collapse and emergence being displayed on rest-frame timescales of $\sim$240-1640 days. These are rapid changes, especially when considering virial black hole mass estimates of $M_{\rm BH} > 10^{9} M_{\odot}$ for all three quasars. Continuum and emission line measurements from the three quasars show changes in the continuum-equivalent width plane with the CLQs seen to be on the edge of the full population distribution, and showing indications of an intrinsic Baldwin effect. We put these observations in context with recent state-change models, and note that even in their observed low-state, the C IV CLQs are generally above $\sim$5\% in Eddington luminosity. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 12 pages, 7 figures, 4 tables. All data, analysis code and text are fully available at: github.com/d80b2t/CIV_CLQs. Comments, questions and suggestions welcome and encouraged

Behaviour of the MgII 2798AA Line Over the Full Range of AGN Variability

Authors: David Homan, Chelsea L. Macleod, Andy Lawrence, Nicholas P. Ross, Alastair Bruce

Abstract: We investigate the responsiveness of the 2798AA Mg II broad emission line in AGN on timescales of several years. Our study is based on a sample of extremely variable AGN as well as a broad population sample. The observed response of the line in previous studies has been mixed. By focussing on extreme variability ($|Δg|>$ 1) we find that Mg II clearly does respond to the continuum. However, the deg… ▽ More We investigate the responsiveness of the 2798AA Mg II broad emission line in AGN on timescales of several years. Our study is based on a sample of extremely variable AGN as well as a broad population sample. The observed response of the line in previous studies has been mixed. By focussing on extreme variability ($|Δg|>$ 1) we find that Mg II clearly does respond to the continuum. However, the degree of responsiveness varies strikingly from one object to another; we see cases of Mg II changing by as much as the continuum, more than the continuum, or very little at all. In 74% of the highly variable sample the behaviour of Mg II corresponds with that of H$β$, with 30% of the objects showing large variations in both lines. We do not detect any change in the line width that would correspond to Broad Line Region `breathing', in accordance with results from literature. Some of the objects in our highly variable sample show a clear asymmetry in the Mg II profile. This skewness can be both to the blue and the red of the line centre. Results from our broad population sample show that highly variable quasars have lower Eddington ratios. This result holds for the variability of the continuum, but the correlation is significantly reduced for the variability of the Mg II line. For the first time, we present an overview of the value of the intrinsic Baldwin Effect for Mg II in a large sample. △ Less

Submitted 13 July, 2020; v1 submitted 24 October, 2019; originally announced October 2019.

Comments: 22 pages, 21 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 496, Issue 1, pp.309-327; published 30 May 2020

Stein's method for the Poisson-Dirichlet distribution and the Ewens Sampling Formula, with applications to Wright-Fisher models

Authors: Han L. Gan, Nathan Ross

Abstract: We provide a general theorem bounding the error in the approximation of a random measure of interest--for example, the empirical population measure of types in a Wright-Fisher model--and a Dirichlet process, which is a measure having Poisson-Dirichlet distributed atoms with i.i.d. labels from a diffuse distribution. The implicit metric of the approximation theorem captures the sizes and locations… ▽ More We provide a general theorem bounding the error in the approximation of a random measure of interest--for example, the empirical population measure of types in a Wright-Fisher model--and a Dirichlet process, which is a measure having Poisson-Dirichlet distributed atoms with i.i.d. labels from a diffuse distribution. The implicit metric of the approximation theorem captures the sizes and locations of the masses, and so also yields bounds on the approximation between the masses of the measure of interest and the Poisson-Dirichlet distribution. We apply the result to bound the error in the approximation of the stationary distribution of types in the finite Wright-Fisher model with infinite-alleles mutation structure (not necessarily parent independent) by the Poisson-Dirichlet distribution. An important consequence of our result is an explicit upper bound on the total variation distance between the random partition generated by sampling from a finite Wright-Fisher stationary distribution, and the Ewens Sampling Formula. The bound is small if the sample size $n$ is much smaller than $N^{1/6}\log(N)^{-1/2}$, where $N$ is the total population size. Our analysis requires a result of separate interest, giving an explicit bound on the second moment of the number of types of a finite Wright-Fisher stationary distribution. The general approximation result follows from a new development of Stein's method for the Dirichlet process, which follows by viewing the Dirichlet process as the stationary distribution of a Fleming-Viot process, and then applying Barbour's generator approach. △ Less

Submitted 5 July, 2020; v1 submitted 11 October, 2019; originally announced October 2019.

Comments: Ver2: 40 pages, minor revision; Ver1: 40 pages

Exponential and Laplace approximation for occupation statistics of branching random walk

Authors: Erol Peköz, Adrian Röllin, Nathan Ross

Abstract: We study occupancy counts for the critical nearest-neighbor branching random walk on the $d$-dimensional lattice, conditioned on non-extinction. For $d\geq 3$, Lalley and Zheng (2011) showed that the properly scaled joint distribution of the number of sites occupied by $j$ generation-$n$ particles, $j=1,2,\ldots$, converges in distribution as $n$ goes to infinity, to a deterministic multiple of a… ▽ More We study occupancy counts for the critical nearest-neighbor branching random walk on the $d$-dimensional lattice, conditioned on non-extinction. For $d\geq 3$, Lalley and Zheng (2011) showed that the properly scaled joint distribution of the number of sites occupied by $j$ generation-$n$ particles, $j=1,2,\ldots$, converges in distribution as $n$ goes to infinity, to a deterministic multiple of a single exponential random variable. The limiting exponential variable can be understood as the classical Yaglom limit of the total population size of generation $n$. Here we study the second order fluctuations around this limit, first, by providing a rate of convergence in the Wasserstein metric that holds for all $d\geq3$, and second, by showing that for $d\geq 7$, the weak limit of the scaled joint differences between the number of occupancy-$j$ sites and appropriate multiples of the total population size converge in the Wasserstein metric to a multivariate symmetric Laplace distribution. We also provide a rate of convergence for this latter result. △ Less

Submitted 23 April, 2020; v1 submitted 4 September, 2019; originally announced September 2019.

Comments: Ver2: 22 pages, minor revision; Ver1: 22 pages

Number-Theoretic Characterizations of Some Restricted Clifford+T Circuits

Authors: Matthew Amy, Andrew N. Glaudell, Neil J. Ross

Abstract: Kliuchnikov, Maslov, and Mosca proved in 2012 that a $2\times 2$ unitary matrix $V$ can be exactly represented by a single-qubit Clifford+$T$ circuit if and only if the entries of $V$ belong to the ring $\mathbb{Z}[1/\sqrt{2},i]$. Later that year, Giles and Selinger showed that the same restriction applies to matrices that can be exactly represented by a multi-qubit Clifford+$T$ circuit. These num… ▽ More Kliuchnikov, Maslov, and Mosca proved in 2012 that a $2\times 2$ unitary matrix $V$ can be exactly represented by a single-qubit Clifford+$T$ circuit if and only if the entries of $V$ belong to the ring $\mathbb{Z}[1/\sqrt{2},i]$. Later that year, Giles and Selinger showed that the same restriction applies to matrices that can be exactly represented by a multi-qubit Clifford+$T$ circuit. These number-theoretic characterizations shed new light upon the structure of Clifford+$T$ circuits and led to remarkable developments in the field of quantum compiling. In the present paper, we provide number-theoretic characterizations for certain restricted Clifford+$T$ circuits by considering unitary matrices over subrings of $\mathbb{Z}[1/\sqrt{2},i]$. We focus on the subrings $\mathbb{Z}[1/2]$, $\mathbb{Z}[1/\sqrt{2}]$, $\mathbb{Z}[1/i\sqrt{2}]$, and $\mathbb{Z}[1/2,i]$, and we prove that unitary matrices with entries in these rings correspond to circuits over well-known universal gate sets. In each case, the desired gate set is obtained by extending the set of classical reversible gates $\{X, CX, CCX\}$ with an analogue of the Hadamard gate and an optional phase gate. △ Less

Submitted 3 April, 2020; v1 submitted 16 August, 2019; originally announced August 2019.

Journal ref: Quantum 4, 252 (2020)

Local limit theorems for occupancy models

Authors: A. D. Barbour, Peter Braunsteins, Nathan Ross

Abstract: We present a rather general method for proving local limit theorems, with a good rate of convergence, for sums of dependent random variables. The method is applicable when a Stein coupling can be exhibited. Our approach involves both Stein's method for distributional approximation and Stein's method for concentration. As applications, we prove local central limit theorems with rate of convergence… ▽ More We present a rather general method for proving local limit theorems, with a good rate of convergence, for sums of dependent random variables. The method is applicable when a Stein coupling can be exhibited. Our approach involves both Stein's method for distributional approximation and Stein's method for concentration. As applications, we prove local central limit theorems with rate of convergence for the number of germs with $d$ neighbours in a germ--grain model, and the number of degree-$d$ vertices in an Erdős--Rényi random graph. In both cases, the error rate is optimal, up to logarithmic factors. △ Less

Submitted 5 July, 2020; v1 submitted 1 August, 2019; originally announced August 2019.

Comments: Ver2: 33 pages, minor revision; Ver 1: 32 pages

Host galaxies of high-redshift extremely red and obscured quasars

Authors: Nadia L. Zakamska, Ai-Lei Sun, Michael A. Strauss, Rachael M. Alexandroff, W. N. Brandt, Marco Chiaberge, Jenny E. Greene, Fred Hamann, Guilin Liu, Serena Perrotta, Nicholas P. Ross, Dominika Wylezalek

Abstract: We present Hubble Space Telescope 1.4-1.6 micron images of the hosts of ten extremely red quasars (ERQs) and six type 2 quasar candidates at z=2-3. ERQs, whose bolometric luminosities range between 10^47 and 10^48 erg/sec, show spectroscopic signs of powerful ionized winds, whereas type 2 quasar candidates are less luminous and show only mild outflows. After performing careful subtraction of the q… ▽ More We present Hubble Space Telescope 1.4-1.6 micron images of the hosts of ten extremely red quasars (ERQs) and six type 2 quasar candidates at z=2-3. ERQs, whose bolometric luminosities range between 10^47 and 10^48 erg/sec, show spectroscopic signs of powerful ionized winds, whereas type 2 quasar candidates are less luminous and show only mild outflows. After performing careful subtraction of the quasar light, we clearly detect almost all host galaxies. The median rest-frame B-band luminosity of the ERQ hosts in our sample is 10^11.2 L_Sun, or 4 L* at this redshift. Two of the ten hosts of ERQs are in ongoing mergers. The hosts of the type 2 quasar candidates are 0.6 dex less luminous, with 2/6 in likely ongoing mergers. Intriguingly, despite some signs of interaction and presence of low-mass companions, our objects do not show nearly as much major merger activity as do high-redshift radio-loud galaxies and quasars. In the absence of an overt connection to major ongoing gas-rich merger activity, our observations are consistent with a model in which the near-Eddington accretion and strong feedback of ERQs are associated with relatively late stages of mergers resulting in early-type remnants. These results are in some tension with theoretical expectations of galaxy formation models, in which rapid black hole growth occurs within a short time of a major merger. Type 2 quasar candidates are less luminous, so they may instead be powered by internal galactic processes. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: Accepted to MNRAS; 23 pages, including 9 figures and 2 tables

Ram-pressure stripping of a kicked Hill sphere: Prompt electromagnetic emission from the merger of stellar mass black holes in an AGN accretion disk

Authors: B. McKernan, K. E. S. Ford, I. Bartos, M. J. Graham, W. Lyra, S. Marka, Z. Marka, N. P. Ross, D. Stern, Y. Yang

Abstract: Accretion disks around supermassive black holes (SMBH) are promising sites for stellar mass black hole (BH) mergers due to mass segregation and merger acceleration by disk gas torques. Here we show that a GW-kick at BH merger causes ram-pressure stripping of gas within the BH Hill sphere. If $R_{H}\geq H$, the disk height, an off-center UV flare at $a_{\rm BH} \sim 10^{3}r_{g}$ emerges within… ▽ More Accretion disks around supermassive black holes (SMBH) are promising sites for stellar mass black hole (BH) mergers due to mass segregation and merger acceleration by disk gas torques. Here we show that a GW-kick at BH merger causes ram-pressure stripping of gas within the BH Hill sphere. If $R_{H}\geq H$, the disk height, an off-center UV flare at $a_{\rm BH} \sim 10^{3}r_{g}$ emerges within $t_{\rm UV} \sim \rm{O}(2 {\rm days})(a_{\rm BH}/10^{3}r_{g})(M_{\rm SMBH}/10^{8}M_{\odot})(v_{\rm kick}/10^{2}\rm{km/s})$ post-merger and lasts O$(R_{H}/v_{\rm kick}) \sim \rm{O}(5 t_{\rm UV}$). The flare emerges with luminosity O($10^{42}{\rm erg/s})(t_{\rm UV}/2{\rm days})^{-1}(M_{\rm Hill}/1M_{\odot})(v_{\rm kick}/10^{2}{\rm km/s})^{2}$. AGN optical/UV photometry alters and asymmetric broad emission line profiles can develop after weeks. If $R_{H}<H$, detectability depends on disk optical depth. Follow-up by large optical sky surveys is optimized for small GW error volumes and for LIGO/Virgo triggers $>50M_{\odot}$. △ Less

Submitted 30 September, 2019; v1 submitted 8 July, 2019; originally announced July 2019.

Comments: 5pages, 2 figures, ApJL (accepted)

Journal ref: 2019 ApJ, 884, L50

The Near and Mid-infrared photometric properties of known redshift $z\geq5$ Quasars

Authors: Nicholas P. Ross, Nicholas J. G. Cross

Abstract: We assemble a catalogue of 488 spectroscopically confirmed very high ($z\geq5.00$) redshift quasars and report their near- ($ZYJHK_{s}/K$) and mid- (WISE W1234) infrared properties. 97\% of the VH$z$Q sample is detected in one or more NIR ($ZYJHK/K_{s}$) band, with lack of coverage rather than lack of depth being the reason for the non-detections. 389 (80\%) of the very high redshift quasars are d… ▽ More We assemble a catalogue of 488 spectroscopically confirmed very high ($z\geq5.00$) redshift quasars and report their near- ($ZYJHK_{s}/K$) and mid- (WISE W1234) infrared properties. 97\% of the VH$z$Q sample is detected in one or more NIR ($ZYJHK/K_{s}$) band, with lack of coverage rather than lack of depth being the reason for the non-detections. 389 (80\%) of the very high redshift quasars are detected at 3.4$μ$m in the W1 band from the unWISE catalog and all of the $z\geq7$ quasars are detected in both unWISE W1 and W2. Using archival WFCAM/UKIRT and VIRCAM/VISTA data we check for photometric variability that might be expected from super-Eddington accretion. We find 28 of the quasars have sufficient NIR measurements and signal-to-noise ratio to look for variability. Weak variability was detected in multiple bands of SDSS J0959+0227, and very marginally in the $Y$-band of MMT J0215-0529. Only one quasar, SDSS J0349+0034, shows significant differences between WFCAM and VISTA magnitudes in one band. With supermassive black hole accretion likely to be redshift invariant up to very high-redshift, further monitoring of these sources is warranted. All the data, analysis codes and plots used and generated here can be found at: https://github.com/d80b2t/VHzQ}{\tt github.com/d80b2t/VHzQ. △ Less

Submitted 7 May, 2020; v1 submitted 17 June, 2019; originally announced June 2019.

Comments: Published in MNRAS, 2020, Vol. 494, pp.789-803. 15 pages; 12 figures; 6 tables. Databases available at WFCAM Science Archive and the VISTA Science Archive. All data, analysis codes and plots can be found at github.com/d80b2t/VHzQ. Play the Five Tones

Journal ref: 2020MNRAS.494..789R

Understanding extreme quasar optical variability with CRTS: II. Changing-state quasars

Authors: Matthew J. Graham, Nicholas P. Ross, Daniel Stern, rew J. Drake, Barry McKernan, K. E. Saavik Ford, S. G. Djorgovski, Ashish Mahabal, Eilat Glikman, Steve Larson, Eric Christensen

Abstract: We present the results of a systematic search for quasars in the Catalina Real-time Transient Survey exhibiting both strong photometric and spectroscopic variability over a decadal baseline. We identify 73 sources with specific patterns of optical and mid-IR photometric behavior and a defined spectroscopic change. These "Changing-State" quasars (CSQs) form a higher luminosity sample to complement… ▽ More We present the results of a systematic search for quasars in the Catalina Real-time Transient Survey exhibiting both strong photometric and spectroscopic variability over a decadal baseline. We identify 73 sources with specific patterns of optical and mid-IR photometric behavior and a defined spectroscopic change. These "Changing-State" quasars (CSQs) form a higher luminosity sample to complement existing sets of "Changing-Look" AGN and quasars in the literature. The CSQs (by selection) exhibit larger photometric variability than the CLQs. The spectroscopic variability is marginally stronger in the CSQs than CLQs as defined by the change in H$β$/[OIII] ratio. We find 36 sources with declining H$β$ flux, 37 sources with increasing H$β$ flux and discover seven sources with $z > 0.8$, further extending the redshift arm. Our CSQ sample compares to the literature CLQ objects in similar distributions of H$β$ flux ratios and differential Eddington ratios between high (bright) and low (dim) states. Taken as a whole, we find that this population of extreme varying quasars is associated with changes in the Eddington ratio and the timescales imply cooling/heating fronts propagating through the disk. △ Less

Submitted 6 May, 2019; originally announced May 2019.

Comments: 43 pages, 22 figures, submitted

The tidal disruption event AT2017eqx: spectroscopic evolution from hydrogen rich to poor suggests an atmosphere and outflow

Authors: M. Nicholl, P. K. Blanchard, E. Berger, S. Gomez, R. Margutti, K. D. Alexander, J. Guillochon, J. Leja, R. Chornock, B. Snios, K. Auchettl, A. G. Bruce, P. Challis, D. J. D'Orazio, M. R. Drout, T. Eftekhari, R. J. Foley, O. Graur, C. D. Kilpatrick, A. Lawrence, A. L. Piro, C. Rojas-Bravo, N. P. Ross, P. Short, S. J. Smartt , et al. (2 additional authors not shown)

Abstract: We present and analyse a new tidal disruption event (TDE), AT2017eqx at redshift z=0.1089, discovered by Pan-STARRS and ATLAS. The position of the transient is consistent with the nucleus of its host galaxy; it peaks at a luminosity of $L \approx 10^{44}$ erg s$^{-1}$; and the spectrum shows a persistent blackbody temperature $T \gtrsim 20,000$ K with broad H I and He II emission. The lines are in… ▽ More We present and analyse a new tidal disruption event (TDE), AT2017eqx at redshift z=0.1089, discovered by Pan-STARRS and ATLAS. The position of the transient is consistent with the nucleus of its host galaxy; it peaks at a luminosity of $L \approx 10^{44}$ erg s$^{-1}$; and the spectrum shows a persistent blackbody temperature $T \gtrsim 20,000$ K with broad H I and He II emission. The lines are initially centered at zero velocity, but by 100 days the H I lines disappear while the He II develops a blueshift of $\gtrsim 5,000$ km s$^{-1}$. Both the early- and late-time morphologies have been seen in other TDEs, but the complete transition between them is unprecedented. The evolution can be explained by combining an extended atmosphere, undergoing slow contraction, with a wind in the polar direction becoming visible at late times. Our observations confirm that a lack of hydrogen a TDE spectrum does not indicate a stripped star, while the proposed model implies that much of the diversity in TDEs may be due to the observer viewing angle. Modelling the light curve suggests AT2017eqx resulted from the complete disruption of a solar-mass star by a black hole of $\sim 10^{6.3} M_\odot$. The host is another quiescent, Balmer-strong galaxy, though fainter and less centrally concentrated than most TDE hosts. Radio limits rule out a relativistic jet, while X-ray limits at 500 days are among the deepest for a TDE at this phase. △ Less

Submitted 2 August, 2019; v1 submitted 23 April, 2019; originally announced April 2019.

Comments: Updated to match published version

Journal ref: MNRAS, 488, 1878 (2019)

Opportunities in Time-Domain Extragalactic Astrophysics with the NASA Near-Earth Object Camera (NEOCam)

Authors: Nicholas P. Ross, Roberto J. Assef, Matthew J. Graham, J. Davy Kirkpatrick

Abstract: This White Pape motivates the time domain extragalactic science case for the NASA Near-Earth Object Camera (NEOCam). NEOCam is a NASA Planetary mission whose goal is to discover and characterize asteroids and comets, to assess the hazard to Earth from near-Earth objects, and to study the origin, evolution, and fate of asteroids and comets. NEOCam will, however, cover 68% of the extragalactic sky a… ▽ More This White Pape motivates the time domain extragalactic science case for the NASA Near-Earth Object Camera (NEOCam). NEOCam is a NASA Planetary mission whose goal is to discover and characterize asteroids and comets, to assess the hazard to Earth from near-Earth objects, and to study the origin, evolution, and fate of asteroids and comets. NEOCam will, however, cover 68% of the extragalactic sky and as the NEOWISE-R mission has recently proved, infrared information is now vital for identifying and characterizing the $\gtrsim$10 million IR bright Active Galactic Nuclei, as well as using the IR light curve to provide deep insights into accretion disk astrophysics. NEOWISE-R data has also been used to discover Super-luminous Supernovae, dust echos in Tidal Disruption Events and detects all of the known $z\geq7$ quasars (and over 80% of the known $z\geq6.70$ quasars). As such, for relatively little additional cost, adding the capacity for additional NEOCam data processing (and/or alerting) would have a massive scientific and legacy impact on extragalactic time domain science. △ Less

Submitted 12 April, 2019; originally announced April 2019.

Comments: 8 pages (5 science case); White paper submissionb to the Astro2020 Decadal Survey. All current information at: https://github.com/d80b2t/Astro2020

AGN (and other) astrophysics with Gravitational Wave Events

Authors: K. E. Saavik Ford, Imre Bartos, Barry McKernan, Zoltan Haiman, Alessandra Corsi, Azadeh Keivani, Szabolcs Marka, Rosalba Perna, Matthew Graham, Nicholas P. Ross, Daniel Stern, Jillian Bellovary, Emanuele Berti, Matthew O'Dowd, Wladimir Lyra, Mordecai-Mark Mac Low, Zsuzsanna Marka

Abstract: The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspec… ▽ More The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {\em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.} △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: 5 pages, 2 figures, Astro2020 Decadal Survey white paper

SDSS-IV eBOSS Spectroscopy of X-ray and WISE AGN in Stripe 82X: Overview of the Demographics of X-ray and Mid-Infrared Selected Active Galactic Nuclei

Authors: Stephanie M. LaMassa, Antonis Georgakakis, M. Vivek, Mara Salvato, Tonima Tasnim Ananna, C. Meg Urry, Chelsea MacLeod, Nicholas Ross

Abstract: We report the results of a Sloan Digital Sky Survey-IV eBOSS program to target X-ray sources and mid-infrared-selected WISE AGN candidates in a 36.8 deg$^2$ region of Stripe 82. About half this survey (15.6 deg$^2$) covers the largest contiguous portion of the Stripe 82 X-ray survey. This program represents the largest spectroscopic survey of AGN candidates selected solely by their WISE colors. We… ▽ More We report the results of a Sloan Digital Sky Survey-IV eBOSS program to target X-ray sources and mid-infrared-selected WISE AGN candidates in a 36.8 deg$^2$ region of Stripe 82. About half this survey (15.6 deg$^2$) covers the largest contiguous portion of the Stripe 82 X-ray survey. This program represents the largest spectroscopic survey of AGN candidates selected solely by their WISE colors. We combine this sample with X-ray and WISE AGN in the field identified via other sources of spectroscopy, producing a catalog of 4847 sources that is 82% complete to $r\sim22$. Based on X-ray luminosities or WISE colors, 4730 of these sources are AGN, with a median sample redshift of $z\sim1$. About 30% of the AGN are optically obscured (i.e., lack broad lines in their optical spectra). BPT analysis, however, indicates that 50% of the WISE AGN at $z<0.5$ have emission line ratios consistent with star-forming galaxies, so whether they are buried AGN or star-forming galaxy contaminants is currently unclear. We find that 61% of X-ray AGN are not selected as MIR AGN, with 22% of X-ray AGN undetected by WISE. Most of these latter AGN have high X-ray luminosities ($L_{\rm x} > 10^{44}$ erg s$^{-1}$), indicating that MIR selection misses a sizable fraction of the highest luminosity AGN, as well as lower luminosity sources where AGN heated dust is not dominating the MIR emission. Conversely, $\sim$58% of WISE AGN are undetected by X-rays, though we do not find that they are preferentially redder than the X-ray detected WISE AGN. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: resubmitted to AAS Journals after addressing referee's comments. 21 pages, 17 figures, 5 Tables. Tables 4 & 5 summarize the demographics of the sample. Figures 15 & 16 compare AGN populations from X-ray and MIR selection. The catalog (in fits format) can be downloaded at http://doi.org/10.5281/zenodo.2577357

Spatially Resolved Spectroscopic Study of nearby Seyfert Galaxies: Implications for a Population of "Missed" Seyferts at High-$\textit{z}$

Authors: Junjie Xia, Matthew A. Malkan, Nathaniel R. Ross, Agnes J. Ancheta

Abstract: We present mosaicked long-slit spectral maps of 18 nearby Active Galactic Nuclei (AGNs), 2 LINERs, and 4 star-forming galaxies. With the resulting data cubes taken using the Kast dual spectrograph on the 3 m Shane telescope of the Lick Observatory, we measure the aperture effects on the spectroscopic classification of AGNs. With more starlight included in a larger aperture, the nuclear spectrum th… ▽ More We present mosaicked long-slit spectral maps of 18 nearby Active Galactic Nuclei (AGNs), 2 LINERs, and 4 star-forming galaxies. With the resulting data cubes taken using the Kast dual spectrograph on the 3 m Shane telescope of the Lick Observatory, we measure the aperture effects on the spectroscopic classification of AGNs. With more starlight included in a larger aperture, the nuclear spectrum that is Seyfert-like may become contaminated. We generated standard spectroscopic classification diagrams in different observing apertures. These show quantitatively how the ensemble of Seyferts migrates toward the H $\scriptsize{\textrm{II}}$ region classification when being observed with increasing aperture sizes. But the effect ranges widely in individual active galaxies. Some of the less luminous Seyferts shfit by a large amount, while some other barely move or even shift in different directions. We find that those Seyfert galaxies with the fraction of nuclear H$α$ emission lower than 0.2 of the host galaxy, 2-10 keV hard X-ray luminosity lower than $10^{43}$ erg s$^{-1}$, and the observed nuclear [O $\scriptsize{\textrm{III}}$] luminosity lower than $10^{40.5}$ erg s$^{-1}$, are more likely to change activity classification type when the entire host galaxy is included. Overall, 4 of our 24 galaxies (18 Seyferts) change their spectral activity classification type when observed with a very large aperture. △ Less

Submitted 19 December, 2018; originally announced December 2018.

Comments: 11 pages with 8 figures

Journal ref: The Astrophysical Journal, Volume 869 (2018) 138

Changing-Look Quasar Candidates: First Results from Follow-up Spectroscopy of Highly Optically Variable Quasars

Authors: Chelsea L. MacLeod, Paul J. Green, Scott F. Anderson, Alastair Bruce, Michael Eracleous, Matthew Graham, David Homan, Andy Lawrence, Amy LeBleu, Nicholas P. Ross, John J. Ruan, Jessie Runnoe, Daniel Stern, William Burgett, Kenneth C. Chambers, Nick Kaiser, Eugene Magnier, Nigel Metcalfe

Abstract: Active galactic nuclei (AGN) that show strong rest-frame optical/UV variability in their blue continuum and broad line emission are classified as "changing-look" AGN, or at higher luminosities changing look quasars (CLQs). These surprisingly large and sometimes rapid transitions challenge accepted models of quasar physics and duty cycles, offer several new avenues for study of quasar host galaxies… ▽ More Active galactic nuclei (AGN) that show strong rest-frame optical/UV variability in their blue continuum and broad line emission are classified as "changing-look" AGN, or at higher luminosities changing look quasars (CLQs). These surprisingly large and sometimes rapid transitions challenge accepted models of quasar physics and duty cycles, offer several new avenues for study of quasar host galaxies, and open a wider interpretation of the cause of differences between broad and narrow line AGN. To better characterize extreme quasar variability, we present follow-up spectroscopy as part of a comprehensive search for CLQs across the full SDSS footprint using spectroscopically confirmed quasars from the SDSS DR7 catalog. Our primary selection requires large-amplitude (|Δg|>1 mag, |Δr|>0.5 mag) variability over any of the available time baselines probed by the SDSS and Pan-STARRS 1 surveys. We employ photometry from the Catalina Sky Survey to verify variability behavior in CLQ candidates where available, and confirm CLQs using optical spectroscopy from the William Herschel, MMT, Magellan, and Palomar telescopes. For our adopted S/N threshold on variability of broad Hβemission, we find 17 new CLQs, yielding a confirmation rate of >~ 20%. These candidates are at lower Eddington ratio relative to the overall quasar population which supports a disk-wind model for the broad line region. Based on our sample, the CLQ fraction increases from 10% to roughly half as the continuum flux ratio between repeat spectra at 3420 Angstroms increases from 1.5 to 6. We release a catalog of over 200 highly variable candidates to facilitate future CLQ searches. △ Less

Submitted 6 February, 2019; v1 submitted 28 September, 2018; originally announced October 2018.

Comments: 26 pages, 9 Figures, 3 Tables. Replaced with accepted version

Spectropolarimetry of High Redshift Obscured and Red Quasars

Authors: Rachael M. Alexandroff, Nadia L. Zakamska, Aaron J. Barth, Fred Hamann, Michael A. Strauss, Julian Krolik, Jenny E. Greene, Isabelle Paris, Nicholas P. Ross

Abstract: Spectropolarimetry is a powerful technique that has provided critical support for the geometric unification model of local active galactic nuclei. In this paper, we present optical (rest-frame UV) Keck spectropolarimetry of five luminous obscured (Type 2) and extremely red quasars (ERQs) at z~2.5. Three objects reach polarization fractions of >10% in the continuum. We propose a model in which dust… ▽ More Spectropolarimetry is a powerful technique that has provided critical support for the geometric unification model of local active galactic nuclei. In this paper, we present optical (rest-frame UV) Keck spectropolarimetry of five luminous obscured (Type 2) and extremely red quasars (ERQs) at z~2.5. Three objects reach polarization fractions of >10% in the continuum. We propose a model in which dust scattering is the dominant scattering and polarization mechanism in our targets, though electron scattering cannot be completely excluded. Emission lines are polarized at a lower level than is the continuum. This suggests that the emission-line region exists on similar spatial scales as the scattering region. In three objects we detect an intriguing 90 degree swing in the polarization position angle as a function of line-of-sight velocity in the emission lines of Ly-alpha, CIV and NV. We interpret this phenomenon in the framework of a geometric model with an equatorial dusty scattering region in which the material is outflowing at several thousand km/sec. Emission lines may also be scattered by dust or resonantly. This model explains several salient features of observations by scattering on scales of a few tens of pc. Our observations provide a tantalizing view of the inner region geometry and kinematics of high-redshift obscured and extremely red quasars. Our data and modeling lend strong support for toroidal obscuration and powerful outflows on the scales of the UV emission-line region, in addition to the larger scale outflows inferred previously from the optical emission-line kinematics. △ Less

Submitted 26 June, 2018; originally announced June 2018.

Comments: 26 pages, MNRAS, in press

A new physical interpretation of optical and infrared variability in quasars

Authors: Nicholas P. Ross, K. E. Saavik Ford, Matthew Graham, Barry McKernan, Daniel Stern, Aaron M. Meisner, Roberto J. Assef, Arjun Dey, Andrew J. Drake, Hyunsung D. Jun, Dustin Lang

Abstract: Changing-look quasars are a recently identified class of active galaxies in which the strong UV continuum and/or broad optical hydrogen emission lines associated with unobscured quasars either appear or disappear on timescales of months to years. The physical processes responsible for this behaviour are still debated, but changes in the black hole accretion rate or accretion disk structure appear… ▽ More Changing-look quasars are a recently identified class of active galaxies in which the strong UV continuum and/or broad optical hydrogen emission lines associated with unobscured quasars either appear or disappear on timescales of months to years. The physical processes responsible for this behaviour are still debated, but changes in the black hole accretion rate or accretion disk structure appear more likely than changes in obscuration. Here we report on four epochs of spectroscopy of SDSS J110057.70-005304.5, a quasar at a redshift of $z=0.378$ whose UV continuum and broad hydrogen emission lines have faded, and then returned over the past $\approx$20 years. The change in this quasar was initially identified in the infrared, and an archival spectrum from 2010 shows an intermediate phase of the transition during which the flux below rest-frame $\approx$3400Å has decreased by close to an order of magnitude. This combination is unique compared to previously published examples of changing-look quasars, and is best explained by dramatic changes in the innermost regions of the accretion disk. The optical continuum has been rising since mid-2016, leading to a prediction of a rise in hydrogen emission line flux in the next year. Increases in the infrared flux are beginning to follow, delayed by a $\sim$3 year observed timescale. If our model is confirmed, the physics of changing-look quasars are governed by processes at the innermost stable circular orbit (ISCO) around the black hole, and the structure of the innermost disk. The easily identifiable and monitored changing-look quasars would then provide a new probe and laboratory of the nuclear central engine. △ Less

Submitted 27 September, 2018; v1 submitted 17 May, 2018; originally announced May 2018.

Comments: 13 pages, 4 figures, 3 tables. Published in MNRAS. All code and data links on GitHub, https://github.com/d80b2t/WISE_LC

Journal ref: 2018MNRAS.480.4468R

A Mid-IR Selected Changing-Look Quasar and Physical Scenarios for Abrupt AGN Fading

Authors: Daniel Stern, Barry McKernan, Matthew J. Graham, K. E. S. Ford, Nicholas P. Ross, Aaron M. Meisner, Roberto J. Assef, Mislav Baloković, Murray Brightman, Arjun Dey, Andrew Drake, S. G. Djorgovski, Peter Eisenhardt, Hyunsung D. Jun

Abstract: We report a new changing-look quasar, WISE~J105203.55+151929.5 at $z=0.303$, found by identifying highly mid-IR variable quasars in the WISE/NEOWISE data stream. Compared to multi-epoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at $3.4$ and $4.6 μ$m since 2009. Swift target-of-opportunity ob… ▽ More We report a new changing-look quasar, WISE~J105203.55+151929.5 at $z=0.303$, found by identifying highly mid-IR variable quasars in the WISE/NEOWISE data stream. Compared to multi-epoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at $3.4$ and $4.6 μ$m since 2009. Swift target-of-opportunity observations in 2017 show even stronger fading in the soft X-rays compared to the ROSAT detection of this source in 1995, with at least a factor of fifteen decrease. We obtained second-epoch spectroscopy with the Palomar telescope in 2017 which, when compared with the 2006 archival SDSS spectrum, reveals that the broad H$β$ emission has vanished and that the quasar has become significantly redder. The two most likely interpretations for this dramatic change are source fading or obscuration, where the latter is strongly disfavored by the mid-IR data. We discuss various physical scenarios that could cause such changes in the quasar luminosity over this timescale, and favor changes in the innermost regions of the accretion disk that occur on the thermal and heating/cooling front timescales. We discuss possible physical triggers that could cause these changes, and predict the multiwavelength signatures that could distinguish these physical scenarios. △ Less

Submitted 17 May, 2018; originally announced May 2018.

Comments: 9 pages, 4 figures, submitted to MNRAS