Showing 1–50 of 121 results for author: Holzmann, M

First principles simulations of dense hydrogen

Authors: Michael Bonitz, Jan Vorberger, Mandy Bethkenhagen, Maximilian Böhme, David Ceperley, Alexey Filinov, Thomas Gawne, Frank Graziani, Gianluca Gregori, Paul Hamann, Stephanie Hansen, Markus Holzmann, S. X. Hu, Hanno Kählert, Valentin Karasiev, Uwe Kleinschmidt, Linda Kordts, Christopher Makait, Burkhard Militzer, Zhandos Moldabekov, Carlo Pierleoni, Martin Preising, Kushal Ramakrishna, Ronald Redmer, Sebastian Schwalbe , et al. (2 additional authors not shown)

Abstract: Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g. planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extre… ▽ More Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g. planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data. rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange-correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principle PIMC and DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Nonrelativistic Limit of Generalized MIT Bag Models and Spectral Inequalities

Authors: Jussi Behrndt, Dale Frymark, Markus Holzmann, Christian Stelzer-Landauer

Abstract: For a family of self-adjoint Dirac operators $-i c (αあるふぁ\cdot \nabla) + \frac{c^2}{2}$ subject to generalized MIT bag boundary conditions on domains in $\mathbb R^3$ it is shown that the nonrelativistic limit in the norm resolvent sense is the Dirichlet Laplacian. This allows to transfer spectral geometry results for Dirichlet Laplacians to Dirac operators for large $c$. For a family of self-adjoint Dirac operators $-i c (αあるふぁ\cdot \nabla) + \frac{c^2}{2}$ subject to generalized MIT bag boundary conditions on domains in $\mathbb R^3$ it is shown that the nonrelativistic limit in the norm resolvent sense is the Dirichlet Laplacian. This allows to transfer spectral geometry results for Dirichlet Laplacians to Dirac operators for large $c$. △ Less

Submitted 26 June, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

Comments: Revised version

On two-dimensional Dirac operators with $δでるた$-shell interactions supported on unbounded curves with straight ends

Authors: Jussi Behrndt, Pavel Exner, Markus Holzmann, Matěj Tušek

Abstract: In this paper we study the self-adjointness and spectral properties of two-dimensional Dirac operators with electrostatic, Lorentz scalar, and anomalous magnetic $δでるた$-shell interactions with constant weights that are supported on a smooth unbounded curve that is straight outside a compact set and whose ends are rays that are not parallel to each other. For all possible combinations of interaction s… ▽ More In this paper we study the self-adjointness and spectral properties of two-dimensional Dirac operators with electrostatic, Lorentz scalar, and anomalous magnetic $δでるた$-shell interactions with constant weights that are supported on a smooth unbounded curve that is straight outside a compact set and whose ends are rays that are not parallel to each other. For all possible combinations of interaction strengths we describe the self-adjoint realizations and compute their essential spectra. Moreover, we prove in different situations the existence of geometrically induced discrete eigenvalues. △ Less

Submitted 30 November, 2023; originally announced December 2023.

Comments: 35 pages

Electronic excitation spectra of molecular hydrogen in Phase I from Quantum Monte Carlo and Many-Body perturbation methods

Authors: Vitaly Gorelov, Markus Holzmann, David M. Ceperley, Carlo Pierleoni

Abstract: We study the electronic excitation spectra in solid molecular hydrogen (phase I) at ambient temperature and 5-90 GPa pressures using Quantum Monte Carlo methods and Many-Body Perturbation Theory. In this range, the system changes from a wide gap molecular insulator to a semiconductor, altering the nature of the excitations from localized to delocalized. Computed gaps and spectra agree with experim… ▽ More We study the electronic excitation spectra in solid molecular hydrogen (phase I) at ambient temperature and 5-90 GPa pressures using Quantum Monte Carlo methods and Many-Body Perturbation Theory. In this range, the system changes from a wide gap molecular insulator to a semiconductor, altering the nature of the excitations from localized to delocalized. Computed gaps and spectra agree with experiments, proving the ability to predict accurately band gaps of many-body systems in presence of nuclear quantum and thermal effects. △ Less

Submitted 20 May, 2024; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: Download the archive to see the Supplemental Material

Training models using forces computed by stochastic electronic structure methods

Authors: David M. Ceperley, Scott Jensen, Yubo Yang, Hongwei Niu, Carlo Pierleoni, Markus Holzmann

Abstract: Quantum Monte Carlo (QMC) can play a very important role in generating accurate data needed for constructing potential energy surfaces. We argue that QMC has advantages in terms of a smaller systematic bias and an ability to cover phase space more completely. The stochastic noise can ease the training of the machine learning model. We discuss how stochastic errors affect the generation of effectiv… ▽ More Quantum Monte Carlo (QMC) can play a very important role in generating accurate data needed for constructing potential energy surfaces. We argue that QMC has advantages in terms of a smaller systematic bias and an ability to cover phase space more completely. The stochastic noise can ease the training of the machine learning model. We discuss how stochastic errors affect the generation of effective models by analyzing the errors within a linear least squares procedure, finding that there is an advantage to having many relatively imprecise data points for constructing models. We then analyze the effect of noise on a model of many-body silicon finding that noise in some situations improves the resulting model. We then study the effect of QMC noise on two machine learning models of dense hydrogen used in a recent study of its phase diagram. The noise enable us to estimate the errors in the model. We conclude with a discussion of future research problems. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Journal ref: Electron. Struct. 6 (2024) 015011

Approximation of Dirac operators with $\boldsymbolδでるた$-shell potentials in the norm resolvent sense

Authors: Jussi Behrndt, Markus Holzmann, Christian Stelzer

Abstract: In this paper the approximation of Dirac operators with general $δでるた$-shell potentials supported on $C^2$-curves in $\mathbb{R}^2$ or $C^2$-surfaces in $\mathbb{R}^3$, which may be bounded or unbounded, is studied. It is shown under suitable conditions on the weight of the $δでるた$-interaction that a family of Dirac operators with regular, squeezed potentials converges in the norm resolvent sense to the… ▽ More In this paper the approximation of Dirac operators with general $δでるた$-shell potentials supported on $C^2$-curves in $\mathbb{R}^2$ or $C^2$-surfaces in $\mathbb{R}^3$, which may be bounded or unbounded, is studied. It is shown under suitable conditions on the weight of the $δでるた$-interaction that a family of Dirac operators with regular, squeezed potentials converges in the norm resolvent sense to the Dirac operator with the $δでるた$-shell interaction. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: 58 pages; comments welcome

Spectral analysis of the Dirac operator with a singular interaction on a broken line

Authors: Dale Frymark, Markus Holzmann, Vladimir Lotoreichik

Abstract: We consider the one-parametric family of self-adjoint realizations of the two-dimensional massive Dirac operator with a Lorentz scalar $δでるた$-shell interaction of strength $τたう\in\mathbb{R}\setminus\{-2,0,2\}$ supported on a broken line of opening angle $2ωおめが$ with $ωおめが\in(0,\fracπぱい{2})$. The essential spectrum of any such self-adjoint realization is symmetric with respect to the origin with a gap around ze… ▽ More We consider the one-parametric family of self-adjoint realizations of the two-dimensional massive Dirac operator with a Lorentz scalar $δでるた$-shell interaction of strength $τたう\in\mathbb{R}\setminus\{-2,0,2\}$ supported on a broken line of opening angle $2ωおめが$ with $ωおめが\in(0,\fracπぱい{2})$. The essential spectrum of any such self-adjoint realization is symmetric with respect to the origin with a gap around zero whose size depends on the mass and, for $τたう< 0$, also on the strength of the interaction, but does not depend on $ωおめが$. As the main result, we prove that for any $N\in\mathbb{N}$ and strength $τたう\in(-\infty,0)\setminus\{-2\}$ the discrete spectrum of any such self-adjoint realization has at least $N$ discrete eigenvalues, with multiplicities taken into account, in the gap of the essential spectrum provided that $ωおめが$ is sufficiently small. Moreover, we obtain an explicit estimate on $ωおめが$ sufficient for this property to hold. For $τたう\in(0,\infty)\setminus\{2\}$, the discrete spectrum consists of at most one simple eigenvalue. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: 28 pages, 3 figures

MSC Class: 47F05; 35P15; 35Q40

Static self energy and effective mass of the homogeneous electron gas from Quantum Monte Carlo calculations

Authors: Markus Holzmann, Francesco Calcavecchia, David M. Ceperley, Valerio Olevano

Abstract: We discuss the methodology of quantum Monte Carlo calculations of the effective mass based on the static self energy, $Σしぐま(k,0)$. We then use variational Monte Carlo calculations of $Σしぐま(k,0)$ of the homogeneous electron gas at various densities to obtain results very close to perturbative $G_0 W_0$ calculations for values of the density parameter $1 \le r_s \le 10$. The obtained values for the effect… ▽ More We discuss the methodology of quantum Monte Carlo calculations of the effective mass based on the static self energy, $Σしぐま(k,0)$. We then use variational Monte Carlo calculations of $Σしぐま(k,0)$ of the homogeneous electron gas at various densities to obtain results very close to perturbative $G_0 W_0$ calculations for values of the density parameter $1 \le r_s \le 10$. The obtained values for the effective mass are close to diagrammatic Monte Carlo results and disagree with previous quantum Monte Carlo calculations based on a heuristic mapping of excitation energies to those of an ideal gas. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Neutral band gap of carbon by quantum Monte Carlo methods

Authors: V. Gorelov, Y. Yang, M. Ruggeri, D. M. Ceperley, C. Pierleoni, M. Holzmann

Abstract: We present a method of calculating the energy gap of a charge-neutral excitation using only ground-state calculations. We report Quantum Monte Carlo calculations of $Γがんま\rightarrowΓがんま$ and $Γがんま\rightarrow X$ particle-hole excitation energies in diamond carbon. We analyze the finite-size effect and find the same $1/L$ decay rate as that in a charged excitation, where $L$ is the linear extension of the su… ▽ More We present a method of calculating the energy gap of a charge-neutral excitation using only ground-state calculations. We report Quantum Monte Carlo calculations of $Γがんま\rightarrowΓがんま$ and $Γがんま\rightarrow X$ particle-hole excitation energies in diamond carbon. We analyze the finite-size effect and find the same $1/L$ decay rate as that in a charged excitation, where $L$ is the linear extension of the supercell. This slow decay is attributed to the delocalized nature of the excitation in supercells too small to accommodate excitonic binding effects. At larger system sizes, the apparent $1/L$ decay crosses over to a $1/L^3$ behavior. Estimation of the scale of exciton binding can be used to correct finite-size effects of neutral gaps. △ Less

Submitted 14 August, 2023; v1 submitted 31 March, 2023; originally announced March 2023.

Comments: 11 pages, 2 figures, 1 table

Journal ref: Condensed Matter Physics, 2023, vol. 26, No. 3, 33701

On the single layer boundary integral operator for the Dirac equation

Authors: Markus Holzmann

Abstract: This paper is devoted to the analysis of the single layer boundary integral operator $\mathcal{C}_z$ for the Dirac equation in the two- and three-dimensional situation. The map $\mathcal{C}_z$ is the strongly singular integral operator having the integral kernel of the resolvent of the free Dirac operator $A_0$ and $z$ belongs to the resolvent set of $A_0$. In the case of smooth boundaries fine ma… ▽ More This paper is devoted to the analysis of the single layer boundary integral operator $\mathcal{C}_z$ for the Dirac equation in the two- and three-dimensional situation. The map $\mathcal{C}_z$ is the strongly singular integral operator having the integral kernel of the resolvent of the free Dirac operator $A_0$ and $z$ belongs to the resolvent set of $A_0$. In the case of smooth boundaries fine mapping properties and a decomposition of $\mathcal{C}_z$ in a 'positive' and 'negative' part are analyzed. The obtained results can be applied in the treatment of Dirac operators with singular electrostatic, Lorentz scalar, and anomalous magnetic interactions that are combined in a critical way. △ Less

Submitted 3 October, 2023; v1 submitted 16 December, 2022; originally announced December 2022.

Comments: 23 pages; to appear in Complex Anal. Oper. Theory

Boundary triples and Weyl functions for Dirac operators with singular interactions

Authors: Jussi Behrndt, Markus Holzmann, Christian Stelzer, Georg Stenzel

Abstract: In this article we develop a systematic approach to treat Dirac operators $A_{ηいーた, τたう, λらむだ}$ with singular electrostatic, Lorentz scalar, and anomalous magnetic interactions of strengths $ηいーた, τたう, λらむだ\in \mathbb{R}$, respectively, supported on points in $\mathbb{R}$, curves in $\mathbb{R}^2$, and surfaces in $\mathbb{R}^3$ that is based on boundary triples and their associated Weyl functions. First, we disc… ▽ More In this article we develop a systematic approach to treat Dirac operators $A_{ηいーた, τたう, λらむだ}$ with singular electrostatic, Lorentz scalar, and anomalous magnetic interactions of strengths $ηいーた, τたう, λらむだ\in \mathbb{R}$, respectively, supported on points in $\mathbb{R}$, curves in $\mathbb{R}^2$, and surfaces in $\mathbb{R}^3$ that is based on boundary triples and their associated Weyl functions. First, we discuss the one-dimensional case which also serves as a motivation for the multidimensional setting. Afterwards, in the two and three-dimensional situation we construct quasi, generalized, and ordinary boundary triples and their Weyl functions, and provide a detailed characterization of the associated Sobolev spaces, trace theorems, and the mapping properties of integral operators which play an important role in the analysis of $A_{ηいーた, τたう, λらむだ}$. We make a substantial step towards more rough interaction supports $Σしぐま$ and consider general compact Lipschitz hypersurfaces. We derive conditions for the interaction strengths such that the operators $A_{ηいーた, τたう, λらむだ}$ are self-adjoint, obtain a Krein-type resolvent formula, and characterize the essential and discrete spectrum. These conditions include purely Lorentz scalar and purely non-critical anomalous magnetic interactions as well as the confinement case, the latter having an important application in the mathematical description of graphene. Using a certain ordinary boundary triple, we show the self-adjointness of $A_{ηいーた, τたう, λらむだ}$ for arbitrary combinations of the interaction strengths (including critical ones) under the condition that $Σしぐま$ is $C^{\infty}$-smooth and derive its spectral properties. In particular, in the critical case, a loss of Sobolev regularity in the operator domain and a possible additional point of the essential spectrum are observed. △ Less

Submitted 18 August, 2023; v1 submitted 9 November, 2022; originally announced November 2022.

Comments: 56 pages; to appear in Reviews in Mathematical Physics

Data Subsampling for Bayesian Neural Networks

Authors: Eiji Kawasaki, Markus Holzmann

Abstract: Markov Chain Monte Carlo (MCMC) algorithms do not scale well for large datasets leading to difficulties in Neural Network posterior sampling. In this paper, we apply a generalization of the Metropolis Hastings algorithm that allows us to restrict the evaluation of the likelihood to small mini-batches in a Bayesian inference context. Since it requires the computation of a so-called "noise penalty"… ▽ More Markov Chain Monte Carlo (MCMC) algorithms do not scale well for large datasets leading to difficulties in Neural Network posterior sampling. In this paper, we apply a generalization of the Metropolis Hastings algorithm that allows us to restrict the evaluation of the likelihood to small mini-batches in a Bayesian inference context. Since it requires the computation of a so-called "noise penalty" determined by the variance of the training loss function over the mini-batches, we refer to this data subsampling strategy as Penalty Bayesian Neural Networks - PBNNs. Its implementation on top of MCMC is straightforward, as the variance of the loss function merely reduces the acceptance probability. Comparing to other samplers, we empirically show that PBNN achieves good predictive performance for a given mini-batch size. Varying the size of the mini-batches enables a natural calibration of the predictive distribution and provides an inbuilt protection against overfitting. We expect PBNN to be particularly suited for cases when data sets are distributed across multiple decentralized devices as typical in federated learning. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Stable solid molecular hydrogen above 900K from a machine-learned potential trained with diffusion Quantum Monte Carlo

Authors: Hongwei Niu, Yubo Yang, Scott Jensen, Markus Holzmann, Carlo Pierleoni, David M. Ceperley

Abstract: We survey the phase diagram of high-pressure molecular hydrogen with path integral molecular dynamics using a machine-learned interatomic potential trained with Quantum Monte Carlo forces and energies. Besides the HCP and C2/c-24 phases, we find two new stable phases both with molecular centers in the Fmmm-4 structure, separated by a molecular orientation transition with temperature. The high temp… ▽ More We survey the phase diagram of high-pressure molecular hydrogen with path integral molecular dynamics using a machine-learned interatomic potential trained with Quantum Monte Carlo forces and energies. Besides the HCP and C2/c-24 phases, we find two new stable phases both with molecular centers in the Fmmm-4 structure, separated by a molecular orientation transition with temperature. The high temperature isotropic Fmmm-4 phase has a reentrant melting line with a maximum at higher temperature (1450K at 150GPa) than previously estimated and crosses the liquid-liquid transition line around 1200K and 200GPa. △ Less

Submitted 14 February, 2023; v1 submitted 1 September, 2022; originally announced September 2022.

Comments: 4 pages, 4 figures, updated acknowledgment

Journal ref: Phys. Rev. Lett. 130, 076102 (2023)

Two-dimensional Dirac operators with general $δでるた$-shell interactions supported on a straight line

Authors: Jussi Behrndt, Markus Holzmann, Matěj Tušek

Abstract: In this paper the two-dimensional Dirac operator with a general hermitian $δでるた$-shell interaction supported on a straight line is introduced as a self-adjoint operator and its spectral properties are investigated in detail. In particular, it is demonstrated that the singularly continuous spectrum is always empty and that by switching a certain $δでるた$-shell interaction on, it is possible to generate an… ▽ More In this paper the two-dimensional Dirac operator with a general hermitian $δでるた$-shell interaction supported on a straight line is introduced as a self-adjoint operator and its spectral properties are investigated in detail. In particular, it is demonstrated that the singularly continuous spectrum is always empty and that by switching a certain $δでるた$-shell interaction on, it is possible to generate an eigenvalue in the gap of the spectrum of the free operator or to partially or even fully close the gap. This suggests that the studied operators may serve as interesting continuum toy-models for Dirac materials. Finally, approximations by Dirac operators with regular potentials are presented. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: 3 figures

MSC Class: 81Q10; 47B93

Schrödinger operators with oblique transmission conditions in $\mathbb{R}^2$

Authors: Jussi Behrndt, Markus Holzmann, Georg Stenzel

Abstract: In this paper we study the spectrum of self-adjoint Schrödinger operators in $L^2(\mathbb{R}^2)$ with a new type of transmission conditions along a smooth closed curve $Σしぐま\subseteq \mathbb{R}^2$. Although these $\textit{oblique}$ transmission conditions are formally similar to $δでるた'$-conditions on $Σしぐま$ (instead of the normal derivative here the Wirtinger derivative is used) the spectral properties are… ▽ More In this paper we study the spectrum of self-adjoint Schrödinger operators in $L^2(\mathbb{R}^2)$ with a new type of transmission conditions along a smooth closed curve $Σしぐま\subseteq \mathbb{R}^2$. Although these $\textit{oblique}$ transmission conditions are formally similar to $δでるた'$-conditions on $Σしぐま$ (instead of the normal derivative here the Wirtinger derivative is used) the spectral properties are significantly different: it turns out that for attractive interaction strengths the discrete spectrum is always unbounded from below. Besides this unexpected spectral effect we also identify the essential spectrum, and we prove a Krein-type resolvent formula and a Birman-Schwinger principle. Furthermore, we show that these Schrödinger operators with oblique transmission conditions arise naturally as non-relativistic limits of Dirac operators with electrostatic and Lorentz scalar $δでるた$-interactions justifying their usage as models in quantum mechanics. △ Less

Submitted 27 March, 2023; v1 submitted 5 July, 2022; originally announced July 2022.

Comments: 16 pages; to appear in Comm. Math. Phys

Wave function Ansatz (but Periodic) Networks and the Homogeneous Electron Gas

Authors: Max Wilson, Saverio Moroni, Markus Holzmann, Nicholas Gao, Filip Wudarski, Tejs Vegge, Arghya Bhowmik

Abstract: We design a neural network Ansatz for variationally finding the ground-state wave function of the Homogeneous Electron Gas, a fundamental model in the physics of extended systems of interacting fermions. We study the spin-polarised and paramagnetic phases with 7, 14 and 19 electrons over a broad range of densities from $r_s=1$ to $r_s=100$, obtaining similar or higher accuracy compared to a state-… ▽ More We design a neural network Ansatz for variationally finding the ground-state wave function of the Homogeneous Electron Gas, a fundamental model in the physics of extended systems of interacting fermions. We study the spin-polarised and paramagnetic phases with 7, 14 and 19 electrons over a broad range of densities from $r_s=1$ to $r_s=100$, obtaining similar or higher accuracy compared to a state-of-the-art iterative backflow baseline even in the challenging regime of very strong correlation. Our work extends previous applications of neural network Ansätze to molecular systems with methods for handling periodic boundary conditions, and makes two notable changes to improve performance: splitting the pairwise streams by spin alignment and generating backflow coordinates for the orbitals from the network. We illustrate the advantage of our high quality wave functions in computing the reduced single particle density matrix. This contribution establishes neural network models as flexible and high precision Ansätze for periodic electronic systems, an important step towards applications to crystalline solids. △ Less

Submitted 23 May, 2023; v1 submitted 2 February, 2022; originally announced February 2022.

Journal ref: Phys. Rev. B 107, 235139 (2023)

Extended Bose-Hubbard model with dipolar excitons

Authors: C. Lagoin, U. Bhattacharya, T. Grass, R. Chhajlany, T. Salamon, K. Baldwin, L. Pfeiffer, M. Lewenstein, M. Holzmann, F. Dubin

Abstract: The Hubbard model constitutes one of the most celebrated theoretical frameworks of condensed-matter physics. It describes strongly correlated phases of interacting quantum particles confined in lattice potentials. For bosons, the Hubbard Hamiltonian has been deeply scrutinised for short-range on-site interactions. However, accessing longer-range couplings has remained elusive experimentally. This… ▽ More The Hubbard model constitutes one of the most celebrated theoretical frameworks of condensed-matter physics. It describes strongly correlated phases of interacting quantum particles confined in lattice potentials. For bosons, the Hubbard Hamiltonian has been deeply scrutinised for short-range on-site interactions. However, accessing longer-range couplings has remained elusive experimentally. This marks the frontier towards the extended Bose-Hubbard Hamiltonian that allows insulating ordered phases at fractional lattice fillings. Here we implement this Hamiltonian by confining semiconductor dipolar excitons in an artificial two-dimensional square lattice. Strong dipolar repulsions between nearest neighbouring lattice sites then stabilise an insulating state at half filling. This characteristic feature of the extended Bose-Hubbard model exhibits signatures theoretically expected for a checkerboard spatial order. Our work thus highlights that dipolar excitons enable controlled implementations of boson-like arrays with strong off-site interactions, in lattices with programmable geometries and over 100 sites. △ Less

Submitted 13 September, 2022; v1 submitted 10 January, 2022; originally announced January 2022.

Comments: 14 pages, 7 figures

Dirac operator spectrum in tubes and layers with a zigzag type boundary

Authors: Pavel Exner, Markus Holzmann

Abstract: We derive a number of spectral results for Dirac operators in geometrically nontrivial regions in $\mathbb{R}^2$ and $\mathbb{R}^3$ of tube or layer shapes with a zigzag type boundary using the corresponding properties of the Dirichlet Laplacian. We derive a number of spectral results for Dirac operators in geometrically nontrivial regions in $\mathbb{R}^2$ and $\mathbb{R}^3$ of tube or layer shapes with a zigzag type boundary using the corresponding properties of the Dirichlet Laplacian. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 16 pages, no figures

MSC Class: 81Q05; 35Q40

Efficient computations of continuous action densities of states for lattice models

Authors: Biagio Lucini, Olmo Francesconi, Markus Holzmann, David Lancaster, Antonio Rago

Abstract: The Logarithmic Linear Relaxation (LLR) algorithm is an efficient method for computing densities of states for systems with a continuous spectrum. A key feature of this method is exponential error reduction, which allows us to evaluate the density of states of a system over hundreds of thousands of orders of magnitude with a fixed level of relative accuracy. As a consequence of exponential error r… ▽ More The Logarithmic Linear Relaxation (LLR) algorithm is an efficient method for computing densities of states for systems with a continuous spectrum. A key feature of this method is exponential error reduction, which allows us to evaluate the density of states of a system over hundreds of thousands of orders of magnitude with a fixed level of relative accuracy. As a consequence of exponential error reduction, the LLR method provides a robust alternative to traditional Monte Carlo calculations in cases in which states suppressed by the Boltzmann weight play nevertheless a relevant role, e.g., as transition regions between dominant configuration sets. After reviewing the algorithm, we will show an application in U(1) Lattice Gauge Theory that has enabled us to obtain the most accurate estimate of the critical coupling with modest computational resources, defeating exponential tunneling times between metastable vacua. As a further showcase, we will then present an application of the LLR method to the decorrelation of the topological charge in SU(3) Lattice Gauge Theory near the continuum limit. Finally, we will review in general applications of the LLR algorithm to systems affected by a strong sign problem and discuss the case of the Bose gas at finite chemical potential. △ Less

Submitted 30 October, 2021; originally announced November 2021.

Comments: 6 pages, 3 figures. Talk presented by B. Lucini at the XXXII IUPAP Conference on Computational Physics 21, Coventry, UK, 1-5 August 2021

The fate of Landau levels under $δでるた$-interactions

Authors: Jussi Behrndt, Markus Holzmann, Vladimir Lotoreichik, Georgi Raikov

Abstract: We consider the self-adjoint Landau Hamiltonian $H_0$ in $L^2(\mathbb{R}^2)$ whose spectrum consists of infinitely degenerate eigenvalues $Λらむだ_q$, $q \in \mathbb{Z}_+$, and the perturbed operator $H_\upsilon = H_0 + \upsilonδでるた_Γがんま$, where $Γがんま\subset \mathbb{R}^2$ is a regular Jordan $C^{1,1}$-curve, and $\upsilon \in L^p(Γがんま;\mathbb{R})$, $p>1$, has a constant sign. We investigate… ▽ More We consider the self-adjoint Landau Hamiltonian $H_0$ in $L^2(\mathbb{R}^2)$ whose spectrum consists of infinitely degenerate eigenvalues $Λらむだ_q$, $q \in \mathbb{Z}_+$, and the perturbed operator $H_\upsilon = H_0 + \upsilonδでるた_Γがんま$, where $Γがんま\subset \mathbb{R}^2$ is a regular Jordan $C^{1,1}$-curve, and $\upsilon \in L^p(Γがんま;\mathbb{R})$, $p>1$, has a constant sign. We investigate ${\rm Ker}(H_\upsilon -Λらむだ_q)$, $q \in \mathbb{Z}_+$, and show that generically $$0 \leq {\rm dim \, Ker}(H_\upsilon -Λらむだ_q) - {\rm dim \, Ker}(T_q(\upsilon δでるた_Γがんま)) < \infty,$$ where $T_q(\upsilon δでるた_Γがんま) = p_q (\upsilon δでるた_Γがんま)p_q$, is an operator of Berezin-Toeplitz type, acting in $p_q L^2(\mathbb{R}^2)$, and $p_q$ is the orthogonal projection on ${\rm Ker}\,(H_0 -Λらむだ_q)$. If $\upsilon \neq 0$ and $q = 0$, we prove that ${\rm Ker}\,(T_0(\upsilon δでるた_Γがんま)) = \{0\}$. If $q \geq 1$, and $Γがんま= \mathcal{C}_r$ is a circle of radius $r$, we show that ${\rm dim \, Ker} (T_q(δでるた_{\mathcal{C}_r})) \leq q$, and the set of $r \in (0,\infty)$ for which ${\rm dim \, Ker}(T_q(δでるた_{\mathcal{C}_r})) \geq 1$, is infinite and discrete. △ Less

Submitted 4 October, 2021; v1 submitted 15 September, 2021; originally announced September 2021.

Comments: 28 pages; to appear in Journal of Spectral Theory

Spectral transition for Dirac operators with electrostatic $δでるた$-shell potentials supported on the straight line

Authors: Jussi Behrndt, Markus Holzmann, Matěj Tušek

Abstract: In this note the two dimensional Dirac operator $A_ηいーた$ with an electrostatic $δでるた$-shell interaction of strength $ηいーた\in\mathbb R$ supported on a straight line is studied. We observe a spectral transition in the sense that for the critical interaction strengths $ηいーた=\pm 2$ the continuous spectrum of $A_ηいーた$ inside the spectral gap of the free Dirac operator $A_0$ collapses abruptly to a single point. In this note the two dimensional Dirac operator $A_ηいーた$ with an electrostatic $δでるた$-shell interaction of strength $ηいーた\in\mathbb R$ supported on a straight line is studied. We observe a spectral transition in the sense that for the critical interaction strengths $ηいーた=\pm 2$ the continuous spectrum of $A_ηいーた$ inside the spectral gap of the free Dirac operator $A_0$ collapses abruptly to a single point. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: 11 pages

Time correlation functions for quantum systems: validating Bayesian approaches for harmonic oscillators and beyond

Authors: Vladislav Efremkin, Jean-Louis Barrat, Stefano Mossa, Markus Holzmann

Abstract: The quantum harmonic oscillator is the fundamental building block to compute thermal properties of virtually any dielectric crystal at low temperatures in terms of phonons, extended further to cases with anharmonic couplings, or even disordered solids. In general, Path Integral Monte Carlo (PIMC) or Molecular Dynamics (PIMD) methods are powerful tools to determine stochastically thermodynamic quan… ▽ More The quantum harmonic oscillator is the fundamental building block to compute thermal properties of virtually any dielectric crystal at low temperatures in terms of phonons, extended further to cases with anharmonic couplings, or even disordered solids. In general, Path Integral Monte Carlo (PIMC) or Molecular Dynamics (PIMD) methods are powerful tools to determine stochastically thermodynamic quantities without systematic bias, not relying on perturbative schemes. Addressing transport properties, for instance calculating thermal conductivity from PIMC, however, is substantially more difficult. Although correlation functions of current operators can be determined by PIMC from analytic continuation on the imaginary-time axis, Bayesian methods are usually employed for the numerical inversion back to real-time response functions. This task not only strongly relies on the accuracy of the PIMC data, but also introduces noticeable dependence on the model used for the inversion. Here, we address both difficulties with care. In particular, we first devise improved estimators for current correlations which substantially reduce the variance of the PIMC data. Next, we provide a neat statistical approach to the inversion problem, blending into a fresh workflow the classical stochastic maximum entropy method together with recent notions borrowed from statistical learning theory. We test our ideas on a single harmonic oscillator and a collection of oscillators with a continuous distribution of frequencies, and provide indications of the performance of our method in the case of a particle in a double well potential. This work establishes solid grounds for an unbiased, fully quantum mechanical calculation of transport properties in solids. △ Less

Submitted 7 September, 2021; v1 submitted 20 May, 2021; originally announced May 2021.

Glacier Calving Front Segmentation Using Attention U-Net

Authors: Michael Holzmann, Amirabbas Davari, Thorsten Seehaus, Matthias Braun, Andreas Maier, Vincent Christlein

Abstract: An essential climate variable to determine the tidewater glacier status is the location of the calving front position and the separation of seasonal variability from long-term trends. Previous studies have proposed deep learning-based methods to semi-automatically delineate the calving fronts of tidewater glaciers. They used U-Net to segment the ice and non-ice regions and extracted the calving fr… ▽ More An essential climate variable to determine the tidewater glacier status is the location of the calving front position and the separation of seasonal variability from long-term trends. Previous studies have proposed deep learning-based methods to semi-automatically delineate the calving fronts of tidewater glaciers. They used U-Net to segment the ice and non-ice regions and extracted the calving fronts in a post-processing step. In this work, we show a method to segment the glacier calving fronts from SAR images in an end-to-end fashion using Attention U-Net. The main objective is to investigate the attention mechanism in this application. Adding attention modules to the state-of-the-art U-Net network lets us analyze the learning process by extracting its attention maps. We use these maps as a tool to search for proper hyperparameters and loss functions in order to generate higher qualitative results. Our proposed attention U-Net performs comparably to the standard U-Net while providing additional insight into those regions on which the network learned to focus more. In the best case, the attention U-Net achieves a 1.5% better Dice score compared to the canonical U-Net with a glacier front line prediction certainty of up to 237.12 meters. △ Less

Submitted 8 January, 2021; originally announced January 2021.

A simplified approach to the repulsive Bose gas from low to high densities and its numerical accuracy

Authors: Eric A. Carlen, Markus Holzmann, Ian Jauslin, Elliott H. Lieb

Abstract: In 1963, a Simple Approach was developed to study the ground state energy of an interacting Bose gas. It consists in the derivation of an Equation, which is not based on perturbation theory, and which gives the exact expansion of the energy at low densities. This Equation is expressed directly in the thermodynamic limit, and only involves functions of $3$ variables, rather than $3N$. Here, we revi… ▽ More In 1963, a Simple Approach was developed to study the ground state energy of an interacting Bose gas. It consists in the derivation of an Equation, which is not based on perturbation theory, and which gives the exact expansion of the energy at low densities. This Equation is expressed directly in the thermodynamic limit, and only involves functions of $3$ variables, rather than $3N$. Here, we revisit this approach, and show that the Equation yields accurate predictions for various observables for all densities. Specifically, in addition to the ground state energy, we have shown that the Simple Approach gives predictions for the condensate fraction, two-point correlation function, and momentum distribution. We have carried out a variety of tests by comparing the predictions of the Equation with Quantum Monte Carlo calculations, and have found remarkable agreement. We thus show that the Simple Approach provides a new theoretical tool to understand the behavior of the many-body Bose gas, not only in the small and large density ranges, which have been studied before, but also in the range of intermediate density, for which little is known. △ Less

Submitted 24 March, 2021; v1 submitted 21 November, 2020; originally announced November 2020.

Comments: v4: change title

Journal ref: Phys. Rev. A 103, 053309 (2021)

Electronic structure and optical properties of quantum crystals from first principles calculations in the Born-Oppenheimer approximation

Authors: Vitaly Gorelov, David M. Ceperley, Markus Holzmann, Carlo Pierleoni

Abstract: We develop a formalism to accurately account for the renormalization of electronic structure due to quantum and thermal nuclear motions within the Born-Oppenheimer approximation. We focus on the fundamental energy gap obtained from electronic addition and removal energies from Quantum Monte Carlo calculations in either the canonical or grand canonical ensembles. The formalism applies as well to ef… ▽ More We develop a formalism to accurately account for the renormalization of electronic structure due to quantum and thermal nuclear motions within the Born-Oppenheimer approximation. We focus on the fundamental energy gap obtained from electronic addition and removal energies from Quantum Monte Carlo calculations in either the canonical or grand canonical ensembles. The formalism applies as well to effective single electron theories such as those based on Density Functional Theory. We show that electronic (Bloch) crystal momentum can be restored by marginalizing the total electron-ion wave function with respect to the nuclear equilibrium distribution, and we describe an explicit procedure to establish the band structure of electronic excitations for quantum crystals within the Born-Oppenheimer approximation. Based on the Kubo-Greenwood equation, we discuss the effects of nuclear motion on optical conductivity. Our methodology applies to the low temperature regime where nuclear motion is quantized and in general differs from the semi-classical approximation. We apply our method to study the electronic structure of C2/c-24 crystalline hydrogen at 200K and 250 GPa and discuss the optical absorption profile of hydrogen crystal at 200K and carbon diamond at 297K. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Bath-induced Zeno localization in driven many-body quantum systems

Authors: Thibaud Maimbourg, Denis M. Basko, Markus Holzmann, Alberto Rosso

Abstract: We study a quantum interacting spin system subject to an external drive and coupled to a thermal bath of spatially localized vibrational modes, serving as a model of Dynamic Nuclear Polarization. We show that even when the many-body eigenstates of the system are ergodic, a sufficiently strong coupling to the bath may effectively localize the spins due to many-body quantum Zeno effect, as manifeste… ▽ More We study a quantum interacting spin system subject to an external drive and coupled to a thermal bath of spatially localized vibrational modes, serving as a model of Dynamic Nuclear Polarization. We show that even when the many-body eigenstates of the system are ergodic, a sufficiently strong coupling to the bath may effectively localize the spins due to many-body quantum Zeno effect, as manifested by the hole-burning shape of the electron paramagnetic resonance spectrum. Our results provide an explanation of the breakdown of the thermal mixing regime experimentally observed above 4 - 5 Kelvin. △ Less

Submitted 23 March, 2021; v1 submitted 24 September, 2020; originally announced September 2020.

Comments: Revised version. 4 + 22 pages

Journal ref: Phys. Rev. Lett. 126, 120603 (2021)

Quasi-condensation of bilayer excitons in a periodic potential

Authors: Camille Lagoin, Stephan Suffit, Kenneth West, Kirk Baldwin, Loren Pfeiffer, Markus Holzmann, Francois Dubin

Abstract: We study two-dimensional excitons confined in a lattice potential, for high fillings of the lattice sites. We show that a quasi-condensate is possibly formed for small values of the lattice depth, but for larger ones the critical phase-space density for quasi-condensation rapidly exceeds our experimental reach, due to the increase of the excitons effective mass. On the other hand, in the regime of… ▽ More We study two-dimensional excitons confined in a lattice potential, for high fillings of the lattice sites. We show that a quasi-condensate is possibly formed for small values of the lattice depth, but for larger ones the critical phase-space density for quasi-condensation rapidly exceeds our experimental reach, due to the increase of the excitons effective mass. On the other hand, in the regime of a deep lattice potential where excitons are strongly localised at the lattice sites, we show that an array of phase-independent quasi-condensates, different from a Mott insulating phase, is realised. △ Less

Submitted 7 December, 2020; v1 submitted 16 September, 2020; originally announced September 2020.

Comments: 5 pages 4 figures

Journal ref: Phys. Rev. Lett. 126, 067404 (2021)

Electronic energy gap closure and metal-insulator transition in dense liquid hydrogen

Authors: Vitaly Gorelov, David M. Ceperley, Markus Holzmann, Carlo Pierleoni

Abstract: Using Quantum Monte Carlo (QMC) calculations, we investigate the insulator-metal transition observed in liquid hydrogen at high pressure. Below the critical temperature of the transition from the molecular to the atomic liquid, the fundamental electronic gap closure occurs abruptly, with a small discontinuity reflecting the weak first-order transition in the thermodynamic equation of state. Above… ▽ More Using Quantum Monte Carlo (QMC) calculations, we investigate the insulator-metal transition observed in liquid hydrogen at high pressure. Below the critical temperature of the transition from the molecular to the atomic liquid, the fundamental electronic gap closure occurs abruptly, with a small discontinuity reflecting the weak first-order transition in the thermodynamic equation of state. Above the critical temperature, molecular dissociation sets in while the gap is still open. When the gap closes, the decay of the off-diagonal reduced density matrix shows that the liquid enters a gapless, but localized phase: there is a cross-over between the insulating and the metallic liquids. Compared to different DFT functionals, our QMC calculations provide larger values for the fundamental gap and the electronic density of states close to the band edges, indicating that optical properties from DFT potentially benefit from error cancellations. △ Less

Submitted 1 September, 2020; originally announced September 2020.

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

Quantum Monte Carlo determination of the principal Hugoniot of deuterium

Authors: Michele Ruggeri, Markus Holzmann, David M. Ceperley, Carlo Pierleoni

Abstract: We present Coupled Electron-Ion Monte Carlo results for the principal Hugoniot of deuterium together with an accurate study of the initial reference state of shock wave experiments. We discuss the influence of nuclear quantum effects, thermal electronic excitations, and the convergence of the energy potential surface by wave function optimization within Variational Monte Carlo and Projection Quant… ▽ More We present Coupled Electron-Ion Monte Carlo results for the principal Hugoniot of deuterium together with an accurate study of the initial reference state of shock wave experiments. We discuss the influence of nuclear quantum effects, thermal electronic excitations, and the convergence of the energy potential surface by wave function optimization within Variational Monte Carlo and Projection Quantum Monte Carlo methods. Compared to a previous study, the new calculations also include low pressure-temperature (P,T) conditions resulting in close agreement with experimental data, while our revised results at higher (P,T) conditions still predict a more compressible Hugoniot than experimentally observed. △ Less

Submitted 1 August, 2020; originally announced August 2020.

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

A note on the three dimensional Dirac operator with zigzag type boundary conditions

Authors: Markus Holzmann

Abstract: In this note the three dimensional Dirac operator $A_m$ with boundary conditions, which are the analogue of the two dimensional zigzag boundary conditions, is investigated. It is shown that $A_m$ is self-adjoint in $L^2(Ωおめが;\mathbb{C}^4)$ for any open set $Ωおめが\subset \mathbb{R}^3$ and its spectrum is described explicitly in terms of the spectrum of the Dirichlet Laplacian in $Ωおめが$. In particular, whenev… ▽ More In this note the three dimensional Dirac operator $A_m$ with boundary conditions, which are the analogue of the two dimensional zigzag boundary conditions, is investigated. It is shown that $A_m$ is self-adjoint in $L^2(Ωおめが;\mathbb{C}^4)$ for any open set $Ωおめが\subset \mathbb{R}^3$ and its spectrum is described explicitly in terms of the spectrum of the Dirichlet Laplacian in $Ωおめが$. In particular, whenever the spectrum of the Dirichlet Laplacian is purely discrete, then also the spectrum of $A_m$ consists of discrete eigenvalues that accumulate at $\pm \infty$ and one additional eigenvalue of infinite multiplicity. △ Less

Submitted 29 January, 2021; v1 submitted 30 June, 2020; originally announced June 2020.

Comments: 13 pages; revised version

Quantum simulation and circuit design for solving multidimensional Poisson equations

Authors: Michael Holzmann, Harald Koestler

Abstract: Many methods solve Poisson equations by using grid techniques which discretize the problem in each dimension. Most of these algorithms are subject to the curse of dimensionality, so that they need exponential runtime. In the paper "Quantum algorithm and circuit design solving the Poisson equation" a quantum algorithm is shown running in polylog time to produce a quantum state representing the solu… ▽ More Many methods solve Poisson equations by using grid techniques which discretize the problem in each dimension. Most of these algorithms are subject to the curse of dimensionality, so that they need exponential runtime. In the paper "Quantum algorithm and circuit design solving the Poisson equation" a quantum algorithm is shown running in polylog time to produce a quantum state representing the solution of the Poisson equation. In this paper a quantum simulation of an extended circuit design based on this algorithm is made on a classical computer. Our purpose is to test an efficient circuit design which can break the curse of dimensionality on a quantum computer. Due to the exponential rise of the Hilbert space this design is optimized on a small number of qubits. We use Microsoft's Quantum Development Kit and its simulator of an ideal quantum computer to validate the correctness of this algorithm. △ Less

Submitted 16 June, 2020; originally announced June 2020.

Observation of Algebraic Time Order for Two-Dimensional Dipolar Excitons

Authors: S. Dang, M. Zamorano, S. Suffit, K. West, K. Baldwin, L. Pfeiffer, M. Holzmann, F. Dubin

Abstract: Emergence of algebraic quasi-long-range order is a key feature of superfluid phase transitions at two dimensions. For this reduced dimensionality interactions prevent Bose-Einstein condensation with true long range order, at any finite temperature. Here, we report the occurence of algebraic order in a strongly interacting quantum liquid formed by dipolar excitons confined in a bilayer semiconducto… ▽ More Emergence of algebraic quasi-long-range order is a key feature of superfluid phase transitions at two dimensions. For this reduced dimensionality interactions prevent Bose-Einstein condensation with true long range order, at any finite temperature. Here, we report the occurence of algebraic order in a strongly interacting quantum liquid formed by dipolar excitons confined in a bilayer semiconductor heterostructure. We observe a transition from exponential to algebraic decay of the excitons temporal coherence, accompanied by a universal scaling behaviour of the equation of state. Our results provide strong evidence for a Berezinskii-Kosterlitz-Thouless (BKT) transition in a multi-component boson-like system governed by strong dipolar interactions. △ Less

Submitted 5 January, 2020; originally announced January 2020.

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

Quantum Monte Carlo Compton profiles of solid and liquid lithium

Authors: Yubo Yang, Nozomu Hiraoka, Kazuhiro Matsuda, Markus Holzmann, David M. Ceperley

Abstract: We computed the Compton profile of solid and liquid lithium using quantum Monte Carlo (QMC) and compared with recent experimental measurements obtaining good agreement. Importantly, we find it crucial to account for proper core-valence orthogonalization and to address density differences when comparing with experiment. To account for disorder effects, we sampled finite-temperature configurations u… ▽ More We computed the Compton profile of solid and liquid lithium using quantum Monte Carlo (QMC) and compared with recent experimental measurements obtaining good agreement. Importantly, we find it crucial to account for proper core-valence orthogonalization and to address density differences when comparing with experiment. To account for disorder effects, we sampled finite-temperature configurations using molecular dynamics (MD), then performed diffusion Monte Carlo (DMC) simulations on each configuration. We used Slater-Jastrow wavefunctions and grand-canonical twist-averaged boundary conditions. A QMC pseudopotential correction, derived from an all-electron DMC simulation of the perfect crystal was also used. Our calculations provide the first all-electron QMC benchmark for the Compton profile of lithium crystal and pseudopotential-corrected QMC Compton profiles for both the liquid and solid. △ Less

Submitted 27 December, 2019; originally announced December 2019.

Comments: 9 pages, 8 figures; supplemental materials 6 pages, 4 figures

Journal ref: Phys. Rev. B 101, 165125 (2020)

Computing general observables in lattice models with complex actions

Authors: Olmo Francesconi, Markus Holzmann, Biagio Lucini, Antonio Rago, Jarno Rantaharju

Abstract: The study of QFTs at finite density is hindered by the presence of the so-called sign problem. The action definition of such systems is, in fact, complex-valued making standard importance sampling Monte Carlo methods ineffective. In this work, we shall review the generalized density of states method for complex action systems and the Linear Logarithmic Relaxation algorithm (LLR). We will focus on… ▽ More The study of QFTs at finite density is hindered by the presence of the so-called sign problem. The action definition of such systems is, in fact, complex-valued making standard importance sampling Monte Carlo methods ineffective. In this work, we shall review the generalized density of states method for complex action systems and the Linear Logarithmic Relaxation algorithm (LLR). We will focus on the recent developments regarding the bias control of the LLR method and the evaluation of general observables in the DoS+LLR framework. Recent results on the well-known relativistic Bose gas will be presented, proving that in our approach the phase factor can be consistently evaluated over hundreds of orders of magnitude. A first exploratory study on the Thirring model in the DoS formalism will be presented as well. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: 7 pages, 4 figures, presented at the 37th International Symposium on Lattice Field Theory - Lattice2019, 16-22 June 2019, Wuhan, China

Energy gap closure of crystalline molecular hydrogen with pressure

Authors: Vitaly Gorelov, Markus Holzmann, David M. Ceperley, Carlo Pierleoni

Abstract: We study the gap closure with pressure of crystalline molecular hydrogen. The gaps are obtained from grand-canonical Quantum Monte Carlo methods properly extended to quantum and thermal crystals, simulated by Coupled Electron Ion Monte Carlo. Nuclear zero point effects cause a large reduction in the gap ($\sim 2eV$). \CP{Depending on the structure,} the fundamental indirect gap closes \CP{between… ▽ More We study the gap closure with pressure of crystalline molecular hydrogen. The gaps are obtained from grand-canonical Quantum Monte Carlo methods properly extended to quantum and thermal crystals, simulated by Coupled Electron Ion Monte Carlo. Nuclear zero point effects cause a large reduction in the gap ($\sim 2eV$). \CP{Depending on the structure,} the fundamental indirect gap closes \CP{between 380GPa and} 530GPa for ideal crystals and 330-380GPa for quantum crystals. Beyond this pressure the system enters into a bad metal phase where the density of states at the Fermi level increases with pressure up to $\sim$450\CP{-500} GPa when the direct gap closes. Our work partially supports the interpretation of recent experiments in high pressure hydrogen. △ Less

Submitted 15 February, 2020; v1 submitted 14 November, 2019; originally announced November 2019.

Comments: 5 pages, 5 figures, to appear on PRL

Journal ref: Phys. Rev. Lett. 124, 116401 (2020)

Superconductor-insulator transition in Josephson junction chains by quantum Monte-Carlo

Authors: D. M. Basko, F. Pfeiffer, P. Adamus, M. Holzmann, F. W. J. Hekking

Abstract: We study the zero-temperature phase diagram of a dissipationless and disorder-free Josephson junction chain. Namely, we determine the critical Josephson energy below which the chain becomes insulating, as a function of the ratio of two capacitances: the capacitance of each Josephson junction and the capacitance between each superconducting island and the ground. We develop an imaginary-time path i… ▽ More We study the zero-temperature phase diagram of a dissipationless and disorder-free Josephson junction chain. Namely, we determine the critical Josephson energy below which the chain becomes insulating, as a function of the ratio of two capacitances: the capacitance of each Josephson junction and the capacitance between each superconducting island and the ground. We develop an imaginary-time path integral Quantum Monte-Carlo algorithm in the charge representation, which enables us to efficiently handle the electrostatic part of the chain Hamiltonian. We find that a large part of the phase diagram is determined by anharmonic corrections which are not captured by the standard Kosterlitz-Thouless renormalization group description of the transition. △ Less

Submitted 7 November, 2019; originally announced November 2019.

Journal ref: Physical Review B 101 , 024518 (2020)

Self-adjoint Dirac operators on domains in $\mathbb{R}^3$

Authors: Jussi Behrndt, Markus Holzmann, Albert Mas

Abstract: In this paper the spectral and scattering properties of a family of self-adjoint Dirac operators in $L^2(Ωおめが; \mathbb{C}^4)$, where $Ωおめが\subset \mathbb{R}^3$ is either a bounded or an unbounded domain with a compact $C^2$-smooth boundary, are studied in a systematic way. These operators can be viewed as the natural relativistic counterpart of Laplacians with Robin boundary conditions. Among the Dirac… ▽ More In this paper the spectral and scattering properties of a family of self-adjoint Dirac operators in $L^2(Ωおめが; \mathbb{C}^4)$, where $Ωおめが\subset \mathbb{R}^3$ is either a bounded or an unbounded domain with a compact $C^2$-smooth boundary, are studied in a systematic way. These operators can be viewed as the natural relativistic counterpart of Laplacians with Robin boundary conditions. Among the Dirac operators treated here is also the so-called MIT bag operator, which has been used by physicists and more recently was discussed in the mathematical literature. Our approach is based on abstract boundary triple techniques from extension theory of symmetric operators and a thorough study of certain classes of (boundary) integral operators, that appear in a Krein-type resolvent formula. The analysis of the perturbation term in this formula leads to a description of the spectrum and a Birman-Schwinger principle, a qualitative understanding of the scattering properties in the case that $Ωおめが$ is unbounded, and corresponding trace formulas. △ Less

Submitted 25 October, 2019; originally announced October 2019.

Comments: 48 pages

Free energy of the self-interacting relativistic lattice Bose gas at finite density

Authors: Olmo Francesconi, Markus Holzmann, Biagio Lucini, Antonio Rago

Abstract: The density of state approach has recently been proposed as a potential route to circumvent the sign problem in systems at finite density. In this study, using the Linear Logarithmic Relaxation (LLR) algorithm, we extract the generalised density of states, which is defined in terms of the imaginary part of the action, for the self-interacting relativistic lattice Bose gas at finite density. After… ▽ More The density of state approach has recently been proposed as a potential route to circumvent the sign problem in systems at finite density. In this study, using the Linear Logarithmic Relaxation (LLR) algorithm, we extract the generalised density of states, which is defined in terms of the imaginary part of the action, for the self-interacting relativistic lattice Bose gas at finite density. After discussing the implementation and testing the reliability of our approach, we focus on the determination of the free energy difference between the full system and its phase-quenched counterpart. Using a set of lattices ranging from $4^4$ to $16^4$ , we show that in the low density phase, this overlap free energy can be reliably extrapolated to the thermodynamic limit. The numerical precision we obtain with the LLR method allows us to determine with sufficient accuracy the expectation value of the phase factor, which is used in the calculation of the overlap free energy, down to values of ${\cal O}(10^{-480})$. When phase factor measurements are extended to the dense phase, a change of behaviour of the overlap free energy is clearly visible as the chemical potential crosses a critical value. Using fits inspired by the approximate validity of mean-field theory, which is confirmed by our simulations, we extract the critical chemical potential as the non-analyticity point in the overlap free energy, obtaining a value that is in agreement with other determinations. Implications of our findings and potential improvements of our methodology are also discussed. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: 21 pages, 15 figures

Journal ref: Phys. Rev. D 101, 014504 (2020)

Electronic band gaps from Quantum Monte Carlo methods

Authors: Yubo Yang, Vitaly Gorelov, Carlo Pierleoni, David M. Ceperley, Markus Holzmann

Abstract: We develop a method for calculating the fundamental electronic gap of semiconductors and insulators using grand canonical Quantum Monte Carlo simulations. We discuss the origin of the bias introduced by supercell calculations of finite size and show how to correct the leading and subleading finite size errors either based on observables accessible in the finite-sized simulations or from DFT calcul… ▽ More We develop a method for calculating the fundamental electronic gap of semiconductors and insulators using grand canonical Quantum Monte Carlo simulations. We discuss the origin of the bias introduced by supercell calculations of finite size and show how to correct the leading and subleading finite size errors either based on observables accessible in the finite-sized simulations or from DFT calculations. Our procedure is applied to solid molecular hydrogen and compared to experiment for carbon and silicon crystals. △ Less

Submitted 2 November, 2019; v1 submitted 16 October, 2019; originally announced October 2019.

Comments: 11 pages, 5 figures; supplemental materials 21 pages

Journal ref: Phys. Rev. B 101, 085115 (2020)

Orbital-dependent backflow wave functions for real-space quantum Monte Carlo

Authors: Markus Holzmann, Saverio Moroni

Abstract: We present and motivate an efficient way to include orbital dependent many--body correlations in trial wave function of real--space Quantum Monte Carlo methods for use in electronic structure calculations. We apply our new orbital--dependent backflow wave function to calculate ground state energies of the first row atoms using variational and diffusion Monte Carlo methods. The systematic overall g… ▽ More We present and motivate an efficient way to include orbital dependent many--body correlations in trial wave function of real--space Quantum Monte Carlo methods for use in electronic structure calculations. We apply our new orbital--dependent backflow wave function to calculate ground state energies of the first row atoms using variational and diffusion Monte Carlo methods. The systematic overall gain of correlation energy with respect to single determinant Jastrow-Slater wave functions is competitive with the best single determinant trial wave functions currently available. The computational cost per Monte Carlo step is comparable to that of simple backflow calculations. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: 5 pages, 2 figures, contains additional reference [18] compared to published version

Journal ref: Phys. Rev. B 99, 085121 (2019)

Itinerant-electron magnetism: the importance of many-body correlations

Authors: Makus Holzmann, Saverio Moroni

Abstract: Do electrons become ferromagnetic just because of their repulisve Coulomb interaction? Our calculations on the three-dimensional electron gas imply that itinerant ferromagnetim of delocalized electrons without lattice and band structure, the most basic model considered by Stoner, is suppressed due to many-body correlations as speculated already by Wigner, and a possible ferromagnetic transition lo… ▽ More Do electrons become ferromagnetic just because of their repulisve Coulomb interaction? Our calculations on the three-dimensional electron gas imply that itinerant ferromagnetim of delocalized electrons without lattice and band structure, the most basic model considered by Stoner, is suppressed due to many-body correlations as speculated already by Wigner, and a possible ferromagnetic transition lowering the density is precluded by the formation of the Wigner crystal. △ Less

Submitted 2 June, 2020; v1 submitted 15 October, 2019; originally announced October 2019.

Comments: published version including supplementary material

Journal ref: Phys. Rev. Lett. 124, 206404 (2020)

Limiting absorption principle and scattering matrix for Dirac operators with $δでるた$-shell interactions

Authors: Jussi Behrndt, Markus Holzmann, Andrea Mantile, Andrea Posilicano

Abstract: We provide a limiting absorption principle for self-adjoint realizations of Dirac operators with electrostatic and Lorentz scalar $δでるた$-shell interactions supported on regular compact surfaces. Then we show completeness of the wave operators and give a representation formula for the scattering matrix. We provide a limiting absorption principle for self-adjoint realizations of Dirac operators with electrostatic and Lorentz scalar $δでるた$-shell interactions supported on regular compact surfaces. Then we show completeness of the wave operators and give a representation formula for the scattering matrix. △ Less

Submitted 12 June, 2020; v1 submitted 3 August, 2019; originally announced August 2019.

Comments: Final version

Journal ref: J. Math. Phys. 61, 033504 (2020)

Two-dimensional Dirac operators with singular interactions supported on closed curves

Authors: Jussi Behrndt, Markus Holzmann, Thomas Ourmières-Bonafos, Konstantin Pankrashkin

Abstract: We study the two-dimensional Dirac operator with an arbitrary combination of electrostatic and Lorentz scalar $δでるた$-interactions of constant strengths supported on a smooth closed curve. For any combination of the coupling constants a rigorous description of the self-adjoint realizations of the operators is given and the qualitative spectral properties are described. The analysis covers also all so-… ▽ More We study the two-dimensional Dirac operator with an arbitrary combination of electrostatic and Lorentz scalar $δでるた$-interactions of constant strengths supported on a smooth closed curve. For any combination of the coupling constants a rigorous description of the self-adjoint realizations of the operators is given and the qualitative spectral properties are described. The analysis covers also all so-called critical combinations of coupling constants, for which there is a loss of regularity in the operator domain. In this case, if the mass is non-zero, the resulting operator has an additional point in the essential spectrum, and the position of this point inside the central gap can be made arbitrary by a suitable choice of the coupling constants. The analysis is based on a combination of the extension theory of symmetric operators with a detailed study of boundary integral operators viewed as periodic pseudodifferential operators. △ Less

Submitted 8 August, 2019; v1 submitted 11 July, 2019; originally announced July 2019.

Comments: 54 pages. Some changes in the structure of the text

Journal ref: J. Funct. Anal. 279 (2020) 108700

Boundary integral formulations of eigenvalue problems for elliptic differential operators with singular interactions and their numerical approximation by boundary element methods

Authors: Markus Holzmann, Gerhard Unger

Abstract: In this paper the discrete eigenvalues of elliptic second order differential operators in $L^2(\mathbb{R}^n)$, $n \in \mathbb{N}$, with singular $δでるた$- and $δでるた'$-interactions are studied. We show the self-adjointness of these operators and derive equivalent formulations for the eigenvalue problems involving boundary integral operators. These formulations are suitable for the numerical computations of… ▽ More In this paper the discrete eigenvalues of elliptic second order differential operators in $L^2(\mathbb{R}^n)$, $n \in \mathbb{N}$, with singular $δでるた$- and $δでるた'$-interactions are studied. We show the self-adjointness of these operators and derive equivalent formulations for the eigenvalue problems involving boundary integral operators. These formulations are suitable for the numerical computations of the discrete eigenvalues and the corresponding eigenfunctions by boundary element methods. We provide convergence results and show numerical examples. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 38 pages, 13 figures

On Dirac operators in $\mathbb{R}^3$ with electrostatic and Lorentz scalar $δでるた$-shell interactions

Authors: Jussi Behrndt, Pavel Exner, Markus Holzmann, Vladimir Lotoreichik

Abstract: In this article Dirac operators $A_{ηいーた, τたう}$ coupled with combinations of electrostatic and Lorentz scalar $δでるた$-shell interactions of constant strength $ηいーた$ and $τたう$, respectively, supported on compact surfaces $Σしぐま\subset \mathbb{R}^3$ are studied. In the rigorous definition of these operators the $δでるた$-potentials are modelled by coupling conditions at $Σしぐま$. In the proof of the self-adjointness of… ▽ More In this article Dirac operators $A_{ηいーた, τたう}$ coupled with combinations of electrostatic and Lorentz scalar $δでるた$-shell interactions of constant strength $ηいーた$ and $τたう$, respectively, supported on compact surfaces $Σしぐま\subset \mathbb{R}^3$ are studied. In the rigorous definition of these operators the $δでるた$-potentials are modelled by coupling conditions at $Σしぐま$. In the proof of the self-adjointness of $A_{ηいーた, τたう}$ a Krein-type resolvent formula and a Birman-Schwinger principle are obtained. With their help a detailed study of the qualitative spectral properties of $A_{ηいーた, τたう}$ is possible. In particular, the essential spectrum of $A_{ηいーた, τたう}$ is determined, it is shown that at most finitely many discrete eigenvalues can appear, and several symmetry relations in the point spectrum are obtained. Moreover, the nonrelativistic limit of $A_{ηいーた, τたう}$ is computed and it is discussed that for some special interaction strengths $A_{ηいーた, τたう}$ is decoupled to two operators acting in the domains with the common boundary $Σしぐま$. △ Less

Submitted 6 March, 2019; v1 submitted 31 January, 2019; originally announced January 2019.

Comments: contribution to the proceedings of the conference "Advances in Operator Theory with Applications to Mathematical Physics" at the Chapman University, November 12-16, 2018; 20 pages

The density of states approach to the sign problem

Authors: Biagio Lucini, Olmo Francesconi, Markus Holzmann, Antonio Rago

Abstract: Approaches to the sign problem based on the density of states have been recently revived by the introduction of the LLR algorithm, which allows us to compute the density of states itself with exponential error reduction. In this work, after a review of the generalities of the method, we show recent results for the Bose gas in four dimensions, focussing on the identification of possible systematic… ▽ More Approaches to the sign problem based on the density of states have been recently revived by the introduction of the LLR algorithm, which allows us to compute the density of states itself with exponential error reduction. In this work, after a review of the generalities of the method, we show recent results for the Bose gas in four dimensions, focussing on the identification of possible systematic errors and on methods of controlling the bias they can introduce in the calculation. △ Less

Submitted 22 January, 2019; originally announced January 2019.

Comments: 1 + 11 pages, 6 figures. Invited talk presented by B. Lucini at the conference "XIIIth Quark Confinement and the Hadron Spectrum", Maynooth University, Ireland, 1-6 August 2018

The Landau Hamiltonian with $δでるた$-potentials supported on curves

Authors: Jussi Behrndt, Pavel Exner, Markus Holzmann, Vladimir Lotoreichik

Abstract: The spectral properties of the singularly perturbed self-adjoint Landau Hamiltonian $A_αあるふぁ=(i \nabla + A)^2 + αあるふぁδでるた$ in $L^2(R^2)$ with a $δでるた$-potential supported on a finite $C^{1,1}$-smooth curve $Σしぐま$ are studied. Here $A = \frac{1}{2} B (-x_2, x_1)^\top$ is the vector potential, $B>0$ is the strength of the homogeneous magnetic field, and $αあるふぁ\in L^\infty(Σしぐま)$ is a position-dependent real coefficient mod… ▽ More The spectral properties of the singularly perturbed self-adjoint Landau Hamiltonian $A_αあるふぁ=(i \nabla + A)^2 + αあるふぁδでるた$ in $L^2(R^2)$ with a $δでるた$-potential supported on a finite $C^{1,1}$-smooth curve $Σしぐま$ are studied. Here $A = \frac{1}{2} B (-x_2, x_1)^\top$ is the vector potential, $B>0$ is the strength of the homogeneous magnetic field, and $αあるふぁ\in L^\infty(Σしぐま)$ is a position-dependent real coefficient modeling the strength of the singular interaction on the curve $Σしぐま$. After a general discussion of the qualitative spectral properties of $A_αあるふぁ$ and its resolvent, one of the main objectives in the present paper is a local spectral analysis of $A_αあるふぁ$ near the Landau levels $B(2q+1)$. Under various conditions on $αあるふぁ$ it is shown that the perturbation smears the Landau levels into eigenvalue clusters, and the accumulation rate of the eigenvalues within these clusters is determined in terms of the capacity of the support of $αあるふぁ$. Furthermore, the use of Landau Hamiltonians with $δでるた$-perturbations as model operators for more realistic quantum systems is justified by showing that $A_αあるふぁ$ can be approximated in the norm resolvent sense by a family of Landau Hamiltonians with suitably scaled regular potentials. △ Less

Submitted 21 December, 2018; originally announced December 2018.

Comments: Merry Christmas

Two-dimensional Mixture of Dipolar Fermions: Equation of State and Magnetic Phases

Authors: Tommaso Comparin, Raul Bombin, Markus Holzmann, Ferran Mazzanti, Jordi Boronat, Stefano Giorgini

Abstract: We study a two-component mixture of fermionic dipoles in two dimensions at zero temperature, interacting via a purely repulsive $1/r^3$ potential. This model can be realized with ultracold atoms or molecules, when their dipole moments are aligned in the confinement direction orthogonal to the plane. We characterize the unpolarized mixture by means of the Diffusion Monte Carlo technique. Computing… ▽ More We study a two-component mixture of fermionic dipoles in two dimensions at zero temperature, interacting via a purely repulsive $1/r^3$ potential. This model can be realized with ultracold atoms or molecules, when their dipole moments are aligned in the confinement direction orthogonal to the plane. We characterize the unpolarized mixture by means of the Diffusion Monte Carlo technique. Computing the equation of state, we identify the regime of validity for a mean-field theory based on a low-density expansion and compare our results with the hard-disk model of repulsive fermions. At high density, we address the possibility of itinerant ferromagnetism, namely whether the ground state can be fully polarized in the fluid phase. Within the fixed-node approximation, we show that the accuracy of Jastrow-Slater trial wave functions, even with the typical two-body backflow correction, is not sufficient to resolve the relevant energy differences. By making use of the iterative-backflow improved trial wave functions, we observe no signature of a fully-polarized ground state up to the freezing density. △ Less

Submitted 9 April, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 12 pages, 6 figures

Journal ref: Phys. Rev. A 99, 043609 (2019)

Defect-Driven Superfluid Crossover for Two-Dimensional Dipolar Excitons Trapped at Thermodynamic Equilibrium

Authors: Suzanne Dang, Romain Anankine, Carmen Gomez, Aristide Lemaître, Markus Holzmann, François Dubin

Abstract: We study ultra-cold dipolar excitons confined in a 10$μみゅー$m trap of a double GaAs quantum well. Based on the local density approximation, we unveil for the first time the equation of state of excitons at pure thermodynamic equilibrium. In this regime we show that, below a critical temperature of about $1$ Kelvin, a superfluid forms in the inner region of the trap at a local exciton density… ▽ More We study ultra-cold dipolar excitons confined in a 10$μみゅー$m trap of a double GaAs quantum well. Based on the local density approximation, we unveil for the first time the equation of state of excitons at pure thermodynamic equilibrium. In this regime we show that, below a critical temperature of about $1$ Kelvin, a superfluid forms in the inner region of the trap at a local exciton density $n \sim 2-3 \, 10^{10} \text{cm}^{-2}$, encircled by a more dilute and normal component in the outer rim of the trap. Remarkably, this spatial arrangement correlates directly with the concentration of defects in the exciton density which exhibits a sudden decrease at the onset of superfluidity, thus pointing towards an underlying Berezinskii-Kosterlitz-Thouless mechanism. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Comments: 8 pages, 8 figures

Journal ref: Phys. Rev. Lett. 122, 117402 (2019)

Electron localization properties in high pressure hydrogen at the liquid-liquid phase transition by Coupled Electron-Ion Monte Carlo

Authors: Carlo Pierleoni, Giovanni Rillo, David M. Ceperley, Markus Holzmann

Abstract: We analyze in detail the electronic properties of high pressure hydrogen around the liquid-liquid phase transition based on Coupled Electron-Ion Monte Carlo calculations. Computing the off-diagonal single particle density matrix and the momentum distribution we discuss localization properties of the electrons. The abrupt changes of these distributions indicate a metal to insulator transition occur… ▽ More We analyze in detail the electronic properties of high pressure hydrogen around the liquid-liquid phase transition based on Coupled Electron-Ion Monte Carlo calculations. Computing the off-diagonal single particle density matrix and the momentum distribution we discuss localization properties of the electrons. The abrupt changes of these distributions indicate a metal to insulator transition occurring together with the structural transition from the atomic to molecular fluid. We further discuss the electron-proton and electron-electron pair correlation functions, which also change abruptly at the transition. △ Less

Submitted 1 December, 2017; originally announced December 2017.

Comments: 10 pages, 6 figures, prepared for the proceedings of CCP2017