Showing 1–50 of 66 results for author: Del Maestro, A

SmoQyDEAC.jl: A differential evolution package for the analytic continuation of imaginary time correlation functions

Authors: James Neuhaus, Nathan S. Nichols, Debshikha Banerjee, Benjamin Cohen-Stead, Thomas A. Maier, Adrian Del Maestro, Steven Johnston

Abstract: We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara… ▽ More We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara frequency axes, and treatment of the covariance error in the input data. This paper presents an overview of the DEAC algorithm and the features implemented in the SmoQyDEAC.jl. It also provides detailed benchmarks of the package's output against the popular maximum entropy and stochastic analytic continuation methods. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDEAC.jl or installed using the Julia package manager. The online documentation, including examples, can be accessed at https://smoqysuite.github.io/SmoQyDEAC.jl/stable/. △ Less

Submitted 30 September, 2024; v1 submitted 5 July, 2024; originally announced July 2024.

Direct Visualization of Disorder Driven Electronic Liquid Crystal Phases in Dirac Nodal Line Semimetal GdSbTe

Authors: Balaji Venkatesan, Syu-You Guan, Jen-Te Chang, Shiang-Bin Chiu, Po-Yuan Yang, Chih-Chuan Su, Tay-Rong Chang, Kalaivanan Raju, Raman Sankar, Somboon Fongchaiya, Ming-Wen Chu, Chia-Seng Chang, Guoqing Chang, Hsin Lin, Adrian Del Maestro, Ying-Jer Kao, Tien-Ming Chuang

Abstract: Electronic liquid crystal (ELC) phases are spontaneous symmetry breaking states believed to arise from strong electron correlation in quantum materials such as cuprates and iron pnictides. Here, we report a direct observation of ELC phases in a Dirac nodal line (DNL) semimetal GdSbxTe2-x. Electronic nanostructures consisting of incommensurate smectic charge modulation and intense local nematic ord… ▽ More Electronic liquid crystal (ELC) phases are spontaneous symmetry breaking states believed to arise from strong electron correlation in quantum materials such as cuprates and iron pnictides. Here, we report a direct observation of ELC phases in a Dirac nodal line (DNL) semimetal GdSbxTe2-x. Electronic nanostructures consisting of incommensurate smectic charge modulation and intense local nematic order are visualized by using spectroscopic imaging - scanning tunneling microscopy. As topological materials with symmetry protected Dirac or Weyl fermions are mostly weakly correlated, the discovery of such ELC phases are anomalous and raise questions on the origin of their emergence. Specifically, we demonstrate how chemical substitution generates these symmetry breaking phases before the system undergoes a charge density wave - orthorhombic structural transition. We further show how dopants can induce nematicity via quasiparticle scattering interference. Our results highlight the importance of impurities in realizing ELC phases and present a new material platform for exploring the interplay among quenched disorder, topology and electron correlation. △ Less

Submitted 7 May, 2024; v1 submitted 29 February, 2024; originally announced February 2024.

Efficient simulations of Hartree--Fock equations by an accelerated gradient descent method

Authors: Y. Ohno, A. Del Maestro, T. I. Lakoba

Abstract: We develop convergence acceleration procedures that enable a gradient descent-type iteration method to efficiently simulate Hartree--Fock equations for atoms interacting both with each other and with an external potential. Our development focuses on three aspects: (i) optimization of a parameter in the preconditioning operator; (ii) adoption of a technique that eliminates the slowest-decaying mode… ▽ More We develop convergence acceleration procedures that enable a gradient descent-type iteration method to efficiently simulate Hartree--Fock equations for atoms interacting both with each other and with an external potential. Our development focuses on three aspects: (i) optimization of a parameter in the preconditioning operator; (ii) adoption of a technique that eliminates the slowest-decaying mode to the case of many equations (describing many atoms); and (iii) a novel extension of the above technique that allows one to eliminate multiple modes simultaneously. We illustrate performance of the numerical method for the 2D model of the first layer of helium atoms above a graphene sheet. We demonstrate that incorporation of aspects (i) and (ii) above into the ``plain" gradient descent method accelerates it by at least two orders of magnitude, and often by much more. Aspect (iii) -- a multiple-mode elimination -- may bring further improvement to the convergence rate compared to aspect (ii), the single-mode elimination. Both single- and multiple-mode elimination techniques are shown to significantly outperform the well-known Anderson Acceleration. We believe that our acceleration techniques can also be gainfully employed by other numerical methods, especially those handling hard-core-type interaction potentials. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: main text (39 pages); supplement appended (7 pages)

Disentangling the Physics of the Attractive Hubbard Model via the Accessible and Symmetry-Resolved Entanglement Entropies

Authors: Tong Shen, Hatem Barghathi, Adrian Del Maestro, Brenda Rubenstein

Abstract: The complicated ways in which electrons interact in many-body systems such as molecules and materials have long been viewed through the lens of local electron correlation and associated correlation functions. However, quantum information science has demonstrated that more global diagnostics of quantum states, like the entanglement entropy, can provide a complementary and clarifying lens on electro… ▽ More The complicated ways in which electrons interact in many-body systems such as molecules and materials have long been viewed through the lens of local electron correlation and associated correlation functions. However, quantum information science has demonstrated that more global diagnostics of quantum states, like the entanglement entropy, can provide a complementary and clarifying lens on electronic behavior. One particularly useful measure that can be used to distinguish between quantum and classical sources of entanglement is the accessible entanglement, the entanglement available as a quantum resource for systems subject to conservation laws, such as fixed particle number, due to superselection rules. In this work, we introduce an algorithm and demonstrate how to compute accessible and symmetry-resolved entanglements for interacting fermion systems. This is accomplished by combining an incremental version of the swap algorithm with a recursive Auxiliary Field Quantum Monte Carlo algorithm recently developed by the authors. We apply these tools to study the pairing and charge density waves exhibited in the paradigmatic attractive Hubbard model via entanglement. We find that the particle and spin symmetry-resolved entanglements and their related full probability distribution functions show very clear - and unique - signatures of the underlying electronic behavior even when those features are less pronounced in more conventional correlation functions. Overall, this work provides a systematic means of characterizing the entanglement within quantum systems that can grant a deeper understanding of the complicated electronic behavior that underlies quantum phase transitions and crossovers in many-body systems. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 19 pages, 9 figures

Journal ref: Phys. Rev. B 109, 195119 (2024)

Effects of substrate corrugation during helium adsorption on graphene in the grand canonical ensemble

Authors: Gage Erwin, Adrian Del Maestro

Abstract: Adsorption of 4He on atomically flat substrates such as graphene provides a route towards the engineering of low dimensional quantum phases including superfluids and strongly interacting insulators. In this study, we explore the effects of graphene corrugation on the helium adsorption process via quantum Monte Carlo simulations in the grand canonical ensemble. We utilize an empirical adsorption po… ▽ More Adsorption of 4He on atomically flat substrates such as graphene provides a route towards the engineering of low dimensional quantum phases including superfluids and strongly interacting insulators. In this study, we explore the effects of graphene corrugation on the helium adsorption process via quantum Monte Carlo simulations in the grand canonical ensemble. We utilize an empirical adsorption potential based on the superposition of individual helium-carbon interactions and systematically control corrugation, from a smooth membrane to the fully rough potential, via the implementation of a cutoff in reciprocal space. The results highlight the importance of using a fully corrugated potential to understand the plethora of commensurate and incommensurate phases in the first adsorbed layer. Surprisingly, some residual effects of corrugation are still present before and during the promotion and onset of a second adsorbed layer, where a smooth adsorption can lead to enhanced particle fluctuations in a helium-interaction dominated regime. △ Less

Submitted 29 April, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: 12 pages, 8 figures (added additional system sizes). Submitted to the JLTP Special Issue: Quantum Fluids and Solids (QFS2023). For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-GrapheneCorrugation

Journal ref: J. Low Temp. Phys. 215, 525 (2024)

Topological and magnetic properties of the interacting Bernevig-Hughes-Zhang model

Authors: Rahul Soni, Harini Radhakrishnan, Bernd Rosenow, Gonzalo Alvarez, Adrian Del Maestro

Abstract: We investigate the effects of electronic correlations on the Bernevig-Hughes-Zhang model using the real-space density matrix renormalization group (DMRG) algorithm. We introduce a method to probe topological phase transitions in systems with strong correlations using DMRG, substantiated by an unsupervised machine learning methodology that analyzes the orbital structure of the real-space edges. Inc… ▽ More We investigate the effects of electronic correlations on the Bernevig-Hughes-Zhang model using the real-space density matrix renormalization group (DMRG) algorithm. We introduce a method to probe topological phase transitions in systems with strong correlations using DMRG, substantiated by an unsupervised machine learning methodology that analyzes the orbital structure of the real-space edges. Including the full multi-orbital Hubbard interaction term, we construct a phase diagram as a function of a gap parameter ($m$) and the Hubbard interaction strength ($U$) via exact DMRG simulations on $N\times 4$ cylinders. Our analysis confirms that the topological phase persists in the presence of interactions, consistent with previous studies, but it also reveals an intriguing phase transition from a paramagnetic to a stripey antiferromagnetic topological insulator. The combination of the magnetic structure factor, strength of magnetic moments, and the orbitally resolved density, provides real-space information on both topology and magnetism in a strongly correlated system. △ Less

Submitted 14 September, 2024; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: Main text: 13 pages, 15 figures. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-interactingBHZmodel

A Scaling Function for the Particle Entanglement Entropy of Fermions

Authors: Harini Radhakrishnan, Matthias Thamm, Hatem Barghathi, Bernd Rosenow, Adrian Del Maestro

Abstract: Entanglement entropy under a particle bipartition provides complementary information to mode entanglement as it is sensitive to interactions and particle statistics at leading order and does not depend on any externally imposed length scale. In this paper, we investigate the particle entanglement entropy in a system of $N$ interacting spinless lattice fermions in one spatial dimension by combining… ▽ More Entanglement entropy under a particle bipartition provides complementary information to mode entanglement as it is sensitive to interactions and particle statistics at leading order and does not depend on any externally imposed length scale. In this paper, we investigate the particle entanglement entropy in a system of $N$ interacting spinless lattice fermions in one spatial dimension by combining bosonization techniques with exact and approximate numerical methods. We introduce a general scaling form for the fermionic particle entanglement entropy captured by a shape function that enters as a extensive interaction induced correction to a known free fermion result. A general asymptotic expansion in the total number of particles demonstrates that its form is robust for different values of the Rényi index and highlights how quantum correlations are encoded in the $n$-particle density matrix of a pure many-body quantum state. △ Less

Submitted 15 August, 2023; v1 submitted 17 February, 2023; originally announced February 2023.

Comments: 33 pages, 12 figures. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-scalingParticleEE

Journal ref: J. Stat. Mech. (2023) 083101

Reduction of Autocorrelation Times in Lattice Path Integral Quantum Monte Carlo via Direct Sampling of the Truncated Exponential Distribution

Authors: Emanuel Casiano-Diaz, Kipton Barros, Ying Wai Li, Adrian Del Maestro

Abstract: In Monte Carlo simulations, proposed configurations are accepted or rejected according to an acceptance ratio, which depends on an underlying probability distribution and an a priori sampling probability. By carefully selecting the probability distribution from which random variates are sampled, simulations can be made more efficient, by virtue of an autocorrelation time reduction. In this paper,… ▽ More In Monte Carlo simulations, proposed configurations are accepted or rejected according to an acceptance ratio, which depends on an underlying probability distribution and an a priori sampling probability. By carefully selecting the probability distribution from which random variates are sampled, simulations can be made more efficient, by virtue of an autocorrelation time reduction. In this paper, we illustrate how to directly sample random variates from a two dimensional truncated exponential distribution. We show that our direct sampling approach converges faster to the target distribution compared to rejection sampling. The direct sampling of one and two dimensional truncated exponential distributions is then applied to a recent Path Integral Monte Carlo (PIMC) algorithm for the simulation of Bose-Hubbard lattice models at zero temperature. The new sampling method results in improved acceptance ratios and reduced autocorrelation times of estimators, providing an effective speed up of the simulation. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: ~8 pages, 8 figures. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-truncExponSampling

Quantum Atomic Matter Near Two-Dimensional Materials in Microgravity

Authors: Adrian Del Maestro, Sang Wook Kim, Nicholas P. Bigelow, Robert J. Thompson, Valeri N. Kotov

Abstract: Novel two-dimensional (2D) atomically flat materials, such as graphene and transition-metal dichalcogenides, exhibit unconventional Dirac electronic spectra. We propose to effectively engineer their interactions with cold atoms in microgravity, leading to a synergy between complex electronic and atomic collective quantum phases and phenomena. Dirac materials are susceptible to manipulation and qua… ▽ More Novel two-dimensional (2D) atomically flat materials, such as graphene and transition-metal dichalcogenides, exhibit unconventional Dirac electronic spectra. We propose to effectively engineer their interactions with cold atoms in microgravity, leading to a synergy between complex electronic and atomic collective quantum phases and phenomena. Dirac materials are susceptible to manipulation and quantum engineering via changes in their electronic properties by application of strain, doping with carriers, adjustment of their dielectric environment, etc. Consequently the interaction of atoms with such materials, namely the van der Waals / Casimir-Polder interaction, can be effectively manipulated, leading to the potential observation of physical effects such as Quantum Reflection off atomically thin materials and confined Bose-Einstein Condensate (BEC) frequency shifts. △ Less

Submitted 18 August, 2023; v1 submitted 1 January, 2023; originally announced January 2023.

Comments: 11 pages, 3 figures; discussion and references added

Journal ref: Quantum Science and Technology 8, 044002 (2023)

A Stable, Recursive Auxiliary Field Quantum Monte Carlo Algorithm in the Canonical Ensemble: Applications to Thermometry and the Hubbard Model

Authors: Tong Shen, Hatem Barghathi, Jiangyong Yu, Adrian Del Maestro, Brenda Rubenstein

Abstract: Many experimentally-accessible, finite-sized interacting quantum systems are most appropriately described by the canonical ensemble of statistical mechanics. Conventional numerical simulation methods either approximate them as being coupled to a particle bath, or use projective algorithms which may suffer from non-optimal scaling with system size or large algorithmic prefactors. In this paper, we… ▽ More Many experimentally-accessible, finite-sized interacting quantum systems are most appropriately described by the canonical ensemble of statistical mechanics. Conventional numerical simulation methods either approximate them as being coupled to a particle bath, or use projective algorithms which may suffer from non-optimal scaling with system size or large algorithmic prefactors. In this paper, we introduce a highly stable, recursive Auxiliary Field Quantum Monte Carlo approach that can directly simulate systems in the canonical ensemble. We apply the method to the fermion Hubbard model in one and two spatial dimensions in a regime known to exhibit a significant "sign" problem and find improved performance over existing approaches including rapid convergence to ground state expectation values. The effects of excitations above the ground state are quantified using an estimator-agnostic approach including studying the temperature dependence of the purity and overlap fidelity of the canonical and grand canonical density matrices. As an important application, we show that thermometry approaches often exploited in ultra-cold atoms that employ an analysis of the velocity distribution in the grand canonical ensemble may be subject to errors leading to an under-estimation of extracted temperatures with respect to the Fermi temperature. △ Less

Submitted 19 December, 2022; v1 submitted 16 December, 2022; originally announced December 2022.

Comments: 14 pages, 7 figures, and Supplemental Information with key derivations. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-CanEnsAFQMC

Journal ref: Phys. Rev. E 107, 055302 (2023)

Strain-induced superfluid transition for atoms on graphene

Authors: Sang Wook Kim, Mohamed Elsayed, Nathan S. Nichols, Taras Lakoba, Juan Vanegas, Carlos Wexler, Valeri N. Kotov, Adrian Del Maestro

Abstract: Bosonic atoms deposited on atomically thin substrates represent a playground for exotic quantum many-body physics due to the highly-tunable, atomic-scale nature of the interaction potentials. The ability to engineer strong interparticle interactions can lead to the emergence of complex collective atomic states of matter, not possible in the context of dilute atomic gases confined in optical lattic… ▽ More Bosonic atoms deposited on atomically thin substrates represent a playground for exotic quantum many-body physics due to the highly-tunable, atomic-scale nature of the interaction potentials. The ability to engineer strong interparticle interactions can lead to the emergence of complex collective atomic states of matter, not possible in the context of dilute atomic gases confined in optical lattices. While it is known that the first layer of adsorbed helium on graphene is permanently locked into a solid phase, we show by a combination of quantum Monte Carlo and mean-field techniques, that simple isotropic graphene lattice expansion effectively unlocks a large variety of two-dimensional ordered commensurate, incommensurate, cluster atomic solid, and superfluid states for adsorbed atoms. It is especially significant that an atomically thin superfluid phase of matter emerges under experimentally feasible strain values, with potentially supersolid phases in close proximity on the phase diagram. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: Main text: 12 pages and 5 figures. Supplementary Information: 3 pages and 3 figures. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-Superfluid4HeStrainGraphene

Journal ref: Phys. Rev. B 109, 064512 (2024)

Robust charge-density wave correlations in the electron-doped single-band Hubbard model

Authors: Peizhi Mai, Nathan S. Nichols, Seher Karakuzu, Feng Bao, Adrian Del Maestro, Thomas A. Maier, Steven Johnston

Abstract: There is growing evidence that the hole-doped single-band Hubbard and $t$-$J$ models do not have a superconducting ground state reflective of the high-temperature cuprate superconductors but instead have striped spin- and charge-ordered ground states. Nevertheless, it is proposed that these models may still provide an effective low-energy model for electron-doped materials. Here we study the finit… ▽ More There is growing evidence that the hole-doped single-band Hubbard and $t$-$J$ models do not have a superconducting ground state reflective of the high-temperature cuprate superconductors but instead have striped spin- and charge-ordered ground states. Nevertheless, it is proposed that these models may still provide an effective low-energy model for electron-doped materials. Here we study the finite temperature spin and charge correlations in the electron-doped Hubbard model using quantum Monte Carlo dynamical cluster approximation calculations and contrast their behavior with those found on the hole-doped side of the phase diagram. We find evidence for a charge modulation with both checkerboard and unidirectional components decoupled from any spin-density modulations. These correlations are inconsistent with a weak-coupling description based on Fermi surface nesting, and their doping dependence agrees qualitatively with resonant inelastic x-ray scattering measurements. Our results provide evidence that the single-band Hubbard model describes the electron-doped cuprates. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Spin Squeezing as a Probe of Emergent Quantum Orders

Authors: Ilija K. Nikolov, Stephen Carr, Adrian G. Del Maestro, Chandrasekhar Ramanathan, Vesna F. Mitrović

Abstract: Nuclear magnetic resonance (NMR) experiments can reveal local properties in materials, but are often limited by the low signal-to-noise ratio. Spin squeezed states have an improved resolution below the Heisenberg limit in one of the spin components, and have been extensively used to improve the sensitivity of atomic clocks, for example. Interacting and entangled spin ensembles with non-linear coup… ▽ More Nuclear magnetic resonance (NMR) experiments can reveal local properties in materials, but are often limited by the low signal-to-noise ratio. Spin squeezed states have an improved resolution below the Heisenberg limit in one of the spin components, and have been extensively used to improve the sensitivity of atomic clocks, for example. Interacting and entangled spin ensembles with non-linear coupling are a natural candidate for implementing squeezing. Here, we propose measurement of the spin-squeezing parameter that itself can act as a local probe of emergent orders in quantum materials. In particular, we demonstrate how to investigate an anisotropic electric field gradient via its coupling to the nuclear quadrupole moment. While squeezed spin states are pure, the squeezing parameter can be estimated for both pure and mixed states. We evaluate the range of fields and temperatures for which a thermal-equilibrium state is sufficient to improve the resolution in an NMR experiment and probe relevant parameters of the quadrupole Hamiltonian, including its anisotropy. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Report number: v38, 011149

Multi-modal spectroscopy of order parameter distributions

Authors: Stephen Carr, Ilija K. Nikolov, Rong Cong, Adrian Del Maestro, Chandrasekhar Ramanathan, V. F. Mitrović

Abstract: We present a multi-modal spectroscopic paradigm that enables independent measurement of charge and spin degrees of freedom (DOF) in strongly correlated materials. This spin-based technique probes symmetry-specific Hamiltonian parameters by analyzing how the time delay between applied pulses ($τ$) affects the response. We demonstrate ways in which charge DOF that couple through the quadrupolar inte… ▽ More We present a multi-modal spectroscopic paradigm that enables independent measurement of charge and spin degrees of freedom (DOF) in strongly correlated materials. This spin-based technique probes symmetry-specific Hamiltonian parameters by analyzing how the time delay between applied pulses ($τ$) affects the response. We demonstrate ways in which charge DOF that couple through the quadrupolar interaction (inversion symmetric) can be independently measured even in the presence of large magnetic noise (inversion asymmetric). The method quantifies both the strength of the interactions and its distribution (noise). We provide protocols to directly and independently measure the distribution of interaction strengths, even when the average value of the interaction is zero. By independently measuring distributions of different forms of disorder, this methodology can elucidate which microscopic symmetry drives a phase transition. We discuss potential applications to study complex phase transitions in strongly interacting quantum materials. △ Less

Submitted 6 March, 2023; v1 submitted 23 August, 2022; originally announced August 2022.

Comments: 14 pages, 6 figures

A Path Integral Ground State Monte Carlo Algorithm for Entanglement of Lattice Bosons

Authors: Emanuel Casiano-Diaz, C. M. Herdman, Adrian Del Maestro

Abstract: A ground state path integral quantum Monte Carlo algorithm is introduced that allows for the study of entanglement in lattice bosons at zero temperature. The Rényi entanglement entropy between spatial subregions is explored across the phase diagram of the one dimensional Bose-Hubbard model for systems consisting of up to $L=256$ sites at unit-filling without any restrictions on site occupancy, far… ▽ More A ground state path integral quantum Monte Carlo algorithm is introduced that allows for the study of entanglement in lattice bosons at zero temperature. The Rényi entanglement entropy between spatial subregions is explored across the phase diagram of the one dimensional Bose-Hubbard model for systems consisting of up to $L=256$ sites at unit-filling without any restrictions on site occupancy, far beyond the reach of exact diagonalization. The favorable scaling of the algorithm is demonstrated through a further measurement of the Rényi entanglement entropy at the two dimensional superfluid-insulator critical point for large system sizes, confirming the existence of the expected entanglement boundary law in the ground state. The Rényi estimator is extended to measure the symmetry resolved entanglement that is operationally accessible as a resource for experimentally relevant lattice gases with fixed total particle number. △ Less

Submitted 17 January, 2023; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: ~50 pages, 21 figures. Updated refs, figures, and discussion. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-pigsfli

Journal ref: SciPost Phys. 14, 054 (2023)

One-particle entanglement for one dimensional spinless fermions after an interaction quantum quench

Authors: Matthias Thamm, Harini Radhakrishnan, Hatem Barghathi, Bernd Rosenow, Adrian Del Maestro

Abstract: Particle entanglement provides information on quantum correlations in systems of indistinguishable particles. Here, we study the one particle entanglement entropy for an integrable model of spinless, interacting fermions both at equilibrium and after an interaction quantum quench. Using both large scale exact diagonalization and time dependent density matrix renormalization group calculations, we… ▽ More Particle entanglement provides information on quantum correlations in systems of indistinguishable particles. Here, we study the one particle entanglement entropy for an integrable model of spinless, interacting fermions both at equilibrium and after an interaction quantum quench. Using both large scale exact diagonalization and time dependent density matrix renormalization group calculations, we numerically compute the one body reduced density matrix for the J-V model, as well as its post-quench dynamics. We include an analysis of the fermionic momentum distribution, showcasing its time evolution after a quantum quench. Our numerical results, extrapolated to the thermodynamic limit, can be compared with field theoretic bosonization in the Tomonaga-Luttinger liquid regime. Excellent agreement is obtained using an interaction cutoff that can be determined uniquely in the ground state. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: 20 pages, 16 figures

Journal ref: Phys. Rev. B 106, 165116 (2022)

Experimental Realization of One Dimensional Helium

Authors: Adrian Del Maestro, Nathan S. Nichols, Timothy R. Prisk, Garfield Warren, Paul E. Sokol

Abstract: The realization of experimental platforms exhibiting one dimensional (1D) quantum phenomena has been elusive, due to their inherent lack of stability, with a few notable exceptions including spin chains, carbon nanotubes and ultracold low-density gasses. The difficulty of such systems in exhibiting long range order is integral to their effective description in terms of the Tomonaga-Luttinger liqui… ▽ More The realization of experimental platforms exhibiting one dimensional (1D) quantum phenomena has been elusive, due to their inherent lack of stability, with a few notable exceptions including spin chains, carbon nanotubes and ultracold low-density gasses. The difficulty of such systems in exhibiting long range order is integral to their effective description in terms of the Tomonaga-Luttinger liquid theory. Recently, it has been proposed that bosonic superfluid 4He could realize a 1D quantum system beyond the Luttinger liquid paradigm. Here we describe an experimental observation of this behavior using nanoengineering by preplating a porous material with a noble gas to enhance dimensional reduction. The resulting excitations of the confined 4He are qualitatively different than 3D and 2D superfluid helium, and can be analyzed in terms of a mobile impurity in a Luttinger liquid allowing for the characterization of the emergent quantum liquid. The confined helium system offers the possibility of tuning via pressure - from weakly interacting, all the way to the super Tonks-Girardeau gas of strongly interacting hard-core particles. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 10 pages, 5 figures. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-preplated-nanopores-scattering

Journal ref: Nature Comm. 13, 3168 (2022)

A Group-Equivariant Autoencoder for Identifying Spontaneously Broken Symmetries

Authors: Devanshu Agrawal, Adrian Del Maestro, Steven Johnston, James Ostrowski

Abstract: We introduce the group-equivariant autoencoder (GE-autoencoder) -- a deep neural network (DNN) method that locates phase boundaries by determining which symmetries of the Hamiltonian have spontaneously broken at each temperature. We use group theory to deduce which symmetries of the system remain intact in all phases, and then use this information to constrain the parameters of the GE-autoencoder… ▽ More We introduce the group-equivariant autoencoder (GE-autoencoder) -- a deep neural network (DNN) method that locates phase boundaries by determining which symmetries of the Hamiltonian have spontaneously broken at each temperature. We use group theory to deduce which symmetries of the system remain intact in all phases, and then use this information to constrain the parameters of the GE-autoencoder such that the encoder learns an order parameter invariant to these ``never-broken'' symmetries. This procedure produces a dramatic reduction in the number of free parameters such that the GE-autoencoder size is independent of the system size. We include symmetry regularization terms in the loss function of the GE-autoencoder so that the learned order parameter is also equivariant to the remaining symmetries of the system. By examining the group representation by which the learned order parameter transforms, we are then able to extract information about the associated spontaneous symmetry breaking. We test the GE-autoencoder on the 2D classical ferromagnetic and antiferromagnetic Ising models, finding that the GE-autoencoder (1) accurately determines which symmetries have spontaneously broken at each temperature; (2) estimates the critical temperature in the thermodynamic limit with greater accuracy, robustness, and time-efficiency than a symmetry-agnostic baseline autoencoder; and (3) detects the presence of an external symmetry-breaking magnetic field with greater sensitivity than the baseline method. Finally, we describe various key implementation details, including a new method for extracting the critical temperature estimate from trained autoencoders and calculations of the DNN initialization and learning rate settings required for fair model comparisons. △ Less

Submitted 31 January, 2023; v1 submitted 13 February, 2022; originally announced February 2022.

Comments: 28 pages, 14 figures. For associated code repository see https://github.com/dagrawa2/ssb_detection_ising

Journal ref: Phys. Rev. E 107, 054104 (2023)

A Parameter-Free Differential Evolution Algorithm for the Analytic Continuation of Imaginary Time Correlation Functions

Authors: Nathan S. Nichols, Paul Sokol, Adrian Del Maestro

Abstract: We report on Differential Evolution for Analytic Continuation (DEAC): a parameter-free evolutionary algorithm to generate the dynamic structure factor from imaginary time correlation functions. Our approach to this long-standing problem in quantum many-body physics achieves enhanced spectral fidelity while using fewer compute (CPU) hours. The need for fine-tuning of algorithmic control parameters… ▽ More We report on Differential Evolution for Analytic Continuation (DEAC): a parameter-free evolutionary algorithm to generate the dynamic structure factor from imaginary time correlation functions. Our approach to this long-standing problem in quantum many-body physics achieves enhanced spectral fidelity while using fewer compute (CPU) hours. The need for fine-tuning of algorithmic control parameters is eliminated by embedding them within the genome to be optimized for this evolutionary computation based algorithm. Benchmarks are presented for models where the dynamic structure factor is known exactly, and experimentally relevant results are included for quantum Monte Carlo simulations of bulk $^4$He below the superfluid transition temperature. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 10 pages, 6 figures. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-DEAC

Measuring post-quench entanglement entropy through density correlations

Authors: Adrian Del Maestro, Hatem Barghathi, Bernd Rosenow

Abstract: Following a sudden change of interactions in an integrable system of one-dimensional fermions, we analyze the dependence of the static structure factor on the observation time after the quantum quench. At small waiting times after the quench, we map the system to non-interacting bosons such that we are able to extract their occupation numbers from the Fourier transform of the density-density corre… ▽ More Following a sudden change of interactions in an integrable system of one-dimensional fermions, we analyze the dependence of the static structure factor on the observation time after the quantum quench. At small waiting times after the quench, we map the system to non-interacting bosons such that we are able to extract their occupation numbers from the Fourier transform of the density-density correlation function, and use these to compute a bosonic entropy from a diagonal ensemble. By comparing this bosonic entropy with the asymptotic steady state entanglement entropy per fermion computed with exact diagonalization we find excellent agreement. △ Less

Submitted 12 November, 2021; originally announced November 2021.

Comments: 7 pages, 4 figures. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-QuenchStructureFactorEntanglement

Journal ref: Phys. Rev. Research 4, L022023 (2022)

Balls and Walls: A Compact Unary Coding for Bosonic States

Authors: Hatem Barghathi, Caleb Usadi, Micah Beck, Adrian Del Maestro

Abstract: We introduce a unary coding of bosonic occupation states based on the famous "balls and walls" counting for the number of configurations of $N$ indistinguishable particles on $L$ distinguishable sites. Each state is represented by an integer with a human readable bit string that has a compositional structure allowing for the efficient application of operators that locally modify the number of boso… ▽ More We introduce a unary coding of bosonic occupation states based on the famous "balls and walls" counting for the number of configurations of $N$ indistinguishable particles on $L$ distinguishable sites. Each state is represented by an integer with a human readable bit string that has a compositional structure allowing for the efficient application of operators that locally modify the number of bosons. By exploiting translational and inversion symmetries, we identify a speedup factor of order $L$ over current methods when generating the basis states of bosonic lattice models. The unary coding is applied to a one-dimensional Bose-Hubbard Hamiltonian with up to $L=N=20$, and the time needed to generate the ground state block is reduced to a fraction of the diagonalization time. For the ground state symmetry resolved entanglement, we demonstrate that variational approaches restricting the local bosonic Hilbert space could result in an error that scales with system size. △ Less

Submitted 15 September, 2021; originally announced September 2021.

Comments: Main text: 7 pages, 3 figures + Supplementary material: 3 pages 1 figure. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-UnaryBosonicBasis

Journal ref: Phys. Rev. B 105, L121116 (2022)

A Perspective on Collective Properties of Atoms on 2D Materials

Authors: Adrian Del Maestro, Carlos Wexler, Juan M. Vanegas, Taras Lakoba, Valeri N. Kotov

Abstract: Atoms deposited on two-dimensional (2D) electronic materials, such as graphene, can exhibit unconventional many-body correlations, not accessible in other settings. All of these are driven by van der Waals forces: between the atoms themselves and atom-material interactions. For example $^4$He atoms on 2D materials can potentially form a variety of exotic quantum states of matter, such as two-dimen… ▽ More Atoms deposited on two-dimensional (2D) electronic materials, such as graphene, can exhibit unconventional many-body correlations, not accessible in other settings. All of these are driven by van der Waals forces: between the atoms themselves and atom-material interactions. For example $^4$He atoms on 2D materials can potentially form a variety of exotic quantum states of matter, such as two-dimensional supersolids and superfluids, in addition to solid phases. For the "most quantum" case of a single helium layer we discuss, from a theoretical perspective, how the effective low-energy (Bose-Hubbard) description can take advantage of the extreme sensitivity of this unique system to the interplay between the atomic (helium) and solid-state (graphene) components. Due to the extraordinary variety and tunability of 2D electronic materials, we envisage that a wide range of correlated atomic phases can be realized under favorable conditions. We also outline exciting possibilities in the opposite extreme of many atomic layers forming a liquid on top of graphene -- in this case a so-called "spinodal de-wetting" pattern can form at the liquid-vapor interface which reflects the presence and electronic properties of graphene underneath. Such patterns could be manipulated by choosing different atoms and materials, with potential technological applications. △ Less

Submitted 18 August, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: 10 pages, 9 figures

Journal ref: Adv. Electron. Mater. 8, 2100607 (2022)

Robustness of Helical Edge States Under Edge Reconstruction

Authors: Niels John, Adrian Del Maestro, Bernd Rosenow

Abstract: The helical edge states of time-reversal invariant two-dimensional topological insulators are protected against backscattering in idealized models. In more realistic scenarios with a shallow confining potential at the sample boundary, additional strongly interacting edge states may arise, that could interfere with the topological protection of edge conduction. We find that interaction effects with… ▽ More The helical edge states of time-reversal invariant two-dimensional topological insulators are protected against backscattering in idealized models. In more realistic scenarios with a shallow confining potential at the sample boundary, additional strongly interacting edge states may arise, that could interfere with the topological protection of edge conduction. We find that interaction effects within the reconstructed edges are well described by the Luttinger liquid model. While interactions between this Luttinger liquid and the helical edge states can in principle give rise to dynamical spin polarization and the breaking of time-reversal symmetry, we demonstrate that random spin-orbit coupling strongly suppresses such dynamical spin polarization, resulting in the persistence of near quantized edge conduction. △ Less

Submitted 30 May, 2022; v1 submitted 31 May, 2021; originally announced May 2021.

Comments: 12 pages, 3 figures; updated figures and extended supplemental material

Journal ref: Europhys Lett. 140 26002 (2022)

Two-Dimensional Bose-Hubbard Model for Helium on Graphene

Authors: Jiangyong Yu, Ethan Lauricella, Mohamed Elsayed, Kenneth Shepherd Jr., Nathan S. Nichols, Todd Lombardi, Sang Wook Kim, Carlos Wexler, Juan M. Vanegas, Taras Lakoba, Valeri N. Kotov, Adrian Del Maestro

Abstract: An exciting development in the field of correlated systems is the possibility of realizing two-dimensional (2D) phases of quantum matter. For a systems of bosons, an example of strong correlations manifesting themselves in a 2D environment is provided by helium adsorbed on graphene. We construct the effective Bose-Hubbard model for this system which involves hard-core bosons $(U\approx\infty)$, re… ▽ More An exciting development in the field of correlated systems is the possibility of realizing two-dimensional (2D) phases of quantum matter. For a systems of bosons, an example of strong correlations manifesting themselves in a 2D environment is provided by helium adsorbed on graphene. We construct the effective Bose-Hubbard model for this system which involves hard-core bosons $(U\approx\infty)$, repulsive nearest-neighbor $(V>0)$ and small attractive $(V'<0)$ next-nearest neighbor interactions. The mapping onto the Bose-Hubbard model is accomplished by a variety of many-body techniques which take into account the strong He-He correlations on the scale of the graphene lattice spacing. Unlike the case of dilute ultracold atoms where interactions are effectively point-like, the detailed microscopic form of the short range electrostatic and long range dispersion interactions in the helium-graphene system are crucial for the emergent Bose-Hubbard description. The result places the ground state of the first layer of $^4$He adsorbed on graphene deep in the commensurate solid phase with $1/3$ of the sites on the dual triangular lattice occupied. Because the parameters of the effective Bose-Hubbard model are very sensitive to the exact lattice structure, this opens up an avenue to tune quantum phase transitions in this solid-state system. △ Less

Submitted 4 June, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: 22 pages, 18 figures. Reorganized text and extended an appendix. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-BoseHubbardModelHeAdsorptionGraphene

Journal ref: Phys. Rev. B 103, 235414 (2021)

Theory of Non-Interacting Fermions and Bosons in the Canonical Ensemble

Authors: Hatem Barghathi, Jiangyong Yu, Adrian Del Maestro

Abstract: We present a self-contained theory for the exact calculation of particle number counting statistics of non-interacting indistinguishable particles in the canonical ensemble. This general framework introduces the concept of auxiliary partition functions, and represents a unification of previous distinct approaches with many known results appearing as direct consequences of the developed mathematica… ▽ More We present a self-contained theory for the exact calculation of particle number counting statistics of non-interacting indistinguishable particles in the canonical ensemble. This general framework introduces the concept of auxiliary partition functions, and represents a unification of previous distinct approaches with many known results appearing as direct consequences of the developed mathematical structure. In addition, we introduce a general decomposition of the correlations between occupation numbers in terms of the occupation numbers of individual energy levels, that is valid for both non-degenerate and degenerate spectra. To demonstrate the applicability of the theory in the presence of degeneracy, we compute energy level correlations up to fourth order in a bosonic ring in the presence of a magnetic field. △ Less

Submitted 8 August, 2020; v1 submitted 30 July, 2020; originally announced July 2020.

Comments: 17 pages, 4 figures. Cross-references added. For associated code and data repository see: https://github.com/DelMaestroGroup/CanonicalEnsembleTheory-paper-code

Journal ref: Phys. Rev. Research 2, 043206 (2020)

Dimensional Reduction of Helium-4 Inside Argon Plated MCM-41 Nanopores

Authors: Nathan Nichols, Timothy R. Prisk, Garfield Warren, Paul Sokol, Adrian Del Maestro

Abstract: The angstrom-scale coherence length describing the superfluid wavefunction of helium-4 at low temperatures has prevented its preparation in a truly one-dimensional geometry. Mesoporous ordered silica-based structures, such as the molecular sieve MCM-41, offer a promising avenue towards physical confinement, but the minimal pore diameters that can be chemically synthesized have proven to be too lar… ▽ More The angstrom-scale coherence length describing the superfluid wavefunction of helium-4 at low temperatures has prevented its preparation in a truly one-dimensional geometry. Mesoporous ordered silica-based structures, such as the molecular sieve MCM-41, offer a promising avenue towards physical confinement, but the minimal pore diameters that can be chemically synthesized have proven to be too large to reach the quasi-one-dimensional limit. We present an active nano-engineering approach to this problem by pre-plating MCM-41 with a single, well controlled layer of Ar gas before filling the pores with helium. The structure inside the pore is investigated via experimental adsorption isotherms and neutron scattering measurements that are in agreement with large scale quantum Monte Carlo simulations. The results demonstrate angstrom and kelvin scale tunability of the effective confinement potential experienced by 4He atoms inside the MCM-41, with the Ar layer reducing the diameter of the confining media into a regime where a number of solid layers surround a one-dimensional quantum liquid. △ Less

Submitted 1 October, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 14 pages, 14 figures; updated Fig. 8 with additional data and corrected typos. For associated data and code repository see: https://github.com/DelMaestroGroup/PlatedHe4Nanopores

Journal ref: Phys. Rev. B 102, 144505 (2020)

Equivalence of Spatial and Particle Entanglement Growth After a Quantum Quench

Authors: Adrian Del Maestro, Hatem Barghathi, Bernd Rosenow

Abstract: We analyze fermions after an interaction quantum quench in one spatial dimension and study the growth of the steady state entanglement entropy density under either a spatial mode or particle bipartition. For integrable lattice models, we find excellent agreement between the increase of spatial and particle entanglement entropy, and for chaotic models, an examination of two further neighbor interac… ▽ More We analyze fermions after an interaction quantum quench in one spatial dimension and study the growth of the steady state entanglement entropy density under either a spatial mode or particle bipartition. For integrable lattice models, we find excellent agreement between the increase of spatial and particle entanglement entropy, and for chaotic models, an examination of two further neighbor interaction strengths suggests similar correspondence. This result highlights the generality of the dynamical conversion of entanglement to thermodynamic entropy under time evolution that underlies our current framework of quantum statistical mechanics. △ Less

Submitted 21 July, 2021; v1 submitted 20 December, 2019; originally announced December 2019.

Comments: 11 pages. Now includes analysis of next-nearest-neighbor interactions. Updated code and data repository see, https://github.com/DelMaestroGroup/papers-code-EntanglementQuantumQuench

Journal ref: Phys. Rev. B. 104, 195101 (2021)

Operationally accessible entanglement of one dimensional spinless fermions

Authors: Hatem Barghathi, Emanuel Casiano-Diaz, Adrian Del Maestro

Abstract: For indistinguishable itinerant particles subject to a superselection rule fixing their total number, a portion of the entanglement entropy under a spatial bipartition of the ground state is due to particle fluctuations between subsystems and thus is inaccessible as a resource for quantum information processing. We quantify the remaining operationally accessible entanglement in a model of interact… ▽ More For indistinguishable itinerant particles subject to a superselection rule fixing their total number, a portion of the entanglement entropy under a spatial bipartition of the ground state is due to particle fluctuations between subsystems and thus is inaccessible as a resource for quantum information processing. We quantify the remaining operationally accessible entanglement in a model of interacting spinless fermions on a one dimensional lattice via exact diagonalization and the density matrix renormalization group. We find that the accessible entanglement exactly vanishes at the first order phase transition between a Tomonaga-Luttinger liquid and phase separated solid for attractive interactions and is maximal at the transition to the charge density wave for repulsive interactions. Throughout the phase diagram, we discuss the connection between the accessible entanglement entropy and the variance of the probability distribution describing intra-subregion particle number fluctuations. △ Less

Submitted 29 May, 2019; v1 submitted 8 May, 2019; originally announced May 2019.

Comments: 16 pages, 8 figures, includes new references

Journal ref: Phys. Rev. A 100, 022324 (2019)

Rényi generalization of the operational entanglement entropy

Authors: Hatem Barghathi, C. M. Herdman, Adrian Del Maestro

Abstract: Operationally accessible entanglement in bipartite systems of indistinguishable particles could be reduced due to restrictions on the allowed local operations as a result of particle number conservation. In order to quantify this effect, Wiseman and Vaccaro [Phys. Rev. Lett. 91, 097902 (2003)] introduced an operational measure of the von Neumann entanglement entropy. Motivated by advances in measu… ▽ More Operationally accessible entanglement in bipartite systems of indistinguishable particles could be reduced due to restrictions on the allowed local operations as a result of particle number conservation. In order to quantify this effect, Wiseman and Vaccaro [Phys. Rev. Lett. 91, 097902 (2003)] introduced an operational measure of the von Neumann entanglement entropy. Motivated by advances in measuring Rényi entropies in quantum many-body systems subject to conservation laws, we derive a generalization of the operational entanglement that is both computationally and experimentally accessible. Using the Widom theorem, we investigate its scaling with the size of a spatial subregion for free fermions and find a logarithmically violated area law scaling, similar to the spatial entanglement entropy, with at most, a double-log leading-order correction. A modification of the correlation matrix method confirms our findings in systems of up to $10^5$ particles. △ Less

Submitted 13 August, 2018; v1 submitted 3 April, 2018; originally announced April 2018.

Comments: Main text: 6 pages, 3 figures + Supplementary material: 6 pages

Journal ref: Phys. Rev. Lett. 121, 150501 (2018)

Theory of liquid film growth and wetting instabilities on graphene

Authors: Sanghita Sengupta, Nathan S. Nichols, Adrian Del Maestro, Valeri N. Kotov

Abstract: We investigate wetting phenomena near graphene within the Dzyaloshinskii-Lifshitz-Pitaevskii theory for light gases composed of hydrogen, helium and nitrogen in three different geometries where graphene is either affixed to an insulating substrate, submerged or suspended. We find that the presence of graphene has a significant effect in all configurations. In a suspended geometry where graphene is… ▽ More We investigate wetting phenomena near graphene within the Dzyaloshinskii-Lifshitz-Pitaevskii theory for light gases composed of hydrogen, helium and nitrogen in three different geometries where graphene is either affixed to an insulating substrate, submerged or suspended. We find that the presence of graphene has a significant effect in all configurations. In a suspended geometry where graphene is able to wet on only one side, liquid film growth becomes arrested at a critical thickness which may trigger surface instabilities and pattern formation analogous to spinodal dewetting. These phenomena are also universally present in other two-dimensional materials. △ Less

Submitted 3 May, 2018; v1 submitted 27 November, 2017; originally announced November 2017.

Comments: 13 pages, includes expanded supplemental material with doping dependence

Journal ref: Phys. Rev. Lett. 120, 236802 (2018)

Pair-breaking quantum phase transition in superconducting nanowires

Authors: Hyunjeong Kim, Fédéric Gay, Adrian Del Maestro, Benjamin Sacépé, Andrey Rogachev

Abstract: A quantum phase transition (QPT) between distinct ground states of matter is a wide-spread phenomenon in nature, yet there are only a few experimentally accessible systems where the microscopic mechanism of the transition can be tested and understood. These cases are unique and form the experimentally established foundation for our understanding of quantum critical phenomena. Here we report the di… ▽ More A quantum phase transition (QPT) between distinct ground states of matter is a wide-spread phenomenon in nature, yet there are only a few experimentally accessible systems where the microscopic mechanism of the transition can be tested and understood. These cases are unique and form the experimentally established foundation for our understanding of quantum critical phenomena. Here we report the discovery that a magnetic-field-driven QPT in superconducting nanowires - a prototypical 1d-system - can be fully explained by the critical theory of pair-breaking transitions characterized by a correlation length exponent $ν\approx 1$ and dynamic critical exponent $z \approx 2$. We find that in the quantum critical regime, the electrical conductivity is in agreement with a theoretically predicted scaling function and, moreover, that the theory quantitatively describes the dependence of conductivity on the critical temperature, field magnitude and orientation, nanowire cross sectional area, and microscopic parameters of the nanowire material. At the critical field, the conductivity follows a $T^{(d-2)/z}$ dependence predicted by phenomenological scaling theories and more recently obtained within a holographic framework. Our work uncovers the microscopic processes governing the transition: The pair-breaking effect of the magnetic field on interacting Cooper pairs overdamped by their coupling to electronic degrees of freedom. It also reveals the universal character of continuous quantum phase transitions. △ Less

Submitted 1 March, 2018; v1 submitted 5 September, 2017; originally announced September 2017.

Comments: 22 pages, 5 figures

Journal ref: Nat. Phys. 14, 912 (2018)

Dissipation in mesoscale superfluids

Authors: Adrian Del Maestro, Bernd Rosenow

Abstract: We investigate the maximum speed at which a driven superfluid can flow through a narrow constriction with a size on the order of the healing length. Considering dissipation via the thermal nucleation of quantized vortices, we calculate the critical velocity for superfluid $^4$He and ultracold atomtronic circuits, identify fundamental length and velocity scales, and are thus able to present results… ▽ More We investigate the maximum speed at which a driven superfluid can flow through a narrow constriction with a size on the order of the healing length. Considering dissipation via the thermal nucleation of quantized vortices, we calculate the critical velocity for superfluid $^4$He and ultracold atomtronic circuits, identify fundamental length and velocity scales, and are thus able to present results obtained in widely different temperature and density ranges in a universal framework. For ultra-narrow channels we predict a drastic reduction in the critical velocity as the energy barrier for flow reducing thermally activated phase slip fluctuations is suppressed. △ Less

Submitted 6 April, 2017; originally announced April 2017.

Comments: 13 pages, 7 figures

Journal ref: Phys. Rev. B 95, 140507(R) (2017)

Particle partition entanglement of one dimensional spinless fermions

Authors: Hatem Barghathi, Emanuel Casiano-Diaz, Adrian Del Maestro

Abstract: We investigate the scaling of the Rényi entanglement entropies for a particle bipartition of interacting spinless fermions in one spatial dimension. In the Tomonaga-Luttinger liquid regime, we calculate the second Rényi entanglement entropy and show that the leading order finite-size scaling is equal to a universal logarithm of the system size plus a non-universal constant. Higher-order correction… ▽ More We investigate the scaling of the Rényi entanglement entropies for a particle bipartition of interacting spinless fermions in one spatial dimension. In the Tomonaga-Luttinger liquid regime, we calculate the second Rényi entanglement entropy and show that the leading order finite-size scaling is equal to a universal logarithm of the system size plus a non-universal constant. Higher-order corrections decay as power-laws in the system size with exponents that depend only on the Luttinger parameter. We confirm the universality of our results by investigating the one dimensional $t-V$ model of interacting spinless fermions via exact-diagonalization techniques. The resulting sensitivity of the particle partition entanglement to boundary conditions and statistics supports its utility as a probe of quantum liquids. △ Less

Submitted 25 August, 2017; v1 submitted 30 March, 2017; originally announced March 2017.

Comments: 19 pages, 6 figures

Journal ref: J. Stat. Mech. (2017) 083108

Absence of superfluidity in a quasi one-dimensional parahydrogen fluid adsorbed inside carbon nanotubes

Authors: Adrian Del Maestro, Massimo Boninsegni

Abstract: A recent claim of superfluid behaviour of parahydrogen adsorbed inside armchair carbon nanotubes [M. Rossi and F. Ancilotto, Phys. Rev. B 94, 100502 (2016)] is disproven by means of first principle computer simulations. Conclusive numerical evidence shows that the low temperature equilibrium thermodynamic phase of parahydrogen adsorbed inside a (10,10) armchair nanotube features a crystalline shel… ▽ More A recent claim of superfluid behaviour of parahydrogen adsorbed inside armchair carbon nanotubes [M. Rossi and F. Ancilotto, Phys. Rev. B 94, 100502 (2016)] is disproven by means of first principle computer simulations. Conclusive numerical evidence shows that the low temperature equilibrium thermodynamic phase of parahydrogen adsorbed inside a (10,10) armchair nanotube features a crystalline shell adsorbed on the inner surface of the tube, and a well separated quasi one-dimensional central column showing no evidence of possible fluid or superfluid behaviour. Rather, the system is quasi-crystalline and non-superfluid; its physical character is qualitatively identical to that observed in similar confined settings extensively studied in the past. The incorrect prediction of superfluidity stems from the erroneous use of an effective one-dimensional interaction potential that does not describe parahydrogen at the center of the nanotube. △ Less

Submitted 27 February, 2017; v1 submitted 11 November, 2016; originally announced November 2016.

Comments: Replaced with published version

Journal ref: Phys. Rev. B 95, 054517 (2017)

Entanglement area law in superfluid $^4$He

Authors: C. M. Herdman, P. -N. Roy, R. G. Melko, A. Del Maestro

Abstract: Area laws were first discovered by Bekenstein and Hawking, who found that the entropy of a black hole grows proportional to its surface area, and not its volume. Entropy area laws have since become a fundamental part of modern physics, from the holographic principle in quantum gravity to ground state wavefunctions of quantum matter, where entanglement entropy is generically found to obey area law… ▽ More Area laws were first discovered by Bekenstein and Hawking, who found that the entropy of a black hole grows proportional to its surface area, and not its volume. Entropy area laws have since become a fundamental part of modern physics, from the holographic principle in quantum gravity to ground state wavefunctions of quantum matter, where entanglement entropy is generically found to obey area law scaling. As no experiments are currently capable of directly probing the entanglement area law in naturally occurring many-body systems, evidence of its existence is based on studies of simplified theories. Using new exact microscopic numerical simulations of superfluid $^4$He, we demonstrate for the first time an area law scaling of entanglement entropy in a real quantum liquid in three dimensions. We validate the fundamental principles underlying its physical origin, and present an "entanglement equation of state" showing how it depends on the density of the superfluid. △ Less

Submitted 26 October, 2016; originally announced October 2016.

Comments: main text: 4 pages, 4 figures + supplementary information: 3 pages, 2 figures

Journal ref: Nature Phys. 13, 556 (2017)

Spatial entanglement entropy in the ground state of the Lieb-Liniger model

Authors: C. M. Herdman, P. -N. Roy, R. G. Melko, A. Del Maestro

Abstract: We consider the entanglement between two spatial subregions in the Lieb-Liniger model of bosons in one spatial dimension interacting via a contact interaction. Using ground state path integral quantum Monte Carlo we numerically compute the Rényi entropy of the reduced density matrix of the subsystem as a measure of entanglement. Our numerical algorithm is based on a replica method previously intro… ▽ More We consider the entanglement between two spatial subregions in the Lieb-Liniger model of bosons in one spatial dimension interacting via a contact interaction. Using ground state path integral quantum Monte Carlo we numerically compute the Rényi entropy of the reduced density matrix of the subsystem as a measure of entanglement. Our numerical algorithm is based on a replica method previously introduced by the authors, which we extend to efficiently study the entanglement of spatial subsystems of itinerant bosons. We confirm a logarithmic scaling of the Rényi entropy with subsystem size that is expected from conformal field theory, and compute the non-universal subleading constant for interaction strengths ranging over two orders of magnitude. In the strongly interacting limit, we find agreement with the known free fermion result. △ Less

Submitted 10 June, 2016; originally announced June 2016.

Comments: 14 pages, 11 figures

Journal ref: Phys. Rev. B 94, 064524 (2016)

Adsorption by design: tuning atom-graphene van der Waals interactions via mechanical strain

Authors: Nathan S. Nichols, Adrian Del Maestro, Carlos Wexler, Valeri N. Kotov

Abstract: We aim to understand how the van der Waals force between neutral adatoms and a graphene layer is modified by uniaxial strain and electron correlation effects. A detailed analysis is presented for three atoms (He, H, and Na) and graphene strain ranging from weak to moderately strong. We show that the van der Waals potential can be significantly enhanced by strain, and present applications of our re… ▽ More We aim to understand how the van der Waals force between neutral adatoms and a graphene layer is modified by uniaxial strain and electron correlation effects. A detailed analysis is presented for three atoms (He, H, and Na) and graphene strain ranging from weak to moderately strong. We show that the van der Waals potential can be significantly enhanced by strain, and present applications of our results to the problem of elastic scattering of atoms from graphene. In particular we find that quantum reflection can be significantly suppressed by strain, meaning that dissipative inelastic effects near the surface become of increased importance. Furthermore we introduce a method to independently estimate the Lennard-Jones parameters used in an effective model of He interacting with graphene, and determine how they depend on strain. At short distances, we find that strain tends to reduce the interaction strength by pushing the location of the adsorption potential minima to higher distances above the deformed graphene sheet. This opens up the exciting possibility of mechanically engineering an adsorption potential, with implications for the formation and observation of anisotropic low dimensional superfluid phases. △ Less

Submitted 12 February, 2016; originally announced February 2016.

Comments: 15 pages, 14 figures

Journal ref: Phys. Rev. B 93, 205412 (2016)

Quasi-one-dimensional electron gas bound to a helium-coated nanotube

Authors: Michael Liebrecht, Adrian Del Maestro, Milton W. Cole

Abstract: A much-studied system is the quasi-2D electron gas in image-potential bound states at the surface of helium and hydrogen. In this paper, we report on an analogous quasi-1D system: electrons bound by image-like polarization forces to the surface of a helium-coated carbon nanotube. The potential is computed from an electron-helium pseudopotential, plus a dynamic image term evaluated from a semi-clas… ▽ More A much-studied system is the quasi-2D electron gas in image-potential bound states at the surface of helium and hydrogen. In this paper, we report on an analogous quasi-1D system: electrons bound by image-like polarization forces to the surface of a helium-coated carbon nanotube. The potential is computed from an electron-helium pseudopotential, plus a dynamic image term evaluated from a semi-classical model of the nanotube's response function. Predictions are made for the bound states and potential many-body properties of this novel electron gas for a specific choice of tube radius and film thickness. △ Less

Submitted 29 December, 2015; originally announced December 2015.

Comments: To appear in the Quantum Fluids and Solids 2015 JLTP conference proceedings

Journal ref: J. Low Temp. Phys., 183, 264 (2016)

Entangling qubit registers via many-body states of ultracold atoms

Authors: R. G. Melko, C. M. Herdman, D. Iouchtchenko, P. -N. Roy, A. Del Maestro

Abstract: Inspired by the experimental measurement of the Renyi entanglement entropy in a lattice of ultracold atoms by Islam et al., [Nature 528, 77 (2015)] we propose a method to entangle two spatially-separated qubits using the quantum many-body state as a resource. Through local operations accessible in an experiment, entanglement is transferred to a qubit register from atoms at the ends of a one-dimens… ▽ More Inspired by the experimental measurement of the Renyi entanglement entropy in a lattice of ultracold atoms by Islam et al., [Nature 528, 77 (2015)] we propose a method to entangle two spatially-separated qubits using the quantum many-body state as a resource. Through local operations accessible in an experiment, entanglement is transferred to a qubit register from atoms at the ends of a one-dimensional chain. We compute the operational entanglement, which bounds the entanglement physically transferable from the many-body resource to the register, and discuss a protocol for its experimental measurement. Finally, we explore measures for the amount of entanglement available in the register after transfer, suitable for use in quantum information applications. △ Less

Submitted 30 December, 2015; v1 submitted 20 December, 2015; originally announced December 2015.

Comments: 5 pages, 4 figures; typos fixed

Journal ref: Phys. Rev. A 93, 042336 (2016)

Moments of the inverse participation ratio for the Laplacian on finite regular graphs

Authors: Timothy B. P. Clark, Adrian Del Maestro

Abstract: We investigate the first and second moments of the inverse participation ratio (IPR) for all eigenvectors of the Laplacian on finite random regular graphs with $n$ vertices and degree $z$. By exactly diagonalizing a large set of $z$-regular graphs, we find that as $n$ becomes large, the mean of the inverse participation ratio on each graph, when averaged over a large ensemble of graphs, approaches… ▽ More We investigate the first and second moments of the inverse participation ratio (IPR) for all eigenvectors of the Laplacian on finite random regular graphs with $n$ vertices and degree $z$. By exactly diagonalizing a large set of $z$-regular graphs, we find that as $n$ becomes large, the mean of the inverse participation ratio on each graph, when averaged over a large ensemble of graphs, approaches the numerical value $3$. This universal number is understood as the large-$n$ limit of the average of the quartic polynomial corresponding to the IPR over an appropriate $(n-2)$-dimensional hypersphere of $\mathbb{R}^n$. For a large, but not exhaustive ensemble of graphs, the mean variance of the inverse participation ratio for all graph Laplacian eigenvectors deviates from its continuous hypersphere average due to large graph-to-graph fluctuations that arise from the existence of highly localized modes. △ Less

Submitted 14 November, 2018; v1 submitted 5 June, 2015; originally announced June 2015.

Comments: 24 pages, 10 figures, fixed typos and included new arguments on graph eigenvector distribution

Journal ref: J. Phys. A: Math. Theor. 51, 495003 (2018)

A Simple Model for Long-Range Interacting Pendula

Authors: Owen Myers, Adrian Del Maestro, Junru Wu, Jeffrey S. Marshall

Abstract: We show that the Hamiltonian mean field (HMF) model describes the equilibrium behavior of a system of long pendula with flat bobs that are coupled through long-range interactions (charged or self gravitating). We solve for the canonical partition function in the coordinate frame of the pendula angles. The Hamiltonian in the angles coordinate frame looks similar to the form of the HMF model but wit… ▽ More We show that the Hamiltonian mean field (HMF) model describes the equilibrium behavior of a system of long pendula with flat bobs that are coupled through long-range interactions (charged or self gravitating). We solve for the canonical partition function in the coordinate frame of the pendula angles. The Hamiltonian in the angles coordinate frame looks similar to the form of the HMF model but with the inclusion of an index dependent phase in the interaction term. We also show interesting non-equilibrium behavior of the pendula angles, namely that a quasistationary clustered state can exist when pendula angles are initially ordered by their index. △ Less

Submitted 29 January, 2015; v1 submitted 16 January, 2015; originally announced January 2015.

Journal ref: J. Appl. Phys. 121, 154501 (2017)

Particle partition entanglement of bosonic Luttinger liquids

Authors: C. M. Herdman, A. Del Maestro

Abstract: We consider the Rényi entanglement entropy of bosonic Luttinger liquids under a particle bipartition and demonstrate that the leading order finite-size scaling is logarithmic in the system size with a prefactor equal to the inverse Luttinger parameter. While higher order corrections involve a microscopic length scale, the leading order scaling depends only on this sole dimensionless parameter whic… ▽ More We consider the Rényi entanglement entropy of bosonic Luttinger liquids under a particle bipartition and demonstrate that the leading order finite-size scaling is logarithmic in the system size with a prefactor equal to the inverse Luttinger parameter. While higher order corrections involve a microscopic length scale, the leading order scaling depends only on this sole dimensionless parameter which characterizes the low energy quantum hydrodynamics. This result contrasts the leading entanglement entropy scaling under a spatial bipartition, for which the coefficient is universal and independent of the Luttinger parameter. Using quantum Monte Carlo calculations, we explicitly confirm the scaling predictions of Luttinger liquid theory for the Lieb-Liniger model of $δ$-function interacting bosons in the one dimensional spatial continuum. △ Less

Submitted 25 February, 2015; v1 submitted 19 December, 2014; originally announced December 2014.

Comments: 10 pages, 8 figures; updated figure 5, added references, revised text, reformatted

Journal ref: Phys. Rev. B 91, 184507 (2015)

Critical Flow and Dissipation in a Quasi-One-Dimensional Superfluid

Authors: P-F Duc, M. Savard, M. Petrescu, B. Rosenow, A. Del Maestro, G. Gervais

Abstract: In one of the most celebrated examples of the theory of universal critical phenomena, the phase transition to the superfluid state of $^{4}$He belongs to the same three dimensional $\mathrm{O}(2)$ universality class as the onset of ferromagnetism in a lattice of classical spins with $XY$ symmetry. Below the transition, the superfluid density $ρ_s$ and superfluid velocity $v_s$ increase as power la… ▽ More In one of the most celebrated examples of the theory of universal critical phenomena, the phase transition to the superfluid state of $^{4}$He belongs to the same three dimensional $\mathrm{O}(2)$ universality class as the onset of ferromagnetism in a lattice of classical spins with $XY$ symmetry. Below the transition, the superfluid density $ρ_s$ and superfluid velocity $v_s$ increase as power laws of temperature described by a universal critical exponent constrained to be equal by scale invariance. As the dimensionality is reduced towards one dimension (1D), it is expected that enhanced thermal and quantum fluctuations preclude long-range order, thereby inhibiting superfluidity. We have measured the flow rate of liquid helium and deduced its superfluid velocity in a capillary flow experiment occurring in single $30~$nm long nanopores with radii ranging down from 20~nm to 3~nm. As the pore size is reduced towards the 1D limit, we observe: {\it i)} a suppression of the pressure dependence of the superfluid velocity; {\it ii)} a temperature dependence of $v_{s}$ that surprisingly can be well-fitted by a powerlaw with a single exponent over a broad range of temperatures; and {\it iii)} decreasing critical velocities as a function of radius for channel sizes below $R \simeq 20$~nm, in stark contrast with what is observed in micron sized channels. We interpret these deviations from bulk behaviour as signaling the crossover to a quasi-1D state whereby the size of a critical topological defect is cut off by the channel radius. △ Less

Submitted 16 December, 2014; originally announced December 2014.

Journal ref: Science Advances 1, e1400222 (2015)

A path integral Monte Carlo method for Rényi entanglement entropies

Authors: C. M. Herdman, Stephen Inglis, P. -N. Roy, R. G. Melko, A. Del Maestro

Abstract: We introduce a quantum Monte Carlo algorithm to measure the Rényi entanglement entropies in systems of interacting bosons in the continuum. This approach is based on a path integral ground state method that can be applied to interacting itinerant bosons in any spatial dimension with direct relevance to experimental systems of quantum fluids. We demonstrate how it may be used to compute spatial mod… ▽ More We introduce a quantum Monte Carlo algorithm to measure the Rényi entanglement entropies in systems of interacting bosons in the continuum. This approach is based on a path integral ground state method that can be applied to interacting itinerant bosons in any spatial dimension with direct relevance to experimental systems of quantum fluids. We demonstrate how it may be used to compute spatial mode entanglement, particle partitioned entanglement, and the entanglement of particles, providing insights into quantum correlations generated by fluctuations, indistinguishability and interactions. We present proof-of-principle calculations, and benchmark against an exactly soluble model of interacting bosons in one spatial dimension. As this algorithm retains the fundamental polynomial scaling of quantum Monte Carlo when applied to sign-problem-free models, future applications should allow for the study of entanglement entropy in large scale many-body systems of interacting bosons. △ Less

Submitted 4 August, 2014; v1 submitted 28 April, 2014; originally announced April 2014.

Comments: 19 pages, 13 figures; updated figures, fixed typos, and fixed notational inconsistencies

Journal ref: Phys. Rev. E 90, 013308 (2014)

Quantum Monte Carlo measurement of the chemical potential of helium-4

Authors: C. M. Herdman, A. Rommal, A. Del Maestro

Abstract: A path integral Monte Carlo method based on the worm algorithm has been developed to compute the chemical potential of interacting bosonic quantum fluids. By applying it to finite-sized systems of helium-4 atoms, we have confirmed that the chemical potential scales inversely with the number of particles to lowest order. The introduction of a simple scaling form allows for the extrapolation of the… ▽ More A path integral Monte Carlo method based on the worm algorithm has been developed to compute the chemical potential of interacting bosonic quantum fluids. By applying it to finite-sized systems of helium-4 atoms, we have confirmed that the chemical potential scales inversely with the number of particles to lowest order. The introduction of a simple scaling form allows for the extrapolation of the chemical potential to the thermodynamic limit, where we observe excellent agreement with known experimental results for helium-4 at saturated vapor pressure. We speculate on future applications of the proposed technique, including its use in studies of confined quantum fluids. △ Less

Submitted 20 December, 2013; originally announced December 2013.

Comments: 7 pages, 4 figures

Journal ref: Phys. Rev. B 89, 224502 (2014)

Particle entanglement in continuum many-body systems via quantum Monte Carlo

Authors: C. M. Herdman, P. -N. Roy, R. G. Melko, A. Del Maestro

Abstract: Entanglement of spatial bipartitions, used to explore lattice models in condensed matter physics, may be insufficient to fully describe itinerant quantum many-body systems in the continuum. We introduce a procedure to measure the Rényi entanglement entropies on a particle bipartition, with general applicability to continuum Hamiltonians via path integral Monte Carlo methods. Via direct simulations… ▽ More Entanglement of spatial bipartitions, used to explore lattice models in condensed matter physics, may be insufficient to fully describe itinerant quantum many-body systems in the continuum. We introduce a procedure to measure the Rényi entanglement entropies on a particle bipartition, with general applicability to continuum Hamiltonians via path integral Monte Carlo methods. Via direct simulations of interacting bosons in one spatial dimension, we confirm a logarithmic scaling of the single-particle entanglement entropy with the number of particles in the system. The coefficient of this logarithmic scaling increases with interaction strength, saturating to unity in the strongly interacting limit. Additionally, we show that the single-particle entanglement entropy is bounded by the condensate fraction, suggesting a practical route towards its measurement in future experiments. △ Less

Submitted 7 April, 2014; v1 submitted 30 October, 2013; originally announced October 2013.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. B 89, 140501 (2014)

Local Superfluidity at the Nanoscale

Authors: B. Kulchytskyy, G. Gervais, A. Del Maestro

Abstract: We have performed quantum Monte Carlo simulations measuring the finite size and temperature superfluid response of helium-4 to the linear and rotational motion of the walls of a nanopore. Within the two-fluid model, the portion of the normal liquid dragged along with the boundaries is dependent on the type of motion and the resulting anisotropic superfluid density saturates far below unity at T=0.… ▽ More We have performed quantum Monte Carlo simulations measuring the finite size and temperature superfluid response of helium-4 to the linear and rotational motion of the walls of a nanopore. Within the two-fluid model, the portion of the normal liquid dragged along with the boundaries is dependent on the type of motion and the resulting anisotropic superfluid density saturates far below unity at T=0.5 K. The origin of the saturation is uncovered by computing the spatial distribution of superfluidity, with only the core of the nanopore exhibiting any evidence of phase coherence. The superfluid core displays scaling behavior consistent with Luttinger liquid theory, thereby providing an experimental test for the emergence of a one dimensional quantum liquid. △ Less

Submitted 2 September, 2013; v1 submitted 27 March, 2013; originally announced March 2013.

Comments: Added a figure and extended discussion

Journal ref: Physical Review B 88, 064512 (2013)

Backscattering Between Helical Edge States via Dynamic Nuclear Polarization

Authors: Adrian Del Maestro, Timo Hyart, Bernd Rosenow

Abstract: We show that that the non-equilibrium spin polarization of one dimensional helical edge states at the boundary of a two dimensional topological insulator can dynamically induce a polarization of nuclei via the hyperfine interaction. When combined with a spatially inhomogeneous Rashba coupling, the steady state polarization of the nuclei produces backscattering between the topologically protected e… ▽ More We show that that the non-equilibrium spin polarization of one dimensional helical edge states at the boundary of a two dimensional topological insulator can dynamically induce a polarization of nuclei via the hyperfine interaction. When combined with a spatially inhomogeneous Rashba coupling, the steady state polarization of the nuclei produces backscattering between the topologically protected edge states leading to a reduction in the conductance which persists to zero temperature. We study these effects in both short and long edges, uncovering deviations from Ohmic transport at finite temperature and a current noise spectrum which may hold the fingerprints for experimental verification of the backscattering mechanism. △ Less

Submitted 31 July, 2012; originally announced July 2012.

Comments: 4+ pages, 4 figures

Journal ref: Phys. Rev. B 87, 165440 (2013)

A Luttinger Liquid Core Inside Helium-4 Filled Nanopores

Authors: Adrian Del Maestro

Abstract: As helium-4 is cooled below 2.17 K it undergoes a phase transition to a fundamentally new state of matter known as a superfluid which supports flow without viscosity. This type of dissipationless transport can be observed by forcing helium to travel through a narrow constriction that the normal liquid could not penetrate. Recent experiments have highlighted the feasibility of fabricating smooth po… ▽ More As helium-4 is cooled below 2.17 K it undergoes a phase transition to a fundamentally new state of matter known as a superfluid which supports flow without viscosity. This type of dissipationless transport can be observed by forcing helium to travel through a narrow constriction that the normal liquid could not penetrate. Recent experiments have highlighted the feasibility of fabricating smooth pores with nanometer radii, that approach the truly one dimensional limit where it is believed that a system of bosons (like helium-4) may have startlingly different behavior than in three dimensions. The one dimensional system is predicted to have a linear hydrodynamic description known as Luttinger liquid theory, where no type of long range order can be sustained. In the limit where the pore radius is small, Luttinger liquid theory would predict that helium inside the channel behaves as a sort of quasi-supersolid with all correlations decaying as power-law functions of distance at zero temperature. We have performed large scale quantum Monte Carlo simulations of helium-4 inside nanopores of varying radii at low temperature with realistic helium-helium and helium-pore interactions. The results indicate that helium inside the nanopore forms concentric cylindrical shells surrounding a core that can be described via Luttinger liquid theory and provides insights into the exciting possibility of the experimental detection of this intriguing low dimensional state of matter. △ Less

Submitted 23 January, 2012; originally announced January 2012.

Comments: 26 pages, 10 figures

Journal ref: Int. J. Mod. Phys. B, 26, 1244002 (2012)

Quantum Criticality for Extended Nodes on a Bethe Lattice in the Large Connectivity Limit

Authors: James M. Murray, Adrian Del Maestro, Zlatko Tesanovic

Abstract: Theoretical description of anisotropic systems, such as layered superconductors and coupled spin chains, is often a challenge due to the different natures of interactions along different directions. As a model of such a system, we present an analytical study of d-dimensional "nodes" arranged as the vertices of a Bethe lattice, where each node has nonzero spatial dimension and is described by an O(… ▽ More Theoretical description of anisotropic systems, such as layered superconductors and coupled spin chains, is often a challenge due to the different natures of interactions along different directions. As a model of such a system, we present an analytical study of d-dimensional "nodes" arranged as the vertices of a Bethe lattice, where each node has nonzero spatial dimension and is described by an O(N) quantum rotor model, and there is hopping between neighboring nodes. In the limit of large connectivity on the Bethe lattice, the hopping can be treated by constructing a self-consistent effective action for a single node. This procedure is akin to dynamical mean field theory, but generalized so that spatial as well as quantum fluctuations are taken into account on each node. The quantum phase transition is studied using this effective action for both infinite and finite N. The importance of the Perron-Frobenius uniform mode on the Bethe lattice is discussed, and its elimination via an "infinite range hopping" term shifts the transition, leading to nontrivial critical behavior. We calculate critical exponents and find that the internode hopping reduces the upper and lower critical dimensions each by one, indicating that--at least for the purposes of quantum criticality--the large number of internode couplings is similar to adding a single extra dimension to the theory describing a single node. △ Less

Submitted 2 April, 2012; v1 submitted 8 November, 2011; originally announced November 2011.

Comments: 8 pages, 3 figures. Published version

Journal ref: Phys. Rev. B 85, 115117 (2012)