BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Tilo Wettig

Tilo Wettig appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-thmath-phmath.MPhep-phphysics.comp-phDistributed, Parallel, and Cluster Computingcond-mat.stat-mechmath.NANumerical Analysis

Catalog footprint

What is connected

37works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

37 published item(s)

preprint2022arXiv

Approximation formula for complex spacing ratios in the Ginibre ensemble

Recently, Sá, Ribeiro and Prosen introduced complex spacing ratios to analyze eigenvalue correlations in non-Hermitian systems. At present there are no analytical results for the probability distribution of these ratios in the limit of large system size. We derive an approximation formula for the Ginibre universality class of random matrix theory which converges exponentially fast to the limit of infinite matrix size. We also give results for moments of the distribution in this limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lattice QCD on a novel vector architecture

The SX-Aurora TSUBASA PCIe accelerator card is the newest model of NEC's SX architecture family. Its multi-core vector processor features a vector length of 16 kbits and interfaces with up to 48 GB of HBM2 memory in the current models, available since 2018. The compute performance is up to 2.45 TFlop/s peak in double precision, and the memory throughput is up to 1.2 TB/s peak. New models with improved performance characteristics are announced for the near future. In this contribution we discuss key aspects of the SX-Aurora and describe how we enabled the architecture in the Grid Lattice QCD framework.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Low-rank tensor methods for Markov chains with applications to tumor progression models

Continuous-time Markov chains describing interacting processes exhibit a state space that grows exponentially in the number of processes. This state-space explosion renders the computation or storage of the time-marginal distribution, which is defined as the solution of a certain linear system, infeasible using classical methods. We consider Markov chains whose transition rates are separable functions, which allows for an efficient low-rank tensor representation of the operator of this linear system. Typically, the right-hand side also has low-rank structure, and thus we can reduce the cost for computation and storage from exponential to linear. Previously known iterative methods also allow for low-rank approximations of the solution but are unable to guarantee that its entries sum up to one as required for a probability distribution. We derive a convergent iterative method using low-rank formats satisfying this condition. We also perform numerical experiments illustrating that the marginal distribution is well approximated with low rank.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Induced QCD I: Theory

We explore an alternative discretization of continuum SU(N_c) Yang-Mills theory on a Euclidean spacetime lattice, originally introduced by Budzcies and Zirnbauer. In this discretization the self-interactions of the gauge field are induced by a path integral over N_b auxiliary boson fields, which are coupled linearly to the gauge field. The main progress compared to earlier approaches is that N_b can be as small as N_c. In the present paper we (i) extend the proof that the continuum limit of the new discretization reproduces Yang-Mills theory in two dimensions from gauge group U(N_c) to SU(N_c), (ii) derive refined bounds on N_b for non-integer values, and (iii) perform a perturbative calculation to match the bare parameter of the induced gauge theory to the standard lattice coupling. In follow-up papers we will present numerical evidence in support of the conjecture that the induced gauge theory reproduces Yang-Mills theory also in three and four dimensions, and explore the possibility to integrate out the gauge fields to arrive at a dual formulation of lattice QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiple right-hand-side setup for the DD-αAMG

The setup cost of a modern solver such as DD-αAMG (Wuppertal Multigrid) is a significant contribution to the total time spent on solving the Dirac equation, and in HMC it can even be dominant. We present an improved implementation of this algorithm with modified computation order in the setup procedure. By processing multiple right-hand sides simultaneously we can alleviate many of the performance issues of the default single right-hand-side setup. The main improvements are as follows: By combining multiple right-hand sides the message size for off-chip communication is larger, which leads to better utilization of the network bandwidth. Many matrix-vector products are replaced by matrix-matrix products, leading to better cache reuse. The synchronization overhead inflicted by on-chip parallelization (threading), which is becoming crucial on many-core architectures such as the Intel Xeon Phi, is effectively reduced. In the parts implemented so far, we observe a speedup of roughly 3x compared to the optimized version of the single right-hand-side setup on realistic lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Adaptive algebraic multigrid on SIMD architectures

We present details of our implementation of the Wuppertal adaptive algebraic multigrid code DD-$α$AMG on SIMD architectures, with particular emphasis on the Intel Xeon Phi processor (KNC) used in QPACE 2. As a smoother, the algorithm uses a domain-decomposition-based solver code previously developed for the KNC in Regensburg. We optimized the remaining parts of the multigrid code and conclude that it is a very good target for SIMD architectures. Some of the remaining bottlenecks can be eliminated by vectorizing over multiple test vectors in the setup, which is discussed in the contribution of Daniel Richtmann.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Induced YM theory with auxiliary bosons

We study pure SU(N) lattice gauge theory with a plaquette weight factor given by an inverse determinant which can be written as an integral over auxiliary bosonic fields (modifying a proposal of Budczies and Zirnbauer). We derive conditions for the existence of a continuum limit and its equivalence to Yang-Mills theory. Furthermore, we perturbatively compute the relation between the coupling constants of the `induced' gauge action and the familiar Wilson gauge action using the background-field technique. The perturbative relation agrees well with numerical results for N=2 in three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Itzykson-Zuber-like Integral and Diffusion for Complex Ordinary and Supermatrices

We compute an analogue of the Itzykson-Zuber integral for the case of arbitrary complex matrices. The calculation is done for both ordinary and supermatrices by transferring the Itzykson-Zuber diffusion equation method to the space of arbitrary complex matrices. The integral is of interest for applications in Quantum Chromodynamics and the theory of two-dimensional Quantum Gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Banks-Casher-type relations for complex Dirac spectra

For theories with a sign problem there is no analog of the Banks-Casher relation. This is true in particular for QCD at nonzero quark chemical potential. However, for QCD-like theories without a sign problem the Banks-Casher relation can be extended to the case of complex Dirac eigenvalues. We derive such extensions for the zero-temperature, high-density limits of two-color QCD, QCD at nonzero isospin chemical potential, and adjoint QCD. In all three cases the density of the complex Dirac eigenvalues at the origin is proportional to the BCS gap squared.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Induced QCD with two auxiliary bosonic fields

Following a proposal of Budczies and Zirnbauer, we investigate an alternative lattice discretization of continuum ${\rm SU}(N_c)$ Yang-Mills theory in which the self-interactions of the gauge field are induced by a path integral over $N_b\ge N_c-1$ auxiliary bosonic fields which are coupled linearly to the gauge field. In two dimensions there exists an analytic proof that the new discretization reproduces Yang-Mills theory in its non-perturbative continuum limit. We provide numerical evidence that this is also the case in three and four dimensions and that, after a suitable matching of the free parameters, the results of the induced theory agree with results from the ordinary plaquette action up to lattice artifacts. The new discretization is ideally suited to change the order of integration in the QCD path integral to arrive at formulations in which the gauge fields have been integrated out. The resulting theories might be amenable to methods previously used in the infinite-coupling limit, and we briefly discuss possibilities to arrive at dual representations of lattice QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lattice QCD with Domain Decomposition on Intel Xeon Phi Co-Processors

The gap between the cost of moving data and the cost of computing continues to grow, making it ever harder to design iterative solvers on extreme-scale architectures. This problem can be alleviated by alternative algorithms that reduce the amount of data movement. We investigate this in the context of Lattice Quantum Chromodynamics and implement such an alternative solver algorithm, based on domain decomposition, on Intel Xeon Phi co-processor (KNC) clusters. We demonstrate close-to-linear on-chip scaling to all 60 cores of the KNC. With a mix of single- and half-precision the domain-decomposition method sustains 400-500 Gflop/s per chip. Compared to an optimized KNC implementation of a standard solver [1], our full multi-node domain-decomposition solver strong-scales to more nodes and reduces the time-to-solution by a factor of 5.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sign problem and subsets in one-dimensional QCD

We present a subset method that solves the sign problem for QCD at nonzero quark chemical potential in 0+1 dimensions. The subsets of gauge configurations are constructed using the center symmetry of the SU(3) group. These subsets completely solve the sign problem for up to five flavors. For a larger number of flavors the sign problem slowly reappears, and we propose an extension of the subsets that also solves the sign problem for these cases. The subset method allows for numerical simulations of the model at nonzero chemical potential. We also present some preliminary results on subsets for QCD in two, three, and four dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stressed Cooper pairing in QCD at high isospin density: effective Lagrangian and random matrix theory

We generalize QCD at asymptotically large isospin chemical potential to an arbitrary even number of flavors. We also allow for small quark chemical potentials, which stress the coincident Fermi surfaces of the paired quarks and lead to a sign problem in Monte Carlo simulations. We derive the corresponding low-energy effective theory in both $p$- and $ε$-expansion and quantify the severity of the sign problem. We construct the random matrix theory describing our physical situation and show that it can be mapped to a known random matrix theory at low baryon density so that new insights can be gained without additional calculations. In particular, we explain the Silver Blaze phenomenon at high isospin density. We also introduce stressed singular values of the Dirac operator and relate them to the pionic condensate. Finally we comment on extensions of our work to two-color QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Chiral Condensate of One-Flavor QCD and the Dirac Spectrum at θ=0

In a sector of fixed topological charge, the chiral condensate has a discontinuity given by the Banks-Casher formula also in the case of one-flavor QCD. However, at fixed θ-angle, the chiral condensate remains constant when the quark mass crosses zero. To reconcile these contradictory observations, we have evaluated the spectral density of one-flavor QCD at θ=0. For negative quark mass, it becomes a strongly oscillating function with a period that scales as the inverse space-time volume and an amplitude that increases exponentially with the space-time volume. As we have learned from QCD at nonzero chemical potential, if this is the case, an alternative to the Banks-Casher formula applies, and as we will demonstrate in this talk, for one-flavor QCD this results in a continuous chiral condensate. A special role is played by the topological zero modes which have to be taken into account exactly in order to get a finite chiral condensate in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Banks-Casher-type relation for the BCS gap at high density

We derive a new Banks-Casher-type relation which relates the density of complex Dirac eigenvalues at the origin to the BCS gap of quarks at high density. Our relation is applicable to QCD and QCD-like theories without a sign problem, such as two-color QCD and adjoint QCD with baryon chemical potential, and QCD with isospin chemical potential. It provides us with a method to measure the BCS gap through the Dirac spectrum on the lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

iDataCool: HPC with Hot-Water Cooling and Energy Reuse

iDataCool is an HPC architecture jointly developed by the University of Regensburg and the IBM Research and Development Lab Böblingen. It is based on IBM's iDataPlex platform, whose air-cooling solution was replaced by a custom water-cooling solution that allows for cooling water temperatures of 70C/158F. The system is coupled to an adsorption chiller by InvenSor that operates efficiently at these temperatures. Thus a significant portion of the energy spent on HPC can be recovered in the form of chilled water, which can then be used to cool other parts of the computing center. We describe the architecture of iDataCool and present benchmarks of the cooling performance and the energy (reuse) efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Subset method for one-dimensional QCD

We present a subset method which solves the sign problem for QCD at nonzero quark chemical potential in 0+1 dimensions. The subsets gather gauge configurations based on the center symmetry of the SU(3) group. We show that the sign problem is solved for one to five quark flavors and that it slowly reappears for a larger number of flavors. We formulate an extension of the center subsets that solves the sign problem for a larger number of flavors as well. We also derive some new analytical results for this toy model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Singular values of the Dirac operator at nonzero density

At nonzero density the eigenvalues of the Dirac operator move into the complex plane, while its singular values remain real and nonnegative. In QCD-like theories, the singular-value spectrum carries information on the diquark (or pionic) condensate. We have constructed low-energy effective theories in different density regimes and derived a number of exact results for the Dirac singular values, including Banks-Casher-type relations for the diquark (or pionic) condensate, Smilga-Stern-type relations for the slope of the singular-value density, and Leutwyler-Smilga-type sum rules for the inverse singular values. We also present a rigorous index theorem for non-Hermitian Dirac operators.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Wigner surmise for mixed symmetry classes in random matrix theory

We consider the nearest-neighbor spacing distributions of mixed random matrix ensembles interpolating between different symmetry classes, or between integrable and non-integrable systems. We derive analytical formulas for the spacing distributions of 2x2 or 4x4 matrices and show numerically that they provide very good approximations for those of random matrices with large dimension. This generalizes the Wigner surmise, which is valid for pure ensembles that are recovered as limits of the mixed ensembles. We show how the coupling parameters of small and large matrices must be matched depending on the local eigenvalue density.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Exact results for two-color QCD at low and high density

We discuss a random matrix theory that was originally constructed to describe two-color QCD at low density in the phase with a nonzero chiral condensate. With a particular choice of a parameter, the same random matrix theory also describes the high-density phase of two-color QCD. In this phase a BCS superfluid of diquark pairs is formed, and the pattern of chiral symmetry breaking is very different from that at low density. Analytical results for the spectral density obtained from this random matrix theory allow for the extraction of the BCS gap from lattice data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Geometry dependence of RMT-based methods to extract the low-energy constants Sigma and F

The lowest-order low-energy constants $Σ$ and $F$ of chiral pertubation theory can be extracted from lattice data using methods based on the equivalence of random matrix theory (RMT) and QCD in the epsilon regime. We discuss how the choice of the lattice geometry affects such methods. In particular, we show how to minimize systematic deviations from RMT by an optimal choice of the lattice geometry in the case of two light quark flavors. We illustrate our findings by determining $Σ$ and $F$ from lattice configurations with two dynamical overlap fermions generated by JLQCD, using two different lattice geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Singular values of the Dirac operator in dense QCD-like theories

We study the singular values of the Dirac operator in dense QCD-like theories at zero temperature. The Dirac singular values are real and nonnegative at any nonzero quark density. The scale of their spectrum is set by the diquark condensate, in contrast to the complex Dirac eigenvalues whose scale is set by the chiral condensate at low density and by the BCS gap at high density. We identify three different low-energy effective theories with diquark sources applicable at low, intermediate, and high density, together with their overlapping domains of validity. We derive a number of exact formulas for the Dirac singular values, including Banks-Casher-type relations for the diquark condensate, Smilga-Stern-type relations for the slope of the singular value density, and Leutwyler-Smilga-type sum rules for the inverse singular values. We construct random matrix theories and determine the form of the microscopic spectral correlation functions of the singular values for all nonzero quark densities. We also derive a rigorous index theorem for non-Hermitian Dirac operators. Our results can in principle be tested in lattice simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The QCD sign problem and dynamical simulations of random matrices

At nonzero quark chemical potential dynamical lattice simulations of QCD are hindered by the sign problem caused by the complex fermion determinant. The severity of the sign problem can be assessed by the average phase of the fermion determinant. In an earlier paper we derived a formula for the microscopic limit of the average phase for general topology using chiral random matrix theory. In the current paper we present an alternative derivation of the same quantity, leading to a simpler expression which is also calculable for finite-sized matrices, away from the microscopic limit. We explicitly prove the equivalence of the old and new results in the microscopic limit. The results for finite-sized matrices illustrate the convergence towards the microscopic limit. We compare the analytical results with dynamical random matrix simulations, where various reweighting methods are used to circumvent the sign problem. We discuss the pros and cons of these reweighting methods.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Chiral random matrix theory for two-color QCD at high density

We identify a non-Hermitian chiral random matrix theory that corresponds to two-color QCD at high density. We show that the partition function of the random matrix theory coincides with the partition function of the finite-volume effective theory at high density, and that the Leutwyler-Smilga-type spectral sum rules of the random matrix theory are identical to those derived from the effective theory. The microscopic Dirac spectrum of the theory is governed by the BCS gap, rather than the conventional chiral condensate. We also show that with a different choice of a parameter the random matrix theory yields the effective partition function at low density.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Random matrix analysis of the QCD sign problem

The severity of the sign problem in lattice QCD at nonzero baryon density is measured by the average phase of the fermion determinant. Motivated by the equivalence of chiral random matrix theory and QCD to leading order in the epsilon regime, we compute the phase of the fermion determinant for general topology in random matrix theory as a function of the quark chemical potential and the quark mass. We find that the sign problem becomes milder with increasing topological charge. The analytic predictions are verified by detailed numerical random matrix simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Short-recurrence Krylov subspace methods for the overlap Dirac operator at nonzero chemical potential

The overlap operator in lattice QCD requires the computation of the sign function of a matrix, which is non-Hermitian in the presence of a quark chemical potential. In previous work we introduced an Arnoldi-based Krylov subspace approximation, which uses long recurrences. Even after the deflation of critical eigenvalues, the low efficiency of the method restricts its application to small lattices. Here we propose new short-recurrence methods which strongly enhance the efficiency of the computational method. Using rational approximations to the sign function we introduce two variants, based on the restarted Arnoldi process and on the two-sided Lanczos method, respectively, which become very efficient when combined with multishift solvers. Alternatively, in the variant based on the two-sided Lanczos method the sign function can be evaluated directly. We present numerical results which compare the efficiencies of a restarted Arnoldi-based method and the direct two-sided Lanczos approximation for various lattice sizes. We also show that our new methods gain substantially when combined with deflation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The epsilon expansion at next-to-next-to-leading order with small imaginary chemical potential

We discuss chiral perturbation theory for two and three quark flavors in the epsilon expansion at next-to-next-to-leading order (NNLO) including a small imaginary chemical potential. We calculate finite-volume corrections to the low-energy constants $Σ$ and $F$ and determine the non-universal modifications of the theory, i.e., modifications that cannot be mapped to random matrix theory (RMT). In the special case of two quark flavors in an asymmetric box we discuss how to minimize the finite-volume corrections and non-universal modifications by an optimal choice of the lattice geometry. Furthermore we provide a detailed calculation of a special version of the massless sunset diagram at finite volume.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Chiral Lagrangian and spectral sum rules for two-color QCD at high density

We report on our analytical study of two-color QCD with an even number of flavors at high baryon density. Based on the pattern of chiral symmetry breaking induced by BCS-type diquark pairing we construct the low-energy effective Lagrangian for the Nambu-Goldstone bosons. We also identify a new epsilon-regime at high baryon density and derive Leutwyler-Smilga-type spectral sum rules for the complex eigenvalues of the Dirac operator in terms of the fermion gap. Our results can in principle be tested in lattice QCD simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Krylov subspace methods and the sign function: multishifts and deflation in the non-Hermitian case

Rational approximations of the matrix sign function lead to multishift methods. For non-Hermitian matrices long recurrences can cause storage problems, which can be circumvented with restarts. Together with deflation we obtain efficient iterative methods, as we show in numerical experiments for the overlap Dirac operator at non-vanishing quark chemical potential for lattices up to size 10^4.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Possible large-N transitions for complex Wilson loop matrices

It is shown that a very simple multiplicative random complex matrix model generalizes the large-N phase structure found in the unitary case: A perturbative regime is joined to a non-perturbative regime at a point where the smoothness of some quantities breaks down. A generic complex Wilson loop matrix in a field theory admitting a 't Hooft planar limit could display a phase transition in that limit as nonlinear effects become dominating over linear ones.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

U(1) staggered Dirac operator and random matrix

We investigate the spectrum of the staggered Dirac operator in 4d quenched U(1) lattice gauge theory and its relationship to random matrix theory. In the confined as well as in the Coulomb phase the nearest-neighbor spacing distribution of the unfolded eigenvalues is well described by the chiral unitary ensemble. The same is true for the distribution of the smallest eigenvalue and the microscopic spectral density in the confined phase. The physical origin of the chiral condensate in this phase deserves further study.

0 citations0 reviewsNo code linkedOpen access