Search

We analyze in detail a sharp transition between the low-energy, low-dimensional eigenstates and the high-energy chaotic bulk of the spectrum for a simple supersymmetric quantum-mechanical model with Hamiltonian $\hat{H}_S = \left(\hat{p}_1^2 + \hat{p}_2^2 + \hat{x}_1^2 \, \hat{x}_2^2 \right) \otimes I + \hat{x}_1 \otimes \sigma_1 + \hat{x}_2 \otimes \sigma_3$, which mimics the structure of the Banks-Fischler-Susskind-Stanford (BFSS) matrix model, the spatially compactified $\mathcal{N} = 1$ super-Yang-Mills theory. We conjecture that this transition might be similar to the transition between the $D0$-brane and $M$-theory regimes in the BFSS model, and find that it does not lead to irregularities in the thermodynamic equation of state. We demonstrate that real-time spectral form-factor for our supersymmetric model exhibits the ``ramp'' behavior typical for quantum chaos. We also analyze the entanglement entropy and the spectrum of the reduced density matrix of the eigenstates of $\hat{H}_S$, considering one of the bosonic degrees of freedom as a subsystem. The entanglement entropy of low-energy eigenstates appears to be practically energy-independent. Exactly at the onset of random-matrix-type level spacing fluctuations, this behavior rapidly changes into a steady growth of entanglement with energy. We demonstrate that the spectrum of the reduced density matrix also exhibits universal level-spacing fluctuations towards its higher end, even for the ground state of the supersymmetric model. Thus even the regularly spaced, non-chaotic eigenstates contain some information about semi-classical chaotic dynamics at high energies., Comment: Proceedings of the 39th International Symposium on Lattice Field Theory (Lattice 2022), 8-12 August 2022, Bonn, Germany. v2: updated references, published version. 17 pages, 6 figures, PoS style

We analyze in detail a sharp transition between the low-energy, low-dimensional eigenstates and the high-energy chaotic bulk of the spectrum for a simple supersymmetric quantum-mechanical model with Hamiltonian $\hat{H}_S = \left(\hat{p}_1^2 + \hat{p}_2^2 + \hat{x}_1^2 \, \hat{x}_2^2 \right) \otimes I + \hat{x}_1 \otimes \sigma_1 + \hat{x}_2 \otimes \sigma_3$, which mimics the structure of the Banks-Fischler-Susskind-Stanford (BFSS) matrix model, the spatially compactified $\mathcal{N} = 1$ super-Yang-Mills theory. We conjecture that this transition might be similar to the transition between the $D0$-brane and $M$-theory regimes in the BFSS model, and find that it does not lead to irregularities in the thermodynamic equation of state. We demonstrate that real-time spectral form-factor for our supersymmetric model exhibits the ``ramp'' behavior typical for quantum chaos. We also analyze the entanglement entropy and the spectrum of the reduced density matrix of the eigenstates of $\hat{H}_S$, considering one of the bosonic degrees of freedom as a subsystem. The entanglement entropy of low-energy eigenstates appears to be practically energy-independent. Exactly at the onset of random-matrix-type level spacing fluctuations, this behavior rapidly changes into a steady growth of entanglement with energy. We demonstrate that the spectrum of the reduced density matrix also exhibits universal level-spacing fluctuations towards its higher end, even for the ground state of the supersymmetric model. Thus even the regularly spaced, non-chaotic eigenstates contain some information about semi-classical chaotic dynamics at high energies., Comment: Proceedings of the 39th International Symposium on Lattice Field Theory (Lattice 2022), 8-12 August 2022, Bonn, Germany. v2: updated references, published version. 17 pages, 6 figures, PoS style

We put the Density-of-States (DoS) approach to Monte-Carlo (MC) simulations under a stress test by applying it to a physical problem with the worst possible sign problem: the real time evolution of a non-integrable quantum spin chain. Benchmarks against numerical exact diagonalisation and stochastic reweighting are presented. Both MC methods, the DoS approach and reweighting, allow for simulations of spin chains as long as $L=40$, far beyond exact diagonalisability, though only for short evolution times $t\lesssim 1$. We identify discontinuities of the density of states as one of the key problems in the MC simulations and propose to calculate some of the dominant contributions analytically, increasing the precision of our simulations by several orders of magnitude. Even after these improvements the density of states is found highly non-smooth and therefore the DoS approach cannot outperform reweighting. We prove this implication theoretically and provide numerical evidence, concluding that the DoS approach is not well suited for quantum real time simulations with discrete degrees of freedom., Comment: 16 + 4 pages, 7 figures; code and data available (DOI: 10.5281/zenodo.7164902)

Wang-Landau simulations offer the possibility to integrate explicitly over a collective coordinate and stochastically over the remainder of configuration space. We propose to choose the so-called "slow mode", which is responsible for large autocorrelation times and thus critical slowing down, for collective integration. We study this proposal for the Ising model and the linear-log-relaxation (LLR) method as simulation algorithm. We firstly demonstrate super critical slowing down in a phase with spontaneously broken symmetry and for the heatbath algorithms, for which autocorrelation times grow exponentially with system size. By contrast, using the magnetisation as collective coordinate, we present evidence that super critical slowing down is absent. We still observe a polynomial increase of the autocorrelation time with volume (critical slowing down), which is however reduced by orders of magnitude when compared to local update techniques., Comment: Author accepted version for publication in Physical Review E, 15 pages, 15 figures

Wang-Landau simulations offer the possibility to integrate explicitly over a collective coordinate and stochastically over the remainder of configuration space. We propose to choose the so-called "slow mode", which is responsible for large autocorrelation times and thus critical slowing down, for collective integration. We study this proposal for the Ising model and the LLR method as simulation algorithm. We firstly show that in a phase with spontaneously broken global symmetry, autocorrelation times grow exponentially with system size for the standard heatbath update. Identifying the magnetisation as collective coordinate, we present evidence that critical slowing down is absent for almost all observables., Comment: 15 pages, 15 figures

Hubbard-type models on the hexagonal lattice are of great interest, as they provide realistic descriptions of graphene and other related materials. Hybrid Monte Carlo simulations offer a first-principles approach to study their phase structure. Here, we review the present status of our work in this direction., Comment: Contributed to the 8th International Conference on New Frontiers in Physics (ICNFP 2019). Revised version, accepted for publication

Hubbard-type models on the hexagonal lattice are of great interest, as they provide realistic descriptions of graphene and other related materials. Hybrid Monte Carlo simulations offer a first-principles approach to study their phase structure. Here, we review the present status of our work in this direction., Comment: Contributed to the 8th International Conference on New Frontiers in Physics (ICNFP 2019). Revised version, accepted for publication

Using state of the art Hybrid-Monte-Carlo (HMC) simulations we carry out an unbiased study of the competition between spin-density wave (SDW) and charge-density wave (CDW) order in suspended graphene. We determine that the realistic inter-electron potential of graphene must be scaled up by a factor of roughly 1.6 to induce a semimetal-SDW phase transition and find no evidence for CDW order. A study of critical properties suggests that the universality class of the three-dimensional chiral Heisenberg Gross-Neveu model with two fermion flavors, predicted by renormalization group studies and strong-coupling expansion, is unlikely to apply to this transition. We propose that our results instead favor an interpretation in terms of a conformal phase transition. In addition, we describe a variant of the HMC algorithm which uses exact fermionic forces during molecular dynamics trajectories and avoids the use of pseudofermions. Compared to standard HMC this allows for a substantial increase of the integrator stepsize while achieving comparable Metropolis acceptance rates and leads to a sizable performance improvement., Comment: revised version, accepted for publication in Phys. Rev. B

We consider the second R\'enyi entropy $S^{(2)}$ in pure lattice gauge theory with $SU(2)$, $SU(3)$ and $SU(4)$ gauge groups, which serves as a first approximation for the entanglement entropy and the entropic $C$-function. We compare the results for different gauge groups using scale setting via the string tension. We confirm that at small distances $l$ our approximation for the entropic $C$-function $C(l)$, calculated for the slab-shaped entangled region of width $l$, scales as $N_c^2 - 1$ in accordance with its interpretation in terms of free gluons. At larger distances $l$ $C(l)$ is found to approach zero for $N_c = 3, 4$, somewhat more rapidly for $N_c = 4$ than for $N_c = 3$. This finding supports the conjectured discontinuity of the entropic $C$-function in the large-$N$ limit, which was found in the context of AdS/CFT correspondence and which can be interpreted as transition between colorful quarks and gluons at small distances and colorless confined states at long distances. On the other hand, for $SU(2)$ gauge group the long-distance behavior of the entropic $C$-function is inconclusive so far. There exists a small region of lattice spacings yielding results consistent with $N_c=3,4$, while results from other lattice spacings deviate without clear systematics. We discuss several possible causes for discrepancies between our results and the behavior of entanglement entropy in holographic models., Comment: 10 pages, 4 figures; v2: journal version, some discussion and references added

We present a non-iterative solver based on the Schur complement method for sparse linear systems of special form which appear in Quantum Monte-Carlo (QMC) simulations of strongly interacting fermions on the lattice. While the number of floating-point operations for this solver scales as the cube of the number of lattice sites, for practically relevant lattice sizes it is still significantly faster than iterative solvers such as the Conjugate Gradient method in the regime of strong inter-fermion interactions, for example, in the vicinity of quantum phase transitions. The speed-up is even more dramatic for the solution of multiple linear systems with different right-hand sides. We present benchmark results for QMC simulations of the tight-binding models on the hexagonal graphene lattice with on-site (Hubbard) and non-local (Coulomb) interactions, and demonstrate the potential for further speed-up using GPU., Comment: 23 pages, 5 figures. Journal version: accepted for publication in Computer Physics Communications

Using first-principle Hybrid-Monte-Carlo (HMC) simulations, we carry out an unbiased study of the competition between spin-density wave (SDW) and charge-density wave (CDW) order in the extended Hubbard model on the two dimensional hexagonal lattice at half filling. We determine the phase diagram in the space of on-site and nearest-neighbor couplings $U$ and $V$ in the region $V