Fermion Model Research Articles

Quantum leap: observing antiferromagnetic transition in a 3D fermionic Hubbard model with ultracold atoms

Quantum leap: observing antiferromagnetic transition in a 3D fermionic Hubbard model with ultracold atoms

Searching accretion-enhanced dark matter annihilation signals in the Galactic Centre

This study reanalyzes the detection prospects of dark matter (DM) annihilation signals in the Galactic Center, focusing on velocity-dependent dynamics within a spike density near the supermassive black hole (Sgr A⋆). We investigate three annihilation processes — p-wave, resonance, and forbidden annihilation — under semi-relativistic velocities, leveraging gamma-ray data from Fermi and DAMPE telescopes. Our analysis integrates a fermionic DM model with an electroweak axion-like particle (ALP) portal, exploring annihilation into two or four photons. Employing a comprehensive six-dimensional integration, we precisely calculate DM-induced gamma-ray fluxes near Sgr A⋆, incorporating velocity and positional dependencies in the annihilation cross-section and photon yield spectra. Our findings highlight scenarios of resonance and forbidden annihilation, where the larger ALP-DM-DM coupling constant Caχχ can affect spike density, potentially yielding detectable gamma-ray line spectra within Fermi and DAMPE energy resolution. We set upper limits for Caχχ across these scenarios, offering insights into the detectability and spectral characteristics of DM annihilation signals from the Galactic Center.

Quantum statistical mechanics of the Sachdev-Ye-Kitaev model and charged black holes

This review is a contribution to a book dedicated to the memory of Michael E. Fisher. The first example of a quantum many body system not expected to have any quasiparticle excitations was the Wilson-Fisher conformal field theory. The absence of quasiparticles can be established in the compressible, metallic state of the Sachdev-Ye-Kitaev model of fermions with random interactions. The solvability of the latter model has enabled numerous computations of the non-quasiparticle dynamics of chaotic many-body states, such as those expected to describe quantum black holes. We review thermodynamic properties of the SYK model, and describe how they have led to an understanding of the universal structure of the low energy density of states of charged black holes without low energy supersymmetry.

Generalization to d-dimensions of a fermionic path integral for exact enumeration of polygons on hypercubic lattices

The generating function for polygons on the square lattice has been known for many decades and is closely related to the path integral formulation of a free fermion model. On the cubic and hypercubic lattices the generating function is still unknown and the problem remains open. It has been conjectured that the three-dimensional (3D) and higher dimensional problems are not solvable—or, to be more precise, that there are no differentiably finite (D-finite) solutions. In this context, very recently a Berezin integral of an exponentiated Grassmann action was found for the polygon generating function on the cubic lattice, making explicit the connection between 3D polygons and a model of interacting fermions. Here we address the problem of how to generalize the 3D result to higher dimensions. We derive a Grassmann representation in terms of a Berezin integral for the generating function of polygons on d-dimensional hypercubic lattices. On the one hand, this new result admittedly brings us no closer to the problem of finding an explicit analytic expression for the desired generating function for polygons. On the other hand, however, the significant advance reported here precisely quantifies the remarkable mathematical difficulty of finding the explicit generating function. Indeed, the non-quadratic functional form of the Grassmann action that we derive here provides a very clear picture of the formidable mathematical obstruction that would need to be overcome. Specifically, in d dimensions, the Grassmann action contains terms of degree 2(d-1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$2(d-1)$$\\end{document}, so the model describes interacting rather than free fermions. It is an open problem whether or not these models of interacting fermions can in principle be free fermionized through some still undiscovered algebraic method, but it is widely believed that this goal is mathematically impossible.

Explanation of puzzling FQHE at the filling fraction 3/4 in a band-hole 2D system in GaAs

A recent experiment revealed an unexpected FQHE at filling fraction 3/4 in a GaAs 2D hole system, which contradicts the composite fermion model prediction and the observation of a compressible Hall metal-type state in a twin 2D electron system in GaAs at the same filling fraction 3/4 at almost same other conditions. This finding challenges conventional effective single-quasiparticle model for FQHE exposing its limitations. We explain this experimental observation within a multiparticle approach based on a topological cyclotron commensurability criterion. This allows to generalize Laughlin function for filling fractions from the complete FQHE hierarchy including observable FQHE states at even denominator fractions. The topological multiparticle approach helps to decipher a structure of composite fermions and provides their generalization for so-called enigmatic states including even denominator filling fractions, and also for quantum fractional Hall-type behavior in Chern topological insulators without a magnetic field.

Prerelaxation in quantum, classical, and quantum-classical two-impurity models

We numerically study the relaxation dynamics of impurity-host systems, focusing on the presence of long-lived metastable states in the nonequilibrium dynamics after an initial excitation of the impurities. In generic systems, an excited impurity coupled to a large bath at zero temperature is expected to relax and approach its ground state over time. However, certain exceptional cases exhibit metastability, where the system remains in an excited state on timescales largely exceeding the typical relaxation time. We study this phenomenon for three prototypical impurity models: a tight-binding quantum model of independent spinless fermions on a lattice with two stub impurities, a classical-spin Heisenberg model with two weakly coupled classical impurity spins, and a tight-binding quantum model of independent electrons with two classical impurity spins. Through numerical integration of the fundamental equations of motion, we find that all three models exhibit similar qualitative behavior: complete relaxation for nearest-neighbor impurities and incomplete or strongly delayed relaxation for next-nearest-neighbor impurities. The underlying mechanisms leading to this behavior differ between models and include impurity-induced bound states, emergent approximately conserved local observables, and exact cancellation of local and nonlocal dissipation effects. Published by the American Physical Society 2024

A quintessential quantum simulator takes a 10 000-fold leap

Experiments on the fermionic Hubbard model can now be made much larger, more uniform, and more quantitative.

Jordan–Wigner transformation constructed for spinful fermions at S=1/2 spins in one dimension

ABSTRACT An exact Jordan–Wigner type of transformation is presented in 1D connecting spin-1/2 operators to spinful canonical Fermi operators. The transformation contains two free parameters allowing a broad interconnection possibility between spin models and fermionic models containing spinful Fermi operators.

Expanding the reach of quantum optimization with fermionic embeddings

Quadratic programming over orthogonal matrices encompasses a broad class of hard optimization problems that do not have an efficient quantum representation. Such problems are instances of the little noncommutative Grothendieck problem (LNCG), a generalization of binary quadratic programs to continuous, noncommutative variables. In this work, we establish a natural embedding for this class of LNCG problems onto a fermionic Hamiltonian, thereby enabling the study of this classical problem with the tools of quantum information. This embedding is accomplished by a new representation of orthogonal matrices as fermionic quantum states, which we achieve through the well-known double covering of the orthogonal group. Correspondingly, the embedded LNCG Hamiltonian is a two-body fermion model. Determining extremal states of this Hamiltonian provides an outer approximation to the original problem, a quantum analogue to classical semidefinite relaxations. In particular, when optimizing over the special orthogonal group our quantum relaxation obeys additional, powerful constraints based on the convex hull of rotation matrices. The classical size of this convex-hull representation is exponential in matrix dimension, whereas our quantum representation requires only a linear number of qubits. Finally, to project the relaxed solution back into the feasible space, we propose rounding procedures which return orthogonal matrices from appropriate measurements of the quantum state. Through numerical experiments we provide evidence that this rounded quantum relaxation can produce high-quality approximations.

Fixed lines in four fermion models in two dimensions

Motivated by conjectures about near-horizon dynamics in quantum gravity, we search for lines of perturbatively accessible fixed points emanating from models of N free fermions. Through two loops we find a new class of models, apart from the well-known Abelian Thirring models. Further study is needed to see whether these can lead to true conformal manifolds, or perhaps a new class of large-N fixed points. Published by the American Physical Society 2024

De Sitter at all loops: the story of the Schwinger model

We consider the two-dimensional Schwinger model of a massless charged fermion coupled to an Abelian gauge field on a fixed de Sitter background. The theory admits an exact solution, first examined by Jayewardena, and can be analyzed efficiently using Euclidean methods. We calculate fully non-perturbative, gauge-invariant correlation functions of the electric field as well as the fermion and analyze these correlators in the late-time limit. We compare these results with the perturbative picture, for example by verifying that the one-loop contribution to the fermion two-point function, as predicted from the exact solution, matches the direct computation of the one-loop Feynman diagram. We demonstrate many features endemic of quantum field theory in de Sitter space, including the appearance of late-time logarithms, their resummation to de Sitter invariant expressions, and Boltzmann suppressed non-perturbative phenomena, with surprising late-time features.

Top-down and bottom-up: Studying the SMEFT beyond leading order in $1/\Lambda^2$

In order to assess the relevance of higher order terms in the Standard Model effective field theory (SMEFT) expansion we consider four new physics models and their impact on the Drell Yan cross section. Of these four, one scalar model has no effect on Drell Yan, a model of fermions while appearing to generate a momentum expansion actually belongs to the vacuum expectation value expansion and so has a nominal effect on the process. The remaining two, a leptoquark and a Z'Z′ model exhibit a momentum expansion. After matching these models to dimension-ten we study the how the inclusion of dimension-eight and dimension-ten operators in hypothetical effective field theory fits to the full ultraviolet models impacts fits. We do this both in the top-down approach, and in a very limited approximation to the bottom up approach of the SMEFT to infer the impact of a fully general fit to the SMEFT. We find that for the more weakly coupled models a strictly dimension-six fit is sufficient. In contrast when stronger interactions or lighter masses are considered the inclusion of dimension-eight operators becomes necessary. However, their Wilson coefficients perform the role of nuisance parameters with best fit values which can differ statistically from the theory prediction. In the most strongly coupled theories considered (which are already ruled out by data) the inclusion of dimension-ten operators allows for the measurement of dimension-eight operator coefficients consistent with theory predictions and the dimension-ten operator coefficients then behave as nuisance parameters. We also study the impact of the inclusion of partial next order results, such as dimension-six squared contributions, and find that in some cases they improve the convergence of the series while in others they hinder it.

Quantum advantage and stability to errors in analogue quantum simulators

Several quantum hardware platforms, while being unable to perform fully fault-tolerant quantum computation, can still be operated as analogue quantum simulators for addressing many-body problems. However, due to the presence of errors, it is not clear to what extent those devices can provide us with an advantage with respect to classical computers. In this work, we make progress on this problem for noisy analogue quantum simulators computing physically relevant properties of many-body systems both in equilibrium and undergoing dynamics. We first formulate a system-size independent notion of stability against extensive errors, which we prove for Gaussian fermion models, as well as for a restricted class of spin systems. Remarkably, for the Gaussian fermion models, our analysis shows the stability of critical models which have long-range correlations. Furthermore, we analyze how this stability may lead to a quantum advantage, for the problem of computing the thermodynamic limit of many-body models, in the presence of a constant error rate and without any explicit error correction.

Survival Probability, Particle Imbalance, and Their Relationship in Quadratic Models.

We argue that the dynamics of particle imbalance in quadratic fermionic models is, for the majority of initial many-body product states in the site occupation basis, virtually indistinguishable from the dynamics of survival probabilities of single-particle states. We then generalize our statement to a similar relationship between the non-equal time and space density correlation functions in many-body states, and the transition probabilities of single-particle states at nonzero distances. Finally, we study the equal-time connected density-density correlation functions in many-body states, which exhibit certain qualitative analogies with the survival and transition probabilities of single-particle states. Our results are numerically tested for two paradigmatic models of single-particle localization: the 3D Anderson model and the 1D Aubry-André model. This work gives an affirmative answer to the question of whether it is possible to measure features of single-particle survival and transition probabilities by the dynamics of observables in many-body states.

Fermionic UV models for neutral triple gauge boson vertices

Searches for anomalous neutral triple gauge boson couplings (NTGCs) provide important tests for the gauge structure of the standard model. In SMEFT (“standard model effective field theory”) NTGCs appear only at the level of dimension-8 operators. While the phenomenology of these operators has been discussed extensively in the literature, renormalizable UV models that can generate these operators are scarce. In this work, we study a variety of extensions of the SM with heavy fermions and calculate their matching to d = 8 NTGC operators. We point out that the complete matching of UV models requires four different CP-conserving d = 8 operators and that the single CPC d = 8 operator, most commonly used by the experimental collaborations, does not describe all possible NTGC form factors. Despite stringent experimental constraints on NTGCs, limits on the scale of UV models are relatively weak, because their contributions are doubly suppressed (being d = 8 and 1-loop). We suggest a series of benchmark UV scenarios suitable for interpreting searches for NTGCs in the upcoming LHC runs, obtain their current limits and provide estimates for the expected sensitivity of the high-luminosity LHC.

Probing Inert Triplet Model at a multi-TeV muon collider via vector boson fusion with forward muon tagging

This study investigates the potential of a multi-TeV Muon Collider (MuC) for probing the Inert Triplet Model (ITM), which introduces a triplet scalar field with hypercharge Y = 0 to the Standard Model. The ITM stands out as a compelling Beyond the Standard Model scenario, featuring a neutral triplet T0 and charged triplets T±. Notably, T0 is posited as a dark matter (DM) candidate, being odd under a Z2 symmetry. Rigorous evaluations against theoretical, collider, and DM experimental constraints corner the triplet scalar mass to a narrow TeV-scale region, within which three benchmark points are identified, with T± masses of 1.21 TeV, 1.68 TeV, and 3.86 TeV, for the collider study. The ITM’s unique TTVV four-point vertex, differing from fermionic DM models, facilitates efficient pair production through Vector Boson Fusion (VBF). This characteristic positions the MuC as an ideal platform for exploring the ITM, particularly due to the enhanced VBF cross-sections at high collision energies. To address the challenge of the soft decay products of T± resulting from the narrow mass gap between T± and T0, we propose using Disappearing Charged Tracks (DCTs) from T± and Forward muons as key signatures. We provide event counts for these signatures at MuC energies of 6 TeV and 10 TeV, with respective luminosities of 4 ab−1 and 10 ab−1. Despite the challenge of beam-induced backgrounds contaminating the signal, we demonstrate that our proposed final states enable the MuC to achieve a 5σ discovery for the identified benchmark points, particularly highlighting the effectiveness of the final state with one DCT and one Forward muon.

Quantum phase transitions in one-dimensional nanostructures: a comparison between DFT and DMRG methodologies.

In the realm of quantum chemistry, the accurate prediction of electronic structure and properties of nanostructures remains a formidable challenge. Density functional theory (DFT) and density matrix renormalization group (DMRG) have emerged as two powerful computational methods for addressing electronic correlation effects in diverse molecular systems. We compare ground-state energies ( ), density profiles ( ), and average entanglement entropies ( ) in metals, insulators and at the transition from metal to insulator, in homogeneous, superlattices, and harmonically confined chains described by the fermionic one-dimensional Hubbard model. While for the homogeneous systems, there is a clear hierarchy between the deviations, , and all the deviations decrease with the chain size; for superlattices and harmonic confinement, the relation among the deviations is less trivial and strongly dependent on the superlattice structure and the confinement strength considered. For the superlattices, in general, increasing the number of impurities in the unit cell represents lower precision in the DFT calculations. For the confined chains, DFT performs better for metallic phases, while the highest deviations appear for the Mott and band-insulator phases. This work provides a comprehensive comparative analysis of these methodologies, shedding light on their respective strengths, limitations, and applications. The DFT calculations were performed using the standard Kohn-Sham scheme within the BALDA approach. It integrated the numerical Bethe-Ansatz (BA) solution of the Hubbard model as the homogeneous density functional within a local-density approximation (LDA) for the exchange-correlation energy. The DMRG algorithms were implemented using the ITensor library, which is based on the matrix product states (MPS) ansatz. The calculations were performed until the energy reaches convergence of at least .

Antiferromagnetic phase transition in a 3D fermionic Hubbard model.

The fermionic Hubbard model (FHM)1 describes a wide range of physical phenomena resulting from strong electron-electron correlations, including conjectured mechanisms for unconventional superconductivity. Resolving its low-temperature physics is, however, challenging theoretically or numerically. Ultracold fermions in optical lattices2,3 provide a clean and well-controlled platform offering a path to simulate the FHM. Doping the antiferromagnetic ground state of a FHM simulator at half-filling is expected to yield various exotic phases, including stripe order4, pseudogap5, and d-wave superfluid6, offering valuable insights into high-temperature superconductivity7-9. Although the observation of antiferromagnetic correlations over short10 and extended distances11 has been obtained, the antiferromagnetic phase has yet to be realized as it requires sufficiently low temperatures in a large and uniform quantum simulator. Here we report the observation of the antiferromagnetic phase transition in a three-dimensional fermionic Hubbard system comprising lithium-6 atoms in a uniform optical lattice with approximately 800,000 sites. When the interaction strength, temperature and doping concentration are finely tuned to approach their respective critical values, a sharp increase in the spin structure factor is observed. These observations can be well described by a power-law divergence, with a critical exponent of 1.396 from the Heisenberg universality class12. At half-filling and with optimal interaction strength, the measured spin structure factor reaches 123(8), signifying the establishment of an antiferromagnetic phase. Our results provide opportunities for exploring the low-temperature phase diagram of the FHM.

Induced motions on Carroll geometries

In this article, we consider some Carrollian dynamical systems as effective models on null hypersurfaces in a Lorentzian spacetime. We show that we can realize Carroll models from more usual ‘relativistic’ theories. In particular, we show how ambient null geodesics imply the classical ʼno Carroll motion’ and, more interestingly, we find that the ambient model of chiral fermions implies Hall motion on null hypersurfaces, in agreement with previous intrinsic Carroll results. We also show how Wigner–Souriau translations imply (apparent) Carroll motion, and how ambient particles with a non vanishing gyromagnetic ratio cannot have a Carrollian description.

Non-unitary Trotter circuits for imaginary time evolution

We propose an imaginary time equivalent of the well-established Pauli gadget primitive for Trotter-decomposed real time evolution, using mid-circuit measurements on a single ancilla qubit. Imaginary time evolution (ITE) is widely used for obtaining the ground state (GS) of a system on classical hardware, computing thermal averages, and as a component of quantum algorithms that perform non-unitary evolution. Near-term implementations on quantum hardware rely on heuristics, compromising their accuracy. As a result, there is growing interest in the development of more natively quantum algorithms. Since it is not possible to implement a non-unitary gate deterministically, we resort to the implementation of probabilistic ITE (PITE) algorithms, which rely on a unitary quantum circuit to simulate a block encoding of the ITE operator—that is, they rely on successful ancillary measurements to evolve the system non-unitarily. Compared with previous PITE proposals, the suggested block encoding in this paper results in shorter circuits and is simpler to implement, requiring only a slight modification of the Pauli gadget primitive. This scheme was tested on the transverse Ising model and the fermionic Hubbard model and is demonstrated to converge to the GS of the system.