Researcher profile

Kristel Michielsen

Kristel Michielsen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
13topics
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

16 published item(s)

preprint2026arXiv

Empirical Quantum Advantage in Constrained Optimization from Encoded Unitary Designs

We introduce the Constraint-Enhanced Quantum Approximate Optimization Algorithm (CE-QAOA), a shallow, constraint-aware ansatz that operates inside the one-hot product space [n]^m, where m is the number of blocks and each block is initialized in an n-qubit W_n state. We give an ancilla-free, depth-optimal encoder that prepares W_n using n-1 two-qubit rotations per block, and a two-local block-XY mixer that preserves the one-hot manifold and has a constant spectral gap on the one-excitation sector. At the level of expressivity, we establish per-block controllability, implying approximate universality per block. At the level of distributional behavior, we show that, after natural block and symbol permutation twirls, shallow CE-QAOA realizes an encoded unitary 1-design and supports approximate second-moment (2-design) behavior; combined with a Paley-Zygmund argument, this yields finite-shot anticoncentration guarantees. Algorithmically, we wrap constant-depth sampling with a deterministic feasibility checker to obtain a polynomial-time hybrid quantum-classical solver (PHQC) that returns the best observed feasible solution in O(S n^2) time, where S is a polynomial shot budget. We obtain two advantages. First, when CE-QAOA fixes r >= 1 locations different from the start city, we achieve a Theta(n^r) reduction in shot complexity even against a classical sampler that draws uniformly from the feasible set. Second, against a classical baseline restricted to raw bitstring sampling, we show an exp(Theta(n^2)) minimax separation. In noiseless circuit simulations of traveling salesman problem instances with n in {4,...,10} locations from the QOPTLib benchmark library, we recover the global optimum at depth p = 1 using polynomial shot budgets and coarse parameter grids defined by the problem size.

preprint2022arXiv

Assessment of the variational quantum eigensolver: application to the Heisenberg model

We present and analyze large-scale simulation results of a hybrid quantum-classical variational method to calculate the ground state energy of the anti-ferromagnetic Heisenberg model. Using a massively parallel universal quantum computer simulator, we observe that a low-depth-circuit ansatz advantageously exploits the efficiently preparable Néel initial state, avoids potential barren plateaus, and works for both one- and two-dimensional lattices. The analysis reflects the decisive ingredients required for a simulation by comparing different ansätze, initial parameters, and gradient-based versus gradient-free optimizers. Extrapolation to the thermodynamic limit accurately yields the analytical value for the ground state energy, given by the Bethe ansatz. We predict that a fully functional quantum computer with 100 qubits can calculate the ground state energy with a relatively small error.

preprint2022arXiv

Benchmarking Advantage and D-Wave 2000Q quantum annealers with exact cover problems

We benchmark the quantum processing units of the largest quantum annealers to date, the 5000+ qubit quantum annealer Advantage and its 2000+ qubit predecessor D-Wave 2000Q, using tail assignment and exact cover problems from aircraft scheduling scenarios. The benchmark set contains small, intermediate, and large problems with both sparsely connected and almost fully connected instances. We find that Advantage outperforms D-Wave 2000Q for almost all problems, with a notable increase in success rate and problem size. In particular, Advantage is also able to solve the largest problems with 120 logical qubits that D-Wave 2000Q cannot solve anymore. Furthermore, problems that can still be solved by D-Wave 2000Q are solved faster by Advantage. We find, however, that D-Wave 2000Q can achieve better success rates for sparsely connected problems that do not require the many new couplers present on Advantage, so improving the connectivity of a quantum annealer does not per se improve its performance.

preprint2022arXiv

Classical, quantum and event-by-event simulation of a Stern-Gerlach experiment with neutrons

We present a comprehensive simulation study of the Newtonian and quantum model of a Stern-Gerlach experiment with cold neutrons.By solving Newton's equation of motion and the time-dependent Pauli equation, for a wide range of uniform magnetic field strengths, we scrutinize the role of the latter for drawing the conclusion that the magnetic moment of the neutron is quantized. We then demonstrate that a marginal modification of the Newtonian model suffices to construct, without invoking any concept of quantum theory, an event-based subquantum model that eliminates the shortcomings of the classical model and yields results that are in qualitative agreement with experiment and quantum theory. In this event-by-event model, the intrinsic angular momentum can take any value on the sphere, yet, for a sufficiently strong uniform magnetic field, the particle beam splits in two, exactly as in experiment and in concert with quantum theory.

preprint2022arXiv

GPU-accelerated simulations of quantum annealing and the quantum approximate optimization algorithm

We study large-scale applications using a GPU-accelerated version of the massively parallel Jülich universal quantum computer simulator (JUQCS--G). First, we benchmark JUWELS Booster, a GPU cluster with 3744 NVIDIA A100 Tensor Core GPUs. Then, we use JUQCS--G to study the relation between quantum annealing (QA) and the quantum approximate optimization algorithm (QAOA). We find that a very coarsely discretized version of QA, termed approximate quantum annealing (AQA), performs surprisingly well in comparison to the QAOA. It can either be used to initialize the QAOA, or to avoid the costly optimization procedure altogether. Furthermore, we study the scaling of the success probability when using AQA for problems with 30 to 40 qubits. We find that the case with the largest discretization error scales most favorably, surpassing the best result obtained from the QAOA.

preprint2022arXiv

Lecture Notes: Programming Quantum Computers

Quantum computing is a new emerging computer technology. Current quantum computing devices are at a development stage where they are gradually becoming suitable for small real-world applications. This lecture is devoted to the practical aspects of programming such quantum computing devices. The first part of these lecture notes focuses on programming gate-based quantum computers, and the second part shows how to program quantum annealers.

preprint2022arXiv

Numerical analysis of effective models for flux-tunable transmon systems

Simulations and analytical calculations that aim to describe flux-tunable transmons are usually based on effective models of the corresponding lumped-element model. However, when a control pulse is applied, in most cases it is not known how much the predictions made with the effective models deviate from the predictions made with the original lumped-element model. In this work we compare the numerical solutions of the time-dependent Schrödinger equation for both the effective and the lumped-element models, for microwave and unimodal control pulses (external fluxes). These control pulses are used to model single-qubit (X) and two-qubit gate (Iswap and Cz) transitions. First, we derive a non-adiabatic effective Hamiltonian for a single flux-tunable transmon and compare the pulse response of this model to the one of the corresponding circuit Hamiltonian. Here we find that both models predict similar outcomes for similar control pulses. Then, we study how different approximations affect single-qubit (X) and two-qubit gate (Iswap and Cz) transitions in two different two-qubit systems. For this purpose we consider three different systems in total: a single flux-tunable transmon and two two-qubit systems. In summary, we find that a series of commonly applied approximations (individually and/or in combination) can change the response of a system substantially, when a control pulse is applied.

preprint2022arXiv

Quantum annealing for hard 2-SAT problems : Distribution and scaling of minimum energy gap and success probability

In recent years, quantum annealing has gained the status of being a promising candidate for solving various optimization problems. Using a set of hard 2-satisfiabilty (2-SAT) problems, consisting of upto 18-variables problems, we analyze the scaling complexity of the quantum annealing algorithm and study the distributions of the minimum energy gap and the success probability. We extend the analysis of the standard quantum annealing Hamiltonian by introducing an additional term, the trigger Hamiltonian, which can be of two types : ferromagnetic and antiferromagnetic. We use these trigger Hamiltonians to study their influence on the success probability for solving the selected 2-SAT problems. We found that although the scaling of the run-time is exponential for the standard and modified quantum annealing Hamiltonians, the scaling constant in case of adding the trigger Hamiltonians can be significantly smaller. Furthermore, certain choices for the trigger Hamiltonian and annealing times can result in a better scaling than that for simulated annealing. Lastly, we also use the quantum annealers of D-Wave Systems Inc. to study their performance in solving the 2-SAT problems and compare it with the simulation results.

preprint2021arXiv

Support vector machines on the D-Wave quantum annealer

Kernel-based support vector machines (SVMs) are supervised machine learning algorithms for classification and regression problems. We introduce a method to train SVMs on a D-Wave 2000Q quantum annealer and study its performance in comparison to SVMs trained on conventional computers. The method is applied to both synthetic data and real data obtained from biology experiments. We find that the quantum annealer produces an ensemble of different solutions that often generalizes better to unseen data than the single global minimum of an SVM trained on a conventional computer, especially in cases where only limited training data is available. For cases with more training data than currently fits on the quantum annealer, we show that a combination of classifiers for subsets of the data almost always produces stronger joint classifiers than the conventional SVM for the same parameters.

preprint2020arXiv

Benchmarking the Quantum Approximate Optimization Algorithm

The performance of the quantum approximate optimization algorithm is evaluated by using three different measures: the probability of finding the ground state, the energy expectation value, and a ratio closely related to the approximation ratio. The set of problem instances studied consists of weighted MaxCut problems and 2-satisfiability problems. The Ising model representations of the latter possess unique ground states and highly-degenerate first excited states. The quantum approximate optimization algorithm is executed on quantum computer simulators and on the IBM Q Experience. Additionally, data obtained from the D-Wave 2000Q quantum annealer is used for comparison, and it is found that the D-Wave machine outperforms the quantum approximate optimization algorithm executed on a simulator. The overall performance of the quantum approximate optimization algorithm is found to strongly depend on the problem instance.

preprint2020arXiv

Discrete-event simulation of an extended Einstein-Podolsky-Rosen-Bohm experiment

We use discrete-event simulation to construct a subquantum model that can reproduce the quantum-theoretical prediction for the statistics of data produced by the Einstein-Podolsky-Rosen-Bohm experiment and an extension thereof. This model satisfies Einstein's criterion of locality and generates data in an event-by-event and cause-and-effect manner. We show that quantum theory can describe the statistics of the simulation data for a certain range of model parameters only.

preprint2020arXiv

Discrete-event simulation of quantum walks

We use discrete-event simulation on a digital computer to study two different models of experimentally realizable quantum walks. The simulation models comply with Einstein locality, are as "realistic" as the one of the simple random walk in that the particles follow well-defined trajectories, are void of concepts such as particle-wave duality and wave-function collapse, and reproduce the quantum-theoretical results by means of a cause-and-effect, event-by-event process. Our simulation model for the quantum walk experiment presented in [C. Robens et al., Phys. Rev. X 5, 011003 (2015)] reproduces the result of that experiment. Therefore, the claim that the result of the experiment "rigorously excludes (i.e., falsifies) any explanation of quantum transport based on classical, well-defined trajectories" needs to be revised.

preprint2020arXiv

Exponential damping induced by random and realistic perturbations

Given a quantum many-body system and the expectation-value dynamics of some operator, we study how this reference dynamics is altered due to a perturbation of the system's Hamiltonian. Based on projection operator techniques, we unveil that if the perturbation exhibits a random-matrix structure in the eigenbasis of the unperturbed Hamiltonian, then this perturbation effectively leads to an exponential damping of the original dynamics. Employing a combination of dynamical quantum typicality and numerical linked cluster expansions, we demonstrate that our theoretical findings for random matrices can, in some cases, be relevant for the dynamics of realistic quantum many-body models as well. Specifically, we study the decay of current autocorrelation functions in spin-$1/2$ ladder systems, where the rungs of the ladder are treated as a perturbation to the otherwise uncoupled legs. We find a convincing agreement between the exact dynamics and the lowest-order prediction over a wide range of interchain couplings.

preprint2020arXiv

Real-time simulation of flux qubits used for quantum annealing

The real-time flux dynamics of up to three superconducting quantum interference devices (SQUIDs) are studied by numerically solving the time-dependent Schrödinger equation. The numerical results are used to scrutinize the mapping of the flux degrees of freedom onto two-level systems (the qubits) as well as the performance of the intermediate SQUID as a tunable coupling element. It is shown that the qubit representation yields a good description of the flux dynamics during quantum annealing and the presence of the tunable coupling element does not have negative effects on the overall performance. Additionally, data obtained from a simulation of the dynamics of two-level systems during quantum annealing are compared to experimental data produced by the D-Wave 2000Q quantum annealer. The effects of finite temperature are incorporated in the simulation by coupling the qubit system to a bath of two-level systems. It is shown that an environment modeled as non-interacting two-level systems coupled to the qubits can produce data which matches the experimental data much better than the simulation data of the qubits without coupling to an environment and better than data obtained from a simulation of an environment modeled as interacting two-level systems coupling to the qubits.

preprint2019arXiv

Finite-Difference Time-Domain simulations of transmission microscopy enable a better interpretation of 3D nerve fiber architectures in the brain

In many laboratories, conventional bright-field transmission microscopes are available to study the structure and organization principles of fibrous tissue samples, but they usually provide only 2D information. To access the third (out-of-plane) dimension, more advanced techniques are employed. An example is 3D Polarized Light Imaging (3D-PLI), which measures the birefringence of histological brain sections to derive the spatial nerve fiber orientations. Here, we show how light scattering in transmission microscopy measurements can be leveraged to gain 3D structural information about fibrous tissue samples like brain tissue. For this purpose, we developed a simulation framework using finite-difference time-domain (FDTD) simulations and high performance computing, which can easily be adapted to other microscopy techniques and tissue types with comparable fibrous structures (e.g., muscle fibers, collagen, or artificial fibers). As conventional bright-field transmission microscopy provides usually only 2D information about tissue structures, a three-dimensional reconstruction of fibers across several sections is difficult. By combining our simulations with experimental studies, we show that the polarization-independent transmitted light intensity (transmittance) contains 3D information: We demonstrate in several experimental studies on brain sections from different species (rodent, monkey, human) that the transmittance decreases significantly (by more than 50%) with the increasing out-of-plane angle of the nerve fibers. Our FDTD simulations show that this decrease is mainly caused by polarization-independent light scattering in combination with the finite numerical aperture of the imaging system. This allows to use standard transmission microscopy techniques to obtain 3D information about the fiber inclination and to detect steep fibers, without need for additional measurements.

preprint2019arXiv

Supplementary information for "Quantum supremacy using a programmable superconducting processor"

This is an updated version of supplementary information to accompany "Quantum supremacy using a programmable superconducting processor", an article published in the October 24, 2019 issue of Nature. The main article is freely available at https://www.nature.com/articles/s41586-019-1666-5. Summary of changes since arXiv:1910.11333v1 (submitted 23 Oct 2019): added URL for qFlex source code; added Erratum section; added Figure S41 comparing statistical and total uncertainty for log and linear XEB; new References [1,65]; miscellaneous updates for clarity and style consistency; miscellaneous typographical and formatting corrections.