Source author record

Dean Lee

Dean Lee 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

58works
23topics
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

58 published item(s)

preprint2026arXiv

The MASK Benchmark: Disentangling Honesty From Accuracy in AI Systems

As large language models (LLMs) become more capable and agentic, the requirement for trust in their outputs grows significantly, yet at the same time concerns have been mounting that models may learn to lie in pursuit of their goals. To address these concerns, a body of work has emerged around the notion of "honesty" in LLMs, along with interventions aimed at mitigating deceptive behaviors. However, some benchmarks claiming to measure honesty in fact simply measure accuracy--the correctness of a model's beliefs--in disguise. Moreover, no benchmarks currently exist for directly measuring whether language models lie. In this work, we introduce a large-scale human-collected dataset for directly measuring lying, allowing us to disentangle accuracy from honesty. Across a diverse set of LLMs, we find that while larger models obtain higher accuracy on our benchmark, they do not become more honest. Surprisingly, most frontier LLMs obtain high scores on truthfulness benchmarks yet exhibit a substantial propensity to lie under pressure, resulting in low honesty scores on our benchmark. We find that simple methods, such as representation engineering interventions, can improve honesty. These results underscore the growing need for robust evaluations and effective interventions to ensure LLMs remain trustworthy.

preprint2022arXiv

Nuclear Forces for Precision Nuclear Physics -- a collection of perspectives

This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges ahead, and new ideas and strategies to make progress in nuclear structure and reaction physics, effective field theory, lattice QCD, quantum information, and quantum computing. The contributed pieces solely reflect the perspectives of the respective authors and do not represent the viewpoints of the Institute for Nuclear theory or the organizers of the program.

preprint2022arXiv

Perturbative quantum Monte Carlo method for nuclear physics

While first order perturbation theory is routinely used in quantum Monte Carlo (QMC) calculations, higher-order terms present significant numerical challenges. We present a new approach for computing perturbative corrections in projection QMC calculations. We demonstrate the method by computing nuclear ground state energies up to second order for a realistic chiral interaction. We calculate the binding energies of several light nuclei up to $^{16}$O by expanding the Hamiltonian around the Wigner SU(4) limit and find good agreement with data. In contrast to the natural ordering of the perturbative series, we find remarkably large second order energy corrections. This occurs because the perturbing interactions break the symmetries of the unperturbed Hamiltonian. Our method is free from the sign problem and can be applied to QMC calculations for many-body systems in nuclear physics, condensed matter physics, ultracold atoms, and quantum chemistry.

preprint2022arXiv

Self-learning Emulators and Eigenvector Continuation

Emulators that can bypass computationally expensive scientific calculations with high accuracy and speed can enable new studies of fundamental science as well as more potential applications. In this work we discuss solving a system of constraint equations efficiently using a self-learning emulator. A self-learning emulator is an active learning protocol that can be used with any emulator that faithfully reproduces the exact solution at selected training points. The key ingredient is a fast estimate of the emulator error that becomes progressively more accurate as the emulator is improved, and the accuracy of the error estimate can be corrected using machine learning. We illustrate with three examples. The first uses cubic spline interpolation to find the solution of a transcendental equation with variable coefficients. The second example compares a spline emulator and a reduced basis method emulator to find solutions of a parameterized differential equation. The third example uses eigenvector continuation to find the eigenvectors and eigenvalues of a large Hamiltonian matrix that depends on several control parameters.

preprint2021arXiv

Alpha-alpha scattering in the Multiverse

We investigate the phase shifts of low-energy alpha-alpha scattering under variations of the fundamental parameters of the Standard Model, namely the light quark mass, the electromagnetic fine-structure constant as well as the QCD theta-angle. As a first step, we recalculate alpha-alpha scattering in our Universe utilizing various improvements in the adiabatic projection method, which leads to an improved, parameter-free prediction of the S- and D-wave phase shifts for laboratory energies below 10 MeV. We find that positive shifts in the pion mass have a small effect on the S-wave phase shift, whereas lowering the pion mass adds some repulsion in the two-alpha system. The effect on the D-wave phase shift turns out to be more pronounced as signaled by the D-wave resonance parameters. Variations of the fine-structure constant have almost no effect on the low-energy alpha-alpha phase shifts. We further show that up-to-and-including next-to-leading order in the chiral expansion, variations of these phase shifts with respect to the QCD theta-angle can be expressed in terms of the theta-dependent pion mass.

preprint2020arXiv

$θ$-dependence of light nuclei and nucleosynthesis

We investigate the impact of the QCD vacuum at nonzero $θ$ on the properties of light nuclei, Big Bang nucleosynthesis, and stellar nucleosynthesis. Our analysis starts with a calculation of the $θ$-dependence of the neutron-proton mass difference and neutron decay using chiral perturbation theory. We then discuss the $θ$-dependence of the nucleon-nucleon interaction using a one-boson-exchange model and compute the properties of the two-nucleon system. Using the universal properties of four-component fermions at large scattering length, we then deduce the binding energies of the three-nucleon and four-nucleon systems. Based on these results, we discuss the implications for primordial abundances of light nuclei, the production of nuclei in stellar environments, and implications for an anthropic view of the universe.

preprint2020arXiv

Convergence of Eigenvector Continuation

Eigenvector continuation is a computational method that finds the extremal eigenvalues and eigenvectors of a Hamiltonian matrix with one or more control parameters. It does this by projection onto a subspace of eigenvectors corresponding to selected training values of the control parameters. The method has proven to be very efficient and accurate for interpolating and extrapolating eigenvectors. However, almost nothing is known about how the method converges, and its rapid convergence properties have remained mysterious. In this letter we present the first study of the convergence of eigenvector continuation. In order to perform the mathematical analysis, we introduce a new variant of eigenvector continuation that we call vector continuation. We first prove that eigenvector continuation and vector continuation have identical convergence properties and then analyze the convergence of vector continuation. Our analysis shows that, in general, eigenvector continuation converges more rapidly than perturbation theory. The faster convergence is achieved by eliminating a phenomenon that we call differential folding, the interference between non-orthogonal vectors appearing at different orders in perturbation theory. From our analysis we can predict how eigenvector continuation converges both inside and outside the radius of convergence of perturbation theory. While eigenvector continuation is a non-perturbative method, we show that its rate of convergence can be deduced from power series expansions of the eigenvectors. Our results also yield new insights into the nature of divergences in perturbation theory.

preprint2020arXiv

Effective interactions between nuclear clusters

The effective interactions between two nuclear clusters, $d+d$, $t+t$, and $α+α$, are investigated within a cluster model using local nucleon-nucleon~($NN$) forces. It is shown that the interaction in the spin-aligned $d+d$ system is repulsive for all inter-cluster distances, whereas the $α+α$ and spin-aligned $t+t$ systems are attractive at intermediate distances. The Pauli blocking between identical-nucleon pairs is responsible for the cluster-cluster repulsion and becomes dominant in the shallow binding limit. We demonstrate that two $d$-clusters could be bound if the $NN$ force has nonzero range and is strong enough to form a deeply bound $d$-cluster, or if the $NN$ force has both even-parity and odd-parity attraction. Effective dimer-dimer interactions for general quantum systems of two-component fermions are also discussed in heavy-light mass limit, where one component is much heavier than the other, and their relation to inter-cluster interactions in nuclear systems are discussed. Our findings provide a conceptual foundation for conclusions obtained numerically in the literature, that increasing the range or strength of the local part of the attractive nucleon-nucleon interaction results in a more attractive cluster-cluster interaction.

preprint2020arXiv

Projected Cooling Algorithm for Quantum Computation

In the current era of noisy quantum devices, there is a need for quantum algorithms that are efficient and robust against noise. Towards this end, we introduce the projected cooling algorithm for quantum computation. The projected cooling algorithm is able to construct the localized ground state of any Hamiltonian with a translationally-invariant kinetic energy and interactions that vanish at large distances. The term "localized" refers to localization in position space. The method can be viewed as the quantum analog of evaporative cooling. We start with an initial state with support over a compact region of a large volume. We then drive the excited quantum states to disperse and measure the remaining portion of the wave function left behind. For the nontrivial examples we consider here, the improvement over other methods is substantial. The only additional resource required is performing the operations in a volume significantly larger than the size of the localized state. These characteristics make the projected cooling algorithm a promising tool for calculations of self-bound systems such as atomic nuclei.

preprint2020arXiv

Report from the A.I. For Nuclear Physics Workshop

This report is an outcome of the workshop "AI for Nuclear Physics" held at Thomas Jefferson National Accelerator Facility on March 4-6, 2020. The workshop brought together 184 scientists to explore opportunities for Nuclear Physics in the area of Artificial Intelligence. The workshop consisted of plenary talks, as well as six working groups. The report includes the workshop deliberations and additional contributions to describe prospects for using AI across Nuclear Physics research.

preprint2020arXiv

Superfluid Condensate Fraction and Pairing Wave Function of the Unitary Fermi Gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gases and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this work we present auxiliary-field lattice Monte Carlo simulations using a novel lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test we compute the ground state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{FG}$, we find $E_0/E_{FG}= 0.369(2), 0.372(2)$, using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is $α= 0.43(2)$. We also extract the pairing wave function and find the pair correlation length to be $ζ_pk_F = 1.8(3) \hbar$, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic P-wave and D-wave superfluids.

preprint2019arXiv

Improved many-body expansions from eigenvector continuation

Quantum many-body theory has witnessed tremendous progress in various fields, ranging from atomic and solid-state physics to quantum chemistry and nuclear structure. Due to the inherent computational burden linked to the ab initio treatment of microscopic fermionic systems, it is desirable to obtain accurate results through low-order perturbation theory. In atomic nuclei however, effects such as strong short-range repulsion between nucleons can spoil the convergence of the expansion and make the reliability of perturbation theory unclear. Mathematicians have devised an extensive machinery to overcome the problem of divergent expansions by making use of so-called resummation methods. In large-scale many-body applications such schemes are often of limited use since no a priori analytical knowledge of the expansion is available. We present here eigenvector continuation as an alternative resummation tool that is both efficient and reliable because it is based on robust and simple mathematical principles.

preprint2016arXiv

Nuclear binding near a quantum phase transition

How do protons and neutrons bind to form nuclei? This is the central question of ab initio nuclear structure theory. While the answer may seem as simple as the fact that nuclear forces are attractive, the full story is more complex and interesting. In this work we present numerical evidence from ab initio lattice simulations showing that nature is near a quantum phase transition, a zero-temperature transition driven by quantum fluctuations. Using lattice effective field theory, we perform Monte Carlo simulations for systems with up to twenty nucleons. For even and equal numbers of protons and neutrons, we discover a first-order transition at zero temperature from a Bose-condensed gas of alpha particles (4He nuclei) to a nuclear liquid. Whether one has an alpha-particle gas or nuclear liquid is determined by the strength of the alpha-alpha interactions, and we show that the alpha-alpha interactions depend on the strength and locality of the nucleon-nucleon interactions. This insight should be useful in improving calculations of nuclear structure and important astrophysical reactions involving alpha capture on nuclei. Our findings also provide a tool to probe the structure of alpha cluster states such as the Hoyle state responsible for the production of carbon in red giant stars and point to a connection between nuclear states and the universal physics of bosons at large scattering length.

preprint2016arXiv

Nucleon-deuteron scattering using the adiabatic projection method

In this paper we discuss the adiabatic projection method, a general framework for scattering and reaction calculations on the lattice. We also introduce several new techniques developed to study nucleus-nucleus scattering and reactions on the lattice. We present technical details of the method for large-scale problems. To estimate the systematic errors of the calculations we consider simple two-particle scattering on the lattice. Then we benchmark the accuracy and efficiency of the numerical methods by applying these to calculate fermion-dimer scattering in lattice effective field theory with and without a long-range Coulomb potential. The fermion-dimer calculations correspond to neutron-deuteron and proton-deuteron scattering in the spin-quartet channel at leading order in the pionless effective field theory.

preprint2016arXiv

Precise determination of lattice phase shifts and mixing angles

We introduce a general and accurate method for determining lattice phase shifts and mixing angles, which is applicable to arbitrary, non-cubic lattices. Our method combines angular momentum projection, spherical wall boundaries and an adjustable auxiliary potential. This allows us to construct radial lattice wave functions and to determine phase shifts at arbitrary energies. For coupled partial waves, we use a complex-valued auxiliary potential that breaks time-reversal invariance. We benchmark our method using a system of two spin-1/2 particles interacting through a finite-range potential with a strong tensor component. We are able to extract phase shifts and mixing angles for all angular momenta and energies, with precision greater than that of extant methods. We discuss a wide range of applications from nuclear lattice simulations to optical lattice experiments.

preprint2015arXiv

Ab initio alpha-alpha scattering

Processes involving alpha particles and alpha-like nuclei comprise a major part of stellar nucleosynthesis and hypothesized mechanisms for thermonuclear supernovae. In an effort towards understanding alpha processes from first principles, we describe in this letter the first ab initio calculation of alpha-alpha scattering. We use lattice effective field theory to describe the low-energy interactions of nucleons and apply a technique called the adiabatic projection method to reduce the eight-body system to an effective two-cluster system. We find good agreement between lattice results and experimental phase shifts for S-wave and D-wave scattering. The computational scaling with particle number suggests that alpha processes involving heavier nuclei are also within reach in the near future.

preprint2015arXiv

Ab initio lattice results for Fermi polarons in two dimensions

We investigate the attractive Fermi polaron problem in two dimensions using non-perturbative Monte Carlo simulations. We introduce a new Monte Carlo algorithm called the impurity lattice Monte Carlo method. This algorithm samples the path integral in a computationally efficient manner and has only small sign oscillations for systems with a single impurity. As a benchmark of the method, we calculate the universal polaron energy in three dimensions in the scale-invariant unitarity limit and find agreement with published results. We then present the first fully non-perturbative calculations of the polaron energy in two dimensions and density correlations between the impurity and majority particles in the limit of zero range interactions. We find evidence for a smooth crossover transition from fermionic quasiparticle to molecular state as a function of interaction strength.

preprint2015arXiv

Breaking and restoration of rotational symmetry for irreducible tensor operators on the lattice

We study the breaking of rotational symmetry on the lattice for irreducible tensor operators and practical methods for suppressing this breaking. We illustrate the features of the general problem using an $α$ cluster model for $^{8}$Be. We focus on the lowest states with non-zero angular momentum and examine the matrix elements of multipole moment operators. We show that the physical reduced matrix element is well reproduced by averaging over all possible orientations of the quantum state, and this is expressed as a sum of matrix elements weighted by the corresponding Clebsch-Gordan coefficients. For our $α$ cluster model we find that the effects of rotational symmetry breaking can be largely eliminated for lattice spacings of $a\leq 1.7$ fm, and we expect similar improvement for actual lattice Monte Carlo calculations.

preprint2015arXiv

Nuclear Lattice Simulations using Symmetry-Sign Extrapolation

Projection Monte Carlo calculations of lattice Chiral Effective Field Theory suffer from sign oscillations to a varying degree dependent on the number of protons and neutrons. Hence, such studies have hitherto been concentrated on nuclei with equal numbers of protons and neutrons, and especially on the alpha nuclei where the sign oscillations are smallest. Here, we introduce the "symmetry-sign extrapolation" method, which allows us to use the approximate Wigner SU(4) symmetry of the nuclear interaction to systematically extend the Projection Monte Carlo calculations to nuclear systems where the sign problem is severe. We benchmark this method by calculating the ground-state energies of the $^{12}$C, $^6$He and $^6$Be nuclei, and discuss its potential for studies of neutron-rich halo nuclei and asymmetric nuclear matter.

preprint2015arXiv

Regularization Methods for Nuclear Lattice Effective Field Theory

We investigate Nuclear Lattice Effective Field Theory for the two-body system for several lattice spacings at lowest order in the pionless as well as in the pionful theory. We discuss issues of regularizations and predictions for the effective range expansion. In the pionless case, a simple Gaussian smearing allows to demonstrate lattice spacing independence over a wide range of lattice spacings. We show that regularization methods known from the continuum formulation are necessary as well as feasible for the pionful approach.

preprint2015arXiv

Scattering cluster wave functions on the lattice using the adiabatic projection method

The adiabatic projection method is a general framework for studying scattering and reactions on the lattice. It provides a low-energy effective theory for clusters which becomes exact in the limit of large Euclidean projection time. Previous studies have used the adiabatic projection method to extract scattering phase shifts from finite periodic-box energy levels using Lüschers method. In this paper we demonstrate that scattering observables can be computed directly from asymptotic cluster wave functions. For a variety of examples in one and three spatial dimensions, we extract elastic phase shifts from asymptotic cluster standing waves corresponding to spherical wall boundary conditions. We find that this approach of extracting scattering wave functions from the adiabatic Hamiltonian to be less sensitive to small stochastic and systematic errors as compared with using periodic-box energy levels.

preprint2014arXiv

Ab initio calculation of the spectrum and structure of $^{16}$O

We present ab initio lattice calculations of the low-energy even-parity states of $^{16}$O using chiral nuclear effective field theory. We find good agreement with the empirical energy spectrum, and with the electromagnetic properties and transition rates. For the ground state, we find that the nucleons are arranged in a tetrahedral configuration of alpha clusters. For the first excited spin-0 state, we find that the predominant structure is a square configuration of alpha clusters, with rotational excitations that include the first spin-2 state.

preprint2014arXiv

Breaking and restoration of rotational symmetry on the lattice for bound state multiplets

We explore the breaking of rotational symmetry on the lattice for bound state energies and practical methods for suppressing this breaking. We demonstrate the general problems associated with lattice discretization errors and finite-volume errors using an $α$ cluster model for $^8$Be and $^{12}$C. We consider the two and three $α$-particle systems and focus on the lowest states with non-zero angular momentum which split into multiplets corresponding to different irreducible representations of the cubic group. We examine the dependence of such splittings on the lattice spacing and box size. We find that lattice spacing errors are closely related to the commensurability of the lattice with the intrinsic length scales of the system. We also show that rotational symmetry breaking effects can be significantly reduced by using improved lattice actions, and that the physical energy levels are accurately reproduced by the weighted average of a given spin multiplets.

preprint2014arXiv

Fermion-dimer scattering using an impurity lattice Monte Carlo approach and the adiabatic projection method

We present lattice Monte Carlo calculations of fermion-dimer scattering in the limit of zero-range interactions using the adiabatic projection method. The adiabatic projection method uses a set of initial cluster states and Euclidean time projection to give a systematically improvable description of the low-lying scattering cluster states in a finite volume. We use Lüscher's finite-volume relations to determine the s-wave, p-wave, and d-wave phase shifts. For comparison, we also compute exact lattice results using Lanczos iteration and continuum results using the Skorniakov-Ter-Martirosian equation. For our Monte Carlo calculations we use a new lattice algorithm called impurity lattice Monte Carlo. This algorithm can be viewed as a hybrid technique which incorporates elements of both worldline and auxiliary-field Monte Carlo simulations.

preprint2014arXiv

Lattice Effective Field Theory for Medium-Mass Nuclei

We extend Nuclear Lattice Effective Field Theory (NLEFT) to medium-mass nuclei, and present results for the ground states of alpha nuclei from $^4$He to $^{28}$Si, calculated up to next-to-next-to-leading order (NNLO) in the EFT expansion. This computational advance is made possible by extrapolations of lattice data using multiple initial and final states. For our soft two-nucleon interaction, we find that the overall contribution from multi-nucleon forces must change sign from attractive to repulsive with increasing nucleon number. This effect is not produced by three-nucleon forces at NNLO, but it can be approximated by an effective four-nucleon interaction. We discuss the convergence of the EFT expansion and the broad significance of our findings for future ab initio calculations.

preprint2014arXiv

The Hoyle state in Nuclear Lattice EFT

We review the calculation of the Hoyle state of $^{12}$C in Nuclear Lattice Effective Field Theory (NLEFT) and its anthropic implications for the nucleosynthesis of $^{12}$C and $^{16}$O in red giant stars. We also review the extension of NLEFT to the regime of medium-mass nuclei, with emphasis on the determination of the ground-state energies of the alpha nuclei $^{16}$O, $^{20}$Ne, $^{24}$Mg and $^{28}$Si by means of Euclidean time projection. Finally, we review recent NLEFT results for the spectrum, electromagnetic properties, and alpha-cluster structure of $^{16}$O.

preprint2014arXiv

Uncertainties of Euclidean Time Extrapolation in Lattice Effective Field Theory

Extrapolations in Euclidean time form a central part of Nuclear Lattice Effective Field Theory (NLEFT) calculations using the Projection Monte Carlo method, as the sign problem in many cases prevents simulations at large Euclidean time. We review the next-to-next-to-leading order NLEFT results for the alpha nuclei up to $^{28}$Si, with emphasis on the Euclidean time extrapolations, their expected accuracy and potential pitfalls. We also discuss possible avenues for improving the reliability of Euclidean time extrapolations in NLEFT.

preprint2013arXiv

Adiabatic projection method for scattering and reactions on the lattice

We demonstrate and test the adiabatic projection method, a general new framework for calculating scattering and reactions on the lattice. The method is based upon calculating a low-energy effective theory for clusters which becomes exact in the limit of large Euclidean projection time. As a detailed example we calculate the adiabatic two-body Hamiltonian for elastic fermion-dimer scattering in lattice effective field theory. Our calculation corresponds to neutron-deuteron scattering in the spin-quartet channel at leading order in pionless effective field theory. We show that the spectrum of the adiabatic Hamiltonian reproduces the spectrum of the original Hamiltonian below the inelastic threshold to arbitrary accuracy. We also show that the calculated s-wave phase shift reproduces the known exact result in the continuum and infinite-volume limits. When extended to more than one scattering channel, the adiabatic projection method can be used to calculate inelastic reactions on the lattice in future work.

preprint2013arXiv

Causality constraints for charged particles

In quantum systems with short-range interactions, causality imposes nontrivial constraints on low-energy scattering parameters. We investigate these causality constraints for systems where a long-range Coulomb potential is present in addition to a short-range interaction. The main result is an upper bound for the Coulomb-modified effective range parameter. We discuss the implications of this bound to the effective field theory for nuclear halo systems. In particular, we consider several examples of proton--nucleus and nucleus--nucleus scattering. For the bound-state regime, we find relations for the asymptotic normalization coefficients (ANCs) of nuclear halo states. As an application of these relations, we extract the ANCs of the excited 2+ and 1- states in O-16 from alpha--C-12 scattering data.

preprint2013arXiv

Causality, universality, and effective field theory for van der Waals interactions

We analyze low-energy scattering for arbitrary short-range interactions plus an attractive 1/r^6 tail. We derive the constraints of causality and unitarity and find that the van der Waals length scale dominates over parameters characterizing the short-distance physics of the interaction. This separation of scales suggests a separate universality class for physics characterizing interactions with an attractive 1/r^6 tail. We argue that a similar universality class exists for any attractive potential 1/r^{alpha} for alpha >= 2. We also discuss the extension to multi-channel systems near a magnetic Feshbach resonance. We discuss the implications for effective field theory with attractive singular power law tails.

preprint2013arXiv

Dependence of the triple-alpha process on the fundamental constants of nature

We present an ab initio calculation of the quark mass dependence of the ground state energies of ^4He, ^8Be and ^{12}C, and of the energy of the Hoyle state in ^{12}C. These investigations are performed within the framework of lattice chiral Effective Field Theory. We address the sensitivity of the production rate of carbon and oxygen in red giant stars to the fundamental constants of nature by considering the impact of variations in the light quark masses and the electromagnetic fine-structure constant on the reaction rate of the triple-alpha process. As carbon and oxygen are essential to life as we know it, we also discuss the implications of our findings for an anthropic view of the Universe. We find strong evidence that the physics of the triple-alpha process is driven by alpha clustering, and that shifts in the fundamental parameters at the \simeq 2 - 3 % level are unlikely to be detrimental to the development of life. Tolerance against much larger changes cannot be ruled out at present, given the relatively limited knowledge of the quark mass dependence of the two-nucleon S-wave scattering parameters. Lattice QCD is expected to provide refined estimates of the scattering parameters in the future.

preprint2013arXiv

Effective Field Theory for Bound State Reflection

Elastic quantum bound-state reflection from a hard-wall boundary provides direct information regarding the structure and compressibility of quantum bound states. We discuss elastic quantum bound-state reflection and derive a general theory for elastic reflection of shallow dimers from hard-wall surfaces using effective field theory. We show that there is a small expansion parameter for analytic calculations of the reflection scattering length. We present a calculation up to second order in the effective Hamiltonian in one, two, and three dimensions. We also provide numerical lattice results for all three cases as a comparison with our effective field theory results. Finally, we provide an analysis of the compressibility of the alpha particle confined to a cubic lattice with vanishing Dirichlet boundaries.

preprint2013arXiv

Finite volume effects in low-energy neutron-deuteron scattering

We present a lattice calculation of neutron-deuteron scattering at very low energies and investigate in detail the impact of the topological finite-volume corrections. Our calculations are carried out in the framework of pionless effective field theory at leading order in the low-energy expansion. Using lattice sizes and a lattice spacing comparable to those employed in nuclear lattice simulations, we find that the topological volume corrections must be taken into account in order to obtain correct results for the neutron-proton S-wave scattering lengths.

preprint2013arXiv

Radiative capture reactions in lattice effective field theory

We outline a general method for computing nuclear capture reactions on the lattice. The method consists of two major parts. In this study we detail the second part which consists of calculating an effective two-body capture reaction on the lattice at finite volume. We solve this problem by calculating the two-point Green's function using an infrared regulator and the capture amplitude to a two-body bound state. We demonstrate the details of this method by calculating on the lattice the leading M1 contribution to the radiative neutron capture on proton at low energies using pionless effective field theory. We find good agreement with exact continuum results.

preprint2013arXiv

Viability of carbon-based life as a function of the light quark mass

The Hoyle state plays a crucial role in the helium burning of stars that have reached the red giant stage. The close proximity of this state to the triple-alpha threshold is needed for the production of carbon, oxygen, and other elements necessary for life. We investigate whether this life-essential condition is robust or delicately fine-tuned by measuring its dependence on the fundamental constants of nature, specifically the light quark mass and the strength of the electromagnetic interaction. We show that there exist strong correlations between the alpha-particle binding energy and the various energies relevant to the triple-alpha process. We derive limits on the variation of these fundamental parameters from the requirement that sufficient amounts of carbon and oxygen be generated in stars. We also discuss the implications of these results for an anthropic view of the universe.

preprint2012arXiv

Benchmark calculations for elastic fermion-dimer scattering

We present continuum and lattice calculations for elastic scattering between a fermion and a bound dimer in the shallow binding limit. For the continuum calculation we use the Skorniakov-Ter-Martirosian (STM) integral equation to determine the scattering length and effective range parameter to high precision. For the lattice calculation we use the finite-volume method of Lüscher. We take into account topological finite-volume corrections to the dimer binding energy which depend on the momentum of the dimer. After subtracting these effects, we find from the lattice calculation kappa a_fd = 1.174(9) and kappa r_fd = -0.029(13). These results agree well with the continuum values kappa a_fd = 1.17907(1) and kappa r_fd = -0.0383(3) obtained from the STM equation. We discuss applications to cold atomic Fermi gases, deuteron-neutron scattering in the spin-quartet channel, and lattice calculations of scattering for nuclei and hadronic molecules at finite volume.

preprint2012arXiv

Causality bounds for neutron-proton scattering

We consider the constraints of causality and unitarity for the low-energy interactions of protons and neutrons. We derive a general theorem that non-vanishing partial-wave mixing cannot be reproduced with zero-range interactions without violating causality or unitarity. We define and calculate interaction length scales which we call the causal range and the Cauchy-Schwarz range for all spin channels up to J = 3. For some channels we find that these length scales are as large as 5 fm. We investigate the origin of these large lengths and discuss their significance for the choice of momentum cutoff scales in effective field theory and universality in many-body Fermi systems.

preprint2012arXiv

Structure and rotations of the Hoyle state

The excited state of the 12C nucleus known as the "Hoyle state" constitutes one of the most interesting, difficult and timely challenges in nuclear physics, as it plays a key role in the production of carbon via fusion of three alpha particles in red giant stars. In this letter, we present ab initio lattice calculations which unravel the structure of the Hoyle state, along with evidence for a low-lying spin-2 rotational excitation. For the 12C ground state and the first excited spin-2 state, we find a compact triangular configuration of alpha clusters. For the Hoyle state and the second excited spin-2 state, we find a "bent-arm" or obtuse triangular configuration of alpha clusters. We also calculate the electromagnetic transition rates between the low-lying states of 12C.

preprint2012arXiv

Weakly bound molecules trapped with discrete scaling symmetries

When the scattering length is proportional to the distance from the center of the system, two particles are shown to be trapped about the center. Furthermore, their spectrum exhibits discrete scale invariance, whose scale factor is controlled by the slope of the scattering length. While this resembles the Efimov effect, our system has a number of advantages when realized with ultracold atoms. We also elucidate how the emergent discrete scaling symmetry is violated for more than two bosons, which may shed new light on Efimov physics. Our system thus serves as a tunable model system to investigate universal physics involving scale invariance, quantum anomaly, and renormalization group limit cycle, which are important in a broad range of quantum physics.

preprint2011arXiv

Ab initio calculation of the Hoyle state

The Hoyle state plays a crucial role in the hydrogen burning of stars heavier than our sun and in the production of carbon and other elements necessary for life. This excited state of the carbon-12 nucleus was postulated by Hoyle [1] as a necessary ingredient for the fusion of three alpha particles to produce carbon at stellar temperatures. Although the Hoyle state was seen experimentally more than a half century ago [2,3], nuclear theorists have not yet uncovered the nature of this state from first principles. In this letter we report the first ab initio calculation of the low-lying states of carbon-12 using supercomputer lattice simulations and a theoretical framework known as effective field theory. In addition to the ground state and excited spin-2 state, we find a resonance at -85(3) MeV with all of the properties of the Hoyle state and in agreement with the experimentally observed energy. These lattice simulations provide insight into the structure of this unique state and new clues as to the amount of fine-tuning needed in nature for the production of carbon in stars.

preprint2011arXiv

How quantum bound states bounce and the structure it reveals

We investigate how quantum bound states bounce from a hard surface. Our analysis has applications to ab initio calculations of nuclear structure and elastic deformation, energy levels of excitons in semiconductor quantum dots and wells, and cold atomic few-body systems on optical lattices with sharp boundaries. We develop the general theory of elastic reflection for a composite body from a hard wall. On the numerical side we present ab initio calculations for the compression of alpha particles and universal results for two-body states. On the analytical side we derive a universal effective potential that gives the reflection scattering length for shallow two-body states.

preprint2011arXiv

Non-relativistic bound states in a finite volume

We derive general results for the mass shift of bound states with angular momentum l >= 1 in a periodic cubic box in two and three spatial dimensions. Our results have applications to lattice simulations of hadronic molecules, halo nuclei, and Feshbach molecules. The sign of the mass shift can be related to the symmetry properties of the state under consideration. We verify our analytical results with explicit numerical calculations. Moreover, we comment on the relations connecting the effective range parameter, the binding momentum of a given state and the asymptotic normalization coefficient of the corresponding wave function. We give explicit expressions for this relation in the shallow binding limit.

preprint2011arXiv

Precision benchmark calculations for four particles at unitarity

The unitarity limit describes interacting particles where the range of the interaction is zero and the scattering length is infinite. We present precision benchmark calculations for two-component fermions at unitarity using three different ab initio methods: Hamiltonian lattice formalism using iterated eigenvector methods, Euclidean lattice formalism with auxiliary-field projection Monte Carlo, and continuum diffusion Monte Carlo with fixed and released nodes. We have calculated the ground state energy of the unpolarized four-particle system in a periodic cube as a dimensionless fraction of the ground state energy for the non-interacting system. We obtain values 0.211(2) and 0.210(2) using two different Hamiltonian lattice representations, 0.206(9) using Euclidean lattice, and an upper bound of 0.212(2) from fixed-node diffusion Monte Carlo. Released-node calculations starting from the fixed-node result yield a decrease of less than 0.002 over a propagation of 0.4/E_F in Euclidean time, where E_F is the Fermi energy. We find good agreement among all three ab initio methods.

preprint2011arXiv

Topological phases for bound states moving in a finite volume

We show that bound states moving in a finite periodic volume have an energy correction which is topological in origin and universal in character. The topological volume corrections contain information about the number and mass of the constituents of the bound states. These results have broad applications to lattice calculations involving nucleons, nuclei, hadronic molecules, and cold atoms. We illustrate and verify the analytical results with several numerical lattice calculations.

preprint2011arXiv

Volume Dependence of Bound States with Angular Momentum

We derive general results for the mass shift of bound states with angular momentum l >= 1 in a finite periodic volume. Our results have direct applications to lattice simulations of hadronic molecules as well as atomic nuclei. While the binding of S-wave bound states increases at finite volume, we show that the binding of P-wave bound states decreases. The mass shift for D-wave bound states as well as higher partial waves depends on the representation of the cubic rotation group. Nevertheless, the multiplet-averaged mass shift for any angular momentum l can be expressed in a simple form, and the sign of the shift alternates for even and odd l. We verify our analytical results with explicit numerical calculations. We also show numerically that similar volume corrections appear in three-body bound states.

preprint2010arXiv

Causality and the effective range expansion

We derive the generalization of Wigner's causality bounds and Bethe's integral formula for the effective range parameter to arbitrary dimension and arbitrary angular momentum. We also discuss the impact of these constraints on the separation of low- and high-momentum scales and universality in low-energy scattering. Some of our results were summarized earlier in a letter publication. In this work, we present full derivations and several detailed examples.

preprint2010arXiv

Lattice calculations for A=3,4,6,12 nuclei using chiral effective field theory

We present lattice calculations for the ground state energies of tritium, helium-3, helium-4, lithium-6, and carbon-12 nuclei. Our results were previously summarized in a letter publication. This paper provides full details of the calculations. We include isospin-breaking, Coulomb effects, and interactions up to next-to-next-to-leading order in chiral effective field theory.

preprint2009arXiv

Ground state energy of dilute neutron matter at next-to-leading order in lattice chiral effective field theory

We present lattice calculations for the ground state energy of dilute neutron matter at next-to-leading order in chiral effective field theory. This study follows a series of recent papers on low-energy nuclear physics using chiral effective field theory on the lattice. In this work we introduce an improved spin- and isospin-projected leading-order action which allows for a perturbative treatment of corrections at next-to-leading order and smaller estimated errors. Using auxiliary fields and Euclidean-time projection Monte Carlo, we compute the ground state of 8, 12, and 16 neutrons in a periodic cube, covering a density range from 2% to 10% of normal nuclear density.

preprint2008arXiv

Lattice simulations for few- and many-body systems

We review the recent literature on lattice simulations for few- and many-body systems. We focus on methods and results that combine the framework of effective field theory with computational lattice methods. Lattice effective field theory is discussed for cold atoms as well as low-energy nucleons with and without pions. A number of different lattice formulations and computational algorithms are considered, and an effort is made to show common themes in studies of cold atoms and low-energy nuclear physics as well as common themes in work by different collaborations.

preprint2008arXiv

Nuclear effective field theory on the lattice

In the low-energy region far below the chiral symmetry breaking scale (which is of the order of 1 GeV) chiral perturbation theory provides a model-independent approach for quantitative description of nuclear processes. In the two- and more-nucleon sector perturbation theory is applicable only at the level of an effective potential which serves as input in the corresponding dynamical equation. To deal with the resulting many-body problem we put chiral effective field theory (EFT) on the lattice. Here we present the results of our lattice EFT study up to next-to-next-to-leading order in the chiral expansion. Accurate description of two-nucleon phase-shifts and ground state energy ratio of dilute neutron matter up to corrections of higher orders shows that lattice EFT is a promising tool for a quantitative description of low-energy few- and many-body systems.

preprint2008arXiv

The symmetric heavy-light ansatz

The symmetric heavy-light ansatz is a method for finding the ground state of any dilute unpolarized system of attractive two-component fermions. Operationally it can be viewed as a generalization of the Kohn-Sham equations in density functional theory applied to N-body density correlations. While the original Hamiltonian has an exact Z_2 symmetry, the heavy-light ansatz breaks this symmetry by skewing the mass ratio of the two components. In the limit where one component is infinitely heavy, the many-body problem can be solved in terms of single-particle orbitals. The original Z_2 symmetry is recovered by enforcing Z_2 symmetry as a constraint on N-body density correlations for the two components. For the 1D, 2D, and 3D attractive Hubbard models the method is in very good agreement with exact Lanczos calculations for few-body systems at arbitrary coupling. For the 3D attractive Hubbard model there is very good agreement with lattice Monte Carlo results for many-body systems in the limit of infinite scattering length.

preprint1999arXiv

The massless Thirring model in spherical field theory

We use the massless Thirring model to demonstrate a new approach to non-perturbative fermion calculations based on the spherical field formalism. The methods we present are free from the problems of fermion doubling and difficulties associated with integrating out massless fermions. Using a non-perturbative regularization, we compute the two-point correlator and find agreement with the known analytic solution.

preprint1998arXiv

Introduction to spherical field theory

Spherical field theory is a new non-perturbative method for studying quantum field theories. It uses the spherical partial wave expansion to reduce a general d-dimensional Euclidean field theory into a set of coupled one-dimensional systems. The coupled one-dimensional systems are then converted to partial differential equations and solved numerically. We demonstrate the methods of spherical field theory by analyzing Euclidean phi^4 theory in two dimensions.