Researcher profile

Gautam Rupak

Gautam Rupak contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
25works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

25 published item(s)

preprint2019arXiv

Bayesian analysis of capture reactions $^3\mathrm{He}(α,γ)^7\mathrm{Be}$ and $^3\mathrm{H}(α,γ)^7\mathrm{Li}$

Bayesian analysis of the radiative capture reactions $^3\mathrm{He}(α,γ)^7\mathrm{Be}$ and $^3\mathrm{H}(α,γ)^7\mathrm{Li}$ are performed to draw inferences about the cross sections at threshold. We do a model comparison of two competing effective field theory power countings for the capture reactions. The two power countings differ in the contribution of two-body electromagnetic currents. In one power counting, two-body currents contribute at leading order, and in the other they contribute at higher orders. The former is favored for $^3\mathrm{He}(α,γ)^7\mathrm{Be}$ if elastic scattering data in the incoming channel is considered in the analysis. Without constraints from elastic scattering data, both the power countings are equally favored. For $^3\mathrm{H}(α,γ)^7\mathrm{Li}$, the first power counting with two-body current contributions at leading order is favored with or without constraints from elastic scattering data.

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.

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

Electromagnetic form factors of one neutron halos with spin 1/2+ ground state

The electromagnetic form factors for single neutron halo nuclei Be-11, C-15 and C-19 are calculated. The calculations are performed in halo effective field theory (EFT) where the halo nuclei are approximated as made of a single neutron and a core. The form factors depend on the single neutron separation energy, the s-wave neutron-core scattering effective range and a two-body current. The EFT expressions are presented to leading order for C-15 and next-to-leading order for Be-11 and C-19.

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.

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

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

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.

preprint2012arXiv

Radiative Neutron Capture on Carbon-14 in Effective Field Theory

The cross section for radiative capture of neutron on carbon-14 is calculated using the model-independent formalism of halo effective field theory. The dominant contribution from E1 transition is considered, and the cross section is expressed in terms of elastic scattering parameters of the effective range expansion. Contributions from both resonant and non-resonant interaction are calculated. Significant interference between these leads to a capture contribution that deviates from simple Breit-Wigner resonance form.

preprint2012arXiv

The r-mode instability in strange stars with a crystalline crust

The r-mode instability, believed to limit the rotation speed of compact stars, can provide empirical confirmation for the existence of stable deconfined phases of quark matter that are predicted by weak coupling calculations in Quantum Chromodynamics. We construct a model for strange quark stars as heavy as 2 solar masses that are made of superconducting quark matter in the bulk and a thin crystalline quark matter crust. This crystalline quark crust is sufficiently robust to withstand r-mode heating and viscous rubbing for realistic mode amplitudes O(10^{-2}), unlike a crust made of neutron-rich nuclei. The dissipation provided by viscous rubbing at the core-crust boundary is both necessary and sufficient to obtain stable rotation speeds that are consistent with the majority of rapidly spinning pulsars in low mass X-ray binaries. Our analysis implies that while bare strange stars are ruled out by the existence of rapidly spinning pulsars, a strange star with a quark matter crust is a distinct possibility.

preprint2011arXiv

Resonance Contribution to Radiative Neutron Capture on Lithium-7

Using halo effective field theory, we provide a model-independent calculation of the radiative neutron capture on lithium-7 over an energy range where the contribution from the 3+ resonance becomes important. One finds that a satisfactory description of the capture reaction, in the present single-particle approximation, suggests the use of a resonance width about three times larger than the experimental value. We also present power counting arguments that establish a hierarchy for electromagnetic one- and two-body currents.

preprint2010arXiv

Constraining phases of quark matter with studies of r-mode damping in neutron stars

The r-mode instability in rotating compact stars is used to constrain the phase of matter at high density. The color-flavor-locked phase with kaon condensation (CFL-K0) and without (CFL) is considered in the temperature range 10^8K < T <10^{11} K. While the bulk viscosity in either phase is only effective at damping the r-mode at temperatures T > 10^{11} K, the shear viscosity in the CFL-K0 phase is the only effective damping agent all the way down to temperatures T > 10^8 K characteristic of cooling neutron stars. However, it cannot keep the star from becoming unstable to gravitational wave emission for rotation frequencies f ~ 56-11 Hz at T ~ 10^8-10^9 K. Stars composed almost entirely of CFL or CFL-K0 matter are ruled out by observation of rapidly rotating neutron stars, indicating that dissipation at the quark-hadron interface or nuclear crust interface must play a key role in damping the instability.

preprint2010arXiv

Radiative Neutron Capture on Lithium-7

The radiative neutron capture on lithium-7 is calculated model independently using a low energy halo effective field theory. The cross section is expressed in terms of scattering parameters directly related to the S-matrix element. The cross section depends on the poorly known p-wave effective range parameter r. This constitutes the leading order uncertainty in traditional model calculations. It is explicitly demonstrated by comparing with potential model calculations. A single parameter fit describes the low energy data extremely well and yields r = -1.47 fm^{-1}.

preprint2008arXiv

Viscous damping of r-mode oscillations in compact stars with quark matter

We determine characteristic timescales for the viscous damping of r-mode oscillations in rapidly rotating compact stars that contain quark matter. We present results for the color-flavor-locked (CFL) phase of dense quark matter, in which the up, down and strange quarks are gapped, as well as the normal (ungapped) quark phase. While the ungapped quark phase supports a temperature window between 10^8 K and 5x10^9 K where the r-mode is damped even for rapid rotation, the r-mode in a rapidly rotating pure CFL star is not damped in the temperature range 10^10 K - 10^11 K. Rotating hybrid stars with quark matter cores display an instability window whose width is determined by the amount of quark matter present, and they can have large spin frequencies outside this window. Except at high temperatures T > 10^10 K, the presence of a quark phase allows for larger critical frequencies and smaller spin-periods compared to rotating neutron stars. If low-mass X-ray binaries contain a large amount of ungapped or CFL quark matter, then our estimates of the r-mode instability suggest that there should be a population of rapidly rotating binaries at frequencies greater than 1000 Hz which have not yet been observed.

preprint2003arXiv

Chiral perturbation theory at O(a^2) for lattice QCD

We construct the chiral effective Lagrangian for two lattice theories: one with Wilson fermions and the other with Wilson sea fermions and Ginsparg-Wilson valence fermions. For each of these theories we construct the Symanzik action through order $a^2$. The chiral Lagrangian is then derived, including terms of order $a^2$, which have not been calculated before. We find that there are only few new terms at this order. Corrections to existing coefficients in the continuum chiral Lagrangian are proportional to $a^2$, and appear in the Lagrangian at order $a^2 p^2$ or higher. Similarly, O(4) symmetry breaking terms enter the Symanzik action at order $a^2$, but contribute to the chiral Lagrangian at order $a^2 p^4$ or higher. We calculate the light meson masses in chiral perturbation theory for both lattice theories. At next-to-leading order, we find that there are no order $a^2$ corrections to the valence-valence meson mass in the mixed theory due to the enhanced chiral symmetry of the valence sector.

preprint2003arXiv

Goldstone Bosons in the 3P2 Superfluid Phase of Neutron Matter and Neutrino Emission

At the high densities present in the interior of neutron stars, the neutrons are condensed into the 3P2 superfluid phase. While this condensation has little impact on the equation of state, it can have an important role in determining the low-temperature energy-momentum transport properties. The spontaneous breaking of baryon number by the condensate gives rise to the familiar Goldstone boson, but in addition, the spontaneous breaking of rotational invariance by the condensate gives rise to three Goldstone bosons, in general, one for each broken generator of rotations. These Goldstone bosons, which couple to the Z0, provide a new mechanism for neutrino emission. Using a low-energy effective field theory to describe the dynamics of these Goldstone bosons we estimate the neutrino emissivity of dense neutron matter and show that their annihilation is the dominant energy-loss mechanism over a range of temperatures.

preprint2002arXiv

Dilute resonating gases and the third virial coefficient

We study dilute gases with short range interactions and large two-body scattering lengths. At temperatures between the condensation temperature and the scale set by therange of the potential there is a high degree of universality. The first two terms in the expansion of thermodynamic functions in powers of the fugacity $z$, which measures the diluteness of the system, are determined by the scattering length only. The term proportional to $z^3$ depends only on one new parameter describing the three-body physics. We compute the third term of the expansion and show that, for many values of this new parameter, the $z^3$ term may be the dominant one.

preprint1999arXiv

Improving the Convergence of NN Effective Field Theory

We study a low-energy effective field theory (EFT) describing the NN system in which all exchanged particles are integrated out. We show that fitting the residue of the 3S1 amplitude at the deuteron pole, rather than the 3S1 effective range, dramatically improves the convergence of deuteron observables in this theory. Reproducing the residue ensures that the tail of the deuteron wave function, which is directly related to NN scattering data via analytic continuation, is correctly reproduced in the EFT at next-to-leading order. The role of multi-nucleon-electroweak operators which produce deviations from effective-range theory can then be explicitly separated from the physics of the wave function tail. Such an operator contributes to the deuteron quadrupole moment, mu_Q, at low order, indicating a sensitivity to short-distance physics. This is consistent with the failure of impulse approximation calculations in NN potential models to reproduce mu_Q. The convergence of NN phase shifts in the EFT is unimpaired by the use of this new expansion.

preprint1997arXiv

Non-1/m_b^n Power Suppressed Contributions to Inclusive b->s l+ l- Decays

We compute non-perturbative contributions to b->s l+ l- that are not explicitly suppressed by powers of the b-quark mass. They are proportional to lambda_2 and arise from an interference between the free-quark amplitude and higher order terms in the matrix element of a four-quark operator. This correction is found to be small over most of the dalitz plot except near the charm threshold. Unfortunately, the perturbative computation we have performed is invalid near charm threshold and we do not except to see the structure found at lowest order reproduced in the data. We conclude that these non-perturbative contributions do not significantly modify the previous analysis of b->s l+ l-.