Source author record

J. Carlson

J. Carlson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

50works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

50 published item(s)

preprint2024arXiv

Equilibration of quantum many-body fast neutrino flavor oscillations

Neutrino gases are expected to form in high density astrophysical environments, and accurately modeling their flavor evolution is critical to understanding such environments. In this work we study a simplified model of such a dense neutrino gas in the regime for which neutrino-neutrino coherent forward scattering is the dominant mechanism contributing to the flavor evolution. We show evidence that the generic potential induced by this effect is non-integrable and that the statistics of its energy level spaces are in good agreement with the Wigner surmise. We also find that individual neutrinos rapidly entangle with all of the others present which results in an equilibration of the flavor content of individual neutrinos. We show that the average neutrino flavor content can be predicted utilizing a thermodynamic partition function. A random phase approximation to the evolution gives a simple picture of this equilibration. In the case of neutrinos and antineutrinos, processes like $ν_e {\barν}_e \leftrightarrows ν_μ{\barν_μ} $ yield a rapid equilibrium satisfying $n( ν_e) n({\bar ν}_e) = n( ν_μ) n({\bar ν}_μ) = n( ν_τ) n({\bar ν}_τ)$ in addition to the standard lepton number conservation in regimes where off-diagonal vacuum oscillations are small compared to $ν-ν$ interactions.

preprint2023arXiv

Many-body neutrino flavor entanglement in a simple dynamic model

Dense neutrino gases form in extreme astrophysical sites, and the flavor content of the neutrinos likely has an important impact on the subsequent dynamical evolution of their environment. Through coherent forward scattering among neutrinos, the flavor content of the gas evolves under a time-dependent potential which can be modeled in a quantum many-body formalism as an all-to-all coupled spin-spin interaction. This two-body potential generically introduces entanglement and greatly complicates the study of these systems. In this work we study the evolution of the quantum many-body problem as well as the typically employed mean-field approximation to it for a small number of neutrinos ($N = 16$). We consider randomly chosen one- and two-body couplings in the Hamiltonian, and the resulting evolution of several initial product states. We subsequently compare many-body and mean-field predictions for one-body observables, and we consider one- and two-body entanglement to assess under what conditions the many-body and mean-field predictions are likely to disagree. Except for a special category of prototypical initial conditions, we find that the typically employed mean-field approximation is insufficient to capture the evolution of one-body operators in the systems we consider. We also observe a loss of coherence in one- and two-body trace-reduced subsystems which suggests that the evolution may be well approximated as a classical mixture of separable states.

preprint2022arXiv

Classical and Quantum Evolution in a Simple Coherent Neutrino Problem

The extraordinary neutrino flux produced in extreme astrophysical environments like the early universe, core-collapse supernovae and neutron star mergers may produce coherent quantum neutrino oscillations on macroscopic length scales. The Hamiltonian describing this evolution can be mapped into quantum spin models with all-to-all couplings arising from neutrino-neutrino forward scattering. To date many studies of these oscillations have been performed in a mean-field limit where the neutrinos time evolve in a product state. In this paper we examine a simple two-beam model evolving from an initial product state and compare the mean-field and many-body evolution. The symmetries in this model allow us to solve the real-time evolution for the quantum many-body system for hundreds or thousands of spins, far beyond what would be possible in a more general case with an exponential number ($2^N$) of quantum states. We compare mean-field and many-body solutions for different initial product states and ratios of one- and two-body couplings, and find that in all cases in the limit of infinite spins the mean-field (product state) and many-body solutions coincide for simple observables. This agreement can be understood as a consequence of the fact that the typical initial condition represents a very local but dense distribution about a mean energy in the spectrum of the Hamiltonian. We explore quantum information measures like entanglement entropy and purity of the many-body solutions, finding intriguing relationships between the quantum information measures and the dynamical behavior of simple physical observables.

preprint2022arXiv

Quality Control of Mass-Produced GEM Detectors for the CMS GE1/1 Muon Upgrade

The series of upgrades to the Large Hadron Collider, culminating in the High Luminosity Large Hadron Collider, will enable a significant expansion of the physics program of the CMS experiment. However, the accelerator upgrades will also make the experimental conditions more challenging, with implications for detector operations, triggering, and data analysis. The luminosity of the proton-proton collisions is expected to exceed $2-3\times10^{34}$~cm$^{-2}$s$^{-1}$ for Run 3 (starting in 2022), and it will be at least $5\times10^{34}$~cm$^{-2}$s$^{-1}$ when the High Luminosity Large Hadron Collider is completed for Run 4. These conditions will affect muon triggering, identification, and measurement, which are critical capabilities of the experiment. To address these challenges, additional muon detectors are being installed in the CMS endcaps, based on Gas Electron Multiplier technology. For this purpose, 161 large triple-Gas Electron Multiplier detectors have been constructed and tested. Installation of these devices began in 2019 with the GE1/1 station and will be followed by two additional stations, GE2/1 and ME0, to be installed in 2023 and 2026, respectively. The assembly and quality control of the GE1/1 detectors were distributed across several production sites around the world. We motivate and discuss the quality control procedures that were developed to standardize the performance of the detectors, and we present the final results of the production. Out of 161 detectors produced, 156 detectors passed all tests, and 144 detectors are now installed in the CMS experiment. The various visual inspections, gas tightness tests, intrinsic noise rate characterizations, and effective gas gain and response uniformity tests allowed the project to achieve this high success rate.

preprint2022arXiv

The $^1$S$_0$ pairing gap in neutron matter

We report ab initio calculations of the S wave pairing gap in neutron matter calculated using realistic nuclear Hamiltonians that include two- and three-body interactions. We use a trial state, properly optimized to capture the essential pairing correlations, from which we extract ground state properties by means of auxiliary field diffusion Monte Carlo simulations. We extrapolate our results to the thermodynamic limit by studying the finite-size effects in the symmetry-restored projected Bardeen-Cooper-Schrieffer (PBCS) theory and compare our results to other ab initio studies done in the past. Our quantum Monte Carlo results for the pairing gap show a modest suppression with respect to the mean-field BCS values. These results can be connected to cold atom experiments, via the unitarity regime where fermionic superfluidity assumes a unified description, and they are important in the prediction of thermal properties and the cooling of neutron stars.

preprint2021arXiv

Fast flavor oscillations in dense neutrino media with collisions

We investigate the impact of the nonzero neutrino splitting and elastic neutrino-nucleon collisions on fast neutrino oscillations. Our calculations confirm that a small neutrino mass splitting and the neutrino mass hierarchy have very little effect on fast oscillation waves. We also demonstrate explicitly that fast oscillations remain largely unaffected for the time/distance scales that are much smaller than the neutrino mean free path but are damped on larger scales. This damping originates from both the direct modification of the dispersion relation of the oscillation waves in the neutrino medium and the flattening of the neutrino angular distributions over time. Our work suggests that fast neutrino oscillation waves produced near the neutrino sphere can propagate essentially unimpeded which may have ramifications in various aspects of the supernova physics.

preprint2021arXiv

Variational approaches to constructing the many-body nuclear ground state for quantum computing

We explore the preparation of specific nuclear states on gate-based quantum hardware using variational algorithms. Large scale classical diagonalization of the nuclear shell model have reached sizes of $10^9 - 10^{10}$ basis states, but are still severely limited by computational resources. Quantum computing can, in principle, solve such systems exactly with exponentially fewer resources than classical computing. Exact solutions for large systems require many qubits and large gate depth, but variational approaches can effectively limit the required gate depth. We use the unitary coupled cluster approach to construct approximations of the ground-state vectors, later to be used in dynamics calculations. The testing ground is the phenomenological shell model space, which allows us to mimic the complexity of the inter-nucleon interactions. We find that often one needs to minimize over a large number of parameters, using a large number of entanglements that makes challenging the application on existing hardware. Prospects for rapid improvements with more capable hardware are, however, very encouraging.

preprint2020arXiv

Ab initio short-range-correlation scaling factors from light to medium-mass nuclei

High-energy scattering processes, such as deep inelastic scattering (DIS) and quasielastic (QE) scattering provide a wealth of information about the structure of atomic nuclei. The remarkable discovery of the empirical linear relationship between the slope of the European Muon Collaboration (EMC) effect in DIS and the short-range-correlation (SRC) scaling factors $a_2$ in QE kinematics is naturally explained in terms of scale separation in effective field theory. This explanation has powerful consequences, allowing us to calculate and predict SRC scaling factors from ab initio low-energy nuclear theory. We present ab initio calculations of SRC scaling factors for a nucleus $A$ relative to the deuteron $a_2(A/d)$ and relative to $^3\rm He$ $a_2(A/^3\rm He)$ in light and medium-mass nuclei. Our framework further predicts that the EMC effect and SRC scaling factors have minimal or negligible isovector corrections.

preprint2020arXiv

Chiral Effective Field Theory Calculations of Weak Transitions in Light Nuclei

We report Quantum Monte Carlo calculations of weak transitions in $A\leq 10$ nuclei, based on the Norfolk two- and three-nucleon chiral interactions, and associated one- and two-body axial currents. We find that the contribution from two-body currents is at the $2$ - $3\%$ level, with the exception of matrix elements entering the rates of $^8$Li, $^8$B, and $^8$He beta decays. These matrix elements are suppressed in impulse approximation based on the (leading order) Gamow Teller transition operator alone; two-body currents provide a $20$ - $30\%$ correction, which is, however, insufficient to bring theory in agreement with experimental data. For the other transitions, the agreement with the data is satisfactory, and the results exhibit a negligible to mild model dependence when different combinations of Norfolk interactions are utilized to construct the nuclear wave functions. We report a complete study of two-body weak transition densities which reveals the expected universal behavior of two-body currents at short distances throughout the range of $A\,$=$\,3$ to $A\,$=$\,10$ systems considered here.

preprint2020arXiv

Clustering of Four-Component Unitary Fermions

Ab initio nuclear physics tackles the problem of strongly interacting four-component fermions. The same setting could foreseeably be probed experimentally in ultracold atomic systems, where two- and three-component experiments have led to major breakthroughs in recent years. Both due to the problem's inherent interest and as a pathway to nuclear physics, in this Letter we study four-component fermions at unitarity via the use of quantum Monte Carlo methods. We explore novel forms of the trial wave function and find one which leads to a ground state of the eight-particle system whose energy is almost equal to that of two four-particle systems. We investigate the clustering properties involved and also extrapolate to the zero-range limit. In addition to being experimentally testable, our results impact the prospects of developing nuclear physics as a perturbation around the unitary limit.

preprint2020arXiv

Nuclear and neutron-star matter from local chiral interactions

We report a quantum Monte Carlo calculation of the equation of state of symmetric nuclear matter using local interactions derived from chiral effective field theory up to next-to-next-to-leading order fit to few-body observables only. The empirical saturation density and energy are well reproduced within statistical and systematic uncertainties. We have also derived the symmetry energy as a function of the density, finding good agreement with available experimentally derived constraints at saturation and twice saturation density. We find that the corresponding pressure is also in excellent agreement with recent constraints extracted from gravitational waves of the neutron-star merger GW170817 by the LIGO-Virgo detection.

preprint2020arXiv

Spectral swaps in a two-dimensional neutrino ring model

Neutrinos emitted deep within a supernova explosion experience a self-induced index of refraction. In the stationary, one-dimensional (1D) supernova "bulb model", this self-induced refraction can lead to a collective flavor transformation which is coherent among different neutrino momentum modes. Such collective oscillations can produce partial swaps of the energy spectra of different neutrino species as the neutrinos stream away from the proto-neutron star. However, it has been demonstrated that the spatial symmetries (such as the spherical symmetry in the bulb model) can be broken spontaneously by collective neutrino oscillations in multi-dimensional models. Using a stationary, 2D neutrino ring model we demonstrate that there exist two limiting scenarios where collective oscillations may occur. In one limit, the collective flavor transformation begins at a radius with relatively high neutrino densities and develops small-scale flavor structures. The loss of the spatial correlation in the neutrino flavor field results in similar (average) energy spectra for the anti-neutrinos of almost all energies and the neutrinos of relatively high energies. In the other limit, the flavor transformation starts at a radius where the neutrino densities are smaller (e.g., due to the suppression of the high matter density near the proto-neutron star). Although the spatial symmetry is broken initially, it is restored as the neutrino densities decrease, and the neutrinos of different flavors partially swap their energy spectra as in the 1D bulb model. This finding may have interesting ramifications in other aspects of supernova physics.

preprint2016arXiv

Chiral Three-Nucleon Interactions in Light Nuclei, Neutron-$α$ Scattering, and Neutron Matter

We present quantum Monte Carlo calculations of light nuclei, neutron-$α$ scattering, and neutron matter using local two- and three-nucleon ($3N$) interactions derived from chiral effective field theory up to next-to-next-to-leading order (N$^2$LO). The two undetermined $3N$ low-energy couplings are fit to the $^4$He binding energy and, for the first time, to the spin-orbit splitting in the neutron-$α$ $P$-wave phase shifts. Furthermore, we investigate different choices of local $3N$-operator structures and find that chiral interactions at N$^2$LO are able to simultaneously reproduce the properties of $A=3,4,5$ systems and of neutron matter, in contrast to commonly used phenomenological $3N$ interactions.

preprint2016arXiv

Electromagnetic response of 12C: a first-principles calculation

The longitudinal and transverse electromagnetic response functions of $^{12}$C are computed in a "first-principles" Green's function Monte Carlo calculation, based on realistic two- and three-nucleon interactions and associated one- and two-body currents. We find excellent agreement between theory and experiment and, in particular, no evidence for the quenching of measured versus calculated longitudinal response. This is further corroborated by a re-analysis of the Coulomb sum rule, in which the contributions from the low-lying $J^π\,$=$\, 2^+$, $0^+_2$ (Hoyle), and $4^+$ states in $^{12}$C are accounted for explicitly in evaluating the total inelastic strength.

preprint2015arXiv

Ab initio calculation of the electromagnetic and neutral-weak response functions of 4He and 12C

Precise measurement of neutrino oscillations, and hence the determination of their masses demands a quantitative understanding of neutrino-nucleus interactions. To this aim, two-body meson-exchange currents have to be accounted for along within realistic models of nuclear dynamics. We summarize our progresses towards the construction of a consistent framework, based on quantum Monte Carlo methods and on the spectral function approach, that can be exploited to accurately describe neutrino interactions with atomic nuclei over the broad kinematical region covered by neutrino experiments.

preprint2015arXiv

Electromagnetic and neutral-weak response functions of 4He and 12C

Ab initio calculations of the quasi-elastic electromagnetic and neutral-weak response functions of 4He and 12C are carried out for the first time. They are based on a realistic approach to nuclear dynamics, in which the strong interactions are described by two- and three-nucleon potentials and the electroweak interactions with external fields include one- and two-body terms. The Green's function Monte Carlo method is used to calculate directly the Laplace transforms of the response functions, and maximum-entropy techniques are employed to invert the resulting imaginary-time correlation functions with associated statistical errors. The theoretical results, confirmed by experiment in the electromagnetic case, show that two-body currents generate excess transverse strength from threshold to the quasi-elastic to the dip region and beyond. These findings challenge the conventional picture of quasi-elastic inclusive scattering as being largely dominated by single-nucleon knockout processes.

preprint2015arXiv

Neutron Matter from Low to High Density

Neutron matter is an intriguing nuclear system with multiple connections to other areas of physics. Considerable progress has been made over the last two decades in exploring the properties of pure neutron fluids. Here we begin by reviewing work done to explore the behavior of very low density neutron matter, which forms a strongly paired superfluid and is thus similar to cold Fermi atoms, though at energy scales differing by many orders of magnitude. We then increase the density, discussing work that ties the study of neutron matter with the determination of the properties of neutron-rich nuclei and neutron-star crusts. After this, we review the impact neutron matter at even higher densities has on the mass-radius relation of neutron stars, thereby making contact with astrophysical observations.

preprint2015arXiv

Quantum Monte Carlo methods for nuclear physics

Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.

preprint2015arXiv

The Intermediate Neutrino Program

The US neutrino community gathered at the Workshop on the Intermediate Neutrino Program (WINP) at Brookhaven National Laboratory February 4-6, 2015 to explore opportunities in neutrino physics over the next five to ten years. Scientists from particle, astroparticle and nuclear physics participated in the workshop. The workshop examined promising opportunities for neutrino physics in the intermediate term, including possible new small to mid-scale experiments, US contributions to large experiments, upgrades to existing experiments, R&D plans and theory. The workshop was organized into two sets of parallel working group sessions, divided by physics topics and technology. Physics working groups covered topics on Sterile Neutrinos, Neutrino Mixing, Neutrino Interactions, Neutrino Properties and Astrophysical Neutrinos. Technology sessions were organized into Theory, Short-Baseline Accelerator Neutrinos, Reactor Neutrinos, Detector R&D and Source, Cyclotron and Meson Decay at Rest sessions.This report summarizes discussion and conclusions from the workshop.

preprint2014arXiv

Exploiting Intrinsic Triangular Geometry in Relativistic He3+Au Collisions to Disentangle Medium Properties

Recent results in d+Au and p+Pb collisions at RHIC and the LHC provide evidence for collective expansion and flow of the created medium. We propose a control set of experiments to directly compare particle emission patterns from p+Au, d+Au, and He3+Au or t+Au collisions at the same sqrt(sNN). Using Monte Carlo Glauber we find that a He3 or triton projectile, with a realistic wavefunction description, induces a significant intrinsic triangular shape to the initial medium and that, even with viscous damping, this survives into a significant third order flow moment v3. By comparing systems with one, two, and three initial hot spots, one can disentangle the effects from the initial spatial distribution of the deposited energy and viscous damping. These are key tools to answering the question of how small a droplet of matter is necessary to form a quark-gluon plasma described by nearly inviscid hydrodynamics.

preprint2014arXiv

From the lightest nuclei to the equation of state of asymmetric nuclear matter with realistic nuclear interactions

We present microscopic calculations of light and medium mass nuclei and the equation of state of symmetric and asymmetric nuclear matter using different nucleon-nucleon forces, including a new Argonne version that has the same spin/isospin structure as local chiral forces at next-to-next-to-leading order (N2LO). The calculations are performed using Auxiliary Field Diffusion Monte Carlo (AFDMC) combined with an improved variational wave function. We show that the AFDMC method can now be used to successfully calculate the energies of very light to medium mass nuclei as well as the energy of isospin-asymmetric nuclear matter, demonstrating microscopically the quadratic dependence of the energy on the symmetry energy.

preprint2014arXiv

Neutral weak current two-body contributions in inclusive scattering from $^{12}$C

An ab initio calculation of the sum rules of the neutral weak response functions in $^{12}$C is reported, based on a realistic Hamiltonian, including two- and three-nucleon potentials, and on realistic currents, consisting of one- and two-body terms. We find that the sum rules of the response functions associated with the longitudinal and transverse components of the (space-like) neutral current are largest and that a significant portion ($\simeq 30$ %) of the calculated strength is due to two-body terms. This fact may have implications for the MiniBooNE and other neutrino quasi-elastic scattering data on nuclei.

preprint2014arXiv

Nuclear Theory and Science of the Facility for Rare Isotope Beams

The Facility for Rare Isotope Beams (FRIB) will be a world-leading laboratory for the study of nuclear structure, reactions and astrophysics. Experiments with intense beams of rare isotopes produced at FRIB will guide us toward a comprehensive description of nuclei, elucidate the origin of the elements in the cosmos, help provide an understanding of matter in neutron stars, and establish the scientific foundation for innovative applications of nuclear science to society. FRIB will be essential for gaining access to key regions of the nuclear chart, where the measured nuclear properties will challenge established concepts, and highlight shortcomings and needed modifications to current theory. Conversely, nuclear theory will play a critical role in providing the intellectual framework for the science at FRIB, and will provide invaluable guidance to FRIB's experimental programs. This article overviews the broad scope of the FRIB theory effort, which reaches beyond the traditional fields of nuclear structure and reactions, and nuclear astrophysics, to explore exciting interdisciplinary boundaries with other areas. \keywords{Nuclear Structure and Reactions. Nuclear Astrophysics. Fundamental Interactions. High Performance Computing. Rare Isotopes. Radioactive Beams.

preprint2014arXiv

Predicting Energies of Small Clusters from the Inhomogeneous Unitary Fermi Gas

We investigate the inhomogeneous unitary Fermi gas and use the long-wavelength properties to predict the energies of small clusters of unitary fermions trapped in harmonic potentials. The large pairing gap and scale invariance place severe restrictions on the form of the density functional. We determine the relevant universal constants needed to constrain the functional from calculations of the bulk in oscillating external potentials. Comparing with exact Quantum Monte Carlo calculations, we find that the same functional correctly predicts the lack of shell closures for small clusters of fermions trapped in harmonic wells as well as their absolute energies. A rapid convergence to the bulk limit in three dimensions, where the surface to volume ratio is quite large, is demonstrated. The resulting functional can be tested experimentally, and is a key ingredient in predicting possible polarized superfluid phases and the properties of the unitary Fermi gas in optical lattices.

preprint2014arXiv

Quantum Monte Carlo Calculations of Light Nuclei Using Chiral Potentials

We present the first Green's function Monte Carlo calculations of light nuclei with nuclear interactions derived from chiral effective field theory up to next-to-next-to-leading order. Up to this order, the interactions can be constructed in a local form and are therefore amenable to quantum Monte Carlo calculations. We demonstrate a systematic improvement with each order for the binding energies of $A=3$ and $A=4$ systems. We also carry out the first few-body tests to study perturbative expansions of chiral potentials at different orders, finding that higher-order corrections are more perturbative for softer interactions. Our results confirm the necessity of a three-body force for correct reproduction of experimental binding energies and radii, and pave the way for studying few- and many-nucleon systems using quantum Monte Carlo methods with chiral interactions.

preprint2013arXiv

Charge form factor and sum rules of electromagnetic response functions in Carbon-12

An "ab initio" calculation of the Carbon-12 elastic form factor, and sum rules of longitudinal and transverse response functions measured in inclusive (e,e') scattering, is reported, based on realistic nuclear potentials and electromagnetic currents. The longitudinal elastic form factor and sum rule are found to be in satisfactory agreement with available experimental data. A direct comparison between theory and experiment is difficult for the transverse sum rule. However, it is shown that the calculated one has large contributions from two-body currents, indicating that these mechanisms lead to a significant enhancement of the quasi-elastic transverse response. This fact may have implications for the anomaly observed in recent neutrino quasi-elastic charge-changing scattering data off Carbon-12.

preprint2013arXiv

Halo Modification of a Supernova Neutronization Neutrino Burst

We give the first self-consistent calculation of the effect of the scattered neutrino halo on flavor evolution in supernovae. Our example case is an O-Ne-Mg core collapse supernova neutronization neutrino burst. We find that the addition of the halo neutrinos produces qualitative and quantitative changes in the final flavor states of neutrinos. We also find that the halo neutrinos produce a novel distortion of the neutrino flavor swap. Our results provide strong motivation for tackling the full multidimensional and composition-dependent aspects of this problem in the future.

preprint2013arXiv

Nucleon and nucleon-pair momentum distributions in A <= 12 nuclei

We report variational Monte Carlo calculations of single-nucleon momentum distributions for A <= 12 nuclei and nucleon-pair and nucleon-cluster momentum distributions for A <= 8. The wave functions have been generated for a Hamiltonian containing the Argonne v18 two-nucleon and Urbana X three-nucleon potentials. The single-nucleon and nucleon-pair momentum distributions exhibit universal features attributable to the one-pion-exchange tensor interaction. The single-nucleon distributions are broken down into proton and neutron components and spin-up and spin-down components where appropriate. The nucleon-pair momentum distributions are given separately for pp and pn pairs. The nucleon-cluster momentum distributions include dp in 3He, tp and dd in 4He, alpha-d in 6Li, alpha-t in 7Li, and alpha-alpha in 8Be. Detailed tables are provided on-line for download.

preprint2013arXiv

Properties of trapped neutrons interacting with realistic nuclear Hamiltonians

We calculate properties of neutron drops in external potentials using both quantum Monte Carlo and no-core full configuration techniques. The properties of the external wells are varied to examine different density profiles. We compare neutron drop results given by a selection of nuclear Hamiltonians, including realistic two-body interactions as well as several three-body forces. We compute a range of properties for the neutron drops: ground-state energies, spin-orbit splittings, excitation energies, radial densities and rms radii. We compare the equations of state for neutron matter for several of these Hamiltonians. Our results can be used as benchmarks to test other many-body techniques, and to constrain properties of energy-density functionals.

preprint2013arXiv

Spin Response and Neutrino Emissivity of Dense Neutron Matter

We study the spin response of cold dense neutron matter in the limit of zero momentum transfer, and show that the frequency dependence of the long-wavelength spin response is well constrained by sum-rules and the asymptotic behavior of the two-particle response at high frequency. The sum-rules are calculated using Auxiliary Field Diffusion Monte Carlo technique and the high frequency two-particle response is calculated for several nucleon-nucleon potentials. At nuclear saturation density, the sum-rules suggest that the strength of the spin response peaks at $ω\simeq$ 40--60 MeV, decays rapidly for $ω\geq $100 MeV, and has a sizable strength below 40 MeV. This strength at relatively low energy may lead to enhanced neutrino production rates in dense neutron-rich matter at temperatures of relevance to core-collapse supernova.

preprint2013arXiv

The equation of state of neutron matter, symmetry energy, and neutron star structure

We review the calculation of the equation of state of pure neutron matter using quantum Monte Carlo (QMC) methods. QMC algorithms permit the study of many-body nuclear systems using realistic two- and three-body forces in a nonperturbative framework. We present the results for the equation of state of neutron matter, and focus on the role of three-neutron forces at supranuclear density. We discuss the correlation between the symmetry energy, the neutron star radius and the symmetry energy. We also combine QMC and theoretical models of the three-nucleon interactions, and recent neutron star observations to constrain the value of the symmetry energy and its density dependence.

preprint2012arXiv

Inclusive neutrino scattering off deuteron from threshold to GeV energies

Background: Neutrino-nucleus quasi-elastic scattering is crucial to interpret the neutrino oscillation results in long baseline neutrino experiments. There are rather large uncertainties in the cross section, due to insufficient knowledge on the role of two-body weak currents. Purpose: Determine the role of two-body weak currents in neutrino-deuteron quasi-elastic scattering up to GeV energies. Methods: Calculate cross sections for inclusive neutrino scattering off deuteron induced by neutral and charge-changing weak currents, from threshold up to GeV energies, using the Argonne $v_{18}$ potential and consistent nuclear electroweak currents with one- and two-body terms. Results: Two-body contributions are found to be small, and increase the cross sections obtained with one-body currents by less than 10% over the whole range of energies. Total cross sections obtained by describing the final two-nucleon states with plane waves differ negligibly, for neutrino energies $\gtrsim 500$ MeV, from those in which interaction effects in these states are fully accounted for. The sensitivity of the calculated cross sections to different models for the two-nucleon potential and/or two-body terms in the weak current is found to be weak. Comparing cross sections to those obtained in a naive model in which the deuteron is taken to consist of a free proton and neutron at rest, nuclear structure effects are illustrated to be non-negligible. Conclusion: Contributions of two-body currents in neutrino-deuteron quasi-elastic scattering up to GeV are found to be smaller than 10%. Finally, it should be stressed that the results reported in this work do not include pion production channels.

preprint2012arXiv

Neutrino scattering and flavor transformation in supernovae

We argue that the small fraction of neutrinos that undergo direction-changing scattering outside of the neutrinosphere could have significant influence on neutrino flavor transformation in core-collapse supernova environments. We show that the standard treatment for collective neutrino flavor transformation is adequate at late times, but could be inadequate in the crucial shock revival/explosion epoch of core-collapse supernovae, where the potentials that govern neutrino flavor evolution are affected by the scattered neutrinos. Taking account of this effect, and the way it couples to entropy and composition, will require a new paradigm in supernova modeling.

preprint2012arXiv

Quantum Monte Carlo approaches to nuclear and atomic physics

Quantum Monte Carlo methods have proven to be valuable in the study of strongly correlated quantum systems, particularly nuclear physics and cold atomic gases. Historically, such ab initio simulations have been used to study properties of light nuclei, including spectra and form factors, low-energy scattering, and high-momentum properties including inclusive scattering and one- and two-body momentum distributions. More recently they have been used to study the properties of homogeneous and inhomogeneous neutron matter and cold atomic gases. There are close analogies between these seemingly diverse systems, including the equation of state, superfluid pairing, and linear response to external probes. In this paper, we compare and contrast results found in nuclear and cold atom physics. We show updated lattice results for the energy of the homogeneous unitary Fermi gas and comparisons with neutron matter, as well as for the dependence of the cold atom energy on the mass ratio between paired particles, which yields insights on the structure of the ground state. We also provide new lattice and continuum results for the harmonically trapped unitary gas, again comparing neutron matter and cold atoms.

preprint2012arXiv

The maximum mass and radius of neutron stars and the nuclear symmetry energy

We calculate the equation of state of neutron matter with realistic two- and three-nucleon interactions using quantum Monte Carlo techniques, and illustrate that the short-range three-neutron interaction determines the correlation between neutron matter energy at nuclear saturation density and higher densities relevant to neutron stars. Our model also makes an experimentally testable prediction for the correlation between the nuclear symmetry energy and its density dependence -- determined solely by the strength of the short-range terms in the three neutron force. The same force provides a significant constraint on the maximum mass and radius of neutron stars.

preprint2011arXiv

Auxiliary Field quantum Monte Carlo for Strongly Paired Fermions

We demonstrate that the inclusion of a BCS importance function dramatically increases the efficiency of the auxiliary field method for strong pairing. We calculate the ground-state energy of an unpolarized fermi gas at unitarity with up to 66 particles and lattices of up to $27^3$ sites. The method has no fermion sign problem, and an accurate result is obtained for the universal parameter $ξ$. Several different forms of the kinetic energy adjusted to the unitary limit but with different effective ranges extrapolate to the same continuum limit within error bars. The finite effective range results for different interactions are consistent with a linear term proportional to the Fermi momentum times the effective range. The new method described herein will have many applications in superfluid cold atom systems and in both electronic and nuclear structures when pairing is important.

preprint2011arXiv

BEC-BCS crossover and universal relations in unitary Fermi gases

The contact parameter in unitary Fermi Gases governs the short-distance, high-momentum, and high-energy properties of the system. We perform accurate quantum Monte Carlo calculations with highly optimized trial functions to precisely determine this parameter at T=0, demonstrate its universal application to a variety of observables, and determine the regions of momentum and energy over which the leading short-range behavior is dominant. We derive Tan's expressions for the contact parameter using just the short-range behavior of the ground-state many-body wave function, and use this behavior to calculate the two-body distribution function, one-body density matrix, and the momentum distribution of unitary Fermi gases; providing a precise value of the contact parameter that can be compared to experiments.

preprint2011arXiv

Cold neutrons trapped in external fields

The properties of inhomogeneous neutron matter are crucial to the physics of neutron-rich nuclei and the crust of neutron stars. Advances in computational techniques now allow us to accurately determine the binding energies and densities of many neutrons interacting via realistic microscopic interactions and confined in external fields. We perform calculations for different external fields and across several shells to place important constraints on inhomogeneous neutron matter, and hence the large isospin limit of the nuclear energy density functionals that are used to predict properties of heavy nuclei and neutron star crusts. We find important differences between microscopic calculations and current density functionals; in particular the isovector gradient terms are significantly more repulsive than in traditional models, and the spin-orbit and pairing forces are comparatively weaker.

preprint2011arXiv

Jastrow functions in double-beta decay

We use simple analytic considerations and a Monte Carlo calculation of nucleons in a box to argue that the use of Jastrow functions as short-range correlators in the commonly employed two-body-cluster approximation causes significant errors in the matrix elements for double-beta decay. The Jastrow approach appears to agree with others, however, if many-body clusters are included. A careful treatment of the charge-changing analog of the nuclear pair density shows, in addition, that differences between Unitary Correlator Operator Method and Brueckner methods for treating short-range correlations in double-beta are less significant than suggested by previous work.

preprint2011arXiv

Terrestrial and Astrophysical Superfluidity: Cold Atoms and Neutron Matter

After a brief historical overview of superfluidity in connection with neutron matter and cold fermionic atoms, we discuss the commonalities between these two systems as well as their relevance to the physics of neutron star crusts. We then present the methodological tools we use to attack the many-body problem, comparing Quantum Monte Carlo to mean-field theory and to the analytic expectations at weak coupling. We review recent results on the equation of state and the pairing gap of cold atoms and neutron matter and contrast them with a variety of calculations by other groups. We conclude by giving a few directions of possible future inquiry.

preprint2010arXiv

Heavy-Light Few Fermion Clusters at Unitarity

We examine the physics of two, three, and four heavy fermions interacting with a single light fermion via short-range interactions. Four-particle bosonic Efimov states have proven important experimentally and also been the subject of significant theoretical effort. Similar fermionic systems are just now being investigated. We find that with some simple interactions the four- and five-particle states collapse to the interaction range at smaller mass ratios than the three-body state, and also before larger clusters can collapse. These states and their excitations can be studied in cold atom experiments, providing unique insights into the role of few-body systems in many-body physics.

preprint2010arXiv

Low-density neutron matter

The properties of low-density neutron matter are important for the understanding of neutron star crusts and the exterior of large neutron-rich nuclei. We examine various properties of dilute neutron matter using quantum Monte Carlo methods, with s- and p-wave terms in the interaction. Our results provide a smooth evolution of the equation of state and pairing gap from extremely small densities, where analytic expressions are available, up to the strongly interacting regime probed experimentally and described theoretically in cold atomic systems, where the Fermi momentum is approximately half an inverse fermi and the pairing gap becomes of the order of magnitude of the Fermi energy. We also present results for the momentum distribution and pair distributions, displaying the same evolution from weak to strong coupling. Combined with previous quantum Monte Carlo and other calculations at moderate densities, these results provide strong constraints on the neutron matter equation of state up to saturation densities.

preprint2004arXiv

Neutron Matter: A Superfluid Gas

We review recent progress in the theory of neutron matter with particular emphasis on its superfluid properties. Results of quantum Monte Carlo calculations of simple and realistic models of uniform superfluid neutron gas are discussed along with those of neutrons interacting in a potential well chosen to approximate neutron-rich oxygen isotopes. The properties of dilute superfluid Fermi gases that may be produced in atom traps, and their relations with neutron matter, are illustrated. The density dependence of the effective interaction between neutrons, used to describe neutron-rich systems with the mean field approximations, is also discussed.

preprint2004arXiv

Parity-Violating Interaction Effects in the np System

We investigate parity-violating observables in the np system, including the longitudinal asymmetry and neutron-spin rotation in np elastic scattering, the photon asymmetry in np radiative capture, and the asymmetries in deuteron photo-disintegration d(gamma,n)p in the threshold region and electro-disintegration d(e,e`)np in quasi-elastic kinematics. To have an estimate of the model dependence for the various predictions, a number of different, latest-generation strong-interaction potentials--Argonne v18, Bonn 2000, and Nijmegen I--are used in combination with a weak-interaction potential consisting of pi-, rho-, and omega-meson exchanges--the model known as DDH. The complete bound and scattering problems in the presence of parity-conserving, including electromagnetic, and parity-violating potentials is solved in both configuration and momentum space. The issue of electromagnetic current conservation is examined carefully. We find large cancellations between the asymmetries induced by the parity-violating interactions and those arising from the associated pion-exchange currents. In the np capture, the model dependence is nevertheless quite small, because of constraints arising through the Siegert evaluation of the relevant E1 matrix elements. In quasi-elastic electron scattering these processes are found to be insignificant compared to the asymmetry produced by gamma-Z interference on individual nucleons.

preprint2002arXiv

Parity Violating Interactions and Currents in the Deuteron

We investigate parity-violating asymmetries in polarized n p radiative capture, and deuteron electro-disintegration in quasi-elastic kinematics, using the DDH model for the parity-violating nucleon-nucleon interaction. We find dramatic cancellations between the asymmetries induced by the parity-violating interaction and those arising from the associated parity-violating pion-exchange currents. In np capture, the model-dependence of the result is nevertheless quite small because of constraints arising through the Siegert evaluation of the relevant E1 matrix element. In quasi-elastic electron scattering these processes are found to be insignificant comared to the asymmetry produced by the gamma-Z interference on individual nucleons. These two experiments, then, provide clean probes of different aspects of weak-interaction physics associated with parity violation in the np system.

preprint2001arXiv

Benchmark Test Calculation of a Four-Nucleon Bound State

In the past, several efficient methods have been developed to solve the Schroedinger equation for four-nucleon bound states accurately. These are the Faddeev-Yakubovsky, the coupled-rearrangement-channel Gaussian-basis variational, the stochastic variational, the hyperspherical variational, the Green's function Monte Carlo, the no-core shell model and the effective interaction hyperspherical harmonic methods. In this article we compare the energy eigenvalue results and some wave function properties using the realistic AV8' NN interaction. The results of all schemes agree very well showing the high accuracy of our present ability to calculate the four-nucleon bound state.

preprint2001arXiv

Longitudinal and Transverse Quasi-Elastic Response Functions of Light Nuclei

The $^3$He and $^4$He longitudinal and transverse response functions are determined from an analysis of the world data on quasi-elastic inclusive electron scattering. The corresponding Euclidean response functions are derived and compared to those calculated with Green's function Monte Carlo methods, using realistic interactions and currents. Large contributions associated with two-body currents are found, particularly in the $^4$He transverse response, in agreement with data. The contributions of two-body charge and current operators in the $^3$He, $^4$He, and $^6$Li response functions are also studied via sum-rule techniques. A semi-quantitative explanation for the observed systematics in the excess of transverse quasi-elastic strength, as function of mass number and momentum transfer, is provided. Finally, a number of model studies with simplified interactions, currents, and wave functions is carried out to elucidate the role played, in the full calculation, by tensor interactions and correlations.

preprint2001arXiv

Parity-Violating Interaction Effects I: the Longitudinal Asymmetry in pp Elastic Scattering

The proton-proton parity-violating longitudinal asymmetry is calculated in the lab-energy range 0--350 MeV, using a number of different, latest-generation strong-interaction potentials--Argonne V18, Bonn-2000, and Nijmegen-I--in combination with a weak-interaction potential consisting of rho- and omega-meson exchanges--the model known as DDH. The complete scattering problem in the presence of parity-conserving, including Coulomb, and parity-violating potentials is solved in both configuration- and momentum-space. The predicted parity-violating asymmetries are found to be only weakly dependent upon the input strong-interaction potential adopted in the calculation. Values for the rho- and omega-meson weak coupling constants $h^{pp}_ρ$ and $h^{pp}_ω$ are determined by reproducing the measured asymmetries at 13.6 MeV, 45 MeV, and 221 MeV.

preprint1995arXiv

A Constrained Path Quantum Monte Carlo Method for Fermion Ground States

We propose a new quantum Monte Carlo algorithm to compute fermion ground-state properties. The ground state is projected from an initial wavefunction by a branching random walk in an over-complete basis space of Slater determinants. By constraining the determinants according to a trial wavefunction $|Ψ_T \rangle$, we remove the exponential decay of signal-to-noise ratio characteristic of the sign problem. The method is variational and is exact if $|Ψ_T\rangle$ is exact. We report results on the two-dimensional Hubbard model up to size $16\times 16$, for various electron fillings and interaction strengths.

preprint1995arXiv

Quantum Monte Carlo Calculations of $A\leq6$ Nuclei

The energies of $^{3}H$, $^{3}He$, and $^{4}He$ ground states, the ${\frac{3}{2}}^{-}$ and ${\frac{1}{2}}^{-}$ scattering states of $^{5}He$, the ground states of $^{6}He$, $^{6}Li$, and $^{6}Be$ and the $3^{+}$ and $0^{+}$ excited states of $^{6}Li$ have been accurately calculated with the Green's function Monte Carlo method using realistic models of two- and three-nucleon interactions. The splitting of the $A=3$ isospin $T=\frac{1}{2}$ and $A=6$ isospin $T=1$, $J^π = 0^{+}$ multiplets is also studied. The observed energies and radii are generally well reproduced, however, some definite differences between theory and experiment can be identified.