New trends in density matrix renormalization

The density matrix renormalization group (DMRG) has become a powerful numerical method that can be applied to low-dimensional strongly correlated fermionic and bosonic systems. It allows for a very precise calculation of static, dynamic and thermodynamic properties. Its field of applicability has now extended beyond condensed matter, and is successfully used in quantum chemistry, statistical mechanics, quantum information theory, and nuclear and high-energy physics as well. In this article, we briefly review the main aspects of the method and present some of the most relevant applications so as to give an overview of the scope and possibilities of DMRG. We focus on the most important extensions of the method such as the calculation of dynamical properties, the application to classical systems, finite-temperature simulations, phonons and disorder, field theory, time-dependent properties and the ab initio calculation of electronic states in molecules. The recent quantum information interpretation, the development of highly accurate time-dependent algorithms and the possibility of using the DMRG as the impurity-solver of the dynamical mean field method (DMFT) give new insights into its present and potential uses. We review the numerous very recent applications of these techniques where the DMRG has shown to be one of the most reliable and versatile methods in modern computational physics.     

基于函数矩阵的一类混沌系统同步

研究了一类混沌系统的同步问题、基于稳定性理论和极点配置技术,设计了两个混沌系统之间的同步方案,实现两个混沌系统之间的同步,通过函数矩阵,实现驱动系统和响应系统的状态变量按给定的函数矩阵同步,同时证明了该方法同样适用于两个混沌系统之间的滞后同步,通过对Lorenz混沌系统和Lorenz超混沌系统的数值模拟,进一步验证了所提方案的有效性。     

Exact analytical results for density profile in Fourier space and elastic scattering function of a rotating harmonically confined ultra-cold Fermi gas

In this paper, the system dealt with consisting of an ultra-cold neutral spin-polarized Fermi gas undergoing rotation (or in the so-called synthetic magnetic field) trapped by an anisotropic harmonic potential in a two and three-dimensional space at zero temperature. Using the so-called Bloch propagator as a tool, we derive exact closed-form expressions for particle density in Fourier space which are valid for an arbitrary particle number confined by a two and three-dimensional rotating anisotropic harmonic trap. Numerical illustrations and discussions are presented. The results can be easily generalized at finite temperatures. The crossover from two-dimensional to the one-dimensional regime is shown to be reflected in the shape of the density distribution in Fourier space at very fast rotating velocity (or at strong synthetic magnetic field). In addition, an exact analytical expression of the elastic scattering factor is found, a quantity of interest used to probe the spatial distribution of the quantum gases.     

Exact steady-state density operator for a collective atomic system in an external field

An exact steady-state density operator is obtained for a model describing the collective behaviour of a system of N two-level atoms driven by a classical field. This is used to obtain the exact steady-state expectation value of the atomic population difference for any N.     

Two-particle density matrix calculation for atom by modified perturbation theory

The two-particle density matrix of hydrogen atoms moving in a solid film is calculated. The probability of surviving a bound state between electron and proton in the course of the atomic motion in the film is discussed. The surviving fraction of the initial beam possesses the increased energy of the Coulomb attraction and lowered energy losses. This fraction can be responsible for the observed number of Rydberg states in the exit beam.     

Classical density functional theory meets Monte-Carlo simulations

ABSTRACT

Typically the quality of an approximate density functional is evaluated by a direct comparison of its predictions in a given test case to exact data obtained by computer simulations. An important example for such an approach is the test of equilibrium structure of a simple fluid as measured by the pair distribution function g(r) or the cavity correlation function y(r). However, the combination of exact density profiles and the analytical structure of density functional theory allows one to determine and potentially improve the quality of a functional in a more sophisticatedway.     

Microscopically-constrained Fock energy density functionals from chiral effective field theory. I. Two-nucleon interactions

The density matrix expansion (DME) of Negele and Vautherin is a convenient tool to map finite-range physics associated with vacuum two- and three-nucleon interactions into the form of a Skyrme-like energy density functional (EDF) with density-dependent couplings. In this work, we apply the improved formulation of the DME proposed recently in arXiv:0910.4979 by Gebremariam et al. to the non-local Fock energy obtained from chiral effective field theory (EFT) two-nucleon (NN) interactions at next-to-next-to-leading-order (N²LO). The structure of the chiral interactions is such that each coupling in the DME Fock functional can be decomposed into a coupling constant arising from zero-range contact interactions and a coupling function of the density arising from the universal long-range pion exchanges. This motivates a new microscopically-guided Skyrme phenomenology where the density-dependent couplings associated with the underlying pion-exchange interactions are added to standard empirical Skyrme functionals, and the density-independent Skyrme parameters subsequently refit to data. A link to a downloadable Mathematica notebook containing the novel density-dependent couplings is provided.     

Transfer matrix spectrum for the finite-width Ising model with adjustable boundary conditions: Exact solution

Using the spinor approach, we calculate exactly the complete spectrum of the transfer matrix for the finite-width, planar Ising model with adjustable boundary conditions. Specifically, in order to control the boundary conditions, we consider an Ising model wrapped around the cylinder, and introduce along the axis a seam of defect bonds of variable strength. Depending on the boundary conditions used, the mass gap is found to vanish algebraically or exponentially with the size of the system. These results are compared with recent numerical simulations, and with random-walk and capillary-wave arguments.     

密度矩阵重正化群的异构并行优化

密度矩阵重正化群方法（DMRG）在求解一维强关联格点模型的基态时可以获得较高的精度，在应用于二维或准二维问题时，要达到类似的精度通常需要较大的计算量与存储空间.本文提出一种新的DMRG异构并行策略，可以同时发挥计算机中央处理器（CPU）和图形处理器（GPU）的计算性能.针对最耗时的哈密顿量对角化部分，实现了数据的分布式存储，并且给出了CPU和GPU之间的负载平衡策略.以费米Hubbard模型为例，测试了异构并行程序在不同DMRG保留状态数下的运行表现，并给出了相应的性能基准.应用于4腿梯子时，观测到了高温超导中常见的电荷密度条纹，此时保留状态数达到10⁴，使用的GPU显存小于12 GB.     

密度矩阵重正化群的异构并行优化

密度矩阵重正化群方法（DMRG）在求解一维强关联格点模型的基态时可以获得较高的精度，在应用于二维或准二维问题时，要达到类似的精度通常需要较大的计算量与存储空间.本文提出一种新的DMRG异构并行策略，可以同时发挥计算机中央处理器（CPU）和图形处理器（GPU）的计算性能.针对最耗时的哈密顿量对角化部分，实现了数据的分布式存储，并且给出了CPU和GPU之间的负载平衡策略.以费米Hubbard模型为例，测试了异构并行程序在不同DMRG保留状态数下的运行表现，并给出了相应的性能基准.应用于4腿梯子时，观测到了高温超导中常见的电荷密度条纹，此时保留状态数达到10⁴，使用的GPU显存小于12 GB.     

密度矩阵重正化群的异构并行优化

密度矩阵重正化群方法（DMRG）在求解一维强关联格点模型的基态时可以获得较高的精度,在应用于二维或准二维问题时,要达到类似的精度通常需要较大的计算量与存储空间.本文提出一种新的DMRG异构并行策略,可以同时发挥计算机中央处理器（CPU）和图形处理器（GPU）的计算性能.针对最耗时的哈密顿量对角化部分,实现了数据的分布式存储,并且给出了CPU和GPU之间的负载平衡策略.以费米Hubbard模型为例,测试了异构并行程序在不同DMRG保留状态数下的运行表现,并给出了相应的性能基准.应用于4腿梯子时,观测到了高温超导中常见的电荷密度条纹,此时保留状态数达到10⁴,使用的GPU显存小于12 GB.     

密度矩阵重正化群的异构并行优化

密度矩阵重正化群方法（DMRG）在求解一维强关联格点模型的基态时可以获得较高的精度，在应用于二维或准二维问题时，要达到类似的精度通常需要较大的计算量与存储空间.本文提出一种新的DMRG异构并行策略，可以同时发挥计算机中央处理器（CPU）和图形处理器（GPU）的计算性能.针对最耗时的哈密顿量对角化部分，实现了数据的分布式存储，并且给出了CPU和GPU之间的负载平衡策略.以费米Hubbard模型为例，测试了异构并行程序在不同DMRG保留状态数下的运行表现，并给出了相应的性能基准.应用于4腿梯子时，观测到了高温超导中常见的电荷密度条纹，此时保留状态数达到10⁴，使用的GPU显存小于12 GB.     

密度矩阵重正化群的异构并行优化

密度矩阵重正化群方法（DMRG）在求解一维强关联格点模型的基态时可以获得较高的精度，在应用于二维或准二维问题时，要达到类似的精度通常需要较大的计算量与存储空间.本文提出一种新的DMRG异构并行策略，可以同时发挥计算机中央处理器（CPU）和图形处理器（GPU）的计算性能.针对最耗时的哈密顿量对角化部分，实现了数据的分布式存储，并且给出了CPU和GPU之间的负载平衡策略.以费米Hubbard模型为例，测试了异构并行程序在不同DMRG保留状态数下的运行表现，并给出了相应的性能基准.应用于4腿梯子时，观测到了高温超导中常见的电荷密度条纹，此时保留状态数达到10⁴，使用的GPU显存小于12 GB.     

密度矩阵重正化群的异构并行优化

密度矩阵重正化群方法（DMRG）在求解一维强关联格点模型的基态时可以获得较高的精度,在应用于二维或准二维问题时,要达到类似的精度通常需要较大的计算量与存储空间.本文提出一种新的DMRG异构并行策略,可以同时发挥计算机中央处理器（CPU）和图形处理器（GPU）的计算性能.针对最耗时的哈密顿量对角化部分,实现了数据的分布式存储,并且给出了CPU和GPU之间的负载平衡策略.以费米Hubbard模型为例,测试了异构并行程序在不同DMRG保留状态数下的运行表现,并给出了相应的性能基准.应用于4腿梯子时,观测到了高温超导中常见的电荷密度条纹,此时保留状态数达到10⁴,使用的GPU显存小于12 GB.     

密度矩阵重正化群的异构并行优化

密度矩阵重正化群方法（DMRG）在求解一维强关联格点模型的基态时可以获得较高的精度,在应用于二维或准二维问题时,要达到类似的精度通常需要较大的计算量与存储空间.本文提出一种新的DMRG异构并行策略,可以同时发挥计算机中央处理器（CPU）和图形处理器（GPU）的计算性能.针对最耗时的哈密顿量对角化部分,实现了数据的分布式存储,并且给出了CPU和GPU之间的负载平衡策略.以费米Hubbard模型为例,测试了异构并行程序在不同DMRG保留状态数下的运行表现,并给出了相应的性能基准.应用于4腿梯子时,观测到了高温超导中常见的电荷密度条纹,此时保留状态数达到10⁴,使用的GPU显存小于12 GB.