Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userComputational efficiency of numerical integration methods for the tangent dynamics of many-body Hamiltonian systems in one and two spatial dimensions
124
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We investigate the computational performance of various numerical methods for the integration of the equations of motion and the variational equations for some typical classical many-body models of condensed matter physics: the Fermi-Pasta-Ulam-Tsingou (FPUT) chain and the one- and two-dimensional disordered, discrete nonlinear Schrodinger equations (DDNLS). In our analysis we consider methods based on Taylor series expansion, Runge-Kutta discretization and symplectic transformations. The latter have the ability to exactly preserve the symplectic structure of Hamiltonian systems, which results in keeping bounded the error of the systems computed total energy. We perform extensive numerical simulations for several initial conditions of the studied models and compare the numerical efficiency of the used integrators by testing their ability to accurately reproduce characteristics of the systems dynamics and quantify their chaoticity through the computation of the maximum Lyapunov exponent. We also report the expressions of the implemented symplectic schemes and provide the explicit forms of the used differential operators. Among the tested numerical schemes the symplectic integrators $ABA864$ and $SRKN^a_{14}$ exhibit the best performance, respectively for moderate and high accuracy levels in the case of the FPUT chain, while for the DDNLS models $s9mathcal{ABC}6$ and $s11mathcal{ABC}6$ (moderate accuracy), along with $s17mathcal{ABC}8$ and $s19mathcal{ABC}8$ (high accuracy) proved to be the most efficient schemes.
rate research
Read More
Many-body localization is a striking mechanism that prevents interacting quantum systems from thermalizing. The absence of thermalization behaviour manifests itself, for example, in a remanence of local particle number configurations, a quantity that is robust over a parameter range. Local particle numbers are directly accessible in programmable quantum simulators, in systems of cold atoms even in two spatial dimensions. Yet, the classical simulation aimed at building trust in quantum simulations is highly challenging. In this work, we present a comprehensive tensor network simulation of a many-body localized systems in two spatial dimensions using a variant of an iPEPS algorithm. The required translational invariance can be restored by implementing the disorder into an auxiliary spin system, providing an exact disorder average under dynamics. We can quantitatively assess signatures of many-body localization for the infinite system: Our methods are powerful enough to provide crude dynamical estimates for the transition between localized and ergodic phases. Interestingly, in this setting of finitely many disorder values, which we also compare with simulations involving non-interacting fermions and for which we discuss the emergent physics, localization emerges in the interacting regime, for which we provide an intuitive argument, while Anderson localization is absent.
A powerful perspective in understanding non-equilibrium quantum dynamics is through the time evolution of its entanglement content. Yet apart from a few guiding principles for the entanglement entropy, to date, not much else is known about the refined characters of entanglement propagation. Here, we unveil signatures of the entanglement evolving and information propagation out-of-equilibrium, from the view of entanglement Hamiltonian. As a prototypical example, we study quantum quench dynamics of a one-dimensional Bose-Hubbard model by means of time-dependent density-matrix renormalization group simulation. Before reaching equilibration, it is found that a current operator emerges in entanglement Hamiltonian, implying that entanglement spreading is carried by particle flow. In the long-time limit subsystem enters a steady phase, evidenced by the dynamic convergence of the entanglement Hamiltonian to the expectation of a thermal ensemble. Importantly, entanglement temperature of steady state is spatially independent, which provides an intuitive trait of equilibrium. We demonstrate that these features are consistent with predictions from conformal field theory. These findings not only provide crucial information on how equilibrium statistical mechanics emerges in many-body dynamics, but also add a tool to exploring quantum dynamics from perspective of entanglement Hamiltonian.
We provide a systematic comparison of two numerical methods to solve the widely used nonlinear Schrodinger equation. The first one is the standard second order split-step (SS2) method based on operator splitting approach. The second one is the Hamiltonian integration method (HIM). It allows the exact conservation of the Hamiltonian at the cost of requiring the implicit time stepping. We found that numerical error for HIM method is systematically smaller than the SS2 solution for the same time step. At the same time, one can take orders of magnitude larger time steps in HIM compared with SS2 still ensuring numerical stability. In contrast, SS2 time step is limited by the numerical stability threshold.
In this Thesis, we report a detailed study of the ground-state properties of a set of quantum few- and many-body systems in one and two dimensions with different types of interactions by using Quantum Monte Carlo methods. Nevertheless, the main focus of this work is the study of the ground-state properties of an ultracold Bose system with dipole-dipole interaction between the particles. We consider the cases where the bosons are confined to a bilayer and multilayer geometries, that consist of equally spaced two-dimensional layers. These layers can be experimentally realized by imposing tight confinement in one direction. We specifically address the study of new quantum phases, their properties, and transitions between them. One expects these systems to have a rich collection of few- and many-body phases because the dipole-dipole interaction is anisotropic and quasi long-range.
We investigate the achievable efficiency of both the time and the space discretisation methods used in Antares for mixed parabolic-hyperbolic problems. We show that the fifth order variant of WENO combined with a second order Runge-Kutta scheme is not only more accurate than standard first and second order schemes, but also more efficient taking the computation time into account. Then, we calculate the error decay rates of WENO with several explicit Runge-Kutta schemes for advective and diffusive problems with smooth and non-smooth initial conditions. With this data, we estimate the computational costs of three-dimensional simulations of stellar surface convection and show that SSP RK(3,2) is the most efficient scheme considered in this comparison.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
