Lattice Fermions Research Articles

Equality of Magnetization and Edge Current for Interacting Lattice Fermions at Positive Temperature

We prove that the magnetization is equal to the edge current in the thermodynamic limit for a large class of models of lattice fermions with finite-range interactions satisfying local indistinguishability of the Gibbs state, a condition known to hold for sufficiently high temperatures. Our result implies that edge currents in such systems are determined by bulk properties and are therefore stable against large perturbations near the boundaries. Moreover, the equality persists also after taking the derivative with respect to the chemical potential. We show that this form of bulk-edge correspondence is essentially a consequence of homogeneity in the bulk and locality of the Gibbs state. An important intermediate result is a new version of Bloch’s theorem for two-dimensional systems, stating that persistent currents vanish in the bulk.

Luttinger liquid tensor network: Sine versus tangent dispersion of massless Dirac fermions

To apply the powerful many-body techniques of tensor networks to massless Dirac fermions, one wants to discretize the p·σ Hamiltonian and construct a matrix-product-operator (MPO) representation. We compare two alternative discretization schemes, one with a sine dispersion, the other with a tangent dispersion, applied to a one-dimensional Luttinger liquid with Hubbard interaction. Both types of lattice fermions allow for an exact MPO representation of low bond dimension, so they are efficiently computable, but only the tangent dispersion gives a power law decay of the propagator in agreement with the continuum limit: The sine dispersion is gapped by the interactions, evidenced by an exponentially decaying propagator. Our construction of a tensor network with an unpaired Dirac cone works around the fermion-doubling obstruction by exploiting the fact that the Hamiltonian of tangent fermions permits a generalized eigenproblem. Published by the American Physical Society 2024

Spectral functions of lattice fermions on the honeycomb lattice with Hubbard and long-range Coulomb interactions

Spectral functions of lattice fermions on the honeycomb lattice with Hubbard and long-range Coulomb interactions

Stout smearing and Wilson flow in lattice perturbation theory

We present the expansion of stout smearing and the Wilson flow in lattice perturbation theory to order g03, which is suitable for one-loop calculations. As the Wilson flow is generated by infinitesimal stout smearing steps, the results are related to each other by taking the appropriate limits. We show how to apply perturbative stout smearing or Wilson flow to the Feynman rules of any lattice fermion action and illustrate them by calculating the self-energy of the clover-improved Wilson fermion. Published by the American Physical Society 2024

Non-Hermitian lattice fermions in the 2D Gross-Neveu-Yukawa model

We work the lattice fermions and non-Hermitian formulation in the 2D Gross-Neveu-Yukawa (GNY) model and demonstrate the numerical implementation for two flavors by the hybrid Monte Carlo. Our approach has a notable advantage in dealing with chiral symmetry on a lattice by avoiding the Nielsen-Ninomiya theorem, due to the nonsymmetrized finite-difference operator. We restore the hypercubic symmetry by averaging over all possible orientations with the proper continuum limit. Our study is the first simulation for the interacting fermion formulated in a non-Hermitian way. We compare the numerical solution with the one-loop resummation. The resummation results matches with the numerical solution in ⟨ϕ⟩, ⟨ϕ2⟩, ⟨Tr(ψ¯1ψ1+ψ¯2ψ2)/2⟩, and ⟨Tr(ψ¯1ψ1+ψ¯2ψ2)ϕ/2⟩. We also used the one-loop resummation to provide the renormalization group flow and asymptotic safety in the 2D GNY model. Published by the American Physical Society 2024

InGaAs quantum dot chains grown by twofold selective area molecular beam epitaxy

Increasing quantum confinement in semiconductor quantum dot (QD) systems is essential to perform robust simulations of many-body physics. By combining molecular beam epitaxy and lithographic techniques, we developed an approach consisting of a twofold selective area growth to build QD chains. Starting from 15 nm-thick and 65 nm-wide in-plane In0.53Ga0.47As nanowires on InP substrates, linear arrays of In0.53Ga0.47As QDs were grown on top, with tunable lengths and separations. Kelvin probe force microscopy performed at room temperature revealed a change of quantum confinement in chains with decreasing QD sizes, which was further emphasized by the spectral shift of quantum levels resolved in the conduction band with low temperature scanning tunneling spectroscopy. This approach, which allows the controlled formation of 25 nm-thick QDs with a minimum length and separation of 30 nm and 22 nm respectively, is suitable for the construction of scalable fermionic quantum lattices.

Low-depth simulations of fermionic systems on square-grid quantum hardware

We present a general strategy for mapping fermionic systems to quantum hardware with square qubit connectivity which yields low-depth quantum circuits, counted in the number of native two-qubit fSIM gates. We achieve this by leveraging novel operator decomposition and circuit compression techniques paired with specifically chosen low-depth fermion-to-qubit mappings and allow for a high degree of gate cancellations and parallelism. Our mappings retain the flexibility to simultaneously optimize for qubit counts or qubit operator weights and can be used to investigate arbitrary fermionic lattice geometries. We showcase our approach by investigating the tight-binding model, the Fermi-Hubbard model as well as the multi-orbital Hubbard-Kanamori model. We report unprecedentedly low circuit depths per single Trotter layer with up to a 70% improvement upon previous state-of-the-art. Our compression technique also results in significant reduction of two-qubit gates. We find the lowest gate-counts when applying the XYZ-formalism to the DK mapping. Additionally, we show that our decomposition and compression formalism produces favourable circuits even when no native parameterized two-qubit gates are available.

Frustration-induced itinerant ferromagnetism of fermions in optical lattices

Frustration-induced itinerant ferromagnetism of fermions in optical lattices

Topologically protected Casimir effect for lattice fermions

The electromagnetic Casimir effect has a fermionic counterpart in topological insulators: Zero-point fluctuations of a massless Dirac fermion field mediate a force between magnetic scatterers. The Casimir force is insensitive to disorder that preserves the topological protection of an unpaired Dirac cone. The protection may be broken if the Dirac equation is discretized, and an exponential suppression of the Casimir effect will result if a gap opens at the Dirac point. Here, we show how this lattice artifact may be avoided, by applying a recently developed local discretization of the Euclidean action that does not suffer from the fermion-doubling obstruction of local discretizations of the Hamiltonian. Published by the American Physical Society 2024

Floquet Insulators and Lattice Fermions.

Floquet insulators are periodically driven quantum systems that can host novel topological phases as a function of the drive parameters. These new phases exhibit features reminiscent of fermion doubling in discrete-time lattice fermion theories. We make this suggestion concrete by mapping the spectrum of a noninteracting (1+1)D Floquet insulator for certain drive parameters onto that of a discrete-time lattice fermion theory with a time-independent Hamiltonian. The resulting Hamiltonian is distinct from the Floquet Hamiltonian that generates stroboscopic dynamics. It can take the form of a discrete-time Su-Schrieffer-Heeger model with half the number of spatial sites of the original model, or of a (1+1)D Wilson-Dirac theory with one quarter of the spatial sites.

Exploring kinetically induced bound states in triangular lattices with ultracold atoms: Spectroscopic approach

Quantum simulations with ultracold fermions in triangular optical lattices have recently emerged as a new platform for studying magnetism in frustrated systems. Experimental realizations of the Fermi Hubbard model revealed striking contrast between magnetism in bipartite and triangular lattices. In bipartite lattices magnetism is strongest at half filling, and doped charge carriers tend to suppress magnetic correlations. In triangular-type lattices for large U/t and t>0, antiferromagnetism (ferromagnetism) gets enhanced by doping away from n=1 with holes (doublons) because kinetic energy of dopants can be lowered through developing magnetic correlations, corresponding to formation of magnetic polarons. Snapshots of many-body states obtained with quantum gas microscopes demonstrated existence of magnetic polarons by revealing the magnetic correlations around dopants at temperatures that considerably exceed superexchange energy scale. In this paper we discuss theoretically that additional insight into properties of magnetic polarons in triangular lattices can be achieved using spectroscopic experiments with ultracold atoms. We consider starting from a spin polarized state with small hole doping and applying a two-photon Raman photoexcitation, which transfers atoms into a different spin state. We show that such magnon injection spectra exhibit a separate peak corresponding to formation of a bound state between a hole and a magnon. This polaron peak is separated from the simple magnon spectrum by energy proportional to single particle tunneling and can be easily resolved with currently available experimental techniques. For some momentum transfer there is an additional peak corresponding to photoexciting a bound state between two holes and a magnon. We point out that in two component Bose mixtures in triangular lattices one can also create dynamical magnetic polarons, with one hole and one magnon forming a repulsive bound state.

A Lattice Formulation of Weyl Fermions on a Single Curved Surface

Abstract In the standard lattice domain-wall fermion formulation, one needs two flat domain-walls where both of the left- and right-handed massless modes appear. In this work we investigate a single domain-wall system with a nontrivial curved background. Specifically we consider a massive fermion on a 3D square lattice, whose domain-wall is a 2D sphere. In the free theory, we find that a single Weyl fermion is localized at the wall and it feels gravity through the induced spin connection. With a topologically nontrivial U(1) link gauge field, however, we find a zero mode with the opposite chirality localized at the center where the gauge field is singular. In the latter case, the low-energy effective theory is not chiral but vectorlike. We discuss how to circumvent this obstacle in formulating lattice chiral gauge theory in the single domain-wall fermion system.

Adiabatic Evolution of Low-Temperature Many-Body Systems

We consider finite-range, many-body fermionic lattice models and we study the evolution of their thermal equilibrium state after introducing a weak and slowly varying time-dependent perturbation. Under suitable assumptions on the external driving, we derive a representation for the average of the evolution of local observables via a convergent expansion in the perturbation, for small enough temperatures. Convergence holds for a range of parameters that is uniform in the size of the system. Under a spectral gap assumption on the unperturbed Hamiltonian, convergence is also uniform in temperature. As an application, our expansion allows us to prove closeness of the time-evolved state to the instantaneous Gibbs state of the perturbed system, in the sense of expectation of local observables, at zero and at small temperatures. As a corollary, we also establish the validity of linear response. Our strategy is based on a rigorous version of the Wick rotation, which allows us to represent the Duhamel expansion for the real-time dynamics in terms of Euclidean correlation functions, for which precise decay estimates are proved using fermionic cluster expansion.

Tilted Dirac superconductor at quantum criticality: restoration of Lorentz symmetry

Lorentz symmetry appears as a quite robust feature of the strongly interacting Dirac materials even though the lattice interactions break such a symmetry. We here demonstrate that the Lorentz symmetry is restored at the quantum-critical point (QCP) separating the tilted Dirac semimetal, breaking this symmetry already at the noninteracting level, from a gapped s-wave superconducting instability. To this end, we employ a one-loop ϵ = (3 − D)-expansion close to the D = 3 upper critical dimension of the corresponding Gross-Neveu-Yukawa field theory. In particular, we show that the tilt parameter is irrelevant and ultimately vanishes at the QCP separating the two phases. In fact, as we argue here, such a Lorentz symmetry restoration may be generic for the strongly interacting tilted Dirac semimetals, irrespective of whether they feature mirror-symmetric or mirror-asymmetric tilting, and is also insensitive to whether the instability represents an insulator or a gapped superconductor. The proposed scenario can be tested in the quantum Monte Carlo simulations of the interacting tilted Dirac fermion lattice models.

DivERGe implements various Exact Renormalization Group examples

We present divERGe, an open source, high-performance C/C++/Python library for functional renormalization group (FRG) calculations on lattice fermions. The versatile model interface is tailored to real materials applications and seamlessly integrates with existing, standard tools from the ab-initio community. The code fully supports multi-site, multi-orbital, and non-SU(2) models in all of the three included FRG variants: TU^22FRG, N-patch FRG, and grid FRG. With this, the divERGe library paves the way for widespread application of FRG as a tool in the study of competing orders in quantum materials.

Codebase release 0.5 for divERGe

We present divERGe, an open source, high-performance C/C++/Python library for functional renormalization group (FRG) calculations on lattice fermions. The versatile model interface is tailored to real materials applications and seamlessly integrates with existing, standard tools from the ab-initio community. The code fully supports multi-site, multi-orbital, and non-SU(2) models in all of the three included FRG variants: TU^22FRG, N-patch FRG, and grid FRG. With this, the divERGe library paves the way for widespread application of FRG as a tool in the study of competing orders in quantum materials.

Interface fluctuations associated with split Fermi seas

We consider the asymptotic behaviour of a family of unidimensional lattice fermion models, which are in exact correspondence with certain probability laws on partitions and on unitary matrices. These models exhibit limit shapes, and in the case where the bulk of these shapes are described by analytic functions, the fluctuations around their interfaces have been shown to follow a universal Tracy–Widom distribution or its higher-order analogue. Non-differentiable bulk limit shape functions arise when a gap appears in some quantum numbers of the model, in other words when the Fermi sea is split. We show that split Fermi seas give rise to new interface fluctuations, governed by integer powers of universal distributions. This breakdown in universality is analogous to the behaviour of a random Hermitian matrix when the support of its limiting eigenvalue distribution has multiple cuts, with oscillations appearing in the limit of the two-point correlation function. We show that when the Fermi sea is split in the lattice fermion model, there are multiple cuts in the eigenvalue support of the corresponding unitary matrix model.

Impact of conditional modelling for a universal autoregressive quantum state

We present a generalized framework to adapt universal quantum state approximators, enabling them to satisfy rigorous normalization and autoregressive properties. We also introduce filters as analogues to convolutional layers in neural networks to incorporate translationally symmetrized correlations in arbitrary quantum states. By applying this framework to the Gaussian process state, we enforce autoregressive and/or filter properties, analyzing the impact of the resulting inductive biases on variational flexibility, symmetries, and conserved quantities. In doing so we bring together different autoregressive states under a unified framework for machine learning-inspired ansätze. Our results provide insights into how the autoregressive construction influences the ability of a variational model to describe correlations in spin and fermionic lattice models, as well as ab initio electronic structure problems where the choice of representation affects accuracy. We conclude that, while enabling efficient and direct sampling, thus avoiding autocorrelation and loss of ergodicity issues in Metropolis sampling, the autoregressive construction materially constrains the expressivity of the model in many systems.

Dynamical chirality production in one dimension

We discuss the quantum computation of dynamical chirality production in lattice gauge theory. Although the chirality of a lattice fermion is complicated in general dimensions, it can be simply formulated on a one-dimensional lattice. The chiral fermion formalism enables us to extract the physical part of the chirality production that would be interpreted as the chiral anomaly in the continuous theory. We demonstrate the computation of the Z2 lattice gauge theory on a classical emulator. Published by the American Physical Society 2024

Calculation of cSW at one-loop order for Brillouin fermions

The Brillouin action is a Wilson-like lattice fermion action with a 81-point stencil, which was found to ameliorate the Wilson action in many respects. The Sheikholeslami-Wohlert coefficient cSW of the clover improvement term has a perturbative expansion cSW=cSW(0)+g02cSW(1)+O(g04). At tree-level cSW(0)=r holds for Wilson and Brillouin fermions alike. We present the Feynman rules for the Brillouin action in lattice perturbation theory, and employ them to calculate the one-loop coefficient cSW(1) with plaquette or Lüscher-Weisz gluons. Numerically its value is found to be about half that of the Wilson action. Published by the American Physical Society 2024