Source author record

Ka-Ming Tam

Ka-Ming Tam appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

21 published item(s)

preprint2025arXiv

Interplay of non-Hermitian skin effect and electronic correlations in the non-Hermitian Hubbard model via Real-space dynamical mean field theory

Non-Hermitian quantum systems, characterized by their ability to model open systems with gain and loss, have unveiled striking phenomena such as the non-Hermitian skin effect (NHSE), where eigenstates localize at boundaries under open boundary conditions. While extensively studied in non-interacting systems, the interplay between NHSE and strong electron correlations remains largely unexplored. Here, we investigate the non-Hermitian Hubbard model with asymmetric hopping, employing real-space dynamical mean-field theory (R-DMFT), a novel extension to such non-Hermitian correlated models-to capture both local correlations and spatial inhomogeneities. By analyzing the end-to-end Green's function as probes of directional amplification, we demonstrate that strong correlations can suppress the skin effect, leading to a crossover from boundary-dominated to correlation-driven dynamics. Our systematic study reveals that correlations dominate at small to intermediate strength of asymmetric hopping, inducing exponential decay in the end-to-end Green's function, but higher strength of asymmetric hopping can restore amplification even at strong interaction. These results illuminate the tunable interplay between correlations and non-Hermitian physics, suggesting avenues for engineering non-reciprocal transport in correlated open quantum systems.

preprint2022arXiv

Ab initio Approaches to High Entropy Alloys: A Comparison of CPA, SQS, and Supercell Methods

We present a comparative study of different modeling approaches to the electronic properties of the $\textrm{Hf}_{0.05}\textrm{Nb}_{0.05}\textrm{Ta}_{0.8}\textrm{Ti}_{0.05}\textrm{Zr}_{0.05}$ high entropy alloy. Common to our modeling is the methodology to compute the one-particle Green's function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e. the locally self-consistent multiple-scatering methods provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the $\textrm{Hf}_{0.05}\textrm{Nb}_{0.05}\textrm{Ta}_{0.8}\textrm{Ti}_{0.05}\textrm{Zr}_{0.05}$ high entropy alloy, instead we find possible superconductivity below $\sim 9$K.

preprint2022arXiv

Beyond Quantum Cluster Theories: Multiscale Approaches for Strongly Correlated Systems

The degrees of freedom that confer to strongly correlated systems their many intriguing properties also render them fairly intractable through typical perturbative treatments. For this reason, the mechanisms responsible for these technologically promising properties remain mostly elusive. Computational approaches have played a major role in efforts to fill this void. In particular, dynamical mean field theory (DMFT) and its cluster extension, the dynamical cluster approximation (DCA) have allowed significant progress. However, despite all the insightful results of these embedding schemes, computational constraints, such as the minus sign problem in Quantum Monte Carlo (QMC), and the exponential growth of the Hilbert space in exact diagonalization (ED) methods, still limit the length scale within which correlations can be treated exactly in the formalism. A recent advance to overcome these difficulties is the development of multiscale many body approaches whereby this challenge is addressed by introducing an intermediate length scale between the short length scale where correlations are treated exactly using a cluster solver such QMC or ED, and the long length scale where correlations are treated in a mean field manner. At this intermediate length scale correlations can be treated perturbatively. This is the essence of multiscale many-body methods. We will review various implementations of these multiscale many-body approaches, the results they have produced, and the outstanding challenges that should be addressed for further advances.

preprint2020arXiv

Deep Learning on the 2-Dimensional Ising Model to Extract the Crossover Region with a Variational Autoencoder

The 2-dimensional Ising model on a square lattice is investigated with a variational autoencoder in the non-vanishing field case for the purpose of extracting the crossover region between the ferromagnetic and paramagnetic phases. The encoded latent variable space is found to provide suitable metrics for tracking the order and disorder in the Ising configurations that extends to the extraction of a crossover region in a way that is consistent with expectations. The extracted results achieve an exceptional prediction for the critical point as well as agreement with previously published results on the configurational magnetizations of the model. The performance of this method provides encouragement for the use of machine learning to extract meaningful structural information from complex physical systems where little a priori data is available.

preprint2020arXiv

Effect of Mitigation Measures on the Spreading of COVID-19 in Hard-Hit States

State government-mandated social distancing measures have helped to slow down the growth of the COVID-19 pandemic in the United States. Current predictive models of the development of COVID-19, especially after mitigation efforts, are largely based on extrapolating the data from other countries. Since most states enforced stay-at-home orders towards the end of March, their effect should be reflected in the death and infection counts at the end of April. Using the data available until April 25th, we investigate the change in the infection rate due to the mitigation efforts, and project death and infection counts until September, 2020, for some of the most heavily impacted states: New York, New Jersey, Michigan, Massachusetts, Illinois and Louisiana. We find that with the current mitigation efforts five of those six states reduce their reproduction number to a value less than one, stopping the exponential growth of the pandemic. We also projected different scenarios after the mitigation is relaxed. Analysis for other states can be found at https://covid19projection.org/.

preprint2020arXiv

Neural Network Solver for Small Quantum Clusters

Machine learning approaches have recently been applied to the study of various problems in physics. Most of the studies are focused on interpreting the data generated by conventional numerical methods or an existing database. An interesting question is whether it is possible to use a machine learning approach, in particular a neural network, for solving the many-body problem. In this paper, we present a solver for interacting quantum problem for small clusters based on the neural network. We study the small quantum cluster which mimics the single impurity Anderson model. We demonstrate that the neural network based solver provides quantitatively accurate results for the spectral function as compared to the exact diagonalization method. This opens the possibility of utilizing the neural network approach as an impurity solver for other many body numerical approaches, such as dynamical mean field theory.

preprint2020arXiv

Non-Fermi Liquid Behaviour in the Three Dimensional Hubbard Model

We present a numerical study on non-Fermi liquid behaviour of a three dimensional system. The Hubbard model in a cubic lattice is simulated by the dynamical cluster approximation, in particular the quasi-particle weight is calculated at finite dopings for a range of temperatures. Near the putative quantum critical point, we find evidence of a separatrix at a finite doping which separates the Fermi liquid from non-Fermi liquid as the doping increases. Our results suggest that a marginal Fermi liquid and possibly a quantum critical point should exist in the three dimensions interacting Fermi system.

preprint2020arXiv

Projected Development of COVID-19 in Louisiana

At the time of writing, Louisiana has the third highest COVID-19 infection per capita in the United States. The state government issued a stay-at-home order effective March 23rd. We analyze the projected spread of COVID-19 in Louisiana without including the effects of the stay-at-home order. We predict that a large fraction of the state population would be infected without the mitigation efforts, and would certainly overwhelm the capacity of Louisiana health care system. We further predict the outcomes with different degrees of reduction in the infection rate. More than 70% of reduction is required to cap the number of infected to under one million.

preprint2015arXiv

Study of multiband disordered systems using the typical medium dynamical cluster approximation

We generalize the typical medium dynamical cluster approximation to multiband disordered systems. Using our extended formalism, we perform a systematic study of the non-local correlation effects induced by disorder on the density of states and the mobility edge of the three-dimensional two-band Anderson model. We include inter-band and intra-band hopping and an intra-band disorder potential. Our results are consistent with the ones obtained by the transfer matrix and the kernel polynomial methods. We apply the method to K$_x$Fe$_{2-y}$Se$_2$ with Fe vacancies. Despite the strong vacancy disorder and anisotropy, we find the material is not an Anderson insulator. Our results demonstrate the application of the typical medium dynamical cluster approximation method to study Anderson localization in real materials.

preprint2014arXiv

Ferromagnetic phase in the polarized two-species bosonic Hubbard Model

We recently studied a doped two-dimensional bosonic Hubbard model with two hard-core species, with different masses, using quantum Monte Carlo simulations [Phys. Rev. B 88, 161101(R) (2013)]. Upon doping away from half-filling, we find several distinct phases, including a phase-separated ferromagnet with Mott behavior for the heavy species and both Mott insulating and superfluid behaviors for the light species. Introducing polarization, an imbalance in the population between species, we find a fully phase-separated ferromagnet. This phase exists for a broad range of temperatures and polarizations. By using finite size scaling of the susceptibility, we find a critical exponent which is consistent with the two-dimensional Ising universality class. Significantly, since the global entropy of this phase is higher than that of the ferromagnetic phase with single species, its experimental observation in cold atoms may be feasible.

preprint2014arXiv

Three Dimensional Edwards-Anderson Spin Glass Model in an External Field

We study the Edwards-Anderson model on a simple cubic lattice with a finite constant external field. We employ an indicator composed of a ratio of susceptibilities at finite wavenumbers, which was recently proposed to avoid the difficulties of a zero momentum quantity, for capturing the spin glass phase transition. Unfortunately, this new indicator is fairly noisy, so a large pool of samples at low temperature and small external field are needed to generate results with sufficiently small statistical error for analysis. We thus implement the Monte Carlo method using graphics processing units to drastically speedup the simulation. We confirm previous findings that conventional indicators for the spin glass transition, including the Binder ratio and the correlation length do not show any indication of a transition for rather low temperatures. However, the ratio of spin glass susceptibilities do show crossing behavior, albeit a systematic analysis is beyond the reach of the present data. This calls for a more thorough study of the three-dimension Edwards-Anderson model in an external field.

preprint2013arXiv

Complex phases in the doped two-species bosonic Hubbard Model

We study a two-dimensional bosonic Hubbard model with two hard-core species away from half filling using Quantum Monte Carlo simulations. The model includes a repulsive interspecies interaction and different nearest-neighbor hopping terms for the two species. By varying the filling we find a total of five distinct phases, including a normal liquid phase at higher temperature, and four different phases at lower temperature. We find an anti-ferromagnetically ordered Mott insulator and a region of coexistent anti-ferromagnetic and superfluid phases near half filling. Further away from half filling the phase diagram displays a superfluid phase and a novel phase inside the superfluid region at even lower temperatures. In this novel phase separated region, the heavy species has a Mott behavior with integer filling, while the lighter species shows phase separated Mott and superfluid behaviors.

preprint2013arXiv

Dominant Superconducting Fluctuations in the One-Dimensional Extended Holstein-Extended Hubbard model

The search for realistic one-dimensional (1D) models that exhibit dominant superconducting (SC) fluctuations effects has a long history. In these 1D systems, the effects of commensurate band fillings--strongest at half-filling--and electronic repulsions typically lead to a finite charge gap and the favoring of insulating density wave ordering over superconductivity. Accordingly, recent proposals suggesting a gapless metallic state in the Holstein-Hubbard (HH) model, possibly superconducting, have generated considerable interest and controversy, with the most recent work demonstrating that the putative dominant superconducting state likely does not exist. In this paper we study a model with non-local electron-phonon interactions, in addition to electron-electron interactions, this model unambiguously possesses dominant superconducting fluctuations at half filling in a large region of parameter space. Using both the numerical multi-scale functional renormalization group for the full model and an analytic conventional renormalization group for a bosonized version of the model, we demonstrate the existence of dominant superconducting (SC) fluctuations. These dominant SC fluctuations arise because the spin-charge coupling at high energy is weakened by the non-local electron-phonon interaction and the charge gap is destroyed by the resultant suppression of the Umklapp process. The existence of the dominant SC pairing instability in this half-filled 1D system suggests that non-local boson-mediated interactions may be important in the superconductivity observed in the organic superconductors.

preprint2013arXiv

Parallel Tempering Simulation of the three-dimensional Edwards-Anderson Model with Compact Asynchronous Multispin Coding on GPU

Monte Carlo simulations of the Ising model play an important role in the field of computational statistical physics, and they have revealed many properties of the model over the past few decades. However, the effect of frustration due to random disorder, in particular the possible spin glass phase, remains a crucial but poorly understood problem. One of the obstacles in the Monte Carlo simulation of random frustrated systems is their long relaxation time making an efficient parallel implementation on state-of-the-art computation platforms highly desirable. The Graphics Processing Unit (GPU) is such a platform that provides an opportunity to significantly enhance the computational performance and thus gain new insight into this problem. In this paper, we present optimization and tuning approaches for the CUDA implementation of the spin glass simulation on GPUs. We discuss the integration of various design alternatives, such as GPU kernel construction with minimal communication, memory tiling, and look-up tables. We present a binary data format, Compact Asynchronous Multispin Coding (CAMSC), which provides an additional $28.4\%$ speedup compared with the traditionally used Asynchronous Multispin Coding (AMSC). Our overall design sustains a performance of 33.5 picoseconds per spin flip attempt for simulating the three-dimensional Edwards-Anderson model with parallel tempering, which significantly improves the performance over existing GPU implementations.

preprint2013arXiv

Phase diagram of the Bose-Hubbard model on a ring-shaped lattice with tunable weak links

Motivated by recent experiments on toroidal Bose-Einstein condensates in all-optical traps with tunable weak links, we study the one-dimensional Bose-Hubbard model on a ring-shaped lattice with a small region of weak hopping integrals using quantum Monte Carlo simulations. Besides the usual Mott insulating and superfluid phases, we find a phase which is compressible but non superfluid with a local Mott region. This `local Mott' phase extends in a large region of the phase diagram. These results suggest that the insulating and conducting phases can be tuned by a local parameter which may provide a new insight to the design of atomtronic devices.

preprint2013arXiv

The Boson-Hubbard Model on a Kagome Lattice with Sextic Ring-Exchange Terms

High order ring-exchange interactions are crucial for the study of quantum fluctuations on highly frustrated systems. We present the first exact quantum Monte Carlo study of a model of hard-core bosons with sixth order ring-exchange interactions on a two-dimensional kagome lattice. By using the Stochastic Green Function algorithm, we show that the system becomes unstable in the limit of large ring-exchange interactions. It undergoes a phase separation at all fillings, except at 1/3 and 2/3 fillings for which the superfluid density vanishes and an unusual mixed valence bond and charge density ordered solid is formed.

preprint2012arXiv

Periodic Anderson model with electron-phonon correlated conduction band

This paper reports dynamical mean field calculations for the periodic Anderson model in which the conduction band is coupled to phonons. Motivated in part by recent attention to the role of phonons in the $γ$-$α$ transition in Ce, this model yields a rich and unexpected phase diagram which is of intrinsic interest. Specifically, above a critical value of the electron-phonon interaction, a first order transition with two coexisting phases develops in the temperature-hybridization plane, which terminates at a second order critical point. The coexisting phases display the familiar Kondo screened and local moment character, yet they also exhibit pronounced polaronic and bipolaronic properties, respectively.

preprint2011arXiv

On the Validity of the Tomonaga Luttinger Liquid Relations for the One-dimensional Holstein Model

For the one-dimensional Holstein model, we show that the relations among the scaling exponents of various correlation functions of the Tomonaga Luttinger liquid (LL), while valid in the thermodynamic limit, are significantly modified by finite size corrections. We obtain analytical expressions for these corrections and find that they decrease very slowly with increasing system size. The interpretation of numerical data on finite size lattices in terms of LL theory must therefore take these corrections into account. As an important example, we re-examine the proposed metallic phase of the zero-temperature, half-filled one-dimensional Holstein model without employing the LL relations. In particular, using quantum Monte Carlo calculations, we study the competition between the singlet pairing and charge ordering. Our results do not support the existence of a dominant singlet pairing state.

preprint2011arXiv

Solving the Parquet Equations for the Hubbard Model beyond Weak Coupling

We find that imposing the crossing symmetry in the iteration process considerably extends the range of convergence for solutions of the parquet equations for the Hubbard model. When the crossing symmetry is not imposed, the convergence of both simple iteration and more complicated continuous loading (homotopy) methods are limited to high temperatures and weak interactions. We modify the algorithm to impose the crossing symmetry without increasing the computational complexity. We also imposed time reversal and a subset of the point group symmetries, but they did not further improve the convergence. We elaborate the details of the latency hiding scheme which can significantly improve the performance in the computational implementation. With these modifications, stable solutions for the parquet equations can be obtained by iteration more quickly even for values of the interaction that are a significant fraction of the bandwidth and for temperatures that are much smaller than the bandwidth. This may represent a crucial step towards the solution of two-particle field theories for correlated electron models.

preprint2010arXiv

Thermodynamic Spin Glass Phase Induced by Weak Random Exchange Disorder in a Classical Spin Liquid: the Case of the Pyrochlore Heisenberg Antiferromagnet

The glassy behavior observed in the pyrochlore magnet Y2Mo2O7, where the magnetic Mo^{4+} ions interact predominantly via isotropic nearest neighbor antiferromagnetic exchange, possibly with additional weak disorder, is a distinct class of spin glass systems where frustration is mostly geometrical. A model proposed to describe such a spin glass behavior is the Heisenberg model on a pyrochlore lattice with random but strictly antiferromagnetic exchange disorder. In this paper, we provide compelling numerical evidence from extensive Monte Carlo simulations which show that the model exhibits a finite temperature spin glass transition and thus is a realization of a spin glass induced by random weak disorder from spin liquid. From our results, we are led to suggest that the spin glass state of Y2Mo2O7 is driven by effective strong disorder.

preprint2009arXiv

Evolution of pairing from weak to strong coupling on a honeycomb lattice

We study the evolution of the pairing from weak to strong coupling on a honeycomb lattice by Quantum Monte Carlo. We show numerical evidence of the BCS-BEC crossover as the coupling strength increases on a honeycomb lattice with small fermi surface by measuring a wide range of observables: double occupancy, spin susceptibility, local pair correlation, and kinetic energy. Although at low energy, the model sustains Dirac fermions, we do not find significant qualitative difference in the BCS-BEC crossover as compared to those with an extended Fermi surface, except at weak coupling, BCS regime.