PIF: Quantum Monte Carlo Methods for Non-Equilibrium Dynamics of Interacting Quantum Many-Body Systems

PIF：相互作用量子多体系统非平衡动力学的量子蒙特卡罗方法

• 批准号: 1211284
• 负责人: Anders Sandvik
• 金额: $ 34.5万
• 依托单位: Trustees of Boston University
• 依托单位国家: 美国
• 项目类别: Continuing Grant
• 财政年份: 2012
• 资助国家: 美国
• 起止时间: 2012-08-15 至 2015-10-31
• 项目状态: 已结题
• 来源: https://www.nsf.gov/awardsearch/showAward?AWD_ID=1211284&HistoricalAwards=false
• 关键词: PIF; Quantum; Monte; Carlo; Methods

In this project, novel computational and theoretical approaches to the dynamics of quantum many-body systems are investigated. The central idea is that quantum quenches using essentially arbitrary evolution protocols in imaginary time can be carried out numerically using modified versions of existing quantum Monte Carlo (QMC) algorithms. We develop two algorithms: (i) In non-equilibrium QMC (NEQMC) simulations, the final state of a system after a Hamiltonian evolution in imaginary time is computed. (ii) In the quasi-adiabatic QMC (QAQMC) method, expectation values along a full time-path are obtained simultaneously, using a product of a large number of evolving Hamiltonians (instead of the standard exponential time evolution operator). Key to the utility of these approaches are theoretical insights into how the imaginary-time information is related to real-times dynamics. Specifically, scaling behaviors as a function of a quench velocity (or a generalized velocity in the case of nonlinear protocols) in the neighborhood of a quantum-critical point are investigated. Susceptibilities characterizing non-equilibrium dynamics are also obtained, including the geometric tensor characterizing the state space. These theoretical ideas are developed in parallel with the QMC algorithms. Several applications to test the methods and demonstrate their uses are considered, including quantum phase transitions in transverse-field Ising models and Heisenberg-type spin models, as well as quantum spin glasses and other disordered systems. Quantum annealing protocols of interest in the context of quantum computations and other applications of adiabatic or quasi-adiabatic evolution of quantum systems are also investigated. In addition, some of the ideas developed within the context of quantum systems are also applied to classical many-body dynamics, e.g., to obtain improved ways of computing the dynamic critical exponent.One of the greatest challenges of contemporary theoretical physics is how to compute (predict) the dynamical evolution of systems of a large (macroscopic) number of interacting microscopical particles. Examples of such systems include the electrons in a solid or confined "clouds" of atoms cooled to ultra-cold temperature. These particles obey the laws of quantum mechanics, which makes computations of their behaviors extremely difficult, especially as regards their evolution in time (quantum dynamics). Even with today's supercomputers, there are no generally applicable numerical algorithms useful for solving this class of problems within reasonable time (the computation time typically scaling exponentially with the number of particles in the system). This project is centered around an idea to partially circumvent these problems for some classes of important quantum systems. Using a mathematical trick of "rotating" the time dimension in the complex plane to the imaginary-time axis, certain well known simulation algorithms for quantum systems in equilibrium can be modified to solve the evolution as a function of imaginary time. In parallel with the development of these computational tools, related theoretical work is conducted in order to obtain precise relationships between real and imaginary time dynamics. The systems under investigation are of interest, e.g., in the fields of quantum magnetism (magnetism at the electronic scale) and ultra-cold atoms, and it can be expected that the results can have impact also broadly beyond these systems (as the quantum dynamics problem is of very general interest in physics), e.g., in quantum computation (currently researched computers based on quantum-mechanical principles). The graduate students involved in the project receive training in cutting-edge scientific computation and theoretical physics.

在这个项目中，我们研究了量子多体系统动力学的新的计算和理论方法。其核心思想是，在虚时间内使用基本上任意演化协议的量子猝灭可以使用现有量子蒙特卡罗(QMC)算法的修改版本进行数值计算。我们提出了两个算法：(I)在非平衡量子力学(NEQMC)模拟中，计算系统在虚时间哈密顿演化后的最终状态。(Ii)在准绝热QMC(QAQMC)方法中，使用大量演化哈密顿量的乘积(而不是标准的指数时间演化算子)同时获得沿整个时间路径的期望值。这些方法的实用性的关键是对想象的时间信息如何与实时动态相关的理论洞察。具体地，研究了在量子临界点附近作为失超速度(或在非线性协议情况下的广义速度)的函数的标度行为。还得到了表征非平衡动力学的极化率，包括表征状态空间的几何张量。这些理论思想是与QMC算法并行发展的。我们考虑了一些应用来测试这些方法并展示它们的用途，包括横场伊辛模型和海森堡型自旋模型中的量子相变，以及量子自旋玻璃和其他无序系统。在量子计算和量子系统的绝热或准绝热演化的其他应用的背景下，我们也研究了感兴趣的量子退火协议。此外，在量子系统背景下发展起来的一些想法也被应用到经典的多体动力学中，例如，获得计算动力学临界指数的改进方法。当代理论物理的最大挑战之一是如何计算(预测)大量(宏观)相互作用的微观粒子系统的动力学演化。这类系统的例子包括固体中的电子或冷却到超低温的受限原子云中的电子。这些粒子遵守量子力学定律，这使得计算它们的行为极其困难，特别是关于它们在时间上的演化(量子动力学)。即使在今天的超级计算机上，也没有普遍适用的数值算法可以在合理的时间内解决这类问题(计算时间通常随系统中粒子的数量呈指数级增长)。这个项目围绕着一个想法，即部分绕过某些类别的重要量子系统的这些问题。利用将复平面中的时间维度旋转到虚时间轴的数学技巧，可以修改某些著名的平衡量子系统的模拟算法，以解决作为虚时间函数的演化问题。在开发这些计算工具的同时，还进行了相关的理论工作，以获得真实和虚构的时间动力学之间的精确关系。正在研究的系统是令人感兴趣的，例如，在量子磁学(电子尺度的磁性)和超冷原子领域，可以预期，结果也可能产生广泛的影响，超出这些系统(因为量子动力学问题在物理学中非常普遍的兴趣)，例如，在量子计算(当前基于量子力学原理研究的计算机)。参与该项目的研究生接受了尖端科学计算和理论物理方面的培训。