-
Robust analytic continuation using sparse modeling approach imposed by semi-positive definiteness for multi-orbital systems
Authors:
Yuichi Motoyama,
Hiroshi Shinaoka,
Junya Otsuki,
Kazuyoshi Yoshimi
Abstract:
Analytic continuation (AC) from imaginary-time Green's function to spectral function is essential in the numerical analysis of dynamical properties in quantum many-body systems. However, this process faces a fundamental challenge: it is an ill-posed problem, leading to unstable spectra against the noise in the Green's function. This instability is further complicated in multi-orbital systems with…
▽ More
Analytic continuation (AC) from imaginary-time Green's function to spectral function is essential in the numerical analysis of dynamical properties in quantum many-body systems. However, this process faces a fundamental challenge: it is an ill-posed problem, leading to unstable spectra against the noise in the Green's function. This instability is further complicated in multi-orbital systems with hybridization between spin-orbitals, where off-diagonal Green's functions yield a spectral matrix with off-diagonal elements, necessitating the matrix's semi-positive definiteness to satisfy the causality. We propose an advanced AC method using sparse modeling for multi-orbital systems, which reduces the effect of noise and ensures the matrix's semi-positive definiteness. We demonstrate the effectiveness of this approach by contrasting it with the conventional sparse modeling method, focusing on handling Green's functions with off-diagonal elements, thereby demonstrating our proposed method's enhanced stability and precision.
△ Less
Submitted 2 September, 2024;
originally announced September 2024.
-
Configuration sampling in multi-component multi-sublattice systems enabled by ab Initio Configuration Sampling Toolkit (abICS)
Authors:
Shusuke Kasamatsu,
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Tatsumi Aoyama
Abstract:
Simulation of the intermediate levels of disorder found in multi-component multi-sublattice systems in various functional materials is a challenging issue, even for state-of-the-art methodologies based on first-principles calculation. Here, we introduce our open-source package ab Initio Configuration Sampling Toolkit (abICS), which combines high-throughput first-principles calculations, machine le…
▽ More
Simulation of the intermediate levels of disorder found in multi-component multi-sublattice systems in various functional materials is a challenging issue, even for state-of-the-art methodologies based on first-principles calculation. Here, we introduce our open-source package ab Initio Configuration Sampling Toolkit (abICS), which combines high-throughput first-principles calculations, machine learning, and parallel extended ensemble sampling in an active learning setting to enable such simulations. The theoretical background is reviewed in some detail followed by brief notes on usage of the software. In addition, our recent applications of abICS to multi-component ionic systems and their interfaces for energy applications are reviewed as demonstration of the power of this approach.
△ Less
Submitted 9 September, 2023;
originally announced September 2023.
-
H-wave -- A Python package for the Hartree-Fock approximation and the random phase approximation
Authors:
Tatsumi Aoyama,
Kazuyoshi Yoshimi,
Kota Ido,
Yuichi Motoyama,
Taiki Kawamura,
Takahiro Misawa,
Takeo Kato,
Akito Kobayashi
Abstract:
H-wave is an open-source software package for performing the Hartree--Fock approximation (HFA) and random phase approximation (RPA) for a wide range of Hamiltonians of interacting fermionic systems. In HFA calculations, H-wave examines the stability of several symmetry-broken phases, such as anti-ferromagnetic and charge-ordered phases, in the given Hamiltonians at zero and finite temperatures. Fu…
▽ More
H-wave is an open-source software package for performing the Hartree--Fock approximation (HFA) and random phase approximation (RPA) for a wide range of Hamiltonians of interacting fermionic systems. In HFA calculations, H-wave examines the stability of several symmetry-broken phases, such as anti-ferromagnetic and charge-ordered phases, in the given Hamiltonians at zero and finite temperatures. Furthermore, H-wave calculates the dynamical susceptibilities using RPA to examine the instability toward the symmetry-broken phases. By preparing a simple input file for specifying the Hamiltonians, users can perform HFA and RPA for standard Hamiltonians in condensed matter physics, such as the Hubbard model and its extensions. Additionally, users can use a Wannier90-like format to specify fermionic Hamiltonians. A Wannier90 format is implemented in RESPACK to derive ab initio Hamiltonians for solids. HFA and RPA for the ab initio Hamiltonians can be easily performed using H-wave. In this paper, we first explain the basis of HFA and RPA, and the basic usage of H-wave, including download and installation. Thereafter, the input file formats implemented in H-wave, including the Wannier90-like format for specifying the interacting fermionic Hamiltonians, are discussed. Finally, we present several examples of H-wave such as zero-temperature HFA calculations for the extended Hubbard model on a square lattice, finite-temperature HFA calculations for the Hubbard model on a cubic lattice, and RPA in the extended Hubbard model on a square lattice.
△ Less
Submitted 20 January, 2024; v1 submitted 1 August, 2023;
originally announced August 2023.
-
Update of $\mathcal{H}Φ$: Newly added functions and methods in versions 2 and 3
Authors:
Kota Ido,
Mitsuaki Kawamura,
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Youhei Yamaji,
Synge Todo,
Naoki Kawashima,
Takahiro Misawa
Abstract:
$\mathcal{H}Φ$ [$aitch$-$phi…
▽ More
$\mathcal{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 $\unicode{x2014}$simplified input format in $\mathcal{H}Φ$$\unicode{x2014}$ to use these functions and methods. We explain the basics of the implemented methods and how to use them.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
MateriApps LIVE! and MateriApps Installer: Environment for starting and scaling up materials science simulations
Authors:
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Takeo Kato,
Synge Todo
Abstract:
In our current era, numerical simulations have become indispensable theoretical and experimental tools for use in daily research activities, particularly in the materials science fields. However, the installation processes for such simulations frequently become problematic because they depend strongly on the device environment, and troubleshooting those processes is a challenging task for beginner…
▽ More
In our current era, numerical simulations have become indispensable theoretical and experimental tools for use in daily research activities, particularly in the materials science fields. However, the installation processes for such simulations frequently become problematic because they depend strongly on the device environment, and troubleshooting those processes is a challenging task for beginners. To minimize such difficulties, we created MateriApps LIVE! and MateriApps Installer, which can solve most of the related issues. Specifically, MateriApps LIVE! offers a virtual environment in which users can quickly try out computational materials science simulations on a personal computer while MateriApps Installer provides a comprehensive set of shell scripts for use when installing software on Unix, Linux, macOS, and supercomputer systems. Herein, we provide detailed descriptions of MateriApps LIVE! and MateriApps Installer together with illustrative examples of their use.
△ Less
Submitted 4 May, 2022;
originally announced May 2022.
-
Reply to Comments on "Universal and Non-Universal Correction Terms of Bose Gases in Dilute Region: A Quantum Monte Carlo Study"
Authors:
Akiko Masaki-Kato,
Yuichi Motoyama,
Naoki Kawashima
Abstract:
We reply to the comments on our previous paper J. Phys. Soc. Jpn. 91, 024001 (2022)
We reply to the comments on our previous paper J. Phys. Soc. Jpn. 91, 024001 (2022)
△ Less
Submitted 21 April, 2022;
originally announced April 2022.
-
Data-analysis software framework 2DMAT and its application to experimental measurements for two-dimensional material structures
Authors:
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Izumi Mochizuki,
Harumichi Iwamoto,
Hayato Ichinose,
Takeo Hoshi
Abstract:
An open-source data-analysis framework 2DMAT has been developed for experimental measurements of two-dimensional material structures. 2DMAT offers five analysis methods: (i) Nelder-Mead optimization, (ii) grid search, (iii) Bayesian optimization, (iv) replica exchange Monte Carlo method, and (v) population-annealing Monte Carlo method. Methods (ii) through (v) are implemented by parallel computati…
▽ More
An open-source data-analysis framework 2DMAT has been developed for experimental measurements of two-dimensional material structures. 2DMAT offers five analysis methods: (i) Nelder-Mead optimization, (ii) grid search, (iii) Bayesian optimization, (iv) replica exchange Monte Carlo method, and (v) population-annealing Monte Carlo method. Methods (ii) through (v) are implemented by parallel computation,which is efficient not only for personal computers but also for supercomputers.The current version of 2DMAT is applicable to total-reflection high-energy positron diffraction (TRHEPD), surface X-ray diffraction (SXRD), and low-energy electron diffraction (LEED) experiments by installing corresponding forward problem solvers that generate diffraction intensity data from a given dataset of the atomic positions. The analysis methods are general and can be applied also to other experiments and problems.
△ Less
Submitted 12 July, 2022; v1 submitted 9 April, 2022;
originally announced April 2022.
-
TeNeS: Tensor Network Solver for Quantum Lattice Systems
Authors:
Yuichi Motoyama,
Tsuyoshi Okubo,
Kazuyoshi Yoshimi,
Satoshi Morita,
Takeo Kato,
Naoki Kawashima
Abstract:
TeNeS (Tensor Network Solver) is a free/libre open-source software program package for calculating two-dimensional many-body quantum states based on the tensor network method and the corner transfer matrix renormalization group (CTMRG) method. This package calculates ground-state wavefunctions for user-defined Hamiltonians and evaluates user-defined physical quantities such as magnetization and co…
▽ More
TeNeS (Tensor Network Solver) is a free/libre open-source software program package for calculating two-dimensional many-body quantum states based on the tensor network method and the corner transfer matrix renormalization group (CTMRG) method. This package calculates ground-state wavefunctions for user-defined Hamiltonians and evaluates user-defined physical quantities such as magnetization and correlation functions. For certain predefined models and lattices, there is a tool that makes it easy to generate input files. TeNeS uses an OpenMP/MPI hybrid parallelized tensor operation library and thus can perform large-scale calculations using massively parallel machines.
△ Less
Submitted 16 June, 2022; v1 submitted 24 December, 2021;
originally announced December 2021.
-
Universal and Non-Universal Correction Terms of Bose Gases in Dilute Region: a Quantum Monte Carlo Study
Authors:
Akiko Masaki-Kato,
Yuichi Motoyama,
Naoki Kawashima
Abstract:
We study dilute gases of interacting bosons at zero-temperature in the region where the system is characterized only by the s-wave scattering length. We carry out quantum Monte Carlo simulation of the Bose Hubbard model and a continuous-space hard-core model. Fitting the extended Lee-Huang-Yang formula to the Monte Carlo results establishes the detailed mapping from the lattice model to the contin…
▽ More
We study dilute gases of interacting bosons at zero-temperature in the region where the system is characterized only by the s-wave scattering length. We carry out quantum Monte Carlo simulation of the Bose Hubbard model and a continuous-space hard-core model. Fitting the extended Lee-Huang-Yang formula to the Monte Carlo results establishes the detailed mapping from the lattice model to the continuous field theory characterized only by the s-wave scattering length. Our estimate of the intrinsic s-wave scattering length $a_s$ of the Bose-Hubbard model is $a_s/a_l = 0.316(2)$ where $a_l$ is the lattice constant. It turned out that inclusion of the universal second correction term of $O(n \log n)$ makes the fitting worse while it is restored by inclusion of the non-universal third correction term, of which the existence was predicted analytically.
△ Less
Submitted 19 December, 2021;
originally announced December 2021.
-
sim-trhepd-rheed -- Open-source simulator of total-reflection high-energy positron diffraction (TRHEPD) and reflection high-energy electron diffraction (RHEED)
Authors:
Takashi Hanada,
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Takeo Hoshi
Abstract:
The present paper reports sim-trhepd-rheed (STR), an open-source simulator of total-reflection high-energy positron diffraction (TRHEPD) and reflection high-energy electron diffraction (RHEED) experiments which are used for atom-scale surface structure determination of a material. The STR simulator is used for the analysis of experimental diffraction data by simulating the rocking curve from a giv…
▽ More
The present paper reports sim-trhepd-rheed (STR), an open-source simulator of total-reflection high-energy positron diffraction (TRHEPD) and reflection high-energy electron diffraction (RHEED) experiments which are used for atom-scale surface structure determination of a material. The STR simulator is used for the analysis of experimental diffraction data by simulating the rocking curve from a given trial surface structure by solving the partial differential equation of the dynamical quantum diffraction theory for positron or electron wavefunctions. Using the obtained surface structure, electronic structure, and other physical quantities can be evaluated through first-principles calculations. For this purpose, a utility software was also developed in order to realize a first principles calculation with the Quantum ESPRESSO suite.
△ Less
Submitted 24 April, 2022; v1 submitted 18 October, 2021;
originally announced October 2021.
-
Bayesian optimization package: PHYSBO
Authors:
Yuichi Motoyama,
Ryo Tamura,
Kazuyoshi Yoshimi,
Kei Terayama,
Tsuyoshi Ueno,
Koji Tsuda
Abstract:
PHYSBO (optimization tools for PHYSics based on Bayesian Optimization) is a Python library for fast and scalable Bayesian optimization. It has been developed mainly for application in the basic sciences such as physics and materials science. Bayesian optimization is used to select an appropriate input for experiments/simulations from candidate inputs listed in advance in order to obtain better out…
▽ More
PHYSBO (optimization tools for PHYSics based on Bayesian Optimization) is a Python library for fast and scalable Bayesian optimization. It has been developed mainly for application in the basic sciences such as physics and materials science. Bayesian optimization is used to select an appropriate input for experiments/simulations from candidate inputs listed in advance in order to obtain better output values with the help of machine learning prediction. PHYSBO can be used to find better solutions for both single and multi-objective optimization problems. At each cycle in the Bayesian optimization, a single proposal or multiple proposals can be obtained for the next experiments/simulations. These proposals can be obtained interactively for use in experiments. PHYSBO is available at https://github.com/issp-center-dev/PHYSBO.
△ Less
Submitted 24 May, 2022; v1 submitted 15 October, 2021;
originally announced October 2021.
-
Data analysis of $ab$ $initio$ effective Hamiltonians in iron-based superconductors $\unicode{x2014}$ Construction of predictors for superconducting critical temperature
Authors:
Kota Ido,
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Takahiro Misawa
Abstract:
High-temperature superconductivity occurs in strongly correlated materials such as copper oxides and iron-based superconductors. Numerous experimental and theoretical works have been done to identify the key parameters that induce high-temperature superconductivity. However, the key parameters governing the high-temperature superconductivity remain still unclear, which hamper the prediction of sup…
▽ More
High-temperature superconductivity occurs in strongly correlated materials such as copper oxides and iron-based superconductors. Numerous experimental and theoretical works have been done to identify the key parameters that induce high-temperature superconductivity. However, the key parameters governing the high-temperature superconductivity remain still unclear, which hamper the prediction of superconducting critical temperatures ($T_\text{c}$s) of strongly correlated materials. Here by using data-science techniques, we clarified how the microscopic parameters in the $ab$ $initio$ effective Hamiltonians correlate with the experimental $T_\text{c}$s in iron-based superconductors. We showed that a combination of microscopic parameters can characterize the compound-dependence of $T_\text{c}$ using the principal component analysis. We also constructed a linear regression model that reproduces the experimental $T_\text{c}$ from the microscopic parameters. Based on the regression model, we showed a way for increasing $T_\text{c}$ by changing the lattice parameters. The developed methodology opens a new field of materials informatics for strongly correlated electron systems.
△ Less
Submitted 7 February, 2023; v1 submitted 19 September, 2021;
originally announced September 2021.
-
Robust analytic continuation combining the advantages of the sparse modeling approach and Padé approximation
Authors:
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Junya Otsuki
Abstract:
Analytic continuation (AC) from the imaginary-time Green's function to the spectral function is a crucial process for numerical studies of the dynamical properties of quantum many-body systems. This process, however, is an ill-posed problem; that is, the obtained spectrum is unstable against the noise of the Green's function. Though several numerical methods have been developed, each of them has i…
▽ More
Analytic continuation (AC) from the imaginary-time Green's function to the spectral function is a crucial process for numerical studies of the dynamical properties of quantum many-body systems. This process, however, is an ill-posed problem; that is, the obtained spectrum is unstable against the noise of the Green's function. Though several numerical methods have been developed, each of them has its own advantages and disadvantages. The sparse modeling (SpM) AC method, for example, is robust against the noise of the Green's function but suffers from unphysical oscillations in the low-energy region. We propose a new method that combines the SpM AC with the Padé approximation. This combination, called SpM-Padé, inherits robustness against noise from SpM and low-energy accuracy from Padé, compensating for the disadvantages of each. We demonstrate that the SpM- Padé method yields low-variance and low-biased results with almost the same computational cost as that of the SpM method.
△ Less
Submitted 25 January, 2022; v1 submitted 17 September, 2021;
originally announced September 2021.
-
Non-monotonic behavior of the Binder Parameter in the discrete spin systems
Authors:
Hiroshi Watanabe,
Yuichi Motoyama,
Satoshi Morita,
Naoki Kawashima
Abstract:
We study a non-monotonic behavior of the Binder parameter, which appears in the discrete spin systems. We show that the Binder parameters of the Potts model are non-monotonic for $q=3$ and $4$, while they are monotonic for the Ising case ($q=2$). Using the Fortuin-Kasteleyn graph representation, we find that the improved estimator of the Binder parameter consists of two terms with values only in h…
▽ More
We study a non-monotonic behavior of the Binder parameter, which appears in the discrete spin systems. We show that the Binder parameters of the Potts model are non-monotonic for $q=3$ and $4$, while they are monotonic for the Ising case ($q=2$). Using the Fortuin-Kasteleyn graph representation, we find that the improved estimator of the Binder parameter consists of two terms with values only in high- and low-temperature regions. The non-monotonic behavior is found to originate from the low-temperature term. With the appropriately defined order parameter, we can reduce the influence of the low-temperature term, and as a result, the non-monotonic behavior can also be reduced. We propose new definitions of the order parameter, which reduces or eliminates the non-monotonic behavior of the Binder parameter in a system for which the improved estimator of the Binder parameter is unknown.
△ Less
Submitted 3 February, 2023; v1 submitted 25 August, 2021;
originally announced August 2021.
-
Facilitating {\it ab initio} configurational sampling of multicomponent solids using an on-lattice neural network model and active learning
Authors:
Shusuke Kasamatsu,
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Ushio Matsumoto,
Akihide Kuwabara,
Takafumi Ogawa
Abstract:
We propose a scheme for {\it ab initio} configurational sampling in multicomponent crystalline solids using Behler-Parinello type neural network potentials (NNPs) in an unconventional way: the NNPs are trained to predict the energies of relaxed structures from the perfect lattice with configurational disorder instead of the usual way of training to predict energies as functions of continuous atom…
▽ More
We propose a scheme for {\it ab initio} configurational sampling in multicomponent crystalline solids using Behler-Parinello type neural network potentials (NNPs) in an unconventional way: the NNPs are trained to predict the energies of relaxed structures from the perfect lattice with configurational disorder instead of the usual way of training to predict energies as functions of continuous atom coordinates. An active learning scheme is employed to obtain a training set containing configurations of thermodynamic relevance. This enables bypassing of the structural relaxation procedure which is necessary when applying conventional NNP approaches to the lattice configuration problem. The idea is demonstrated on the calculation of the temperature dependence of the degree of A/B site inversion in three spinel oxides, MgAl$_2$O$_4$, ZnAl$_2$O$_4$, and MgGa$_2$O$_4$. The present scheme may serve as an alternative to cluster expansion for `difficult' systems, e.g., complex bulk or interface systems with many components and sublattices that are relevant to many technological applications today.
△ Less
Submitted 20 April, 2022; v1 submitted 6 August, 2020;
originally announced August 2020.
-
DSQSS: Discrete Space Quantum Systems Solver
Authors:
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Akiko Masaki-Kato,
Takeo Kato,
Naoki Kawashima
Abstract:
The Discrete Space Quantum Systems Solver (DSQSS) is a program package for solving quantum many-body problems defined on lattices. The DSQSS is based on the quantum Monte Carlo method in Feynman's path integral representation and covers a broad range of problems using flexible input files that define arbitrary unit cells in arbitrary dimensions and arbitrary matrix elements representing the intera…
▽ More
The Discrete Space Quantum Systems Solver (DSQSS) is a program package for solving quantum many-body problems defined on lattices. The DSQSS is based on the quantum Monte Carlo method in Feynman's path integral representation and covers a broad range of problems using flexible input files that define arbitrary unit cells in arbitrary dimensions and arbitrary matrix elements representing the interactions among an arbitrary number of degrees of freedom. Finite temperature calculations of quantum spin and the Bose- Hubbard models can be performed by specifying parameters such as the number of dimensions, the lattice size, coupling constants, and temperature. The present paper details the use of DSQSS and presents a number of applications thereof.
△ Less
Submitted 23 March, 2021; v1 submitted 22 July, 2020;
originally announced July 2020.
-
K$ω$ -- Open-source library for the shifted Krylov subspace method of the form $(zI-H)x=b$
Authors:
Takeo Hoshi,
Mitsuaki Kawamura,
Kazuyoshi Yoshimi,
Yuichi Motoyama,
Takahiro Misawa,
Youhei Yamaji,
Synge Todo,
Naoki Kawashima,
Tomohiro Sogabe
Abstract:
We develop K$ω$, an open-source linear algebra library for the shifted Krylov subspace methods. The methods solve a set of shifted linear equations $(z_k I-H)x^{(k)}=b\, (k=0,1,2,...)$ for a given matrix $H$ and a vector $b$, simultaneously. The leading order of the operational cost is the same as that for a single equation. The shift invariance of the Krylov subspace is the mathematical foundatio…
▽ More
We develop K$ω$, an open-source linear algebra library for the shifted Krylov subspace methods. The methods solve a set of shifted linear equations $(z_k I-H)x^{(k)}=b\, (k=0,1,2,...)$ for a given matrix $H$ and a vector $b$, simultaneously. The leading order of the operational cost is the same as that for a single equation. The shift invariance of the Krylov subspace is the mathematical foundation of the shifted Krylov subspace methods. Applications in materials science are presented to demonstrate the advantages of the algorithm over the standard Krylov subspace methods such as the Lanczos method. We introduce benchmark calculations of (i) an excited (optical) spectrum and (ii) intermediate eigenvalues by the contour integral on the complex plane. In combination with the quantum lattice solver $\mathcal{H} Φ$, K$ω$ can realize parallel computation of excitation spectra and intermediate eigenvalues for various quantum lattice models.
△ Less
Submitted 28 May, 2020; v1 submitted 23 January, 2020;
originally announced January 2020.
-
RESPACK: An ab initio tool for derivation of effective low-energy model of material
Authors:
Kazuma Nakamura,
Yoshihide Yoshimoto,
Yusuke Nomura,
Terumasa Tadano,
Mitsuaki Kawamura,
Taichi Kosugi,
Kazuyoshi Yoshimi,
Takahiro Misawa,
Yuichi Motoyama
Abstract:
RESPACK is a first-principles calculation software for evaluating the interaction parameters of materials and is able to calculate maximally localized Wannier functions, response functions based on the random phase approximation and related optical properties, and frequency-dependent electronic interaction parameters. RESPACK receives its input data from a band-calculation code using norm-conservi…
▽ More
RESPACK is a first-principles calculation software for evaluating the interaction parameters of materials and is able to calculate maximally localized Wannier functions, response functions based on the random phase approximation and related optical properties, and frequency-dependent electronic interaction parameters. RESPACK receives its input data from a band-calculation code using norm-conserving pseudopotentials with plane-wave basis sets. Automatic generation scripts that convert the band-structure results to the RESPACK inputs are prepared for xTAPP and Quantum ESPRESSO. An input file for specifying the RESPACK calculation conditions is designed pursuing simplicity and is given in the Fortran namelist format. RESPACK supports hybrid parallelization using OpenMP and MPI and can treat large systems including a few hundred atoms in the calculation cell.
△ Less
Submitted 7 January, 2020;
originally announced January 2020.
-
SpM: Sparse modeling tool for analytic continuation of imaginary-time Green's function
Authors:
Kazuyoshi Yoshimi,
Junya Otsuki,
Yuichi Motoyama,
Masayuki Ohzeki,
Hiroshi Shinaoka
Abstract:
We present SpM, a sparse modeling tool for the analytic continuation of imaginary-time Green's function, licensed under GNU General Public License version 3. In quantum Monte Carlo simulation, dynamic physical quantities such as single-particle and magnetic excitation spectra can be obtained by applying analytic continuation to imaginary-time data. However, analytic continuation is an ill-conditio…
▽ More
We present SpM, a sparse modeling tool for the analytic continuation of imaginary-time Green's function, licensed under GNU General Public License version 3. In quantum Monte Carlo simulation, dynamic physical quantities such as single-particle and magnetic excitation spectra can be obtained by applying analytic continuation to imaginary-time data. However, analytic continuation is an ill-conditioned inverse problem and thus sensitive to noise and statistical errors. SpM provides stable analytic continuation against noise by means of a modern regularization technique, which automatically selects bases that contain relevant information unaffected by noise. This paper details the use of this program and shows some applications.
△ Less
Submitted 5 April, 2019;
originally announced April 2019.
-
Asymmetric melting of one-third plateau in kagome quantum antiferromagnets
Authors:
Takahiro Misawa,
Yuichi Motoyama,
Youhei Yamaji
Abstract:
Asymmetric destruction of the one-third magnetization plateau upon heating is found in the spin-1/2 kagome Heisenberg antiferromagnets using the typical pure quantum state approach. The asymmetry originates from larger density of states of low-lying excited sates of $N_{\rm s}$ spin systems with magnetization $(1/3-2/N_{\rm s})$ than that of low-lying sates with magnetization $(1/3+2/N_{\rm s})$.…
▽ More
Asymmetric destruction of the one-third magnetization plateau upon heating is found in the spin-1/2 kagome Heisenberg antiferromagnets using the typical pure quantum state approach. The asymmetry originates from larger density of states of low-lying excited sates of $N_{\rm s}$ spin systems with magnetization $(1/3-2/N_{\rm s})$ than that of low-lying sates with magnetization $(1/3+2/N_{\rm s})$. The enhanced specific heat and entropy that reflect the larger density of states in the lower-field side of the plateau are detectable in candidate materials of the kagome antiferromagnets. We discuss that the asymmetry originates from the unprecedented preservation of the ice rule around the plateau.
△ Less
Submitted 9 September, 2020; v1 submitted 22 January, 2018;
originally announced January 2018.
-
mVMC - Open-source software for many-variable variational Monte Carlo method
Authors:
Takahiro Misawa,
Satoshi Morita,
Kazuyoshi Yoshimi,
Mitsuaki Kawamura,
Yuichi Motoyama,
Kota Ido,
Takahiro Ohgoe,
Masatoshi Imada,
Takeo Kato
Abstract:
mVMC (many-variable Variational Monte Carlo) is an open-source software 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…
▽ More
mVMC (many-variable Variational Monte Carlo) is an open-source software 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.
△ Less
Submitted 2 September, 2018; v1 submitted 30 November, 2017;
originally announced November 2017.
-
Dimer-Mott and charge-ordered insulating states in the quasi-one-dimensional organic conductors $δ'_{P}$- and $δ'_{C}$-(BPDT-TTF)$_2$ICl$_2$
Authors:
R. Kobayashi,
K. Hashimoto,
N. Yoneyama,
K. Yoshimi,
Y. Motoyama,
S. Iguchi,
Y. Ikemoto,
T. Moriwaki,
H. Taniguchi,
T. Sasaki
Abstract:
We investigated the electronic states of the quasi-one-dimensional organic conductors $δ'_{P}$-(BPDT-TTF)$_2$ICl$_2$ and $δ'_{C}$-(BPDT-TTF)$_2$ICl$_2$, both of which are insulating at room temperature owing to strong electron correlations. Through measurements of electrical resistivity, optical conductivity, and magnetic susceptibility, as well as band-structure calculations, we have revealed tha…
▽ More
We investigated the electronic states of the quasi-one-dimensional organic conductors $δ'_{P}$-(BPDT-TTF)$_2$ICl$_2$ and $δ'_{C}$-(BPDT-TTF)$_2$ICl$_2$, both of which are insulating at room temperature owing to strong electron correlations. Through measurements of electrical resistivity, optical conductivity, and magnetic susceptibility, as well as band-structure calculations, we have revealed that the two materials possess completely different ground states, even though they have the same chemical composition and stacking configuration of the donor molecules. We have found that the $δ_P'$-type salt with an effective half-filled band behaves as a dimer-Mott (DM) insulator and exhibits a nonmagnetic transition at 25 K, whereas the $δ'_C$-type salt with a 3/4-filled band shows a charge ordering (CO) transition just above room temperature and becomes nonmagnetic below 20 K. The optical spectra of the $δ_P'$-type salt are composed of two characteristic bands due to intra- and interdimer charge transfers, supporting the DM insulating behavior arising from the strong on-site Coulomb interaction. By contrast, in the $δ'_C$-type salt, a single band characterizing the formation of CO arising from the off-site Coulomb interactions is observed. Upon lowering the temperature, the shape of the optical spectra in the $δ_C'$-type salt becomes asymmetric and shifts to much lower frequencies, suggesting the emergence of domain-wall excitations with fractional charges expected in a one-dimensional CO chain. The temperature dependence of the magnetic susceptibility of the $δ_P'$-type salt is well described by a 2D spin-1/2 Heisenberg AFM model on an anisotropic square lattice in the dimerized picture, while in the $δ_C'$-type salt, it can be explained by a 2D spin-1/2 Heisenberg AFM model on an anisotropic honeycomb lattice formed in the CO state.
△ Less
Submitted 10 September, 2017; v1 submitted 30 May, 2017;
originally announced May 2017.
-
Z_N Berry phase and symmetry protected topological phases of SU(N) antiferromagnetic Heisenberg chain
Authors:
Yuichi Motoyama,
Synge Todo
Abstract:
The local Z_N quantized Berry phase for the SU(N) antiferromagnetic Heisenberg spin model is formulated. This quantity, which is a generalization of the local Z_2 Berry phase for SU(2) symmetry, has a direct correspondence to the number of singlet pairs spanning on a particular bond, and is effective as a tool to characterize and classify the various symmetry protected topological phases of one-di…
▽ More
The local Z_N quantized Berry phase for the SU(N) antiferromagnetic Heisenberg spin model is formulated. This quantity, which is a generalization of the local Z_2 Berry phase for SU(2) symmetry, has a direct correspondence to the number of singlet pairs spanning on a particular bond, and is effective as a tool to characterize and classify the various symmetry protected topological phases of one-dimensional SU(N) spin systems. We extend the path-integral quantum Monte Carlo method for the Z_2 Berry phase in order to calculate the ZN Berry phase numerically. We demonstrate our method by calculating the Z_4 Berry phase for the bond-alternating SU(4) antiferromagnetic Heisenberg chain, which is represented by the Young diagram of four columns.
△ Less
Submitted 4 August, 2015;
originally announced August 2015.
-
Path-integral Monte Carlo for the Gauge-fixed Berry Connection and the Local Z_2 Berry Phase
Authors:
Yuichi Motoyama,
Synge Todo
Abstract:
We present a general formula of the gauge-fixed Berry connection which can be evaluated by path- integral Monte Carlo method. We also propose that the gauge-fixed local Berry connection can be an- other effective tool to estimate precisely the quantum critical point. For a demonstration, we calculate the gauge-fixed Berry connection and the local Z2 Berry phase of the antiferromagnetic Heisenberg…
▽ More
We present a general formula of the gauge-fixed Berry connection which can be evaluated by path- integral Monte Carlo method. We also propose that the gauge-fixed local Berry connection can be an- other effective tool to estimate precisely the quantum critical point. For a demonstration, we calculate the gauge-fixed Berry connection and the local Z2 Berry phase of the antiferromagnetic Heisenberg model on a staggered bond-alternating ladder, and estimated quantum critical point is consistent with other methods.
△ Less
Submitted 1 October, 2013;
originally announced October 2013.
-
Path-integral Monte Carlo method for the local Z_2 Berry phase
Authors:
Yuichi Motoyama,
Synge Todo
Abstract:
We present a loop cluster algorithm Monte Carlo method for calculating the local Z_2 Berry phase of the quantum spin models. The Berry connection, which is given as the inner product of two ground states with different local twist angles, is expressed as a Monte Carlo average on the worldlines with fixed spin configurations at the imaginary-time boundaries. The "complex weight problem" caused by t…
▽ More
We present a loop cluster algorithm Monte Carlo method for calculating the local Z_2 Berry phase of the quantum spin models. The Berry connection, which is given as the inner product of two ground states with different local twist angles, is expressed as a Monte Carlo average on the worldlines with fixed spin configurations at the imaginary-time boundaries. The "complex weight problem" caused by the local twist is solved by adopting the meron cluster algorithm. We present the results of simulation on the antiferromagnetic Heisenberg model on an out-phase bond-alternating ladder to demonstrate that our method successfully detects the change in the valence bond pattern at the quantum phase transition point. We also propose that the gauge-fixed local Berry connection can be an effective tool to estimate pricisely the quantum critical point.
△ Less
Submitted 7 February, 2013; v1 submitted 14 November, 2012;
originally announced November 2012.
-
Screening and Fabrication of Half-Heusler phases for thermoelectric applications
Authors:
Wilfried Wunderlich,
Yuichiro Motoyama
Abstract:
Half-Heusler phases have gained recently much interest as thermoelectric materials. Screening of possible systems was performed by ab-initio simulation using VASP-software. The energy-versus-Volume (E(V)) curves were calculated and calibrated. For TiCoSb, NbNiSn, FeMoSb the stability of Half-Heusler phase against concurrent crystal structures like TiNiSi, ZrCoAl, ZrBeSi, FeSiV, ZrNiP and Full He…
▽ More
Half-Heusler phases have gained recently much interest as thermoelectric materials. Screening of possible systems was performed by ab-initio simulation using VASP-software. The energy-versus-Volume (E(V)) curves were calculated and calibrated. For TiCoSb, NbNiSn, FeMoSb the stability of Half-Heusler phase against concurrent crystal structures like TiNiSi, ZrCoAl, ZrBeSi, FeSiV, ZrNiP and Full Heusler was confirmed. However, the thermo-dynamical driving force as calculated from the difference in lattice energies is less than 0.1eV/atom. Hence, the fabrication of Half Heusler phases is a challenge and requires three steps, surface activation of the raw material by ball milling, arc-melting of pressed pellets and finally long-term annealing treatment in a vacuum furnace. On doped TiCoSb specimens, Seebeck coefficients up to 0.1 mV/K, on NiNbSn 0.16 mV/K were measured, although the microstructure was not yet optimized.
△ Less
Submitted 20 February, 2009; v1 submitted 11 January, 2009;
originally announced January 2009.
