Efficient low temperature simulations for fermionic reservoirs with the hierarchical equations of motion method: Application to the Anderson impurity model

• URL: http://arxiv.org/abs/2211.04089v2
• Date: Tue, 9 May 2023 11:40:17 GMT
• Title: Efficient low temperature simulations for fermionic reservoirs with the hierarchical equations of motion method: Application to the Anderson impurity model
• Authors: Xiaohan Dan (1 and 2), Meng Xu (3), J. T. Stockburger (3), J. Ankerhold (3), Qiang Shi (1 and 2) ((1) Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Beijing, China, (2) University of Chinese Academy of Sciences, Beijing, China, (3) Institute for Complex Quantum Systems and IQST, Ulm University, Ulm, Germany)
• Abstract summary: In this work, we employ the barycentric representation to approximate the Fermi function and hybridization functions in the frequency domain. The new method, by approxing these functions with optimized rational decomposition, greatly reduces the number of basis functions in decomposing the reservoir correlation functions. We demonstrate the efficiency, accuracy, and long-time stability of the new decomposition scheme by applying it to the Anderson impurity model (AIM) in the low-temperature regime.
• Score: 0.0
• License: http://creativecommons.org/licenses/by/4.0/
• Abstract: The hierarchical equations of motion (HEOM) approach is an accurate method to simulate open system quantum dynamics, which allows for systematic convergence to numerically exact results. To represent the effects of the bath, the reservoir correlation functions are usually decomposed into the summation of multiple exponential terms in the HEOM method. Since the reservoir correlation functions become highly non-Markovian at low temperatures or when the bath has complex band structures, a present challenge is to obtain accurate exponential decompositions that allow efficient simulation with the HEOM. In this work, we employ the barycentric representation to approximate the Fermi function and hybridization functions in the frequency domain. The new method, by approximating these functions with optimized rational decomposition, greatly reduces the number of basis functions in decomposing the reservoir correlation functions, which further allows the HEOM method to be applied to ultra-low temperature and general bath structures. We demonstrate the efficiency, accuracy, and long-time stability of the new decomposition scheme by applying it to the Anderson impurity model (AIM) in the low-temperature regime with the Lorentzian and tight-binding hybridization functions.

Related papers

• Generalized system-bath entanglement theorem for Gaussian environments [10.07019992236203]
  A system-bath entanglement theorem (SBET) with Gaussian environments was established previously in J. Chem. Phys. 152, 034102 ( 2020) In this work, we generalize it to correlation functions. Numerical demonstrations are carried out for the evaluation of the solvation free energy of an electron transfer system.
  arXiv  Detail & Related papers  (2023-12-21T07:11:23Z)
• Stochastic Schr\"odinger equation approach to real-time dynamics of Anderson-Holstein impurities: an open quantum system perspective [3.105656247358225]
  We develop a framework to simulate dynamics of Anderson-Holstein impurities coupled to a continuous fermionic bath. We show that such an SSE treatment provides a middle ground between numerically expansive microscopic simulations and master equations.
  arXiv  Detail & Related papers  (2023-09-16T06:03:54Z)
• An Optimization-based Deep Equilibrium Model for Hyperspectral Image Deconvolution with Convergence Guarantees [71.57324258813675]
  We propose a novel methodology for addressing the hyperspectral image deconvolution problem. A new optimization problem is formulated, leveraging a learnable regularizer in the form of a neural network. The derived iterative solver is then expressed as a fixed-point calculation problem within the Deep Equilibrium framework.
  arXiv  Detail & Related papers  (2023-06-10T08:25:16Z)
• Tree tensor network state approach for solving hierarchical equations of motion [0.0]
  The hierarchical equations of motion (HEOM) method is a numerically exact open quantum system dynamics approach. We show that, the proposed HEOM+TTNS approach yields consistent results with that of the conventional HEOM method. Besides, the simulation with a genuine TTNS is four times faster than a one-dimensional matrix product state decomposition scheme.
  arXiv  Detail & Related papers  (2023-04-11T11:40:15Z)
• D4FT: A Deep Learning Approach to Kohn-Sham Density Functional Theory [79.50644650795012]
  We propose a deep learning approach to solve Kohn-Sham Density Functional Theory (KS-DFT) We prove that such an approach has the same expressivity as the SCF method, yet reduces the computational complexity. In addition, we show that our approach enables us to explore more complex neural-based wave functions.
  arXiv  Detail & Related papers  (2023-03-01T10:38:10Z)
• Monte Carlo Neural PDE Solver for Learning PDEs via Probabilistic Representation [59.45669299295436]
  We propose a Monte Carlo PDE solver for training unsupervised neural solvers. We use the PDEs' probabilistic representation, which regards macroscopic phenomena as ensembles of random particles. Our experiments on convection-diffusion, Allen-Cahn, and Navier-Stokes equations demonstrate significant improvements in accuracy and efficiency.
  arXiv  Detail & Related papers  (2023-02-10T08:05:19Z)
• Sampling with Mollified Interaction Energy Descent [57.00583139477843]
  We present a new optimization-based method for sampling called mollified interaction energy descent (MIED) MIED minimizes a new class of energies on probability measures called mollified interaction energies (MIEs) We show experimentally that for unconstrained sampling problems our algorithm performs on par with existing particle-based algorithms like SVGD.
  arXiv  Detail & Related papers  (2022-10-24T16:54:18Z)
• Adaptive LASSO estimation for functional hidden dynamic geostatistical model [69.10717733870575]
  We propose a novel model selection algorithm based on a penalized maximum likelihood estimator (PMLE) for functional hiddenstatistical models (f-HD) The algorithm is based on iterative optimisation and uses an adaptive least absolute shrinkage and selector operator (GMSOLAS) penalty function, wherein the weights are obtained by the unpenalised f-HD maximum-likelihood estimators.
  arXiv  Detail & Related papers  (2022-08-10T19:17:45Z)
• Universal Prony fitting decomposition for optimized hierarchical quantum master equations [9.385683681088087]
  In this work, we propose the Prony fitting decomposition (PFD) as an accurate and efficient exponential series method. The resulting hierarchical equations of motion (HEOM) formalism is greatly optimized, especially in extremely low temperature regimes. For demonstration, we calibrate the present PFD against the celebrated Pad'e spectrum decomposition method, followed by converged HEOM evaluations on the single-impurity Anderson model system.
  arXiv  Detail & Related papers  (2022-04-14T10:55:39Z)
• Taming Quantum Noise for Efficient Low Temperature Simulations of Open Quantum Systems [4.866728358750297]
  We introduce an effective treatment of quantum noise in frequency space by systematically clustering higher order Matsubara poles equivalent to an optimized rational decomposition. This leads to an elegant extension of the HEOM to arbitrary temperatures and very general reservoirs in combination with efficiency, high accuracy and long-time stability. As one highly non-trivial application, for the sub-ohmic spin-boson model at vanishing temperature the Shiba relation is quantitatively verified which predicts the long-time decay of correlation functions.
  arXiv  Detail & Related papers  (2022-02-08T18:46:11Z)
• A Near-Optimal Gradient Flow for Learning Neural Energy-Based Models [93.24030378630175]
  We propose a novel numerical scheme to optimize the gradient flows for learning energy-based models (EBMs) We derive a second-order Wasserstein gradient flow of the global relative entropy from Fokker-Planck equation. Compared with existing schemes, Wasserstein gradient flow is a smoother and near-optimal numerical scheme to approximate real data densities.
  arXiv  Detail & Related papers  (2019-10-31T02:26:20Z)

This list is automatically generated from the titles and abstracts of the papers in this site.

This site does not guarantee the quality of this site (including all information) and is not responsible for any consequences.