Advisor Feed
https://feeds.library.caltech.edu/people/Koonin-S-E/advisor.rss
A Caltech Library Repository Feedhttp://www.rssboard.org/rss-specificationpython-feedgenenTue, 16 Apr 2024 15:31:04 +0000Application of TDHF Methods to Nuclear Physics
https://resolver.caltech.edu/CaltechTHESIS:03022018-134310889
Authors: {'items': [{'id': 'Flanders-Bradley-A', 'name': {'family': 'Flanders', 'given': 'Bradley A.'}}]}
Year: 1981
DOI: 10.7907/vakw-ws15
<p>Part I presents results for TDHF calculations of realistic heavy-ion reactions. Results for a separable approximation, neglecting motion normal to the scattering plane, agree very well with the results of the full 3-D calculation. Results for the fusion cross section of two systems leading to the compound nucleus <sup>56</sup>Ni are compared with experimental data. The <sup>16</sup>O + <sup>40</sup>Ca results agree quite well with the experimental data, and the <sup>28</sup>Si + <sup>28</sup>Si results agree quite well with experimental data for the similar <sup>32</sup>S + <sup>27</sup>Al system. Results of calculations with the separable approximation for <sup>86</sup>Kr + <sup>139</sup>La are compared with both axially symmetric calculations and experimental results. All three show substantial agreement.</p>
<p>Part II presents a stability criterion for the validity of TDHF solutions based on a time-dependent generalization of RPA theory. Results are tested in an exactly soluble model, the SU(3) generalization of the Lipkin model. Unfortunately, the exact solution could not be computed for a large enough number of particles to permit quantitative testing of this criterion. However, it does seem to correctly indicate the stability or instability of the TDHF path.</p>https://thesis.library.caltech.edu/id/eprint/10747A Proximity Formulation of Nuclear Dynamics
https://resolver.caltech.edu/CaltechTHESIS:12082010-110205851
Authors: {'items': [{'email': 'gjball@gmail.com', 'id': 'Ball-Gregory-John', 'name': {'family': 'Ball', 'given': 'Gregory John'}, 'show_email': 'NO'}]}
Year: 1981
DOI: 10.7907/0H2W-T919
<p>The nuclear potential, the transfer-induced dissipation, and the mass diffusion coefficient in heavy-ion collisions are investigated in a proximity formulation. An energy-dependent nuclear potential is calculated in the frozen wave function approximation using two slabs of symmetric nuclear matter, each described by Hartree-Fock single-particle wave functions. Corrections to the inertia parameter are also evaluated from this potential. The flux entering the window formula for the friction between two heavy ions is calculated in a simple barrier penetration model. The classically forbidden flux is found to make a significant contribution. The transfer flux arising from both the relative motion and finite temperature of the nuclei is calculated and the latter is used to estimate the mass diffusion coefficient. Using the mean trajectories from time-dependent Hartree-Fock calculations the charge variance is calculated for the reaction <sup>84</sup>Kr(712 MeV) + <sup>209</sup>Bi and is found to be in agreement with experiment.</p>https://thesis.library.caltech.edu/id/eprint/6198Realistic Calculations of Excitations in Nuclear Matter
https://resolver.caltech.edu/CaltechTHESIS:10182019-162452942
Authors: {'items': [{'id': 'Kwong-Nai-Hang', 'name': {'family': 'Kwong', 'given': 'Nai-Hang'}, 'show_email': 'NO'}]}
Year: 1983
DOI: 10.7907/wcb1-mh63
<p>A numerical method has been developed to solve the RPA equation, exchange term included, in nuclear matter. The dynamic form factor S(q, ω) is extracted for several v4 and v6 phenomenological potentials, including the d1-potential of Gogny et al. The limits of validity of the long-wavelength (Landau) approximation and the often adopted local-kernel approximation are discussed. Substantial disagreements with the exact results are found for the latter. The method is then applied to solve a Jastrow-correlated extension of the RPA equation, using the hardcore OMY potential. Results of calculations performed in two-body cluster approximation and Fermi-Hypernetted-Chain (FHNC) approximation are compared. The two-body results predict an instability against density fluctuations, which disappears at the FHNC level. The validity and consequences of employing the FHNC effective potential within the self-consistent HF/RPA framework are discussed. Future developments include applying the method to other Fermi systems such as liquid <sup>3</sup>He and the microscopic calculation of Landau parameters.</p>https://thesis.library.caltech.edu/id/eprint/11828Semiclassical Quantization in Many Dimensions
https://resolver.caltech.edu/CaltechTHESIS:10312019-103914805
Authors: {'items': [{'id': 'Williams-Roy-David', 'name': {'family': 'Williams', 'given': 'Roy David'}, 'show_email': 'NO'}]}
Year: 1983
DOI: 10.7907/ezwt-k047
<p>We examine the semiclassical limit of the quantum energy spectrum in many dimensions: by means of a WKB-like ansatz leading to Einstein-Brillouin-Keller (EBK) quantization, by means of a path integral, hence associating a bound state with a particular classical periodic trajectory, and by the Birkhoff-Gustavson (BG) transformation to action-angle variables. We extend the EBK method to many-fermion systems using coherent states; and apply both EBK using surfaces of section, and the BG transformation to an SU(3) schematic nuclear shell model. We describe a new algorithm for finding periodic trajectories of a Lagrangian system with polynomial potential. It is applied to the Henon-Heiles system with good results, and these trajectories are used to quantize the system. The EBK and BG methods have some success, while periodic trajectory quantization fails. We discuss possible reasons for this failure and future approaches to these problems.</p>https://thesis.library.caltech.edu/id/eprint/11876Monte Carlo Calculation of Three-Body Scattering
https://resolver.caltech.edu/CaltechTHESIS:01252019-102522205
Authors: {'items': [{'id': 'Mahale-Narayan-Krishna', 'name': {'family': 'Mahale', 'given': 'Narayan Krishna'}, 'show_email': 'NO'}]}
Year: 1985
DOI: 10.7907/09kg-d808
<p>We construct an eigenvalue problem by confining many-body system to a bounded domain with the boundary condition that the wave function vanishes. By changing the boundary, however, the eigenvalues of the energy can be varied continuously. The D-matrix is defined for a series of bounded problems with the same value for the ground state energy. The D-matrix is related to the S-matrix, enabling us to calculate the S-matrix at a given energy. The Schrodinger equation for the system is transformed to a diffusion equation by regarding time as imaginary. Initial ensemble, representing an approximate wave function, is evolved, through Monte Carlo simulation of random walks and branching, to the ground state ensemble. The limitations of investigation are: 1. Ingoing and outgoing channels have two fragments. 2. The interaction between the fragments is negligible outside the boundary mentioned above. 3. The particles are bosons or we know the zeros of the wave function.</p>
<p>First we consider the scattering of a particle by a potential, which is equivalent to the two-body problem, in one dimension. Here we use the Poschl-Teller potential for which the exact solution is known. We use this case to investigate a new sampling method and study of various parameters. Next we consider three particles in one dimension. Here we take interaction to be a potential well, where at least one of the interactions is attractive so that a two-body bound state is possible.</p>https://thesis.library.caltech.edu/id/eprint/11362Auxiliary Field Monte-Carlo for Quantum Many-Body Systems
https://resolver.caltech.edu/CaltechTHESIS:01072019-110818820
Authors: {'items': [{'id': 'Sugiyama-Gayle-A', 'name': {'family': 'Sugiyama', 'given': 'Gayle A.'}, 'show_email': 'NO'}]}
Year: 1985
DOI: 10.7907/t9kp-qe57
<p>An algorithm is developed for determining the exact ground state properties of quantum many-body systems which is equally applicable to bosons and fermions. The Schroedinger eigenvalue equation for the ground state energy is recast into the form of a many-dimensional integral through the use of the Hubbard-Stratonovitch representation of the imaginary time many- body evolution operator. The resulting functional integral is then evaluated stochastically. The algorithm is tested for an exactly soluble boson system and is then extended to include fermions and repulsive potentials. Importance sampling is crucial to the success of the method, particularly for more complex systems. Improved computational efficiency is attained by performing the calculations in momentum space.</p>https://thesis.library.caltech.edu/id/eprint/11335Subbarier Fusion of the Oxygen Isotopes
https://resolver.caltech.edu/CaltechTHESIS:06302010-080153741
Authors: {'items': [{'id': 'Pönisch-Volker', 'name': {'family': 'Pönisch', 'given': 'Volker'}, 'show_email': 'NO'}]}
Year: 1986
DOI: 10.7907/KZX1-P789
<p>The subbarrier fusion process is studied for systems involving oxygen isotopes. A one-channel incoming wave boundary condition (IWBC) calculation gives an excellent fit to fusion cross section data for <sup>16</sup>O+<sup>16</sup>O. An IWBC coupled channels calculation for <sup>17</sup>O+<sup>16</sup>O that includes inelastic excitations as well as one-neutron transfer with formfactors calculated in a consistent single-particle framework reproduces the subbarrier enhancement down to four fifths the barrier height, but not below that. The calculation does not invoke the adiabatic approximation, which would create non-unitarity in the coupled channels equations. The measured subbarrier fusion cross section for <sup>18</sup>O+<sup>16</sup>O is well reproduced by an IWBC coupled channels calculation with two-neutron transfer, but the calculation disagrees with the above-barrier data.</p>
https://thesis.library.caltech.edu/id/eprint/5966Quantum Chaos: Spectral Fluctuations and Overlap Distributions of the Three Level Lipkin-Meshkov-Glick Model
https://resolver.caltech.edu/CaltechTHESIS:02202019-120226199
Authors: {'items': [{'id': 'Meredith-Dawn-Christine', 'name': {'family': 'Meredith', 'given': 'Dawn Christine'}, 'show_email': 'NO'}]}
Year: 1987
DOI: 10.7907/7y5k-ex17
<p>We test the prediction that quantum systems with chaotic classical analogs have spectral fluctuations and overlap distributions equal to those of the Gaussian Orthogonal Ensemble (GOE). The subject of our study is the three level Lipkin-Meshkov-Glick model of nuclear physics. This model differs from previously investigated systems because the quantum basis and classical phase space are compact, and the classical Hamiltonian has quartic momentum dependence. We investigate the dynamics of the classical analog to identify values of coupling strength and energy ranges for which the motion is chaotic, quasi-chaotic, and quasi-integrable. We then analyze the fluctuation properties of the eigenvalues for those same energy ranges and coupling strength, and we find that the chaotic eigenvalues are in good agreement with GOE fluctuations, while the quasi-integrable and quasichaotic levels fluctuations are closer to the Poisson fluctuations that are predicted for integrable systems. We also study the distribution of the overlap of a chaotic eigenvector with a basis vector, and find that in some cases it is a Gaussian random variable as predicted by GOE. This result, however, is not universal. </p>
https://thesis.library.caltech.edu/id/eprint/11406Hydrodynamics of Ultra-Relativistic Heavy-Ion Collisions
https://resolver.caltech.edu/CaltechTHESIS:08302019-142642675
Authors: {'items': [{'id': 'Chu-Ming-chung', 'name': {'family': 'Chu', 'given': 'Ming-chung'}, 'show_email': 'NO'}]}
Year: 1987
DOI: 10.7907/kq9w-6f74
<p>Relativistic hydrodynamic calculations are presented to describe the dynamics
of ultra-relativistic heavy-ion collisions. In contrast to the "standard picture" of
the field, our calculations do not assume scaling symmetry, and in fact we find
large scaling violations near the fragmentation regions. In our 1+1-dimensional
calculations, we find that while the hydrodynamic evolution is very sensitive to the
formation and thermalization time and to the models of the source terms, the effects
of changing the viscosity and the equation of state are small. Our 2+1-dimensional
calculations show that transverse expansion is not important in the central rapidity
region. We also present a brief review of the proposed signatures of the formation of
quark-gluon plasma in high energy heavy-ion collisions, as examples of applications
of hydrodynamics.</p>https://thesis.library.caltech.edu/id/eprint/11780Two Problems in Many Body Physics: I. Photon Production: A Probe of Heavy Ion Collisions. II. Structure of Matter in Strong Magnetic Fields
https://resolver.caltech.edu/CaltechTHESIS:10222019-115047089
Authors: {'items': [{'id': 'Neuhauser-Daniel', 'name': {'family': 'Neuhauser', 'given': 'Daniel'}, 'show_email': 'NO'}]}
Year: 1987
DOI: 10.7907/dn5s-g690
<p>In Part 1 we examine the theoretical framework for the use of photon spectra to probe heavy ion collisions. We first calculate single photon emission spectra from nuclear matter in the incoherent limit and combine them with the simplified participant-spectator model and with the semiclassical VUU model, to predict photon production cross sections in heavy ion collisions. The spectra differ from previous estimates based on a classical soft-photon approximation and lead to good agreement with experiment, except for an overall normalization factor of order (2-5), which we interpret as direct evidence for medium effects.</p>
<p>We then proceed to examine the Hanbury-Brown-Twiss correlation of high-energy photons emitted from heavy ion collisions. We find that both the polarization average and a possible coherent component complicate the extraction of the size and lifetime of the emitting source from the correlation function.</p>
<p>In Part 2 we calculate the binding energies of atoms and molecular chains in 10<sup>12</sup> G magnetic fields using the Hartree-Fock method. Our calculations are the first self-consistent ones treating exchange properly for atoms heavier than helium in high fields. For <i>Z</i> > 2 at 10<sup>12</sup>G and <i>Z</i> > 4 at 5 x 10<sup>12</sup>G the isolated atom is energetically favored over the molecular chains.</p>https://thesis.library.caltech.edu/id/eprint/11840Scaling Violations in Quasielastic Electron Scattering
https://resolver.caltech.edu/CaltechTHESIS:01172013-133508860
Authors: {'items': [{'id': 'Butler-Malcolm-Nicholas', 'name': {'family': 'Butler', 'given': 'Malcolm Nicholas'}, 'show_email': 'NO'}]}
Year: 1988
DOI: 10.7907/pp3m-h071
<p>We study the phenomenon of y-scaling in inclusive quasielastic electron scattering. Emphasis is placed on the approach to scaling at finite four-momentum transfers, and the effects of final state interactions. Brueckner-Goldstone perturbation theory for nuclear matter is used to perform a detailed, microscopic calculation of the dynamic structure function of nuclear matter. This is compared to the naive prediction of the Impulse approximation, and we find that the approach to scaling is quite different. The Brueckner-Goldstone approach reproduces the trends seen experimentally (Impulse approximation does not); however, there is still not good quantitative agreement with the data. We take this to be a possible sign of problems in the conventional nucleon-nucleon interactions studied.</p>https://thesis.library.caltech.edu/id/eprint/7404A Path Integral Monte Carlo Method for the Quasielastic Response
https://resolver.caltech.edu/CaltechETD:etd-06072007-080221
Authors: {'items': [{'email': 'Carraro@berkeley.edu', 'id': 'Carraro-Carlo', 'name': {'family': 'Carraro', 'given': 'Carlo'}, 'orcid': '0000-0002-3482-9226', 'show_email': 'NO'}]}
Year: 1990
DOI: 10.7907/mkyz-4067
<p>We formulate the quasielastic response of a non-relativistic many-body system at zero temperature in terms of ground state density matrix elements and real time path integrals that embody the final state interactions. While the former provide the weight for a conventional Monte Carlo calculation, the latter require a more sophisticated treatment. We argue that the recently developed Stationary Phase Monte Carlo technique can be used to study the approach to "Y-scaling." We perform calculations for a particle in a potential well in one and three dimensions and compare them to the exact results available for these models. We then derive an eikonal approximation to the Path Integrals. This method is suitably generalized to treat strongly repulsive interactions, and allows comparison to Silver's theory of final state interactions in a straightforward way. We also give an exact prescription to calculate the scaling limit for potentials comprising a hard core. Finally, we study the approach to scaling in a model ⁴He nucleus, and find good agreement with experimental data.</p>https://thesis.library.caltech.edu/id/eprint/2498Relativistic mean field theory: methods and applications
https://resolver.caltech.edu/CaltechTHESIS:03182015-112613348
Authors: {'items': [{'id': 'Wasson-D-A', 'name': {'family': 'Wasson', 'given': 'David Allan'}, 'show_email': 'NO'}]}
Year: 1990
DOI: 10.7907/gs6h-m244
<p>We develop a method for performing one-loop calculations in finite systems that is based on using the WKB approximation for the high energy states. This approximation allows us to absorb all the counterterms analytically and thereby avoids the need for extreme numerical precision that was required by previous methods. In addition, the local approximation makes this method well suited for self-consistent calculations. We then discuss the application of relativistic mean field methods to the atomic nucleus. Self-consistent, one loop calculations in the Walecka model are performed and the role of the vacuum in this model is analyzed. This model predicts that vacuum polarization effects are responsible for up to five percent of the local nucleon density. Within this framework the possible role of strangeness degrees of freedom is studied. We find that strangeness polarization can increase the kaon-nucleus scattering cross section by ten percent. By introducing a cutoff into the model, the dependence of the model on short-distance physics, where its validity is doubtful, is calculated. The model is very sensitive to cutoffs around one GeV.</p>
https://thesis.library.caltech.edu/id/eprint/8785Calculations of quasielastic electron scattering from nuclei at high momentum transfers
https://resolver.caltech.edu/CaltechTHESIS:09022011-110852198
Authors: {'items': [{'id': 'Allder-Anne', 'name': {'family': 'Allder', 'given': 'Anne'}, 'show_email': 'NO'}]}
Year: 1992
DOI: 10.7907/d08g-aa33
I perform calculations of high momentum-transfer,
quasielastic electron scattering from nuclei using a
relativistic, microscopic model of the reaction that includes the effects of final-state interactions via an eikonal approximation and the effects of color transparency via a cross-section that varies along the exit path of the proton. The results of this microscopic calculation agree to 5% with a semiclassical estimate based on the mean free path of the proton in the nuclear medium. The overall uncertainty in the calculation is on the order of 10%; this figure is dominated by uncertainties in the potential that describes the final state interaction. Color transparency effects an enhancement on the order of a few percent (in total cross-section) for Q^2=1 GeV^2, and ten percent for Q^2=10 GeV^2. The predictions of various specific transparency models differ by approximately 20%, and thus are probably not experimentally distinguishable for
Q^2<10 GeV^2.
https://thesis.library.caltech.edu/id/eprint/6648Auxiliary-field Monte Carlo methods for interacting fermions : application to the nuclear shell model
https://resolver.caltech.edu/CaltechTHESIS:10222009-160215418
Authors: {'items': [{'id': 'Lang-G-H-W', 'name': {'family': 'Lang', 'given': 'Gladys Hau-Wan'}, 'show_email': 'NO'}]}
Year: 1993
DOI: 10.7907/xh84-k642
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.
https://thesis.library.caltech.edu/id/eprint/5323Nuclear Hartree-Fock calculations on parallel computers
https://resolver.caltech.edu/CaltechTHESIS:04262013-131130293
Authors: {'items': [{'id': 'Chen-B-Q', 'name': {'family': 'Chen', 'given': 'Bing-Qing'}, 'show_email': 'NO'}]}
Year: 1994
DOI: 10.7907/kp30-db55
<p>Hartree-Fock (HF) calculations have had remarkable success in describing large nuclei at
high spin, temperature and deformation. To allow full range of possible deformations,
the Skyrme HF equations can be discretized on a three-dimensional mesh. However, such
calculations are currently limited by the computational resources provided by traditional
supercomputers. To take advantage of recent developments in massively parallel computing
technology, we have implemented the LLNL Skyrme-force static and rotational
HF codes on Intel's DELTA and GAMMA systems at Caltech.</p>
<p>We decomposed the HF code by assigning a portion of the mesh to each node, with
nearest neighbor meshes assigned to nodes connected by communication· channels. This
kind of decomposition is well-suited for the DELTA and the GAMMA architecture because
the only non-local operations are wave function orthogonalization and the boundary
conditions of the Poisson equation for the Coulomb field.</p>
<p>Our first application of the HF code on parallel computers has been the study of
identical superdeformed (SD) rotational bands in the Hg region. In the last ten years,
many SD rotational bands have been found experimentally. One very surprising feature
found in these SD rotational bands is that many pairs of bands in nuclei that differ
by one or two mass units have nearly identical deexcitation gamma-ray energies. Our
calculations of the five rotational bands in ^(192)Hg and ^(194)Pb show that the filling of
specific orbitals can lead to bands with deexcitation gamma-ray energies differing by at
most 2 keV in nuclei differing by two mass units and over a range of angular momenta comparable to that observed experimentally. Our calculations of SD rotational bands
in the Dy region also show that twinning can be achieved by filling or emptying some specific orbitals.</p>
<p>The interpretation of future precise experiments on atomic parity nonconservation
(PNC) in terms of parameters of the Standard Model could be hampered by uncertainties
in the atomic and nuclear structure. As a further application of the massively parallel
HF calculations, we calculated the proton and neutron densities of the Cesium isotopes
from A = 125 to A = 139. Based on our good agreement with experimental charge
radii, binding energies, and ground state spins, we conclude that the uncertainties in
the ratios of weak charges are less than 10^(-3), comfortably smaller than the anticipated experimental error.</p>https://thesis.library.caltech.edu/id/eprint/7635Shell Model Monte Carlo for Gamow-Teller Strengths and Two-Neutrino Double Beta Decay
https://resolver.caltech.edu/CaltechTHESIS:10312019-172947012
Authors: {'items': [{'id': 'Bahukutumbi-Radha-Pillapakkam', 'name': {'family': 'Bahukutumbi', 'given': 'Radha Pillapakkam'}, 'show_email': 'NO'}]}
Year: 1996
DOI: 10.7907/xg95-wy83
<p>In this thesis, a method to calculate two-neutrino double beta decay matrix elements employing the Shell Model Monte Carlo is presented. This method is validated against direct-diagonalization for the decay of <sup>48</sup>Ca. The first realistic calculation of the nuclear matrix element within the shell model for <sup>76</sup>Ge is performed; the result is in reasonable agreement with experiment.</p>
<p>The sensitivity of the shell model results to the nuclear Hamiltonian has been studied for the case of <sup>48</sup>Ca where the Hamiltonian used is known to be an optimal one. While one cannot make the nuclear matrix element arbitrarily small, the uncertainty in certain pieces of the Hamiltonian such as the monopole isovector pairing, provides room for at least a factor of two in the matrix element (and hence a factor of four in the half-life) from such calculations.</p>
<p>A Maximum Entropy method to obtain realistic strength functions from imaginary time response functions has been applied to Gamow-Teller response functions calculated using the Shell Model Monte Carlo and the results are validated against direct-diagonalization and experiment.</p>
<p>Future prospects for double beta decay calculations and astrophysical applications of the Gamow-Teller strength functions are briefly discussed.</p>https://thesis.library.caltech.edu/id/eprint/11879Fermionic quantum systems. Part I: Phase transitions in quantum dots. Part II: Nuclear matter on a lattice
https://resolver.caltech.edu/CaltechTHESIS:12042013-111115672
Authors: {'items': [{'id': 'Müller-H-M', 'name': {'family': 'Müller', 'given': 'Hans-Michael'}, 'show_email': 'NO'}]}
Year: 1999
DOI: 10.7907/yznz-2b69
<p>In the first part I perform Hartree-Fock calculations to show that quantum dots (i.e.,
two-dimensional systems of up to twenty interacting electrons in an external parabolic
potential) undergo a gradual transition to a spin-polarized Wigner crystal with increasing
magnetic field strength. The phase diagram and ground state energies have
been determined. I tried to improve the ground state of the Wigner crystal by introducing
a Jastrow ansatz for the wave function and performing a variational Monte
Carlo calculation. The existence of so called magic numbers was also investigated.
Finally, I also calculated the heat capacity associated with the rotational degree of
freedom of deformed many-body states and suggest an experimental method to detect
Wigner crystals.</p>
<p>The second part of the thesis investigates infinite nuclear matter on a cubic lattice.
The exact thermal formalism describes nucleons with a Hamiltonian that accommodates
on-site and next-neighbor parts of the central, spin-exchange and isospin-exchange
interaction. Using auxiliary field Monte Carlo methods, I show that energy
and basic saturation properties of nuclear matter can be reproduced. A first order
phase transition from an uncorrelated Fermi gas to a clustered system is observed
by computing mechanical and thermodynamical quantities such as compressibility,
heat capacity, entropy and grand potential. The structure of the clusters is investigated
with the help two-body correlations. I compare symmetry energy and first
sound velocities with literature and find reasonable agreement. I also calculate the
energy of pure neutron matter and search for a similar phase transition, but the survey
is restricted by the infamous Monte Carlo sign problem. Also, a regularization
scheme to extract potential parameters from scattering lengths and effective ranges
is investigated.</p>https://thesis.library.caltech.edu/id/eprint/8037