skip to main content


Title: An updated LLVM-based quantum research compiler with further OpenQASM support
Abstract

Quantum computing is a rapidly growing field with the potential to change how we solve previously intractable problems. Emerging hardware is approaching a complexity that requires increasingly sophisticated programming and control. Scaffold is an older quantum programming language that was originally designed for resource estimation for far-future, large quantum machines, and ScaffCC is the corresponding LLVM-based compiler. For the first time, we provide a full and complete overview of the language itself, the compiler as well as its pass structure. While previous works Abhariet al(2015Parallel Comput.452–17), Abhariet al(2012 Scaffold: quantum programming languagehttps://cs.princeton.edu/research/techreps/TR-934-12), have piecemeal descriptions of different portions of this toolchain, we provide a more full and complete description in this paper. We also introduce updates to ScaffCC including conditional measurement and multidimensional qubit arrays designed to keep in step with modern quantum assembly languages, as well as an alternate toolchain targeted at maintaining correctness and low resource count for noisy-intermediate scale quantum (NISQ) machines, and compatibility with current versions of LLVM and Clang. Our goal is to provide the research community with a functional LLVM framework for quantum program analysis, optimization, and generation of executable code.

 
more » « less
NSF-PAR ID:
10303668
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
Quantum Science and Technology
Volume:
5
Issue:
3
ISSN:
2058-9565
Page Range / eLocation ID:
Article No. 034013
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Motivated by the challenge of sampling Gibbs measures with nonconvex potentials, we study a continuum birth–death dynamics. We improve results in previous works (Liuet al2023Appl. Math. Optim.8748; Luet al2019 arXiv:1905.09863) and provide weaker hypotheses under which the probability density of the birth–death governed by Kullback–Leibler divergence or byχ2divergence converge exponentially fast to the Gibbs equilibrium measure, with a universal rate that is independent of the potential barrier. To build a practical numerical sampler based on the pure birth–death dynamics, we consider an interacting particle system, which is inspired by the gradient flow structure and the classical Fokker–Planck equation and relies on kernel-based approximations of the measure. Using the technique of Γ-convergence of gradient flows, we show that on the torus, smooth and bounded positive solutions of the kernelised dynamics converge on finite time intervals, to the pure birth–death dynamics as the kernel bandwidth shrinks to zero. Moreover we provide quantitative estimates on the bias of minimisers of the energy corresponding to the kernelised dynamics. Finally we prove the long-time asymptotic results on the convergence of the asymptotic states of the kernelised dynamics towards the Gibbs measure.

     
    more » « less
  2. Abstract

    We study nonlinear optimization problems with a stochastic objective and deterministic equality and inequality constraints, which emerge in numerous applications including finance, manufacturing, power systems and, recently, deep neural networks. We propose an active-set stochastic sequential quadratic programming (StoSQP) algorithm that utilizes a differentiable exact augmented Lagrangian as the merit function. The algorithm adaptively selects the penalty parameters of the augmented Lagrangian, and performs a stochastic line search to decide the stepsize. The global convergence is established: for any initialization, the KKT residuals converge to zeroalmost surely. Our algorithm and analysis further develop the prior work of Na et al. (Math Program, 2022.https://doi.org/10.1007/s10107-022-01846-z). Specifically, we allow nonlinear inequality constraintswithoutrequiring the strict complementary condition; refine some of designs in Na et al. (2022) such as the feasibility error condition and the monotonically increasing sample size; strengthen the global convergence guarantee; and improve the sample complexity on the objective Hessian. We demonstrate the performance of the designed algorithm on a subset of nonlinear problems collected in CUTEst test set and on constrained logistic regression problems.

     
    more » « less
  3. Abstract

    To quantitatively convert upper mantle seismic wave speeds measured into temperature, density, composition, and corresponding and uncertainty, we introduce theWhole‐rockInterpretativeSeismicToolboxForUltramaficLithologies (WISTFUL). WISTFUL is underpinned by a database of 4,485 ultramafic whole‐rock compositions, their calculated mineral modes, elastic moduli, and seismic wave speeds over a range of pressure (P) and temperature (T) (P = 0.5–6 GPa,T = 200–1,600°C) using the Gibbs free energy minimization routine Perple_X. These data are interpreted with a toolbox of MATLAB® functions, scripts, and three general user interfaces:WISTFUL_relations, which plots relationships between calculated parameters and/or composition;WISTFUL_geotherms, which calculates seismic wave speeds along geotherms; andWISTFUL_inversion, which inverts seismic wave speeds for best‐fit temperature, composition, and density. To evaluate our methodology and quantify the forward calculation error, we estimate two dominant sources of uncertainty: (a) the predicted mineral modes and compositions, and (b) the elastic properties and mixing equations. To constrain the first source of uncertainty, we compiled 122 well‐studied ultramafic xenoliths with known whole‐rock compositions, mineral modes, and estimatedPTconditions. We compared the observed mineral modes with modes predicted using five different thermodynamic solid solution models. The Holland et al. (2018,https://doi.org/10.1093/petrology/egy048) solution models best reproduce phase assemblages (∼12 vol. % phase root‐mean‐square error [RMSE]) and estimated wave speeds. To assess the second source of uncertainty, we compared wave speed measurements of 40 ultramafic rocks with calculated wave speeds, finding excellent agreement (VpRMSE = 0.11 km/s). WISTFUL easily analyzes seismic datasets, integrates into modeling, and acts as an educational tool.

     
    more » « less
  4. Abstract

    “Classical shadows” are estimators of an unknown quantum state, constructed from suitably distributed random measurements on copies of that state (Huang et al. in Nat Phys 16:1050, 2020,https://doi.org/10.1038/s41567-020-0932-7). In this paper, we analyze classical shadows obtained using random matchgate circuits, which correspond to fermionic Gaussian unitaries. We prove that the first three moments of the Haar distribution over thecontinuousgroup of matchgate circuits are equal to those of thediscreteuniform distribution over only the matchgate circuits that are also Clifford unitaries; thus, the latter forms a “matchgate 3-design.” This implies that the classical shadows resulting from the two ensembles are functionally equivalent. We show how one can use these matchgate shadows to efficiently estimate inner products between an arbitrary quantum state and fermionic Gaussian states, as well as the expectation values of local fermionic operators and various other quantities, thus surpassing the capabilities of prior work. As a concrete application, this enables us to apply wavefunction constraints that control the fermion sign problem in the quantum-classical auxiliary-field quantum Monte Carlo algorithm (QC-AFQMC) (Huggins et al. in Nature 603:416, 2022,https://doi.org/10.1038/s41586-021-04351-z), without the exponential post-processing cost incurred by the original approach.

     
    more » « less
  5. Abstract

    We mapped the Soi crater region at 1:800,000 scale and produced a geomorphological map using methodology presented by Malaska, Lopes, Williams, et al. (2016),https://doi.org/10.1016/j.icarus.2016.02.021and Schoenfeld et al. (2021),https://doi.org/10.1016/j.icarus.2021.114516. This region spans longitude 110° to 180°W and latitude 0° to 60°N and is representative of the transition between the equatorial, mid‐latitude, and high‐latitude northern regions of Titan. We used Cassini Synthetic Aperture Radar (SAR) as our primary mapping data set. For areas where SAR was not available, we used lower resolution data from the Imaging Science Subsystem, the Visible and Infrared Mapping Spectrometer, radiometry, and high‐altitude SAR for complete mapping coverage of the region. We identified 22 geomorphological units, 3 of which have been discussed in existing literature but have not yet been incorporated into our mapping investigations. These include sharp‐edged depressions (bse), ramparts (brh), and bright gradational plains (pgh). All six major terrain classes are represented in this region: Craters, Labyrinth, Hummocky/mountainous, Plains, Dunes, and Basin and Lakes. We find that plains dominate the surface of the Soi crater region, comprising ∼73% of the mapped area, followed by dunes (∼14%), mountains/hummocky terrains (∼12%), basin and lakes (∼0.7%), labyrinth terrains (∼0.5%), and crater terrains (∼0.4%). We also observe empty lakes as far south as 40°N. The Soi crater region largely has the same collection and proportion of geomorphological units to other mapped regions on Titan. These results further support the hypothesis that surface processes are, broadly speaking, the same across Titan's middle and equatorial latitudes, with the exception of Xanadu.

     
    more » « less