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

Bryan K. Clark

Bryan K. Clark 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-elquant-phphysics.comp-phcond-mat.dis-nncond-mat.quant-gascond-mat.othercond-mat.stat-mechhep-latcond-mat.mtrl-scicond-mat.mes-hallcond-mat.supr-conDistributed, Parallel, and Cluster Computinghep-phhep-thMachine Learningmath.NAphysics.atom-phphysics.chem-ph

Catalog footprint

What is connected

27works
18topics
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

27 published item(s)

preprint2022arXiv

Gauge equivariant neural networks for quantum lattice gauge theories

Gauge symmetries play a key role in physics appearing in areas such as quantum field theories of the fundamental particles and emergent degrees of freedom in quantum materials. Motivated by the desire to efficiently simulate many-body quantum systems with exact local gauge invariance, gauge equivariant neural-network quantum states are introduced, which exactly satisfy the local Hilbert space constraints necessary for the description of quantum lattice gauge theory with Zd gauge group on different geometries. Focusing on the special case of Z2 gauge group on a periodically identified square lattice, the equivariant architecture is analytically shown to contain the loop-gas solution as a special case. Gauge equivariant neural-network quantum states are used in combination with variational quantum Monte Carlo to obtain compact descriptions of the ground state wavefunction for the Z2 theory away from the exactly solvable limit, and to demonstrate the confining/deconfining phase transition of the Wilson loop order parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Simulating Quantum Mechanics with a $θ$-term and an 't Hooft Anomaly on a Synthetic Dimension

A topological $θ$-term in gauge theories, including quantum chromodynamics in 3+1 dimensions, gives rise to a sign problem that makes classical Monte Carlo simulations impractical. Quantum simulations are not subject to such sign problems and are a promising approach to studying these theories in the future. In the near term, it is interesting to study simpler models that retain some of the physical phenomena of interest and their implementation on quantum hardware. For example, dimensionally-reducing gauge theories on small spatial tori produces quantum mechanical models which, despite being relatively simple to solve, retain interesting vacuum and symmetry structures from the parent gauge theories. Here we consider quantum mechanical particle-on-a-circle models, related by dimensional reduction to the 1+1d Schwinger model, that possess a $θ$-term and realize an 't Hooft anomaly or global inconsistency at $θ= π$. These models also exhibit the related phenomena of spontaneous symmetry breaking and instanton-anti-instanton interference in real time. We propose an experimental scheme for the real-time simulation of a particle on a circle with a $θ$-term and a $\mathbb{Z}_n$ potential using a synthetic dimension encoded in a Rydberg atom. Simulating the Rydberg atom with realistic experimental parameters, we demonstrate that the essential physics can be well-captured by the experiment, with expected behavior in the tunneling rate as a function of $θ$. Similar phenomena and observables can also arise in more complex quantum mechanical models connected to higher-dimensional nonabelian gauge theories by dimensional reduction.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Unitary Block Optimization for Variational Quantum Algorithms

Variational quantum algorithms are a promising hybrid framework for solving chemistry and physics problems with broad applicability to optimization as well. They are particularly well suited for noisy intermediate scale quantum (NISQ) computers. In this paper, we describe the unitary block optimization scheme (UBOS) and apply it to two variational quantum algorithms: the variational quantum eigensolver (VQE) and variational time evolution. The goal of VQE is to optimize a classically intractable parameterized quantum wave function to target a physical state of a Hamiltonian or solve an optimization problem. UBOS is an alternative to other VQE optimization schemes with a number of advantages including fast convergence, less sensitivity to barren plateaus, the ability to tunnel through some local minima and no hyperparameters to tune. We additionally describe how UBOS applies to real and imaginary time-evolution (TUBOS).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

(1+1)-d U(1) Quantum link models from effective Hamiltonians of dipolar molecules

We study the promising idea of using dipolar molecular systems as analog quantum simulators for quantum link models, which are discrete versions of lattice gauge theories. In a quantum link model the link variables have a finite number of degrees of freedom and discrete values. We construct the effective Hamiltonian of a system of dipolar molecules with electric dipole-dipole interactions, where we use the tunable parameters of the system to match it to the target Hamiltonian describing a U(1) quantum link model in 1+1 dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Characterizing the many-body localization transition through correlations

Closed, interacting, quantum systems have the potential to transition to a many-body localized (MBL) phase under the presence of sufficiently strong disorder, hence breaking ergodicity and failing to thermalize. In this work we study the distribution of correlations throughout the ergodic-MBL phase diagram. We find the typical correlations in the MBL phase decay as a stretched exponential with range $r$ eventually crossing over to an exponential decay deep in the MBL phase. At the transition, the stretched exponential goes as $e^{-A\sqrt{r}}$, a decay that is reminiscent of the random singlet phase. While the standard deviation of the $\log(QMI)$ has a range dependence, the $\log(QMI)$ converges to a range-invariant distribution on all other moments (i.e., the skewness and higher) at the transition. The universal nature of these distributions provides distinct phenomenology of the transition different from both the ergodic and MBL phenomenologies. In addition to the typical correlations, we study the extreme correlations in the system, finding that the probability of strong long-range correlations is maximal at the transition, suggesting the proliferation of resonances there. Finally, we analyze the probability that a single bit of information is shared across two halves of a system, finding that this probability is non-zero deep in the MBL phase but vanishes at moderate disorder well above the transition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deep Learning Enabled Strain Mapping of Single-Atom Defects in 2D Transition Metal Dichalcogenides with Sub-picometer Precision

2D materials offer an ideal platform to study the strain fields induced by individual atomic defects, yet challenges associated with radiation damage have so-far limited electron microscopy methods to probe these atomic-scale strain fields. Here, we demonstrate an approach to probe single-atom defects with sub-picometer precision in a monolayer 2D transition metal dichalcogenide, WSe$_{2-2x}$Te$_{2x}$. We utilize deep learning to mine large datasets of aberration-corrected scanning transmission electron microscopy images to locate and classify point defects. By combining hundreds of images of nominally identical defects, we generate high signal-to-noise class-averages which allow us to measure 2D atomic coordinates with up to 0.3 pm precision. Our methods reveal that Se vacancies introduce complex, oscillating strain fields in the WSe$_{2-2x}$Te$_{2x}$ lattice which cannot be explained by continuum elastic theory. These results indicate the potential impact of computer vision for the development of high-precision electron microscopy methods for beam-sensitive materials.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Distributed-Memory DMRG via Sparse and Dense Parallel Tensor Contractions

The Density Matrix Renormalization Group (DMRG) algorithm is a powerful tool for solving eigenvalue problems to model quantum systems. DMRG relies on tensor contractions and dense linear algebra to compute properties of condensed matter physics systems. However, its efficient parallel implementation is challenging due to limited concurrency, large memory footprint, and tensor sparsity. We mitigate these problems by implementing two new parallel approaches that handle block sparsity arising in DMRG, via Cyclops, a distributed memory tensor contraction library. We benchmark their performance on two physical systems using the Blue Waters and Stampede2 supercomputers. Our DMRG performance is improved by up to 5.9X in runtime and 99X in processing rate over ITensor, at roughly comparable computational resource use. This enables higher accuracy calculations via larger tensors for quantum state approximation. We demonstrate that despite having limited concurrency, DMRG is weakly scalable with the use of efficient parallel tensor contraction mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Eigenstates hybridize on all length scales at the many-body localization transition

An interacting quantum system can transition from an ergodic to a many-body localized (MBL) phase under the presence of sufficiently large disorder. Both phases are radically different in their dynamical properties, which are characterized by highly excited eigenstates of the Hamiltonian. Each eigenstate can be characterized by the set of quantum numbers over the set of (local, in the MBL phase) integrals of motion of the system. In this work we study the evolution of the eigenstates of the disordered Heisenberg model as the disorder strength, $W$, is varied adiabatically. We focus on the probability that two `colliding' eigenstates hybridize as a function of both the range $R$ at which they differ as well as the strength of their hybridization. We find, in the MBL phase, that the probability of a colliding eigenstate hybridizing strongly at range $R$ decays as $Pr(R)\propto \exp [-R/η]$, with a length scale $η(W) = 1 / (B \log(W / W_c) )$ which diverges at the critical disorder strength $W_c$. This leads to range-invariance at the transition, suggesting the formation of resonating cat states at all ranges. This range invariance does not survive to the ergodic phase, where hybridization is exponentially more likely at large range, a fact that can be understood with simple combinatorial arguments. In fact, compensating for these combinatorial effects allows us to define an additional correlation length $ξ$ in the MBL phase which is in excellent agreement with previous works and which takes the critical value $1 / \log(2)$ at the transition, found in previous works to destabilize the MBL phase. Finally, we show that deep in the MBL phase hybridization is dominated by two-level collisions of eigenstates close in energy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Topology and the one-dimensional Kondo-Heisenberg model

The Kondo-Heinsberg chain is an interesting model of a strongly correlated system which has a broad superconducting state with pair-density wave (PDW) order. Some of us have recently proposed that this PDW state is a symmetry-protected topological (SPT) state, and the gapped spin sector of the model supports Majorana zero modes. In this work, we reexamine this problem using a combination of numeric and analytic methods. In extensive density matrix renormalization group calculations, we find no evidence of a topological ground state degeneracy or the previously proposed Majorana zero modes in the PDW phase of this model. This result motivated us to reexamine the original arguments for the existence of the Majorana zero modes. A careful analysis of the effective continuum field theory of the model shows that the Hilbert space of the spin sector of the theory does not contain any single Majorana fermion excitations. This analysis shows that the PDW state of the doped 1D Kondo-Heisenberg model is not an SPT with Majorana zero modes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Engineering Topological Models with a General-Purpose Symmetry-to-Hamiltonian Approach

Symmetry is at the heart of modern physics. Phases of matter are classified by symmetry breaking, topological phases are characterized by non-local symmetries, and point group symmetries are critical to our understanding of crystalline materials. Symmetries could then be used as a criterion to engineer quantum systems with targeted properties. Toward that end, we have developed a novel approach, the symmetric Hamiltonian construction (SHC), that takes as input symmetries, specified by integrals of motion or discrete symmetry transformations, and produces as output all local Hamiltonians consistent with these symmetries (see github.com/ClarkResearchGroup/qosy for our open-source code). This approach builds on the slow operator method [PRE 92, 012128]. We use our new approach to construct new Hamiltonians for topological phases of matter. Topological phases of matter are exotic quantum phases with potential applications in quantum computation. In this work, we focus on two types of topological phases of matter: superconductors with Majorana zero modes and $Z_2$ quantum spin liquids. In our first application of the SHC approach, we analytically construct a large and highly tunable class of superconducting Hamiltonians with Majorana zero modes with a given targeted spatial distribution. This result lays the foundation for potential new experimental routes to realizing Majorana fermions. In our second application, we find new $Z_2$ spin liquid Hamiltonians on the square and kagome lattices. These new Hamiltonians are not sums of commuting operators nor frustration-free and, when perturbed appropriately (in a way that preserves their $Z_2$ spin liquid behavior), exhibit level-spacing statistics that suggest non-integrability. This result demonstrates how our approach can automatically generate new spin liquid Hamiltonians with interesting properties not often seen in solvable models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bimodal entanglement entropy distribution in the many-body localization transition

We introduce the cut averaged entanglement entropy in disordered periodic spin chains and prove it to be a concave function of subsystem size for individual eigenstates. This allows us to identify the entanglement scaling as a function of subsystem size for individual states in inhomogeneous systems. Using this quantity, we probe the critical region between the many-body localized (MBL) and ergodic phases in finite systems. In the middle of the spectrum, we show evidence for bimodality of the entanglement distribution in the MBL critical region, finding both volume law and area law eigenstates over disorder realizations as well as within \emph{single disorder realizations}. The disorder averaged entanglement entropy in this region then scales as a volume law with a coefficient below its thermal value. We discover in the critical region, as we approach the thermodynamic limit, that the cut averaged entanglement entropy density falls on a one-parameter family of curves. Finally, we also show that without averaging over cuts the slope of the entanglement entropy \vs subsystem size can be negative at intermediate and strong disorder, caused by rare localized regions in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical evidence for a chiral spin liquid in the XXZ antiferromagnetic Heisenberg model on the kagome lattice at $m=\frac{2}{3}$ magnetization

We perform an exact diagonalization study of the spin-$1/2$ XXZ Heisenberg antiferromagnet on the kagome lattice at finite magnetization $m = \frac{2}{3}$ with an emphasis on the XY point ($J_z = 0$), and in the presence of a small chiral term. Recent analytic work by Kumar, Sun and Fradkin [Phys. Rev. B 90, 174409 (2014)] on the same model, using a newly developed flux attachment transformation, predicts a plateau at this value of the magnetization described by a chiral spin liquid (CSL) with a spin Hall conductance of $σ_{xy} = \frac{1}{2}$. Such a state is topological in nature, has a ground state degeneracy and exhibits fractional excitations. We analyze the degeneracy structure in the low energy manifold, identify the candidate topological states and use them to compute the modular matrices and Chern numbers all of which strongly agree with expected theoretical behavior for the $σ_{xy} = \frac{1}{2}$ CSL. We argue that the evidence suggests the CSL is robust even in the limit of zero external chirality.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Variational identification of a fractional Chern insulator in an extended Bose-Hubbard model

We study the extended Bose-Hubbard model on the square lattice at half filling as a function of next-nearest neighbor hopping amplitude and interaction strength. To variationally map out the phase diagram of this model, we develop a two-parameter family of wave-functions based on the parton construction which can describe both topological and broken symmetry phases on equal footing. In addition, our wave-functions resolve long standing issues with more conventional short-range Jastrow wave-functions. Using this variational ansatz, we show that a spontaneous time-reversal symmetry breaking fractional Chern insulator is energetically favored over a critical region between two superfluid phases. In verifying the properties of these parton wave-functions we exemplify a more robust way to identify topology through the Hall conductance.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Many-body Localization Transition in Rokhsar-Kivelson-type wave functions

We construct a family of many-body wave functions to study the many-body localization phase transition. The wave functions have a Rokhsar-Kivelson form, in which the weight for the configurations are chosen from the Gibbs weights of a classical spin glass model, known as the Random Energy Model, multiplied by a random sign structure to represent a highly excited state. These wave functions show a phase transition into an MBL phase. In addition, we see three regimes of entanglement scaling with subsystem size: scaling with entanglement corresponding to an infinite temperature thermal phase, constant scaling, and a sub-extensive scaling between these limits. Near the phase transition point, the fluctuations of the Rényi entropies are non-Gaussian. We find that Rényi entropies with different Rényi index transition into the MBL phase at different points and have different scaling behavior, suggesting a multifractal behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solving strongly correlated electron models on a quantum computer

One of the main applications of future quantum computers will be the simulation of quantum models. While the evolution of a quantum state under a Hamiltonian is straightforward (if sometimes expensive), using quantum computers to determine the ground state phase diagram of a quantum model and the properties of its phases is more involved. Using the Hubbard model as a prototypical example, we here show all the steps necessary to determine its phase diagram and ground state properties on a quantum computer. In particular, we discuss strategies for efficiently determining and preparing the ground state of the Hubbard model starting from various mean-field states with broken symmetry. We present an efficient procedure to prepare arbitrary Slater determinants as initial states and present the complete set of quantum circuits needed to evolve from these to the ground state of the Hubbard model. We show that, using efficient nesting of the various terms each time step in the evolution can be performed with just $\mathcal{O}(N)$ gates and $\mathcal{O}(\log N)$ circuit depth. We give explicit circuits to measure arbitrary local observables and static and dynamic correlation functions, both in the time and frequency domain. We further present efficient non-destructive approaches to measurement that avoid the need to re-prepare the ground state after each measurement and that quadratically reduce the measurement error.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gate count estimates for performing quantum chemistry on small quantum computers

As quantum computing technology improves and quantum computers with a small but non-trivial number of N > 100 qubits appear feasible in the near future the question of possible applications of small quantum computers gains importance. One frequently mentioned application is Feynman's original proposal of simulating quantum systems, and in particular the electronic structure of molecules and materials. In this paper, we analyze the computational requirements for one of the standard algorithms to perform quantum chemistry on a quantum computer. We focus on the quantum resources required to find the ground state of a molecule twice as large as what current classical computers can solve exactly. We find that while such a problem requires about a ten-fold increase in the number of qubits over current technology, the required increase in the number of gates that can be coherently executed is many orders of magnitude larger. This suggests that for quantum computation to become useful for quantum chemistry problems, drastic algorithmic improvements will be needed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stochastically Projecting Tensor Networks

We apply a series of projection techniques on top of tensor networks to compute energies of ground state wave functions with higher accuracy than tensor networks alone with minimal additional cost. We consider both matrix product states as well as tree tensor networks in this work. Building on top of these approaches, we apply fixed-node quantum Monte Carlo, Lanczos steps, and exact projection. We demonstrate these improvements for the triangular lattice Heisenberg model, where we capture up to 57 percent of the remaining energy not captured by the tensor network alone. We conclude by discussing further ways to improve our approach.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Searching for Topological Degeneracy in the Hubbard Model with Quantum Monte Carlo

Z2 spin liquids have topological order. One manifestation of this is that a Z2 spin liquid on a torus exhibits a four-fold degeneracy. Recent numerical evidence has argued for the existence of a spin liquid ground state in the Hubbard model on a honeycomb lattice near U approximately 4. The evidence for this claim involves the presence of a gapped state that lacks any identifiable order. This argument relies on being able to distinguish small order from no order which is notoriously difficult. In this paper we demonstrate an approach which uses quantum Monte Carlo to search for one of the key features that positively identify the topological spin liquid: the topological degeneracy. For any finite system, this topological degeneracy is split where the splitting decays exponentially with system size. We search for low lying states in the energy spectrum that could be identified as these topologically degenerate states. We show that, for system sizes up to N=162 sites, there is no evidence for these states being significantly below the first excited state giving evidence against the existence of a topological phase. We discuss the possible options for the Hubbard model on the honeycomb lattice in the absence of such degeneracy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effect of quantization on the FCIQMC sign problem

The sign problem in Full Configuration Interaction Quantum Monte Carlo (FCIQMC) without annihilation can be understood as an instability of the psi-particle population to the ground state of the matrix obtained by making all off-diagonal elements of the Hamiltonian negative. Such a matrix, and hence the sign problem, is basis dependent. In this paper we discuss the properties of a physically important basis choice: first versus second quantization. For a given choice of single-particle orbitals, we identify the conditions under which the fermion sign problem in the second quantized basis of antisymmetric Slater determinants is identical to the sign problem in the first quantized basis of unsymmetrized Hartree products. We also show that, when the two differ, the fermion sign problem is always less severe in the second quantized basis. This supports the idea that FCIQMC, even in the absence of annihilation, improves the sign problem relative to first quantized methods. Finally, we point out some theoretically interesting classes of Hamiltonians where first and second quantized sign problems differ, and others where they do not.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A second-quantized red herring in full configuration-interaction Monte Carlo

This paper deals with the sign problem in full configuration-interaction quantum Monte Carlo. After putting this article on the arxiv, it was pointed out to us that our argument applies to a number of model Hamiltonians we considered numerically but not to the most generic case. Please see arxiv: 1209.3044, where we have a corrected comprehensive discussion of the necessary and sufficient conditions for when the sign problem for a given Hamiltonian (and basis) differ between first and second quantization.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

FCI-QMC approach to the Fermi polaron

Finding the ground state of a fermionic Hamiltonian using quantum Monte Carlo is a very difficult problem, due to the Fermi sign problem. While still scaling exponentially, full configuration-interaction Monte Carlo (FCI-QMC) mitigates some of the exponential variance by allowing annihilation of noise -- whenever two walkers arrive at the same configuration with opposite signs, they are removed from the simulation. While FCI-QMC has been quite successful for quantum chemistry problems, its application to problems in condensed systems has been limited. In this paper, we apply FCI-QMC to the Fermi polaron problem, which provides an ideal test-bed for improving the algorithm. In its simplest form, FCI-QMC is unstable for even a fairly small system sizes. However, with a series of algorithmic improvements, we are able to significantly increase its effectiveness. We modify fixed node QMC to work in these systems, introduce a well chosen importance sampled trial wave function, a partial node approximation, and a variant of released node. Finally, we develop a way to perform FCI-QMC directly in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-equilibrium dynamic critical scaling of the quantum Ising chain

We solve for the time-dependent finite-size scaling functions of the 1D transverse-field Ising chain during a linear-in-time ramp of the field through the quantum critical point. We then simulate Mott-insulating bosons in a tilted potential, an experimentally-studied system in the same equilibrium universality class, and demonstrate that universality holds for the dynamics as well. We find qualitatively athermal features of the scaling functions, such as negative spin correlations, and show that they should be robustly observable within present cold atom experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Path Integral Monte Carlo Simulation of the Warm-Dense Homogeneous Electron Gas

We perform calculations of the {3D} finite-temperature homogeneous electron gas (HEG) in the warm-dense regime ({r_{s} \equiv (3/4πn)^{1/3}a_{B}^{- 1} = 1.0- 40.0} and {Θ\equiv T/T_{F} = 0.0625- 8.0}) using restricted path integral Monte Carlo (RPIMC). Precise energies, pair correlation functions, and structure factors are obtained. For all densities, we find a significant discrepancy between the ground state parameterized local density approximation (LDA) and our results around {T_{F}} . These results can be used as a benchmark for improved functionals, as well as input for orbital-free DFT formulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Computing the energy of a water molecule using MultiDeterminants: A simple, efficient algorithm

Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave-functions are critical to ascertaining new physics. One such wave function is the multiSlater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wavefunctions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally we implement this method and use it to compute the ground state energy of a water molecule.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Improved Scaling for Quantum Monte Carlo on Insulators

Quantum Monte Carlo (QMC) methods are often used to calculate properties of many body quantum systems. The main cost of many QMC methods, for example the variational Monte Carlo (VMC) method, is in constructing a sequence of Slater matrices and computing the ratios of determinants for successive Slater matrices. Recent work has improved the scaling of constructing Slater matrices for insulators so that the cost of constructing Slater matrices in these systems is now linear in the number of particles, whereas computing determinant ratios remains cubic in the number of particles. With the long term aim of simulating much larger systems, we improve the scaling of computing the determinant ratios in the VMC method for simulating insulators by using preconditioned iterative solvers. The main contribution of this paper is the development of a method to efficiently compute for the Slater matrices a sequence of preconditioners that make the iterative solver converge rapidly. This involves cheap preconditioner updates, an effective reordering strategy, and a cheap method to monitor instability of ILUTP preconditioners. Using the resulting preconditioned iterative solvers to compute determinant ratios of consecutive Slater matrices reduces the scaling of QMC algorithms from O(n^3) per sweep to roughly O(n^2), where n is the number of particles, and a sweep is a sequence of n steps, each attempting to move a distinct particle. We demonstrate experimentally that we can achieve the improved scaling without increasing statistical errors. Our results show that preconditioned iterative solvers can dramatically reduce the cost of VMC for large(r) systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-equilibrium dynamics of Bosonic Mott insulators in an electric field

We study the non-equilibrium dynamics of one-dimensional Mott insulating bosons in the presence of a tunable effective electric field E which takes the system across a quantum critical point (QCP) separating a disordered and a translation symmetry broken ordered phase. We provide an exact numerical computation of the residual energy Q, the log-fidelity F, the excess defect density D, and the order parameter correlation function for a linear-in-time variation of E with a rate v. We discuss the temporal and spatial variation of these quantities for a range of v and for finite system sizes as relevant to realistic experimental setups [J. Simon et al., Nature 472, 307 (2011)]. We show that in finite-sized systems Q, F, and D obey Kibble-Zurek scaling, and suggest further experiments within this setup to test our theory.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hexatic and mesoscopic phases in the 2D quantum Coulomb system

We study the Wigner crystal melting in a two dimensional quantum system of particles interacting via the 1/r Coulomb potential. We use quantum Monte Carlo methods to calculate its phase diagram, locate the Wigner crystal region, and analyze its instabilities towards the liquid phase. We discuss the role of quantum effects in the critical behavior of the system, and compare our numerical results with the classical theory of melting, and the microemulsion theory of frustrated Coulomb systems. We find a Pomeranchuk effect much larger then in solid helium. In addition, we find that the exponent for the algebraic decay of the hexatic phase differs significantly from the Kosterilitz-Thouless theory of melting. We search for the existence of mesoscopic phases and find evidence of metastable bubbles but no mesoscopic phase that is stable in equilibrium.

0 citations0 reviewsNo code linkedOpen access