BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

M. Troyer

M. Troyer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.quant-gascond-mat.stat-mechquant-phcond-mat.dis-nncond-mat.mes-hallcond-mat.supr-conphysics.comp-phcond-mat.mtrl-scicond-mat.otherComputational Engineering, Finance, and Sciencehep-th

Catalog footprint

What is connected

41works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2022arXiv

Propagation of errors and quantitative quantum simulation with quantum advantage

The rapid development in hardware for quantum computing and simulation has led to much interest in problems where these devices can exceed the capabilities of existing classical computers and known methods. Approaching this for problems that go beyond testing the performance of a quantum device is an important step, and quantum simulation of many-body quench dynamics is one of the most promising candidates for early practical quantum advantage. We analyse the requirements for quantitatively reliable quantum simulation beyond the capabilities of existing classical methods for analogue quantum simulators with neutral atoms in optical lattices and trapped ions. Considering the primary sources of error in analogue devices and how they propagate after a quench in studies of the Hubbard or long-range transverse field Ising model, we identify the level of error expected in quantities we extract from experiments. We conclude for models that are directly implementable that regimes of practical quantum advantage are attained in current experiments with analogue simulators. We also identify the hardware requirements to reach the same level of accuracy with future fault-tolerant digital quantum simulation. Verification techniques are already available to test the assumptions we make here, and demonstrating these in experiments will be an important next step.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Towards quantum computing for high-energy excited states in molecular systems: quantum phase estimations of core-level states

This paper explores the utility of the quantum phase estimation (QPE) in calculating high-energy excited states characterized by promotions of electrons occupying inner energy shells. These states have been intensively studied over the last few decades especially in supporting the experimental effort at light sources. Results obtained with the QPE are compared with various high-accuracy many-body techniques developed to describe core-level states. The feasibility of the quantum phase estimator in identifying classes of challenging shake-up states characterized by the presence of higher-order excitation effects is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Autonomous tuning and charge state detection of gate defined quantum dots

Defining quantum dots in semiconductor based heterostructures is an essential step in initializing solid-state qubits. With growing device complexity and increasing number of functional devices required for measurements, a manual approach to finding suitable gate voltages to confine electrons electrostatically is impractical. Here, we implement a two-stage device characterization and dot-tuning process which first determines whether devices are functional and then attempts to tune the functional devices to the single or double quantum dot regime. We show that automating well established manual tuning procedures and replacing the experimenter's decisions by supervised machine learning is sufficient to tune double quantum dots in multiple devices without pre-measured input or manual intervention. The quality of measurement results and charge states are assessed by four binary classifiers trained with experimental data, reflecting real device behaviour. We compare and optimize eight models and different data preprocessing techniques for each of the classifiers to achieve reliable autonomous tuning, an essential step towards scalable quantum systems in quantum dot based qubit architectures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hybrid quantum-classical approach to correlated materials

Recent improvements in control of quantum systems make it seem feasible to finally build a quantum computer within a decade. While it has been shown that such a quantum computer can in principle solve certain small electronic structure problems and idealized model Hamiltonians, the highly relevant problem of directly solving a complex correlated material appears to require a prohibitive amount of resources. Here, we show that by using a hybrid quantum-classical algorithm that incorporates the power of a small quantum computer into a framework of classical embedding algorithms, the electronic structure of complex correlated materials can be efficiently tackled using a quantum computer. In our approach, the quantum computer solves a small effective quantum impurity problem that is self-consistently determined via a feedback loop between the quantum and classical computation. Use of a quantum computer enables much larger and more accurate simulations than with any known classical algorithm, and will allow many open questions in quantum materials to be resolved once a small quantum computer with around one hundred logical qubits becomes available.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Robust Type-II Weyl Semimetal Phase in Transition Metal Diphosphides XP$_2$ (X = Mo, W)

The recently discovered type-II Weyl points appear at the boundary between electron and hole pockets. Type-II Weyl semimetals that host such points are predicted to exhibit a new type of chiral anomaly and possess thermodynamic properties very different from their type-I counterparts. In this Letter, we describe the prediction of a type-II Weyl semimetal phase in the transition metal diphosphides MoP$_2$ and WP$_2$. These materials are characterized by relatively simple band structures with four pairs of type-II Weyl points. Neighboring Weyl points have the same chirality, which makes the predicted topological phase robust with respect to small perturbations of the crystalline lattice. In addition, this peculiar arrangement of the Weyl points results in long topological Fermi arcs, thus making them readily accessible in angle-resolved photoemission spectroscopy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Training A Quantum Optimizer

We study a variant of the quantum approximate optimization algorithm [ E. Farhi, J. Goldstone, and S. Gutmann, arXiv:1411.4028] with slightly different parametrization and different objective: rather than looking for a state which approximately solves an optimization problem, our goal is to find a quantum algorithm that, given an instance of MAX-2-SAT, will produce a state with high overlap with the optimal state. Using a machine learning approach, we chose a "training set" of instances and optimized the parameters to produce large overlap for the training set. We then tested these optimized parameters on a larger instance set. As a training set, we used a subset of the hard instances studied by E. Crosson, E. Farhi, C. Yen-Yu Lin, H.-H. Lin, and P. Shor (CFLLS) [arXiv:1401.7320]. When tested on the full set, the parameters that we find produce significantly larger overlap than the optimized annealing times of CFLLS. Testing on other random instances from $20$ to $28$ bits continues to show improvement over annealing, with the improvement being most notable on the hardest instances. Further tests on instances of MAX-3-SAT also showed improvement on the hardest instances. This algorithm may be a possible application for near-term quantum computers with limited coherence times.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Bespoke Single-Band Hubbard Model Material

The Hubbard model, which augments independent-electron band theory with a single parameter to describe electron-electron correlations, is widely regarded to be the `standard model' of condensed matter physics. The model has been remarkably successful at addressing a range of correlation effects in solids, but beyond one dimension its solution is intractable. Much current research aims, therefore, at finding appropriate approximations to the Hubbard model phase diagram. Here we take the new approach of using ab initio electronic structure methods to design a material whose Hamiltonian is that of the single-band Hubbard model. Solution of the Hubbard model will then be available through measurement of the material's properties. After identifying an appropriate crystal class and several appropriate chemistries, we use density functional theory and dynamical mean-field theory to screen for the desired electronic band structure and metal-insulator transition. We then explore the most promising candidates for structural stability and suitability for doping and propose specific materials for subsequent synthesis. Finally, we identify a regime -- that should manifest in our bespoke material -- in which the single-band Hubbard model on a triangular lattice exhibits exotic d-wave superconductivity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimized simulated annealing for Ising spin glasses

We present several efficient implementations of the simulated annealing algorithm for Ising spin glasses on sparse graphs. In particular, we provide a generic code for any choice of couplings, an optimized code for bipartite graphs, and highly optimized implementations using multi-spin coding for graphs with small maximum degree and discrete couplings with a finite range. The latter codes achieve up to 50 spin flips per nanosecond on modern Intel CPUs. We also compare the performance of the codes to that of the special purpose D-Wave devices built for solving such Ising spin glass problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Towards Practical Quantum Variational Algorithms

The preparation of quantum states using short quantum circuits is one of the most promising near-term applications of small quantum computers, especially if the circuit is short enough and the fidelity of gates high enough that it can be executed without quantum error correction. Such quantum state preparation can be used in variational approaches, optimizing parameters in the circuit to minimize the energy of the constructed quantum state for a given problem Hamiltonian. For this purpose we propose a simple-to-implement class of quantum states motivated by adiabatic state preparation. We test its accuracy and determine the required circuit depth for a Hubbard model on ladders with up to 12 sites (24 spin-orbitals), and for small molecules. We find that this ansatz converges faster than previously proposed schemes based on unitary coupled clusters. While the required number of measurements is astronomically large for quantum chemistry applications to molecules, applying the variational approach to the Hubbard model (and related models) is found to be far less demanding and potentially practical on small quantum computers. We also discuss another application of quantum state preparation using short quantum circuits, to prepare trial ground states of models faster than using adiabatic state preparation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two-dimensional epitaxial superconductor-semiconductor heterostructures: A platform for topological superconducting networks

Progress in the emergent field of topological superconductivity relies on synthesis of new material combinations, combining superconductivity, low density, and spin-orbit coupling (SOC). For example, theory [1-4] indicates that the interface between a one-dimensional (1D) semiconductor (Sm) with strong SOC and a superconductor (S) hosts Majorana modes with nontrivial topological properties [5-8]. Recently, epitaxial growth of Al on InAs nanowires was shown to yield a high quality S-Sm system with uniformly transparent interfaces [9] and a hard induced gap, indicted by strongly suppressed sub gap tunneling conductance [10]. Here we report the realization of a two-dimensional (2D) InAs/InGaAs heterostructure with epitaxial Al, yielding a planar S-Sm system with structural and transport characteristics as good as the epitaxial wires. The realization of 2D epitaxial S-Sm systems represent a significant advance over wires, allowing extended networks via top-down processing. Among numerous potential applications, this new material system can serve as a platform for complex networks of topological superconductors with gate-controlled Majorana zero modes [1-4]. We demonstrate gateable Josephson junctions and a highly transparent 2D S-Sm interface based on the product of excess current and normal state resistance.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A model project for reproducible papers: critical temperature for the Ising model on a square lattice

In this paper we present a simple, yet typical simulation in statistical physics, consisting of large scale Monte Carlo simulations followed by an involved statistical analysis of the results. The purpose is to provide an example publication to explore tools for writing reproducible papers. The simulation estimates the critical temperature where the Ising model on the square lattice becomes magnetic to be Tc /J = 2.26934(6) using a finite size scaling analysis of the crossing points of Binder cumulants. We provide a virtual machine which can be used to reproduce all figures and results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ferromagnetism of a Repulsive Atomic Fermi Gas in an Optical Lattice: a Quantum Monte Carlo Study

Using continuous-space quantum Monte Carlo methods we investigate the zero-temperature ferromagnetic behavior of a two-component repulsive Fermi gas under the influence of periodic potentials that describe the effect of a simple-cubic optical lattice. Simulations are performed with balanced and with imbalanced components, including the case of a single impurity immersed in a polarized Fermi sea (repulsive polaron). For an intermediate density below half filling, we locate the transitions between the paramagnetic, and the partially and the fully ferromagnetic phases. As the intensity of the optical lattice increases, the ferromagnetic instability takes place at weaker interactions, indicating a possible route to observe ferromagnetism in experiments performed with ultracold atoms. We compare our findings with previous predictions based on the standard computational method used in material science, namely density functional theory, and with results based on tight-binding models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Improving Quantum Algorithms for Quantum Chemistry

We present several improvements to the standard Trotter-Suzuki based algorithms used in the simulation of quantum chemistry on a quantum computer. First, we modify how Jordan-Wigner transformations are implemented to reduce their cost from linear or logarithmic in the number of orbitals to a constant. Our modification does not require additional ancilla qubits. Then, we demonstrate how many operations can be parallelized, leading to a further linear decrease in the parallel depth of the circuit, at the cost of a small constant factor increase in number of qubits required. Thirdly, we modify the term order in the Trotter-Suzuki decomposition, significantly reducing the error at given Trotter-Suzuki timestep. A final improvement modifies the Hamiltonian to reduce errors introduced by the non-zero Trotter-Suzuki timestep. All of these techniques are validated using numerical simulation and detailed gate counts are given for realistic molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-local quantum fluctuations and fermionic superfluidity in the imbalanced attractive Hubbard model

We study fermionic superfluidity in strongly anisotropic optical lattices with attractive interactions utilizing the cluster DMFT method, and focusing in particular on the role of non-local quantum fluctuations. We show that non-local quantum fluctuations impact the BCS superfluid transition dramatically. Moreover, we show that exotic superfluid states with delicate order parameter structure, such as the Fulde-Ferrell-Larkin-Ovchinnikov phase driven by spin population imbalance, can emerge even in the presence of such strong fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Critical Temperature of Interacting Bose Gases in Periodic Potentials

The superfluid transition of a repulsive Bose gas in the presence of a sinusoidal potential which represents a simple-cubic optical lattice is investigate using quantum Monte Carlo simulations. At the average filling of one particle per well the critical temperature has a nonmonotonic dependence on the interaction strength, with an initial sharp increase and a rapid suppression at strong interactions in the vicinity of the Mott transition. In an optical lattice the positive shift of the transition is strongly enhanced compared to the homogenous gas. By varying the lattice filling we find a crossover from a regime where the optical lattice has the dominant effect to a regime where interactions dominate and the presence of the lattice potential becomes almost irrelevant.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bosonic superfluid-insulator transition in continuous space

We investigate the zero-temperature phase diagram of interacting Bose gases in the presence of a simple cubic optical lattice, going beyond the regime where the mapping to the single-band Bose-Hubbard model is reliable. Our computational approach is a new hybrid quantum Monte Carlo method which combines algorithms used to simulate homogeneous quantum fluids in continuous space with those used for discrete lattice models of strongly correlated systems. We determine the critical interaction strength and optical lattice intensity where the superfluid-to-insulator transition takes place, considering also the regime of shallow optical lattices and strong inter-atomic interactions. The implications of our findings for the super-solid state of matter are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dipolar Dynamics for Interacting Ultracold Fermions in a Trapped Optical Lattice

Using the time-dependent density matrix renormalization group method, we calculate transport properties of an interacting Fermi gas in an optical lattice with a confining trap after a sudden displacement of the trap center. In the regime of attractive interactions, the dipolar motion after the displacement can be classified into underdamped oscillations and overdamped relaxations, depending on the interaction strength. These numerical calculations are consistent with experimental results. In the regime of repulsive interactions, we predict a revival of the oscillations of the center of mass when the interaction strength is increased. This unique feature can be considered as a dynamical signature for the emergence of a Mott plateau for an interacting trapped Fermi gas in an optical lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Magnetic susceptibility of Ce: an LDA+DMFT study

The magnetic properties of Ce in the $α$ and $γ$ phase are calculated within the LDA+DMFT approach. The magnetic susceptibility in these two phases shows a similar behavior over a wide temperature range: a Curie-Weiss law at high temperatures, indicating the presence of local moments, followed by a maximum in a crossover regime, and a saturation characteristic of a state with screened local moments at low temperature. The difference in experimentally observable magnetic properties is caused by the shift of the susceptibility to higher temperatures in the $α-$phase compared to the $γ-$phase, due to the larger hybridization between localized $4f-$states and conductive $spd-$electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multigrid Algorithms for Tensor Network States

The widely used density matrix renormalization group (DRMG) method often fails to converge in systems with multiple length scales, such as lattice discretizations of continuum models and dilute or weakly doped lattice models. The local optimization employed by DMRG to optimize the wave function is ineffective in updating large-scale features. Here we present a multigrid algorithm that solves these convergence problems by optimizing the wave function at different spatial resolutions. We demonstrate its effectiveness by simulating bosons in continuous space, and study non-adiabaticity when ramping up the amplitude of an optical lattice. The algorithm can be generalized to tensor network methods, and be combined with the contractor renormalization group (CORE) method to study dilute and weakly doped lattice models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Strong-disorder renormalization for interacting non-Abelian anyon systems in two dimensions

We consider the effect of quenched spatial disorder on systems of interacting, pinned non-Abelian anyons as might arise in disordered Hall samples at filling fractions ν=5/2 or ν=12/5. In one spatial dimension, such disordered anyon models have previously been shown to exhibit a hierarchy of infinite randomness phases. Here, we address systems in two spatial dimensions and report on the behavior of Ising and Fibonacci anyons under the numerical strong-disorder renormalization group (SDRG). In order to manage the topology-dependent interactions generated during the flow, we introduce a planar approximation to the SDRG treatment. We characterize this planar approximation by studying the flow of disordered hard-core bosons and the transverse field Ising model, where it successfully reproduces the known infinite randomness critical point with exponent ψ~ 0.43. Our main conclusion for disordered anyon models in two spatial dimensions is that systems of Ising anyons as well as systems of Fibonacci anyons do not realize infinite randomness phases, but flow back to weaker disorder under the numerical SDRG treatment.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Diagrammatic Quantum Monte Carlo solution of the two-dimensional Cooperon-Fermion model

We investigate the two-dimensional cooperon-fermion model in the correlated regime with a new continuous-time diagrammatic determinant quantum Monte Carlo (DDQMC) algorithm. We estimate the transition temperature $T_{c}$, examine the effectively reduced band gap and cooperon mass, and find that delocalization of the cooperons enhances the diamagnetism. When applied to diamagnetism of the pseudogap phase in high-$T_{c}$ cuprates, we obtain results in a qualitative agreement with recent torque magnetization measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamics at and near conformal quantum critical points

We explore the dynamical behavior at and near a special class of two-dimensional quantum critical points. Each is a conformal quantum critical point (CQCP), where in the scaling limit the equal-time correlators are those of a two-dimensional conformal field theory. The critical theories include the square-lattice quantum dimer model, the quantum Lifshitz theory, and a deformed toric code model. We show that under generic perturbation the latter flows toward the ordinary Lorentz-invariant (2+1) dimensional Ising critical point, illustrating that CQCPs are generically unstable. We exploit a correspondence between the classical and quantum dynamical behavior in such systems to perform an extensive numerical study of two lines of CQCPs in a quantum eight-vertex model, or equivalently, two coupled deformed toric codes. We find that the dynamical critical exponent z remains 2 along the U(1)-symmetric quantum Lifshitz line, while it continuously varies along the line with only Z_2 symmetry. This illustrates how two CQCPs can have very different dynamical properties, despite identical equal-time ground-state correlators. Our results equally apply to the dynamics of the corresponding purely classical models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Identifying Quantum Topological Phases Through Statistical Correlation

We theoretically examine the use of a statistical distance measure, the indistinguishability, as a generic tool for the identification of topological order. We apply this measure to the toric code and two fractional quantum Hall models. We find that topologically ordered states can be identified with the indistinguishability for both models. Calculations with the indistinguishability also underscore a key distinction between symmetries that underlie topological order in the toric code and quantum Hall models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Implementing global Abelian symmetries in projected entangled-pair state algorithms

Due to the unfavorable scaling of tensor network methods with the refinement parameter M, new approaches are necessary to improve the efficiency of numerical simulations based on such states in particular for gapless, strongly entangled systems. In one-dimensional DMRG, the use of Abelian symmetries has lead to large computational gain. In higher-dimensional tensor networks, this is associated with significant technical efforts and additional approximations. We explain a formalism to implement such symmetries in two-dimensional tensor network states and present benchmark results that confirm the validity of these approximations in the context of projected entangled-pair state algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stripes in the two-dimensional t-J model with infinite projected entangled-pair states

We simulate the $t$ $J$ model in two dimensions by means of infinite projected entangled-pair states (iPEPS) generalized to arbitrary unit cells, finding results similar to those previously obtained by the density-matrix renormalization group (DMRG) for wide ladders. In particular, we show that states exhibiting stripes, {\it i.e}. a unidirectional modulation of hole-density and antiferromagnetic order with a $π$-phase shift between adjacent stripes, have a lower variational energy than uniform phases predicted by variational and fixed-node Monte Carlo simulations. For a fixed unit-cell size the energy per hole is minimized for a hole density $ρ_l\sim 0.5$ per unit length of a stripe. The superconducting order parameter is maximal around $ρ_l\sim 0.75-0.8$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The ALPS project release 2.0: Open source software for strongly correlated systems

We present release 2.0 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. The code development is centered on common XML and HDF5 data formats, libraries to simplify and speed up code development, common evaluation and plotting tools, and simulation programs. The programs enable non-experts to start carrying out serial or parallel numerical simulations by providing basic implementations of the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended ensemble simulations, exact and full diagonalization (ED), the density matrix renormalization group (DMRG) both in a static version and a dynamic time-evolving block decimation (TEBD) code, and quantum Monte Carlo solvers for dynamical mean field theory (DMFT). The ALPS libraries provide a powerful framework for programers to develop their own applications, which, for instance, greatly simplify the steps of porting a serial code onto a parallel, distributed memory machine. Major changes in release 2.0 include the use of HDF5 for binary data, evaluation tools in Python, support for the Windows operating system, the use of CMake as build system and binary installation packages for Mac OS X and Windows, and integration with the VisTrails workflow provenance tool. The software is available from our web server at http://alps.comp-phys.org/.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Trapped Ultracold Bosons in Periodically Modulated Lattices

Motivated by the recent rapid development of the field of quantum gases in optical lattices, we present a comprehensive study of the spectrum of ultracold atoms in a one-dimensional optical lattice subjected to a periodic lattice modulation. Using the time-dependent density-matrix renormalization group method, we study the dynamical response due to lattice modulations in different quantum phases of the system with varying density. For the Mott insulating state, we identify several excitation processes, which provide important information about the density profile of the gases. For the superfluid, the dynamical response can be well described in a local density approximation. This simplification can be valuable in understanding the strong-correlated superfluid in a slow-varying harmonic potential. All these spectroscopic features of an inhomogeneous system can be used as a test for the validity of the Bose-Hubbard model in a parabolic trapping potential.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Diagrammatic Monte Carlo for Correlated Fermions

We show that Monte Carlo sampling of the Feynman diagrammatic series (DiagMC) can be used for tackling hard fermionic quantum many-body problems in the thermodynamic limit by presenting accurate results for the repulsive Hubbard model in the correlated Fermi liquid regime. Sampling Feynman's diagrammatic series for the single-particle self-energy we can study moderate values of the on-site repulsion ($U/t \sim 4$) and temperatures down to $T/t=1/40$. We compare our results with high temperature series expansion and with single-site and cluster dynamical mean-field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Itinerant ferromagnetism of a repulsive atomic Fermi gas: a quantum Monte Carlo study

We investigate the phase diagram of a two-component repulsive Fermi gas at T=0 by means of quantum Monte Carlo simulations. For a given value of the positive s-wave scattering length, both purely repulsive and purely attractive model potentials are considered in order to analyze the limits of the universal regime where the details of interatomic forces can be neglected. The equation of state of both balanced and unbalanced systems is calculated as a function of the interaction strength and the critical density for the onset of ferromagnetism is determined. The energy per particle of the strongly polarized gas is calculated and parametrized in terms of the physical properties of repulsive polarons, which are relevant for the stability of the fully magnetized ferromagnetic state. Finally, we analyze the phase diagram in the polarization/interaction plane under the assumption that only phases with homogeneous magnetization can be produced.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermodynamics of the 3D Hubbard model on approach to the Neel transition

We study the thermodynamic properties of the 3D Hubbard model for temperatures down to the Neel temperature using cluster dynamical mean-field theory. In particular we calculate the energy, entropy, density, double occupancy and nearest-neighbor spin correlations as a function of chemical potential, temperature and repulsion strength. To make contact with cold-gas experiments, we also compute properties of the system subject to an external trap in the local density approximation. We find that an entropy per particle $S/N \approx 0.65(6)$ at $U/t=8$ is sufficient to achieve a Neel state in the center of the trap, substantially higher than the entropy required in a homogeneous system. Precursors to antiferromagnetism can clearly be observed in nearest-neighbor spin correlators.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Absence of a Direct Superfluid to Mott Insulator Transition in Disordered Bose Systems

We prove the absence of a direct quantum phase transition between a superfluid and a Mott insulator in a bosonic system with generic, bounded disorder. We also prove compressibility of the system on the superfluid--insulator critical line and in its neighborhood. These conclusions follow from a general {\it theorem of inclusions} which states that for any transition in a disordered system one can always find rare regions of the competing phase on either side of the transition line. Quantum Monte Carlo simulations for the disordered Bose-Hubbard model show an even stronger result, important for the nature of the Mott insulator to Bose glass phase transition: The critical disorder bound, $Δ_c$, corresponding to the onset of disorder-induced superfluidity, satisfies the relation $Δ_c > E_{\rm g/2}$, with $E_{\rm g/2}$ the half-width of the Mott gap in the pure system.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Beliaev technique for a weakly interacting Bose gas

Aiming for simplicity of explicit equations and at the same time controllable accuracy of the theory we present results for all thermodynamic quantities and correlation functions for the weakly interacting Bose gas at short-to-intermediate distances obtained within an improved version of Beliaev's diagrammatic technique. With a small symmetry breaking term Beliaev's diagrammatic technique becomes regular in the infrared limit. Up to higher-order terms (for which we present order-of-magnitude estimates), the partition function and entropy of the system formally correspond to those of a non-interacting bosonic (pseudo-)Hamiltonian with a temperature dependent Bogoliubov-type dispersion relation. Away from the fluctuation region, this approach provides the most accurate--in fact, the best possible within the Bogoliubov-type pseudo-Hamiltonian framework--description of the system with controlled accuracy. It produces accurate answers for the off-diagonal correlation functions up to distances where the behaviour of correlators is controlled by generic hydrodynamic relations, and thus can be accurately extrapolated to arbitrarily large scales. In the fluctuation region, the non-perturbative contributions are given by universal (for all weakly interacting U(1) systems) constants and scaling functions, which can be obtained separately--by simulating classical U(1) models--and then used to extend the description of the weakly interacting Bose gas to the fluctuation region. The theory works in all spatial dimensions and we explicitly check its validity against first-principle Monte Carlo simulations for various thermodynamic properties and the single-particle density matrix.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Distinguishing Phases with Ansatz Wavefunctions

We propose an indistinguishability measure for assessment of ansatz wavefunctions with numerically determined wavefunctions. The measure efficiently compares all correlation functions of two states and can therefore be used to distinguish phases by defining correlator classes for ansatz wavefunctions. It also allows identification of quantum critical points. We demonstrate the approach for the transverse Ising model, using the matrix product state formalism with the time evolving block decimation algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Phase diagram of the disordered Bose-Hubbard model

We establish the phase diagram of the disordered three-dimensional Bose-Hubbard model at unity filling, which has been controversial for many years. The theorem of inclusions, proven in Ref. [1], states that the Bose glass phase always intervenes between the Mott insulating and superfluid phases. Here, we note that assumptions on which the theorem is based exclude phase transitions between gapped (Mott insulator) and gapless phases (Bose glass). The apparent paradox is resolved through a unique mechanism: such transitions have to be of the Griffiths type when the vanishing of the gap at the critical point is due to a zero concentration of rare regions where extreme fluctuations of disorder mimic a {\it regular} gapless system. An exactly solvable random transverse field Ising model in one dimension is used to illustrate the point. A highly non-trivial overall shape of the phase diagram is revealed with the worm algorithm. The phase diagram features a long superfluid finger at strong disorder and on-site interaction. Moreover, bosonic superfluidity is extremely robust against disorder in a broad range of interaction parameters; it persists in random potentials nearly 50 (!) times larger than the particle half-bandwidth. Finally, we comment on the feasibility of obtaining this phase diagram in cold-atom experiments, which work with trapped systems at finite temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Supersolid phase with cold polar molecules on a triangular lattice

We study a system of heteronuclear molecules on a triangular lattice and analyze the potential of this system for the experimental realization of a supersolid phase. The ground state phase diagram contains superfluid, solid and supersolid phases. At finite temperatures and strong interactions there is an additional emulsion region, in contrast to similar models with short-range interactions. We derive the maximal critical temperature $T_c$ and the corresponding entropy $S/N = 0.04(1)$ for supersolidity and find feasible experimental conditions for its realization.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Suppression of the critical temperature for superfluidity near the Mott transition: validating a quantum simulator

Ultracold atomic gases in optical lattices have proven to be a controllable, tunable and clean implementation of strongly interacting quantum many-body systems. An essential prospect for such quantum simulators is their ability to map out the phase diagram of fundamental many-body model Hamiltonians. However, the results need to be validated first for representative benchmark problems via state-of-the-art numerical methods of quantum many-body theory. Here we present the first ab-initio comparison between experiments and quantum Monte Carlo simulations for strongly interacting Bose gases on a lattice for large systems (up to N = 3e5 particles). The comparison has enabled us to perform thermometry for the interacting quantum gas and to experimentally determine the finite temperature phase diagram for bosonic superfluids in an optical lattice. Our results reveal a downshift of the critical temperature as the transition to the Mott insulator is approached.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Local stress and superfluid properties of solid Helium-4

More than half a century ago Penrose asked: are the superfluid and solid state of matter mutually exclusive or do there exist "supersolid" materials where the atoms form a regular lattice and simultaneously flow without friction? Recent experiments provide evidence that supersolid behavior indeed exists in Helium-4 -- the most quantum material known in Nature. In this paper we show that large local strain in the vicinity of crystalline defects is the origin of supersolidity in Helium-4. Although ideal crystals of Helium-4 are not supersolid, the gap for vacancy creation closes when applying a moderate stress. While a homogeneous system simply becomes unstable at this point, the stressed core of crystalline defects (dislocations and grain boundaries) undergoes a radical transformation and can become superfluid.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Quantum Monte Carlo study of a one-dimensional phase-fluctuating condensate

We study numerically the low temperature behavior of a one-dimensional Bose gas trapped in an optical lattice. For a sufficient number of particles and weak repulsive interactions, we find a clear regime of temperatures where density fluctuations are negligible but phase fluctuations are considerable, i.e., a quasicondensate. In the weakly interacting limit, our results are in very good agreement with those obtained using a mean-field approximation. In coupling regimes beyond the validity of mean-field approaches, a phase-fluctuating condensate also appears, but the phase-correlation properties are qualitatively different. It is shown that quantum depletion plays an important role.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Superfluidity of Grain Boundaries in Solid Helium-4

By large-scale quantum Monte Carlo simulations we show that grain boundaries in Helium-4 crystals are generically superfluid at low temperature, with a transition temperature of the order of ~0.5K at the melting pressure; insulating grain boundaries are found only for special orientations of the grains. We also find that close vicinity to the melting line is not a necessary condition for superfluid grain boundaries, and a grain boundary in direct contact with the superfluid liquid at the melting curve is found to be mechanically stable and the grain boundary superfluidity observed by Sasaki et al. [Science 313, 1098 (2006)] is not just a crack filled with superfluid.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Bosons in optical lattices - from the Mott transition to the Tonks-Girardeau gas

We present results from quantum Monte Carlo simulations of trapped bosons in optical lattices, focusing on the crossover from a gas of softcore bosons to a Tonks-Girardeau gas in a one-dimensional optical lattice. We find that depending on the quantity being measured, the behavior found in the Tonks-Girardeau regime is observed already at relatively small values of the interaction strength. A finite critical value for entering the Tonks-Girardeau regime does not exist. Furthermore, we discuss the computational efficiency of two quantum Monte Carlo methods to simulate large scale trapped bosonic systems: directed loops in stochastic series expansions and the worm algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Thermodynamics of Spin S = 1/2 Antiferromagnetic Uniform and Alternating-Exchange Heisenberg Chains

The magnetic susceptibility chi and specific heat C versus temperature T of the spin-1/2 antiferromagnetic alternating-exchange (J1 and J2) Heisenberg chain are studied for the entire range 0 \leq alpha \leq 1 of the alternation parameter alpha = J2/J1. For the uniform chain (alpha = 1), detailed comparisons of the high-accuracy chi(T) and C(T) Bethe ansatz data of Kluemper and Johnston are made with the asymptotically exact low-T field theory predictions of Lukyanov. QMC simulations and TMRG calculations of chi(alpha,T) are presented. From the low-T TMRG data, the spin gap Delta(alpha)/J1 is extracted for 0.8 \leq alpha \leq 0.995. High accuracy fits to all of the above numerical data are obtained. We examine in detail the theoretical predictions of Bulaevskii for chi(alpha,T) and compare them with our results. Our experimental chi(T) and C(T) data for NaV2O5 single crystals are modeled in detail. The chi(T) data above the spin dimerization temperature Tc = 34 K are not in agreement with the prediction for the uniform Heisenberg chain, but can be explained if there is a moderate ferromagnetic interchain coupling and/or if J changes with T. By fitting the chi(T) data, we obtain Delta(T = 0) = 103(2) K, alternation parameter delta(0) = (1 - alpha)/(1 + alpha) = 0.034(6) and average exchange constant J(0) = 640(80) K. The delta(T) and Delta(T) are derived from the data. A spin pseudogap with a large magnitude \approx 0.4 Delta(0) is consistently found just above Tc, which decreases with increasing T. Analysis of our C(T) data indicates that at Tc, at least 77% of the entropy change due to the transition at Tc and associated order parameter fluctuations arise from the lattice and/or charge degrees of freedom and less than 23% from the spin degrees of freedom.

0 citations0 reviewsNo code linkedOpen access