Professor Alarcon’s lab is producing proton beam detectors, and this project is focused on informing the decision as to which layout of detector is more effective at producing an accurate backprojection for an equal number of data channels. The comparison is between “square pad” detectors and “wire pad” detectors. The square pad detector consists of a grid of square pads all of identical size, that each collect their own data. The wire pad detector consists of large rectangular pads that span the entire detector in one direction, with 2 additional layers of identical pads each rotated by 60° from the previous. In order to test each design Python was used to simulate Gaussian beams of varying amplitudes, position and size and integrate them in each of the two methods. They were then backprojected and fit to a Gaussian function and the error between the backprojected parameters and the original parameters of the beam were measured.

In this dissertation two kinds of strongly interacting fermionic systems were studied: cold atomic gases and nucleon systems. In the first part I report T=0 diffusion Monte Carlo results for the ground-state and vortex excitation of unpolarized spin-1/2 fermions in a two-dimensional disk. I investigate how vortex core structure properties behave over the BEC-BCS crossover. The vortex excitation energy, density profiles, and vortex core properties related to the current are calculated. A density suppression at the vortex core on the BCS side of the crossover and a depleted core on the BEC limit is found. Size-effect dependencies in the disk geometry were carefully studied. In the second part of this dissertation I turn my attention to a very interesting problem in nuclear physics. In most simulations of nonrelativistic nuclear systems, the wave functions are found by solving the many-body Schrödinger equations, and they describe the quantum-mechanical amplitudes of the nucleonic degrees of freedom. In those simulations the pionic contributions are encoded in nuclear potentials and electroweak currents, and they determine the low-momentum behavior. By contrast, in this work I present a novel quantum Monte Carlo formalism in which both relativistic pions and nonrelativistic nucleons are explicitly included in the quantum-mechanical states of the system. I report the renormalization of the nucleon mass as a function of the momentum cutoff, an Euclidean time density correlation function that deals with the short-time nucleon diffusion, and the pion cloud density and momentum distributions. In the two nucleon sector the interaction of two static nucleons at large distances reduces to the one-pion exchange potential, and I fit the low-energy constants of the contact interactions to reproduce the binding energy of the deuteron and two neutrons in finite volumes. I conclude by showing that the method can be readily applied to light-nuclei.

Observations of four times ionized iron and nickel (Fe V & Ni V) in the G191-B2B white dwarf spectrum have been used to test for variations in the fine structure constant, α, in the presence of strong gravitational fields. The laboratory wavelengths for these ions were thought to be the cause of inconsistent conclusions regarding the
variation of α as observed through the white dwarf spectrum. This thesis presents 129 revised Fe V wavelengths (1200 Å to 1600 Å) and 161 revised Ni V wavelengths (1200 Å to 1400 Å) with uncertainties of approximately 3 mÅ. A systematic calibration error
is identified in the previous Ni V wavelengths and is corrected in this work. The evaluation of the fine structure variation is significantly improved with the results
found in this thesis.

Our group has constructed a ring-imaging Cherenkov (RICH) detector with the goal of testing the performance of aerogel tiles in charged particle detectors. In previous work, tiles produced by Aspen Aerogels were tested as radiators in Cherenkov threshold counters and compared to commercial-grade samples. As an extension of this work we built a counter of the RICH type, which is used in practice to extract more particle identification information than threshold counters, and we have studied the images resulting from various aerogel samples.
The detector was designed for use in table-top experiments in which our particle source would be cosmic rays. Due to the vast energy range of cosmic rays, the window in which we can discriminate velocities is relatively small. Since the particles we do detect generally have β≈1, the relativistic limit β→1 motivates imaging by the Focusing Aerogel RICH (FARICH) technique, in which photons from multiple tiles are focused together at a detection plane.
Our detection plane is an array of flat-panel, multi-anode photomultiplier tubes (PMTs). Readout consists of multiplexing the anode outputs, recording the digitized signal, and converting this into a matrix of integrated charge values. The charge distribution in that matrix should directly imply the particle's speed; however, in practice, final recorded images are the influenced by many intermediate processes, so we have attempted to make meaningful measurements by averaging over numerous events.
For a given configuration and data collection, we produce the spatial distribution of observed signals relative to the cosmic ray's point of impact. These distributions have the expected form of a ring and their characteristics compare favorably with the predictions of geometric optics. Our confidence in the images is increased by verifying that changes to the configuration are reflected by the changes in the rings. We find that FARICH improves the sharpness of our ring images, but tiles must be used individually for actual aerogel analysis. So far we have shown that the Aspen tiles behave as one would expect for the purposes of RICH. Their images do resemble those produced by commercial-grade tiles, but we do not have tiles sufficiently similar for side-by-side comparison. A method of quantifying tile performance has proven difficult and is the only remaining task for our group.

A series of experiments using a polarized beam incident on a polarized frozen spin target

(FROST) was conducted at Jefferson Lab in 2010. Results presented here were taken

during the second running period with the FROST target using the CEBAF Large Acceptance

Spectrometer (CLAS) detector at Jefferson Lab, which used transversely-polarized

protons in a butanol target and a circularly-polarized incident tagged photon beam with

energies between 0.62 and 2.93 GeV. Data are presented for the F and T polarization observables

for h meson photoproduction on the proton from W = 1.55 GeV to 1.80 GeV.

The data presented here will improve the world database and refine theoretical approaches

of nucleon structure.

The muon problem of flavor physics presents a rich opportunity to study beyond standard model physics. The as yet undiscovered bound state (μ+μ-), called true muonium, presents a unique opportunity to investigate the muon problem. The near-future experimental searches for true muonium will produce it relativistically, preventing the easy application of non-relativistic quantum mechanics. In this thesis, quantum field theory methods based on light-front quantization are used to solve an effective Hamiltonian for true muonium in the Fock space of |μ+μ-> , |μ+μ-γ> , |e+e->, |e+e-γ>, |τ+τ-> , and |τ+τ-γ> . To facilitate these calculations a new parallel code, True Muonium Solver With Front-Form Techniques (TMSWIFT), has been developed. Using this code, numerical results for the wave functions, energy levels, and decay constants of true muonium have been obtained for a range of coupling constants α. Work is also presented for deriving the effective interaction arising from the |γγ sector’s inclusion into the model.

Proton and fluorine diffusivity and ionic conductivity of 2-fluoropyridinium triflate (2-FPTf) and proton and fluorine diffusivity, ionic conductivity, and viscosity of trifluoromethanesulfonic acid (TFMSA) monohydrate have been measured over a wide range of temperatures. Diffusivities were measured using the pulsed-gradient spin-echo (PGSE) technique on a 300 MHz NMR spectrometer. Conductivities were measured using electrochemical impedance spectroscopy (EIS) on standard equipment and viscosities were determined using a Cannon-Ubbelohde viscometer. For 2-FPTF, the diffusivity of mobile protons increased from 1.84+/-0.06 x 10(-11) m2/s at 55 degC to 1.64+/-0.05 x 10(-10) m2/s at 115 degC while the diffusivity of 2-fluoropyridine fluorines increased from 2.22+/-0.07 x 10(-11) m2/s at 70 degC to 9.40+/-0.28 x 10(-11) m2/s at 115 degC. For TFMSA monohydrate, the diffusivity of protons increased from 7.67+/-0.23 x 10(-11) m2/s at 40 degC to 3.92+/-0.12 x 10(-10) m2/s at 110 degC while the diffusivity of fluorines increased from 4.63+/-0.14 x 10(-11) m2/s at 40 degC to 2.79+/-0.08 x 10(-10) m2/s at 110 degC, protons diffusing approximately 1.5 times faster than fluorines over the entire temperature range. NMR spectra indicate that proton diffusion occurs via direct hopping from TFMSA molecule to molecule. The conductivity of 2-FPTf varied from 0.85+/-0.03 mS/cm to 35.9+/-1.08 mS/cm between 25 and 110 degC. The conductivity of TFMSA monohydrate varied from 6.60+/-0.2 mS/cm to 84.6+/-2.5 mS/cm between 23 and 139 degC and its viscosity varied from 27.5+/-1.4 mPa.s to 4.38+/-0.22 mPa.s between 49 and 121.5 degC, in good agreement with literature values. Temperature dependences of the measured properties showed Arrhenius behavior with activation energies for proton diffusion, fluorine diffusion and ionic conduction for 2-FPTf above the melting point of 16.9+/-0.8 kJ/mol, 48.0+/-2.4 kJ/mol and 27.8+/-1.4 kJ/mol respectively. Activation energies for proton diffusion, fluorine diffusion, ionic conduction and viscosity for TFMSA monohydrate were 23.4+/-1.2 kJ/mol, 26.0+/-1.3 kJ/mol, 22.1+/-1.1 kJ/mol, and 26.9+/-1.3 kJ/mol respectively. The degree of dissociation of the charged species, calculated using the Nernst-Einstein relation, varied from 13 to 24% for 2-FPTf and from 25 to 29% for TFMSA monohydrate over the temperature range.

With the discovery of the Higgs Boson in 2012, particle physics has decidedly moved beyond the Standard Model into a new epoch. Though the Standard Model particle content is now completely accounted for, there remain many theoretical issues about the structure of the theory in need of resolution. Among these is the hierarchy problem: since the renormalized Higgs mass receives quadratic corrections from a higher cutoff scale, what keeps the Higgs boson light? Many possible solutions to this problem have been advanced, such as supersymmetry, Randall-Sundrum models, or sub-millimeter corrections to gravity. One such solution has been advanced by the Lee-Wick Standard Model. In this theory, higher-derivative operators are added to the Lagrangian for each Standard Model field, which result in propagators that possess two physical poles and fall off more rapidly in the ultraviolet regime. It can be shown by an auxiliary field transformation that the higher-derivative theory is identical to positing a second, manifestly renormalizable theory in which new fields with opposite-sign kinetic and mass terms are found. These so-called Lee-Wick fields have opposite-sign propagators, and famously cancel off the quadratic divergences that plague the renormalized Higgs mass. The states in the Hilbert space corresponding to Lee-Wick particles have negative norm, and implications for causality and unitarity are examined.

This dissertation explores a variant of the theory called the N = 3 Lee-Wick

Standard Model. The Lagrangian of this theory features a yet-higher derivative operator, which produces a propagator with three physical poles and possesses even better high-energy behavior than the minimal Lee-Wick theory. An analogous auxiliary field transformation takes this higher-derivative theory into a renormalizable theory with states of alternating positive, negative, and positive norm. The phenomenology of this theory is examined in detail, with particular emphasis on the collider signatures of Lee-Wick particles, electroweak precision constraints on the masses that the new particles can take on, and scenarios in early-universe cosmology in which Lee-Wick particles can play a significant role.

Spin-orbit interactions are important in determining nuclear structure. They lead to a shift in the energy levels in the nuclear shell model, which could explain the sequence of magic numbers in nuclei. Also in nucleon-nucleon scattering, the large nucleon polarization observed perpendicular to the plane of scattering needs to be explained by adding the spin-orbit interactions in the potential. Their effects change the equation of state and other properties of nuclear matter. Therefore, the simulation of spin-orbit interactions is necessary in nuclear matter.

The auxiliary field diffusion Monte Carlo is an effective and accurate method for calculating the ground state and low-lying exited states in nuclei and nuclear matter. It has successfully employed the Argonne v6' two-body potential to calculate the equation of state in nuclear matter, and has been applied to light nuclei with reasonable agreement with experimental results. However, the spin-orbit interactions were not included in the previous simulations, because the isospin-dependent spin-orbit potential is difficult in the quantum Monte Carlo method. This work develops a new method using extra auxiliary fields to break up the interactions between nucleons, so that the spin-orbit interaction with isospin can be included in the Hamiltonian, and ground-state energy and other properties can be obtained.

For this project, the diffuse supernova neutrino background (DSNB) has been calculated based on the recent direct supernova rate measurements and neutrino spectrum from SN1987A. The estimated diffuse electron antineutrino flux is ∼ 0.10 – 0.59 /cm2/s at 99% confidence level, which is 5 times lower than the Super-Kamiokande 2012 upper limit of 3.0 /cm2/s, above energy threshold of 17.3 MeV. With a Megaton scale water detector, 40 events could be detected above the threshold per year. In addition, the detectability of neutrino bursts from direct black hole forming collapses (failed supernovae) at Megaton detectors is calculated. These neutrino bursts are energetic and with short time duration, ∼ 1s. They could be identified by the time coincidence of N ≥2 or N ≥3 events within 1s time window from nearby (4 – 5 Mpc) failed supernovae. The detection rate of these neutrino bursts could get up to one per decade. This is a realistic way to detect a failed supernova and gives a promising method for studying the physics of direct black hole formation mechanism. Finally, the absorption of ultra high energy (UHE) neutrinos by the cosmic neutrino background, with full inclusion of the effect of the thermal distribution of the background on the resonant annihilation channel, is discussed. Results are applied to serval models of UHE neutrino sources. Suppression effects are strong for sources that extend beyond z ∼ 10. This provides a fascinating probe of the physics of the relic neutrino background in the unexplored redshift interval z ∼ 10 – 100. Ultimately this research will examine the detectability of DSNB, neutrino bursts from failed supernovae and absorption effects in the neutrino spectrum.