Recently proposed two methods such as Quantum Monte Carlo Diagonalization method in nuclear structure physics and path-integral renormalization group method in condensed matter physics have a similarity as quantum many-body approach, especially in description beyond mean-field approximation. As these two methods have been developed in different fields, there are significant differences as quantum-number projection technique in former method and extrapolation technique in latter method. As these techniques are independent of details of considered system, these can be applied to other fields. In this paper, we discuss how to apply such extrapolation technique to shell model calculations and discuss how to apply quantum-number projection to path-integral renormalization group method.
Recently proposed two methods such as Quantum Monte Carlo Diagonalization method in nuclear structure physics and path-integral renormalization group method in condensed matter physics have a similarity as quantum many-body approach, especially in description beyond mean-field approximation. As these two methods have been developed in different fields, there are significant differences as quantum-number projection technique in former method and extrapolation technique in latter method. As these techniques are independent of details of considered system, these can be applied to other fields. In this paper, we discuss how to apply such extrapolation technique to shell model calculations and discuss how to apply quantum-number projection to path-integral renormalization group method.
The formalism for the auxiliary-field Monte Carlo approach to the nuclear shell model is presented. The method is based on a linearization of the two-body part of the Hamiltonian in an imaginary-time propagator using the Hubbard-Stratonovich transformation. The foundation of the method, as applied to the nuclear many-body problem, is discussed. Topics presented in detail include: (1) the density-density formulation of the method, (2) computation of the overlaps, (3) the sign of the Monte Carlo weight function, (4) techniques for performing Monte Carlo sampling, and (5) the reconstruction of response functions from an imaginary-time auto-correlation function using MaxEnt techniques. Results obtained using schematic interactions, which have no sign problem, are presented to demonstrate the feasibility of the method, while an extrapolation method for realistic Hamiltonians is presented. In addition, applications at finite temperature are outlined.
The shell model provides a powerful framework for nuclear structure calculations. The nucleons beyond an inert magic core are confined in a valence shell and interact through an effective two-body potential generally determined from the G-matrix method. However, the applicability of the shell model remains limited by the exponential growth of the many-body space with the number of valence nucleons and with the size of the one-body space. The quantum Monte Carlo (QMC) methods may thus be considered as a potentially attractive alternative to the direct diagonalization of the Hamiltonian matrix. The many-body problem is indeed reduced to a set of stochastic one-body problems describing independent particles in a fluctuating external field. Up to date, the QMC approaches of the shell model lead to the ground-state and the thermodynamic properties, but a detailed spectroscopy of nuclei is not reachable. Furthermore, they suffer from a common pathology of fermionic QMC methods, namely the sign/phase problem, which strongly contaminates the convergence of the calculation. In this context, we will present a new QMC method, providing the "yrast spectroscopy" of nuclei, and in which the sign/phase problem is managed through the constrained path approximation that is widely used in the ab initio nuclear calculations and in condensed matter physics.
The formalism for the auxiliary-field Monte Carlo approach to the nuclear shell model is presented.The method is based on a linearization of the two-body part of the Hamiltonian in an imaginary-time propagator using the Hubbard-Stratonovich transformation.The foundation of the method, as applied to the nuclear many-body problem, is discussed.Topics presented in detail include: (1) the density-density formulation of the method, (2) computation of the overlaps, (3) the sign of the Monte Carlo weight function, (4) techniques for performing Monte Carlo sampling, and ( 5) the reconstruction of response functions from an imaginary-time auto-correlation function using MaxEnt techniques.Results obtained using schematic interactions, which have no sign problem, are presented to demonstrate the feasibility of the method, while an extrapolation method for realistic Hamiltonians is presented.In addition, applications at finite temperature are outlined.
We present a quantum-number projection technique which enables us to exactly treat spin, momentum, and other symmetries embedded in the Hubbard model. By combining this projection technique, we extend the path-integral renormalization-group method to improve the efficiency of numerical computations. By taking numerical calculations for the standard Hubbard model and the Hubbard model with next-nearest-neighbor transfer, we show that the present extended method can extremely enhance numerical accuracy and that it can handle excited states, in addition to the ground state.
The formulation and recent applications of the Quantum Monte Carlo diagonalization (QMCD) method are reported. The QMCD has been proposed for solving the quantum many-body inteacting systems, providing us with energy eigenvalues, transition matrix elements and wave functions. Its application to the nuclear shell model is referred to as the Monte Carlo Shell Model. By the Monte Carlo Shell Model calculations, the level structure of low-lying states can be studied with realistic interactions, providing a useful tool for nuclear spectroscopy. The Monte Carlo Shell Model has been applied to the study of a variety of nuclei, and can be characterized as the importance truncation scheme to the full diagonalization which is infeasible in many cases due to extremely large dimensions.
This thesis presents the path-integral formulation of the nuclear shell model using the Hubbard-Stratonovich transformation, which linearizes the two-body interaction by auxiliary fields. The path-integral was evaluated via Monte Carlo. The method scales favorably with valence-nucleon number and shell-model basis: full-basis calculations can be done up to the rare-earth region, which cannot be treated by other methods. Observables are calculated for the ground state and in a thermal ensemble. Dynamical correlations are obtained, from which strength functions are extracted through the Maximum Entropy method. Examples in the s-d shell, where exact diagonalization can be carried out, compare well with exact results. The "sign problem", which is generic to fermion Monte Carlo calculations, is proved to be absent in a wide class of interactions including the attractive pairing-plus-multipole interactions. The formulation is general for interacting fermion systems and is well suited for parallel computation. The method has been implemented on the Intel Touchstone Delta System, achieving better than 99% parallelization.