We study interacting Majorana fermions in two dimensions as a low-energy effective model of a vortex lattice in two-dimensional time-reversal-invariant topological superconductors. For that purpose, we implement ab-initio quantum Monte Carlo simulation to the Majorana fermion system in which the path-integral measure is given by a semi-positive Pfaffian. We discuss spontaneous breaking of time-reversal symmetry at finite temperature.
We report large scale determinant Quantum Monte Carlo calculations of the effective bandwidth, momentum distribution, and magnetic correlations of the square lattice fermion Hubbard Hamiltonian at half-filling. The sharp Fermi surface of the non-interacting limit is significantly broadened by the electronic correlations, but retains signatures of the approach to the edges of the first Brillouin zone as the density increases. Finite size scaling of simulations on large lattices allows us to extract the interaction dependence of the antiferromagnetic order parameter, exhibiting its evolution from weak coupling to the strong coupling Heisenberg limit. Our lattices provide improved resolution of the Greens function in momentum space, allowing a more quantitative comparison with time-of-flight optical lattice experiments.
We present different methods to increase the performance of Hybrid Monte Carlo simulations of the Hubbard model in two-dimensions. Our simulations concentrate on a hexagonal lattice, though can be easily generalized to other lattices. It is found that best results can be achieved using a flexible GMRES solver for matrix
The interplay between lattice gauge theories and fermionic matter accounts for fundamental physical phenomena ranging from the deconfinement of quarks in particle physics to quantum spin liquid with fractionalized anyons and emergent gauge structures in condensed matter physics. However, except for certain limits (for instance large number of flavors of matter fields), analytical methods can provide few concrete results. Here we show that the problem of compact $U(1)$ lattice gauge theory coupled to fermionic matter in $(2+1)$D is possible to access via sign-problem-free quantum Monte Carlo simulations. One can hence map out the phase diagram as a function of fermion flavors and the strength of gauge fluctuations. By increasing the coupling constant of the gauge field, gauge confinement in the form of various spontaneous symmetry breaking phases such as valence bond solid (VBS) and Neel antiferromagnet emerge. Deconfined phases with algebraic spin and VBS correlation functions are also observed. Such deconfined phases are an incarnation of exotic states of matter, $i.e.$ the algebraic spin liquid, which is generally viewed as the parent state of various quantum phases. The phase transitions between deconfined and confined phases, as well as that between the different confined phases provide various manifestations of deconfined quantum criticality. In particular, for four flavors, $N_f = 4$, our data suggests a continuous quantum phase transition between the VBS and N{e}el order. We also provide preliminary theoretical analysis for these quantum phase transitions.
We investigate the interaction potential of superconducting vortices at the full quantum level. We formulate the interaction potential in a constrained path integral and calculate it by the quantum Monte Carlo simulation. The vortex-vortex potential is attractive (type-I), repulsive (type-II), and flat (critical) depending on a coupling constant. The vortex-antivortex potential also depends on the coupling constant at long range but is always attractive at short range.
We have performed numerical studies of the Hubbard-Holstein model in two dimensions using determinant quantum Monte Carlo (DQMC). Here we present details of the method, emphasizing the treatment of the lattice degrees of freedom, and then study the filling and behavior of the fermion sign as a function of model parameters. We find a region of parameter space with large Holstein coupling where the fermion sign recovers despite large values of the Hubbard interaction. This indicates that studies of correlated polarons at finite carrier concentrations are likely accessible to DQMC simulations. We then restrict ourselves to the half-filled model and examine the evolution of the antiferromagnetic structure factor, other metrics for antiferromagnetic and charge-density-wave order, and energetics of the electronic and lattice degrees of freedom as a function of electron-phonon coupling. From this we find further evidence for a competition between charge-density-wave and antiferromagnetic order at half-filling.