Source author record

Paulo F. Bedaque

Paulo F. Bedaque 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

39works
17topics
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

39 published item(s)

preprint2021arXiv

The spectrum of qubitized QCD: glueballs in a $S(1080)$ gauge theory

Quantum simulations of QCD require digitization of the infinite-dimensional gluon field. Schemes for doing this with the minimum amount of qubits are desirable. We present a practical digitization for $SU(3)$ gauge theories via its discrete subgroup $S(1080)$. Using a modified action that allows classical simulations down to $a\approx 0.08$ fm, the low-lying glueball spectrum is computed with percent-level precision at multiple lattice spacings and shown to extrapolate to the continuum limit $SU(3)$ results. This suggests that this digitization scheme is sufficient for precision quantum simulations of QCD.

preprint2021arXiv

Universality of a truncated sigma-model

Bosonic quantum field theories, even when regularized using a finite lattice, possess an infinite dimensional Hilbert space and, therefore, cannot be simulated in quantum computers with a finite number of qubits. A truncation of the Hilbert space is then needed and the physical results are obtained after a double limit: one to remove the truncation and another to remove the regulator (the continuum limit). A simpler alternative is to find a model with a finite dimensional Hilbert space belonging to the same universality class as the continuum model (a "qubitization"), so only the space continuum limit is required. A qubitization of the $1+1$ dimensional asymptotically free $O(3)$ nonlinear $σ$-model based on ideas of non-commutative geometry was previously proposed arXiv:1903.06577 and, in this paper, we provide evidence that it reproduces the physics of the $σ$-model both in the infrared and the ultraviolet regimes.

preprint2020arXiv

Complex Paths Around The Sign Problem

The Monte Carlo evaluation of path integrals is one of a few general purpose methods to approach strongly coupled systems. It is used in all branches of Physics, from QCD/nuclear physics to the correlated electron systems. However, many systems of great importance (dense matter inside neutron stars, the repulsive Hubbard model away from half-filling, dynamical and non-equilibrium observables) are not amenable to the Monte Carlo method as it currently stands due to the so-called "sign-problem". We review a new set of ideas recently developed to tackle the sign problem based on the complexification of field space and the Picard-Lefshetz theory accompanying it. The mathematical ideas underpinning this approach, as well as the algorithms so far developed, are described together with non-trivial examples where the method has already been proved successful. Directions of future work, including the burgeoning use of machine learning techniques, are delineated.

preprint2020arXiv

Quantum algorithms for disordered physics

We show how a quantum computer may efficiently simulate a disordered Hamiltonian, by incorporating a pseudo-random number generator directly into the time evolution circuit. This technique is applied to quantum simulation of few-body disordered systems in the large volume limit; in particular, Anderson localization. The method requires a number of (error corrected) qubits proportional to the logarithm of the volume of the system, and each time evolution step requires a number of gates polylogarithmic in the volume. We simulate the method to observe the metal-insulator transition on a three-dimensional lattice. Additionally, we demonstrate the algorithm on a one-dimensional lattice, using physical quantum processors.

preprint2019arXiv

Structure Factors of The Unitary Gas Under Supernova Conditions

We compute with lattice field theory the vector and axial static structure factors of the unitary gas for arbitrary temperature above the superfluid transition and for fugacities 0.1 < z < 1.0. Using the lattice formulation, we calculate beyond the validity of the virial expansion, a commonly used technique in many-body physics. We find qualitative differences in the behavior of the structure factors at high fugacity compared to the predictions of the virial expansion. Due to the large scattering length of neutrons, we expect the unitary gas structure factors to approximate the structure factors of hot neutron gases, and we therefore expect our calculations to be useful in supernova simulations, where neutron gas structure factors are needed to compute in-medium neutrino-neutron scattering rates.

preprint2016arXiv

Fast Estimator of jacobians in Monte Carlo Integration on Lefschetz Thimbles

A solution to the sign problem is the so-called "Lefschetz thimble approach" where the domain of integration for field variables in the path integral is deformed from the real axis to a sub-manifold in the complex space. For properly chosen sub-manifolds ("thimbles") the sign problem disappears or is drastically alleviated. The parametrization of the thimble by real coordinates require the calculation of a jacobian with a computational cost of order O(V^3), where V is proportional to the spacetime volume. In this note we propose two estimators for this jacobian with a computational cost of order O(V). We discuss analytically the regimes where we expect the estimator to work and show numerical examples in two different models.

preprint2016arXiv

Monte Carlo study of real time dynamics

Monte Carlo studies involving real time dynamics are severely restricted by the sign problem that emerges from highly oscillatory phase of the path integral. In this letter, we present a new method to compute real time quantities on the lattice using the Schwinger-Keldysh formalism via Monte Carlo simulations. The key idea is to deform the path integration domain to a complex manifold where the phase oscillations are mild and the sign problem is manageable. We use the previously introduced "contraction algorithm" to create a Markov chain on this alternative manifold. We substantiate our approach by analyzing the quantum mechanical anharmonic oscillator. Our results are in agreement with the exact ones obtained by diagonalization of the Hamiltonian. The method we introduce is generic and in principle applicable to quantum field theory albeit very slow. We discuss some possible improvements that should speed up the algorithm.

preprint2016arXiv

Sign problem and Monte Carlo calculations beyond Lefschetz thimbles

We point out that Monte Carlo simulations of theories with severe sign problems can be profitably performed over manifolds in complex space different from the one with fixed imaginary part of the action. We describe a family of such manifolds that interpolate between the tangent space at one critical point, where the sign problem is milder compared to the real plane but in some cases still severe, and the union of relevant thimbles, where the sign problem is mild but a multimodal distribution function complicates the Monte Carlo sampling. We exemplify this approach using a simple 0 + 1 dimensional fermion model previously used on sign problem studies and show that it can solve the model for some parameter values where a solution using Lefshetz thimbles was elusive.

preprint2015arXiv

An Upper Bound on Neutron Star Masses from Models of Short Gamma-ray Bursts

The discovery of two neutron stars with gravitational masses $\approx 2~M_\odot$ has placed a strong lower limit on the maximum mass of nonrotating neutron stars, and with it a strong constraint on the properties of cold matter beyond nuclear density. Current upper mass limits are much looser. Here we note that, if most short gamma-ray bursts are produced by the coalescence of two neutron stars, and if the merger remnant collapses quickly, then the upper mass limit is constrained tightly. If the rotation of the merger remnant is limited only by mass-shedding (which seems probable based on numerical studies), then the maximum gravitational mass of a nonrotating neutron star is $\approx 2-2.2~M_\odot$ if the masses of neutron stars that coalesce to produce gamma-ray bursts are in the range seen in Galactic double neutron star systems. These limits would be increased by $\sim 4$% in the probably unrealistic case that the remnants rotate at $\sim 30$% below mass-shedding, and by $\sim 15$% in the extreme case that the remnants do not rotate at all. Future coincident detection of short gamma-ray bursts with gravitational waves will strengthen these arguments because they will produce tight bounds on the masses of the components for individual events. If these limits are accurate then a reasonable fraction of double neutron star mergers might not produce gamma-ray bursts. In that case, or in the case that many short bursts are produced instead by the mergers of neutron stars with black holes, the implied rate of gravitational wave detections will be increased.

preprint2015arXiv

Sound velocity bound and neutron stars

It has been conjectured that the velocity of sound in any medium is smaller than the velocity of light in vacuum divided by $\sqrt{3}$. Simple arguments support this bound in non-relativistic and/or weakly coupled theories. The bound has been demonstrated in several classes of strongly coupled theories with gravity duals and is saturated only in conformal theories. We point out that the existence of neutron stars with masses around two solar masses combined with the knowledge of the equation of state of hadronic matter at "low" densities is in strong tension with this bound.

preprint2014arXiv

Hypernuclei and the hyperon problem in neutron stars

The likely presence of $Λ$ baryons in dense hadronic matter tends to soften the equation of state to an extend that the observed heaviest neutron stars are difficult to explain. We analyze this "hyperon problem" with a phenomenological approach. First, we review what can be learned about the interaction of $Λ$ particle with dense matter from the observed hypernuclei and extend this phenomenological analysis to asymmetric matter. We add to this the current knowledge on non-strange dense matter, including its uncertainties, to conclude that the interaction between $Λ$s and dense matter has to become repulsive at densities below three times the nuclear saturation density.

preprint2013arXiv

A new crystalline phase in magnetar crusts

We show that ions at the low densities and high magnetic fields relevant to the outer crust of magnetars form a novel crystalline phase where ions are strongly coupled along the magnetic field and loosely coupled in the transverse direction. The underlying cause is the anisotropic screening of the Coulomb force by electrons in the presence of a strongly quantizing magnetic field which leads to Friedel oscillations in the ion-ion potential. In particular, the Friedel oscillations are much longer-ranged in the direction of the magnetic field than is the case in the absence of magnetic fields, a factor that has been neglected in previous studies. These "Friedel crystals" have very anisotropic elastic moduli, with potentially interesting implications for the Quasi-periodic Oscillations seen in the X-ray flux of magnetars during their giant flares. We find the minimum energy configuration of ions taking into account these anisotropic effects and find that, depending on the density, temperature and magnetic field strength, different lattice structures (fcc, hcp or bcc oriented in different ways in relation to the magnetic field) are favored.

preprint2013arXiv

Friedel crystals and the outer crust of magnetars

The strong magnetic fields found on the surface of magnetars are known to have profound effects in the physics of atoms in magnetar envelopes. We argue that the Friedel oscillations in the Coulomb force between the ions due to electron shielding can, for certain values of the parameters, be the dominant effect determining the crystal structure in the outer crust of magnetars. We estimate the densities and magnetic fields for which this occurs, compute some of the elastic moduli and lattice phonon dispersion relations in this "Friedel crystal" phase.

preprint2013arXiv

Goldstone modes in the neutron star core

We formulate a theory of Goldstone bosons and their interactions in the superfluid and superconducting phase of dense nucleonic matter at densities of relevance to the neutron star core. For typical neutron star temperatures in the range T = 10^6 to 10^9 K, the Goldstone mode associated with rotational symmetry, called angulons, couple weakly to each other and to electrons. Consequently, these modes have anomalously large mean free paths and can contribute to both diffusive and ballistic transport of heat and momentum. In contrast, the two Goldstone bosons associated with density oscillations of the neutron and electron + proton fluids, called superfluid phonons, mix and couple strongly to electrons. They have shorter mean free paths, and their contribution to transport is negligible. Long-wavelength superfluid phonons and angulons can play a role in neutron star seismology, and lead to interesting phenomenology as angulons couple to magnetic fields and have anisotropic dispersion relations.

preprint2013arXiv

The low lying modes of triplet-condensed neutron matter and their effective theory

The condensation of neutrons into a 3P2 superfluid phase occurs at densities relevant for the interior of neutron stars. The triplet pairing breaks rotational symmetry spontaneously and leads to the existence of gapless modes (angulons) that are relevant for many transport coefficients and to the star's cooling properties. We derive the leading terms of the low energy effective field theory, including the leading coupling to electroweak currents, valid for a variety of possible 3P2 phases.

preprint2012arXiv

Neutrino Emission from Helium White Dwarfs with Condensed Cores

The possibility that ions in a helium white dwarf star are in a Bose-Einstein condensed state has been explored recently. In particular, it has been argued that the resulting novel quantum liquid has a new kind of quasiparticle excitation with a phonon-like dispersion relation. We investigate the neutrino emission rate due to this gapless state and the resulting impact on the total luminosity of helium white dwarf stars, as a possible observable way of detecting this exotic phase. If the condensation temperature for the quantum liquid state, which is currently not known very precisely, turns out to be high enough, our calculations indicate that neutrino emission due to the gapless mode would make a large contribution to the total luminosity of the helium white dwarf stars.

preprint2012arXiv

Thermodynamics of nuclear condensates and phase transitions in white dwarfs

We study the thermodynamics of helium at densities relevant for white dwarf physics. We find evidence that, as the temperature is increased, there is first a first order transition between two superconducting phases followed by a second order transition to the normal state. These transitions occur, for realistic densities, at temperatures below the crystallization temperature and the crystalline state is likely to remain as the true ground state of the system. The calculations are performed with a screening but non-dynamical electron background and we comment on the impact of this and other approximations to our result.

preprint2011arXiv

Nuclear condensate and helium white dwarfs

We consider a high density region of the helium phase diagram, where the nuclei form a Bose-Einstein condensate rather than a classical plasma or a crystal. Helium in this phase may be present in helium-core white dwarfs. We show that in this regime there is a new gapless quasiparticle not previously noticed, arising when the constraints imposed by gauge symmetry are taken into account. The contribution of this quasiparticle to the specific heat of a white dwarf core turns out to be comparable in a range of temperatures to the contribution from the particle-hole excitations of the degenerate electrons. The specific heat in the condensed phase is two orders of magnitude smaller than in the uncondensed plasma phase, which is the ground state at higher temperatures, and four orders of magnitude smaller than the specific heat that an ion lattice would provide, if formed. Since the specific heat of the core is an important input for setting the rate of cooling of a white dwarf star, it may turn out that such a change in the thermal properties of the cores of helium white dwarfs has observable implications.

preprint2011arXiv

Quark mass variation constraints from Big Bang nucleosynthesis

We study the impact on the primordial abundances of light elements created by a variation of the quark masses at the time of Big Bang nucleosynthesis (BBN). In order to navigate through the particle and nuclear physics required to connect quark masses to binding energies and reaction rates in a model-independent way, we use lattice QCD data and a hierarchy of effective field theories. We find that the measured Helium-4 abundances put a bound of -1 % <~ d m_q/m_q <~ 0.7 % on a possible variation of quark masses. The effect of quark mass variations on the deuterium abundances can be largely compensated by changes of the baryon-to-photon ratio eta. Including bounds on the variation of eta coming from WMAP results and adding some additional assumptions further narrows the range of allowed values of d m_q/m_q.

preprint2011arXiv

Vortons in dense quark matter

At large baryon number density, it is likely that the ground state of QCD is a color-flavor-locked phase with a K0 condensate. The CFL+K0 phase is known to support superconducting vortex strings, and it has been previously suggested that it may also support vortons, which are superconducting vortex rings. We reexamine the question of the stability of vortons, taking into account electromagnetic effects, which make leading-order contributions to vorton dynamics but were not investigated in previous work. We find that current-carrying and electrically charged vortons can be stabilized either by their angular momentum, by Coulomb repulsion, or by a combination of both effects.

preprint2010arXiv

The phases of deuterium at extreme densities

We consider deuterium compressed to higher than atomic, but lower than nuclear densities. At such densities deuterium is a superconducting quantum liquid. Generically, two superconducting phases compete, a "ferromagnetic" and a "nematic" one. We provide a power counting argument suggesting that the dominant interactions in the deuteron liquid are perturbative (but screened) Coulomb interactions. At very high densities the ground state is determined by very small nuclear interaction effects that probably favor the ferromagnetic phase. At lower densities the symmetry of the theory is effectively enhanced to SU(3), and the quantum liquid enters a novel phase, neither ferromagnetic nor nematic. Our results can serve as a starting point for investigations of the phase dynamics of deuteron liquids, as well as exploration of the stability and dynamics of the rich variety of topological objects that may occur in phases of the deuteron quantum liquid, which range from Alice strings to spin skyrmions to Z_2 vortices.

preprint2009arXiv

Can fermions save large N dimensional reduction?

This paper explores whether Eguchi-Kawai reduction for gauge theories with adjoint fermions is valid. The Eguchi-Kawai reduction relates gauge theories in different numbers of dimensions in the large $N$ limit provided that certain conditions are met. In principle, this relation opens up the possibility of learning about the dynamics of 4D gauge theories through techniques only available in lower dimensions. Dimensional reduction can be understood as a special case of large $N$ equivalence between theories related by an orbifold projection. In this work, we focus on the simplest case of dimensional reduction, relating a 4D gauge theory to a 3D gauge theory via an orbifold projection. A necessary condition for the large N equivalence between the 4D and 3D theories to hold is that certain discrete symmetries in the two theories must not be broken spontaneously. In pure 4D Yang-Mills theory, these symmetries break spontaneously as the size of one of the spacetime dimensions shrinks. An analysis of the effect of adjoint fermions on the relevant symmetries of the 4D theory shows that the fermions help stabilize the symmetries. We consider the same problem from the point of view of the lower dimensional 3D theory and find that, surprisingly, adjoint fermions are not generally enough to stabilize the necessary symmetries of the 3D theory. In fact, a rich phase diagram arises, with a complicated pattern of symmetry breaking. We discuss the possible causes and consequences of this finding.

preprint2009arXiv

Meson-Baryon Scattering Parameters from Lattice QCD with an Isospin Chemical Potential

The extraction of meson-baryon scattering parameters from lattice QCD is complicated by the necessity, in some channels, of including annihilation diagrams. We consider a strategy to avoid the need for these extremely costly and noisy contributions. The strategy is based on simulations with an isospin chemical potential which, contrary to a baryon chemical potential, has no sign problem. When the isospin chemical potential is larger than a critical value, a charged pion condensate forms. Baryons propagating in the pion condensate will have their masses modified by pion-baryon scattering parameters. Consequently these parameters can be extracted from lattice QCD simulations in an isospin chemical potential, and we detail precisely which low-energy constants using baryon chiral perturbation theory.

preprint2006arXiv

Finite volume corrections to pi-pi scattering

Lattice QCD studies of hadron-hadron interactions are performed by computing the energy levels of the system in a finite box. The shifts in energy levels proportional to inverse powers of the volume are related to scattering parameters in a model independent way. In addition, there are non-universal exponentially suppressed corrections that distort this relation. These terms are proportional to exp(-m_pi L) and become relevant as the chiral limit is approached. In this paper we report on a one-loop chiral perturbation theory calculation of the leading exponential corrections in the case of I=2 pi-pi scattering near threshold.

preprint2006arXiv

Superfluid phases of the three-species fermion gas

We discuss the zero temperature phase diagram of a dilute gas with three fermionic species. We make use of solvable limits to conjecture the behavior of the system in the "unitary" regions. The physics of the Thomas-Efimov effect plays a role in these considerations. We find a rich phase diagram with superfluid, gapless superfluid and inhomogeneous phases with different symmetry breaking patterns. We then discuss one particular possible experimental implementation in a system of ^6Li atoms and the possible phases arising in this system as an external magnetic field is varied across three overlaping Feshbach resonances. We also suggest how to experimentally distinguish the different phases.

preprint2005arXiv

Twisted valence quarks and hadron interactions on the lattice

We consider QCD with valence and sea quarks obeying different boundary conditions. We point out that the energy of low lying two hadron states do not depend on the boundary condition of the sea quarks (up to exponentially small corrections). Thus, the advantages in using twisted boundary conditions on the lattice QCD extraction of nucleon-nucleon phase shifts can be gained without the need of new gauge configurations, even in fully unquenched calculations.

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.

preprint2000arXiv

From Hadrons to Nuclei: Crossing the Border

The study of nuclei predates by many years the theory of quantum chromodynamics. More recently, effective field theories have been used in nuclear physics to ``cross the border'' from QCD to a nuclear theory. We are now entering the second decade of efforts to develop a perturbative theory of nuclear interactions using effective field theory. This work describes the current status of these efforts.

preprint1999arXiv

Color Superconductivity in Asymmetric Matter

The influence of different chemical potential for different flavors on color superconductivity is analyzed. It is found that there is a first order transition as the asymmetry grows. This transition proceeds through the formation of bubbles of low density, flavor asymmetric normal phase inside a high density, superconducting phase with a gap {\it larger} than the one found in the symmetric case. For small fixed asymmetries the system is normal at low densities and superconducting only above some critical density. For larger asymmetries the two massless quarks system stays in the mixed state for arbitrarily high densities.

preprint1995arXiv

Baryon masses at second order in large-$N$ chiral perturbation theory

We consider flavor breaking in the the octet and decuplet baryon masses at second order in large-$N$ chiral perturbation theory, where $N$ is the number of QCD colors. We assume that $1/N \sim 1/N_F \sim m_s / Λ\gg m_{u,d}/Λ, α_{EM}$, where $N_F$ is the number of light quark flavors, and $m_{u,d,s} / Λ$ are the parameters controlling $SU(N_F)$ flavor breaking in chiral perturbation theory. We consistently include non-analytic contributions to the baryon masses at orders $m_q^{3/2}$, $m_q^2 \ln m_q$, and $(m_q \ln m_q) / N$. The $m_q^{3/2}$ corrections are small for the relations that follow from $SU(N_F)$ symmetry alone, but the corrections to the large-$N$ relations are large and have the wrong sign. Chiral power-counting and large-$N$ consistency allow a 2-loop contribution at order $m_q^2 \ln m_q$, and a non-trivial explicit calculation is required to show that this contribution vanishes. At second order in the expansion, there are eight relations that are non-trivial consequences of the $1/N$ expansion, all of which are well satisfied within the experimental errors. The average deviation at this order is $7 \MeV$ for the $\De I = 0$ mass differences and $0.35 \MeV$ for the $\De I \ne 0$ mass differences, consistent with the expectation that the error is of order $1/N^2 \sim 10\%$.

preprint1995arXiv

Chiral perturbation theory analysis of baryon temperature mass shifts

We compute the finite temperature pole mass shifts of the octet and decuplet baryons using heavy baryon chiral perturbation theory and the 1/N_c expansion, where N_c is the number of QCD colors. We consider temperatures of the order of the pion mass m_π, and truncate the chiral and 1/N_c expansions assuming that m_π\sim 1/N_c. There are three scales in the problem: the temperature T, the pion mass m_π, and the octet--decuplet mass difference. Therefore, the result is not simply a power series in T. We find that the nucleon and Δtemperature mass shifts are opposite in sign, and that their mass difference changes by 20% in the temperature range 90 MeV < T < 130 MeV, that is the range where the freeze out in relativistic heavy ion collisions is expected to occur. We argue that our results are insensitive to the neglect of 1/N_c- supressed effects; the main purpose of the 1/N_c expansion in this work is to justify our treatment of the decuplet states.