-
Hybrid Oscillator-Qubit Quantum Processors: Simulating Fermions, Bosons, and Gauge Fields
Authors:
Eleanor Crane,
Kevin C. Smith,
Teague Tomesh,
Alec Eickbusch,
John M. Martyn,
Stefan Kühn,
Lena Funcke,
Michael Austin DeMarco,
Isaac L. Chuang,
Nathan Wiebe,
Alexander Schuckert,
Steven M. Girvin
Abstract:
We develop a hybrid oscillator-qubit processor framework for quantum simulation of strongly correlated fermions and bosons that avoids the boson-to-qubit mapping overhead encountered in qubit hardware. This framework gives exact decompositions of particle interactions such as density-density terms and gauge-invariant hopping, as well as approximate methods based on the Baker-Campbell Hausdorff for…
▽ More
We develop a hybrid oscillator-qubit processor framework for quantum simulation of strongly correlated fermions and bosons that avoids the boson-to-qubit mapping overhead encountered in qubit hardware. This framework gives exact decompositions of particle interactions such as density-density terms and gauge-invariant hopping, as well as approximate methods based on the Baker-Campbell Hausdorff formulas including the magnetic field term for the $U(1)$ quantum link model in $(2+1)$D. We use this framework to show how to simulate dynamics using Trotterisation, perform ancilla-free partial error detection using Gauss's law, measure non-local observables, estimate ground state energies using a oscillator-qubit variational quantum eigensolver as well as quantum signal processing, and we numerically study the influence of hardware errors in circuit QED experiments. To show the advantages over all-qubit hardware, we perform an end-to-end comparison of the gate complexity for the gauge-invariant hopping term and find an improvement of the asymptotic scaling with the boson number cutoff $S$ from $\mathcal{O}(\log(S)^2)$ to $\mathcal{O}(1)$ in our framework as well as, for bosonic matter, a constant factor improvement of better than $10^4$. We also find an improvement from $\mathcal{O}(\log(S))$ to $\mathcal{O}(1)$ for the $U(1)$ magnetic field term. While our work focusses on an implementation in superconducting hardware, our framework can also be used in trapped ion, and neutral atom hardware. This work establishes digital quantum simulation with hybrid oscillator-qubit hardware as a viable and advantageous method for the study of qubit-boson models in materials science, chemistry, and high-energy physics.
△ Less
Submitted 5 September, 2024;
originally announced September 2024.
-
Hybrid Oscillator-Qubit Quantum Processors: Instruction Set Architectures, Abstract Machine Models, and Applications
Authors:
Yuan Liu,
Shraddha Singh,
Kevin C. Smith,
Eleanor Crane,
John M. Martyn,
Alec Eickbusch,
Alexander Schuckert,
Richard D. Li,
Jasmine Sinanan-Singh,
Micheline B. Soley,
Takahiro Tsunoda,
Isaac L. Chuang,
Nathan Wiebe,
Steven M. Girvin
Abstract:
Quantum computing with discrete variable (DV, qubit) hardware is approaching the large scales necessary for computations beyond the reach of classical computers. However, important use cases such as quantum simulations of physical models containing bosonic modes, and quantum error correction are challenging for DV-only systems. Separately, hardware containing native continuous-variable (CV, oscill…
▽ More
Quantum computing with discrete variable (DV, qubit) hardware is approaching the large scales necessary for computations beyond the reach of classical computers. However, important use cases such as quantum simulations of physical models containing bosonic modes, and quantum error correction are challenging for DV-only systems. Separately, hardware containing native continuous-variable (CV, oscillator) systems has received attention as an alternative approach, yet the universal control of such systems is non-trivial. In this work, we show that hybrid CV-DV hardware offers a great advantage in meeting these challenges, offering a powerful computational paradigm that inherits the strengths of both DV and CV processors. We provide a pedagogical introduction to CV-DV systems and the multiple abstraction layers needed to produce a full software stack connecting applications to hardware. We present a variety of new hybrid CV-DV compilation techniques, algorithms, and applications, including the extension of quantum signal processing concepts to CV-DV systems and strategies to simulate systems of interacting spins, fermions, and bosons. To facilitate the development of hybrid CV-DV processor systems, we introduce formal Abstract Machine Models and Instruction Set Architectures -- essential abstractions that enable developers to formulate applications, compile algorithms, and explore the potential of current and future hardware for realizing fault-tolerant circuits, modules, and processors. Hybrid CV-DV quantum computations are beginning to be performed in superconducting, trapped ion, and neutral atom platforms, and large-scale experiments are set to be demonstrated in the near future. We present a timely and comprehensive guide to this relatively unexplored yet promising approach to quantum computation and providing an architectural backbone to guide future development.
△ Less
Submitted 5 August, 2024; v1 submitted 14 July, 2024;
originally announced July 2024.
-
Analog Quantum Simulator of a Quantum Field Theory with Fermion-Spin Systems in Silicon
Authors:
Ali Rad,
Alexander Schuckert,
Eleanor Crane,
Gautam Nambiar,
Fan Fei,
Jonathan Wyrick,
Richard M. Silver,
Mohammad Hafezi,
Zohreh Davoudi,
Michael J. Gullans
Abstract:
Simulating fermions coupled to spin degrees of freedom, relevant for a range of quantum field theories, represents a promising application for quantum simulators. Mapping fermions to qubits is challenging in $2+1$ and higher spacetime dimensions, and mapping bosons demands substantial quantum-computational overhead. These features complicate the realization of mixed fermion-boson quantum systems i…
▽ More
Simulating fermions coupled to spin degrees of freedom, relevant for a range of quantum field theories, represents a promising application for quantum simulators. Mapping fermions to qubits is challenging in $2+1$ and higher spacetime dimensions, and mapping bosons demands substantial quantum-computational overhead. These features complicate the realization of mixed fermion-boson quantum systems in digital quantum computers. We propose a native fermion-(large-)spin analog quantum simulator by utilizing dopant arrays in silicon. Specifically, we show how to use a dynamical lattice of coupled nuclear spins and conduction-band electrons to realize a quantum field theory: an extended Jackiw-Rebbi model involving coupled fermions and quantum rotors. We demonstrate the feasibility of observing dynamical mass generation and a confinement-deconfinement quantum phase transition in 1+1 dimensions on this platform, even in the presence of strong long-range Coulomb interactions. Furthermore, we employ finite-temperature Hartree-Fock-Bogoliubov simulations to investigate the dynamics of mass generation in two-dimensional square and honeycomb arrays, showing that this phenomenon can be simulated with realistic experimental parameters. Our findings reveal two distinct phases, and demonstrate robustness against the addition of Coulomb interactions. Finally, we discuss experimental signatures of the phases through transport and local charge sensing in dopant arrays. This study lays the foundation for quantum simulations of quantum field theories exhibiting fermions coupled to spin degrees of freedom using donors in silicon.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
Gaussian consensus processes and their Lyapunov exponents
Authors:
Edward Crane,
Stanislav Volkov
Abstract:
We introduce a simple dynamic model of opinion formation, in which a finite population of individuals hold vector-valued opinions. At each time step, each individual's opinion moves towards the mean opinion but is then perturbed independently by a centred multivariate Gaussian random variable, with covariance proportional to the covariance matrix of the opinions of the population. We establish pre…
▽ More
We introduce a simple dynamic model of opinion formation, in which a finite population of individuals hold vector-valued opinions. At each time step, each individual's opinion moves towards the mean opinion but is then perturbed independently by a centred multivariate Gaussian random variable, with covariance proportional to the covariance matrix of the opinions of the population. We establish precise necessary and sufficient conditions on the parameters of the model, under which all opinions converge to a common limiting value. Asymptotically perfect correlation emerges between opinions on different topics. Our results are rigorous and based on properties of the partial products of an i.i.d. sequence of random matrices. Each matrix is a fixed linear combination of the identity matrix and a real Ginibre matrix. We derive an analytic expression for the maximal Lyapunov exponent of this product sequence. We also analyze a continuous-time analogue of our model.
△ Less
Submitted 8 May, 2024;
originally announced May 2024.
-
Observation of a finite-energy phase transition in a one-dimensional quantum simulator
Authors:
Alexander Schuckert,
Or Katz,
Lei Feng,
Eleanor Crane,
Arinjoy De,
Mohammad Hafezi,
Alexey V. Gorshkov,
Christopher Monroe
Abstract:
One of the most striking many-body phenomena in nature is the sudden change of macroscopic properties as the temperature or energy reaches a critical value. Such equilibrium transitions have been predicted and observed in two and three spatial dimensions, but have long been thought not to exist in one-dimensional (1D) systems. Fifty years ago, Dyson and Thouless pointed out that a phase transition…
▽ More
One of the most striking many-body phenomena in nature is the sudden change of macroscopic properties as the temperature or energy reaches a critical value. Such equilibrium transitions have been predicted and observed in two and three spatial dimensions, but have long been thought not to exist in one-dimensional (1D) systems. Fifty years ago, Dyson and Thouless pointed out that a phase transition in 1D can occur in the presence of long-range interactions, but an experimental realization has so far not been achieved due to the requirement to both prepare equilibrium states and realize sufficiently long-range interactions. Here we report on the first experimental demonstration of a finite-energy phase transition in 1D. We use the simple observation that finite-energy states can be prepared by time-evolving product initial states and letting them thermalize under the dynamics of a many-body Hamiltonian. By preparing initial states with different energies in a 1D trapped-ion quantum simulator, we study the finite-energy phase diagram of a long-range interacting quantum system. We observe a ferromagnetic equilibrium phase transition as well as a crossover from a low-energy polarized paramagnet to a high-energy unpolarized paramagnet in a system of up to $23$ spins, in excellent agreement with numerical simulations. Our work demonstrates the ability of quantum simulators to realize and study previously inaccessible phases at finite energy density.
△ Less
Submitted 30 October, 2023;
originally announced October 2023.
-
Momentum-space imaging of ultra-thin electron liquids in delta-doped silicon
Authors:
Procopios Constantinou,
Taylor J. Z. Stock,
Eleanor Crane,
Alexander Kölker,
Marcel van Loon,
Juerong Li,
Sarah Fearn,
Henric Bornemann,
Nicolò D'Anna,
Andrew J. Fisher,
Vladimir N. Strocov,
Gabriel Aeppli,
Neil J. Curson,
Steven R. Schofield
Abstract:
Two-dimensional dopant layers ($δ$-layers) in semiconductors provide the high-mobility electron liquids (2DELs) needed for nanoscale quantum-electronic devices. Key parameters such as carrier densities, effective masses, and confinement thicknesses for 2DELs have traditionally been extracted from quantum magnetotransport. In principle, the parameters are immediately readable from the one-electron…
▽ More
Two-dimensional dopant layers ($δ$-layers) in semiconductors provide the high-mobility electron liquids (2DELs) needed for nanoscale quantum-electronic devices. Key parameters such as carrier densities, effective masses, and confinement thicknesses for 2DELs have traditionally been extracted from quantum magnetotransport. In principle, the parameters are immediately readable from the one-electron spectral function that can be measured by angle-resolved photoemission spectroscopy (ARPES). Here, buried 2DEL $δ$-layers in silicon are measured with soft X-ray (SX) ARPES to obtain detailed information about their filled conduction bands and extract device-relevant properties. This study takes advantage of the larger probing depth and photon energy range of SX-ARPES relative to vacuum ultraviolet (VUV) ARPES to accurately measure the $δ$-layer electronic confinement. The measurements are made on ambient-exposed samples and yield extremely thin ($\approx 1$ $nm$) and dense ($\approx$ $10^{14}$ $cm^2$) 2DELs. Critically, this method is used to show that $δ$-layers of arsenic exhibit better electronic confinement than $δ$-layers of phosphorus fabricated under identical conditions.
△ Less
Submitted 29 September, 2023;
originally announced September 2023.
-
Measuring the Loschmidt amplitude for finite-energy properties of the Fermi-Hubbard model on an ion-trap quantum computer
Authors:
Kévin Hémery,
Khaldoon Ghanem,
Eleanor Crane,
Sara L. Campbell,
Joan M. Dreiling,
Caroline Figgatt,
Cameron Foltz,
John P. Gaebler,
Jacob Johansen,
Michael Mills,
Steven A. Moses,
Juan M. Pino,
Anthony Ransford,
Mary Rowe,
Peter Siegfried,
Russell P. Stutz,
Henrik Dreyer,
Alexander Schuckert,
Ramil Nigmatullin
Abstract:
Calculating the equilibrium properties of condensed matter systems is one of the promising applications of near-term quantum computing. Recently, hybrid quantum-classical time-series algorithms have been proposed to efficiently extract these properties from a measurement of the Loschmidt amplitude $\langle ψ| e^{-i \hat H t}|ψ\rangle$ from initial states $|ψ\rangle$ and a time evolution under the…
▽ More
Calculating the equilibrium properties of condensed matter systems is one of the promising applications of near-term quantum computing. Recently, hybrid quantum-classical time-series algorithms have been proposed to efficiently extract these properties from a measurement of the Loschmidt amplitude $\langle ψ| e^{-i \hat H t}|ψ\rangle$ from initial states $|ψ\rangle$ and a time evolution under the Hamiltonian $\hat H$ up to short times $t$. In this work, we study the operation of this algorithm on a present-day quantum computer. Specifically, we measure the Loschmidt amplitude for the Fermi-Hubbard model on a $16$-site ladder geometry (32 orbitals) on the Quantinuum H2-1 trapped-ion device. We assess the effect of noise on the Loschmidt amplitude and implement algorithm-specific error mitigation techniques. By using a thus-motivated error model, we numerically analyze the influence of noise on the full operation of the quantum-classical algorithm by measuring expectation values of local observables at finite energies. Finally, we estimate the resources needed for scaling up the algorithm.
△ Less
Submitted 22 September, 2023; v1 submitted 19 September, 2023;
originally announced September 2023.
-
Leveraging Hamiltonian Simulation Techniques to Compile Operations on Bosonic Devices
Authors:
Christopher Kang,
Micheline B. Soley,
Eleanor Crane,
S. M. Girvin,
Nathan Wiebe
Abstract:
Circuit QED enables the combined use of qubits and oscillator modes. Despite a variety of available gate sets, many hybrid qubit-boson (i.e., oscillator) operations are realizable only through optimal control theory (OCT) which is oftentimes intractable and uninterpretable. We introduce an analytic approach with rigorously proven error bounds for realizing specific classes of operations via two ma…
▽ More
Circuit QED enables the combined use of qubits and oscillator modes. Despite a variety of available gate sets, many hybrid qubit-boson (i.e., oscillator) operations are realizable only through optimal control theory (OCT) which is oftentimes intractable and uninterpretable. We introduce an analytic approach with rigorously proven error bounds for realizing specific classes of operations via two matrix product formulas commonly used in Hamiltonian simulation, the Lie--Trotter and Baker--Campbell--Hausdorff product formulas. We show how this technique can be used to realize a number of operations of interest, including polynomials of annihilation and creation operators, i.e., $a^p {a^\dagger}^q$ for integer $p, q$. We show examples of this paradigm including: obtaining universal control within a subspace of the entire Fock space of an oscillator, state preparation of a fixed photon number in the cavity, simulation of the Jaynes--Cummings Hamiltonian, simulation of the Hong-Ou-Mandel effect and more. This work demonstrates how techniques from Hamiltonian simulation can be applied to better control hybrid boson-qubit devices.
△ Less
Submitted 10 August, 2024; v1 submitted 27 March, 2023;
originally announced March 2023.
-
The limit point in the Jante's law process has an absolutely continuous distribution
Authors:
Edward Crane,
Stanislav Volkov
Abstract:
We study a stochastic model of consensus formation, introduced in 2015 by Grinfeld, Volkov and Wade, who called it a multidimensional randomized Keynesian beauty contest. The model was generalized by Kennerberg and Volkov, who called their generalization the Jante's law process. We consider a version of the model where the space of possible opinions is a convex body $\mathcal{B}$ in…
▽ More
We study a stochastic model of consensus formation, introduced in 2015 by Grinfeld, Volkov and Wade, who called it a multidimensional randomized Keynesian beauty contest. The model was generalized by Kennerberg and Volkov, who called their generalization the Jante's law process. We consider a version of the model where the space of possible opinions is a convex body $\mathcal{B}$ in $\mathbb{R}^d$. $N$ individuals in a population each hold a (multidimensional) opinion in $\mathcal{B}$. Repeatedly, the individual whose opinion is furthest from the center of mass of the $N$ current opinions chooses a new opinion, sampled uniformly at random from $\mathcal{B}$. Kennerberg and Volkov showed that the set of opinions that are not furthest from the center of mass converges to a random limit point. We show that the distribution of the limit opinion is continuous, thus proving the conjecture made after Proposition 3.2 in Grinfeld et al.
△ Less
Submitted 8 November, 2023; v1 submitted 28 November, 2022;
originally announced November 2022.
-
Deterministic constant-depth preparation of the AKLT state on a quantum processor using fusion measurements
Authors:
Kevin C. Smith,
Eleanor Crane,
Nathan Wiebe,
S. M. Girvin
Abstract:
The ground state of the spin-1 Affleck, Kennedy, Lieb and Tasaki (AKLT) model is a paradigmatic example of both a matrix product state and a symmetry-protected topological phase, and additionally holds promise as a resource state for measurement-based quantum computation. Having a nonzero correlation length, the AKLT state cannot be exactly prepared by a constant-depth unitary circuit composed of…
▽ More
The ground state of the spin-1 Affleck, Kennedy, Lieb and Tasaki (AKLT) model is a paradigmatic example of both a matrix product state and a symmetry-protected topological phase, and additionally holds promise as a resource state for measurement-based quantum computation. Having a nonzero correlation length, the AKLT state cannot be exactly prepared by a constant-depth unitary circuit composed of local gates. In this work, we demonstrate that this no-go limit can be evaded by augmenting a constant-depth circuit with fusion measurements, such that the total preparation time is independent of system size and entirely deterministic. We elucidate our preparation scheme using the language of tensor networks, and furthermore show that the $\mathbb{Z}_2\times\mathbb{Z}_2$ symmetry of the AKLT state directly affords this speed-up over previously known preparation methods. To demonstrate the practical advantage of measurement-assisted preparation on noisy intermediate-scale quantum (NISQ) devices, we carry out our protocol on an IBM Quantum processor. We measure both the string order and entanglement spectrum of prepared AKLT chains and, employing these as metrics, find improved results over the known (purely unitary) sequential preparation approach. We conclude with a demonstration of quantum teleportation using the AKLT state prepared by our measurement-assisted scheme. This work thus serves to provide an efficient strategy to prepare a specific resource in the form of the AKLT state and, more broadly, experimentally demonstrates the possibility for realizable improvement in state preparation afforded by measurement-based circuit depth reduction strategies on NISQ-era devices.
△ Less
Submitted 10 April, 2023; v1 submitted 31 October, 2022;
originally announced October 2022.
-
Bosonic Qiskit
Authors:
Timothy J Stavenger,
Eleanor Crane,
Kevin Smith,
Christopher T Kang,
Steven M Girvin,
Nathan Wiebe
Abstract:
The practical benefits of hybrid quantum information processing hardware that contains continuous-variable objects (bosonic modes such as mechanical or electromagnetic oscillators) in addition to traditional (discrete-variable) qubits have recently been demonstrated by experiments with bosonic codes that reach the break-even point for quantum error correction and by efficient Gaussian boson sampli…
▽ More
The practical benefits of hybrid quantum information processing hardware that contains continuous-variable objects (bosonic modes such as mechanical or electromagnetic oscillators) in addition to traditional (discrete-variable) qubits have recently been demonstrated by experiments with bosonic codes that reach the break-even point for quantum error correction and by efficient Gaussian boson sampling simulation of the Franck-Condon spectra of triatomic molecules that is well beyond the capabilities of current qubit-only hardware. The goal of this Co-design Center for Quantum Advantage (C2QA) project is to develop an instruction set architecture (ISA) for hybrid qubit/bosonic mode systems that contains an inventory of the fundamental operations and measurements that are possible in such hardware. The corresponding abstract machine model (AMM) would also contain a description of the appropriate error models associated with the gates, measurements and time evolution of the hardware. This information has been implemented as an extension of Qiskit. Qiskit is an opensource software development toolkit (SDK) for simulating the quantum state of a quantum circuit on a system with Python 3.7+ and for running the same circuits on prototype hardware within the IBM Quantum Lab. We introduce the Bosonic Qiskit software to enable the simulation of hybrid qubit/bosonic systems using the existing Qiskit software development kit. This implementation can be used for simulating new hybrid systems, verifying proposed physical systems, and modeling systems larger than can currently be constructed. We also cover tutorials and example use cases included within the software to study Jaynes- Cummings models, bosonic Hubbard models, plotting Wigner functions and animations, and calculating maximum likelihood estimations using Wigner functions.
△ Less
Submitted 2 December, 2022; v1 submitted 22 September, 2022;
originally announced September 2022.
-
Probing finite-temperature observables in quantum simulators of spin systems with short-time dynamics
Authors:
Alexander Schuckert,
Annabelle Bohrdt,
Eleanor Crane,
Michael Knap
Abstract:
Preparing finite temperature states in quantum simulators of spin systems, such as trapped ions or Rydberg atoms in optical tweezers, is challenging due to their almost perfect isolation from the environment. Here, we show how finite-temperature observables can be obtained with an algorithm motivated from the Jarzynski equality and equivalent to the one in Lu, Banuls and Cirac, PRX Quantum 2, 0203…
▽ More
Preparing finite temperature states in quantum simulators of spin systems, such as trapped ions or Rydberg atoms in optical tweezers, is challenging due to their almost perfect isolation from the environment. Here, we show how finite-temperature observables can be obtained with an algorithm motivated from the Jarzynski equality and equivalent to the one in Lu, Banuls and Cirac, PRX Quantum 2, 020321 (2021). It consists of classical importance sampling of initial states and a measurement of the Loschmidt echo with a quantum simulator. We use the method as a quantum-inspired classical algorithm and simulate the protocol with matrix product states to analyze the requirements on a quantum simulator. This way, we show that a finite temperature phase transition in the long-range transverse field Ising model can be characterized in trapped ion quantum simulators. We propose a concrete measurement protocol for the Loschmidt echo and discuss the influence of measurement noise, dephasing, as well as state preparation and measurement errors. We argue that the algorithm is robust against those imperfections under realistic conditions.
△ Less
Submitted 4 May, 2023; v1 submitted 3 June, 2022;
originally announced June 2022.
-
Visualizing spinon Fermi surfaces with time-dependent spectroscopy
Authors:
Alexander Schuckert,
Annabelle Bohrdt,
Eleanor Crane,
Fabian Grusdt
Abstract:
Quantum simulation experiments have started to explore regimes that are not accessible with exact numerical methods. In order to probe these systems and enable new physical insights, the need for measurement protocols arises that can bridge the gap to solid state experiments, and at the same time make optimal use of the capabilities of quantum simulation experiments. Here we propose applying time-…
▽ More
Quantum simulation experiments have started to explore regimes that are not accessible with exact numerical methods. In order to probe these systems and enable new physical insights, the need for measurement protocols arises that can bridge the gap to solid state experiments, and at the same time make optimal use of the capabilities of quantum simulation experiments. Here we propose applying time-dependent photo-emission spectroscopy, an established tool in solid state systems, in cold atom quantum simulators. Concretely, we suggest combining the method with large magnetic field gradients, unattainable in experiments on real materials, to drive Bloch oscillations of spinons, the emergent quasiparticles of spin liquids. We show in exact diagonalization simulations of the one-dimensional $t-J$ model that the spinons start to populate previously unoccupied states in an effective band structure, thus allowing to visualize states invisible in the equilibrium spectrum. The dependence of the spectral function on the time after the pump pulse reveals collective interactions among spinons. In numerical simulations of small two-dimensional systems, spectral weight appears at the ground state energy at momentum $\mathbf{q} = (π,π)$, where the equilibrium spectral response is strongly suppressed up to higher energies, indicating a possible route towards solving the mystery of the Fermi arcs in the cuprate materials.
△ Less
Submitted 27 May, 2021;
originally announced May 2021.
-
Rydberg Entangling Gates in Silicon
Authors:
Eleanor Crane,
Alexander Schuckert,
Nguyen H. Le,
Andrew J. Fisher
Abstract:
In this paper, we propose a new Rydberg entangling gate scheme which we demonstrate theoretically to have an order of magnitude improvement in fidelities and speed over existing protocols. We find that applying this gate to donors in silicon would help overcome the strenuous requirements on atomic precision donor placement and substantial gate tuning, which so far has hampered scaling. We calculat…
▽ More
In this paper, we propose a new Rydberg entangling gate scheme which we demonstrate theoretically to have an order of magnitude improvement in fidelities and speed over existing protocols. We find that applying this gate to donors in silicon would help overcome the strenuous requirements on atomic precision donor placement and substantial gate tuning, which so far has hampered scaling. We calculate multivalley Rydberg interactions for several donor species using the Finite Element Method, and show that induced electric dipole and Van der Waals interactions, calculated here for the first time, are important even for low-lying excited states. We show that Rydberg gate operation is possible within the lifetime of donor excited states with 99.9% fidelity for the creation of a Bell state in the presence of decoherence.
△ Less
Submitted 24 November, 2020; v1 submitted 26 August, 2020;
originally announced August 2020.
-
Well-posedness of the mean field forest fire age evolution equation
Authors:
Edward Crane
Abstract:
We prove the well-posedness of a differential equation that describes the evolution of the large-system limit of the empirical age measure in the mean field forest fire model of Ráth and Tóth (arXiv:0808.2116). This forest fire model is a random graph process on $n$ vertices, whose dynamics combine the Erdős-Rényi dynamics with a Poisson rain of lightning strikes. All edges in any connected compon…
▽ More
We prove the well-posedness of a differential equation that describes the evolution of the large-system limit of the empirical age measure in the mean field forest fire model of Ráth and Tóth (arXiv:0808.2116). This forest fire model is a random graph process on $n$ vertices, whose dynamics combine the Erdős-Rényi dynamics with a Poisson rain of lightning strikes. All edges in any connected component are deleted as soon as any of its vertices is struck by lightning. Each vertex has an age, which increases at rate $1$ but is reset to $0$ each time it burns. We consider the asymptotic lightning regime in which the model displays self-organized criticality. Crane, Ráth and Yeo (arXiv:1811.07981) take the initial state to be an inhomogeneous random graph whose edge probabilities depend on the ages of the vertices. They show that as $n \to \infty$ the empirical age distribution converges as a process to the solution of a deterministic autonomous differential equation. It is a nonlinear age-dependent population dynamics model whose age-specific mortality modulus involves the leading eigenfunction of the branching operator of an associated multitype branching process. The differential equation displays self-organized criticality in the sense that the leading eigenvalue of the branching operator is held at $1$ without this being imposed as a boundary condition.
△ Less
Submitted 11 July, 2020;
originally announced July 2020.
-
Atomic-Scale Patterning of Arsenic in Silicon by Scanning Tunneling Microscopy
Authors:
Taylor J. Z. Stock,
Oliver Warschkow,
Procopios C. Constantinou,
Juerong Li,
Sarah Fearn,
Eleanor Crane,
Emily V. S. Hofmann,
Alexander Kölker,
David R. McKenzie,
Steven R. Schofield,
Neil J. Curson
Abstract:
Over the last two decades, prototype devices for future classical and quantum computing technologies have been fabricated, by using scanning tunneling microscopy and hydrogen resist lithography to position phosphorus atoms in silicon with atomic-scale precision. Despite these successes, phosphine remains the only donor precursor molecule to have been demonstrated as compatible with the hydrogen re…
▽ More
Over the last two decades, prototype devices for future classical and quantum computing technologies have been fabricated, by using scanning tunneling microscopy and hydrogen resist lithography to position phosphorus atoms in silicon with atomic-scale precision. Despite these successes, phosphine remains the only donor precursor molecule to have been demonstrated as compatible with the hydrogen resist lithography technique. The potential benefits of atomic-scale placement of alternative dopant species have, until now, remained unexplored. In this work, we demonstrate successful fabrication of atomic-scale structures of arsenic-in-silicon. Using a scanning tunneling microscope tip, we pattern a monolayer hydrogen mask to selectively place arsenic atoms on the Si(001) surface using arsine as the precursor molecule. We fully elucidate the surface chemistry and reaction pathways of arsine on Si(001), revealing significant differences to phosphine. We explain how these differences result in enhanced surface immobilization and in-plane confinement of arsenic compared to phosphorus, and a dose-rate independent arsenic saturation density of $0.24{\pm}0.04$ ML. We demonstrate the successful encapsulation of arsenic delta-layers using silicon molecular beam epitaxy, and find electrical characteristics that are competitive with equivalent structures fabricated with phosphorus. Arsenic delta-layers are also found to offer improvement in out-of-plane confinement compared to similarly prepared phosphorus layers, while still retaining >80% carrier activation and sheet resistances of $<2 kΩ/{\square}$. These excellent characteristics of arsenic represent opportunities to enhance existing capabilities of atomic-scale fabrication of dopant structures in silicon, and are particularly important for three-dimensional devices, where vertical control of the position of device components is critical.
△ Less
Submitted 15 October, 2019;
originally announced October 2019.
-
Optically Controlled Entangling Gates in Randomly Doped Silicon
Authors:
Eleanor Crane,
Thomas Crane,
Nguyen H. Le,
Alexander Schuckert,
Andrew J. Fisher
Abstract:
Randomly-doped silicon has many competitive advantages for quantum computation; not only is it fast to fabricate but it could naturally contain high numbers of qubits and logic gates as a function of doping densities. We determine the densities of entangling gates in randomly doped silicon comprising two different dopant species. First, we define conditions and plot maps of the relative locations…
▽ More
Randomly-doped silicon has many competitive advantages for quantum computation; not only is it fast to fabricate but it could naturally contain high numbers of qubits and logic gates as a function of doping densities. We determine the densities of entangling gates in randomly doped silicon comprising two different dopant species. First, we define conditions and plot maps of the relative locations of the dopants necessary for them to form exchange interaction mediated entangling gates. Second, using nearest neighbour Poisson point process theory, we calculate the doping densities necessary for maximal densities of single and dual-species gates. We find agreement of our results with a Monte Carlo simulation, for which we present the algorithms, which handles multiple donor structures and scales optimally with the number of dopants and use it to extract donor structures not captured by our Poisson point process theory. Third, using the moving average cluster expansion technique, we make predictions for a proof of principle experiment demonstrating the control of one species by the orbital excitation of another. These combined approaches to density optimization in random distributions may be useful for other condensed matter systems as well as applications outside physics.
△ Less
Submitted 1 February, 2019;
originally announced February 2019.
-
Age evolution in the mean field forest fire model via multitype branching processes
Authors:
Edward Crane,
Balazs Rath,
Dominic Yeo
Abstract:
We study the distribution of ages in the mean field forest fire model introduced by Ráth and Tóth. This model is an evolving random graph whose dynamics combine Erdős-Rényi edge-addition with a Poisson rain of lightning strikes. All edges in a connected component are deleted when any of its vertices is struck by lightning. We consider the asymptotic regime of lightning rates for which the model di…
▽ More
We study the distribution of ages in the mean field forest fire model introduced by Ráth and Tóth. This model is an evolving random graph whose dynamics combine Erdős-Rényi edge-addition with a Poisson rain of lightning strikes. All edges in a connected component are deleted when any of its vertices is struck by lightning. We consider the asymptotic regime of lightning rates for which the model displays self-organized criticality. The age of a vertex increases at unit rate, but it is reset to zero at each burning time. We show that the empirical age distribution converges as a process to a deterministic solution of an autonomous measure-valued differential equation. The main technique is to observe that, conditioned on the vertex ages, the graph is an inhomogeneous random graph in the sense of Bollobás, Janson and Riordan. We then study the evolution of the ages via the multitype Galton--Watson trees that arise as the limit in law of the component of an identified vertex at any fixed time. These trees are critical from the gelation time onwards.
△ Less
Submitted 9 December, 2020; v1 submitted 19 November, 2018;
originally announced November 2018.
-
Steady state clusters and the Rath-Toth mean field forest fire model
Authors:
Edward Crane
Abstract:
We introduce a random finite rooted tree $\mathcal{C}$, the steady state cluster, characterized by a recursive description: $\mathcal{C}$ is a singleton with probability $1/2$ and otherwise is obtained by joining by an edge the roots of two independent trees $\mathcal{C}'$ and $\mathcal{C}''$, each having the law of $\mathcal{C}$, then re-rooting the resulting tree at a uniform random vertex.
We…
▽ More
We introduce a random finite rooted tree $\mathcal{C}$, the steady state cluster, characterized by a recursive description: $\mathcal{C}$ is a singleton with probability $1/2$ and otherwise is obtained by joining by an edge the roots of two independent trees $\mathcal{C}'$ and $\mathcal{C}''$, each having the law of $\mathcal{C}$, then re-rooting the resulting tree at a uniform random vertex.
We construct a stationary regenerative stochastic process $\mathcal{C}(t)$, the steady state cluster growth process. It is characterized by a simple fixed-point property. Its stationary distribution is the law of the steady state cluster $\mathcal{C}$. We conjecture that $\mathcal{C}(t)$ is the local limit of the evolution of the cluster of a tagged vertex in the stationary state of the mean field forest fire model of Ráth and Tóth. We describe its explosions in terms of a Lévy subordinator, using a state-dependent time change.
The steady state cluster is also a multitype Galton-Watson tree with a continuum of types. The steady state cluster conditioned on its size is a random weighted spanning tree of the complete graph equipped random edge weights with a simple explicit joint distribution. The time-reversal of the steady state cluster growth process is realised as the component of a `uniform' vertex in a logging process of a critical multitype Galton-Watson tree conditioned to be infinite.
We construct a stationary forest fire model on the infinite rooted tree $\mathbb{Z}^*$ with the property that the evolution of the cluster of the root is a version of the steady state cluster growth process. This model is similar in spirit to Aldous' frozen percolation model on the rooted infinite binary tree. We conjecture that it is the local weak limit of the stationary Ráth-Tóth model.
△ Less
Submitted 10 September, 2018;
originally announced September 2018.
-
Diffusion and superdiffusion in lattice models of colliding particles with stored momentum
Authors:
Edward Crane,
Sean Ledger,
Balint Toth
Abstract:
We introduce two discrete models of a collection of colliding particles with stored momentum and study the asymptotic growth of the mean-square displacement of an active particle. We prove that the models are superdiffusive in one dimension (with power law correction) and diffusive in three and higher dimensions. In two dimensions we demonstrate superdiffusivity (with logarithmic correction) for c…
▽ More
We introduce two discrete models of a collection of colliding particles with stored momentum and study the asymptotic growth of the mean-square displacement of an active particle. We prove that the models are superdiffusive in one dimension (with power law correction) and diffusive in three and higher dimensions. In two dimensions we demonstrate superdiffusivity (with logarithmic correction) for certain anisotropic initial conditions.
△ Less
Submitted 10 September, 2018;
originally announced September 2018.
-
Functional Large Deviations for Cox Processes and $Cox/G/\infty$ Queues, with a Biological Application
Authors:
Justin Dean,
Ayalvadi Ganesh,
Edward Crane
Abstract:
We consider an infinite-server queue into which customers arrive according to a Cox process and have independent service times with a general distribution. We prove a functional large deviations principle for the equilibrium queue length process. The model is motivated by a linear feed-forward gene regulatory network, in which the rate of protein synthesis is modulated by the number of RNA molecul…
▽ More
We consider an infinite-server queue into which customers arrive according to a Cox process and have independent service times with a general distribution. We prove a functional large deviations principle for the equilibrium queue length process. The model is motivated by a linear feed-forward gene regulatory network, in which the rate of protein synthesis is modulated by the number of RNA molecules present in a cell. The system can be modelled as a tandem of infinite-server queues, in which the number of customers present in a queue modulates the arrival rate into the next queue in the tandem. We establish large deviation principles for this queueing system in the asymptotic regime in which the arrival process is sped up, while the service process is not scaled.
△ Less
Submitted 27 March, 2020; v1 submitted 13 August, 2018;
originally announced August 2018.
-
Circle Packings with Generalized Branching
Authors:
James Ashe,
Edward Crane,
Kenneth Stephenson
Abstract:
Attempts to build a discrete theory for rational maps on the sphere via circle packing have foundered on discretization effects in locating branch points. The authors remove this impediment by introducing generalized branch points. A generalized branch point need no longer be attached to an individual circle, but with the help of chaperones and other devices, can be positioned anywhere that the ge…
▽ More
Attempts to build a discrete theory for rational maps on the sphere via circle packing have foundered on discretization effects in locating branch points. The authors remove this impediment by introducing generalized branch points. A generalized branch point need no longer be attached to an individual circle, but with the help of chaperones and other devices, can be positioned anywhere that the geometry requires. The effects will be illustrated in images and videos as we use generalized branching to fix flaws in discrete Ahlfors and Weierstrasse functions.
△ Less
Submitted 10 July, 2016;
originally announced July 2016.
-
Cluster growth in the dynamical Erdős-Rényi process with forest fires
Authors:
Edward Crane,
Nic Freeman,
Bálint Tóth
Abstract:
We investigate the growth of clusters within the forest fire model of Ráth and Tóth [22]. The model is a continuous-time Markov process, similar to the dynamical Erdős-Rényi random graph but with the addition of so-called fires. A vertex may catch fire at any moment and, when it does so, causes all edges within its connected cluster to burn, meaning that they instantaneously disappear. Each burned…
▽ More
We investigate the growth of clusters within the forest fire model of Ráth and Tóth [22]. The model is a continuous-time Markov process, similar to the dynamical Erdős-Rényi random graph but with the addition of so-called fires. A vertex may catch fire at any moment and, when it does so, causes all edges within its connected cluster to burn, meaning that they instantaneously disappear. Each burned edge may later reappear.
We give a precise description of the process $C_t$ of the size of the cluster of a tagged vertex, in the limit as the number of vertices in the model tends to infinity. We show that $C_t$ is an explosive branching process with a time-inhomogeneous offspring distribution and instantaneous return to $1$ on each explosion. Additionally, we show that the characteristic curves used to analyse the Smoluchowski-type coagulation equations associated to the model have a probabilistic interpretation in terms of the process $C_t$.
△ Less
Submitted 12 December, 2014; v1 submitted 20 May, 2014;
originally announced May 2014.
-
Intrinsic circle domains
Authors:
Edward Crane
Abstract:
Using quasiconformal mappings, we prove that any Riemann surface of finite connectivity and finite genus is conformally equivalent to an intrinsic circle domain U in a compact Riemann surface S. This means that each connected component B of S \ U is either a point or a closed geometric disc with respect to the complete constant curvature conformal metric of the Riemann surface (U union B). Moreove…
▽ More
Using quasiconformal mappings, we prove that any Riemann surface of finite connectivity and finite genus is conformally equivalent to an intrinsic circle domain U in a compact Riemann surface S. This means that each connected component B of S \ U is either a point or a closed geometric disc with respect to the complete constant curvature conformal metric of the Riemann surface (U union B). Moreover the pair (U,S) is unique up to conformal isomorphisms. We give a generalization to countably infinite connectivity. Finally we show how one can compute numerical approximations to intrinsic circle domains using circle packings and conformal welding.
△ Less
Submitted 4 November, 2013; v1 submitted 27 March, 2013;
originally announced March 2013.
-
Rigidity of configurations of balls and points in the $N$-sphere
Authors:
Edward Crane,
Ian Short
Abstract:
We answer two questions of Beardon and Minda that arose from their study of the conformal symmetries of circular regions in the complex plane. We show that a configuration of closed balls in the $N$-sphere is determined up to Möbius transformations by the signed inversive distances between pairs of its elements, except when the boundaries of the balls have a point in common, and that a configura…
▽ More
We answer two questions of Beardon and Minda that arose from their study of the conformal symmetries of circular regions in the complex plane. We show that a configuration of closed balls in the $N$-sphere is determined up to Möbius transformations by the signed inversive distances between pairs of its elements, except when the boundaries of the balls have a point in common, and that a configuration of points in the $N$-sphere is determined by the absolute cross-ratios of 4-tuples of its elements. The proofs use the hyperboloid model of hyperbolic $(N+1)$-space.
△ Less
Submitted 10 December, 2009;
originally announced December 2009.
-
The simple harmonic urn
Authors:
Edward Crane,
Nicholas Georgiou,
Stanislav Volkov,
Andrew R. Wade,
Robert J. Waters
Abstract:
We study a generalized Pólya urn model with two types of ball. If the drawn ball is red, it is replaced together with a black ball, but if the drawn ball is black it is replaced and a red ball is thrown out of the urn. When only black balls remain, the roles of the colors are swapped and the process restarts. We prove that the resulting Markov chain is transient but that if we throw out a ball eve…
▽ More
We study a generalized Pólya urn model with two types of ball. If the drawn ball is red, it is replaced together with a black ball, but if the drawn ball is black it is replaced and a red ball is thrown out of the urn. When only black balls remain, the roles of the colors are swapped and the process restarts. We prove that the resulting Markov chain is transient but that if we throw out a ball every time the colors swap, the process is recurrent. We show that the embedded process obtained by observing the number of balls in the urn at the swapping times has a scaling limit that is essentially the square of a Bessel diffusion. We consider an oriented percolation model naturally associated with the urn process, and obtain detailed information about its structure, showing that the open subgraph is an infinite tree with a single end. We also study a natural continuous-time embedding of the urn process that demonstrates the relation to the simple harmonic oscillator; in this setting, our transience result addresses an open problem in the recurrence theory of two-dimensional linear birth and death processes due to Kesten and Hutton. We obtain results on the area swept out by the process. We make use of connections between the urn process and birth--death processes, a uniform renewal process, the Eulerian numbers, and Lamperti's problem on processes with asymptotically small drifts; we prove some new results on some of these classical objects that may be of independent interest. For instance, we give sharp new asymptotics for the first two moments of the counting function of the uniform renewal process. Finally, we discuss some related models of independent interest, including a "Poisson earthquakes" Markov chain on the homeomorphisms of the plane.
△ Less
Submitted 16 January, 2012; v1 submitted 2 November, 2009;
originally announced November 2009.
-
Conical limit sets and continued fractions
Authors:
Edward Crane,
Ian Short
Abstract:
Inspired by questions of convergence in continued fraction theory, Erdős, Piranian and Thron studied the possible sets of divergence for arbitrary sequences of Möbius maps acting on the Riemann sphere, $S^2$. By identifying $S^2$ with the boundary of three-dimensional hyperbolic space, $H^3$, we show that these sets of divergence are precisely the sets that arise as conical limit sets of subsets…
▽ More
Inspired by questions of convergence in continued fraction theory, Erdős, Piranian and Thron studied the possible sets of divergence for arbitrary sequences of Möbius maps acting on the Riemann sphere, $S^2$. By identifying $S^2$ with the boundary of three-dimensional hyperbolic space, $H^3$, we show that these sets of divergence are precisely the sets that arise as conical limit sets of subsets of $H^3$. Using hyperbolic geometry, we give simple geometric proofs of the theorems of Erdős, Piranian and Thron that generalise to arbitrary dimensions. New results are also obtained about the class of conical limit sets, for example, that it is closed under locally quasisymmetric homeomorphisms. Applications are given to continued fractions.
△ Less
Submitted 13 August, 2007;
originally announced August 2007.
-
A note on the Hayman-Wu theorem
Authors:
Edward Crane
Abstract:
The Hayman-Wu theorem states that the preimage of a line or circle L under a conformal mapping from the unit disc to a simply-connected domain U has total Euclidean length bounded by an absolute constant. The best possible constant is known to lie in the interval [pi^2, 4 pi), thanks to work of Øyma and Rohde. Earlier, Brown Flinn showed that the total length is at most pi^2 in the special case…
▽ More
The Hayman-Wu theorem states that the preimage of a line or circle L under a conformal mapping from the unit disc to a simply-connected domain U has total Euclidean length bounded by an absolute constant. The best possible constant is known to lie in the interval [pi^2, 4 pi), thanks to work of Øyma and Rohde. Earlier, Brown Flinn showed that the total length is at most pi^2 in the special case in which U contains L. Let r be the anti-Möbius map that fixes L pointwise. In this note we extend the sharp bound pi^2 to the case where each connected component of the intersection of U with r(U) is bounded by one arc of U and its image under r. We also strengthen the bounds slightly by replacing Euclidean length with the strictly larger spherical length restricted to the unit disc.
△ Less
Submitted 12 July, 2007;
originally announced July 2007.
-
Mean Value Conjectures for Rational Maps
Authors:
Edward Crane
Abstract:
Let p be a polynomial in one complex variable. Smale's mean value conjecture estimates |p'(z)| in terms of the gradient of a chord from (z, p(z)) to some stationary point on the graph of $p$. The conjecture does not immediately generalise to rational maps since its formulation is invariant under the group of affine maps, not the full Mobius group. Here we give two possible generalisations to rat…
▽ More
Let p be a polynomial in one complex variable. Smale's mean value conjecture estimates |p'(z)| in terms of the gradient of a chord from (z, p(z)) to some stationary point on the graph of $p$. The conjecture does not immediately generalise to rational maps since its formulation is invariant under the group of affine maps, not the full Mobius group. Here we give two possible generalisations to rational maps, both of which are Mobius invariant. In both cases we prove a version with a weaker constant, in parallel to the situation for Smale's mean value conjecture. Finally, we discuss some candidate extremal rational maps, namely rational maps all of whose critical points are fixed points.
△ Less
Submitted 26 November, 2004;
originally announced November 2004.
-
The Area of Polynomial Images and Preimages
Authors:
Edward Crane
Abstract:
Let p be a monic polynomial in one complex variable and K a measurable subset of the complex plane. In terms of the area of K, we give an upper bound on the area of the preimage of K under p and a lower bound on the area of the image of K under p, (counted with multiplicity). Both bounds are sharp. The former extends an inequality of Polya. The proof uses Carleman's isoperimetric inequality for…
▽ More
Let p be a monic polynomial in one complex variable and K a measurable subset of the complex plane. In terms of the area of K, we give an upper bound on the area of the preimage of K under p and a lower bound on the area of the image of K under p, (counted with multiplicity). Both bounds are sharp. The former extends an inequality of Polya. The proof uses Carleman's isoperimetric inequality for plane condensers. We include a summary of the necessary potential theory.
△ Less
Submitted 17 February, 2003;
originally announced February 2003.