Fermionic molecular dynamics

Size: px
Start display at page:

Download "Fermionic molecular dynamics"

Transcription

1 Fermionic molecular dynamics Can it cook nuclear pasta? K. Vantournhout H. Feldmeier GSI Helmholtzzentrum für Schwerionenforschung GmbH, Darmstadt, Germany NuSTAR seminar April 25, 2012 Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

2 Outline 1 Introduction 2 Recipes for pasta 3 Molecular dynamics for fermions 4 Periodic boundary conditions and FMD 5 Conclusion and outlook Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

3 Introduction : preface Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

4 Introduction : preface Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

5 Introduction : preface Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

6 Introduction : preface Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

7 Introduction : preface Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

8 Introduction : preface Eta Carinae and the bipolar Homunculus Nebula which surrounds the star. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

9 Introduction : preface March 16, 2012 SN2012aw in M95 Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

10 Introduction : preface Crab pulsar Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

11 Introduction : composition of a neutron star Coulomb Crystal Relativistic electron gas (e Z) Mantel Outer Crust Inner Crust Outer Core n+p+e+µ pasta phases Coulomb Crystal Relativistic electron gas dripped neutrons (e Z n)? Inner Core hyperons? -condensate? K-condensate? quark matter? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

12 Introduction : crustal matter and nuclear pastas meatballs spaghetti lasagne ziti Swiss cheese Inner Crust sauce Outer Core Starting from the inner crust, with increasing density three-dimensional spherical nuclear clusters (meatballs) two-dimensional cylindrical tubes of dense matter (spaghetti) one-dimensional slabs interlaid with planar voids (lasagne) two-dimensional cylindrical neutron liquids within the nuclear matter (bucatini) three-dimensional spherical neutron-liquid-bubbles embedded in the nuclear matter (Swiss cheese) transition to the neutron star core, uniform neutron matter (sauce) Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

13 Introduction : what is bucatini? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

14 Introduction : crustal matter and nuclear pastas meatballs spaghetti lasagne ziti Swiss cheese Inner Crust sauce Outer Core Where does the crustal and subnuclear density matter play a role? Supernova physics Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

15 Introduction : crustal matter and nuclear pastas meatballs spaghetti lasagne ziti Swiss cheese Inner Crust sauce Outer Core Where does the crustal and subnuclear density matter play a role? Supernova physics Neutron star cooling Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

16 Introduction : crustal matter and nuclear pastas meatballs spaghetti lasagne ziti Swiss cheese Inner Crust sauce Outer Core Where does the crustal and subnuclear density matter play a role? Supernova physics Neutron star cooling The r-process Neutron star glitches Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

17 Introduction : crustal matter and nuclear pastas meatballs spaghetti lasagne ziti Swiss cheese Inner Crust sauce Outer Core Where does the crustal and subnuclear density matter play a role? Supernova physics Neutron star cooling The r-process Neutron star glitches Gravitational waves Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

18 pastas : is there a relation to Turing instabilities? Tomographic reconstruction of the Belousov-Zhabotinsky reaction incorporated in an water-in-oil aerosol micro-emulsion. Turing patterns form in a reaction-diffusion medium under the condition of a long-ranged reaction inhibition and a short-ranged activation. T. Jr. Bánsági et al., Science 331, pp (2011) Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

19 pasta recipe : pasta alla liquid drop The liquid-drop model K. Nakazato et al., Phys. Rev. Lett. 103, (2009) Semi-empirical model mainly studies the canonical pasta phases does not allow dynamics the hard-coded geometry is studied computationally fairly simple Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

20 pasta recipe : pasta in scatola alla Hartree-Fock The Hartree-Fock model W.G. Newton et al., Phys. Rev. C 79, (2009) mean-field model studies the canonical pasta phases does not allow dynamics the studied geometry strongly depends on the box-geometry computationally very intensive (large configuration space) Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

21 pasta recipe : pasta classica con dinamica molecolare The molecular-dynamics technique G. Watanabe et al., Phys. Rev. Lett. 103, (2009) a quantum-dressed classical many-body technique studies the canonical and intermediate pasta phases allows dynamics the studied geometry is independent of the box-geometry computationally fairly intensive Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

22 pasta recipe : crosta croccante con dinamica molecolare The molecular-dynamics technique is awesome! C. Horowitz et al., Phys. Rev. Lett. 102, (2009) allows laboratory like simulations (e.g. viscosity or stress tensors) allows to study the effect of external probes time-dependent, out-of-equilibrium, thermodynamics, phase transitions,... lacks quantum features (mimics fermion behaviour by a Pauli potential) How can we use MD with correct quantum features? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

23 What is classical molecular dynamics I. Newton Molecular dynamics studies systems of interacting particles by numerically solving the equations of motion. These particles can be point-like, spherical or non-spherical rigid bodies with extra internal degrees of freedom. Time-dependent or static studies are done by solving Newton s or Hamilton s equations of motion (here for point-particles) : W.R. Hamilton q = H p, ṗ = H q. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

24 Fermionic molecular dynamics : the wave function J.C.F. Gauß The single particle fermion wave function is a parametrised quantum-dressed Gaussian wave packet q p = c kp A kp b kp χ kp ζ kp, k x Ab = exp 1 2 (x b) A 1 (x b) The antisymmetrised many-body fermion wave function Q is a Slater determinant of non-orthogonal fermion wave-packets : x 1,..., x A Q = 1 x 1 q 1... x 1 q A det. A! x A q 1... x A q A J.C. Slater Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

25 Fermionic molecular dynamics : the math Adam Ries FMD has a matrix formalism. All relevant fermion physics is contained in the overlap matrix n = [ q p q q ] A pq=1 and its inverse o = n 1. Expectation values are evaluated as : I = Q ˆB I Q Q Q = A q p ˆB I q q o qp pq=1 Time-dependent or static studies are done by means of variational principles on the variables z of trail state Q(z) : Lord Rayleigh ic ż = Q Ĥ Q z Q Q, Q Ĥ ˆT C M Q z = 0. Q Q H. Feldmeier et al., Rev. Mod. Phys. 72, pp (2000) Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

26 Fermionic molecular dynamics : its fascinating! Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

27 How to simulate bulk systems? Finite-sized simulations introduce unwanted surface effects. The more particles, the heavier the computation (N 2(4) ) The less particles, the larger the influence of surface effects. (488/1000 particles) Can we use a small number of particles and avoid surface effects? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

28 How to simulate bulk systems? Finite-sized simulations introduce unwanted surface effects. The more particles, the heavier the computation (N 2(4) ) The less particles, the larger the influence of surface effects. (488/1000 particles) Can we use a small number of particles and avoid surface effects? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

29 How to simulate bulk systems? Finite-sized simulations introduce unwanted surface effects. The more particles, the heavier the computation (N 2(4) ) The less particles, the larger the influence of surface effects. (488/1000 particles) Can we use a small number of particles and avoid surface effects? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

30 Fermionic molecular dynamics : bulk systems Periodic boundary conditions One unit cell contains A single particle states The unit cell is periodically replicated in all directions Which shape of unit cell? What kind of periodicity? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

31 Periodic boundary conditions : Bravais lattices A Bravais lattice is an infinite array of points periodically placed throughout space. a 2 a 1 The lattice appears identical from any point. Each lattice point is referred to through a Bravais lattice vector : R = n 1 a 1 + n 2 a 2 +. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

32 Periodic boundary conditions : Primitive cell A primitive cell tessellates space perfectly without overlap or voids. A primitive cell contains only a single lattice point. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

33 Periodic boundary conditions : Wigner-Seitz cell A Wigner-Seitz cell is a primitive cell preserving the lattice symmetry. It is a region of space that is closer to one lattice point than to any other. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

34 Periodic boundary conditions : Reciprocal lattice The reciprocal lattice is a dual lattice in the Fourier space of the Bravais lattice. Its lattice vectors K satisfy e i K R = 1. Its Wigner-Seitz cell is referred to as the first Brillouin zone. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

35 Periodic Boundary Conditions : Common lattices Square lattice Hexagonal lattice Simple cubic Face-centred cubic Body-centred cubic Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

36 Fermionic molecular dynamics : bulk systems Periodic boundary conditions One unit cell contains A single particle states The unit cell is periodically replicated in all directions Which Bravais lattice? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

37 Fermionic molecular dynamics : bulk systems What about FMD and PBC? Fermions require antisymmetrisation Antisymmetry does not stop at the border of the unit cell FMD must be applied to the infinite periodic system Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

38 Fermionic molecular dynamics : the bulk overlap matrix The matrix formalism of FMD becomes infinite dimensional but has a multi-index block-toeplitz structure that can be exploited. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

39 Fermionic molecular dynamics : the bulk overlap matrix The matrix formalism of FMD becomes infinite dimensional but has a multi-index block-toeplitz structure that can be exploited. N = Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

40 Fermionic molecular dynamics : the bulk overlap matrix The matrix formalism of FMD becomes infinite dimensional but has a multi-index block-toeplitz structure that can be exploited. N = Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

41 Fermionic molecular dynamics : the bulk overlap matrix The matrix formalism of FMD becomes infinite dimensional but has a multi-index block-toeplitz structure that can be exploited. N = Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

42 Fermionic molecular dynamics : the bulk overlap matrix The matrix formalism of FMD becomes infinite dimensional but has a multi-index block-toeplitz structure that can be exploited. N = Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

43 Fermionic molecular dynamics : the bulk overlap matrix The matrix formalism of FMD becomes infinite dimensional but has a multi-index block-toeplitz structure that can be exploited. N = Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

44 Fermionic molecular dynamics : the bulk overlap matrix The matrix formalism of FMD becomes infinite dimensional but has a multi-index block-toeplitz structure that can be exploited. N = Expectation values become finite-dimensional ρ,i = 1 A I,pq (k) qp (k)dk, V BZ BZ pq=1 = n P Q,pq = q p ˆT (Q P) q q where a volume-integral reflects the infinite system. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

45 ` ` Results : the one-dimensional free Fermi gas 4 (a) a/ ` = (b) a/ ` = ρ x 2 ρ kkf (c) a/ ` = (d) a/ ` = 1.0 k F = π/ ` ρ x 1.0 ρ kkf scaled position x/ ` scaled momentum k / k F For small overlap : the particles behave as distinguishable particles For large overlap : the particles reproduce free Fermi gas behaviour Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

46 Results : understanding antisymmetry and PBC Coordinate density ρv W Z Momentum density ρv BZ y/l k y /k l x/l k x /k l 0.0 Gaussian width : a = l/5 without periodicity without antisymmetry Densities in momentum and coordinate space are summed. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

47 Results : understanding antisymmetry and PBC Coordinate density ρv W Z 0.5 Momentum density ρv BZ y/l k y /k l x/l k x /k l 0.0 Gaussian width : a = l/5 without periodicity with antisymmetry Antisymmetry enables the Pauliexclusion principle. It seems like particles are pushed away. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

48 Results : understanding antisymmetry and PBC Coordinate density ρv W Z Momentum density ρv BZ y/l k y /k l x/l k x /k l 0.0 Gaussian width : a = l/5 without periodicity without antisymmetry Densities in momentum and coordinate space are summed. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

49 Results : understanding antisymmetry and PBC Coordinate density ρv W Z Momentum density ρv BZ y/l k y /k l x/l k x /k l 0.0 Gaussian width : a = l/5 with periodicity without antisymmetry Densities in momentum and coordinate space are summed and periodic corrections added in the spatial density. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

50 Results : understanding antisymmetry and PBC Coordinate density ρv W Z 0.5 Momentum density ρv BZ y/l k y /k l x/l k x /k l 0.0 Gaussian width : a = l/5 with periodicity with antisymmetry Antisymmetry is a long-range manybody correlation. Both densities tend to that of a free fermion system. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

51 Results : the two-dimensional free Fermi gas Minimal kinetic energy for N fermions per unit-cell. ρ = 1/fm 2. With increasing overlap, the system moves towards the results of the free Fermi gas, however there is a strong particle dependence. What is the origin of this seemingly erratic behaviour? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

52 Where is Wally? Simulations often introduce constraints to overcome certain problems. But how do these constraints influence the physics of the system under study? Periodic boundary conditions are such a restriction. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

53 Where is Wally? Where is my box? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

54 The effect of periodic boundary conditions Coordinate space Periodic boundary conditions introduce a discrete translation symmetry into the system as a result of the lattice (box). By construction, the simulation cannot leave this symmetry. To uphold the lattice symmetry, Bragg planes and Brillouin zones define the Fermi surface and replace the spherical Fermi distribution. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

55 The effect of periodic boundary conditions Momentum space Periodic boundary conditions introduce a discrete translation symmetry into the system as a result of the lattice (box). By construction, the simulation cannot leave this symmetry. To uphold the lattice symmetry, Bragg planes and Brillouin zones define the Fermi surface and replace the spherical Fermi distribution. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

56 The effect of periodic boundary conditions Momentum space Periodic boundary conditions introduce a discrete translation symmetry into the system as a result of the lattice (box). By construction, the simulation cannot leave this symmetry. To uphold the lattice symmetry, Bragg planes and Brillouin zones define the Fermi surface and replace the spherical Fermi distribution. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

57 The effect of periodic boundary conditions Momentum space Periodic boundary conditions introduce a discrete translation symmetry into the system as a result of the lattice (box). By construction, the simulation cannot leave this symmetry. To uphold the lattice symmetry, Bragg planes and Brillouin zones define the Fermi surface and replace the spherical Fermi distribution. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

58 The effect of periodic boundary conditions Momentum space Periodic boundary conditions introduce a discrete translation symmetry into the system as a result of the lattice (box). By construction, the simulation cannot leave this symmetry. To uphold the lattice symmetry, Bragg planes and Brillouin zones define the Fermi surface and replace the spherical Fermi distribution. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

59 The effect of periodic boundary conditions Momentum space Periodic boundary conditions introduce a discrete translation symmetry into the system as a result of the lattice (box). By construction, the simulation cannot leave this symmetry. To uphold the lattice symmetry, Bragg planes and Brillouin zones define the Fermi surface and replace the spherical Fermi distribution. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

60 The effect of periodic boundary conditions Momentum space Periodic boundary conditions introduce a discrete translation symmetry into the system as a result of the lattice (box). By construction, the simulation cannot leave this symmetry. To uphold the lattice symmetry, Bragg planes and Brillouin zones define the Fermi surface and replace the spherical Fermi distribution. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

61 The effect of periodic boundary conditions Momentum space Periodic boundary conditions introduce a discrete translation symmetry into the system as a result of the lattice (box). By construction, the simulation cannot leave this symmetry. To uphold the lattice symmetry, Bragg planes and Brillouin zones define the Fermi surface and replace the spherical Fermi distribution. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

62 The effect of periodic boundary conditions Momentum space Periodic boundary conditions introduce a discrete translation symmetry into the system as a result of the lattice (box). By construction, the simulation cannot leave this symmetry. To uphold the lattice symmetry, Bragg planes and Brillouin zones define the Fermi surface and replace the spherical Fermi distribution. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

63 The effect of periodic boundary conditions Momentum space Periodic boundary conditions introduce a discrete translation symmetry into the system as a result of the lattice (box). By construction, the simulation cannot leave this symmetry. To uphold the lattice symmetry, Bragg planes and Brillouin zones define the Fermi surface and replace the spherical Fermi distribution. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

64 Results : the two-dimensional free Fermi gas Spatial and momentum distribution for minimal kinetic energy y/l 0.0 k y /k l y/l 0.0 k y /k l y/l x/l ρ = 1/fm 2, square lattice x/l k y /k l k x /k l k x /k l Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

65 Results : the two-dimensional free Fermi gas Spatial and momentum distribution for minimal kinetic energy y/l y/l y/l x/l ρ = 1/fm 2, hexagonal lattice x/l k y /k l k y /k l k y /k l k x /k l k x /k l Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

66 Results : the two-dimensional free Fermi gas Brillouin zones are the limit for large overlap, and become more spherical with increasing N. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

67 Conclusion and outlook Conclusion : FMD can study bulk matter using PBC It reproduces free fermion behaviour The effect of PBC is well understood Outlook : We are finally on the stage where we can investigate the neutron star s crust. Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

68 THE ANTISYMMETRY OF A WAVE FUNCTION IS A LONG-RANGE EFFECT, TREAT IT WITH RESPECT! THE GEOMETRY OF THE LATTICE INFLUENCES THE PHYSICS AND SHOULD BE TREATED AS AN ACTIVE PARAMETER IN SIMULATIONS.

69 Introduction : composition of a neutron star Coulomb Crystal Relativistic electron gas (e Z) Mantel Outer Crust Inner Crust Outer Core n+p+e+µ pasta phases Coulomb Crystal Relativistic electron gas dripped neutrons (e Z n)? Inner Core hyperons? -condensate? K-condensate? quark matter? Klaas Vantournhout (GSI Darmstadt) Can we cook pasta with FMD? NuSTAR / 33

70 THIS SLIDE INTENTIONALLY LEFT BLANK

Structure of light hypernuclei in the framework of Fermionic Molecular Dynamics

Structure of light hypernuclei in the framework of Fermionic Molecular Dynamics 1 Structure of light hypernuclei in the framework of Fermionic Molecular Dynamics Martin Schäfer, Jiří Mareš Nuclear Physics Institute, Řež, Czech Republic H. Feldmeier, T. Neff GSI Helmholtzzentrum für

More information

Impact of Terrestrial Facilities on the Structure of the Neutron Star Crust

Impact of Terrestrial Facilities on the Structure of the Neutron Star Crust Impact of Terrestrial Facilities on the Structure of the Neutron Star Crust Jorge Piekarewicz Florida State University The Neutron Star Crust and Surface (INT - June, 2007) My Collaborators: C.J. Horowitz,

More information

Nuclear Structure for the Crust of Neutron Stars

Nuclear Structure for the Crust of Neutron Stars Nuclear Structure for the Crust of Neutron Stars Peter Gögelein with Prof. H. Müther Institut for Theoretical Physics University of Tübingen, Germany September 11th, 2007 Outline Neutron Stars Pasta in

More information

Clusters in Dense Matter and the Equation of State

Clusters in Dense Matter and the Equation of State Clusters in Dense Matter and the Equation of State Excellence Cluster Universe, Technische Universität München GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt in collaboration with Gerd Röpke

More information

E. Fermi: Notes on Thermodynamics and Statistics (1953))

E. Fermi: Notes on Thermodynamics and Statistics (1953)) E. Fermi: Notes on Thermodynamics and Statistics (1953)) Neutron stars below the surface Surface is liquid. Expect primarily 56 Fe with some 4 He T» 10 7 K ' 1 KeV >> T melting ( 56 Fe) Ionization: r Thomas-Fermi

More information

Nuclear equation of state with realistic nuclear forces

Nuclear equation of state with realistic nuclear forces Nuclear equation of state with realistic nuclear forces Hajime Togashi (RIKEN) Collaborators: M. Takano, K. Nakazato, Y. Takehara, S. Yamamuro, K. Sumiyoshi, H. Suzuki, E. Hiyama 1:Introduction Outline

More information

Self-Consistent Equation of State for Hot Dense Matter: A Work in Progress

Self-Consistent Equation of State for Hot Dense Matter: A Work in Progress Self-Consistent Equation of State for Hot Dense Matter: A Work in Progress W.G.Newton 1, J.R.Stone 1,2 1 University of Oxford, UK 2 Physics Division, ORNL, Oak Ridge, TN Outline Aim Self-consistent EOS

More information

Lecture 11 - Phonons II - Thermal Prop. Continued

Lecture 11 - Phonons II - Thermal Prop. Continued Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon

More information

New Frontiers in Nuclear Structure Theory

New Frontiers in Nuclear Structure Theory New Frontiers in Nuclear Structure Theory From Realistic Interactions to the Nuclear Chart Robert Roth Institut für Kernphysik Technical University Darmstadt Overview Motivation Nucleon-Nucleon Interactions

More information

Translation Symmetry, Space Groups, Bloch functions, Fermi energy

Translation Symmetry, Space Groups, Bloch functions, Fermi energy Translation Symmetry, Space Groups, Bloch functions, Fermi energy Roberto Orlando and Silvia Casassa Università degli Studi di Torino July 20, 2015 School Ab initio Modelling of Solids (UniTo) Symmetry

More information

User s Guide for Neutron Star Matter EOS

User s Guide for Neutron Star Matter EOS User s Guide for Neutron Star Matter EOS CQMC model within RHF approximation and Thomas-Fermi model Tsuyoshi Miyatsu (Tokyo Univ. of Sci.) Ken ichiro Nakazato (Kyushu University) May 1 2016 Abstract This

More information

Quantum simula+ons of nuclear pasta

Quantum simula+ons of nuclear pasta Quantum simula+ons of nuclear pasta William Newton, Sarah Cantu, Mike Gearheart, Farrukh Fa=oyev, Bao-An Li Texas A&M University-Commerce Jirina Rikovska Stone, Helena Pais, Alex Kaltenborn University

More information

Symmetry energy of dilute warm nuclear matter

Symmetry energy of dilute warm nuclear matter Symmetry energy of dilute warm nuclear matter J. B. Natowitz, G. Röpke, 1 S. Typel, 2,3 D. Blaschke, 4, 5 A. Bonasera, 6 K. Hagel, T. Klähn, 4, 7 S. Kowalski, L. Qin, S. Shlomo, R. Wada, and H. H. Wolter

More information

Minimal Surfaces in Nuclear Pasta with the Time-Dependent Hartree-Fock Approach

Minimal Surfaces in Nuclear Pasta with the Time-Dependent Hartree-Fock Approach Minimal Surfaces in Nuclear Pasta with the Time-Dependent Hartree-Fock Approach Bastian Schuetrumpf Institut für Theoretische Physik, Universität Frankfurt, D-60438 Frankfurt, Germany E-mail: schuetrumpf@th.physik.uni-frankfurt.de

More information

Electronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch

Electronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3

More information

Crust-core transitions in neutron stars revisited

Crust-core transitions in neutron stars revisited Crust-core transitions in neutron stars revisited X. Viñas a, C. González-Boquera a, B.K. Sharma a,b M. Centelles a a Departament de Física Quàntica i Astrofísica and Institut de Ciències del Cosmos, Universitat

More information

Calculation of realistic electrostatic potentials in star crusts

Calculation of realistic electrostatic potentials in star crusts Calculation of realistic electrostatic potentials in star crusts Claudio Ebel 1 2 November 3, 212 1 Strong electric fields induced on a sharp stellar boundary I.N.Mishustin, C.Ebel, and W.Greiner, J.Phys.G

More information

Solid State Physics Lecture 3 Diffraction and the Reciprocal Lattice (Kittel Ch. 2)

Solid State Physics Lecture 3 Diffraction and the Reciprocal Lattice (Kittel Ch. 2) Solid State Physics 460 - Lecture 3 Diffraction and the Reciprocal Lattice (Kittel Ch. 2) Diffraction (Bragg Scattering) from a powder of crystallites - real example of image at right from http://www.uni-wuerzburg.de/mineralogie/crystal/teaching/pow.html

More information

Equation-of-State of Nuclear Matter with Light Clusters

Equation-of-State of Nuclear Matter with Light Clusters Equation-of-State of Nuclear Matter with Light Clusters rmann Wolter Faculty of Physics, University of Munich, D-878 Garching, Germany E-mail: hermann.wolter@lmu.de The nuclear equation-of-state (EoS)

More information

Introduction to Dense Matter. C. J. Pethick (U. of Copenhagen and NORDITA)

Introduction to Dense Matter. C. J. Pethick (U. of Copenhagen and NORDITA) Introduction to Dense Matter C. J. Pethick (U. of Copenhagen and NORDITA) Astro-Solids, Dense Matter, and Gravitational Waves INT, Seattle, April 16, 2018 Bottom lines Exciting time for neutron star studies:

More information

Electrons in a weak periodic potential

Electrons in a weak periodic potential Electrons in a weak periodic potential Assumptions: 1. Static defect-free lattice perfectly periodic potential. 2. Weak potential perturbative effect on the free electron states. Perfect periodicity of

More information

Neutron star structure explored with a family of unified equations of state of neutron star matter

Neutron star structure explored with a family of unified equations of state of neutron star matter Neutron star structure explored with a family of unified equations of state of neutron star matter Department of Human Informatics, ichi Shukutoku University, 2-9 Katahira, Nagakute, 48-1197, Japan E-mail:

More information

Thermodynamic properties of nuclear pasta in neutron star crusts

Thermodynamic properties of nuclear pasta in neutron star crusts Thermodynamic properties of nuclear pasta in neutron star crusts Gentaro Watanabe a, Kei Iida a,b, Katsuhiko Sato a,c a Department of Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033,

More information

Graphene and Planar Dirac Equation

Graphene and Planar Dirac Equation Graphene and Planar Dirac Equation Marina de la Torre Mayado 2016 Marina de la Torre Mayado Graphene and Planar Dirac Equation June 2016 1 / 48 Outline 1 Introduction 2 The Dirac Model Tight-binding model

More information

Phase Transitions and the Casimir Effect in Neutron Stars

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

More information

CHEM6085: Density Functional Theory

CHEM6085: Density Functional Theory Lecture 11 CHEM6085: Density Functional Theory DFT for periodic crystalline solids C.-K. Skylaris 1 Electron in a one-dimensional periodic box (in atomic units) Schrödinger equation Energy eigenvalues

More information

equals the chemical potential µ at T = 0. All the lowest energy states are occupied. Highest occupied state has energy µ. For particles in a box:

equals the chemical potential µ at T = 0. All the lowest energy states are occupied. Highest occupied state has energy µ. For particles in a box: 5 The Ideal Fermi Gas at Low Temperatures M 5, BS 3-4, KK p83-84) Applications: - Electrons in metal and semi-conductors - Liquid helium 3 - Gas of Potassium 4 atoms at T = 3µk - Electrons in a White Dwarf

More information

Density Functional Theory: from theory to Applications

Density Functional Theory: from theory to Applications Density Functional Theory: from theory to Applications Uni Mainz May 14, 2012 All electrons vs pseudopotentials Classes of Basis-set Condensed phase: Bloch s th and PBC Hamann-Schlüter-Chiang pseudopotentials

More information

arxiv: v2 [nucl-th] 13 Apr 2012

arxiv: v2 [nucl-th] 13 Apr 2012 Nuclear pasta in supernovae and neutron stars 1 Chapter 1 arxiv:1109.3511v2 [nucl-th] 13 Apr 2012 NUCLEAR PASTA IN SUPERNOVAE AND NEUTRON STARS Gentaro Watanabe 1,2,3 and Toshiki Maruyama 4 1 Asia Pacific

More information

Modeling Materials. Continuum, Atomistic and Multiscale Techniques. gg CAMBRIDGE ^0 TADMOR ELLAD B. HHHHM. University of Minnesota, USA

Modeling Materials. Continuum, Atomistic and Multiscale Techniques. gg CAMBRIDGE ^0 TADMOR ELLAD B. HHHHM. University of Minnesota, USA HHHHM Modeling Materials Continuum, Atomistic and Multiscale Techniques ELLAD B. TADMOR University of Minnesota, USA RONALD E. MILLER Carleton University, Canada gg CAMBRIDGE ^0 UNIVERSITY PRESS Preface

More information

2. Diffraction as a means to determine crystal structure

2. Diffraction as a means to determine crystal structure 2. Diffraction as a means to determine crystal structure Recall de Broglie matter waves: He atoms: [E (ev)] 1/2 = 0.14 / (Å) E 1Å = 0.0196 ev Neutrons: [E (ev)] 1/2 = 0.28 / (Å) E 1Å = 0.0784 ev Electrons:

More information

Author(s) Nakazato, Ken ichiro; Oyamatsu, Kaz. Citation Physical Review Letters (2009), 103.

Author(s) Nakazato, Ken ichiro; Oyamatsu, Kaz. Citation Physical Review Letters (2009), 103. TitleGyroid Phase in Nuclear Pasta Author(s) Nakazato, Ken ichiro; Oyamatsu, Kaz Citation Physical Review Letters (2009), 103 Issue Date 2009-09 URL http://hdl.handle.net/2433/87304 Rightc 2009 The American

More information

Chapter 7 Neutron Stars

Chapter 7 Neutron Stars Chapter 7 Neutron Stars 7.1 White dwarfs We consider an old star, below the mass necessary for a supernova, that exhausts its fuel and begins to cool and contract. At a sufficiently low temperature the

More information

An EOS implementation for astrophyisical simulations

An EOS implementation for astrophyisical simulations Introduction Formalism Neutron Stars CCSN An EOS implementation for astrophyisical simulations A S Schneider 1, L F Roberts 2, C D Ott 1 1 TAPIR, Caltech, Pasadena, CA 2 NSCL, MSU, East Lansing, MI East

More information

Superfluid Heat Conduction in the Neutron Star Crust

Superfluid Heat Conduction in the Neutron Star Crust Superfluid Heat Conduction in the Neutron Star Crust Sanjay Reddy Los Alamos National Lab Collaborators : Deborah Aguilera Vincenzo Cirigliano Jose Pons Rishi Sharma arxiv:0807.4754 Thermal Conduction

More information

Phys 412 Solid State Physics. Lecturer: Réka Albert

Phys 412 Solid State Physics. Lecturer: Réka Albert Phys 412 Solid State Physics Lecturer: Réka Albert What is a solid? A material that keeps its shape Can be deformed by stress Returns to original shape if it is not strained too much Solid structure

More information

Wavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts.

Wavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts. Wavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts. Silvia Casassa Università degli Studi di Torino July 12, 2017 Minnesota Workshop on ab initio MSC Symmetry

More information

2. Diffraction as a means to determine crystal structure

2. Diffraction as a means to determine crystal structure Page 1 of 22 2. Diffraction as a means to determine crystal structure Recall de Broglie matter waves: 2 p h E = where p = 2m λ h 1 E = ( ) 2m λ hc E = hυ = ( photons) λ ( matter wave) He atoms: [E (ev)]

More information

Outline. 1 Relativistic field theory with variable space-time. 3 Extended Hamiltonians in field theory. 4 Extended canonical transformations

Outline. 1 Relativistic field theory with variable space-time. 3 Extended Hamiltonians in field theory. 4 Extended canonical transformations Outline General Relativity from Basic Principles General Relativity as an Extended Canonical Gauge Theory Jürgen Struckmeier GSI Helmholtzzentrum für Schwerionenforschung GmbH, Darmstadt, Germany j.struckmeier@gsi.de,

More information

Chapter 1. Crystal structure. 1.1 Crystal lattices

Chapter 1. Crystal structure. 1.1 Crystal lattices Chapter 1 Crystal structure 1.1 Crystal lattices We will concentrate as stated in the introduction, on perfect crystals, i.e. on arrays of atoms, where a given arrangement is repeated forming a periodic

More information

Nearly Free Electron Gas model - I

Nearly Free Electron Gas model - I Nearly Free Electron Gas model - I Contents 1 Free electron gas model summary 1 2 Electron effective mass 3 2.1 FEG model for sodium...................... 4 3 Nearly free electron model 5 3.1 Primitive

More information

PRINCIPLES OF PHYSICS. \Hp. Ni Jun TSINGHUA. Physics. From Quantum Field Theory. to Classical Mechanics. World Scientific. Vol.2. Report and Review in

PRINCIPLES OF PHYSICS. \Hp. Ni Jun TSINGHUA. Physics. From Quantum Field Theory. to Classical Mechanics. World Scientific. Vol.2. Report and Review in LONDON BEIJING HONG TSINGHUA Report and Review in Physics Vol2 PRINCIPLES OF PHYSICS From Quantum Field Theory to Classical Mechanics Ni Jun Tsinghua University, China NEW JERSEY \Hp SINGAPORE World Scientific

More information

Structure and propermes of nuclear ma3er in compact stars

Structure and propermes of nuclear ma3er in compact stars Structure and propermes of nuclear ma3er in compact stars Toshiki Maruyama (JAEA) Nobutoshi Yasutake (Chiba Inst. of Tech.) Minoru Okamoto (Univ. of Tsukuba & JAEA) Toshitaka Tatsumi (Kyoto Univ.) Low-

More information

Structure and Dynamics : An Atomic View of Materials

Structure and Dynamics : An Atomic View of Materials Structure and Dynamics : An Atomic View of Materials MARTIN T. DOVE Department ofearth Sciences University of Cambridge OXFORD UNIVERSITY PRESS Contents 1 Introduction 1 1.1 Observations 1 1.1.1 Microscopic

More information

Solid State Physics 460- Lecture 5 Diffraction and the Reciprocal Lattice Continued (Kittel Ch. 2)

Solid State Physics 460- Lecture 5 Diffraction and the Reciprocal Lattice Continued (Kittel Ch. 2) Solid State Physics 460- Lecture 5 Diffraction and the Reciprocal Lattice Continued (Kittel Ch. 2) Ewald Construction 2θ k out k in G Physics 460 F 2006 Lect 5 1 Recall from previous lectures Definition

More information

Electronic Structure of Crystalline Solids

Electronic Structure of Crystalline Solids Electronic Structure of Crystalline Solids Computing the electronic structure of electrons in solid materials (insulators, conductors, semiconductors, superconductors) is in general a very difficult problem

More information

Clusters in Low-Density Nuclear Matter

Clusters in Low-Density Nuclear Matter Clusters in Low-Density Nuclear Matter Excellence Cluster Universe, Technische Universität München GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt Helmholtz International Summer School: Nuclear

More information

Mean-field concept. (Ref: Isotope Science Facility at Michigan State University, MSUCL-1345, p. 41, Nov. 2006) 1/5/16 Volker Oberacker, Vanderbilt 1

Mean-field concept. (Ref: Isotope Science Facility at Michigan State University, MSUCL-1345, p. 41, Nov. 2006) 1/5/16 Volker Oberacker, Vanderbilt 1 Mean-field concept (Ref: Isotope Science Facility at Michigan State University, MSUCL-1345, p. 41, Nov. 2006) 1/5/16 Volker Oberacker, Vanderbilt 1 Static Hartree-Fock (HF) theory Fundamental puzzle: The

More information

Neutron-star properties with unified equations of state

Neutron-star properties with unified equations of state Neutron-star properties with unified equations of state Nicolas Chamel in collaboration with J. M. Pearson, S. Goriely, A. F. Fantina Institute of Astronomy and Astrophysics Université Libre de Bruxelles,

More information

Symmetry Energy within the Brueckner-Hartree-Fock approximation

Symmetry Energy within the Brueckner-Hartree-Fock approximation Symmetry Energy within the Brueckner-Hartree-Fock approximation Isaac Vidaña CFC, University of Coimbra International Symposium on Nuclear Symmetry Energy Smith College, Northampton ( Massachusetts) June

More information

The crust-core transition and the stellar matter equation of state

The crust-core transition and the stellar matter equation of state The crust-core transition and the stellar matter equation of state Helena Pais CFisUC, University of Coimbra, Portugal Nuclear Physics, Compact Stars, and Compact Star Mergers YITP, Kyoto, Japan, October

More information

Neutron vs. Quark Stars. Igor Shovkovy

Neutron vs. Quark Stars. Igor Shovkovy Neutron vs. Quark Stars Igor Shovkovy Neutron stars Radius: R 10 km Mass: 1.25M M 2M Period: 1.6 ms P 12 s? Surface magnetic field: 10 8 G B 10 14 G Core temperature: 10 kev T 10 MeV April 21, 2009 Arizona

More information

Electronic Structure Methodology 1

Electronic Structure Methodology 1 Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E

More information

Clusters in Nuclear Matter

Clusters in Nuclear Matter Clusters in Nuclear Matter Excellence Cluster Universe, Technische Universität München GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt in collaboration with Gerd Röpke (Universität Rostock) Thomas

More information

Equation of state for supernovae and neutron stars

Equation of state for supernovae and neutron stars Equation of state for supernovae and neutron stars H. Shen Nankai University, Tianjin, China 申虹南開大学天津中国 In collaboration with H. Toki RCNP, Osaka University, Japan K. Sumiyoshi Numazu College of Technology,

More information

Nuclear equation of state for supernovae and neutron stars

Nuclear equation of state for supernovae and neutron stars Nuclear equation of state for supernovae and neutron stars H. Shen 申虹 In collaboration with Nankai University, Tianjin, China 南開大学 天津 中国 H. Toki RCNP, Osaka University, Japan K. Sumiyoshi Numazu College

More information

Dipole Response of Exotic Nuclei and Symmetry Energy Experiments at the LAND R 3 B Setup

Dipole Response of Exotic Nuclei and Symmetry Energy Experiments at the LAND R 3 B Setup Dipole Response of Exotic Nuclei and Symmetry Energy Experiments at the LAND R 3 B Setup Dominic Rossi for the LAND collaboration GSI Helmholtzzentrum für Schwerionenforschung GmbH D 64291 Darmstadt, Germany

More information

Nuclear symmetry energy and Neutron star cooling

Nuclear symmetry energy and Neutron star cooling Nuclear symmetry energy and Neutron star cooling Yeunhwan Lim 1 1 Daegu University. July 26, 2013 In Collaboration with J.M. Lattimer (SBU), C.H. Hyun (Daegu), C-H Lee (PNU), and T-S Park (SKKU) NuSYM13

More information

Review of lattice EFT methods and connections to lattice QCD

Review of lattice EFT methods and connections to lattice QCD Review of lattice EFT methods and connections to lattice QCD Dean Lee Michigan State University uclear Lattice EFT Collaboration Multi-Hadron Systems from Lattice QCD Institute for uclear Theory Feburary

More information

Quantum Condensed Matter Physics Lecture 5

Quantum Condensed Matter Physics Lecture 5 Quantum Condensed Matter Physics Lecture 5 detector sample X-ray source monochromator David Ritchie http://www.sp.phy.cam.ac.uk/drp2/home QCMP Lent/Easter 2019 5.1 Quantum Condensed Matter Physics 1. Classical

More information

Isovector And Pairing Properties. Of The Gogny Force In The. Context Of Neutron Stars

Isovector And Pairing Properties. Of The Gogny Force In The. Context Of Neutron Stars Isovector And Pairing Properties Of The Gogny Force In The Context Of Neutron Stars Rosh Sellahewa Department of Physics University of Surrey Submitted for the degree of Doctor of Philosophy 17 th September

More information

The Equation of State for Neutron Stars from Fermi Gas to Interacting Baryonic Matter. Laura Tolós

The Equation of State for Neutron Stars from Fermi Gas to Interacting Baryonic Matter. Laura Tolós The Equation of State for Neutron Stars from Fermi Gas to Interacting Baryonic Matter Laura Tolós Outline Outline Neutron Star (I) first observations by the Chinese in 1054 A.D. and prediction by Landau

More information

Quantum Theory of Many-Particle Systems, Phys. 540

Quantum Theory of Many-Particle Systems, Phys. 540 Quantum Theory of Many-Particle Systems, Phys. 540 Questions about organization Second quantization Questions about last class? Comments? Similar strategy N-particles Consider Two-body operators in Fock

More information

Road map (Where are we headed?)

Road map (Where are we headed?) Road map (Where are we headed?) oal: Fairly high level understanding of carrier transport and optical transitions in semiconductors Necessary Ingredients Crystal Structure Lattice Vibrations Free Electron

More information

Section 10 Metals: Electron Dynamics and Fermi Surfaces

Section 10 Metals: Electron Dynamics and Fermi Surfaces Electron dynamics Section 10 Metals: Electron Dynamics and Fermi Surfaces The next important subject we address is electron dynamics in metals. Our consideration will be based on a semiclassical model.

More information

PHYSICS-PH (PH) Courses. Physics-PH (PH) 1

PHYSICS-PH (PH) Courses. Physics-PH (PH) 1 Physics-PH (PH) 1 PHYSICS-PH (PH) Courses PH 110 Physics of Everyday Phenomena (GT-SC2) Credits: 3 (3-0-0) Fundamental concepts of physics and elementary quantitative reasoning applied to phenomena in

More information

NUCLEAR STRUCTURE AB INITIO

NUCLEAR STRUCTURE AB INITIO December, 6:8 WSPC/Trim Size: 9in x 6in for Proceedings master NUCLEAR STRUCTURE AB INITIO H. FELDMEIER AND T. NEFF Gesellschaft für Schwerionenforschung mbh Planckstr., D-69 Darmstadt, Germany E-mail:

More information

Experimental Determination of Crystal Structure

Experimental Determination of Crystal Structure Experimental Determination of Crystal Structure Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, U.S.A. PHYS 624: Introduction to Solid State Physics http://www.physics.udel.edu/~bnikolic/teaching/phys624/phys624.html

More information

Simulating Quantum High Energy Density Plasmas: Approximations According to the Time-Dependent Variational Principle

Simulating Quantum High Energy Density Plasmas: Approximations According to the Time-Dependent Variational Principle Simulating Quantum High Energy Density Plasmas: Approximations According to the Time-Dependent Variational Principle (LANL) with contributions from Andreas Markmann (Yale), Mike Surh (LLNL), Michael S.

More information

Relativistic Viscous Hydrodynamics for Multi-Component Systems with Multiple Conserved Currents

Relativistic Viscous Hydrodynamics for Multi-Component Systems with Multiple Conserved Currents Reference: AM and T. Hirano, arxiv:1003:3087 Relativistic Viscous Hydrodynamics for Multi-Component Systems with Multiple Conserved Currents Akihiko Monnai Department of Physics, The University of Tokyo

More information

Density Functional Theory: from theory to Applications

Density Functional Theory: from theory to Applications Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation

More information

Basics of DFT applications to solids and surfaces

Basics of DFT applications to solids and surfaces Basics of DFT applications to solids and surfaces Peter Kratzer Physics Department, University Duisburg-Essen, Duisburg, Germany E-mail: Peter.Kratzer@uni-duisburg-essen.de Periodicity in real space and

More information

Nuclear structure I: Introduction and nuclear interactions

Nuclear structure I: Introduction and nuclear interactions Nuclear structure I: Introduction and nuclear interactions Stefano Gandolfi Los Alamos National Laboratory (LANL) National Nuclear Physics Summer School Massachusetts Institute of Technology (MIT) July

More information

Contents. 1.1 Prerequisites and textbooks Physical phenomena and theoretical tools The path integrals... 9

Contents. 1.1 Prerequisites and textbooks Physical phenomena and theoretical tools The path integrals... 9 Preface v Chapter 1 Introduction 1 1.1 Prerequisites and textbooks......................... 1 1.2 Physical phenomena and theoretical tools................. 5 1.3 The path integrals..............................

More information

Possibility of hadron-quark coexistence in massive neutron stars

Possibility of hadron-quark coexistence in massive neutron stars Possibility of hadron-quark coexistence in massive neutron stars Tsuyoshi Miyatsu Department of Physics, Soongsil University, Korea July 17, 2015 Nuclear-Astrophysics: Theory and Experiments on 2015 2nd

More information

An effective potential for electron-nucleus scattering in neutrino-pair bremsstrahlung in neutron star crust

An effective potential for electron-nucleus scattering in neutrino-pair bremsstrahlung in neutron star crust Journal of Physics: Conference Series PAPER OPEN ACCESS An effective potential for electron-nucleus scattering in neutrino-pair bremsstrahlung in neutron star crust To cite this article: D D Ofengeim et

More information

Phase transitions in dilute stellar matter. Francesca Gulminelli & Adriana Raduta

Phase transitions in dilute stellar matter. Francesca Gulminelli & Adriana Raduta Phase transitions in dilute stellar matter Francesca Gulminelli & Adriana Raduta LPC Caen, France IFIN Bucharest Supernova remnant and neutron star in Puppis A (ROSAT x-ray) χ 1/2 Τ 10 12 Κ Motivation:

More information

neutron star basics λ. λ

neutron star basics λ. λ neutron star basics A solar mass consists of ~ 10 57 nucleons. If they are separated by typical inter-nucleon distances, what would the radius of the volume containing them be? neutron star basics ( )

More information

α Particle Condensation in Nuclear systems

α Particle Condensation in Nuclear systems Particle Condensation in Nuclear systems A. Tohsaki, H. Horiuchi, G. Röpke, P. Sch. T. Yamada and Y. Funaki -condensation in matter 8 Be and Hoyle state in 12 C -condensate wave function Effective GPE

More information

Dense Matter and Neutrinos. J. Carlson - LANL

Dense Matter and Neutrinos. J. Carlson - LANL Dense Matter and Neutrinos J. Carlson - LANL Neutron Stars and QCD phase diagram Nuclear Interactions Quantum Monte Carlo Low-Density Equation of State High-Density Equation of State Neutron Star Matter

More information

arxiv:astro-ph/ v2 24 Apr 2001

arxiv:astro-ph/ v2 24 Apr 2001 Neutron Star Structure and the Neutron Radius of 208 Pb C. J. Horowitz Nuclear Theory Center and Dept. of Physics, Indiana University, Bloomington, IN 47405 J. Piekarewicz Department of Physics Florida

More information

1 Introduction. 2 The hadronic many body problem

1 Introduction. 2 The hadronic many body problem Models Lecture 18 1 Introduction In the next series of lectures we discuss various models, in particluar models that are used to describe strong interaction problems. We introduce this by discussing the

More information

A Tour to the Stars. Some Compelling Questions

A Tour to the Stars. Some Compelling Questions Some Compelling Questions A Tour to the Stars The chemical elements essential for life: How did they form? Where did they form? Stars: How are they born? How do they live? How do they die? The stellar

More information

First, we need a rapid look at the fundamental structure of superfluid 3 He. and then see how similar it is to the structure of the Universe.

First, we need a rapid look at the fundamental structure of superfluid 3 He. and then see how similar it is to the structure of the Universe. Outline of my talk: First, we need a rapid look at the fundamental structure of superfluid 3 He and then see how similar it is to the structure of the Universe. Then we will look at our latest ideas on

More information

The nature of superfluidity in the cold atomic unitary Fermi gas

The nature of superfluidity in the cold atomic unitary Fermi gas The nature of superfluidity in the cold atomic unitary Fermi gas Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo (AFMC) method The trapped unitary Fermi gas

More information

Molecular Aggregation

Molecular Aggregation Molecular Aggregation Structure Analysis and Molecular Simulation of Crystals and Liquids ANGELO GAVEZZOTTI University of Milano OXFORD UNIVERSITY PRESS Contents PART I FUNDAMENTALS 1 The molecule: structure,

More information

Computational Physics. J. M. Thijssen

Computational Physics. J. M. Thijssen Computational Physics J. M. Thijssen Delft University of Technology CAMBRIDGE UNIVERSITY PRESS Contents Preface xi 1 Introduction 1 1.1 Physics and computational physics 1 1.2 Classical mechanics and statistical

More information

PoS(NIC XII)250. A new equation of state with abundances of all nuclei in core collapse simulations of massive stars

PoS(NIC XII)250. A new equation of state with abundances of all nuclei in core collapse simulations of massive stars A new equation of state with abundances of all nuclei in core collapse simulations of massive stars 1, Kohsuke Sumiyoshi 2, Shoichi Yamada 1,3, Hideyuki Suzuki 4 1 Department of Science and Engineering,

More information

文德华 Department of Physics, South China Univ. of Tech. ( 华南理工大学物理系 )

文德华 Department of Physics, South China Univ. of Tech. ( 华南理工大学物理系 ) Investigation on the oscillation modes of neutron stars 文德华 Department of Physics, South China Univ. of Tech. ( 华南理工大学物理系 ) collaborators Bao-An Li, William Newton, Plamen Krastev Department of Physics

More information

Chapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set

Chapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless

More information

Nuclear structure Anatoli Afanasjev Mississippi State University

Nuclear structure Anatoli Afanasjev Mississippi State University Nuclear structure Anatoli Afanasjev Mississippi State University 1. Nuclear theory selection of starting point 2. What can be done exactly (ab-initio calculations) and why we cannot do that systematically?

More information

Nuclear Physics from Lattice Effective Field Theory

Nuclear Physics from Lattice Effective Field Theory Nuclear Physics from Lattice Effective Field Theory Dean Lee (NCSU/Bonn) work done in collaboration with Evgeny Epelbaum (Bochum) Hermann Krebs (Bochum) Ulf-G. Meißner (Bonn/Jülich) Buḡra Borasoy (now

More information

arxiv:nucl-th/ v1 21 Mar 2001

arxiv:nucl-th/ v1 21 Mar 2001 Relativistic Hartree-Bogoliubov Calculation of Specific Heat of the Inner Crust of Neutron Stars arxiv:nucl-th/5v Mar akuya Nakano and Masayuki Matsuzaki Department of Physics, Kyushu University, Fukuoka

More information

Relativistic Astrophysics Neutron Stars, Black Holes & Grav. W. ... A brief description of the course

Relativistic Astrophysics Neutron Stars, Black Holes & Grav. W. ... A brief description of the course Relativistic Astrophysics Neutron Stars, Black Holes & Grav. Waves... A brief description of the course May 2, 2009 Structure of the Course Introduction to General Theory of Relativity (2-3 weeks) Gravitational

More information

Keble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering

Keble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering Tomi Johnson Keble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering Please leave your work in the Clarendon laboratory s J pigeon hole by 5pm on Monday of

More information

Quantum Theory of Many-Particle Systems, Phys. 540

Quantum Theory of Many-Particle Systems, Phys. 540 Quantum Theory of Many-Particle Systems, Phys. 540 IPM? Atoms? Nuclei: more now Other questions about last class? Assignment for next week Wednesday ---> Comments? Nuclear shell structure Ground-state

More information

Phys 460 Describing and Classifying Crystal Lattices

Phys 460 Describing and Classifying Crystal Lattices Phys 460 Describing and Classifying Crystal Lattices What is a material? ^ crystalline Regular lattice of atoms Each atom has a positively charged nucleus surrounded by negative electrons Electrons are

More information

Theory of mean motion resonances.

Theory of mean motion resonances. Theory of mean motion resonances. Mean motion resonances are ubiquitous in space. They can be found between planets and asteroids, planets and rings in gaseous disks or satellites and planetary rings.

More information

2. Surface geometric and electronic structure: a primer

2. Surface geometric and electronic structure: a primer 2. Surface geometric and electronic structure: a primer 2.1 Surface crystallography 2.1.1. Crystal structures - A crystal structure is made up of two basic elements: lattice + basis Basis: Lattice: simplest

More information