Source author record

Ning-Hua Tong

Ning-Hua Tong 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

17works
7topics
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

17 published item(s)

preprint2022arXiv

Interacting spinless fermions on the square lattice: Charge order, phase separation, and superconductivity

We investigate the phase diagram of spinless fermions on a square lattice with nearest-neighbor interaction, using the recently developed projective truncation approximation in Green's function equation of motion. For attractive interaction, the ground state is in an homogeneous p + ip superconducting (SC) phase at high or low electron densities. Near half filling is a phase separation (PS) between the SC phases. Allowing inhomogeneous solution, we obtain p-wave SC domains with positive interface energy. As temperature increases, the SC phases transit into normal phases above Tsc, generating an homogeneous normal phase (far away from n = 1/2), or a PS between normal phases with different densities (close to n = 1/2). Further increasing temperature to Tps, the PS disappears and the particle-hole symmetry of the Hamiltonian is recovered. For repulsive interaction, depending on electron filling, the ground state is in charge-ordered phase (half filling), charge-disordered phase (large hole/electron doping), or PS between them (weak doping). At finite temperature, the regime of charge order phase moves to finite V and extends to finite doping regime.

preprint2022arXiv

Projective-truncation-approximation study of the one-dimensional $ϕ^4$ lattice model

In this paper, we first develop the projective truncation approximation (PTA) in the Green's function equation of motion (EOM) formalism for classical statistical models. To implement PTA for a given Hamiltonian, we choose a set of basis variables and projectively truncate the hierarchical EOM. We apply PTA to the one-dimensional $ϕ^4$ lattice model. Phonon dispersion and static correlation functions are studied in detail. Using one- and two-dimensional bases, we obtain results identical to and beyond the quadratic variational approximation, respectively. In particular, we analyze the power-law temperature dependence of the static averages in the low- and high-temperature limits, and we give exact exponents.

preprint2020arXiv

Equilibrium Dynamics of the Sub-Ohmic Spin-boson Model At Finite Temperature

We use the full-density matrix (FDM) numerical renormalization group (NRG) method to calculate the equilibrium dynamical correlation function $C(ω)$ of the spin operator $σ_z$ at finite temperature for the sub-Ohmic spin-boson model. A peak is observed at the frequency $ω_{T}\sim T$ in the curve of $C(ω)$. The curve merges with the zero temperature $C(ω)$ in $ω\gg ω_{T}$ and deviate significantly from the power-law form $ω^{\pm s}$ of the zero temperature curve in $ω\llω_{T}$.

preprint2019arXiv

Improved strong-coupling perturbation theory of the symmetric Anderson impurity model

In a previous work (N. H. Tong, Phys. Rev. B 92, 165126 (2015)), an equation-of-motion based series expansion formalism was used to do the second-order strong-coupling expansion for the single-particle Green function of the Anderson impurity model. In this paper, we improve this theory in two aspects. We first use a more accurate scheme to self-consistently calculate the averages that appear in G1. In the resummation process, we use updated coefficients for the continued fraction, guided by the formally exact continued fraction from the Mori-Zwanzig theory. These changes lead to more accurate impurity spin response to the magnetic bias of the bath. Combined with the dynamical mean-field theory, our theory gives improved description for the antiferromagnetism of Hubbard model at half filling.

preprint2018arXiv

Controllable Precision of the Projective Truncation Approximation for Green's Functions

Recently, we developed the projective truncation approximation for the equation of motion of two-time Green's functions (P. Fan et al., Phys. Rev. B 97, 165140 (2018)). In that approximation, the precision of results depends on the selection of operator basis. Here, for three successively larger operator bases, we calculate the local static averages and the impurity density of states of the single-band Anderson impurity model. The results converge systematically towards those of numerical renormalization group as the basis size is enlarged. We also propose a quantitative gauge of the truncation error within this method and demonstrate its usefulness using the Hubbard-I basis. We thus confirm that the projective truncation approximation is a method of controllable precision for quantum many-body systems.

preprint2015arXiv

A New Monte Carlo Algorithm for Free Energy Calculation

We propose a new Monte Carlo algorithm for the free energy calculation based on configuration space sampling. We implement this algorithm for Ising model. Comparison with the exact free energy shows an excellent agreement. We analyse the properties of this algorithm and compares it with Wang-Landau algorithm which samples in energy space. This method is applicable to classical statistical models. The possibility of extending it to quantum systems is discussed.

preprint2015arXiv

A Standard Basis Operator Equation of Motion Impurity Solver for Dynamical Mean Field Theory

We present an efficient impurity solver for the dynamical mean-field theory (DMFT). It is based on the separation of bath degrees of freedom into the low energy and the high energy parts. The former is solved exactly using exact diagonalization and the latter is treated approximately using Green's function equation of motion decoupling approximation. The two parts are combined coherently under the standard basis operator formalism. The impurity solver is applied to the Anderson impurity model and, combined with DMFT, to the one-band Hubbard model. Qualitative agreement is found with other well established methods. Some promising features and possible improvements of the present solver are discussed.

preprint2015arXiv

Equation of Motion Series Expansion of Double Time Green's Functions

Based on the Green's function (GF) equation-of-motion formalism, we develop a method to expand the double time Green's function into Taylor series of the parameter $λ$ in the Hamiltonian $H=H_0 + λH_1$. Here $H_0$ is the exactly solvable part and $H_1$ is regarded as the perturbation. To restore the analytical structure of GF, we use the continued fraction to do resummation for the obtained series. The problem of zero-temperature divergence is identified and remedied by the self-consistent series expansion. To demonstrate the implementation of this method, we carry out the weak- as well as the strong-coupling expansion for the Anderson impurity model to order $λ^2$. Improved result for the local density of states is obtained by self-consistent second-order strong-coupling expansion and continued fraction resummation.

preprint2015arXiv

Spin susceptibility of Anderson impurities in arbitrary conduction bands

Spin susceptibility of Anderson impurities is a key quantity in understanding the physics of Kondo screening. Traditional numerical renormalization group (NRG) calculation of the impurity contribution $χ_{\textrm{imp}}$ to susceptibility, defined originally by Wilson in a flat wide band, has been generalized before to structured conduction bands. The results brought about non-Fermi-liquid and diamagnetic Kondo behaviors in $χ_{\textrm{imp}}$, even when the bands are not gapped at the Fermi energy. Here, we use the full density-matrix (FDM) NRG to present high-quality data for the local susceptibility $χ_{\textrm{loc}}$ and to compare them with $χ_{\textrm{imp}}$ obtained by the traditional NRG. Our results indicate that those exotic behaviors observed in $χ_{\textrm{imp}}$ are unphysical. Instead, the low-energy excitations of the impurity in arbitrary bands only without gap at the Fermi energy are still a Fermi liquid and paramagnetic. We also demonstrate that unlike the traditional NRG yielding $χ_{\textrm{loc}}$ less accurate than $χ_{\textrm{imp}}$, the FDM method allows a high-precision dynamical calculation of $χ_{\textrm{loc}}$ at much reduced computational cost, with an accuracy at least one order higher than $χ_{\textrm{imp}}$. Moreover, artifacts in the FDM algorithm to $χ_{\textrm{imp}}$, and origins of the spurious non-Fermi-liquid and diamagnetic features are clarified. Our work provides an efficient high-precision algorithm to calculate the spin susceptibility of impurity for arbitrary structured bands, while negating the applicability of Wilson's definition to such cases.

preprint2013arXiv

Cluster mean-field theory study of $J_{1}-J_{2}$ Heisenberg model on a square lattice

We study the spin-1/2 $J_{1}$-$J_{2}$ Heisenberg model on a square lattice using the cluster mean-field theory. We find a rapid convergence of phase boundaries with increasing cluster size. By extrapolating the cluster size $L$ to infinity, we obtain accurate phase boundaries $J_{2}^{c1} \approx 0.42$ (between the N$\acute{e}$el antiferromagnetic phase and nonmagnetic phase), and $J_{2}^{c2} \approx 0.59$ (between nonmagnetic phase and the collinear antiferromagnetic phase). The transitions are identified unambiguously as second order at $J_{2}^{c1}$ and first order at $J_{2}^{c2}$. At finite temperature, we present a complete phase diagram with stable, meta-stable and unstable states near $J_{2}^{c2}$, being relevant to that of the anisotropic $J_1-J_2$ model. The uniform as well as staggered magnetic susceptibilities are also discussed.

preprint2013arXiv

Local Entanglement Entropy at the Mott Metal-Insulator Transition in Infinite Dimensions

We study the critical behavior of the single-site entanglement entropy S at the Mott metal-insulator transition in infinite-dimensional Hubbard model. For this model, the entanglement between a single site and rest of the lattice can be evaluated exactly, using the dynamical mean-field theory (DMFT). Both the numerical solution using exact diagonalization and the analytical one using two-site DMFT gives S-Sc \propto α\log_{2}[(1/2-Dc)/Dc](U-Uc), with Dc the double occupancy at Uc and α< 0 being different on two sides of the transition.

preprint2012arXiv

Scaling Analysis in the Numerical Renormalization Group Study of the Sub-Ohmic Spin-Boson Model

The spin-boson model has nontrivial quantum phase transitions in the sub-Ohmic regime. For the bath spectra exponent $0 \leqslant s<1/2$, the bosonic numerical renormalization group (BNRG) study of the exponents $β$ and $δ$ are hampered by the boson state truncation which leads to artificial interacting exponents instead of the correct Gaussian ones. In this paper, guided by a mean-field calculation, we study the order parameter function $m(τ=α-α_c, ε, Δ)$ using BNRG. Scaling analysis with respect to the boson state truncation $N_{b}$, the logarithmic discretization parameter $Λ$, and the tunneling strength $Δ$ are carried out. Truncation-induced multiple-power behaviors are observed close to the critical point, with artificial values of $β$ and $δ$. They cross over to classical behaviors with exponents $β=1/2$ and $δ=3$ on the intermediate scales of $τ$ and $ε$, respectively. We also find $τ/Δ^{1-s}$ and $ε/Δ$ scalings in the function $m(τ, ε, Δ)$. The role of boson state truncation as a scaling variable in the BNRG result for $0 \leqslant s<1/2$ is identified and its interplay with the logarithmic discretization revealed. Relevance to the validity of quantum-to-classical mapping in other impurity models is discussed.

preprint2010arXiv

Criticality of the Mean-Field Spin-Boson Model: Boson State Truncation and Its Scaling Analysis

The spin-boson model has nontrivial quantum phase transitions at zero temperature induced by the spin-boson coupling. The bosonic numerical renormalization group (BNRG) study of the critical exponents $β$ and $δ$ of this model is hampered by the effects of boson Hilbert space truncation. Here we analyze the mean-field spin boson model to figure out the scaling behavior of magnetization under the cutoff of boson states $N_{b}$. We find that the truncation is a strong relevant operator with respect to the Gaussian fixed point in $0<s<1/2$ and incurs the deviation of the exponents from the classical values. The magnetization at zero bias near the critical point is described by a generalized homogeneous function (GHF) of two variables $τ=α-α_{c}$ and $x=1/N_{b}$. The universal function has a double-power form and the powers are obtained analytically as well as numerically. Similarly, $m(α=α_{c})$ is found to be a GHF of $ε$ and $x$. In the regime $s>1/2$, the truncation produces no effect. Implications of these findings to the BNRG study are discussed.

preprint2010arXiv

Interacting Dirac Fermions on Honeycomb Lattice

We investigate the interacting Dirac fermions on honeycomb lattice by cluster dynamical mean-field theory (CDMFT) combined with continuous time quantum Monte Carlo simulation (CTQMC). A novel scenario for the semimetal-Mott insulator transition of the interacting Dirac fermions is found beyond the previous DMFT studies. We demonstrate that the non-local spatial correlations play a vital role in the Mott transition on the honeycomb lattice. We also elaborate the experimental protocol to observe this phase transition by the ultracold atoms on optical honeycomb lattice.

preprint2010arXiv

Simulation of the spin-boson model with superconducting phase qubit coupled to a transmission line

Based on the rapid experimental developments of circuit QED, we propose a feasible scheme to simulate a spin-boson model with the superconducting circuits, which can be used to detect quantum Kosterlitz-Thouless (KT) phase transition. We design the spin-boson model by using a superconducting phase qubit coupled with a semi-infinite transmission line, which is regarded as bosonic reservoir with a continuum spectrum. By tuning the bias current or the coupling capacitance, the quantum KT transition can be directly detected through tomography measurement on the states of the phase qubit. We also estimate the experimental parameters using numerical renormalization group method.

preprint2009arXiv

Dynamical Mean Field Theory for the Bose-Hubbard Model

The dynamical mean field theory (DMFT), which is successful in the study of strongly correlated fermions, was recently extended to boson systems [Phys. Rev. B {\textbf 77}, 235106 (2008)]. In this paper, we employ the bosonic DMFT to study the Bose-Hubbard model which describes on-site interacting bosons in a lattice. Using exact diagonalization as the impurity solver, we get the DMFT solutions for the Green's function, the occupation density, as well as the condensate fraction on a Bethe lattice. Various phases are identified: the Mott insulator, the Bose-Einstein condensed (BEC) phase, and the normal phase. At finite temperatures, we obtain the crossover between the Mott-like regime and the normal phase, as well as the BEC-to-normal phase transition. Phase diagrams on the $μ/U-\tilde{t}/U$ plane and on the $T/U-\tilde{t}/U$ plane are produced ($\tilde{t}$ is the scaled hopping amplitude). We compare our results with the previous ones, and discuss the implication of these results to experiments.

preprint2005arXiv

Extended Variational Cluster Approximation

The variational cluster approximation (VCA) proposed by M. Potthoff {\it et al.} [Phys. Rev. Lett. {\bf 91}, 206402 (2003)] is extended to electron or spin systems with nonlocal interactions. By introducing more than one source field in the action and employing the Legendre transformation, we derive a generalized self-energy functional with stationary properties. Applying this functional to a proper reference system, we construct the extended VCA (EVCA). In the limit of continuous degrees of freedom for the reference system, EVCA can recover the cluster extension of the extended dynamical mean-field theory (EDMFT). For a system with correlated hopping, the EVCA recovers the cluster extension of the dynamical mean-field theory for correlated hopping. Using a discrete reference system composed of decoupled three-site single impurities, we test the theory for the extended Hubbard model. Quantitatively good results as compared with EDMFT are obtained. We also propose VCA (EVCA) based on clusters with periodic boundary conditions. It has the (extended) dynamical cluster approximation as the continuous limit. A number of related issues are discussed.