Researcher profile

Karl Jansen

Karl Jansen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

27 published item(s)

preprint2026arXiv

Computing Statistical Properties of Velocity Fields on Current Quantum Hardware

Quantum algorithms are gaining attention in Computational Fluid Dynamics (CFD) for their favorable scaling, as encoding physical fields into quantum probability amplitudes enables representation of two to the power of n spatial points with only n qubits. A key challenge in Quantum CFD is the efficient readout of simulation results, a topic that has received limited attention in literature. This work presents methods to extract statistical properties of spatial velocity fields, such as central moments and structure functions, directly from parameterized ansatz circuits, avoiding full quantum state tomography. As a proof of concept, we implement our approach for 1D velocity fields, encoding 16 spatial points with 4 qubits, and analyze both a sine wave signal and four snapshots from Burgers' equation evolution. Using Qedma's error mitigation software QESEM, we demonstrate that such computations achieve high accuracy on current quantum devices, specifically IBMQ's Heron2 system ibm_fez.

preprint2025arXiv

Processing through encoding: Quantum circuit approaches for point-wise multiplication and convolution

This paper introduces quantum circuit methodologies for pointwise multiplication and convolution of complex functions, conceptualized as "processing through encoding". Leveraging known techniques, we describe an approach where multiple complex functions are encoded onto auxiliary qubits. Applying the proposed scheme for two functions $f$ and $g$, their pointwise product $f(x)g(x)$ is shown to naturally form as the coefficients of part of the resulting quantum state. Adhering to the convolution theorem, we then demonstrate how the convolution $f*g$ can be constructed. Similarly to related work, this involves the encoding of the Fourier coefficients $\mathcal{F}[f]$ and $\mathcal{F}[g]$, which facilitates their pointwise multiplication, followed by the inverse Quantum Fourier Transform. We discuss the simulation of these techniques, their integration into an extended \verb|quantumaudio| package for audio signal processing, and present initial experimental validations. This work offers a promising avenue for quantum signal processing, with potential applications in areas such as quantum-enhanced audio manipulation and synthesis.

preprint2022arXiv

Digitising SU(2) Gauge Fields and the Freezing Transition

Efficient discretisations of gauge groups are crucial with the long term perspective of using tensor networks or quantum computers for lattice gauge theory simulations. For any Lie group other than U$(1)$, however, there is no class of asymptotically dense discrete subgroups. Therefore, discretisations limited to subgroups are bound to lead to a freezing of Monte Carlo simulations at weak couplings, necessitating alternative partitionings without a group structure. In this work we provide a comprehensive analysis of this freezing for all discrete subgroups of SU$(2)$ and different classes of asymptotically dense subsets. We find that an appropriate choice of the subset allows unfrozen simulations for arbitrary couplings, though one has to be careful with varying weights of unevenly distributed points. A generalised version of the Fibonacci spiral appears to be particularly efficient and close to optimal.

preprint2022arXiv

Lattice QCD Study of Transverse-Momentum Dependent Soft Function

In this work, we perform a lattice QCD study of the intrinsic, rapidity-independent soft function within the framework of large momentum effective theory. The computation is carried out using a gauge ensemble of $N_f=2+1+1$ clover-improved twisted mass fermion. After applying an appropriate renormalization procedure and the removal of significant higher-twist contamination, we obtain the intrinsic soft function that is comparable to the one-loop perturbative result at large external momentum. The determination of the nonperturbative soft function from first principles is crucial to sharpen our understanding of the processes with small transverse momentum such as the Drell-Yan production and the semi-inclusive deep inelastic scattering. Additionally, we calculate the Collins-Soper evolution kernel using the quasi-transverse-momentum-dependent wave function as input.

preprint2022arXiv

Measurement Error Mitigation in Quantum Computers Through Classical Bit-Flip Correction

We develop a classical bit-flip correction method to mitigate measurement errors on quantum computers. This method can be applied to any operator, any number of qubits, and any realistic bit-flip probability. We first demonstrate the successful performance of this method by correcting the noisy measurements of the ground-state energy of the longitudinal Ising model. We then generalize our results to arbitrary operators and test our method both numerically and experimentally on IBM quantum hardware. As a result, our correction method reduces the measurement error on the quantum hardware by up to one order of magnitude. We finally discuss how to pre-process the method and extend it to other errors sources beyond measurement errors. For local Hamiltonians, the overhead costs are polynomial in the number of qubits, even if multi-qubit correlations are included.

preprint2022arXiv

New Directions in Quantum Music: concepts for a quantum keyboard and the sound of the Ising model

We explore ideas for generating sounds and eventually music by using quantum devices in the NISQ era using quantum circuits. In particular, we first consider a concept for a "qeyboard", i.e. a quantum keyboard, where the real-time behaviour of expectation values using a time evolving quantum circuit can be associated to sound features like intensity, frequency and tone. Then, we examine how these properties can be extracted from physical quantum systems, taking the Ising model as an example. This can be realized by measuring physical quantities of the quantum states of the system, e.g. the energies and the magnetization obtained via variational quantum simulation techniques.

preprint2022arXiv

Phase structure of the CP(1) model in the presence of a topological $θ$-term

We numerically study the phase structure of the CP(1) model in the presence of a topological $θ$-term, a regime afflicted by the sign problem for conventional lattice Monte Carlo simulations. Using a bond-weighted tensor renormalization group method, we compute the free energy for inverse couplings ranging from $0\leq β\leq 1.1$ and find a CP-violating, first-order phase transition at $θ=π$. In contrast to previous findings, our numerical results provide no evidence for a critical coupling $β_c<1.1$ above which a second-order phase transition emerges at $θ=π$ and/or the first-order transition line bifurcates at $θ\neqπ$. If such a critical coupling exists, as suggested by Haldane&#39;s conjecture, our study indicates that is larger than $β_c>1.1$.

preprint2022arXiv

Pion-pole contribution to HLbL from twisted mass lattice QCD at the physical point

We report on our computation of the pion transition form factor ${\cal F}_{P\rightarrow γ^*γ^*}$ from twisted mass lattice QCD in order to determine the numerically dominant light pseudoscalar pole contribution in the hadronic light-by-light scattering contribution to the anomalous magnetic moment of the muon $a_μ=(g-2)_μ$. The pion transition form factor is computed directly at the physical point. We present first results for our estimate of the pion-pole contribution with kinematic setup for the pion at rest.

preprint2022arXiv

Proton generalized parton distributions from lattice QCD

Momentum and spatial distributions of quarks and gluons inside hadrons are typically encoded in the so-called generalized parton distributions (GPDs). GPDs are multi-dimensional quantities that are very challenging to extract, both experimentally and within lattice QCD. We present the first lattice results on the $x$-dependence of isovector unpolarized, helicity and transversity GPDs of the proton, obtained from lattice QCD using an ensemble of $N_f=2+1+1$ maximally twisted mass fermions, with pion mass $M_π=260$ MeV and lattice spacing $a\simeq 0.093$ fm. Our calculations use the quasi-distribution formalism and the final distributions are presented in the MS-bar scheme at a renormalization scale of 2 GeV.

preprint2022arXiv

Snowmass White Paper: Quantum Computing Systems and Software for High-energy Physics Research

Quantum computing offers a new paradigm for advancing high-energy physics research by enabling novel methods for representing and reasoning about fundamental quantum mechanical phenomena. Realizing these ideals will require the development of novel computational tools for modeling and simulation, detection and classification, data analysis, and forecasting of high-energy physics (HEP) experiments. While the emerging hardware, software, and applications of quantum computing are exciting opportunities, significant gaps remain in integrating such techniques into the HEP community research programs. Here we identify both the challenges and opportunities for developing quantum computing systems and software to advance HEP discovery science. We describe opportunities for the focused development of algorithms, applications, software, hardware, and infrastructure to support both practical and theoretical applications of quantum computing to HEP problems within the next 10 years.

preprint2022arXiv

Strategies for the Determination of the Running Coupling of $(2+1)$-dimensional QED with Quantum Computing

We propose to utilize NISQ-era quantum devices to compute short distance quantities in $(2+1)$-dimensional QED and to combine them with large volume Monte Carlo simulations and perturbation theory. On the quantum computing side, we perform a calculation of the mass gap in the small and intermediate regime, demonstrating, in the latter case, that it can be resolved reliably. The so obtained mass gap can be used to match corresponding results from Monte Carlo simulations, which can be used eventually to set the physical scale. In this paper we provide the setup for the quantum computation and show results for the mass gap and the plaquette expectation value. In addition, we discuss some ideas that can be applied to the computation of the running coupling. Since the theory is asymptotically free, it would serve as a training ground for future studies of QCD in $(3+1)$-dimensions on quantum computers.

preprint2022arXiv

Transversity GPDs of the proton from lattice QCD

We present the first calculation of the $x$-dependence of the isovector transversity generalized parton distributions (GPDs) for the proton within lattice QCD. We compute the matrix elements with non-local operators containing a Wilson line. The calculation implements the Breit symmetric frame. The proton momenta are chosen as $0.83,\,1.25,\,1.67$ GeV, and the values of the momentum transfer squared are $0.69,\,1.02$ GeV$^2$. These combinations include cases with zero and nonzero skewness. The calculation is performed using one ensemble of two degenerate-mass light, a strange and a charm quark of maximally twisted mass fermions with a clover term. The lattice results are renormalized non-perturbatively and finally matched to the light-cone GPDs using one-loop perturbation theory within the framework of large momentum effective theory. The final GPDs are given in the $\overline{\rm MS}$ scheme at a scale of 2 GeV. In addition to the individual GPDs, we form the combination of the transversity GPDs that is related to the transverse spin structure of the proton. Finally, we extract the lowest two moments of GPDs and draw a number of important qualitative conclusions.

preprint2022arXiv

Twisted mass gauge ensembles at physical values of the light, strange and charm quark masses

Lattice QCD simulations directly at physical masses of dynamical light, strange and charm quarks are highly desirable especially to remove systematic errors due to chiral extrapolations. However such simulations are still challenging. We discuss the adaption of efficient algorithms, like multi-grid methods or higher order integrators, within the molecular dynamic steps of the Hybrid Monte Carlo algorithm, that are enabling simulations of a new set of gauge ensembles by the Extended Twisted Mass collaboration (ETMC). We present the status of the on-going ETMC simulation effort that aims to enabling studies of finite size and discretization effects. We work within the twisted mass discretization which is free of odd-discretization effects at maximal twist and present our tuning procedure.

preprint2022arXiv

Using classical bit-flip correction for error mitigation including 2-qubit correlations

We present an error mitigation scheme which corrects readout errors on Noisy Intermediate-Scale Quantum (NISQ) computers [1,2]. After a short review of applying the method to one qubit, we proceed to discuss the case when correlations between different qubits occur. We demonstrate how the readout error can be mitigated in this case. By performing experiments on IBMQ hardware, we show that such correlations do not have a strong effect on the results, justifying to neglect them.

preprint2022arXiv

Zeta-regularized Lattice Field Theory with Lorentzian background metrics

Lattice field theory is a very powerful tool to study Feynman&#39;s path integral non-perturbatively. However, it usually requires Euclidean background metrics to be well-defined. On the other hand, a recently developed regularization scheme based on Fourier integral operator $ζ$-functions can treat Feynman&#39;s path integral non-pertubatively in Lorentzian background metrics. In this article, we formally $ζ$-regularize lattice theories with Lorentzian backgrounds and identify conditions for the Fourier integral operator $ζ$-function regularization to be applicable. Furthermore, we show that the classical limit of the $ζ$-regularized theory is independent of the regularization. Finally, we consider the harmonic oscillator as an explicit example. We discuss multiple options for the regularization and analytically show that they all reproduce the correct ground state energy on the lattice and in the continuum limit. Additionally, we solve the harmonic oscillator on the lattice in Minkowski background numerically.

preprint2021arXiv

A resource efficient approach for quantum and classical simulations of gauge theories in particle physics

Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present limitations of MCMC techniques may be overcome by Hamiltonian-based simulations on classical or quantum devices, which further provide the potential to address questions that lay beyond the capabilities of the current approaches. However, for continuous gauge groups, Hamiltonian-based formulations involve infinite-dimensional gauge degrees of freedom that can solely be handled by truncation. Current truncation schemes require dramatically increasing computational resources at small values of the bare couplings, where magnetic field effects become important. Such limitation precludes one from `taking the continuous limit&#39; while working with finite resources. To overcome this limitation, we provide a resource-efficient protocol to simulate LGTs with continuous gauge groups in the Hamiltonian formulation. Our new method allows for calculations at arbitrary values of the bare coupling and lattice spacing. The approach consists of the combination of a Hilbert space truncation with a regularization of the gauge group, which permits an efficient description of the magnetically-dominated regime. We focus here on Abelian gauge theories and use $2+1$ dimensional quantum electrodynamics as a benchmark example to demonstrate this efficient framework to achieve the continuum limit in LGTs. This possibility is a key requirement to make quantitative predictions at the field theory level and offers the long-term perspective to utilise quantum simulations to compute physically meaningful quantities in regimes that are precluded to quantum Monte Carlo.

preprint2021arXiv

Estimation of Thermodynamic Observables in Lattice Field Theories with Deep Generative Models

In this work, we demonstrate that applying deep generative machine learning models for lattice field theory is a promising route for solving problems where Markov Chain Monte Carlo (MCMC) methods are problematic. More specifically, we show that generative models can be used to estimate the absolute value of the free energy, which is in contrast to existing MCMC-based methods which are limited to only estimate free energy differences. We demonstrate the effectiveness of the proposed method for two-dimensional $ϕ^4$ theory and compare it to MCMC-based methods in detailed numerical experiments.

preprint2021arXiv

Ruling Out the Massless Up-Quark Solution to the Strong CP Problem by Computing the Topological Mass Contribution with Lattice QCD

The infamous strong CP problem in particle physics can in principle be solved by a massless up quark. In particular, it was hypothesized that topological effects could substantially contribute to the observed nonzero up-quark mass without reintroducing CP violation. Alternatively to previous work using fits to chiral perturbation theory, in this Letter, we bound the strength of the topological mass contribution with direct lattice QCD simulations, by computing the dependence of the pion mass on the dynamical strange-quark mass. We find that the size of the topological mass contribution is inconsistent with the massless up-quark solution to the strong CP problem.

preprint2021arXiv

Towards Quantum Simulations in Particle Physics and Beyond on Noisy Intermediate-Scale Quantum Devices

We review two algorithmic advances that bring us closer to reliable quantum simulations of model systems in high energy physics and beyond on noisy intermediate-scale quantum (NISQ) devices. The first method is the dimensional expressivity analysis of quantum circuits, which allows for constructing minimal but maximally expressive quantum circuits. The second method is an efficient mitigation of readout errors on quantum devices. Both methods can lead to significant improvements in quantum simulations, e.g., when variational quantum eigensolvers are used.

preprint2020arXiv

Avoiding the sign-problem in lattice field theory

In lattice field theory, the interactions of elementary particles can be computed via high-dimensional integrals. Markov-chain Monte Carlo (MCMC) methods based on importance sampling are normally efficient to solve most of these integrals. But these methods give large errors for oscillatory integrands, exhibiting the so-called sign-problem. We developed new quadrature rules using the symmetry of the considered systems to avoid the sign-problem in physical one-dimensional models for the resulting high-dimensional integrals. This article gives a short introduction to integrals used in lattice QCD where the interactions of gluon and quark elementary particles are investigated, explains the alternative integration methods we developed and shows results of applying them to models with one physical dimension. The new quadrature rules avoid the sign-problem and can therefore be used to perform simulations at until now not reachable regions in parameter space, where the MCMC errors are too big for affordable sample sizes. However, it is still a challenge to develop these techniques further for applications with physical higher-dimensional systems.

preprint2020arXiv

Comparison of topological charge definitions in Lattice QCD

In this paper, we show a comparison of different definitions of the topological charge on the lattice. We concentrate on one small-volume ensemble with 2 flavours of dynamical, maximally twisted mass fermions and use three more ensembles to analyze the approach to the continuum limit. We investigate several fermionic and gluonic definitions. The former include the index of the overlap Dirac operator, the spectral flow of the Wilson--Dirac operator and the spectral projectors. For the latter, we take into account different discretizations of the topological charge operator and various smoothing schemes to filter out ultraviolet fluctuations: the gradient flow, stout smearing, APE smearing, HYP smearing and cooling. We show that it is possible to perturbatively match different smoothing schemes and provide a well-defined smoothing scale. We relate the smoothing parameters for cooling, stout and APE smearing to the gradient flow time $τ$. In the case of hypercubic smearing the matching is performed numerically. We investigate which conditions have to be met to obtain a valid definition of the topological charge and susceptibility and we argue that all valid definitions are highly correlated and allow good control over topology on the lattice.

preprint2020arXiv

From spin chains to real-time thermal field theory using tensor networks

One of the most interesting directions in theoretical high-energy and condensed-matter physics is understanding dynamical properties of collective states of quantum field theories. The most elementary tool in this quest is retarded equilibrium correlators governing the linear response theory. In this article we examine tensor networks as a way of determining them in a fully ab initio way in a class of (1+1)-dimensional quantum field theories arising as infrared descriptions of quantum Ising chains.We show that, complemented with signal analysis using the Prony method, tensor network calculations for intermediate times provide a powerful way to explore the structure of singularities of the correlator in the complex frequency plane and to make predictions about the thermal response to perturbations in a class of nonintegrable interacting quantum field theories.

preprint2020arXiv

Improvement, generalization, and scheme conversion of Wilson-line operators on the lattice in the auxiliary field approach

Nonlocal quark bilinear operators connected by link paths are used for studying parton distribution functions (PDFs) and transverse momentum-dependent PDFs of hadrons using lattice QCD. The nonlocality makes it difficult to understand the renormalization and improvement of these operators using standard methods. In previous work, we showed that by introducing an auxiliary field on the lattice, one can understand an on-axis Wilson-line operator as the product of two local operators in an extended theory. In this paper, we provide details about the calculation in perturbation theory of the factor for conversion from our lattice-suitable renormalization scheme to the MS-bar scheme. Extending our work, we study Symanzik improvement of the extended theory to understand the pattern of discretization effects linear in the lattice spacing, $a$, which are present even if the lattice fermion action exactly preserves chiral symmetry. This provides a prospect for an eventual $O(a)$ improvement of lattice calculations of PDFs. We also generalize our approach to apply to Wilson lines along lattice diagonals and to piecewise-straight link paths.

preprint2020arXiv

Multilevel Monte Carlo for quantum mechanics on a lattice

Monte Carlo simulations of quantum field theories on a lattice become increasingly expensive as the continuum limit is approached since the cost per independent sample grows with a high power of the inverse lattice spacing. Simulations on fine lattices suffer from critical slowdown, the rapid growth of autocorrelations in the Markov chain. This causes a strong increase in the number of lattice configurations that have to be generated to obtain statistically significant results. This paper discusses hierarchical sampling methods to tame the growth in autocorrelations. Combined with multilevel variance reduction, this significantly reduces the computational cost of simulations for given tolerances $ε_{\text{disc}}$ on the discretisation error and $ε_{\text{stat}}$ on the statistical error. For observables with lattice errors of order $α$ and integrated autocorrelation times that grow like $τ_{\mathrm{int}}\propto a^{-z}$, multilevel Monte Carlo (MLMC) reduces the cost from $\mathcal{O}(ε_{\text{stat}}^{-2}ε_{\text{disc}}^{-(1+z)/α})$ to $\mathcal{O}(ε_{\text{stat}}^{-2}\vert\log ε_{\text{disc}} \vert^2+ε_{\text{disc}}^{-1/α})$ or $\mathcal{O}(ε_{\text{stat}}^{-2}+ε_{\text{disc}}^{-1/α})$. Higher gains are expected for simulations of quantum field theories in $D$ dimensions. The efficiency of the approach is demonstrated on two model systems, including a topological oscillator that is badly affected by critical slowdown from topological charge freezing. On fine lattices, the new methods are orders of magnitude faster than standard Hybrid Monte Carlo sampling. For high resolutions, MLMC can be used to accelerate even the cluster algorithm for the topological oscillator. Performance is further improved through perturbative matching which guarantees efficient coupling of theories on the multilevel hierarchy.

preprint2020arXiv

Parton distribution functions of $Δ^+$ on the lattice

We perform a first calculation for the unpolarized parton distribution function of the $Δ^+$ baryon using lattice QCD simulations within the framework of Large Momentum Effective Theory. Two ensembles of $N_f=2+1+1$ twisted mass fermions are utilized with a pion mass of 270 MeV and 360 MeV, respectively. The baryon, which is treated as a stable single-particle state, is boosted with momentum $P_3$ with values $\{0.42,0.83,1.25\}$ GeV, and we utilize momentum smearing to improve the signal. The unpolarized parton distribution function of $Δ^+$ is obtained using a non-perturbative renormalization and a one-loop formula for the matching, with encouraging precision. In particular, we compute the $\overline{d}(x)-\overline{u}(x)$ asymmetry and compare it with the same quantity in the nucleon, in a first attempt towards resolving the physical mechanism responsible for generating such asymmetry.

preprint2020arXiv

Topological vacuum structure of the Schwinger model with matrix product states

We numerically study the single-flavor Schwinger model with a topological $θ$-term, which is practically inaccessible by standard lattice Monte Carlo simulations due to the sign problem. By using numerical methods based on tensor networks, especially the one-dimensional matrix product states, we explore the non-trivial $θ$-dependence of several lattice and continuum quantities in the Hamiltonian formulation. In particular, we compute the ground-state energy, the electric field, the chiral fermion condensate, and the topological vacuum susceptibility for positive, zero, and even negative fermion mass. In the chiral limit, we demonstrate that the continuum model becomes independent of the vacuum angle $θ$, thus respecting CP invariance, while lattice artifacts still depend on $θ$. We also confirm that negative masses can be mapped to positive masses by shifting $θ\rightarrow θ+π$ due to the axial anomaly in the continuum, while lattice artifacts non-trivially distort this mapping. This mass regime is particularly interesting for the (3+1)-dimensional QCD analog of the Schwinger model, the sign problem of which requires the development and testing of new numerical techniques beyond the conventional Monte Carlo approach.

preprint2020arXiv

Unpolarized and helicity generalized parton distributions of the proton within lattice QCD

We present the first calculation of the $x$-dependence of the proton generalized parton distributions (GPDs) within lattice QCD. Results are obtained for the isovector unpolarized and helicity GPDs. We compute the appropriate matrix elements of fast-moving protons coupled to non-local operators containing a Wilson line. We present results for proton momenta $0.83,\,1.25,\,1.67$ GeV, and momentum transfer squared $0.69,\,1.38$ GeV$^2$. These combinations include cases with zero and nonzero skewness. The calculation is performed using one ensemble of two degenerate mass light, a strange and a charm quark of maximally twisted mass fermions with a clover term. The lattice results are matched to the light-cone GPDs using one-loop perturbation theory within the framework of large momentum effective theory. The final GPDs are given in the $\overline{\rm MS}$ scheme at a scale of 2 GeV.