Related papers: Simulating the Generalized Gibbs Ensemble (GGE): a…
We develop Hybrid Monte Carlo (HMC) algorithms for constrained Hamiltonian systems of gauge- Higgs models and introduce a new observable for the constraint effective Higgs potential. We use an extension of the so-called Rattle algorithm to…
The numerical simulation of quantum many-body dynamics is typically limited by the linear growth of entanglement with time. Recently numerical studies have shown, however, that for 1D Bethe-integrable models the simulation of local…
The Gibbs sampler (a.k.a. Glauber dynamics and heat-bath algorithm) is a popular Markov Chain Monte Carlo algorithm which iteratively samples from the conditional distributions of a probability measure $\pi$ of interest. Under the…
After a quench, observables in an integrable system may not relax to the standard thermal values, but can relax to the ones predicted by the generalized Gibbs ensemble (GGE) [M. Rigol et al., Phys. Rev. Lett. 98, 050405 (2007)]. The GGE has…
The widespread popularity of replica exchange and expanded ensemble algorithms for simulating complex molecular systems in chemistry and biophysics has generated much interest in enhancing phase space mixing of these protocols, thus…
We study a sequential Monte Carlo algorithm to sample from the Gibbs measure with a non-convex energy function at a low temperature. We use the practical and popular geometric annealing schedule, and use a Langevin diffusion at each…
The presence of erratic or unstable paths in standard kinetic Monte Carlo simulations significantly undermines the accurate simulation and sampling of transition pathways. While typically reliable methods, such as the Gillespie algorithm,…
We study the approximation of expectations $\E(f(X))$ for Gaussian random elements $X$ with values in a separable Hilbert space $H$ and Lipschitz continuous functionals $f \colon H \to \R$. We consider restricted Monte Carlo algorithms,…
We investigate off-diagonal matrix elements of local operators in integrable spin chains, focusing on the isotropic spin-$1/2$ Heisenberg chain ($XXX$ chain). We employ state-of-the-art Algebraic Bethe Ansatz results, which allow us to…
It has been recently observed for a particular quantum quench in the XXZ spin chain that local observables do not equilibrate to the predictions of the Generalized Gibbs Ensemble (GGE). In this work we argue that the breakdown of the GGE…
Monte Carlo simulation provides a powerful tool for understanding and exploring thermodynamic phase equilibria in many-particle interacting systems. Among the most physically intuitive simulation methods is Gibbs ensemble Monte Carlo…
We numerically study the unitary time evolution of a nonintegrable model of hard-core bosons with an extensive number of local Z2 symmetries. We find that the expectation values of local observables in the stationary state are described…
The importance-sampling Monte Carlo algorithm appears to be the universally optimal solution to the problem of sampling the state space of statistical mechanical systems according to the relative importance of configurations for the…
Using a random-matrix approach and Monte-Carlo simulations, we generate scattering matrices and cross sections for compound-nucleus reactions. In the absence of direct reactions we compare the average cross sections with the analytic…
Coulomb and log-gases are exchangeable singular Boltzmann-Gibbs measures appearing in mathematical physics at many places, in particular in random matrix theory. We explore experimentally an efficient numerical method for simulating such…
Monte Carlo methods are essential tools for Bayesian inference. Gibbs sampling is a well-known Markov chain Monte Carlo (MCMC) algorithm, extensively used in signal processing, machine learning, and statistics, employed to draw samples from…
We implement fully the algebraic Bethe ansatz for the XXX Heisenberg spin chain in the case when both boundary matrices can be brought to the upper-triangular form. We define the Bethe vectors which yield the strikingly simple expression…
We propose a generalization of the algebraic Bethe ansatz to obtain the eigenvectors of the Heisenberg spin chain with general boundaries associated to the eigenvalues and the Bethe equations found recently by Cao et al. The ansatz takes…
A quantum Monte Carlo algorithm is constructed starting from the standard perturbation expansion in the interaction representation. The resulting configuration space is strongly related to that of the Stochastic Series Expansion (SSE)…
We introduce and study a class of two-dimensional integrable quantum field theories that carry an internal $\mathbb{Z}_n$ structure. These models extend factorised scattering beyond the conventional framework, featuring both the usual…