1 / 19

Ab initio Coupled Cluster Calculations of Nuclei: Status and Outlook

and. Ab initio Coupled Cluster Calculations of Nuclei: Status and Outlook. UNEDF: J. Holt (UT), T. Papenbrock (UT/ORNL) CS support: H. A. Nam (ORNL)

long
Download Presentation

Ab initio Coupled Cluster Calculations of Nuclei: Status and Outlook

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. and Ab initio Coupled Cluster Calculations of Nuclei: Status and Outlook UNEDF: J. Holt (UT), T. Papenbrock (UT/ORNL) CS support: H. A. Nam (ORNL) External collaborators: D. J. Dean (ORNL), M. Hjorth-Jensen (Oslo), G. Baardsen (Oslo), G. Jansen (Oslo), O. Jensen (Bergen), T. Lesinski (Seattle) Gaute Hagen • Status and main accomplishments since last meeting • Status regarding high-performance computing • Roadmap for remainder of year 5 and outlook Annual UNEDF collaboration meeting East Lansing, MI, June 19-24 2011 Research partly funded by the US Department of Energy

  2. Personnel • UNEDF collaborators: G. Hagen, J. Holt, T. Papenbrock, T. Lesinski (UW) • CS support: Hai Ah Nam (ORNL) • Development of coupled-cluster method with relevance for UNEDF (by students supervised by M. Hjorth-Jensen and G. Hagen) • Øyvind Jensen (student at University of Bergen – Graduated July 2011; one-nucleon overlap functions and spectroscopic factors in J-coupled scheme) • Gustav Jansen (student at University of Oslo; closed shell ±2 nucleons) • Gustav Baardsen (student at University of Oslo; coupled-cluster theory for infinite nuclear matter)

  3. Publication Summary Published / accepted since July 2010 Quenching of spectroscopic factors for proton removal in oxygen isotopes, Ø. Jensen, G. Hagen, M. Hjorth-Jensen, B. Alex Brown, and A. Gade, Phys. Rev. Lett., in press (2011), arXiv:1104.1552 Ab initio coupled-cluster approach to nuclear structure with modern nucleon-nucleon interactions, G. Hagen, T. Papenbrock, D. J. Dean, and M. Hjorth-Jensen, Phys. Rev. C 82, 034330 (2010). Fermionic-Bosonic Couplings in a Weakly-Deformed Odd-Mass Nucleus, 93Nb, J. N. Orce, J. D. Holt, A. Linnemann, C. J. McKay, et al. Phys. Rev. C 82, 044317 (2010) Computation of spectroscopic factors with the coupled-cluster method, O. Jensen, G. Hagen, T. Papenbrock, D. J. Dean, and J. S. Vaagen, Phys. Rev. C 82, 014310 (2010). Closed-shell properties of 24O with ab initio coupled-cluster theory, Ø. Jensen, G. Hagen, M. Hjorth-Jensen, and J. S. Vaagen, Phys. Rev. C(R), 83, 021305 (2011). Towards open-shell nuclei with coupled-cluster theory, G. R. Jansen, M. Hjorth-Jensen, G. Hagen, and T. Papenbrock Phys. Rev. C 83, 054306 (2011) Submitted articles Ab initio computation of circular quantum dots, M. Pedersen Lohne, G. Hagen, M. Hjorth-Jensen, S. Kvaal, and F. Pederiva, Phys. Rev. B (2011) submitted, arXiv:1009.4833v1. Three-Body Forces and Shell Structure in Calcium Isotopes, J. D. Holt, T. Otsuka, A. Schwenk, and T. Suzuki arXiv:1009.5984, submitted to Phys. Rev. Lett. Chiral Three-Nucleon Forces and Bound Excited States in Neutron-Rich Oxygen Isotopes, J. D. Holt, A. Schwenk, T. Otsuka, and T. Suzuki,Submitted to Phys. Rev. C Role of many-body processes and three-nucleon forces in nuclear pairing, J. D. Holt, J. Menendez, and A. Schwenk, Submitted to Phys. Rev. C

  4. “Coupled-cluster approach to medium-mass and neutron-rich nuclei”, Nuclear Structure 2010 (NS2010), Clark-Kerr Campus, U. C. Berkeley, CA, August 8th – August 13th, 2010. (GH) “Coupled-cluster approach to medium-mass and neutron-rich nuclei”, Frontiers In NUclear STructure, Astrophysics and Reactions (FINUSTAR 3),Rhodes, Greece, August 23 to August 27, 2010. (GH) “Probing the driplines with ab‐initio Coupled‐Cluster theory”, ECT Workshop on Limits of existence of Light nuclei, ECT, Trento, Italy, October 25, 2010. (GH) “The nuclear many-body problem”, Second Uio-MSU-ORNL-UT School on Topics in Nuclear Physics: Many-Body Theory, Connections to Experiment and Nuclear Astrophysics (five lectures), January 3-7, 2011, Michigan State University. (GH) “Towards ab initio structure and reactions in neutronrich nuclei with coupled‐cluster theory”, ECT workshop on Effective theories and the nuclear many-body problem, ECT, Trento, Italy, March 7-11, 2011. (GH) “Towards ab initio structure and reactions in neutronrich nuclei with coupled‐cluster theory”, invited seminar talk at CEA, CEA, Saclay, May 13, 2011. (GH) “Towards a unified description of nuclear structure and reactions with coupled-cluster theory”, FUSTIPEN Topical Meeting «Structure and reactions at the drip lines», May 24-25, 2011, GANIL, Caen, France. (GH) “Ab initio calculations of weakly bound nuclei”, The first International Conference on Advances in Radioactive Isotope Science (ARIS - 2011), Leuven, Belgium, June 1, 2011. (GH) "Nuclear structure at the neutron drip-line", Gordon Research Conferences, Intersections Between Structure and Reactions: Pushing the Frontiers of Nuclear Science, Colby-Sawyer College, New London, June 14th 2011. (GH) “Ab initio Coupled Cluster Calculations of Nuclei: Status and Outlook”, annual UNEDF meeting, Henry Center, East Lansing, June 20th, 2011. “Model independent computation of atomic nuclei: status and perspectives,” “The INT at 20 years: The Future of Nuclear Physics and its Intersections,” Seattle, WA, July 1-2, 2010. (TP) “Ab initio nuclear structure calculations,” EMMI workshop “Strongly coupled systems,” GSI, Darmstadt, Germany, November 15-17, 2010. (TP) “Frontiers in nuclear structure theory,” Physics Colloquium, Technische Universität Darmstadt, Darmstadt, Germany, December 3, 2010. (TP) “Toward model-independent nuclear structure computations,” Physics seminar, Universität Frankfurt, Frankfurt, Germany, February 24, 2011. (TP) “Model-independent approaches to nuclear structure,” Annual NuSTAR meeting, GSI, Darmstadt, Germany, March 2-4, 2011. (TP) “Model-independent nuclear structure computations,” HISKP Colloquium, Universität Bonn, Bonn, Germany, May 12 2011 (TP) “Toward model-independent nuclear structure computations,” ECT* workshop on “Nuclear many-body open quantum systems,” Trento, Italy, June 7 2011. (TP) “Three-Nucleon Forces for Neutron-Rich, Medium-Mass Nuclei”, International Nuclear Physics Conference Vancouver, BC, Canada; July 4, 2010 (JH) “Three-Nucleon Forces for Medium Mass Nuclei Towards the Driplines“, Pan-American Advanced Studies Institute on Rare Isotopes Joao Pessoa, Brazil; August 1, 2010 (JH) “Three-Nucleon Forces and the Evolution of Nuclear Structure Towards the Driplines”, ECT* Workshop: “The Limits of Existence of Light Nuclei” Trento, Italy; October 27, 2010 (JH) “Three-Nucleon Forces and the Evolution of Nuclear Structure in Exotic Nuclei”, Nuclear Theory / EMMI Seminar; Technical University Darmstadt, Darmstadt, Germany; November 29, 2010.(JH) “Exploring Neutron-Rich Exotic Nuclei with Three-Nucleon Forces”, Nuclear Physics Seminar Instituto de Fisica Teorica Universidad Autonoma de Madrid, Madrid, Spain; February 2, 2011 (JH) “Three-Body Forces and Shell Structure in Exotic Nuclei”, ECT* Workshop: “Effective Theories and The Nuclear Many-Body Problem” Trento, Italy; March 7, 2011 (JH) “Three-Nucleon Forces and the Evolution of Nuclear Structure in Exotic Nuclei”, The Fifth LACM-EFES-JUSTIPEN Workshop Oak Ridge, TN; March 17, 2011. (JH) Talks

  5. Main accomplishments since last meeting • Science: • Long article on nuclear structure with NN interactions from chiral EFT • Quenching of spectroscopic factors for proton removal in neutron oxygen isotopes (Phys. Rev. Lett. 2011). • Computation of one-nucleon overlap functions for nuclear reactions and scattering (unpublished) • Computation of intrinsic densities with application to oxygen isotopes (unpublished) • Towards open shell nuclei with CCM extension to A±2: 18O and 26F • Coupled-cluster theory for nuclear matter. Implemented exact Pauli operator and HF + MBPT(2). • Constraining EDFs from ab-initio CC theory (See J. Holt’s talk) • Computing / algorithm developments: • Implementation of J-coupled code for computing one-nucleon overlaps and spectroscopic factors • Implementation of two-particle attached EOM-CCSD in J-coupled scheme.

  6. Saturation of N3LO (NN only) in medium mass nuclei Benchmarks in light nuclei: Coupled-cluster meets few-body benchmarks for 4He. Recent IT-NCSM and UMOA calculations of 16O agree with CCM. R. Roth et al, arXiv:1105.3173 (2011) Fujii et al, PRL 103, 182501(2009) G. Hagen, T. Papenbrock, D. J. Dean, M. Hjorth-Jensen, Phys. Rev. C 82, 034330 (2010).

  7. Overlap functions and spectroscopic factors Microscopic definition of Spectroscopic Factors (SF): Elastic scattering, capture and transfer of a nucleon on/to a target nucleus with mass A is given by the overlap function The norm of the overlap functions defines the spectroscopic factor: Asymptotic properties of the one-nucleon overlap function: Outside the range of the interaction the overlap function is proportional to a single-particle wave function.

  8. Quenching of SFs for proton removal in the oxygens • SF is a useful tool to study the role correlations toward the dripline. • We find a significant quenching of the SFs due to enhanced correlations coming from coupling to the scattering continuum. • Ø. Jensen, G. Hagen and M. Hjorth-Jensen, Phys. Rev. C(R) 83, 021305 (2011) • Ø. Jensen, G. Hagen, Hjorth-Jensen, Brown, Gade, in press Phys. Rev. Lett. (2011) arXiv:1104.1552 • For threshold effects in SFs see N. Michel et al Phys. Rev. C(R) 75, 031301 (2007)

  9. Asymptotic normalization coefficients Overlap function for a bound A + 1 nucleus: We use an SRG evolved interaction with cutoff 2.66 fm-1. The CCSD ground state energy for 16O in N = 11 major shells is 140.52 MeV. One neutron overlap functions for the bound J = 1/2+ state in 17O with the ground state of 16O. For 17O we get E(1/2+) = -3.83 MeV.

  10. Elastic scattering of neutrons on 16O with CCM(Preliminary) Overlap functions provides a simple and intuitive picture of reactions Left figure: One neutron overlap functions for various J = 1=2+ scattering states in 17O using SRG evolved interaction with cutoff = 2.66fm-1 Right figure: By matching the known asymptotic form of the overlap functions to scattering solutions we can extract low-energy elastic scattering phase shifts.

  11. Intrinsic densities from microscopic CCM Computing intrinsic densities The one-body density matrix for the A and A-1 nuclei are computed in the coupled-cluster formalism and the local laboratory density is CC wave function factorizes and the CoM wave function is a Gaussian with fixed width (G. Hagen et al Phys. Rev. Lett. 103, 062503 (2009)) We can get the internal density by a deconvolution of the local laboratory density (B. G. Giraud, Phys. Rev. C 77, 014311 (2008))

  12. Matter and charge radii in the 21,22,23,24Oisotopes Binding energies of 21,22,23,24O with respect to the ground state of 21O. Point matter and charge radii for computed from intrinsic densities For the nucleon-nucleon (NN) interaction we used the SRG evolved N3LO nucleon-nucleon interaction with momentum cutoffs = 3.6, 3.8, 4.0fm-1.

  13. Matter and charge radii in the 21,22,23,24Oisotopes Binding energies of 21,22,23,24O with respect to the ground state of 21O. Point matter and charge radii for computed from intrinsic densities For the nucleon-nucleon (NN) interaction we used the SRG evolved N3LO nucleon-nucleon interaction with momentum cutoffs = 3.6, 3.8, 4.0fm-1.

  14. Going beyond closed-shell nuclei. Low-lying states in 18O and 26F (Preliminary) Two-particle attached coupled-cluster works very well for low-lying states in open-shell nuclei like 18O and 26F. Our results for 26F seem to suggest a more compressed spectrum as compared to USDA/USDB calculations. G. Jansen, M. Hjorth-Jensen, G. Hagen, T. Papenbrock, Phys. Rev. C 83, 054306 (2011).

  15. Low-lying states in neutron rich Potassium isotopes Proton separation energies 40,48,52,54,60Ca. Low lying states in Potassium isotopes calculated using PA/PR-EOMCCSD with “bare" chiral interactions. Model space consists of 15 major harmonic oscillator shells with fixed oscillator frequency hw = 30MeV.

  16. Coupled-cluster theory for nuclear matter(Preliminary) Implemented exact Pauli operator in relative-center of mass coordinates Implemented Hartree-Fock and MBPT(2) for nuclear matter in relative and Center-of-mass coordinates. Code has been validated and verified using Argonne-V18 and vlow-k. Particle-particle and hole-hole ladders are in progress. Particle-hole channels will be implemented using exact and angle average Pauli operator

  17. Progress for years 4-5 Project is reasonably well on track

  18. Computational status • Data organization • Interaction spread across processors • cluster amplitudes stored locally • oxygen-16 in 20 shells ~ 1500 processor hours per model space • oxygen-28 with continuum: total cost of about 100,000 CPU hours on Jaguar. • Challenge • Number of j-coupled matrix elements of interactions fluctuate strongly for given sets of quantum numbers • load balancing non-trivial: Calculation of estimated computational cost  distribution of data • scaling (load-balancing) up to few thousand processors • Working on optimizing code with OpenMP.

  19. Future plans (remainder of year 5 and outlook) • Interface with DFT, and toward ab initio reactions • Compute phaseshifts for elastic scattering of neutrons on selected closed shell nuclei using one-nucleon overlap functions. • Constraining EDFs from ab-initio coupled-cluster theory • Develop spherical EOM-CCSD code for general spin and compute dipole response function using Lorentz integral transform. • Role of three-nucleon forces in medium-mass nuclei • Outlook / SciDAC-3 ? • Ab initio computation of 48Ca with 3NF as possible highlight • Validation and verification of nuclear forces from Chiral EFT with uncertainty quantification. • Exploring Magic Numbers at the extreme using nuclear coupled-cluster • One- and two-nucleon overlap functions from coupled-cluster theory • Time dependent coupled-cluster approach to reactions • Nuclear matter (PhD student Gustav Baardsen from UiO)

More Related