Towards first-principle description of electromagnetic reactions in medium-mass nuclei

Similar documents
A path to lepton-nucleus reactions from first principles

Electromagnetic reactions from few to many-body systems Giuseppina Orlandini

Electromagentic Reactions and Structure of Light Nuclei

Understanding the proton radius puzzle: Nuclear structure corrections in light muonic atoms

Nucleon Nucleon Forces and Mesons

Electron Sca ering off 4 He

Recent developments in the Lorentz. Transform (LIT)) method

Recirculating Electron Beam Photo-converter for Rare Isotope Production CYC2016

The Lorentz Integral Transform (LIT) method and its connections with other approaches

Coupled-cluster computations of weak decays in nuclei

Theory Group Plans and Highlights

Coupled-cluster theory for medium-mass nuclei

Weak decays from coupled cluster computations

(Todays) Progress in coupled cluster compuations of atomic nuclei

Coupled-cluster theory for nuclei

Study of Photo-Production of Radioactive Ion Beams

Charge State Breeding for the. at TRIUMF

Weakly-Bound Systems in Atomic and Nuclear Physics March 2010

Three-nucleon forces and shell structure of neutron-rich Ca isotopes

Current status and challenges of ab-initio computations of nuclei

Probing the Nuclear Symmetry Energy and Neutron Skin from Collective Excitations. N. Paar

Recent Developments and Applications of Integral Transforms in Few- and Many-Body Physics

Evolution Of Shell Structure, Shapes & Collective Modes. Dario Vretenar

Light Nuclei from chiral EFT interactions

Ab initio description of the unbound 7 He

Physics of neutron-rich nuclei

Nuclear structure from chiral-perturbation-theory two- plus three-nucleon interactions

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

Three-cluster dynamics within an ab initio framework

Coupled-cluster computations of neutron-rich nuclei

From few-body to many-body systems

Shell evolution and pairing in calcium isotopes with two- and three-body forces

Few Body Methods in Nuclear Physics - Lecture I

Benchmarks from high-precision mass measurements at TITAN

Nuclear structure and the proton radius puzzle

Shell evolution in neutron rich nuclei

Ab Initio Theory for All Medium-Mass Nuclei

Frontiers in Ab Initio Nuclear Structure Theory. Robert Roth

Nuclear electric dipole moment in the Gaussian expansion method

Electromagnetic Response of Light Nuclei with Integral Transforms

Pushing the boundaries of Penning trap mass spectrometry

Beyond mean-field study on collective vibrations and beta-decay

New Trends in the Nuclear Shell Structure O. Sorlin GANIL Caen

The Nuclear Many-Body Problem. Lecture 2

Theory of neutron-rich nuclei and nuclear radii Witold Nazarewicz (with Paul-Gerhard Reinhard) PREX Workshop, JLab, August 17-19, 2008

New Horizons in Ab Initio Nuclear Structure Theory. Robert Roth

Constraining the symmetry energy based on relativistic point coupling interactions and excitations in finite nuclei

An Electron Linac Photo-Fission Driver for the Rare Isotope Program at TRIUMF

Radiative-capture reactions

Shell model description of dipole strength at low energy

Effective Field Theory for light nuclear systems

Laser spectroscopy studies of neutron-rich nuclei

Inves&ga&on of neutron- rich isotopes: Understanding the r process

Charge exchange reactions and photo-nuclear reactions

Modeling the Atomic Nucleus. Theoretical bag of tricks

Neutron stars at JLAB and the Pb Radius Experiment

Neutron-rich matter and neutrino-matter interactions based on chiral effective field theory

ab-initio alpha-alpha scattering

Ab Initio Electromagnetic Transitions with the IMSRG

arxiv: v1 [nucl-th] 13 Nov 2017

The No Core Shell Model: Its Formulation, Application and Extensions. Bruce R. Barrett University of Arizona,

PREX and CREX. R N from Electroweak Asymmetry in Elastic Electron-Nucleus Scattering. Neutron Skin.

Central density. Consider nuclear charge density. Frois & Papanicolas, Ann. Rev. Nucl. Part. Sci. 37, 133 (1987) QMPT 540

Coupling of giant resonances to soft E1 and E2 modes in 8 B

The No-Core Shell Model

The Nuclear Equation of State

Application of Equation of Motion Phonon Method to Nuclear and Exotic Nuclear Systems

Pb, 120 Sn and 48 Ca from High-Resolution Proton Scattering MSU 2016

Joint ICTP-IAEA Workshop on Nuclear Structure Decay Data: Theory and Evaluation August Introduction to Nuclear Physics - 1

Effective Field Theory for Nuclear Physics! Akshay Vaghani! Mississippi State University!

Ab initio many-body calculations of nuclear scattering and reactions

Three-Nucleon Forces and Masses of Neutron-Rich Nuclei Jason D. Holt

Continuum States in Drip-line Oxygen isotopes

Linking nuclear reactions and nuclear structure to on the way to the drip lines

hg: Chiral Structure of Few-Nucleon Systems

in covariant density functional theory.

Photo Reactions in Few-Body Systems - From Nuclei to Cold-Atoms

Ab initio effective interactions and operators from IM-SRG

Semi-Classical perturbation theory Coulomb only First-order most used

Description of (p, pn) reactions induced by Borromean nuclei

Pairing Correlations in Nuclei on Neutron Drip Line

Can the shell model be truly ab initio? and other questions

Solar Fusion Cross Sections for the pp chain and CNO cycles

Coupled-Cluster Theory. Nuclear Structure

Beyond Schiff: Atomic EDMs from two-photon exchange Satoru Inoue (work with Wick Haxton & Michael Ramsey-Musolf) ACFI Workshop, November 6, 2014

Part III: The Nuclear Many-Body Problem

Introduction to Nuclear Physics

The Shell Model: An Unified Description of the Structure of th

Ab Initio Nuclear Structure Theory with Chiral NN+3N Interactions. Robert Roth

Three-nucleon forces and neutron-rich nuclei

Quantum Theory of Many-Particle Systems, Phys. 540

Neutron Halo in Deformed Nuclei

From EFTs to Nuclei. Thomas Papenbrock. and. CANHP 2015 Research partly funded by the US Department of Energy

In-medium Similarity Renormalization Group for nuclear many-body systems

New Frontiers in Nuclear Structure Theory

Discerning the symmetry energy and neutron star properties from nuclear collective excitations

Direct reactions methodologies for use at fragmentation beam energies

Nature of low-energy dipole states in exotic nuclei

SciDAC project NUCLE lead PI: Joe Carlson (LA PetaApps award lead PI: Jerry Draayer (L. lead PI: James P Vary (I NERSC

Ab initio rotational bands in medium and heavy nuclei

Transcription:

Canada s National Laboratory for Particle and Nuclear Physics Laboratoire national canadien pour la recherche en physique nucléaire et en physique des particules Towards first-principle description of electromagnetic reactions in medium-mass nuclei Nuclear Halo Sonia Bacca Theory Group TRIUMF In collaboration with: N.Barnea, G.Hagen, M.Miorelli, G.Orlandini, T.Papenbrock Elba XIII Workshop on Electron-nucleus scattering, June 23-27, 214, Marciana Marina Electromagnetic Reactions Owned and operated as a joint venture by a consortium of Canadian universities via a contribution through the National Research Council Canada Propriété d un consortium d universités canadiennes, géré en co-entreprise à partir d une contribution administrée par le Conseil national de recherches Canada

Motivations Electromagnetic probes (coupling constant <<1) With the electro-magnetic probe, we can immediately relate the cross section to the transition matrix element of the current operator, thus to the structure of the target itself [De Forest-Walecka, Ann. Phys. 1966] f J µ 2! = q k k µ µ!, q µ µ µ q = k k µ q = ( ω, q) µ P f µ P June 25th 214 Sonia Bacca 2

Motivations Electromagnetic probes (coupling constant <<1) With the electro-magnetic probe, we can immediately relate the cross section to the transition matrix element of the current operator, thus to the structure of the target itself [De Forest-Walecka, Ann. Phys. 1966] f J µ 2! = q k k µ µ!, q µ µ µ q = k k µ q = ( ω, q) µ P f µ P For few-nucleons one can perform exact calculations both for bound and scattering states test the nuclear theory on light nuclei and extend it to heavier mass number June 25th 214 Sonia Bacca 2

Motivations Electromagnetic probes (coupling constant <<1) With the electro-magnetic probe, we can immediately relate the cross section to the transition matrix element of the current operator, thus to the structure of the target itself [De Forest-Walecka, Ann. Phys. 1966] f J µ 2! = q k k µ µ!, q µ µ µ q = k k µ q = ( ω, q) µ P f µ P For few-nucleons one can perform exact calculations both for bound and scattering states test the nuclear theory on light nuclei and extend it to heavier mass number Provide important informations in other fields of physics, where nuclear physics plays a crucial role: - Astrophysics: interactions with nucleonic matter, radiative capture reactions, interactions with nucleonic matter (vector current as em) - Atomic physics (nuclear corrections to atomic levels, etc.) PRL 111, 14342 (213) arxiv:143.7651, to appear on PRC - Particle physics (neutrino-interaction with nuclei) arxiv:146.523 June 25th 214 Sonia Bacca 2

Electromagnetic Reactions Photonuclear Reactions Electric-Dipole Polarizability... 3

Electromagnetic Reactions Photo-nuclear Reactions Reactions resulting from the interaction of a photon with the nucleus For photon energy 15-25 MeV stable nuclei across the periodic table show a wide and large peak 4 35 3 16 O Ahrens et al. σ γ [mb] 25 2 15 1 5 1 2 3 4 5 E γ [MeV] 4

Electromagnetic Reactions Photo-nuclear Reactions Reactions resulting from the interaction of a photon with the nucleus For photon energy 15-25 MeV stable nuclei across the periodic table show a wide and large peak 4 35 3 16 O Ahrens et al. σ γ [mb] 25 2 15 1 5 1 2 3 4 5 E γ [MeV] Giant Dipole Resonance 5

Electromagnetic Reactions Photo-nuclear Reactions Reactions resulting from the interaction of a photon with the nucleus For photon energy 15-25 MeV stable nuclei across the periodic table show a wide and large peak 4 35 3 16 O Ahrens et al. σ γ [mb] 25 2 15 1 5 1 2 3 4 5 E γ [MeV] Giant Dipole Resonance 5

Electromagnetic Reactions 4 35 Photo-nuclear Reactions Reactions resulting from the interaction of a photon with the nucleus For photon energy 15-25 MeV stable nuclei across the periodic table show a wide and large peak 16 O Ahrens et al. Coulomb excitations Inelastic scattering between two charged particles. Can use unstable nuclei as projectiles. Neutron-rich nuclei show fragmented low-lying strength Leistenschneider et al. 3 σ γ [mb] 25 2 15 1 5 1 2 3 4 5 E γ [MeV] Giant Dipole Resonance 5

Electromagnetic Reactions 4 35 Photo-nuclear Reactions Reactions resulting from the interaction of a photon with the nucleus For photon energy 15-25 MeV stable nuclei across the periodic table show a wide and large peak 16 O Ahrens et al. Coulomb excitations Inelastic scattering between two charged particles. Can use unstable nuclei as projectiles. Neutron-rich nuclei show fragmented low-lying strength Leistenschneider et al. 3 σ γ [mb] 25 2 15 1 5 1 2 3 4 5 E γ [MeV] Giant Dipole Resonance Can we give a microscopic explanation of these observations? 5

Ab-initio Theory Tools s 1 r 2 s 2... r 1 s A r A H i = E i i H = T + V NN + V 3N +... High precision two-nucleon potentials: well constraint on NN phase shifts Three nucleon forces: less known, constraint on A>2 observables 6

Ab-initio Theory Tools s 1 r 2 s 2... J µ = J µ N + J µ NN +... r 1 s A m + m r A H i = E i i H = T + V NN + V 3N +... High precision two-nucleon potentials: well constraint on NN phase shifts N N N N two-body currents (or MEC) subnuclear d.o.f. J µ consistent with V r J = i[v, ] Three nucleon forces: less known, constraint on A>2 observables 6

Ab-initio Theory Tools s 1 r 2 s 2... Traditional Nuclear Physics J µ = J µ N + J µ NN +... r 1 s A AV18+UIX,..., J2 m + m r A H i = E i i H = T + V NN + V 3N +... High precision two-nucleon potentials: well constraint on NN phase shifts Three nucleon forces: less known, constraint on A>2 observables Effective Field Theory N 2 LO, N 3 LO... N N N N two-body currents (or MEC) subnuclear d.o.f. J µ consistent with V r J = i[v, ] 6

Ab-initio Theory Tools s 1 r 2 s 2... Traditional Nuclear Physics J µ = J µ N + J µ NN +... r 1 s A AV18+UIX,..., J2 m + m r A H i = E i i H = T + V NN + V 3N +... High precision two-nucleon potentials: well constraint on NN phase shifts Three nucleon forces: less known, constraint on A>2 observables Effective Field Theory N 2 LO, N 3 LO... N N N N two-body currents (or MEC) subnuclear d.o.f. J µ consistent with V r J = i[v, ] f J µ 2 Exact Initial state & Final state in the continuum at different energies and for different A 6

Efros, et al., JPG.: Nucl.Part.Phys. 34 (27) R459 R( ) = L(, )= d where Lorentz Integral Transform Method Reduce the continuum problem to a bound-state problem f 2 f Ô (Ef E ) J µ R( ) ( ) 2 + is obtained solving = 2 < 1 (H E + i ) = J Ô µ Due to imaginary part the solution is unique h i Since is finite, has bound state asymptotic behaviour June 25th 214 Sonia Bacca 7

Efros, et al., JPG.: Nucl.Part.Phys. 34 (27) R459 R( ) = L(, )= d where Lorentz Integral Transform Method Reduce the continuum problem to a bound-state problem f 2 f Ô (Ef E ) J µ R( ) ( ) 2 + is obtained solving = 2 < 1 (H E + i ) = J Ô µ Due to imaginary part the solution is unique h i Since is finite, has bound state asymptotic behaviour inversion L(, ) R( ) The exact final state interaction is included in the continuum rigorously! June 25th 214 Sonia Bacca 7

Efros, et al., JPG.: Nucl.Part.Phys. 34 (27) R459 R( ) = L(, )= d where Lorentz Integral Transform Method Reduce the continuum problem to a bound-state problem f 2 f Ô (Ef E ) J µ R( ) ( ) 2 + is obtained solving = 2 < 1 (H E + i ) = J Ô µ Due to imaginary part the solution is unique h i Since is finite, has bound state asymptotic behaviour inversion L(, ) R( ) The exact final state interaction is included in the continuum rigorously! Solved for A=3,4,6,7 with hyper-spherical harmonics expansions and for A=4 with NCSM June 25th 214 Sonia Bacca 7

Photoabsorption of A=6 with Hyperspherical Harmonics E1 = 4 2 3 RE1 ( ) = ZX (z i Z cm ) i 5 4 6 Li 6 He stable unstable Aumann et al. 6 He at GSI T 1/2 = 86 ms σ γ [mb] 3 2 6 Li 1 6 He S.B. et al. PRL 89 5252 (22) 1 2 3 4 5 6 ω [MeV] No cluster assumption was used in this calculation. Potential: AV4 8

Hyperspherical Harmonics A basis set, mostly used for A=3,4. Challenge to go up to A=7,8 b 1e+5 Ground States 1 # HH states 1 1 A=3 A=4 1 A=6 A=8 1 2 4 6 8 1 12 14 16 18 2 22 24 26 K max (H E + i ) = J Ô µ J µ i For the reactions expanding increase of dimension by at least one order of magnitude if angular momentum is changed 9

Hyperspherical Harmonics A basis set, mostly used for A=3,4. Challenge to go up to A=7,8 b 1e+5 Ground States 1 # HH states 1 1 A=3 A=4 1 A=6 A=8 1 2 4 6 8 1 12 14 16 18 2 22 24 26 K max Total number of states: #HH X #Hyper-radial states 1 6 dense matrix (H E + i ) = J Ô µ J µ i For the reactions expanding increase of dimension by at least one order of magnitude if angular momentum is changed 9

Extension to medium-mass nuclei Develop new many-body methods that can extend the frontiers to heavier and neutron nuclei Coupled Cluster Theory CC future aims CC theory now CC is optimal for closed shell nuclei ( ± 1, ± 2) Uses particle coordinates ( r 1, r 2,..., r A ) = e T ( r 1, r 2,..., r A ) T = T (A) cluster expansion reference SD with any sp states T 1 = T 2 = 1 4 ia X ij,ab t a i a aa i t ab ij a aa b a ja i... T 1 T 2 T 3 CCSD CCSDT CCSD Computational scaling n 2 on 4 u June 25th 214 Sonia Bacca

Extension to medium-mass nuclei Develop new many-body methods that can extend the frontiers to heavier and neutron nuclei Coupled Cluster Theory CC future aims CC theory now CC is optimal for closed shell nuclei ( ± 1, ± 2) Uses particle coordinates ( r 1, r 2,..., r A ) = e T ( r 1, r 2,..., r A ) T = T (A) cluster expansion reference SD with any sp states... CC is a very mature theory for g.s., see e.g. Hagen et al. PRL 11, 9252 (28),PRC 82, 3433 (21) PRL 18, 24251 (212), PRL 19, 3252 (212) experiment R. Roth et al., Phys. Rev. Lett. 19, 5251 (212) S. Binder et al, arxiv:1312.5685 June 25th 214 Sonia Bacca 11

Extension to medium-mass nuclei Develop new many-body methods that can extend the frontiers to heavier and neutron nuclei Coupled Cluster Theory CC future aims CC theory now CC is optimal for closed shell nuclei ( ± 1, ± 2) Uses particle coordinates ( r 1, r 2,..., r A ) = e T ( r 1, r 2,..., r A ) T = T (A) cluster expansion reference SD with any sp states... CC is a very mature theory for g.s., see e.g. Hagen et al. PRL 11, 9252 (28),PRC 82, 3433 (21) PRL 18, 24251 (212), PRL 19, 3252 (212) What about electromagnetic reactions? experiment R. Roth et al., Phys. Rev. Lett. 19, 5251 (212) S. Binder et al, arxiv:1312.5685 June 25th 214 Sonia Bacca 11

LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass D R( ) L(, )= d = h i L RE = h ( ) 2 + 2 ˆL(z) R(z ˆ ) i (H z ) i = J µ i with z = E + + i ( H z ) Ri = i H = e T He T =e T e T ( H z ) ˆR(z ) i = i EoM with source 12

LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass D R( ) L(, )= d = h i L RE = h ( ) 2 + 2 ˆL(z) R(z ˆ ) i (H z ) i = J µ i with z = E + + i ( H z ) Ri = i H = e T He T =e T e T ( H z ) ˆR(z ) i = i EoM with source 12

LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass D R( ) L(, )= d = h i L RE = h ( ) 2 + 2 ˆL(z) R(z ˆ ) i (H z ) i = J µ i with z = E + + i ( H z ) Ri = i H = e T He T =e T e T ( H z ) ˆR(z ) i = i EoM with source 12

LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass D R( ) L(, )= d = h i L RE = h ( ) 2 + 2 ˆL(z) R(z ˆ ) i (H z ) i = J µ i with z = E + + i ( H z ) Ri = i H = e T He T =e T e T ( H z ) ˆR(z ) i = i EoM with source In the CCSD scheme T = T 1 + T 2 ˆR = R + X ia R a i ĉ aĉ i + 1 4 X ˆ ˆ Rˆij ab ĉ aĉ bĉjĉ i +... ijab 12

L(, )= d LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass R( ) ( ) 2 + 2 = h i D L RE = h L(z) R(z ) i ˆ ˆ (H z ) i = J µ i with z = E + + i ( H z ) Ri = i H = e T He T =e T e T ( H z ) ˆR(z ) i = i EoM with source In the CCSD scheme T = T 1 + T 2 ˆR = R + X ia R a i ĉ aĉ i + 1 4 X ˆ ˆ Rˆij ab ĉ aĉ bĉjĉ i +... ijab Validation for 4 He Dipole Response Functions with NN forces from EFT (N 3 LO) R(ω) [mb/mev].5.4.3.2 4 He EIHH Exact HH CCSD.1 2 4 6 8 1 ω [MeV] 12

LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass Extension to Dipole Response Function in 16 O with NN forces derived from EFT (N 3 LO) Convergence in the model space expansion L [fm 2 MeV -2 1-2 ] 4 3 2 1 Γ=1 MeV N max 18 16 14 12 1 8 hω=26 MeV -2 2 4 6 8 1 12 σ [MeV] Good convergence! 13

LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass Extension to Dipole Response Function in 16 O with NN forces derived from EFT (N 3 LO) Convergence in the model space expansion L [fm 2 MeV -2 1-2 ] 4 3 2 1 Γ=1 MeV N max 18 16 14 12 1 8 hω=26 MeV L [fm 2 MeV -2 1-2 ] 4 3 2 1 Γ=1 MeV N max =18 hω=26 MeV hω=2 MeV -2 2 4 6 8 1 12 σ [MeV] Good convergence! -2 2 4 6 8 1 12 σ [MeV] Small HO dependence: use it as error bar 13

LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass Extension to Dipole Response Function in 16 O with NN forces derived from EFT (N 3 LO) Comparison to the experiment L [fm 2 MeV -2 1-2 ] 4 3 2 1 = 1 MeV LIT of data CCSD -2 2 4 6 8 1 12 σ [MeV] 14

LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass Extension to Dipole Response Function in 16 O with NN forces derived from EFT (N 3 LO) Comparison to the experiment 4 = 1 MeV LIT of data 5 16 O Ahrens et al. L [fm 2 MeV -2 1-2 ] 3 2 1 CCSD σ γ (ω)/4π 2 αω [mb/mev] 4 3 2 1 Ishkanov et al. CCSD -2 2 4 6 8 1 12 σ [MeV] 2 4 6 8 1 ω[mev] 15

LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass Extension to Dipole Response Function in 16 O with NN forces derived from EFT (N 3 LO) Comparison to the experiment 4 = 1 MeV LIT of data 5 16 O Ahrens et al. L [fm 2 MeV -2 1-2 ] 3 2 1 CCSD σ γ (ω)/4π 2 αω [mb/mev] 4 3 2 1 Ishkanov et al. CCSD -2 2 4 6 8 1 12 σ [MeV] 2 4 6 8 1 ω[mev] The GDR of 16 O is described from first principles for the first time! 15

LIT with Coupled Cluster Theory PRL 111, 12252 (213) New theoretical method aimed at extending ab-initio calculations towards medium mass Extension to Dipole Response Function in 16 O with NN forces derived from EFT (N 3 LO) Comparison to the experiment L [fm 2 MeV -2 1-2 ] 4 3 2 1 = 1 MeV LIT of data CCSD σ γ (ω)/4π 2 αω [mb/mev] 5 4 3 2 1 16 O Ahrens et al. Ishkanov et al. CCSD PRL 19 9252 (212) -2 2 4 6 8 1 12 σ [MeV] 2 4 6 8 1 ω[mev] The GDR of 16 O is described from first principles for the first time! 15

LIT with Coupled Cluster Theory Extension to Dipole Response Function in 22 O with NN forces derived from EFT (N 3 LO) L [fm 2 MeV -2 1-2 ] Work in Progress 6 22 5 O 4 3 2 1 = 1 MeV N max 16 14 12 1 8 hω=24 MeV -2 2 4 6 8 1 12 σ [MeV] Convergence pattern similar to 16 O 16

LIT with Coupled Cluster Theory L [fm 2 MeV -2 1-2 ] Extension to Dipole Response Function in 22 O with NN forces derived from EFT (N 3 LO) Work in Progress 6 22 5 O 4 3 2 1 = 1 MeV N max 16 14 12 1 8 hω=24 MeV -2 2 4 6 8 1 12 σ [MeV] L [fm 2 MeV -2 1-2 ] 5 4 3 2 Bremsstrahlung sum rule / NZ A 1 22 O 16 O 2 R PN R PN Preliminary Convergence pattern similar to 16 O -2 2 4 6 8 1 12 σ [MeV] 16

LIT with Coupled Cluster Theory L [fm 2 MeV -2 1-2 ] Extension to Dipole Response Function in 22 O with NN forces derived from EFT (N 3 LO) Work in Progress 6 22 5 O 4 3 2 1 = 1 MeV N max 16 14 12 1 8 hω=24 MeV -2 2 4 6 8 1 12 σ [MeV] L [fm 2 MeV -2 1-2 ] 5 4 3 2 Bremsstrahlung sum rule / NZ A 1 22 O 16 O 2 R PN R PN Preliminary Convergence pattern similar to 16 O -2 2 4 6 8 1 12 σ [MeV] The centroid of 22 O seems ~ 3 MeV lower than for 16 O 16

LIT with Coupled Cluster Theory with M.Miorelli Calcium isotopes with NN(N 3 LO) L [fm 2 MeV -2 ].1.5 4 Ca N max 8 1 12 14 16 18 ~ = 24 MeV Preliminary -2 2 4 6 8 1 12 σ[mev] Theoretical GDR peaks are located at higher energies with respect to experiment: nuclei are over bound with smaller radii and thus one needs higher energy to excited the GDR L(ω ) [fm 2 MeV -2 ].1.5 48 Ca N max 8 1 12 14 16 ~ = 24 MeV Preliminary -2 2 4 6 8 1 12 ω [MeV] 17

LIT with Coupled Cluster Theory with M.Miorelli Electric Dipole Polarizability D = 1 E 2 2 with NN(N 3 LO) Z 1 d! (!)! th! 2 calculated with the Lanczos algorithm α E [fm 3 ].1.1 4 He EIHH hω=2 MeV.9 hω=24 MeV.9.8.8 EXPERIMENT hω=26 MeV.7.7.6.6 6 8 1 12 14 16 18 N max.8.8 16 O.7.7.6.6 EXPERIMENT hω = 2MeV hω = 24MeV hω = 26MeV.5.5.4.4.3.3.2.2 6 8 1 12 14 16 18 N max α D [fm 3 ] 18

LIT with Coupled Cluster Theory with M.Miorelli Electric Dipole Polarizability D = 1 E 2 2 with NN(N 3 LO) Z 1 d! (!)! th! 2 calculated with the Lanczos algorithm α E [fm 3 ] Towards an ab-initio theory for 48 Ca 2.5 2.5 2.4 2.4 2.3 DFT prediction 2.3 48 Ca 2.2 2.2 2.1 2.1 2 2 1.9 1.9 1.8 1.8 1.7 1.7 1.6 1.6 1.5 1.5 1.4 hω = 2MeV 1.4 1.3 hω = 24MeV 1.3 hω = 26MeV 1.2 hω = 28MeV 1.2 Preliminary 1.1 1.1 1 1 6 8 1 12 14 16 N max 48 Ca E being measured at RCNP 19

LIT with Coupled Cluster Theory with M.Miorelli Electric Dipole Polarizability D = 1 E 2 2 Z 1 d! (!)! th! 2 Energy Density Functional Theory Phys. Rev. C 85, 4132 (212) very correlated to the neutron-skin radius Do we also see such a correlation with ab-initio methods? E 2

LIT with Coupled Cluster Theory with M.Miorelli Correlation of polarizability vs radius α E [fm 3 ] n2lo_opt_bare n2lo_opt38 (N=12).6 cdbonn_vlowk2.6 cdbonn_vlowk25 cdbonn_vlowk3 cdbonn_srg35 cdbonn_srg4.5 n3lo5em_srg2.5 n3lo5em_srg25 n3lo5em_srg3 n3lo5em_srg35 n3lo5em_bare.4 av18_vlowk25.4 av18_vlowk3 n3lo6em_srg25 n3lo6em_srg3 n3lo6em_srg35 16 O.3.3 Preliminary Exp +3NF.2.2 C AB =.961.1.1 2.1 2.2 2.3 2.4 2.5 2.6 2.7 r charge [fm] 21

LIT with Coupled Cluster Theory with M.Miorelli Correlation of polarizability vs radius α E [fm 3 ] n2lo_opt_bare n2lo_opt38 (N=12).6 cdbonn_vlowk2.6 cdbonn_vlowk25 cdbonn_vlowk3 cdbonn_srg35 cdbonn_srg4.5 n3lo5em_srg2.5 n3lo5em_srg25 n3lo5em_srg3 n3lo5em_srg35 n3lo5em_bare.4 av18_vlowk25.4 av18_vlowk3 n3lo6em_srg25 n3lo6em_srg3 n3lo6em_srg35 16 O Exp +3NF.3.3 Preliminary Optimized EFT interaction at N 2 LO from Ekström et al..2.2 C AB =.961.1.1 2.1 2.2 2.3 2.4 2.5 2.6 2.7 r charge [fm] 21

LIT with Coupled Cluster Theory with M.Miorelli Correlation of polarizability vs radius 2.6 2.6 α E [fm 3 ] 2.4 n2lo_opt_bare 2.4 n2lo_opt38 (N=12) 2.2 cdbonn_vlowk2 2.2 cdbonn_vlowk25 cdbonn_vlowk3 cdbonn_srg35 2 2 1.8 cdbonn_srg4 1.8 +3NF n3lo5em_srg2 +3NF 1.6 1.6 n3lo5em_srg25 n3lo5em_srg3 n3lo5em_srg35 1.4 1.4 1.2 n3lo6em_srg2 1.2 n3lo6em_srg25 n3lo6em_srg3 4 Ca Exp 1 1.8.8.6.6 Preliminary Optimized EFT interaction at N 2 LO from Ekström et al. C AB =.989.4.4.2.2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 r charge [fm] 22

LIT with Coupled Cluster Theory with M.Miorelli Correlation of polarizability vs skin-radius α E [fm 3 ] 48 Ca 2.6 2.6 2.4 DFT predictions 2.4 2.2 2.2 n2lo_opt38 (N=12) n2lo_opt_bare 2 cdbonn_vlowk3 2 cdbonn_vlowk25 cdbonn_vlowk2 +3NF 1.8 cdbonn_srg4 1.8 cdbonn_srg35 av18_vlowk2 1.6 1.6 av18_srg35 av18_srg3 1.4 n3lo5em_srg2 1.4 n3lo5em_srg25 n3lo6em_srg2 1.2 1.2 n3lo6em_srg25 n3lo6em_srg35 1 1.8.8 Preliminary Optimized EFT interaction at N 2 LO from Ekström et al..6.6 C AB =.974.4.4.13.14.15.16.17.18.19.2.21.22 r skin [fm] 48 Ca E being measured at RCNP 48 Ca parity violating electron scattering at JLab, CREX 23

Conclusions and Outlook Electromagnetic break up reactions are very rich observables to test our understanding of nuclear forces Extending these calculations to medium mass nuclei is possible and very exciting, with hopefully more applications/impact on future experiments. June 25th 214 Sonia Bacca 24

Conclusions and Outlook Electromagnetic break up reactions are very rich observables to test our understanding of nuclear forces Extending these calculations to medium mass nuclei is possible and very exciting, with hopefully more applications/impact on future experiments. Perspectives Dipole response function of neutron-rich Oxygen isotope Other multipole excitation (quadrupole or monopole) of medium mass nuclei need extension to two-body operators Add triples and three-nucleon forces June 25th 214 Sonia Bacca 24

Conclusions and Outlook Electromagnetic break up reactions are very rich observables to test our understanding of nuclear forces Extending these calculations to medium mass nuclei is possible and very exciting, with hopefully more applications/impact on future experiments. Perspectives Dipole response function of neutron-rich Oxygen isotope Other multipole excitation (quadrupole or monopole) of medium mass nuclei need extension to two-body operators Add triples and three-nucleon forces Thank you! June 25th 214 Sonia Bacca 24