Pairing gaps in Hartree-Fock Bogoliubov theory with the Gogny D1S interaction

As part of a program to study odd-A nuclei in the Hartree-Fock-Bogoliubov (HFB) theory, we have developed a new calculational tool to find the HFB minima of odd-A nuclei based on the gradient method and using interactions of Gogny's form. The HFB minimization includes both time-even and time-odd fields in the energy functional, avoiding the commonly used "filling approximation". Here we apply the method to calculate neutron pairing gaps in some representative isotope chains of spherical and deformed nuclei, namely the Z=8,50 and 82 spherical chains and the Z=62 and 92 deformed chains. We find that the gradient method is quite robust, permitting us to carry out systematic surveys involving many nuclei. We find that the time-odd field does not have large effect on the pairing gaps calculated with the Gogny D1S interaction. Typically, adding the T-odd field as a perturbation increases the pairing gap by ~100 keV, but the re-minimization brings the gap back down. This outcome is very similar to results reported for the Skyrme family of nuclear energy density functionals. Comparing the calculated gaps with the experimental ones, we find that the theoretical errors have both signs implying that the D1S interaction has a reasonable overall strength. However, we find some systematic deficiencies comparing spherical and deformed chains and comparing the lighter chains with the heavier ones. The gaps for heavy spherical nuclei are too high, while those for deformed nuclei tend to be too low. The calculated gaps of spherical nuclei show hardly any A-dependence, contrary to the data. Inclusion of the T-odd component of the interaction does not change these qualitative findings

Finding high-order analytic post-Newtonian parameters from a high-precision numerical self-force calculation

We present a novel analytic extraction of high-order post-Newtonian (pN) parameters that govern quasi-circular binary systems. Coefficients in the pN expansion of the energy of a binary system can be found from corresponding coefficients in an extreme-mass-ratio inspiral (EMRI) computation of the change $\Delta U$ in the redshift factor of a circular orbit at fixed angular velocity. Remarkably, by computing this essentially gauge-invariant quantity to accuracy greater than one part in $10^{225}$, and by assuming that a subset of pN coefficients are rational numbers or products of $\pi$ and a rational, we obtain the exact analytic coefficients. We find the previously unexpected result that the post-Newtonian expansion of $\Delta U$ (and of the change $\Delta\Omega$ in the angular velocity at fixed redshift factor) have conservative terms at half-integral pN order beginning with a 5.5 pN term. This implies the existence of a corresponding 5.5 pN term in the expansion of the energy of a binary system. Coefficients in the pN series that do not belong to the subset just described are obtained to accuracy better than 1 part in $10^{265-23n}$ at $n$th pN order. We work in a radiation gauge, finding the radiative part of the metric perturbation from the gauge-invariant Weyl scalar $\psi_0$ via a Hertz potential. We use mode-sum renormalization, and find high-order renormalization coefficients by matching a series in $L=\ell+1/2$ to the large-$L$ behavior of the expression for $\Delta U$. The non-radiative parts of the perturbed metric associated with changes in mass and angular momentum are calculated in the Schwarzschild gauge

The box diagram in Yukawa theory

We present a light-front calculation of the box diagram in Yukawa theory. The covariant box diagram is finite for the case of spin-1/2 constituents exchanging spin-0 particles. In light-front dynamics, however, individual time-ordered diagrams are divergent. We analyze the corresponding light-front singularities and show the equivalence between the light-front and covariant results by taming the singularities.Comment: 21 pages, 17 figures. submittes to Phys. Rev.

Topics in Chiral Perturbation Theory

I consider some selected topics in chiral perturbation theory (CHPT). For the meson sector, emphasis is put on processes involving pions in the isospin zero S-wave which require multi-loop calculations. The advantages and shortcomings of heavy baryon CHPT are discussed. Some recent results on the structure of the baryons are also presented.Comment: 30 pp, TeX, Review talk, Third Workshop on High Energy Particle Physics (WHEPP III), Madras, India, January 1994. 7 figures available upon request. CRN--94/0

Nonlocal Scalar Quantum Field Theory: Functional Integration, Basis Functions Representation and Strong Coupling Expansion

Nonlocal QFT of one-component scalar field $\varphi$ in $D$-dimensional Euclidean spacetime is considered. The generating functional (GF) of complete Green functions $\mathcal{Z}$ as a functional of external source $j$, coupling constant $g$, and spatial measure $d\mu$ is studied. An expression for GF $\mathcal{Z}$ in terms of the abstract integral over the primary field $\varphi$ is given. An expression for GF $\mathcal{Z}$ in terms of integrals over the primary field and separable Hilbert space (HS) is obtained by means of a separable expansion of the free theory inverse propagator $\hat{L}$ over the separable HS basis. The classification of functional integration measures $\mathcal{D}\left[\varphi\right]$ is formulated, according to which trivial and two nontrivial versions of GF $\mathcal{Z}$ are obtained. Nontrivial versions of GF $\mathcal{Z}$ are expressed in terms of $1$-norm and $0$-norm, respectively. The definition of the $0$-norm generator $\varPsi$ is suggested. Simple cases of sharp and smooth generators are considered. Expressions for GF $\mathcal{Z}$ in terms of integrals over the separable HS with new integrands are obtained. For polynomial theories $\varphi^{2n},\, n=2,3,4,\ldots,$ and for the nonpolynomial theory $\sinh^{4}\varphi$, integrals over the separable HS in terms of a power series over the inverse coupling constant $1/\sqrt{g}$ for both norms ($1$-norm and $0$-norm) are calculated. Critical values of model parameters when a phase transition occurs are found numerically. A generalization of the theory to the case of the uncountable integral over HS is formulated. A comparison of two GFs $\mathcal{Z}$, one in the case of uncountable HS integral and one obtained using the Parseval-Plancherel identity, is given.Comment: 26 pages, 2 figures; v2: significant additions in the text; prepared for the special issue "QCD and Hadron Structure" of the journal Particles; v3: minimal corrections; v4: paragraphs added related to Reviewer comment

Towards an understanding of isospin violation in pion-nucleon scattering

We investigate isospin breaking in low-energy pion-nucleon scattering in the framework of chiral perturbation theory. This work extends the systematic analysis of [1] to the energy range above threshold. Various relations, which identically vanish in the limit of isospin symmetry, are used to quantify isospin breaking effects. We study the energy dependence of the S- and P-wave projections of these ratios and find dramatic effects in the S-waves of those two relations which are given in terms of isoscalar quantities only. This effect drops rather quickly with growing center-of-mass energy.Comment: 12 pp, REVTeX, 8 figs, FZJ-IKP(TH)-2000-2

Improved Moving Puncture Gauge Conditions for Compact Binary Evolutions

Robust gauge conditions are critically important to the stability and accuracy of numerical relativity (NR) simulations involving compact objects. Most of the NR community use the highly robust---though decade-old---moving-puncture (MP) gauge conditions for such simulations. It has been argued that in binary black hole (BBH) evolutions adopting this gauge, noise generated near adaptive-mesh-refinement (AMR) boundaries does not converge away cleanly with increasing resolution, severely limiting gravitational waveform accuracy at computationally feasible resolutions. We link this noise to a sharp (short-wavelength), initial outgoing gauge wave crossing into progressively lower resolution AMR grids, and present improvements to the standard MP gauge conditions that focus on stretching, smoothing, and more rapidly settling this outgoing wave. Our best gauge choice greatly reduces gravitational waveform noise during inspiral, yielding less fluctuation in convergence order and $\sim 40$ lower waveform phase and amplitude errors at typical resolutions. Noise in other physical quantities of interest is also reduced, and constraint violations drop by more than an order of magnitude. We expect these improvements will carry over to simulations of all types of compact binary systems, as well as other $N$+1 formulations of gravity for which MP-like gauge conditions can be chosen.Comment: 25 pages, 16 figures, 2 tables. Matches published versio

Applications of Partially Quenched Chiral Perturbation Theory

Partially quenched theories are theories in which the valence- and sea-quark masses are different. In this paper we calculate the nonanalytic one-loop corrections of some physical quantities: the chiral condensate, weak decay constants, Goldstone boson masses, B_K and the K+ to pi+ pi0 decay amplitude, using partially quenched chiral perturbation theory. Our results for weak decay constants and masses agree with, and generalize, results of previous work by Sharpe. We compare B_K and the K+ decay amplitude with their real-world values in some examples. For the latter quantity, two other systematic effects that plague lattice computations, namely, finite-volume effects and unphysical values of the quark masses and pion external momenta are also considered. We find that typical one-loop corrections can be substantial.Comment: 22 pages, TeX, refs. added, minor other changes, version to appear in Phys. Rev.