Computational Approaches to the Nuclear Many-Body Problem

<p>The nuclear many-body problem is conceptualized in an infinite-dimensional Hilbert space, but computationally solved in a finite one. Thus, the predictive power of microscopic calculations relies on the truncated representation of the infinite-dimensional space as well as leveraging advanced computational methods. My dissertation research focuses on three problems related to computational nuclear physics: exploring aspects of nuclear structure, efficiently solving the large sparse matrix-eigenvalue problem and improving the construction of the many-body basis in the no-core configuration-interaction (NCCI) framework.</p> <p>I) Elliott’s rotational SU(3) model, and its later extension, the symplectic Sp(3, <em>R</em>) model, both played a foundational role in improving the description of nuclear rotational spectral bands. My study of several beryllium isotopes, and 20Ne uses the decomposition of no-core shell-model wavefunctions into symmetry defined subspaces to show the Sp(3, <em>R</em>) picture provides a more consistent description of rotational band structure.</p> <p>II) Solving the non-relativistic many-body Schrödinger equation is often cast as a large sparse Hamiltonian eigenvalue problem. State-of-the-art NCCI calculation dimensions can exceed several billion and typically require supercomputers and thousands of core hours to compute small numbers of low-lying eigenstates. Thus, there is strong motivation for ways to reduce computational costs. In this research, I augment the block Lanczos algorithm using a <em>bootstrapped</em> pivot to significantly reduce the number of Hamiltonian-matrix multiplications typically dominating the algorithm’s total time-to-solution. My results demonstrate significant speedup in time-to-solution, often by a factor of two or more, and up to ten, can be achieved through the use of bootstrapping.</p> <p>III) In NCCI, the many-body basis is historically constructed from antisymmeterized products of harmonic oscillator (HO) single-particle wavefunctions. However, one often needs many HO antisymmeterized products states to produce accurate theoretical predictions of the properties of low-lying nuclear states. Alternative choices of single-particle basis which provide better descriptions of nuclear observables relative to the problem dimension and underlying basis parameters motivate continued explorations. In this research, I explore the use of a natural orbital (NO) single-particle basis, that is, one which diagonalizes the one-body density matrix of a reference many-body state, as a means of improving the description of energy, electromagnetic transitions and radii calculations relative to the problem dimension for select <em>sd</em>-shell nuclei.</p>

Name thesis:degree_name

Computational Science Joint PhD with San Diego State University, PhD

Level thesis:degree_level

Open Access Dissertation

Discipline thesis:degree_discipline

Institute of Mathematical Sciences

Year dc:date.available

2024

Author dc:creator

• Zbikowski, Ryan M

Contributors dc:contributor

• Fridolin Weber
• Peter Blomgren
• Marina Chugunova & Ali Nadim

• Computational Physics
• Computational Science
• Lanczos Algorithm
• Nuclear Structure

Repository record dc:identifier

https://scholarship.claremont.edu/cgu_etd/816

OAI identifier oai:identifier

oai:scholarship.claremont.edu:cgu_etd-1838