HΦ [aitch-phi] is an open-source software package of numerically exact and stochastic calculations for a wide range of quantum many-body systems. In this paper, we present the newly added functions and the implemented methods in vers. 2 and 3. In ver. 2, we implement spectrum calculations by the shifted Krylov method, and low-energy excited state calculations by the locally optimal blocking preconditioned conjugate gradient (LOBPCG) method. In ver. 3, we implement the full diagonalization method using ScaLAPACK and GPGPU computing via MAGMA. We also implement a real-time evolution method and the canonical thermal pure quantum (cTPQ) state method for finite-temperature calculations. The Wannier90 format for specifying the Hamiltonians is also implemented. Using the Wannier90 format, it is possible to perform the calculations for the ab initio low-energy effective Hamiltonians of solids obtained by the open-source software RESPACK. We also update Standard mode—simplified input format in HΦ—to use these functions and methods. We explain the basics of the implemented methods and how to use them. Program Title: HΦ [aitch-phi] CPC Library link to program files: https://doi.org/10.17632/vnfthtyctm.2 Developer's repository link: https://github.com/issp-center-dev/HPhi Licensing provisions: GPLv3 Programming language: C, Fortran Journal reference of previous version: Comput. Phys. Commun. 217 (2017) 180–192 Does the new version supersede the previous version?: Yes. The latest version has compatibility with the old versions. Although most functions are available in newer versions, some redundant/unnecessary functions were abolished in newer versions. Reasons for the new version: Implementation of new functions and methods, development of utilities, and bug fixes. Summary of revisions: We added new functions to obtain excited spectrum, low-energy excited states, and time-dependent physical quantities. We also implemented the full diagonalization method by using MAGMA on GPGPUs and finite-temperature simulations by using the canonical thermal pure quantum state method. In addition, we developed utilities for connection with RESPACK and a submodule for a generator of input files. Nature of problem: Physical properties in quantum lattice models with finite system sizes Solution method: In HΦ, we implemented several numerical methods such as the Lanczos method, the full diagonalization method, the locally optimal block preconditioned conjugate gradient (LOBPCG) method, the real-time evolution method based on the Taylor expansion, the shifted Krylov method, and the microcanonical/canonical thermal pure quantum state method.
HΦ [aitch-phi] is an open-source software package of numerically exact and stochastic calculations for a wide range of quantum many-body systems. In this paper, we present the newly added functions and the implemented methods in vers. 2 and 3. In ver. 2, we implement spectrum calculations by the shifted Krylov method, and low-energy excited state calculations by the locally optimal blocking preconditioned conjugate gradient (LOBPCG) method. In ver. 3, we implement the full diagonalization method using ScaLAPACK and GPGPU computing via MAGMA. We also implement a real-time evolution method and the canonical thermal pure quantum (cTPQ) state method for finite-temperature calculations. The Wannier90 format for specifying the Hamiltonians is also implemented. Using the Wannier90 format, it is possible to perform the calculations for the ab initio low-energy effective Hamiltonians of solids obtained by the open-source software RESPACK. We also update Standard mode—simplified input format in HΦ—to use these functions and methods. We explain the basics of the implemented methods and how to use them. Program Title: HΦ [aitch-phi] CPC Library link to program files: https://doi.org/10.17632/vnfthtyctm.2 Developer's repository link: https://github.com/issp-center-dev/HPhi Licensing provisions: GPLv3 Programming language: C, Fortran Journal reference of previous version: Comput. Phys. Commun. 217 (2017) 180–192 Does the new version supersede the previous version?: Yes. The latest version has compatibility with the old versions. Although most functions are available in newer versions, some redundant/unnecessary functions were abolished in newer versions. Reasons for the new version: Implementation of new functions and methods, development of utilities, and bug fixes. Summary of revisions: We added new functions to obtain excited spectrum, low-energy excited states, and time-dependent physical quantities. We also implemented the full diagonalization method by using MAGMA on GPGPUs and finite-temperature simulations by using the canonical thermal pure quantum state method. In addition, we developed utilities for connection with RESPACK and a submodule for a generator of input files. Nature of problem: Physical properties in quantum lattice models with finite system sizes Solution method: In HΦ, we implemented several numerical methods such as the Lanczos method, the full diagonalization method, the locally optimal block preconditioned conjugate gradient (LOBPCG) method, the real-time evolution method based on the Taylor expansion, the shifted Krylov method, and the microcanonical/canonical thermal pure quantum state method.
H Φ [ aitch-phi ] is a program package based on the Lanczos-type eigenvalue solution applicable to a broad range of quantum lattice models, i.e., arbitrary quantum lattice models with two-body interactions, including the Heisenberg model, the Kitaev model, the Hubbard model and the Kondo-lattice model. While it works well on PCs and PC-clusters, H Φ also runs efficiently on massively parallel computers, which considerably extends the tractable range of the system size. In addition, unlike most existing packages, H Φ supports finite-temperature calculations through the method of thermal pure quantum (TPQ) states. In this paper, we explain theoretical background and user-interface of H Φ . We also show the benchmark results of H Φ on supercomputers such as the K computer at RIKEN Advanced Institute for Computational Science (AICS) and SGI ICE XA (Sekirei) at the Institute for the Solid State Physics (ISSP). Program summary Program Title: H Φ Program Files doi: http://dx.doi.org/10.17632/vnfthtyctm.1 Licensing provisions: GNU General Public License, version 3 or later Programming language: C External routines/libraries: MPI, BLAS, LAPACK Nature of problem: Physical properties (such as the magnetic moment, the specific heat, the spin susceptibility) of strongly correlated electrons at zero- and finite-temperature. Solution method: Application software based on the full diagonalization method, the exact diagonalization using the Lanczos method, and the microcanonical thermal pure quantum state for quantum lattice model such as the Hubbard model, the Heisenberg model and the Kondo model. Restrictions: System size less than about 20 sites for an itinerant electronic system and 40 site for a local spin system. Unusual features: Finite-temperature calculation of the strongly correlated electronic system based on the iterative scheme to construct the thermal pure quantum state. This method is efficient for highly frustrated system which is difficult to treat with other methods such as the unbiased quantum Monte Carlo.
mVMC (many-variable Variational Monte Carlo) is an open-source software package based on the variational Monte Carlo method applicable for a wide range of Hamiltonians for interacting fermion systems. In mVMC, we introduce more than ten thousands variational parameters and simultaneously optimize them by using the stochastic reconfiguration (SR) method. In this paper, we explain basics and user interfaces of mVMC. By using mVMC, users can perform the calculation by preparing only one input file of about ten lines for widely studied quantum lattice models, and can also perform it for general Hamiltonians by preparing several additional input files. We show the benchmark results of mVMC for the Hubbard model, the Heisenberg model, and the Kondo-lattice model. These benchmark results demonstrate that mVMC provides ground-state and low-energy-excited-state wave functions for interacting fermion systems with high accuracy. Program title: mVMC Program Files doi: http://dx.doi.org/10.17632/xhgyp6ncvt.1 Licensing provisions: GNU General Public License version 3 Programming language: C External routines/libraries: MPI, BLAS, LAPACK, Pfapack, ScaLAPACK (optional) Nature of problem: Physical properties (such as the charge/spin structure factors) of strongly correlated electrons at zero temperature. Solution method: Application software based on the variational Monte Carlo method for quantum lattice model such as the Hubbard model, the Heisenberg model and the Kondo model. Unusual features: It is possible to perform the highly-accurate calculations for ground states in a wide range of theoretical Hamiltonians in quantum many-body systems. In addition to the conventional orders such as magnetic and/or charge orders, user can treat the anisotropic superconductivities within the same framework. This flexibility is the main advantage of mVMC.
Quantum many-body systems are challenging targets for computational physics due to their large number of degrees of freedom. The tensor networks, particularly Tensor Product States (TPS) and Projected Entangled Pair States (PEPS), effectively represent these systems on two-dimensional lattices. However, the technical complexity of TPS/PEPS-based coding can be challenging for many researchers to manage effectively. To reduce this problem, we developed TeNeS (Tensor Network Solver). This paper introduces TeNeS -v2, which extends TeNeS with real-time and finite temperature simulations, providing deeper insights into quantum many-body systems. We detail the new algorithms, input/output design, and application examples, demonstrating TeNeS -v2's applicability to various quantum spin and Bose models on two-dimensional lattices. New version program summary Program Title: TeNeS CPC Library link to program files: https://doi.org/10.17632/psm26xxbvd.2 Developer's repository link:: https://github.com/issp-center-dev/TeNeS Licensing provisions: GPLv3 Programming language: C++11 Journal reference of previous version: Comput. Phys. Commun. (2022) 279, 108437, doi:10.1016/j.cpc.2022.108437 Reasons for the new version: To extend the capabilities of TeNeS to include real-time and finite temperature simulations, enabling deeper insights into quantum many-body systems beyond ground states. Summary of revisions: TeNeS -v2 introduces real-time evolution and finite temperature simulation capabilities, expanding the range of quantum many-body phenomena that can be studied. Nature of problem: Quantum many-body systems are extremely difficult to simulate because of the huge dimension of the Hilbert space. Conventional methods have difficulty in accurately representing these systems, especially for properties beyond small lattices and ground states. Advanced computational techniques are required to efficiently capture the dynamics and thermal properties of these systems. Solution method: TeNeS employs tensor networks, specifically Tensor Product States (TPS) and Projected Entangled Pair States (PEPS), to efficiently represent quantum many-body states on two-dimensional lattices. The real-time evolution is handled through a straightforward extension of imaginary time evolution methods, allowing the study of dynamical properties. Finite temperature simulations are conducted using an imaginary time evolution approach starting from the infinite-temperature mixed state. These methods enable the accurate calculation of various physical properties in different quantum states. Additional comments including restrictions and unusual features: TeNeS -v2 maintains ease of use by providing flexible input file configurations and supporting a variety of two-dimensional lattice models. Increasing the bond dimension to improve accuracy requires more computational resources. Finite temperature calculations may exhibit unphysical behavior due to limitations of the tensor network approximation method. Users should be aware of these potential problems and interpret the results following the suggestions made in the text.
We present an approach to the calculation of arbitrary spectral, thermal, and excited state properties within the full configuration interaction quzantum Monte Carlo framework. This is achieved via an unbiased projection of the Hamiltonian eigenvalue problem into a space of stochastically sampled Krylov vectors, thus, enabling the calculation of real-frequency spectral and thermal properties and avoiding explicit analytic continuation. We use this approach to calculate temperature-dependent properties and one- and two-body spectral functions for various Hubbard models, as well as isolated excited states in ab initio systems.
We present release 3.0 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an open-source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. As in previous releases, development is centered on common data formats, on libraries to simplify and speed up code development, and on full-featured simulation programs that let non-experts carry out serial or parallel numerical simulations using the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended-ensemble simulations, exact and full diagonalization (ED), the density matrix renormalization group (DMRG), and continuous-time QMC solvers for dynamical mean-field theory (DMFT). Major changes in release 3.0 include distribution of the pyalps binary through the Python Package Index (pip install pyalps) and through Spack for HPC systems; migration of development to GitHub with continuous integration and automated testing; relicensing of the package under the permissive MIT license; a rebuilt documentation and tutorial website, including a set of Jupyter-notebook tutorials and localized content; a broad modernization of the C++ codebase (C++17 compliance, Boost and NumPy 2.0 compatibility, and warning and dead-code cleanup) together with a major DMRG update and associated reliability and build-compatibility fixes; the removal of legacy components (the VisTrails provenance integration and the TEBD, MPS, and directed-worm-algorithm application codes); archival of release 3.0.0 with a Zenodo DOI; and a formal governance and sustainability model developed under the US National Science Foundation (NSF) POSE program. The software is available at https://alps.comp-phys.org.
We present a general-purpose, open-source software suite, SQUIRREL (Streamlined Quantum Unified Interface for Researching Real-time Excitations with Light), for propagating the time-dependent Schrödinger equation on complex geometries in the presence of time-dependent electric and/or magnetic fields. To handle large systems that can be executed on a conventional desktop computer, the SQUIRREL software suite uses a suite of efficient propagation methods for various quantum dynamics applications, including a new perturbation-based element-dropping algorithm that improves computational performance with minimal loss of accuracy. We analyze the efficacy of these optimizations for Crank-Nicolson, scaled Taylor series approximation, and split-operator propagation methods and discuss the range of their applicability to a variety of quantum dynamics problems. In addition, we provide several examples of time-dependent dynamics calculations and extensive documentation for generating custom geometries, potentials, and time-propagation approaches. Our numerical benchmarks and results demonstrate the versatility of the SQUIRREL software suite for efficiently calculating quantum dynamics in complex nanoscale geometries, particularly in the presence of time-dependent magnetic fields, which have received less attention in previous quantum dynamics studies. Program summary Program Title: SQUIRREL CPC Library link to program files: https://doi.org/10.17632/kfvs5s88sj.1 Licensing provisions: GNU General Public License 3 Programming language: MATLAB Supplementary material: Additional details on calculations with an anisotropic effective mass, determination of timesteps, propagation timings, and element-dropping timings. Nature of problem: The SQUIRREL software suite solves the time-dependent Schrödinger equation for electronic states in the presence of time-dependent electric and/or magnetic fields. The code is well-suited for calculating electron dynamics of complex nanostructures in the effective mass approximation, using various propagation and sparse-matrix techniques to reduce the computational complexity of these problems. The program provides a user-friendly interface for testing these optimizations and generating custom geometries and potentials for systems of interest. Solution method: The SQUIRREL software suite uses Crank-Nicolson, split-operator, Padé approximant, and scaled-Taylor-series-based propagation methods to calculate electron dynamics on a finite-element basis. The code includes a novel perturbation-based element-dropping method, which enables sparse representations of the typically dense finite element Hamiltonian matrix. This capability enables a highly efficient and flexible approach for calculating driven dynamics in quantum systems with complex geometries and time-dependent electric/magnetic fields.