Achieving high-fidelity single-qubit gates, two-qubit gates, and qubit readout is critical for building scalable, error-corrected quantum computers. However, device parameters that enhance one operation often degrade the others, making simultaneous optimization challenging. Here, we demonstrate that careful tuning of qubit-coupler coupling strengths in a superconducting circuit with two transmon qubits coupled via a tunable coupler enables high-fidelity single- and two-qubit gates, without compromising readout performance. Furthermore, we introduce a new calibration protocol for diabatic controlled-Z gates, phased-averaged leakage error amplification (PALEA), which enables efficient suppression of coherent gate errors and leakage to noncomputational states. Using PALEA in combination with careful tuning of qubit-coupler coupling strengths, we achieve a 40 h-averaged CZ gate fidelity of 99.93 %, simultaneous single-qubit gate fidelities of 99.98 %, and readout fidelities over 99.94 % in a single device. Our results demonstrate a viable path toward scaling up superconducting quantum processors while maintaining consistently high gate and readout fidelities.

We experimentally demonstrate the fast generation of thermal states of a transmon using a single-junction quantum-circuit refrigerator (QCR) as an in-situ-tunable environment. Through single-shot readout, we monitor the transmon up to its third-excited state, assessing population distributions controlled by QCR drive pulses. Whereas cooling can be achieved in the weak-drive regime, high-amplitude pulses can generate Boltzmann-distributed populations from a temperature of 110 mK up to 500 mK within 100 ns. As we propose in our work, this fast and efficient temperature control provides an appealing opportunity to demonstrate a quantum heat engine. Our results also pave the way for efficient dissipative state preparation and for reducing the circuit depth in thermally assisted quantum algorithms and quantum annealing.

We present a differential geometry-based model for linear wave equations in (2+1)-dimensional spacetime. This model encompasses acoustic, elastic, and electromagnetic waves and is also applicable in quantum mechanical simulations. For discretization, we introduce a spacetime extension of discrete exterior calculus, resulting in a leapfrog-style time evolution. The scheme further supports numerical simulations of moving and deforming domains. The numerical tests presented in this paper demonstrate the method’s stability limits and computational efficiency.

We experimentally demonstrate a recently proposed single-junction quantum-circuit refrigerator (QCR) as an in-situ-tunable low-temperature environment for a superconducting 4.7-GHz resonator. With the help of a transmon qubit, we measure the populations of the different resonator Fock states, thus providing reliable access to the temperature of the engineered electromagnetic environment and its effect on the resonator. We demonstrate coherent and thermal resonator states and that the on-demand dissipation provided by the QCR can drive these to a small fraction of a photon on average, even if starting above 1 K. We observe that the QCR can be operated either with a dc bias voltage or a gigahertz rf drive, or a combination of these. The bandwidth of the rf drive is not limited by the circuit itself and consequently, we show that 2.9-GHz continuous and 10-ns-pulsed drives lead to identical desired refrigeration of the resonator. These observations answer to the shortcomings of previous works where the Fock states were not resolvable and the QCR exhibited slow charging dynamics. Thus this work introduces a versatile tool to study open quantum systems, quantum thermodynamics, and to quickly reset superconducting qubits.

Quantum computing has tremendous potential to overcome some of the fundamental limitations present in classical information processing. Yet, today's technological limitations in the quality and scaling prevent exploiting its full potential. Quantum computing based on superconducting quantum processing units (QPUs) is among the most promising approaches towards practical quantum advantage. In this article the basic technological approach of IQM Quantum Computers is described covering both the QPU and the rest of the full-stack quantum computer. In particular, the focus is on a 20-qubit quantum computer featuring the Garnet QPU and its architecture, which we will scale up to 150 qubits. We also present QPU and system-level benchmarks, including a median 2-qubit gate fidelity of 99.5 Greenberger-Horne-Zeilinger (GHZ) state.

In this paper, we apply the exact controllability concept for time-harmonic electromagnetic scattering. The problem is presented in terms of the differential forms, and the discrete exterior calculus is utilized for spatial discretization. Accordingly, the physical properties of the problem are maintained. Despite we consider time-harmonic problems, we concentrate on transient wave equations treated by the exact controllability approach. Essentially, we use a controlled variation of the asymptotic approach with periodic constraints, in which the time-dependent equation is simulated in time, until the time-harmonic solution is reached.

We study the scattering properties of a cloud of particles. The particles are spherical, close to the incident wavelength in size, have a high albedo, and are randomly packed to 20% volume density. We show, using both numerically exact methods for solving the Maxwell equations and radiative-transfer-approximation methods, that the scattering properties of the cloud converge after about ten million particles in the system. After that, the backward-scattered properties of the system should estimate the properties of a macroscopic, practically infinite system. (C) 2021 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/)

Finite difference kind of schemes are popular in approximating wave propagation problems in finite dimensional spaces. While Yee’s original paper on the finite difference method is already from the sixties, mathematically there still remains questions which are not yet satisfactorily covered. In this paper, we address two issues of this kind. Firstly, in the literature Yee’s scheme is constructed separately for each particular type of wave problem. Here, we explicitly generalize the Yee scheme to a class of wave problems that covers at large physics field theories. For this we introduce Yee’s scheme for all problems of a class characterised on a Minkowski manifold by (i) a pair of first order partial differential equations and by (ii) a constitutive relation that couple the differential equations with a Hodge relation. In addition, we introduce a strategy to systematically exploit higher order Whitney elements in Yee-like approaches. This makes higher order interpolation possible both in time and space. For this, we show that Yee-like schemes preserve the local character of the Hodge relation, which is to say, the constitutive laws become imposed on a finite set of points instead of on all ordinary points of space. As a result, the usage of higher order Whitney forms does not compel to change the actual solution process at all. This is demonstrated with a simple example.

Ill-posed inversion problems are one of the major challenges when there is a need to combine measurements with the theory and numerical model. In this study, we demonstrate the use of wave propagation simulations to train a convolutional neural network (CNN) for retrieving sub-wavelength thickness profiles of thin film coatings from hyperspectral images. The simulations are produced by solving numerically one-dimensional wave equation with a method based on Discrete Exterior Calculus (DEC). This approach provides a powerful tool to produce large sets of training data for the neural network. CNN was verified by simulated verification sets and measured reflectance spectra, both of which showed strong correlations. A hyperspectral image that cover a region of sample provides sufficient number of spectra for reliable thickness analysis, but at the same time allows the use of a small detection spots to solve non-uniformity problems. The non-uniformity of film thickness is characterized and the results are promising. The approach introduced in this study provides a potential solution to the challenges of thin film analytics in the field of sub-wavelength thickness and its non-uniformity.

A general conservation law that defines a class of physical field theories is constructed. First, the notion of a general field is introduced as a formal sum of differential forms on a Minkowski manifold. By the action principle the conservation law is defined for such a general field. By construction, particular field notions of physics, e.g., magnetic flux, electric field strength, stress, strain etc. become instances of the general field. Hence, the differential equations that constitute physical field theories become also instances of the general conservation law. The general field and the general conservation law together correspond to a large class of relativistic hyperbolic physical field models. The parabolic and elliptic models can thereafter be derived by adding constraints. The approach creates solid foundations for developing software systems for scientific computing; the unifying structure shared by the class of field models makes it possible to implement software systems which are not restricted to certain predefined problems. The versatility of the proposed approach is demonstrated by numerical experiments with moving and deforming domains.

Light scattering by particles large compared to the wavelength of incident light is traditionally solved using ray optics which considers absorption inside the particle approximately, along the ray paths. To study the effects rising from this simplification, we have updated the ray-optics code SIRIS to take into account the propagation of light as inhomogeneous plane waves inside an absorbing particle. We investigate the impact of this correction on traditional ray-optics computations in the example case of light scattering by ice crystals through the extended near-infrared (NIR) wavelength regime. In this spectral range, ice changes from nearly transparent to opaque, and therefore provides an interesting test case with direct connection and applicability to atmospheric remote-sensing measurements at NIR wavelengths. We find that the correction for inhomogeneous waves systematically increases the single-scattering albedo throughout the NIR spectrum for both randomly-oriented, column-like hexagonal crystals and ice crystals shaped like Gaussian random spheres. The largest increase in the single-scattering albedo is 0.042 for hexagonal crystals and 0.044 for Gaussian random spheres, both at lambda = 2.725 mu m. Although the effects on the 4 x 4 scattering-matrix elements are generally small, the largest differences are seen at 2.0 mu m and 3.969 mu m wavelengths where the correction for inhomogeneous waves affects mostly the backscattering hemisphere of the depolarization-connected P-22/P-11, P-33/P-11, and P-44/P-11. We evaluated the correction for inhomogeneous waves through comparisons against the discrete exterior calculus (DEC) method. We computed scattering by hexagonal ice crystals using the DEC, a traditional ray-optics code (SIRIS3), and a ray-optics code with inhomogeneous waves (SIRIS4). Comparisons of the scattering-matrix elements from SIRIS3 and SIRIS4 against those from the DEC suggest that consideration of the inhomogeneous waves brings the ray-optics solution generally closer to the exact result and, therefore, should be taken into account in scattering by absorbing particles large compared to the wavelength of incident light. (C) 2018 The Authors. Published by Elsevier Ltd.

We introduce a general class of second-order boundary value problems unifying application areas such as acoustics, electromagnetism, elastodynamics, quantum mechanics, and so on, into a single framework. This also enables us to solve wave propagation problems very efficiently with a single software system. The solution method precisely follows the conservation laws in finite-dimensional systems, whereas the constitutive relations are imposed approximately. We employ discrete exterior calculus for the spatial discretization, use natural crystal structures for three-dimensional meshing, and derive a “discrete Hodge” adapted to harmonic wave. The numerical experiments indicate that the cumulative pollution error can be practically eliminated in the case of harmonic wave problems. The restrictions following from the CFL condition can be bypassed with a local time-stepping scheme. The computational savings are at least one order of magnitude.

Jukka Räbinä,1, ∗ Pekko Kuopanportti,2 Markus Kivioja,1 Mikko Möttönen,3 and Tuomo Rossi1 University of Jyvaskyla, Faculty of Information Technology, P.O. Box 35, FI-40014 University of Jyvaskyla, Finland Department of Physics, University of Helsinki, P.O. Box 43, FI-00014 Helsinki, Finland QCD Labs, QTF Centre of Excellence, Department of Applied Physics, Aalto University, P.O. Box 13500, FI-00076 AALTO, Finland. (Dated: October 5, 2018)

We study the splitting dynamics of giant vortices in dilute Bose-Einstein condensates by numerically integrating the three-dimensional Gross-Pitaevskii equation in time. By taking advantage of tetrahedral tiling in the spatial discretization, we decrease the error and increase the reliability of the numerical method. An extensive survey of vortex splitting patterns is presented for different aspect ratios of the harmonic trapping potential. The discrete rotational symmetries of the splitting patterns that emerge in the time evolution are in good agreement with predictions obtained by solving the prevailing dynamical instabilities from the Bogoliubov equations. Furthermore, we observe intertwining of the split vortices in prolate condensates and a split-and-revival phenomenon in a spherical condensate.

We model meteoroids entering the Earth׳s atmosphere as objects surrounded by non-magnetized plasma, and consider efficient numerical simulation of radar reflections from meteors in the time domain. Instead of the widely used finite difference time domain method (FDTD), we use more generalized finite differences by applying the discrete exterior calculus (DEC) and non-uniform leapfrog-style time discretization. The computational domain is presented by convex polyhedral elements. The convergence of the time integration is accelerated by the exact controllability method. The numerical experiments show that our code is efficiently parallelized. The DEC approach is compared to the volume integral equation (VIE) method by numerical experiments. The result is that both methods are competitive in modelling non-magnetized plasma scattering. For demonstrating the simulation capabilities of the DEC approach, we present numerical experiments of radar reflections and vary parameters in a wide range.

We apply high-quality discretizations to simulate electromagnetic microwaves. Instead of the vector field presentations, we focus on differential forms and discretize the model in the spatial domain using the discrete exterior calculus. At the discrete level, both the Hodge operators and the time discretization are optimized for time-harmonic simulations. Non-uniform spatial and temporal discretization are applied in problems in which the wavelength is highly-variable and geometry contains sub-wavelength structures.

We consider the computationally efficient time integration of Maxwell's equations using discrete exterior calculus (DEC) as the computational framework. With the theory of DEC, we associate the degrees of freedom of the electric and magnetic fields with primal and dual mesh structures, respectively. We concentrate on mesh constructions that imitate the geometry of the close packing in crystal lattices that is typical of elemental metals and intermetallic compounds. This class of computational grids has not been used previously in electromagnetics. For the simulation of wave propagation driven by time-harmonic source terms, we provide an optimized Hodge operator and a novel time discretization scheme with nonuniform time step size. The numerical experiments show a significant improvement in accuracy and a decrease in computing time compared to simulations with well-known variants of the finite difference time domain method.