![2D-Hubbard-model](https://upload.wikimedia.org/wikipedia/commons/a/a8/2D-Hubbard-model.png)
The Hubbard model is an
approximate model used to describe the transition between
conducting
Conducting is the art of directing a musical performance, such as an orchestral or choral concert. It has been defined as "the art of directing the simultaneous performance of several players or singers by the use of gesture." The primary duti ...
and
insulating systems.
[
] It is particularly useful in
solid-state physics
Solid-state physics is the study of rigid matter, or solids, through methods such as quantum mechanics, crystallography, electromagnetism, and metallurgy. It is the largest branch of condensed matter physics. Solid-state physics studies how the l ...
. The model is named for
John Hubbard.
The Hubbard model states that each electron experiences competing forces: one pushes it to tunnel to neighboring atoms, while the other pushes it away from its neighbors.
Its
Hamiltonian
Hamiltonian may refer to:
* Hamiltonian mechanics, a function that represents the total energy of a system
* Hamiltonian (quantum mechanics), an operator corresponding to the total energy of that system
** Dyall Hamiltonian, a modified Hamiltonian ...
thus has two terms: a kinetic term allowing for
tunneling ("hopping") of particles between lattice sites and a potential term reflecting on-site interaction. The particles can either be
fermions
In particle physics, a fermion is a particle that follows Fermi–Dirac statistics. Generally, it has a half-odd-integer spin: spin , spin , etc. In addition, these particles obey the Pauli exclusion principle. Fermions include all quarks and ...
, as in Hubbard's original work, or
bosons
In particle physics, a boson ( ) is a subatomic particle whose spin quantum number has an integer value (0,1,2 ...). Bosons form one of the two fundamental classes of subatomic particle, the other being fermions, which have odd half-integer ...
, in which case the model is referred to as the "
Bose–Hubbard model
The Bose–Hubbard model gives a description of the physics of interacting spinless bosons on a lattice. It is closely related to the Hubbard model that originated in solid-state physics as an approximate description of superconducting systems an ...
".
The Hubbard model is a useful approximation for particles in a periodic potential at sufficiently low temperatures, where all the particles may be assumed to be in the lowest
Bloch band, and long-range interactions between the particles can be ignored. If interactions between particles at different sites of the lattice are included, the model is often referred to as the "extended Hubbard model". In particular, the Hubbard term, most commonly denoted by ''U'', is applied in first principles based simulations using
Density Functional Theory
Density-functional theory (DFT) is a computational quantum mechanical modelling method used in physics, chemistry and materials science to investigate the electronic structure (or nuclear structure) (principally the ground state) of many-body ...
, DFT. The inclusion of the Hubbard term in DFT simulations is important as this improves the prediction of electron localisation and thus it prevents the incorrect prediction of metallic conduction in insulating systems.
The Hubbard model introduces short-range interactions between electrons to the
tight-binding model
In solid-state physics, the tight-binding model (or TB model) is an approach to the calculation of electronic band structure using an approximate set of wave functions based upon Quantum superposition, superposition of wave functions for isolated ...
, which only includes kinetic energy (a "hopping" term) and interactions with the atoms of the lattice (an "atomic" potential). When the interaction between electrons is strong, the behavior of the Hubbard model can be qualitatively different from a tight-binding model. For example, the Hubbard model correctly predicts the existence of
Mott insulator
Mott insulators are a class of materials that are expected to conduct electricity according to conventional band theories, but turn out to be insulators (particularly at low temperatures). These insulators fail to be correctly described by band ...
s: materials that are insulating due to the strong repulsion between electrons, even though they satisfy the usual criteria for conductors, such as having an odd number of electrons per unit cell.
History
The model was originally proposed in 1963 to describe electrons in solids. Hubbard,
Martin Gutzwiller Martin may refer to:
Places
* Martin City (disambiguation)
* Martin County (disambiguation)
* Martin Township (disambiguation)
Antarctica
* Martin Peninsula, Marie Byrd Land
* Port Martin, Adelie Land
* Point Martin, South Orkney Islands
Austral ...
and Junjiro Kanamori each independently proposed it.
Since then, it has been applied to the study of
high-temperature superconductivity
High-temperature superconductors (abbreviated high-c or HTS) are defined as materials that behave as superconductors at temperatures above , the boiling point of liquid nitrogen. The adjective "high temperature" is only in respect to previ ...
, quantum magnetism, and charge density waves.
Narrow energy band theory
The Hubbard model is based on the
tight-binding
In solid-state physics, the tight-binding model (or TB model) is an approach to the calculation of electronic band structure using an approximate set of wave functions based upon superposition of wave functions for isolated atoms located at each ...
approximation from
solid-state physics
Solid-state physics is the study of rigid matter, or solids, through methods such as quantum mechanics, crystallography, electromagnetism, and metallurgy. It is the largest branch of condensed matter physics. Solid-state physics studies how the l ...
, which describes particles moving in a periodic potential, typically referred to as a
lattice
Lattice may refer to:
Arts and design
* Latticework, an ornamental criss-crossed framework, an arrangement of crossing laths or other thin strips of material
* Lattice (music), an organized grid model of pitch ratios
* Lattice (pastry), an orna ...
. For real materials, each lattice site might correspond with an ionic core, and the particles would be the valence electrons of these ions. In the tight-binding approximation, the Hamiltonian is written in terms of
Wannier states, which are localized states centered on each lattice site. Wannier states on neighboring lattice sites are coupled, allowing particles on one site to "hop" to another. Mathematically, the strength of this coupling is given by a "hopping integral", or "transfer integral", between nearby sites. The system is said to be in the tight-binding limit when the strength of the hopping integrals falls off rapidly with distance. This coupling allows states associated with each lattice site to hybridize, and the
eigenstates
In quantum physics, a quantum state is a mathematical entity that provides a probability distribution for the outcomes of each possible measurement on a system. Knowledge of the quantum state together with the rules for the system's evolution i ...
of such a
crystal
A crystal or crystalline solid is a solid material whose constituents (such as atoms, molecules, or ions) are arranged in a highly ordered microscopic structure, forming a crystal lattice that extends in all directions. In addition, macro ...
line system are
Bloch's functions, with the energy levels divided into separated
energy bands. The width of the bands depends upon the value of the hopping integral.
The Hubbard model introduces a contact interaction between particles of opposite spin on each site of the lattice. When the Hubbard model is used to describe electron systems, these interactions are expected to be repulsive, stemming from the
screened Coulomb interaction. However, attractive interactions have also been frequently considered. The physics of the Hubbard model is determined by competition between the strength of the hopping integral, which characterizes the system's
kinetic energy
In physics, the kinetic energy of an object is the energy that it possesses due to its motion.
It is defined as the work needed to accelerate a body of a given mass from rest to its stated velocity. Having gained this energy during its accele ...
, and the strength of the interaction term. The Hubbard model can therefore explain the transition from metal to insulator in certain interacting systems. For example, it has been used to describe
metal
A metal (from Greek μέταλλον ''métallon'', "mine, quarry, metal") is a material that, when freshly prepared, polished, or fractured, shows a lustrous appearance, and conducts electricity and heat relatively well. Metals are typicall ...
oxides as they are heated, where the corresponding increase in nearest-neighbor spacing reduces the hopping integral to the point where the on-site potential is dominant. Similarly, the Hubbard model can explain the transition from conductor to insulator in systems such as
rare-earth
The rare-earth elements (REE), also called the rare-earth metals or (in context) rare-earth oxides or sometimes the lanthanides (yttrium and scandium are usually included as rare earths), are a set of 17 nearly-indistinguishable lustrous silve ...
pyrochlore
Pyrochlore () is a mineral group of the niobium end member of the pyrochlore supergroup.
The general formula, (where A and B are metals), represent a family of phases isostructural to the mineral pyrochlore.
Pyrochlores are an important class of ...
s as the
atomic number
The atomic number or nuclear charge number (symbol ''Z'') of a chemical element is the charge number of an atomic nucleus. For ordinary nuclei, this is equal to the proton number (''n''p) or the number of protons found in the nucleus of every ...
of the rare-earth metal increases, because the
lattice parameter
A lattice constant or lattice parameter is one of the physical dimensions and angles that determine the geometry of the unit cells in a crystal lattice, and is proportional to the distance between atoms in the crystal. A simple cubic crystal has o ...
increases (or the angle between atoms can also change) as the rare-earth element atomic number increases, thus changing the relative importance of the hopping integral compared to the on-site repulsion.
Example: one dimensional hydrogen atom chain
The
hydrogen atom has one electron, in the so-called ''s'' orbital, which can either be spin up (
) or spin down (
). This orbital can be occupied by at most two electrons, one with
spin up and one down (see
Pauli exclusion principle
In quantum mechanics, the Pauli exclusion principle states that two or more identical particles with half-integer spins (i.e. fermions) cannot occupy the same quantum state within a quantum system simultaneously. This principle was formulated ...
).
Under
band theory, for a 1D chain of hydrogen atoms, the 1s orbital forms a continuous band, which would be exactly half-full. The 1D chain of hydrogen atoms is thus predicted to be a conductor under conventional band theory. This 1D string is the only configuration simple enough to be solved directly.
But in the case where the spacing between the hydrogen atoms is gradually increased, at some point the chain must become an insulator.
Expressed using the Hubbard model, the Hamiltonian is made up of two terms. The first term describes the kinetic energy of the system, parameterized by the hopping integral,
. The second term is the on-site interaction of strength
that represents the electron repulsion. Written out in
second quantization
Second quantization, also referred to as occupation number representation, is a formalism used to describe and analyze quantum many-body systems. In quantum field theory, it is known as canonical quantization, in which the fields (typically as ...
notation, the Hubbard
Hamiltonian
Hamiltonian may refer to:
* Hamiltonian mechanics, a function that represents the total energy of a system
* Hamiltonian (quantum mechanics), an operator corresponding to the total energy of that system
** Dyall Hamiltonian, a modified Hamiltonian ...
then takes the form
:
where
is the spin-density operator for spin
on the
-th site. The density operator is
and occupation of
-th site for the wavefunction
is
. Typically ''t'' is taken to be positive, and ''U'' may be either positive or negative, but is assumed to be positive when considering electronic systems.
Without the contribution of the second term, the Hamiltonian resolves to the
tight binding
In solid-state physics, the tight-binding model (or TB model) is an approach to the calculation of electronic band structure using an approximate set of wave functions based upon superposition of wave functions for isolated atoms located at each ...
formula from regular band theory.
Including the second term yields a realistic model that also predicts a transition from conductor to insulator as the ratio of interaction to hopping,
, is varied. This ratio can be modified by, for example, increasing the inter-atomic spacing, which would decrease the magnitude of
without affecting
. In the limit where
, the chain simply resolves into a set of isolated
magnetic moment
In electromagnetism, the magnetic moment is the magnetic strength and orientation of a magnet or other object that produces a magnetic field. Examples of objects that have magnetic moments include loops of electric current (such as electromagne ...
s. If
is not too large, the overlap integral provides for
superexchange interactions between neighboring magnetic moments, which may lead to a variety of interesting magnetic correlations, such as ferromagnetic, antiferromagnetic, etc. depending on the model parameters. The one-dimensional Hubbard model was solved by
Lieb and Wu using the
Bethe ansatz. Essential progress was achieved in the 1990s: a
hidden symmetry was discovered, and the
scattering matrix
In physics, the ''S''-matrix or scattering matrix relates the initial state and the final state of a physical system undergoing a scattering process. It is used in quantum mechanics, scattering theory and quantum field theory (QFT).
More forma ...
,
correlation functions,
thermodynamic
Thermodynamics is a branch of physics that deals with heat, work, and temperature, and their relation to energy, entropy, and the physical properties of matter and radiation. The behavior of these quantities is governed by the four laws of the ...
and
quantum entanglement
Quantum entanglement is the phenomenon that occurs when a group of particles are generated, interact, or share spatial proximity in a way such that the quantum state of each particle of the group cannot be described independently of the state of ...
were evaluated.
[
]
More complex systems
Although Hubbard is useful in describing systems such as a 1D chain of hydrogen atoms, it is important to note that more complex systems may experience other effects that the Hubbard model does not consider. In general, insulators can be divided into
Mott–Hubbard insulators and
charge-transfer insulators.
A Mott–Hubbard insulator can be described as
:
This can be seen as analogous to the Hubbard model for hydrogen chains, where conduction between unit cells can be described by a transfer integral.
However, it is possible for the electrons to exhibit another kind of behavior:
:
This is known as charge transfer and results in
charge-transfer insulators. Unlike Mott–Hubbard insulators electron transfer happens only within a unit cell.
Both of these effects may be present and compete in complex ionic systems.
Numerical treatment
The fact that the Hubbard model has not been solved analytically in arbitrary dimensions has led to intense research into numerical methods for these strongly correlated electron systems.
[
][
] One major goal of this research is to determine the low-temperature phase diagram of this model, particularly in two-dimensions. Approximate numerical treatment of the Hubbard model on finite systems is possible via various methods.
One such method, the
Lanczos algorithm
The Lanczos algorithm is an iterative method devised by Cornelius Lanczos that is an adaptation of power methods to find the m "most useful" (tending towards extreme highest/lowest) eigenvalues and eigenvectors of an n \times n Hermitian matri ...
, can produce static and dynamic properties of the system. Ground state calculations using this method require the storage of three vectors of the size of the number of states. The number of states scales exponentially with the size of the system, which limits the number of sites in the lattice to about 20 on 21st century hardware. With projector and finite-temperature
auxiliary-field Monte Carlo
Auxiliary-field Monte Carlo is a method that allows the calculation, by use of Monte Carlo techniques, of averages of operators in many-body quantum mechanical (Blankenbecler 1981, Ceperley 1977) or classical problems (Baeurle 2004, Baeurle 2003, ...
, two statistical methods exist that can obtain certain properties of the system. For low temperatures, convergence problems appear that lead to an exponential computational effort with decreasing temperature due to the so-called fermion
sign problem In applied mathematics, the numerical sign problem is the problem of numerically evaluating the integral of a highly oscillatory function of a large number of variables. Numerical methods fail because of the near-cancellation of the positive and ne ...
.
The Hubbard model can be studied within
dynamical mean-field theory Dynamical mean-field theory (DMFT) is a method to determine the electronic structure of strongly correlated materials. In such materials, the approximation of independent electrons, which is used in density functional theory and usual band structure ...
(DMFT). This scheme maps the Hubbard Hamiltonian onto a
single-site impurity model, a mapping that is formally exact only in infinite dimensions and in finite dimensions corresponds to the exact treatment of all purely local correlations only. DMFT allows one to compute the local
Green's function
In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions.
This means that if \operatorname is the linear differenti ...
of the Hubbard model for a given
and a given temperature. Within DMFT, the evolution of the
spectral function
The power spectrum S_(f) of a time series x(t) describes the distribution of power into frequency components composing that signal. According to Fourier analysis, any physical signal can be decomposed into a number of discrete frequencies, ...
can be computed and the appearance of the upper and lower Hubbard bands can be observed as correlations increase.
Simulator
Stacks of heterogeneous 2-dimensional
transition metal dichalcogenides (TMD) have been used to simulate geometries in more than one dimension.
Tungsten diselenide
Tungsten diselenide is an inorganic compound with the formula WSe2. The compound adopts a hexagonal crystalline structure similar to molybdenum disulfide. Every tungsten atom is covalently bonded to six selenium ligands in a trigonal prismatic c ...
and tungsten sulfide were stacked. This created a
moiré superlattice consisting of hexagonal
supercells
A supercell is a thunderstorm characterized by the presence of a mesocyclone: a deep, persistently rotating vertical draft, updraft. Due to this, these storms are sometimes referred to as rotating thunderstorms. Of the four classifications of ...
(repetition units defined by the relationship of the two materials). Each supercell then behaves as though it were a single atom. The distance between supercells is roughly 100x that of the atoms within them. This larger distance drastically reduces electron tunneling across supercells.
They can be used to form
Wigner crystal
A Wigner crystal is the solid (crystalline) phase of electrons first predicted by Eugene Wigner in 1934. A gas of electrons moving in a uniform, inert, neutralizing background (i.e. Jellium Model) will crystallize and form a lattice if the electr ...
s.
Electrodes
An electrode is an electrical conductor used to make contact with a nonmetallic part of a circuit (e.g. a semiconductor, an electrolyte, a vacuum or air). Electrodes are essential parts of batteries that can consist of a variety of materials dep ...
can be attached to regulate an
electric field. The electric field controls how many electrons fill each supercell. The number of electrons per supercell effectively determines which "atom" the lattice simulates. One electron/cell behaves like hydrogen, two/cell like helium, etc. As of 2022, supercells with up to eight electrons (
oxygen
Oxygen is the chemical element with the symbol O and atomic number 8. It is a member of the chalcogen group in the periodic table, a highly reactive nonmetal, and an oxidizing agent that readily forms oxides with most elements as wel ...
) could be simulated. One result of the simulation showed that the difference between metal and insulator is a continuous function of the electric field strength.
A "backwards" stacking regime allows the creation of a Chern insulator via the
anomalous quantum Hall effect (with the edges of the device acting as a conductor while the interior acted as an insulator.) The device functioned at a temperature of 5
Kelvin
The kelvin, symbol K, is the primary unit of temperature in the International System of Units (SI), used alongside its prefixed forms and the degree Celsius. It is named after the Belfast-born and University of Glasgow-based engineer and phy ...
s, far above the temperature at which the effect had first been observed.
See also
*
Anderson impurity model
The Anderson impurity model, named after Philip Warren Anderson, is a Hamiltonian that is used to describe magnetic impurities embedded in metals.
It is often applied to the description of Kondo effect-type problems, such as heavy fermion sys ...
*
Bloch's theorem
In condensed matter physics, Bloch's theorem states that solutions to the Schrödinger equation in a periodic potential take the form of a plane wave modulated by a periodic function. The theorem is named after the physicist Felix Bloch, who d ...
*
Electronic band structure
In solid-state physics, the electronic band structure (or simply band structure) of a solid describes the range of energy levels that electrons may have within it, as well as the ranges of energy that they may not have (called ''band gaps'' or ' ...
*
Solid-state physics
Solid-state physics is the study of rigid matter, or solids, through methods such as quantum mechanics, crystallography, electromagnetism, and metallurgy. It is the largest branch of condensed matter physics. Solid-state physics studies how the l ...
*
Bose–Hubbard model
*
t-J model
*
Heisenberg model (quantum)
The quantum Heisenberg model, developed by Werner Heisenberg, is a statistical mechanical model used in the study of critical points and phase transitions of magnetic systems, in which the spins of the magnetic systems are treated quantum mec ...
*
Dynamical mean-field theory Dynamical mean-field theory (DMFT) is a method to determine the electronic structure of strongly correlated materials. In such materials, the approximation of independent electrons, which is used in density functional theory and usual band structure ...
*
Stoner criterion
References
Further reading
*
*
*
*
*
{{DEFAULTSORT:Hubbard Model
Correlated electrons
Condensed matter physics
Quantum chemistry
Lattice models
Quantum lattice models