Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNew Generation of the Monte Carlo Shell Model for the K Computer Era
128
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a newly enhanced version of the Monte Carlo Shell Model method by incorporating the conjugate gradient method and energy-variance extrapolation. This new method enables us to perform large-scale shell-model calculations that the direct diagonalization method cannot reach. This new generation framework of the MCSM provides us with a powerful tool to perform most-advanced large-scale shell-model calculations on current massively parallel computers such as the K computer. We discuss the validity of this method in ab initio calculations of light nuclei, and propose a new method to describe the intrinsic wave function in terms of the shell-model picture. We also apply this new MCSM to the study of neutron-rich Cr and Ni isotopes using the conventional shell-model calculations with an inert 40Ca core and discuss how the magicity of N = 28, 40, 50 remains or is broken.
rate research
Read More
Rotational motion of heated 72-Ge is studied within the microscopic Shell Model Monte Carlo approach. We investigate the the angular momentum alignment and nuclear pairing correlations associated with J-pi Cooper pairs as a function of the rotational frequency and temperature. The reentrance of pairing correlations with temperature is predicted at high rotational frequencies. It manifests itself through the anomalous behavior of specific heat and level density.
In the last years auxiliary field diffusion Monte Carlo has been used to assess the properties of hypernuclear systems, from light- to medium-heavy hypernuclei and hyper-neutron matter. One of the main findings is the key role played by the three-body hyperon-nucleon-nucleon interaction in the determination of the hyperon separation energy of hypernuclei and as a possible solution to the hyperon puzzle. However, there are still aspects of the employed hypernuclear potential that remain to be carefully investigated. In this paper we show that the isospin dependence of the Lambda-NN force, which is crucial in determining the NS structure, is poorly constrained by the available experimental data.
We consider the problem of including $Lambda$ hyperons into the ab initio framework of nuclear lattice effective field theory. In order to avoid large sign oscillations in Monte Carlo simulations, we make use of the fact that the number of hyperons is typically small compared to the number of nucleons in the hypernuclei of interest. This allows us to use the impurity lattice Monte Carlo method, where the minority species of fermions in the full nuclear Hamiltonian is integrated out and treated as a worldline in Euclidean projection time. The majority fermions (nucleons) are treated as explicit degrees of freedom, with their mutual interactions described by auxiliary fields. This is the first application of the impurity lattice Monte Carlo method to systems where the majority particles are interacting. Here, we show how the impurity Monte Carlo method can be applied to compute the binding energy of the light hypernuclei. In this exploratory work we use spin-independent nucleon-nucleon and hyperon-nucleon interactions to test the computational power of the method. We find that the computational effort scales approximately linearly in the number of nucleons. The results are very promising for future studies of larger hypernuclear systems using chiral effective field theory and realistic hyperon-nucleon interactions, as well as applications to other quantum many-body systems.
The configuration interaction method, which is well-known as the shell-model calculation in the nuclear physics community, plays a key role in elucidating various properties of nuclei. In general, these studies require a huge number of shell-model calculations to be repeated for parameter calibration and quantifying uncertainties. To reduce these computational costs, we propose a new workflow of shell-model calculations using a method called eigenvector continuation (EC). It enables us to efficiently approximate the eigenpairs under a given Hamiltonian by previously sampled eigenvectors. We demonstrate the validity of EC as an emulator of the shell-model calculations for a valence space, where the dimension of parameters is relatively large compared to the previous studies using EC. We also discuss its possible applications to the quantification of theoretical uncertainty, using an example of Markov chain Monte Carlo sampling for a simplified problem. Furthermore, we propose a new usage of EC: preprocessing, in which we start the Lanczos iterations from the approximate eigenvectors, and demonstrate that this can accelerate the shell-model calculations and the subsequent research cycles. With the aid of the eigenvector continuation, the eigenvectors obtained during the parameter optimization are not necessarily to be discarded, even if their eigenvalues are far from the experimental data. Those eigenvectors can become accumulated knowledge. In order to enable efficient sampling of shell-model results and to demonstrate the usefulness of the methodology described above, we developed a new shell-model code, ShellModel.jl. This code is written in Julia language and therefore flexible to add extensions for the users purposes.
We report $J^pi = 0^+$ ground-state energies and point-proton radii of $^4$He, $^8$Be, $^{12}$C, $^{16}$O and $^{20}$Ne nuclei calculated by the {it ab initio} no-core Monte Carlo shell model with the JISP16 and Daejeon16 nonlocal $NN$ interactions. Ground-state energies are obtained in the basis spaces up to 7 oscillator shells ($N_{rm shell} = 7$) with several oscillator energies ($hbar omega$) around the optimal oscillator energy for the convergence of ground-state energies. These energy eigenvalues are extrapolated to obtain estimates of converged ground state energies in each basis space using energy variances of computed energy eigenvalues. We further extrapolate these energy-variance-extrapolated energies obtained in the finite basis spaces to infinite basis-space results with an empirical exponential form. This form features a dependence on the basis-space size but is independent of the $hbaromega$ used for the harmonic-oscillator basis functions. Point-proton radii for these states of atomic nuclei are also calculated following techniques employed for the energies. From these results, it is found that the Daejeon16 $NN$ interaction provides good agreement with experimental data up to approximately $^{16}$O, while the JISP16 $NN$ interaction provides good agreement with experimental data up to approximately $^{12}$C. Beyond these nuclei, the interactions produce overbinding accompanied by radii that are too small. These findings suggest and encourage further revisions of nonlocal $NN$ interactions towards the investigation of nuclear structure in heavier-mass regions.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
