Researcher profile

N. Anders Petersson

N. Anders Petersson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2022arXiv

Binary Optimal Control Of Single-Flux-Quantum Pulse Sequences

We introduce a binary, relaxed gradient, trust-region method for optimizing pulse sequences for single flux quanta (SFQ) control of a quantum computer. The pulse sequences are optimized with the goal of realizing unitary gate transformations. Each pulse has a fixed amplitude and duration. We model this process as an binary optimal control problem, constrained by Schrödinger's equation, where the binary variables indicate whether each pulse is on or off. We introduce a first-order trust-region method, which takes advantage of a relaxed gradient to determine an optimal pulse sequence that minimizes the gate infidelity, while also suppressing leakage to higher energy levels. The proposed algorithm has a computational complexity of ${\cal O}(p\log(p)$, where $p$ is the number of pulses in the sequence. We present numerical results for the H and X gates, where the optimized pulse sequences give gate fidelity's better than $99.9\%$, in $\approx 25$ trust-region iterations.

preprint2022arXiv

Optimal Control of Closed Quantum Systems via B-Splines with Carrier Waves

We consider the optimal control problem of determining electromagnetic pulses for implementing logical gates in a closed quantum system, where the Hamiltonian models the dynamics of coupled superconducting qudits. The quantum state is governed by Schrödinger's equation, which we formulate in terms of the real and imaginary parts of the state vector and solve by the Störmer-Verlet scheme, which is a symplectic partitioned Runge-Kutta method. A novel parameterization of the control functions based on B-splines with carrier waves is introduced. The carrier waves are used to trigger the resonant frequencies in the system Hamiltonian, and the B-spline functions specify their amplitude and phase. This approach allows the number of control parameters to be independent of, and significantly smaller than, the number of time steps for integrating Schrödinger's equation. We present numerical examples of how the proposed technique can be combined with an interior point L-BFGS algorithm for realizing quantum gates, and generalize our approach to calculate risk-neutral controls that are resilient to noise in the Hamiltonian model. The proposed method is also shown to compare favorably with QuTiP/pulse\_optim and Grape-Tensorflow.

preprint2020arXiv

Energy conservative SBP discretizations of the acoustic wave equation in covariant form on staggered curvilinear grids

We develop a numerical method for solving the acoustic wave equation in covariant form on staggered curvilinear grids in an energy conserving manner. The use of a covariant basis decomposition leads to a rotationally invariant scheme that outperforms a Cartesian basis decomposition on rotated grids. The discretization is based on high order Summation-By-Parts (SBP) operators and preserves both symmetry and positive definiteness of the contravariant metric tensor. To improve accuracy and decrease computational cost, we also derive a modified discretization of the metric tensor that leads to a conditionally stable discretization. Bounds are derived that yield a point-wise condition that can be evaluated to check for stability of the modified discretization. This condition shows that the interpolation operators should be constructed such that their norm is close to one.

preprint2011arXiv

Boundary estimates for the elastic wave equation in almost incompressible materials

We study the half-plane problem for the elastic wave equation subject to a free surface boundary condition, with particular emphasis on almost incompressible materials. A normal mode analysis is developed to estimate the solution in terms of the boundary data, showing that the problem is boundary stable. The dependence on the material properties, which is difficult to analyze by the energy method, is made transparent by our estimates. The normal mode technique is used to analyze the influence of truncation errors in a finite difference approximation. Our analysis explains why the number of grid points per wave length must be increased when the shear modulus ($μ$) becomes small, that is, for almost incompressible materials. To obtain a fixed error in the phase velocity of Rayleigh surface waves as $μ\to 0$, our analysis predicts that the grid size must be proportional to $μ^{1/2}$ for a second order method. For a fourth order method, the grid size can be proportional to $μ^{1/4}$. Numerical experiments confirm this scaling and illustrate the superior efficiency of the fourth order method.

preprint2010arXiv

Initial-boundary value problems for second order systems of partial differential equations

We develop a well-posedness theory for second order systems in bounded domains where boundary phenomena like glancing and surface waves play an important role. Attempts have previously been made to write a second order system consisting of n equations as a larger first order system. Unfortunately, the resulting first order system consists, in general, of more than 2n equations which leads to many complications, such as side conditions which must be satisfied by the solution of the larger first order system. Here we will use the theory of pseudo-differential operators combined with mode analysis. There are many desirable properties of this approach: 1) The reduction to first order systems of pseudo-differential equations poses no difficulty and always gives a system of 2n equations. 2) We can localize the problem, i.e., it is only necessary to study the Cauchy problem and halfplane problems with constant coefficients. 3) The class of problems we can treat is much larger than previous approaches based on "integration by parts". 4) The relation between boundary conditions and boundary phenomena becomes transparent.