Researcher profile

Aaram J. Kim

Aaram J. Kim contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

8 published item(s)

preprint2022arXiv

Misleading convergence of the skeleton diagrammatic technique: when the correct solution can be found

Convergence of the self-consistent skeleton diagrammatic technique (SDT) -- in which the full Green's function is determined through summation of Feynman diagrams in terms of itself -- to the wrong answer has been associated with the existence of non-perturbative branches of the Luttinger-Ward functional. Although it has been possible to detect misleading convergence without the knowledge of the exact result, the SDT has remained inapplicable in the regimes where this happens. We show that misleading convergence does not always preclude recovering the exact solution. In addition to the established mechanism, convergence of the SDT to the wrong answer can stem from divergence of the inherent diagrammatic series, which allows us to recover the exact solution by a modified SDT protocol based on controlled analytic continuation. We illustrate this approach by its application to the analytically solvable (0+0)d Hubbard model, the Hubbard atom, and the 2d Hubbard model in a challenging strong-coupling regime, for which the SDT is solved with controlled accuracy by the diagrammatic Monte Carlo (DiagMC) technique.

preprint2022arXiv

Pseudoparticle vertex solver for quantum impurity models

We present a quantum impurity solver based on a pseudo-particle framework, which combines diagrammatic resummations for a three-point vertex with diagrammatic Monte Carlo sampling of a four-point vertex. This recently proposed approach [A. J. Kim et al., arXiv:2112.15549] is generalized here to fermionic impurity problems and we discuss the technical details of the implementation, including the time-stepping approach, the Monte Carlo updates, and the routines for checking the two-particle irreducibility of the four-point vertex. We also explain how the vertex information can be efficiently stored using a Dubiner basis representation. The convergence properties of the algorithm are demonstrated with applications to exactly solvable impurity models and dynamical mean field theory simulations of the single-orbital Hubbard model. It is furthermore shown that the algorithm can handle a two-orbital problem with off-diagonal hybridizations, which would cause a severe sign problem in standard hybridization-expansion Monte Carlo simulations. Since the vertex-based algorithm successfully handles sign oscillating integrals in equilibrium and samples only connected diagrams, it may be a promising approach for real-time simulations.

preprint2021arXiv

Evaluating second-order phase transitions with Diagrammatic Monte Carlo: Néel Transition in the doped three-dimensional Hubbard model

Diagrammatic Monte Carlo -- the technique for numerically exact summation of all Feynman diagrams to high orders -- offers a unique unbiased probe of continuous phase transitions. Being formulated directly in the thermodynamic limit, the diagrammatic series is bound to diverge and is not resummable at the transition due to the non-analyticity of physical observables. This enables the detection of the transition with controlled error bars from an analysis of the series coefficients alone, avoiding the challenge of evaluating physical observables near the transition. We demonstrate this technique by the example of the Néel transition in the $3d$ Hubbard model. At half-filling and higher temperatures, the method matches the accuracy of state-of-the-art finite-size techniques, but surpasses it at low temperatures and allows us to map the phase diagram in the doped regime, where finite-size techniques struggle from the fermion sign problem. At low temperatures and sufficient doping, the transition to an incommensurate spin density wave state is observed.

preprint2021arXiv

Vertex-based Diagrammatic Treatment of Light-Matter-Coupled Systems

We propose a diagrammatic Monte Carlo approach for general spin-boson models, which can be regarded as a generalization of the strong-coupling expansion for fermionic impurity models. The algorithm is based on a self-consistently computed three-point vertex and a stochastically sampled four-point vertex, and achieves convergence to the numerically exact result in a wide parameter regime. The performance of the algorithm is demonstrated with applications to a spin-boson model representing an emitter in a waveguide. As a function of the coupling strength, the spin exhibits a delocalization-localization crossover at low temperatures, signaling a qualitative change in the real-time relaxation. In certain parameter regimes, the response functions of the emitter coupled to the electromagnetic continuum can be described by an effective Rabi model with appropriately defined parameters. We also discuss the spatial distribution of the photon density around the emitter.

preprint2020arXiv

Alleviating the Sign Problem in Quantum Monte Carlo Simulations of Spin-Orbit-Coupled Multi-Orbital Hubbard Models

We present a strategy to alleviate the sign problem in continuous-time quantum Monte Carlo (CTQMC) simulations of the dynamical-mean-field-theory (DMFT) equations for the spin-orbit-coupled multiorbital Hubbard model. We first identify the combinations of rotationally invariant Hund coupling terms present in the relativistic basis which lead to a severe sign problem. Exploiting the fact that the average sign in CTQMC depends on the choice of single-particle basis, we propose a bonding-antibonding basis $V_{j3/2\mathrm{BA}}$ which shows an improved average sign compared to the widely used relativistic basis for most parameter sets investigated. We then generalize this procedure by introducing a stochastic optimization algorithm that exploits the space of single-particle bases and show that $V_{j3/2\mathrm{BA}}$ is very close to optimal within the parameter space investigated. Our findings enable more efficient DMFT simulations of materials with strong spin-orbit coupling.

preprint2020arXiv

Extended crossover from Fermi liquid to quasi-antiferromagnet in the half-filled 2D Hubbard model

The ground state of the Hubbard model with nearest-neighbor hopping on the square lattice at half filling is known to be that of an antiferromagnetic (AFM) band insulator for any on-site repulsion. At finite temperature, the absence of long-range order makes the question of how the interaction-driven insulator is realized nontrivial. We address this problem with controlled accuracy in the thermodynamic limit using self-energy diagrammatic determinant Monte Carlo and dynamical cluster approximation methods and show that development of long-range AFM correlations drives an extended crossover from Fermi liquid to insulating behavior in the parameter regime that precludes a metal-to-insulator transition. The intermediate crossover state is best described as a non-Fermi liquid with a partially gapped Fermi surface.

preprint2020arXiv

Spin and charge correlations across the metal-to-insulator crossover in the half-filled $2d$ Hubbard model

The $2d$ Hubbard model with nearest-neighbour hopping on the square lattice and an average of one electron per site is known to undergo an extended crossover from metallic to insulating behavior driven by proliferating antiferromagnetic correlations. We study signatures of this crossover in spin and charge correlation functions and present results obtained with controlled accuracy using diagrammatic Monte Carlo in the range of parameters amenable to experimental verification with ultracold atoms in optical lattices. The qualitative changes in charge and spin correlations associated with the crossover are observed at well-separated temperature scales, which encase the intermediary regime of non-Fermi-liquid character, where local magnetic moments are formed and non-local fluctuations in both channels are essential.

preprint2018arXiv

Multi-valuedness of the Luttinger-Ward functional in the Fermionic and Bosonic System with Replicas

We study the properties of the Luttinger-Ward functional (LWF) in a simplified Hubbard-type model without time or spatial dimensions, but with $N$ identical replicas located on a single site. The simplicity of this $(0+0)d$ model permits an exact solution for all $N$ and for both bosonic and fermionic statistics. We show that fermionic statistics are directly linked to the fact that multiple values of the noninteracting Green function $G_0$ map to the same value of the interacting Green function $G$, i.e. the mapping $G_0 \mapsto G$ is non-injective. This implies that with fermionic statistics the $(0+0)d$ model has a multiply-valued LWF. The number of LWF values in the fermionic model increases proportionally to the number of replicas $N$, while in the bosonic model the LWF has a single value regardless of $N$. We also discuss the formal connection between the $(0+0)d$ model and the $(0+1)d$ model which was used in previous studies of LWF multivaluedness.