We introduce a model of interacting bosons exhibiting an infinite collection of fractal symmetries-termed "Pascal's triangle symmetries"-which provides a natural U(1) generalization of a spin-(1/2) system with Sierpinski triangle fractal symmetries introduced in Newman et al., [Phys. Rev. E 60, 5068 (1999).PLEEE81063-651X10.1103/PhysRevE.60.5068]. The Pascal's triangle symmetry gives rise to exact degeneracies, as well as a manifold of low-energy states which are absent in the Sierpinski triangle model. Breaking the U(1) symmetry of this model to Z_{p}, with prime integer p, yields a lattice model with a unique fractal symmetry which is generated by an operator supported on a fractal subsystem with Hausdorff dimension d_{H}=ln(p(p+1)/2)/lnp. The Hausdorff dimension of the fractal can be probed through correlation functions at finite temperature. The phase diagram of these models at zero temperature in the presence of quantum fluctuations, as well as the potential physical construction of the U(1) model, is discussed.
We study the quantum Ising model on the Sierpiński triangle, whose Hausdorff dimension is
log3/log2≈1.585, and demonstrate that it undergoes second-order phase transition with scaling relations satisfied precisely. We also study the quantum
3-state Potts model on the Sierpiński triangle and find first-order phase transition, which is consistent with a prediction from
ε-expansion that the transition becomes first-order for
D>1.3. We then compute critical exponents of the Ising model on higher-dimensional Sierpiński pyramids with various Hausdorff dimension via Monte-Carlo simulations and real-space RG analysis for
D∈[1,3]. We find that only the correlation length exponent
ν interpolates the values of integer-dimensional models. This implies that, contrary to a generally held belief, the universality class of quantum phase transition may not be uniquely determined by symmetry and spatial dimension of the system. This work initiates studies on quantum critical phenomena on graphs and networks which may be of significant importance in the context of quantum networks and communication.
We consider an extended Hubbard model of interacting fermions on a lattice. The fermion kinetic energy corresponds to a tight-binding Hamiltonian with nearest-neighbor (-t) and next-nearest-neighbor (
t∣IH) hopping matrix elements. In addition to the on-site Hubbard interaction (U) we also consider a nearest-neighbor repulsion (V). We obtain the zero-temperature phase diagram of our model within the Hartree-Fock approximation. We consider ground states having charge- and spin-density-wave ordering as well as states with orbital antiferromagnetism or spin nematic order. The latter two states correspond to particle-hole binding with ${\mathrm{d}}_{{\mathrm{x}}^{2}\mathrm{\ensuremath{-}}{\mathrm{y}}^{2}}$ symmetry in the charge and spin channels, respectively. For
t∣IH=0, only the charge-density-wave and spin-density-wave states are energetically stable. For nonzero
t∣IH, we find that orbital antiferromagnetism (or spin nematic) order is stable over a finite portion of the phase diagram at weak coupling. This region of stability is seen to grow with increasing values of
t∣IH.
Abstract We study many-body localization in a hardcore boson model in the presence of random disorder on finite generation fractal lattices with different Hausdorff dimensions and different local lattice structures. In particular, we consider the Vicsek, T-shaped, Sierpinski gasket, and modified Koch-curve fractal lattices. In the single-particle case, these systems display Anderson localization for arbitrary disorder strength if they are large enough. In the many-body case, the systems available to exact diagonalization exhibit a transition between a delocalized and localized regime, visible in the spectral and entanglement properties of these systems. The position of this transition depends on the Hausdorff dimension of the given fractal, as well as on its local structure.
Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of
U/t from 2 to 8, for a range of densities including half-filling and
n=0.3,0.5,0.6,0.75, and
0.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.
The noninteracting electronic structures of tight-binding models on bipartite lattices with unequal numbers of sites in the two sublattices have a number of unique features, including the presence of spatially localized eigenstates and flat bands. When a uniform on-site Hubbard interaction
U is turned on, Lieb proved rigorously that at half-filling ($\ensuremath{\rho}=1$) the ground state has a nonzero spin. In this paper we consider a ``
CuO2 lattice'' (also known as ``Lieb lattice,'' or as a decorated square lattice), in which ``
d orbitals'' occupy the vertices of the squares, while ``
p orbitals'' lie halfway between two
d orbitals; both
d and
p orbitals can accommodate only up to two electrons. We use exact determinant quantum Monte Carlo (DQMC) simulations to quantify the nature of magnetic order through the behavior of correlation functions and sublattice magnetizations in the different orbitals as a function of
U and temperature; we have also calculated the projected density of states, and the compressibility. We study both the homogeneous (H) case,
Ud=Up, originally considered by Lieb, and the inhomogeneous (IH) case, ${U}_{d}\ensuremath{\ne}{U}_{p}$. For the H case at half-filling, we found that the global magnetization rises sharply at weak coupling, and then stabilizes towards the strong-coupling (Heisenberg) value, as a result of the interplay between the ferromagnetism of like sites and the antiferromagnetism between unlike sites; we verified that the system is an insulator for all
U. For the IH system at half-filling, we argue that the case ${U}_{p}\ensuremath{\ne}{U}_{d}$ falls under Lieb's theorem, provided they are positive definite, so we used DQMC to probe the cases
Up=0,Ud=U and
Up=U,Ud=0. We found that the different environments of
d and
p sites lead to a ferromagnetic insulator when
Ud=0; by contrast,
Up=0 leads to to a metal without any magnetic ordering. In addition, we have also established that at density $\ensuremath{\rho}=1/3$, strong antiferromagnetic correlations set in, caused by the presence of one fermion on each
d site.