ab initio no-core shell model/ resonating group method › conferences › elba › 2014_talks ›...

58
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 Canada’s national laboratory for particle and nuclear physics Laboratoire national canadien pour la recherche en physique nucléaire et en physique des particules Accelerating Science for Canada Un accélérateur de la démarche scientifique canadienne Carolina Romero-Redondo Petr Navrátil, Sofia Quaglioni and Guillaume Hupin Ab initio no-core shell model/ resonating group method

Upload: others

Post on 07-Feb-2021

1 views

Category:

Documents


0 download

TRANSCRIPT

  • 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

    Canada’s national laboratory for particle and nuclear physics

    Laboratoire national canadien pour la recherche en physique nucléaire

    et en physique des particules

    Accelerating Science for Canada

    Un accélérateur de la démarche scientifique canadienne

    Carolina Romero-Redondo

    Petr Navrátil, Sofia Quaglioni and Guillaume Hupin

    Ab initio no-core shell model/resonating group method

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Introduction: NCSM/RGM

    Binary clusters

    Three-body cluster states

    6He: 4He+n+n

    NCSMC

    Summary and outlook

    Outline

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Introduction

    * Assumes nucleons as the effective degrees of freedom

    * Uses realistic interactions

    The goal is to achieve a predictive theory for light nuclear systems to study:

    - Exotic nuclei

    - Reactions important in nuclear astrophysics

    - Reactions important for energy production projects

    Ab initio in nuclear physicsnucleus

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM

    No-core shell model (NCSM)Is an ab initio method capable of studying light bound nuclei from an accurate Hamiltonian.

    Traditionally, has not been able to deal with continuum states.

    Resonances

    ReactionsHalo nuclei wf

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM

    No-core shell model (NCSM)Is an ab initio method capable of studying light bound nuclei from an accurate Hamiltonian.

    Traditionally, has not been able to deal with continuum states.

    Resonating group method (RGM)Microscopic cluster approach.

    Permits studying the scattering of clusters

    Non-realistic Hamiltonian

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    No-core shell model (NCSM)Is an ab initio method capable of studying light bound nuclei from an accurate Hamiltonian.

    Is not able to deal with continuum states and therefore is not applicable to

    reactions.

    Resonating group method (RGM)Microscopic cluster approach.

    Permits studying the scattering of clusters

    Non-realistic Hamiltonian

    NCSM/RGM

    Combines NCSM and RGM to obtain anab initio formalism which uses an accurate

    nuclear Hamiltonian and is capable of studying both structure and scattering

    problems in light nuclear systems

    NCSM/RGM

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-2B

    S. Quaglioni and P. Navrátil- PRL 101, 092501 (2008)- PRC 79, 044606 (2009)

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    n-4He

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-2B

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    R-matrix

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    R-matrix

    R matrix:

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    R-matrix

    R matrix:

    Scattering Matrix

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    A-nucleon system

    S. Quaglioni and P. Navratil

    - PRL 101, 092501 (2008)- PRC 79, 044606 (2009)

    NCSM/RGM

    Formalism presented for the first time. Nucleon-nucleous scattering: n+4He, n+3Hen+10Be, p+3He and p+4He

    Applications

    Radiative Capture: 7Be(p,γ)8B S-factor -PLB 704 (2011) 379 In collaboration with R. Roth

    -PRC 83, 044609 (2011) Deuteron projectile: scattering (d+4He), 6Li states-PRL 108, 042503 (2012) Fusion reactions: 3He(d,p)4He and 3H(d,n)4He

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    C. Romero-Redondo

    S. Quaglioni, P. Navrátil, G. Hupin

    Extension to three-body cluster

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    Why?

    Extension to three-body cluster

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    Extension to three-body cluster

    Bound and resonant states:2n Halo nuclei

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    Extension to three-body cluster

    9Li

    n

    n4He n

    n

    Bound and resonant states:2n Halo nuclei

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    Extension to three-body cluster

    9Li

    n

    n4He n

    n

    3-body continuum states:Transfer reactions

    Bound and resonant states:2n Halo nuclei

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    Extension to three-body cluster

    +3H 3H 4He n

    n

    n

    n

    9Li

    n

    n4He n

    n

    3-body continuum states:Transfer reactions

    Bound and resonant states:2n Halo nuclei

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    Basis

    36

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    Basis

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    NCSM wave functions

    Basis

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    Schrödinger equation

    Norm kernelHamiltonian Kernel

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    Schrödinger equation

    Norm kernelHamiltonian Kernel Relative motion wavefunction

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Orthogonalization

    NCSM/RGM-3B

    Schrödinger equation

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSM/RGM-3B

    Hyperspherical coordinates:

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Hyperspherical coordinates:

    NCSM/RGM-3B

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Hyperspherical coordinates:

    NCSM/RGM-3B

    After changing to hyperspherical coordinates and integrating in α,α':

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Hyperspherical coordinates:

    After changing to hyperspherical coordinates and integrating in α,α':

    Coupled-channel microscopic R-matrix method on a Lagrange mesh

    NCSM/RGM-3B

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    R-matrix

    External region: known asymptotic behaviour (ρ > a)

    Internal region: expansion on a basis (ρ < a)

    * Bound state:

    * Continuum state:

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    4He+n+n

    4He n

    n6He: 2 neutron halo (4He-n-n)

    42

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    4He+n+n

    4He n

    n6He: 2 neutron halo (4He-n-n)

    Input:

    - Accurate soft NN potential(SRG-evolved chiral N3LO with λ=1.5fm-1)

    - NSCM 4He wave function

    - Fully antisymmetric channel states

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Norm Kernel

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Norm Kernel

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Hamiltonian Kernel

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    HO model space Eg.s. (4He) [MeV](NCSM)

    Eg.s. (6He) [MeV](NCSM)

    Eg.s. (6He) [MeV](NCSM/RGM)

    Nmax = 8 -28.173 -28.952 -28.616Nmax = 10 -28.215 -29.452 -28.696Nmax = 12 -28.224 -29.658 -28.697

    4He+n+n

    6He ground state, NCSMNCSM/RGM results- n+n+4He, Nmax= 11, ħΩ = 14 MeV

    - SRG-NN chiral with λ=1.5 fm-1

    Comparison with NCSM:

    - ~ 1MeV difference in Egs due to excitations

    of the 4He core only included in the NCSM

    calculation.

    -Contrary to NCSM, NCSM/RGM n+n+4He

    w.f. Has the appropriate asymptotic behaviour

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    4He+n+n

    6He ground state wave function

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    4He+n+n

    6He ground state wave functionProbability distribution

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    4He+n+n

    6He ground state wave functionProbability distribution

    Small

    Small

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    4He+n+n

    6He ground state wave functionProbability distribution

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    4He+n+n

    6He ground state wave functionProbability distribution

    Di-neutron

    Cigar

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    4He+n+n continuum

    Continuum

    Calculation of phase shifts for different channels and study the stability of the results

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    4He+n+n continuumDependence in Nmax

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Dependence in Kmax

    4He+n+n continuum

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    4He+n+n continuum

    SpectrumComparison with experiment

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    In progress

    No-core shell model with continuum (NCSMC) for three-cluster basis

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    NCSMC

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    +

    For a 3-cluster system like 6He:

    NCSMC

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Improvements in 6He calculationin order to obtain a complete ab initio calculation

    Introduce 4He core excitations by coupling the 4He-n-n basis to 6He NCSM eigenstates (NCSMC)

    Run calculations with 3N force

    4He+n+n

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    Transfer reactions, i.e, 3H(3H,2n)4He

    Derive and calculate couplings between two and three body clusters

    Other systems: 5H, 11Li, 12Be

    Outlook

    Approach is very useful for studying different types of nuclear systems

  • 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

    Canada’s national laboratory for particle and nuclear physics

    Laboratoire national canadien pour la recherche en physique nucléaire

    et en physique des particules

    Thank you!

    TRIUMF: Alberta | British Columbia | Calgary Carleton | Guelph | Manitoba | McMaster Montréal | Northern British Columbia | Queen’s Regina | Saint Mary’s | Simon Fraser | Toronto Victoria | Winnipeg | York

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    7He: 6He+n

  • C. Romero Redondo, Elba XIII Workshop. June 25th, 2014

    R-matrix

    Slide 1Slide 2Slide 3Slide 4Slide 5Slide 6Slide 7Slide 8Slide 9Slide 10Slide 11Slide 12Slide 13Slide 14Slide 15Slide 16Slide 17Slide 18Slide 19Slide 20Slide 21Slide 22Slide 23Slide 24Slide 25Slide 26Slide 27Slide 28Slide 29Slide 30Slide 31Slide 32Slide 33Slide 34Slide 35Slide 36Slide 37Slide 38Slide 39Slide 40Slide 41Slide 42Slide 43Slide 44Slide 45Slide 46Slide 47Slide 48Slide 49Slide 50Slide 51Slide 52Slide 53Slide 54Slide 55Slide 56Slide 57Slide 58