Source author record

William J. Porter

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

5works
7topics
2close collaborators

Actions

Connect this record

Log in to claim

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

5 published item(s)

preprint2016arXiv

Entanglement, noise, and the cumulant expansion

We put forward a simpler and improved variation of a recently proposed method to overcome the signal-to-noise problem found in Monte Carlo calculations of the entanglement entropy of interacting fermions. The present method takes advantage of the approximate lognormal distributions that characterize the signal-to-noise properties of other approaches. In addition, we show that a simple rewriting of the formalism allows circumvention of the inversion of the restricted one-body density matrix in the calculation of the $n$-th Rényi entanglement entropy for $n>2$. We test our technique by implementing it in combination with the hybrid Monte Carlo algorithm and calculating the $n=2,3,4, \dots, 10$ Rényi entropies of the 1D attractive Hubbard model. We use that data to extrapolate to the von Neumann ($n=1$) and $n\to\infty$ cases.

preprint2016arXiv

The entanglement spectrum and Rényi entropies of non-relativistic conformal fermions

We characterize non-perturbatively the Rényi entropies of degree n=2,3,4, and 5 of three-dimensional, strongly coupled many-fermion systems in the scale-invariant regime of short interaction range and large scattering length, i.e. in the unitary limit. We carry out our calculations using lattice methods devised recently by us. Our results show the effect of strong pairing correlations on the entanglement entropy, which modify the sub-leading behavior for large subsystem sizes (as characterized by the dimensionless parameter x=kF L_A, where kF is the Fermi momentum and L_A the linear subsystem size), but leave the leading order unchanged relative to the non-interacting case. Moreover, we find that the onset of the sub-leading asymptotic regime is at surprisingly small x=2-4. We provide further insight into the entanglement properties of this system by analyzing the spectrum of the entanglement Hamiltonian of the two-body problem from weak to strong coupling. The low-lying entanglement spectrum displays clear features as the strength of the coupling is varied, such as eigenvalue crossing, a sharp change in the Schmidt gap, and scale invariance at unitarity. Beyond the low-lying component, the spectrum appears as a quasi-continuum distribution, for which we present a statistical characterization; we find, in particular, that the mean shifts to infinity as the coupling is turned off, which indicates that that part of the spectrum represents non-perturbative contributions to the entanglement Hamiltonian. In contrast, the low-lying entanglement spectrum evolves to finite values in the noninteracting limit. The scale invariance of the unitary regime guarantees that our results are universal features intrinsic to 3D quantum mechanics and represent a well-defined prediction for ultracold atom experiments, which were recently shown to have direct access to the entanglement entropy.

preprint2015arXiv

A hybrid Monte Carlo approach to the entanglement entropy of interacting fermions

The Monte Carlo calculation of Rényi entanglement entropies $S^{}_n$ of interacting fermions suffers from a well-known signal-to-noise problem, even for a large number of situations in which the infamous sign problem is absent. A few methods have been proposed to overcome this issue, such as ensemble switching and the use of auxiliary partition-function ratios. Here, we present an approach that builds on the recently proposed free-fermion decomposition method; it incorporates entanglement in the probability measure in a natural way; it takes advantage of the hybrid Monte Carlo algorithm (an essential tool in lattice quantum chromodynamics and other gauge theories with dynamical fermions); and it does not suffer from noise problems. This method displays no sign problem for the same cases as other approaches and is therefore useful for a wide variety of systems. As a proof of principle, we calculate $S_2^{}$ for the one-dimensional, half-filled Hubbard model and compare with results from exact diagonalization and the free-fermion decomposition method.

preprint2015arXiv

Convexity of the entanglement energy of SU($2N$)-symmetric fermions with attractive interactions

The positivity of the probability measure of attractively interacting systems of $2N$-component fermions enables the derivation of an exact convexity property for the ground-state energy of such systems. Using analogous arguments, applied to path-integral expressions for the $n$-th Rényi entanglement entropy $S^{}_n$ derived recently, we prove non-perturbative analytic relations for the entanglement energies $\mathcal E^{}_n$ of those systems defined via $\mathcal{E}^{}_n \equiv \frac{n-1}{n}T S^{}_n + F$ where $β= 1/T$ is the extent of the imaginary time direction and $-βF = \ln \mathcal Z$ where $\mathcal Z$ is the partition sum appropriate to the temperature. These relations are valid for all sub-system sizes, particle numbers and dimensions, and in arbitrary external trapping potentials.

preprint2015arXiv

Hard-Wall and Non-Uniform Lattice Monte Carlo Approaches to One-Dimensional Fermi Gases in a Harmonic Trap

We present in detail two variants of the lattice Monte Carlo method aimed at tackling systems in external trapping potentials: a uniform-lattice approach with hard-wall boundary conditions, and a non-uniform Gauss-Hermite lattice approach. Using those two methods, we compute the ground-state energy and spatial density profile for systems of N=4 - 8 harmonically trapped fermions in one dimension. From the favorable comparison of both energies and density profiles (particularly in regions of low density), we conclude that the trapping potential is properly resolved by the hard-wall basis. Our work paves the way to higher dimensions and finite temperature analyses, as calculations with the hard-wall basis can be accelerated via fast Fourier transforms, the cost of unaccelerated methods is otherwise prohibitive due to the unfavorable scaling with system size.