Related papers: Lanczos, the transfer matrix, and the signal-to-no…
One proposal to compute parton distributions from first principles is the large momentum effective theory (LaMET), which requires the Fourier transform of matrix elements computed non-perturbatively. Lattice quantum chromodynamics (QCD)…
Two numerical algorithms for the computation of eigenvalues of Dirac operators in lattice gauge theories are described: one is an accelerated conjugate gradient method, the other one a standard Lanczos method. Results obtained by Cullum's…
The runtime for Kernel Partial Least Squares (KPLS) to compute the fit is quadratic in the number of examples. However, the necessity of obtaining sensitivity measures as degrees of freedom for model selection or confidence intervals for…
The one-hole spectral weight for two chains and two dimensional lattices is studied numerically using a new method of analysis of the spectral function within the Lanczos iteration scheme: the Lanczos spectra decoding method. This technique…
A deflated restarted Lanczos algorithm is given for both solving symmetric linear equations and computing eigenvalues and eigenvectors. The restarting limits the storage so that finding eigenvectors is practical. Meanwhile, the deflating…
We introduce a new approach for estimating the number of spikes in a general class of spiked covariance models without directly computing the eigenvalues of the sample covariance matrix. This approach is based on the Lanczos algorithm and…
The variational optimization of high-dimensional neural network models, such as those used in neural quantum states (NQS), presents a significant challenge in machine intelligence. Conventional first-order stochastic methods (e.g., Adam)…
This work considers large-scale Lyapunov matrix equations of the form $AX + XA = \boldsymbol{c}\boldsymbol{c}^T$, where $A$ is a symmetric positive definite matrix and $\boldsymbol{c}$ is a vector. Motivated by the need to solve such…
The time-ordered exponential is defined as the function that solves a system of coupled first-order linear differential equations with generally non-constant coefficients. In spite of being at the heart of much system dynamics, control…
In this work we generalize the concept of modular spread complexity to the cases where the reduced density matrix is non-Hermitian. This notion of complexity and associated Lanczos coefficients contain richer information than the…
The Levenberg-Marquardt (LM) method is commonly used for inverting models used to describe geothermal, groundwater, or oil and gas reservoirs. In previous studies LM parameter updates have been made tractable for highly parameterized…
Lanczos-type algorithms are well known for their inherent instability. They typically breakdown when relevant orthogonal polynomials do not exist. Current approaches to avoiding breakdown rely on jumping over the non-existent polynomials to…
We present a completely unbiased and controlled numerical method to solve quantum impurity problems in d-dimensional lattices. This approach is based on a canonical transformation, of the Lanczos form, where the complete lattice Hamiltonian…
To understand the dynamics of quantum many-body systems, it is essential to study excited eigenstates. While tensor network states have become a standard tool for computing ground states in computational many-body physics, obtaining…
In this paper we compare two methods for finding extremal eigenvalues and eigenvectors: the restarted Lanczos method and momentum accelerated power iterations. The convergence of both methods is based on ratios of Chebyshev polynomials…
A numerical simulation algorithm for lattice QCD is described, in which the short- and long-distance effects of the sea quarks are treated separately. The algorithm can be regarded, to some extent, as an implementation at the quantum level…
An outstanding goal of physics is to find solutions that describe hadrons in the theory of strong interactions, Quantum Chromodynamics (QCD). For this goal, the light-front Hamiltonian formulation of QCD (LFQCD) is a complementary approach…
We study the equilibration times $T_\text{eq}$ of local observables in quantum chaotic systems by considering their auto-correlation functions. Based on the recursion method, we suggest a scheme to estimate $T_\text{eq}$ from the…
In the context of Monte Carlo sampling for lattice models, the complexity of the energy landscape often leads to Markov chains being trapped in local optima, thereby increasing the correlation between samples and reducing sampling…
The eigenvalue spectrum $\rho(\lambda)$ of the Dirac operator is numerically calculated in lattice QCD with 2+1 flavors of dynamical domain-wall fermions. In the high-energy regime, the discretization effects become significant. We subtract…