Analytic continuation is a central step in the simulation of finite-temperature field theories in which numerically obtained Matsubara data are continued to the real frequency axis for a physical interpretation. Numerical analytic continuation is considered to be an ill-posed problem where uncertainties on the Matsubara axis are amplified exponentially. Here, we present a systematic and controlled procedure that approximates any Matsubara function by a minimal pole representation to within a predefined precision. We then show a systematic convergence to the exact spectral function on the real axis as a function of our control parameter for a range of physically relevant setups. Our methodology is robust to noise and paves the way towards reliable analytic continuation in many-body theory and, by providing access to the analytic structure of the functions, a direct theoretical interpretation of physical properties.
more »
« less
Analytic continuation of noisy data using Adams Bashforth residual neural network
We propose a data-driven learning framework for the analytic continuation problem in numerical quantum many-body physics. Designing an accurate and efficient framework for the analytic continuation of imaginary time using computational data is a grand challenge that has hindered meaningful links with experimental data. The standard Maximum Entropy (MaxEnt)-based method is limited by the quality of the computational data and the availability of prior information. Also, the MaxEnt is not able to solve the inversion problem under high level of noise in the data. Here we introduce a novel learning model for the analytic continuation problem using a Adams-Bashforth residual neural network (AB-ResNet). The advantage of this deep learning network is that it is model independent and, therefore, does not require prior information concerning the quantity of interest given by the spectral function. More importantly, the ResNet-based model achieves higher accuracy than MaxEnt for data with higher level of noise. Finally, numerical examples show that the developed AB-ResNet is able to recover the spectral function with accuracy comparable to MaxEnt where the noise level is relatively small.
more »
« less
- Award ID(s):
- 1720222
- PAR ID:
- 10323046
- Date Published:
- Journal Name:
- Discrete & Continuous Dynamical Systems - S
- Volume:
- 15
- Issue:
- 4
- ISSN:
- 1937-1632
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We introduce a Julia implementation of the recently proposed Nevanlinna analytic continuation method. The method is based on Nevanlinna interpolants and inherently preserves the causality of a response function due to its construction. For theoretical calculations without statistical noise, this continuation method is a powerful tool to extract real-frequency information from numerical input data on the Matsubara axis. This method has been applied to first-principles calculations of correlated materials. This paper presents its efficient and full-featured open-source implementation of the method including the Hamburger moment problem and smoothing.more » « less
-
We present the combination of a complex-time tensor-network impurity solver with an analytic continuation scheme based on exponential fitting as an efficient framework for single and multiorbital dynamical mean-field calculations. By performing time evolution along a complex-time contour, the approach balances computational cost with the difficulty of spectral recovery, offering greater flexibility than methods confined to the real or imaginary axis. By complementing the complex-time evolution with an exponential fitting scheme, we faithfully extract real-time information at negligible cost. The resulting method obtains high-resolution spectra at a significantly lower computational cost than real-time evolution, offering a promising tool for studies of strongly correlated materials.more » « less
-
In this work, we propose Mutual Information Maximization Knowledge Distillation (MIMKD). Our method uses a contrastive objective to simultaneously estimate and maximize a lower bound on the mutual information of local and global feature representations between a teacher and a student network. We demonstrate through extensive experiments that this can be used to improve the performance of low capacity models by transferring knowledge from more performant but computationally expensive models. This can be used to produce better models that can be run on devices with low computational resources. Our method is flexible, we can distill knowledge from teachers with arbitrary network architectures to arbitrary student networks. Our empirical results show that MIMKD outperforms competing approaches across a wide range of student-teacher pairs with different capacities, with different architectures, and when student networks are with extremely low capacity. We are able to obtain 74.55% accuracy on CIFAR100 with a ShufflenetV2 from a baseline accuracy of 69.8% by distilling knowledge from ResNet-50. On Imagenet we improve a ResNet-18 network from 68.88% to 70.32% accuracy (1.44%+) using a ResNet-34 teacher network.more » « less
-
We study a class of spectral design problems in which a prior positive semidefinite information matrix is updated by a sum of rank-one matrices constructed from chosen design vectors subject to a bound on their Euclidean norm. The objective of a spectral design problem is any symmetric convex function of the eigenvalues of the updated information matrix. Thisframeworkunifiesclassicaloptimalexperimentaldesigncriteria, includingA-, D-, andE-optimality. Italsoarises in model-basedderivative-freeoptimization, wheresampling directions determine the conditioning and accuracy of regression models. Although the objective is symmetric and convex in the eigenvalues, the optimization problem with design vectors/matrix as decision variables is nonconvex, and optimal solutions of their convex relaxations may not be feasible for the spectral design problem. We use tight eigenvalue relaxations to obtain a convex reformulation, and we apply the Schur–Horn theorem to construct a simple polynomial-time algorithm for solving the spectral design problem. We illustrate the optimal spectral designs computed by our algorithm. Moreover, a small set of numerical experiments shows the potential of spectral designs for derivative-free optimization.more » « less
An official website of the United States government

