The Crusts of Accreting Neutron Stars 0.01 Edward Brown, Michigan State University 1

Total Page:16

File Type:pdf, Size:1020Kb

The Crusts of Accreting Neutron Stars 0.01 Edward Brown, Michigan State University 1 depth (m) envelope The Crusts of Accreting Neutron Stars 0.01 Edward Brown, Michigan State University 1 0.1 Neutron stars should have a km-thick crust composed of nuclei, electrons, and 1 free neutrons. unstable H/He 5 10 burning/rp-process Accretion pushes matter 10 through this crust and unstable 12C 8 induces nuclear reactions. 10 burning electron Observing the response of 100 captures 1011 the star to these reactions neutron drip allows us to infer the nuclear pasta properties of matter in the 1000 deep crust and core. ≈1014 -3 core A. Piro, Carnegie Obs. ρ (g cm ) SYMMETRY ENERGY AT SUBNUCLEAR DENSITIES AND ... PHYSICAL REVIEW C 75,015801(2007) 2 2 the nucleon part nw mnc nn mpc np [see Eq. (1)] and 0.20 + + the electron part (15). By comparing the resultant six energy inner crust densities, we can determine the equilibrium phase. 0.15 average proton fraction G H 0.10 I D III. EQUILIBRIUM SIZE AND SHAPE OF NUCLEI E proton fraction B A We proceed to show the results for the equilibrium nuclear 0.05 matter configuration obtained for various sets of the EOS F C 0.00 parameters L and K0 as shown in Fig. 1.Theseparametersare 5 6 2 3 4 5 6 2 3 4 5 6 still uncertain because they are little constrained from the mass 0.001 0.01 0.1 n (fm-3) and radius data for stable nuclei [2]. As we shall see, the charge b number of spherical nuclei and the density region containing FIG. 4. (Color online) The average proton fraction as a function bubbles and nonspherical nuclei have a strong correlation of nb,calculatedfortheEOSmodelsA–I. with L. We first focus on spherical nuclei, which constitute an obtained in earlier investigations [1,7,10], will be discussed equilibrium state in the low-density region. We calculate the below in terms of fission instability. charge number of the equilibrium nuclide as a function of n for b The average proton fraction, which is the charge number the EOS models A–I as depicted in Fig. 2. Note that the recent divided by the total nucleon number in the cell, is plotted in GFMC calculations of the energy of neutron matter based on Fig. 4.Weobservethatthedependenceoftheaverageproton the Argonne v8’ potential [15]areclosetothebehaviorof fraction on the EOS models is similar to that of Z.Wealso the model E. Hereafter we will thus call the model E as a find that the average proton fraction basically decreases with typical one. The result is shown in Fig. 3.Fordensitiesbelow baryon density. This is a feature coming from the fact that as 0.01 fm 3,thecalculateddensitydependenceofthecharge − the baryon density increases, the electron chemical potential number∼ Z is almost flat, a feature consistent with the results increases under charge neutrality and then the nuclei become in earlier investigations [1]. More important, the calculated more neutron-rich under weak equilibrium. charge number is larger for the EOS models having smaller L, We next consider the density region where bubbles and and this difference in Z is more remarkable at higher densities. nonspherical nuclei appear in equilibrium, i.e., the density As we shall see later in this section, this property of region of the “pasta” phases. We start with such a density Z is related to the tendency that with increasing L,the region calculated for the EOS models A–I. The results are nuclear density decreases while the density of the neutron plotted in Fig. 5.ExceptforthemodelC,weobtainthe gas increases. Note that Z is, within a liquid-drop model [1], successive first order transitions with increasing density: determined by the size equilibrium condition relating the sphere cylinder slab cylindrical hole spherical Coulomb and surface energies in such a way that Z increases hole →uniform matter.→ A marked→ correlation of the→ upper end with increasing surface tension. Because the Thomas-FermiPasta and crust-core transition of the→ density region with the parameter L can be observed by model adopted here can be mapped onto a compressible liquid- drop model [2], the present results may well be interpreted in Oyamatsu & Iida ’07 terms of the liquid-drop model. In fact we shall estimate the 0.14 G surface tension from the Thomas-Fermi model as a function of H L and discuss how the surface tension depends on the nuclear I ) pasta nuclei density and the neutron sea density. -3 0.12 fission instability We also note that the density at which the phase with proton clustering spherical nuclei ceases to be in the ground state is between 0.10 3 D 0.05 and 0.07 fm− . This result, consistent with the results E A 0.08 F 100 B C G 80 inner crust H 0.06 D nucleon number density (fm proton number Z I 60 E 0.04 40 0 20 40 60 80 100 120 140 160 A proton number 20 B L (MeV) C F FIG. 5. (Color online) The density region containing bubbles and 05 6 2 3 4 5 6 2 3 4 5 6 0.001 0.01 0.1 nonspherical nuclei as a function of L,calculatedfortheEOSmodels -3 nb (fm ) A–I. For comparison, the density corresponding to u 1/8 in the phase with spherical nuclei and the onset density, n(Q=), of proton FIG. 3. (Color online) The charge number of spherical nuclei as clustering in uniform nuclear matter, which will be discussed in afunctionofnb,calculatedfortheEOSmodelsA–I. Sec. IV, are also plotted by circles and crosses, respectively. 015801-5 schematic of Z, A in (outer) crust / (/ ;) ; #(;, /)=B7" B4" B" − B$ start with a simple liquid drop − − " − "/ µ = µ µ enforce beta-equilibrium, F O − Q / charge neutrality: Ye decreases µF B$" − :F + linearly with rising electron ≈ − B B ! " " ! " " chemical potential compute neutron separation ∂# / 4O = = B7 B" : B$: " energy Sn and set it zero to fnd ∂/ − − F − F ! "; density of neutron drip: # $ 4 = µ B .F7 ndrip ≈ 0.001 n0 O → ⇒ F ≈ 7 ≈ ionic lattice. Horowitz et al. (2007) nuclear pasta. Horowitz et al. (2004) µ kT, m c2 e e Z2e2 kT a ionic lattice. Horowitz et al. (2007) nuclear pasta. Horowitz et al. (2004) http://news.nationalgeographic.com µ kT, m c2 e e Z2e2 kT a structure of crust Rcore = 11 km, Mcore = 1.6 Msun, accreted composition outer inner 924 BROWN & BILDSTEN Vol. 496 5. OHMIC DIFFUSION IN THE DEEP OCEAN AND CRUST reactions occurring deep in the crust. As in that section, we assume that temperature is a function of only depth y, as at Present uncertainties in the composition of matter after these depths the accreted matter will have spread around hydrogen/helium burning prohibit a calculation of the sub- the star. We plot in Figure 9 (Q 1.0) the local Ohmic sequent chemical evolution of the ocean for hydrogen-rich \ di†usion time (solid lines) for accretion rates of 0.5, 1.0, and accretion. However, even though the composition is not 5.0 times Eddington. We also show the Ñow time over a well known, we can still use the thermal pro⇡les as estimates scale height, t y/m5 (dashed lines). A few conclusions are of the crust temperatures. This is important to the evolution 4 immediate. First,fl where the ions vibrate classically (T Z ), of the magnetic ⇡eld, as the accretion-induced heating of the # the ratio /t is nearly independent of depth until near crust reduces its conductivity and hastens the Ohmic di†u- q neutron dripdiff ; flmoreover, for accretion rates m5 /m5 [ sion of crust magnetic ⇡elds (Geppert & Urpin 1994). This 0.23(A/2Z)2(T /5 ] 108 K), the di†usion time is alwaysEdd heating also increases the mass of the ocean. These e†ects greater than the time for matter to Ñow through one scale have been considered (Romani 1990; Geppert & Urpin height. Second, impurity scattering is unimportant through- 1994; Pethick & Sahrling 1995; Urpin & Geppert 1995; out the crust for Q [ 1. Because we placed the heat sources Urpin & Geppert 1996; Konar & Bhattacharya 1997) for 904 G. Ushomirsky, C. Cutler and L. Bildsten at a ⇡xed depth y \ y , the thermal gradient changes sign stars accreting at M 10 M yr . w 0 [ ~9 _ ~1 there (see ° 4.3). Electron captures remove pressure support yields There are, however, a few neutron stars accreting globally crust | observables and therefore decrease the pressure scale height, causing the at or near the Eddington rate. There are two X-ray pulsars 3 2 2 abrupt decrease in the Ohmic timescale (solid line, Fig. 9). c Vd ha Ngw ; 8 (LMC X-4 and SMC X-1) and the six bright ““ Z ÏÏ sources Once neutron pressure dominates the equation of state, the G (Sco X-1, GX 5[1, GX 349]2, GX 17]2, GX 340]0, Cyg scale height again increases with depth. In this region, the which, when we write it in terms of Q quiescent(using equation 3)thermal gives emission from transients X-2). The accreted material will have spread over the sur- Ñow timescale is always longer than the di†usion timescale 22 faces of these star for column densities Z1014 g cm~2 (see The Radius of Neutron Stars 9 for . 1=2 2 (from Guillot et al. ’13) m5 [ m5 16 p GQ22V ° 6), so that a spherically symmetrical approach is war- Edd ha : 9 ranted for this calculation. We thus use our solutions for the magnetic field evolution 5 3 dc4 thermal pro⇡le of the deep crust at accretion rates m5 D m5 (from Brown & Bildsten 98) When the angular momentum1 loss by gravitational radiation Edd − 0.01 (see ° 4.3) to estimate the Ohmic di†usion timescales in the balances angular momentum gain by accretion, deep crust of these neutron stars.
Recommended publications
  • Astromaterial Science and Nuclear Pasta
    Astromaterial Science and Nuclear Pasta M. E. Caplan∗ and C. J. Horowitzy Center for Exploration of Energy and Matter and Department of Physics, Indiana University, Bloomington, IN 47405, USA (Dated: June 28, 2017) We define `astromaterial science' as the study of materials in astronomical objects that are qualitatively denser than materials on earth. Astromaterials can have unique prop- erties related to their large density, though they may be organized in ways similar to more conventional materials. By analogy to terrestrial materials, we divide our study of astromaterials into hard and soft and discuss one example of each. The hard as- tromaterial discussed here is a crystalline lattice, such as the Coulomb crystals in the interior of cold white dwarfs and in the crust of neutron stars, while the soft astro- material is nuclear pasta found in the inner crusts of neutron stars. In particular, we discuss how molecular dynamics simulations have been used to calculate the properties of astromaterials to interpret observations of white dwarfs and neutron stars. Coulomb crystals are studied to understand how compact stars freeze. Their incredible strength may make crust \mountains" on rotating neutron stars a source for gravitational waves that the Laser Interferometer Gravitational-Wave Observatory (LIGO) may detect. Nu- clear pasta is expected near the base of the neutron star crust at densities of 1014 g/cm3. Competition between nuclear attraction and Coulomb repulsion rearranges neutrons and protons into complex non-spherical shapes such as sheets (lasagna) or tubes (spaghetti). Semi-classical molecular dynamics simulations of nuclear pasta have been used to study these phases and calculate their transport properties such as neutrino opacity, thermal conductivity, and electrical conductivity.
    [Show full text]
  • Phase Transitions and the Casimir Effect in Neutron Stars
    University of Tennessee, Knoxville TRACE: Tennessee Research and Creative Exchange Masters Theses Graduate School 12-2017 Phase Transitions and the Casimir Effect in Neutron Stars William Patrick Moffitt University of Tennessee, Knoxville, [email protected] Follow this and additional works at: https://trace.tennessee.edu/utk_gradthes Part of the Other Physics Commons Recommended Citation Moffitt, William Patrick, "Phase Transitions and the Casimir Effect in Neutron Stars. " Master's Thesis, University of Tennessee, 2017. https://trace.tennessee.edu/utk_gradthes/4956 This Thesis is brought to you for free and open access by the Graduate School at TRACE: Tennessee Research and Creative Exchange. It has been accepted for inclusion in Masters Theses by an authorized administrator of TRACE: Tennessee Research and Creative Exchange. For more information, please contact [email protected]. To the Graduate Council: I am submitting herewith a thesis written by William Patrick Moffitt entitled "Phaser T ansitions and the Casimir Effect in Neutron Stars." I have examined the final electronic copy of this thesis for form and content and recommend that it be accepted in partial fulfillment of the requirements for the degree of Master of Science, with a major in Physics. Andrew W. Steiner, Major Professor We have read this thesis and recommend its acceptance: Marianne Breinig, Steve Johnston Accepted for the Council: Dixie L. Thompson Vice Provost and Dean of the Graduate School (Original signatures are on file with official studentecor r ds.) Phase Transitions and the Casimir Effect in Neutron Stars A Thesis Presented for the Master of Science Degree The University of Tennessee, Knoxville William Patrick Moffitt December 2017 Abstract What lies at the core of a neutron star is still a highly debated topic, with both the composition and the physical interactions in question.
    [Show full text]
  • Arxiv:1707.04966V3 [Astro-Ph.HE] 17 Dec 2017
    From hadrons to quarks in neutron stars: a review Gordon Baym,1, 2, 3 Tetsuo Hatsuda,2, 4, 5 Toru Kojo,6, 1 Philip D. Powell,1, 7 Yifan Song,1 and Tatsuyuki Takatsuka5, 8 1Department of Physics, University of Illinois at Urbana-Champaign, 1110 W. Green Street, Urbana, Illinois 61801, USA 2iTHES Research Group, RIKEN, Wako, Saitama 351-0198, Japan 3The Niels Bohr International Academy, The Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, DK-2100 Copenhagen Ø, Denmark 4iTHEMS Program, RIKEN, Wako, Saitama 351-0198, Japan 5Theoretical Research Division, Nishina Center, RIKEN, Wako 351-0198, Japan 6Key Laboratory of Quark and Lepton Physics (MOE) and Institute of Particle Physics, Central China Normal University, Wuhan 430079, China 7Lawrence Livermore National Laboratory, 7000 East Ave., Livermore, CA 94550 8Iwate University, Morioka 020-8550, Japan (Dated: December 19, 2017) In recent years our understanding of neutron stars has advanced remarkably, thanks to research converging from many directions. The importance of understanding neutron star behavior and structure has been underlined by the recent direct detection of gravitational radiation from merging neutron stars. The clean identification of several heavy neutron stars, of order two solar masses, challenges our current understanding of how dense matter can be sufficiently stiff to support such a mass against gravitational collapse. Programs underway to determine simultaneously the mass and radius of neutron stars will continue to constrain and inform theories of neutron star interiors. At the same time, an emerging understanding in quantum chromodynamics (QCD) of how nuclear matter can evolve into deconfined quark matter at high baryon densities is leading to advances in understanding the equation of state of the matter under the extreme conditions in neutron star interiors.
    [Show full text]
  • The Nuclear Physics of Neutron Stars
    The Nuclear Physics of Neutron Stars Exotic Beam Summer School 2015 Florida State University Tallahassee, FL - August 2015 J. Piekarewicz (FSU) Neutron Stars EBSS 2015 1 / 21 My FSU Collaborators My Outside Collaborators Genaro Toledo-Sanchez B. Agrawal (Saha Inst.) Karim Hasnaoui M. Centelles (U. Barcelona) Bonnie Todd-Rutel G. Colò (U. Milano) Brad Futch C.J. Horowitz (Indiana U.) Jutri Taruna W. Nazarewicz (MSU) Farrukh Fattoyev N. Paar (U. Zagreb) Wei-Chia Chen M.A. Pérez-Garcia (U. Raditya Utama Salamanca) P.G.- Reinhard (U. Erlangen-Nürnberg) X. Roca-Maza (U. Milano) D. Vretenar (U. Zagreb) J. Piekarewicz (FSU) Neutron Stars EBSS 2015 2 / 21 S. Chandrasekhar and X-Ray Chandra White dwarfs resist gravitational collapse through electron degeneracy pressure rather than thermal pressure (Dirac and R.H. Fowler 1926) During his travel to graduate school at Cambridge under Fowler, Chandra works out the physics of the relativistic degenerate electron gas in white dwarf stars (at the age of 19!) For masses in excess of M =1:4 M electrons becomes relativistic and the degeneracy pressure is insufficient to balance the star’s gravitational attraction (P n5=3 n4=3) ∼ ! “For a star of small mass the white-dwarf stage is an initial step towards complete extinction. A star of large mass cannot pass into the white-dwarf stage and one is left speculating on other possibilities” (S. Chandrasekhar 1931) Arthur Eddington (1919 bending of light) publicly ridiculed Chandra’s on his discovery Awarded the Nobel Prize in Physics (in 1983 with W.A. Fowler) In 1999, NASA lunches “Chandra” the premier USA X-ray observatory J.
    [Show full text]
  • Beyond Nuclear Pasta: Phase Transitions and Neutrino Opacity of Non-Traditional Pasta
    Beyond Nuclear Pasta: Phase Transitions and Neutrino Opacity of Non-Traditional Pasta P. N. Alcain, P. A. Gim´enezMolinelli and C. O. Dorso Departamento de F´ısica, FCEyN, UBA and IFIBA, Conicet, Pabell´on1, Ciudad Universitaria, 1428 Buenos Aires, Argentina and IFIBA-CONICET (Dated: October 8, 2018) In this work, we focus on different length scales within the dynamics of nucleons in conditions according to the neutron star crust, with a semiclassical molecular dynamics model, studying isospin symmetric matter at subsaturation densities. While varying the temperature, we find that a solid- liquid phase transition exists, that can be also characterized with a morphology transition. For higher temperatures, above this phase transition, we study the neutrino opacity, and find that in the liquid phase, the scattering of low momenta neutrinos remain high, even though the morphology of the structures differ significatively from those of the traditional nuclear pasta. PACS numbers: PACS 24.10.Lx, 02.70.Ns, 26.60.Gj, 21.30.Fe I. INTRODUCTION static structure factor of nuclear pasta [5]: dσ dσ A neutron star has a radius of approximately 10 km = × S(q): dΩ dΩ uniform and a mass of a few solar masses. According to current models [1], the structure neutron star can be roughly de- And since neutron stars cooling is associated with neu- scribed as composed in two parts: the crust, of 1.5 km trino emission from the core, the interaction between the thick and a density of up to half the normal nuclear den- neutrinos and the particular structure of the crust would dramatically affect the thermal history of young neutron sity ρ0, with the structures known as nuclear pasta [2]; and the core, where the structure is still unknown and stars.
    [Show full text]
  • Nuclear Pasta with a Touch of Quantum Towards Fully Antisymmetrised Dynamics for Bulk Fermion Systems
    Nuclear Pasta with a touch of Quantum Towards fully antisymmetrised dynamics for bulk fermion systems Klaas Vantournhout Proefschrift ingediend tot het behalen van de academische graad van Doctor in de Wetenschappen: Fysica Academiejaar 2009-2010 Promotor: Prof. Dr. Jan Ryckebusch Promotor: Prof. Dr. Natalie Jachowicz THIS PAGE INTENTIONALLY LEFT BLANK voor mijn ouders Philosophy is written in this grand book – I mean the universe – which stands continually open to our gaze, but it cannot be understood unless one first learns to comprehend the language and interpret the characters in which it is written. It is written in the language of mathematics, and its characters are triangles, circles, and other geometrical figures, without which it is humanly impossible to understand a single word of it; without these, one is wandering around in a dark labyrinth. Galileo Galilei, Il Saggiatore (The Assayer, 1623) THIS PAGE INTENTIONALLY LEFT BLANK DANKWOORD Mijn proefschrift, 160 bladzijden! Een pad doorheen de werelden van fysica, wiskunde en sterrenkunde; geplaveid met continue, binaire en kwantummechanische tegels. Een zes jaar durende queeste naar vergelijkingen. Bloed, zweet en tranen heeft het gekost. Lezen, schrijven, herschrijven, beschrijven, verschrijven, verwijderen, vergis- sen, verfrommelen, ontfrommelen, ontrafelen, ontwarren, verwarren, vertalen, her- talen en ook een beetje euforie . de eerste print. Een eindwerk schrijf je niet alleen. Velen hebben bijgedragen aan dit proces. Veel dank ben ik verschuldigd aan mijn beide promotoren, professoren Jan Ryckebusch en Natalie Jachowicz. De wetenschappelijke vrijheid die ik kreeg, jullie steun, vertrouwen en inzicht waren voor deze scriptie heel belangrijk. I especially would like to express my gratitude to Professor Hans Feldmeier and Doctor Thomas Neff.
    [Show full text]
  • Defects in Nuclear Pasta David Reagan David
    Defects in Nuclear Pasta David Reagan David C. J. Horowitz, Indiana University! Computational Challenges … , Stockholm, Sep. 2014 Nuclear Pasta • Nuclear matter, at somewhat below !0, forms complex shapes because of competition between short range nuclear attraction and long range Coulomb repulsion —> “Coulomb frustration”.! • Gravitational collapse during SN increases density: nuclei -> nuclear pasta -> uniform nuclear matter -> neutron star. ! • Nuclear pasta expected in neutron stars at base of crust about 1 km below surface at ~1/3ρ0.! • During NS merger tidal excitation decreases density: uniform nuclear matter -> nuclear pasta -> nuclei + n -> r-process.! • Quantum density functional calculations. [Sergey Postnikov, Irina Sagert]! • Semiclassical molecular dynamics model: -r2/" -r2/2" -r/# v(r)=a e + bij e + eiej e /r Parameters of short range interaction fit to binding E and density of nuclear matter. [Andre Schneider] NS#Born#in#Core#Collapse#Supernovae# Core#of#massive# star#collapses#to# Envelope form#proto5 neutron#star.######νs# form#neutron#star# energizes#shock# ν" that#ejects#outer# 90%#of#star.# July 5, 1054 Shock Crab nebula Neutrino Crab#Pulsar# Sphere Proto-neutron star: hot, e rich Hubble#ST# Audio: Jordal Bank Neutron stars • Mass ~1.4 Msun, Radius ~10 km • Solid crust ~ 1km thick over liquid (outer) core of neutron rich matter. • Possible exotic phase in center: de-confined quark matter, strange matter, meson condensates, color superconductor... • Structure determined by Equation of State (pressure vs density) of
    [Show full text]
  • The Crust of Merging Neutron Stars
    The Crust of Merging Neutron Stars 208Pb • C. J. Horowitz, Indiana University • INT, July 2014 1 The crust of merging neutron stars • PREX uses parity violating electron scattering to accurately measure the neutron radius of 208Pb. • Merger nucleosynthesys and nuclear pasta. • MD simulations of nuclear pasta. • MD simulations of ions and 208Pb strength of neutron star (outer) Brian Alder crust. 2 Parity Violation Isolates Neutrons • In Standard Model Z0 boson • Apv from interference of photon couples to the weak charge. and Z0 exchange. In Born approximation • Proton weak charge is small: G Q2 F (Q2) p − 2Θ ≈ ! F √ W ! QW =1 4sin W 0.05 Apv = 2 2πα 2 Fch(Q ) • Neutron weak charge is big: ! ! sin(Qr) ! Qn = −1 F (Q2)= d3r ρ (r) W W ! Qr W • Weak interactions, at low Q2, ! probe neutrons. • Model independently map out distribution of weak charge in a • Parity violating asymmetry Apv is nucleus. cross section difference for positive and negative helicity • Electroweak reaction free electrons from most strong dσ/dΩ+ − dσ/dΩ− interaction uncertainties. Apv = ! dσ/dΩ dσ/dΩ− + + – Donnelly, Dubach, Sick first suggested PV to measure neutrons. 3 208 Pb PREX Spokespersons K. Kumar R. Michaels K. Paschke P. Souder G. Urciuoli PREX measures how much neutrons stick out past protons (neutron skin). • 4 PREX in Hall A at JLab Spectometers Lead Foil Target Pol. Source Hall A CEBAF R. Michaels 5 MSU Seminar, Nov 2008 First PREX result and future plans •At Jefferson Laboratory, 1.05 GeV electrons elastically scattered from thick 208Pb foil. PRL 108, 112502, PRC 85, 032501 • A PV=0.66 ±0.06(stat) ±0.014(sym) ppm •Neutron skin thickness: +0.16 Rn-Rp=0.33 -0.18 fm •Experiment achieved systematic error goals.
    [Show full text]
  • Unified Equation of State for Neutron Stars on a Microscopic Basis
    A&A 584, A103 (2015) Astronomy DOI: 10.1051/0004-6361/201526642 & c ESO 2015 Astrophysics Unified equation of state for neutron stars on a microscopic basis B. K. Sharma1,2, M. Centelles1,X.Viñas1,M.Baldo3, and G. F. Burgio3 1 Departament d’Estructura i Constituents de la Matèria and Institut de Ciències del Cosmos (ICC), Facultat de Física, Universitat de Barcelona, Diagonal 645, 08028 Barcelona, Spain 2 Department of Sciences, Amrita Vishwa Vidyapeetham Ettimadai, 641112 Coimbatore, India 3 INFN Sezione di Catania, via Santa Sofia 64, 95123 Catania, Italy e-mail: [email protected] Received 30 May 2015 / Accepted 11 September 2015 ABSTRACT We derive a new equation of state (EoS) for neutron stars (NS) from the outer crust to the core based on modern microscopic calculations using the Argonne v18 potential plus three-body forces computed with the Urbana model. To deal with the inhomogeneous structures of matter in the NS crust, we use a recent nuclear energy density functional that is directly based on the same microscopic calculations, and which is able to reproduce the ground-state properties of nuclei along the periodic table. The EoS of the outer crust requires the masses of neutron-rich nuclei, which are obtained through Hartree-Fock-Bogoliubov calculations with the new functional when they are unknown experimentally. To compute the inner crust, Thomas-Fermi calculations in Wigner-Seitz cells are performed with the same functional. Existence of nuclear pasta is predicted in a range of average baryon densities between 0.067 fm−3 and 0.0825 fm−3, where the transition to the core takes place.
    [Show full text]
  • Pasta-Phase Transitions in the Inner Crust of Neutron Stars∗
    Vol. 10 (2017) Acta Physica Polonica B Proceedings Supplement No 1 PASTA-PHASE TRANSITIONS IN THE INNER CRUST OF NEUTRON STARS∗ X. Viñasa, C. Gonzalez-Boqueraa, B.K. Sharmab, M. Centellesa aDepartament de Física Quàntica i Astrofísica and Institut de Ciències del Cosmos (ICCUB), Facultat de Física Universitat de Barcelona, Martí i Franqués 1, 08028 Barcelona, Spain bDepartment of Sciences, Amrita School of Engineering Amrita Vishwa Vidyapeetham, Amrita University, 641112 Coimbatore, India (Received January 10, 2017) We perform calculations of nuclear pasta phases in the inner crust of neutron stars with the Thomas–Fermi method and the Compressible Liquid Drop Model using the Barcelona–Catania–Paris–Madrid (BCPM) energy density functional and several Skyrme forces. We compare the crust–core transition density estimated from the crust side with the predictions ob- tained from the core by using the thermodynamical and dynamical meth- ods. Finally, the correlation between the crust–core transition density and the slope of the symmetry energy at saturation is briefly analyzed. DOI:10.5506/APhysPolBSupp.10.259 1. Introduction Neutron stars (NS) are among the most fascinating objects of the uni- verse. They are ultradense and compact remnants formed at the last stage of life of massive stars [1,2]. An NS is composed by neutrons, protons, leptons and eventually other exotic particles distributed in a solid crust encompass- ing a dense core in a liquid phase. The full NS is electrically neutral and the distribution of nucleons and leptons is such that the system is maintained locally in β-equilibrium. The structure of the crust consists of clusters of positive charge distributed in a solid lattice embedded in a bath of neu- trons and free electrons.
    [Show full text]
  • Neutron Stars from Crust to Core Within the Quark-Meson Coupling Model
    EPJ Web of Conferences 232, 03001 (2020) https://doi.org/10.1051/epjconf/202023203001 HIAS 2019 Neutron stars from crust to core within the Quark-meson coupling model S. Antic´1;∗, J. R. Stone2;3, and A. W. Thomas1 1CSSM, Department of Physics, University of Adelaide, SA 5005 Australia 2Department of Physics (Astro), University of Oxford, OX1 3RH United Kingdom 3Department of Physics and Astronomy, University of Tennessee, TN 37996 USA Abstract. Recent years continue to be an exciting time for the neutron star physics, providing many new observations and insights to these natural ‘laboratories’ of cold dense matter. To describe them, there are many models on the market but still none that would reproduce all observed and experimental data. The quark-meson coupling model stands out with its natural inclusion of hyperons as dense matter building blocks, and fewer parameters necessary to obtain the nuclear matter equation of state. The latest advances of the QMC model and its application to the neutron star physics will be presented, within which we build the neutron star’s outer crust from finite nuclei up to the neutron drip line. The appearance of different elements and their position in the crust of a neutron star is explored and compared to the predictions of various models, giving the same quality of the results for the QMC model as for the models when the nucleon structure is not taken into account. 1 Introduction an average NS with M ∼ 1:4M the nearly pure neu- tron matter with admixture of protons and electrons is all The neutron stars (NS) are in general very complex ob- there is in a NS core.
    [Show full text]
  • Relativistic Density Functional Theory for Finite Nuclei and Neutron Stars
    Relativistic density functional theory for finite nuclei and neutron stars1 J. Piekarewicz Department of Physics Florida State University Tallahassee, FL 32306 February 6, 2015 Abstract: The main goal of the present contribution is a pedagogical in- troduction to the fascinating world of neutron stars by relying on relativistic density functional theory. Density functional theory provides a powerful|and perhaps unique|framework for the calculation of both the properties of finite nuclei and neutron stars. Given the enormous densities that may be reached in the core of neutron stars, it is essential that such theoretical framework incorporates from the outset the basic principles of Lorentz covariance and special relativity. After a brief historical perspective, we present the necessary details required to compute the equation of state of dense, neutron-rich mat- ter. As the equation of state is all that is needed to compute the structure of neutron stars, we discuss how nuclear physics|particularly certain kind of laboratory experiments|can provide significant constrains on the behavior of neutron-rich matter. arXiv:1502.01559v1 [nucl-th] 5 Feb 2015 1 Contributing chapter to the book \Relativistic Density Functional for Nuclear Structure"; World Scientific Publishing Company (Sin- gapore); Editor Prof. Jie Meng. 1 I. INTRODUCTION The birth of a star is marked by the conversion of hydrogen into helium nuclei (α particles) in their hot dense cores. This thermonuclear reaction is the main source of energy generation during the main stage of stellar evolution and provides the pressure support against gravitational collapse. Once the hydrogen in the stellar core is exhausted, thermonuclear fusion stops and the star contracts.
    [Show full text]