Showing 1–50 of 65 results for author: Papakonstantinou, P

Nuclear mass table in deformed relativistic Hartree-Bogoliubov theory in continuum, II: Even-$Z$ nuclei

Authors: DRHBc Mass Table Collaboration, Peng Guo, Xiaojie Cao, Kangmin Chen, Zhihui Chen, Myung-Ki Cheoun, Yong-Beom Choi, Pak Chung Lam, Wenmin Deng, Jianmin Dong, Pengxiang Du, Xiaokai Du, Kangda Duan, Xiaohua Fan, Wei Gao, Lisheng Geng, Eunja Ha, Xiao-Tao He, Jinniu Hu, Jingke Huang, Kun Huang, Yanan Huang, Zidan Huang, Kim Da Hyung, Hoi Yat Chan , et al. (58 additional authors not shown)

Abstract: The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-ne… ▽ More The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-neutron separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. A total of 4829 even-$Z$ nuclei are predicted to be bound, with an rms deviation of 1.477 MeV from the 1244 mass data. Good agreement with the available experimental odd-even mass differences, $αあるふぁ$ decay energies, and charge radii is also achieved. The description accuracy for nuclear masses and nucleon separation energies as well as the prediction for drip lines is compared with the results obtained from other relativistic and nonrelativistic density functional. The comparison shows that the DRHBc theory with PC-PK1 provides an excellent microscopic description for the masses of even-$Z$ nuclei. The systematics of the nucleon separation energies, odd-even mass differences, pairing energies, two-nucleon gaps, $αあるふぁ$ decay energies, rms radii, quadrupole deformations, potential energy curves, neutron density distributions, and neutron mean-field potentials are discussed. △ Less

Submitted 10 June, 2024; v1 submitted 5 February, 2024; originally announced February 2024.

Comments: 394 pages, 17 figures, 2 tables, published in Atomic Data and Nuclear Data Tables, data file in the TXT form is available for download under "Ancillary files"

Journal ref: Peng Guo, et. al. (DRHBc Mass Table Collaboration), Atomic Data and Nuclear Data Tables 158 (2024) 101661

Landscape of nuclear deformation softness with spherical quasi-particle random phase approximation

Authors: Le-Anh Nguyen, Minh-Loc Bui, Panagiota Papakonstantinou, Naftali Auerbach

Abstract: We investigate the stability and softness of nuclei against quadrupole, octupole, and hexadecapole deformation. By applying the spherical Skyrme-force Hartree-Fock Bardeen-Cooper-Schrieffer quasi-particle random phase approximation, we diagnose ground-state deformation when imaginary solutions are obtained, i.e., the spherical ground state {\em collapses}. We also calculate the multipole polarizab… ▽ More We investigate the stability and softness of nuclei against quadrupole, octupole, and hexadecapole deformation. By applying the spherical Skyrme-force Hartree-Fock Bardeen-Cooper-Schrieffer quasi-particle random phase approximation, we diagnose ground-state deformation when imaginary solutions are obtained, i.e., the spherical ground state {\em collapses}. We also calculate the multipole polarizability in spherical nuclei with no collapse, as a measure of softness. This numerically light and theoretically sound method is found able to capture deformation patterns across the nuclide chart. The connection between the intrinsic shape of nuclei and the dynamics of their low-lying collective states is established and the role of shell structure is discussed. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Comments: 8 pages, 8 figures, accepted for publication in Physical Review C

Variational and parquet-diagram calculations for neutron matter. V. Triplet pairing

Authors: E. Krotscheck\', P. Papakonstantinou, J. Wang

Abstract: We apply a large-scale summation of Feynman diagrams, including the class of parquet-diagrams {\em plus} important contributions outside the parquet class, for calculating effective pairing interactions and subsequently the superfluid gap in P-wave pairing in neutron matter. We employ realistic nucleon-nucleon interactions of the $v_8$ type and perform calculations up to a Fermi momentum of… ▽ More We apply a large-scale summation of Feynman diagrams, including the class of parquet-diagrams {\em plus} important contributions outside the parquet class, for calculating effective pairing interactions and subsequently the superfluid gap in P-wave pairing in neutron matter. We employ realistic nucleon-nucleon interactions of the $v_8$ type and perform calculations up to a Fermi momentum of $1.8\,$fm$^{-1}$. We find that many-body correlations lead to a strong reduction of the spin-orbit interaction, and, therefore, to an almost complete suppression of the $^3$P$_2$ and $^3$P$_2$-$^3$F$_2$ gaps. We also find pairing in $^3$P$_0$ states; the strength of the pairing gap depends sensitively on the potential model employed. Our results for triplet pairing are relevant for assessing superfluidity in neutron star interiors, whose presence can affect the cooling of neutron stars. △ Less

Submitted 12 November, 2023; v1 submitted 1 August, 2023; originally announced August 2023.

Missed prediction of the neutron halo in $^{37}$Mg

Authors: K. Y. Zhang, S. Q. Yang, J. L. An, S. S. Zhang, P. Papakonstantinou, M. -H. Mun, Y. Kim, H. Yan

Abstract: Halo phenomena have long been an important frontier in both experimental and theoretical nuclear physics. $^{37}$Mg was identified as a halo nucleus in 2014 and remains the heaviest nuclear halo system to date. While the halo phenomenon in $^{37}$Mg was not predicted before the discovery, its description has been still challenging afterwards. In this Letter, we report a microscopic and self-consis… ▽ More Halo phenomena have long been an important frontier in both experimental and theoretical nuclear physics. $^{37}$Mg was identified as a halo nucleus in 2014 and remains the heaviest nuclear halo system to date. While the halo phenomenon in $^{37}$Mg was not predicted before the discovery, its description has been still challenging afterwards. In this Letter, we report a microscopic and self-consistent description of the neutron halo in $^{37}$Mg using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) that was developed in 2010. The experimental neutron separation energies and empirical matter radii of neutron-rich magnesium isotopes as well as the deformed $p$-wave halo characteristics of $^{37}$Mg are well reproduced without any free parameters. In particular, the orbital occupied by the halo neutron in $^{37}$Mg, exhibiting $p$-wave components comparable to those suggested in experiments, remains consistent across various employed density functionals including PC-F1, PC-PK1, NL3*, and PK1. The DRHBc theory investigated only even-even magnesium isotopes in previous works and for that reason missed predicting $^{37}$Mg as a halo nucleus before 2014. Although the core and the halo of $^{37}$Mg are both prolate, higher-order shape decoupling on the hexadecapole and hexacontatetrapole levels is predicted. △ Less

Submitted 9 August, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: updated references and figure 3; matches the published version; 8 pages, 4 figures, 1 table

Journal ref: Phys. Lett. B 844 (2023) 138112

Triplet paring in neutron matter

Authors: Eckhard Krotscheck, Panagiota Papakonstantinou, Jiawei Wang

Abstract: The presence of superfluidity in neutron star interiors can affect the cooling of neutron stars in intricate ways, enhancing certain mechanisms and suppressing others. Model calculations employing realistic nuclear potentials in Bardeen-Cooper-Schrieffer theory generally suggest the development of a $^3$P$_2$-$^3$F$_2$ pairing gap, and therefore the presence of superfluidity in dense neutron-star… ▽ More The presence of superfluidity in neutron star interiors can affect the cooling of neutron stars in intricate ways, enhancing certain mechanisms and suppressing others. Model calculations employing realistic nuclear potentials in Bardeen-Cooper-Schrieffer theory generally suggest the development of a $^3$P$_2$-$^3$F$_2$ pairing gap, and therefore the presence of superfluidity in dense neutron-star matter.Improved models that go beyond conventional mean-field calculations by including polarization effects suggest a suppression of the triplet gap. We have evaluated the pairing interaction by summing the ``parquet'' Feynman diagrams which include both ladder and ring diagrams systematically, {\em plus\} a set of important non-parquet diagrams, making this the most comprehensive diagram-based approach presently available. Our results suggest a radical suppression of the $^3$P$_2$-$^3$F$_2$ triplet pairing gap and an enhancement of $^3$P$_0$ pairing. △ Less

Submitted 31 July, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Collapse of the $N=28$ shell closure in the newly discovered $^{39}$Na and the development of deformed halos towards the neutron dripline

Authors: K. Y. Zhang, P. Papakonstantinou, M. -H. Mun, Y. Kim, H. Yan, X. -X. Sun

Abstract: Halos and changes of nuclear magicities have been extensively investigated in exotic nuclei during past decades. The newly discovered $^{39}$Na with the neutron number $N=28$ provides a new platform to explore such novel phenomena near the neutron dripline of the sodium isotopic chain. We study the shell property and the possible halo structure in $^{39}$Na within the deformed relativistic Hartree… ▽ More Halos and changes of nuclear magicities have been extensively investigated in exotic nuclei during past decades. The newly discovered $^{39}$Na with the neutron number $N=28$ provides a new platform to explore such novel phenomena near the neutron dripline of the sodium isotopic chain. We study the shell property and the possible halo structure in $^{39}$Na within the deformed relativistic Hartree-Bogoliubov theory in continuum. It is found that the lowering of $2p$ orbitals in the spherical limit results in the collapse of the $N=28$ shell closure in $^{39}$Na, and a well deformed ground state is established. The pairing correlations and the mixing of $pf$ components driven by deformation lead to the occupation of weakly bound or continuum $p$-wave neutron orbitals. An oblate halo is therefore formed around the prolate core in $^{39,41}$Na, making $^{39}$Na a single nucleus with the coexistence of several exotic structures, including the quenched $N=28$ shell closure, Borromean structure, deformed halo, and shape decoupling. The microscopic mechanisms behind the shape decoupling phenomenon and the development of halos towards dripline are revealed. △ Less

Submitted 4 May, 2023; originally announced May 2023.

Comments: 8 pages, 5 figures

Journal ref: Phys. Rev. C 107, L041303 (2023)

Origin of octupole deformation softness in atomic nuclei

Authors: Minh-Loc Bui, Le-Anh Nguyen, Panagiota Papakonstantinou, Naftali Auerbach

Abstract: Recent high-energy heavy ion collision experiments have revealed that some atomic nuclei exhibit unusual softness and significant shape fluctuations. In this work, we use the fully self-consistent mean-field theory to identify all even-even nuclei that are unstable or soft against octupole deformation. All exceptional cases of enhanced octupole transition strengths in stable even-even nuclei throu… ▽ More Recent high-energy heavy ion collision experiments have revealed that some atomic nuclei exhibit unusual softness and significant shape fluctuations. In this work, we use the fully self-consistent mean-field theory to identify all even-even nuclei that are unstable or soft against octupole deformation. All exceptional cases of enhanced octupole transition strengths in stable even-even nuclei throughout the nuclide chart are resolved and the origin is found in basic shell structure. The presence of atomic nuclei exhibiting significant softness to quadrupole-octupole deformation is suggested. These results represent a significant advance in our understanding of the underlying mechanisms of nuclear octupole deformation and have implications for further experimental and theoretical studies. △ Less

Submitted 25 July, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: 7 pages, 2 figures, accepted for publication in Physical Review C

Dense Nuclear Matter Equation of State from Heavy-Ion Collisions

Authors: Agnieszka Sorensen, Kshitij Agarwal, Kyle W. Brown, Zbigniew Chajęcki, Paweł Danielewicz, Christian Drischler, Stefano Gandolfi, Jeremy W. Holt, Matthias Kaminski, Che-Ming Ko, Rohit Kumar, Bao-An Li, William G. Lynch, Alan B. McIntosh, William G. Newton, Scott Pratt, Oleh Savchuk, Maria Stefaniak, Ingo Tews, ManYee Betty Tsang, Ramona Vogt, Hermann Wolter, Hanna Zbroszczyk, Navid Abbasi, Jörg Aichelin , et al. (111 additional authors not shown)

Abstract: The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of mu… ▽ More The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. This white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS. △ Less

Submitted 25 January, 2024; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: White paper prepared for the 2023 Long Range Plan. v3: Updated version as published in Progress in Particle and Nuclear Physics. Note: the published version does not include the executive summary; in the updated arXiv version, the executive summary is included as an appendix. v4: Corrected list of authors

Report number: INT-PUB-23-001, LA-UR-23-20514, LLNL-TR-844629

Journal ref: Prog.Part.Nucl.Phys. 134 (2024) 104080

Bayesian inference of neutron-star observables based on effective nuclear interactions

Authors: Jia Zhou, Jun Xu, Panagiota Papakonstantinou

Abstract: Based on the Skyrme-Hartree-Fock model (SHF) as well as its extension (the Korea-IBS-Daegu-SKKU (KIDS) model) and the relativistic mean-field (RMF) model, we have studied the constraints on the parameters of the nuclear matter equation of state (EOS) from adopted astrophysical observables using a Bayesian approach. While the masses and radii of neutron stars generally favors a stiff isoscalar EOS… ▽ More Based on the Skyrme-Hartree-Fock model (SHF) as well as its extension (the Korea-IBS-Daegu-SKKU (KIDS) model) and the relativistic mean-field (RMF) model, we have studied the constraints on the parameters of the nuclear matter equation of state (EOS) from adopted astrophysical observables using a Bayesian approach. While the masses and radii of neutron stars generally favors a stiff isoscalar EOS and a moderately soft nuclear symmetry energy, model dependence on the constraints is observed and mostly originates from the incorporation of higher-order EOS parameters and difference between relativistic and non-relativistic models. At twice saturation density, the value of the symmetry energy is constrained to be $48^{+15}_{-11}$ MeV in the standard SHF model, $48^{+8}_{-15}$ MeV in the KIDS model, and $48^{+5}_{-6}$ MeV in the RMF model, around their maximum {\it a posteriori} values within $68\%$ confidence intervals. Our study helps to obtain a robust constraint on nuclear matter EOS, and meanwhile, to understand the model dependence of the results. △ Less

Submitted 5 May, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: 16 pages, 8 figures

Journal ref: Physical Review C 107, 055803 (2023)

Comparative analysis of formalisms and performances of three different beyond mean-field approaches

Authors: František Knapp, Panagiota Papakonstantinou, Petr Veselý, Giovanni De Gregorio, Jakub Herko, Nicola Lo Iudice

Abstract: We investigate the differences and analogies between the equation of motion phonon method (EMPM) and second Tamm-Dancoff and random-phase approximations (STDA and SRPA) paying special attention to the problem of spurious center-of-mass (c.m.) admixtures. In order to compare them on an equal footing, we perform self-consistent calculations of the multipole strength distributions in selected doubly… ▽ More We investigate the differences and analogies between the equation of motion phonon method (EMPM) and second Tamm-Dancoff and random-phase approximations (STDA and SRPA) paying special attention to the problem of spurious center-of-mass (c.m.) admixtures. In order to compare them on an equal footing, we perform self-consistent calculations of the multipole strength distributions in selected doubly magic nuclei within a space including up to two-particle-two-hole (2p-2h) basis states using the UCOM two-body intrinsic Hamiltonian and we explore the tools each approach supplies for removing the spurious c.m. admixtures. We find that the EMPM and STDA yield exactly the same results when the same intrinsic Hamiltonian is used and the coupling of the Hartree-Fock state with the 2p-2h space is neglected, but, unlike STDA and SRPA, the EMPM offers the possibility to completely remove c.m. admixtures. △ Less

Submitted 6 November, 2022; originally announced November 2022.

Comments: 13 pages, 17 figures

Nuclear symmetry energy and the PREX-CREX neutron skin puzzle within the KIDS framework

Authors: Panagiota Papakonstantinou

Abstract: I briefly review the KIDS theoretical framework for the nuclear equation of state (EoS) and energy density functional (EDF), I discuss recent results for the curvature parameter of the symmetry energy, and I address the PREX-CREX puzzle. I show that it is possible to obtain EDF models which can reproduce both PREX-II and CREX results each within its respective error bars. Such EDFs correspond to E… ▽ More I briefly review the KIDS theoretical framework for the nuclear equation of state (EoS) and energy density functional (EDF), I discuss recent results for the curvature parameter of the symmetry energy, and I address the PREX-CREX puzzle. I show that it is possible to obtain EDF models which can reproduce both PREX-II and CREX results each within its respective error bars. Such EDFs correspond to EoSs which soften towards low densities, as could be attributed to clusterization. Before such a scenario is considered viable, the dipole polarizability should also be examined. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: 10 pages incl. 1 figure; pilot study; prepared for the Proceedings of the 39th International Workshop on Nuclear Theory, Rila Mountains, Bulgaria, July 3-9, 2022

Deformed relativistic Hartree-Bogoliubov theory in continuum with a point-coupling functional. II. Examples of odd Nd isotopes

Authors: DRHBc Mass Table Collaboration, Cong Pan, Myung-Ki Cheoun, Yong-Beom Choi, Jianmin Dong, Xiaokai Du, Xiao-Hua Fan, Wei Gao, Lisheng Geng, Eunja Ha, Xiao-Tao He, Jinke Huang, Kun Huang, Seonghyun Kim, Youngman Kim, Chang-Hwan Lee, Jenny Lee, Zhipan Li, Zhi-Rui Liu, Yiming Ma, Jie Meng, Myeong-Hwan Mun, Zhongming Niu, Panagiota Papakonstantinou, Xinle Shang , et al. (16 additional authors not shown)

Abstract: The aim of this work is to extend the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) based on the point-coupling density functionals to odd-$A$ and odd-odd nuclei and examine its applicability by taking odd-$A$ Nd isotopes as examples. In the DRHBc theory, the densities and potentials with axial deformation are expanded in terms of Legendre polynomials, and the relativistic H… ▽ More The aim of this work is to extend the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) based on the point-coupling density functionals to odd-$A$ and odd-odd nuclei and examine its applicability by taking odd-$A$ Nd isotopes as examples. In the DRHBc theory, the densities and potentials with axial deformation are expanded in terms of Legendre polynomials, and the relativistic Hartree-Bogoliubov equations are solved in a Dirac Woods-Saxon basis to include the continuum effects. For an odd-$A$ or odd-odd nucleus, the blocking effect of unpaired nucleon(s) is taken into account with the equal filling approximation. To determine its ground state, an automatic blocking procedure is adopted, in which the orbital with the lowest quasiparticle energy is blocked during the iteration. This procedure is justified by comparing with the results from the orbital-fixed blocking calculations, in which the blocked orbital near the Fermi surface is fixed during the iteration. The ground states for both light and heavy nuclei can be provided by the automatic blocking procedure as the orbital-fixed blocking procedure, but with considerably reduced computational cost. The numerical details for even-even nuclei are found to be valid for odd-$A$ and odd-odd nuclei as well. Taking Nd isotopes including both even-even and odd-$A$ ones as examples, the calculated ground-state properties with PC-PK1 are in good agreement with the available experimental data. This work paves the way to construct the DRHBc mass table including all even-even, odd-$A$ and odd-odd nuclei in the nuclear chart. △ Less

Submitted 24 August, 2023; v1 submitted 3 May, 2022; originally announced May 2022.

Comments: 39 pages, 13 figures, and 1 table

Journal ref: Phys. Rev. C 106, 014316 (2022)

Damping of the isovector giant dipole resonance in $^{40,48}$Ca

Authors: J. Carter, L. M. Donaldson, H. Fujita, Y. Fujita, M. Jingo, C. O. Kureba, M. B. Latif, E. Litvinova, F. Nemulodi, P. von Neumann-Cosel, R. Neveling, P. Papakonstantinou, P. Papka, L. Pellegri, V. Yu. Ponomarev, A. Richter, R. Roth, E. Sideras-Haddad, F. D. Smit, J. A. Swartz, A. Tamii, R. Trippel, I. T. Usman, H. Wibowo

Abstract: The fine structure of the IsoVector Giant Dipole Resonance (IVGDR) in the doubly-magic nuclei $^{40,48}$Ca observed in inelastic proton scattering experiments under $0^\circ$ is used to investigate the role of different mechanisms contributing to the IVGDR decay width. Characteristic energy scales are extracted from the fine structure by means of wavelet analysis. The experimental scales are compa… ▽ More The fine structure of the IsoVector Giant Dipole Resonance (IVGDR) in the doubly-magic nuclei $^{40,48}$Ca observed in inelastic proton scattering experiments under $0^\circ$ is used to investigate the role of different mechanisms contributing to the IVGDR decay width. Characteristic energy scales are extracted from the fine structure by means of wavelet analysis. The experimental scales are compared to different theoretical approaches allowing for the inclusion of complex configurations beyond the mean-field level. Calculations are performed in the framework of RPA and beyond-RPA in a relativistic approach based on an effective meson-exchange interaction, with the UCOM effective interaction and, for the first time, with realistic two- plus three-nucleon interactions from chiral effective field theory employing the in-medium similarity renormalization group. All models highlight the role of Landau fragmentation for the damping of the IVGDR, while the differences in the coupling strength between one particle-one hole (1p-1h) and two particle-two hole (2p-2h) correlated (relativistic) and non-correlated (non-relativistic) configurations lead to very different pictures of the importance of the spreading width resulting in wavelet scales being a sensitive measure of their interplay. The relativistic approach with particle-vibration coupling, in particular, shows impressive agreement with the number and absolute values of the scales extracted from the experimental data. △ Less

Submitted 11 July, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: 10 pages, 4 figures

Bayesian inference of finite-nuclei observables based on the KIDS model

Authors: Jun Xu, Panagiota Papakonstantinou

Abstract: Bayesian analyses on both isoscalar and isovector nuclear interaction parameters are carried out based on the Korea-IBS-Daegu-SKKU (KIDS) model under the constraints of nuclear structure data of $^{208}$Pb and $^{120}$Sn. Under the constraint of the neutron-skin thickness, it is found that incorporating the curvature parameter $K_{sym}$ of nuclear symmetry energy as an independent variable signifi… ▽ More Bayesian analyses on both isoscalar and isovector nuclear interaction parameters are carried out based on the Korea-IBS-Daegu-SKKU (KIDS) model under the constraints of nuclear structure data of $^{208}$Pb and $^{120}$Sn. Under the constraint of the neutron-skin thickness, it is found that incorporating the curvature parameter $K_{sym}$ of nuclear symmetry energy as an independent variable significantly broadens the posterior probability distribution function (PDF) of the slope parameter $L$, and affects the related correlations. Typically, the anticorrelation between $L$ and the symmetry energy at saturation density disappears, while a positive correlation between $L$ and $K_{sym}$ is observed. Under the constraint of the isoscalar giant monopole resonance (ISGMR), incorporating the skewness parameter as an independent variable also significantly broadens the posterior PDF of the nuclear matter incompressibility $K_0$. Even with the broad uncertainties of higher-order parameters of the equation of state (EOS), robust constraints of $L<90$ MeV and $K_0<270$ MeV are obtained. Our study quantifies the consistency between the constraints on $L$ from the neutron-skin data of PREXII and isovector giant dipole resonance (IVGDR) data, and the constraints on $K_0$ from the data of ISGMR in $^{208}$Pb and $^{120}$Sn. △ Less

Submitted 19 February, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: 18 pages, 17 figures

Journal ref: Physical Review C 105, 044305 (2022)

Nuclear mass table in deformed relativistic Hartree-Bogoliubov theory in continuum: I. even-even nuclei

Authors: DRHBc Mass Table Collaboration, Kaiyuan Zhang, Myung-Ki Cheoun, Yong-Beom Choi, Pooi Seong Chong, Jianmin Dong, Zihao Dong, Xiaokai Du, Lisheng Geng, Eunja Ha, Xiao-Tao He, Chan Heo, Meng Chit Ho, Eun Jin In, Seonghyun Kim, Youngman Kim, Chang-Hwan Lee, Jenny Lee, Hexuan Li, Zhipan Li, Tianpeng Luo, Jie Meng, Myeong-Hwan Mun, Zhongming Niu, Cong Pan , et al. (22 additional authors not shown)

Abstract: Ground-state properties of even-even nuclei with $8\le Z\le120$ from the proton drip line to the neutron drip line have been investigated using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the density functional PC-PK1. With the effects of deformation and continuum included simultaneously, 2583 even-even nuclei are predicted to be bound. The calculated binding ener… ▽ More Ground-state properties of even-even nuclei with $8\le Z\le120$ from the proton drip line to the neutron drip line have been investigated using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the density functional PC-PK1. With the effects of deformation and continuum included simultaneously, 2583 even-even nuclei are predicted to be bound. The calculated binding energies, two-nucleon separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. The rms deviation from the 637 mass data is 1.518 MeV, providing one of the best microscopic descriptions for nuclear masses. The drip lines obtained from DRHBc calculations are compared with other calculations, including the spherical relativistic continuum Hartree-Bogoliubov (RCHB) and triaxial relativistic Hartree-Bogoliubov (TRHB) calculations with PC-PK1. The deformation and continuum effects on the limits of the nuclear landscape are discussed. Possible peninsulas consisting of bound nuclei beyond the two-neutron drip line are predicted. The systematics of the two-nucleon separation energies, two-nucleon gaps, rms radii, quadrupole deformations, potential energy curves, neutron densities, neutron mean-field potentials, and pairing energies in the DRHBc calculations are also discussed. In addition, the $αあるふぁ$ decay energies extracted are in good agreement with available data. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 217 pages, 15 figures, 2 tables, accepted for publication in Atomic Data and Nuclear Data Tables, data file in the TXT form is available for download under "Ancillary files"

Density dependence of the nuclear symmetry energy and neutron skin thickness in the KIDS framework

Authors: Panagiota Papakonstantinou

Abstract: The KIDS framework for the nuclear equation of state (EoS) and energy density functional (EDF) offers the possibility to explore systematically the effect of EoS parameters on predictions for a variety of observables. The EoS parameters can be varied independently of each other and independently of assumptions regarding the in-medium nucleon effective mass. Here I present a pilot study of the neut… ▽ More The KIDS framework for the nuclear equation of state (EoS) and energy density functional (EDF) offers the possibility to explore systematically the effect of EoS parameters on predictions for a variety of observables. The EoS parameters can be varied independently of each other and independently of assumptions regarding the in-medium nucleon effective mass. Here I present a pilot study of the neutron skin thickness (NST) in nuclei of current interest. The results indicate that variations of the symmetry energy slope parameter L by roughly 10 MeV and variations of the droplet-model counterpart of the curvature parameter $K_τたう$ by roughly 20 MeV affect predictions by comparable amounts. However, structural details may also have sizable effects on predictions, notably in the cases of $^{68}$Ni and $^{208}$Pb. This work is part of a systematic investigation of the NST within the KIDS framework and of a broader effort to constrain the density dependence of the nuclear symmetry energy. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: 6 pages incl. figure, 2 tables; prepared for HNPS Adv. Nucl. Phys. (Proc. 29th Symposium of the HNPS, September 2021); a complementary KIDS study with explicit neutron-star constraints can be found in arXiv:2112.00996

Constraints on the curvature of nuclear symmetry energy from recent astronomical data within the KIDS framework

Authors: Hana Gil, Panagiota Papakonstantinou, Chang Ho Hyun

Abstract: We investigate the density dependence of the nuclear symmetry energy $S(ρろー) $ in the KIDS (Korea-IBS-Daegu-SKKU) framework for the nuclear equation of state (EoS) and energy-density functional (EDF). The aim is to constrain the value of the curvature parameter ($K_{\rm sym}$) based on recent astronomical data. First, assuming a standard saturation point, we calculate bulk nuclear properties within… ▽ More We investigate the density dependence of the nuclear symmetry energy $S(ρろー) $ in the KIDS (Korea-IBS-Daegu-SKKU) framework for the nuclear equation of state (EoS) and energy-density functional (EDF). The aim is to constrain the value of the curvature parameter ($K_{\rm sym}$) based on recent astronomical data. First, assuming a standard saturation point, we calculate bulk nuclear properties within KIDS-EDF for different values of the compression modulus of symmetric nuclear matter ($K_0$) and of the leading-order symmetry energy parameters, i.e., the symmetry energy ($J$) and slope ($L$) at saturation density, each within a broadly accepted range, as well as $K_{\rm sym}$. All of the above EoS parameters are varied independently of each other. The skewness parameter ($Q_{\rm sym}$) is presently kept fixed at 650 MeV. For all EoS parameter sets which describe the selected nuclear data within better than $0.3\%$, we calculate the neutron-star equation of state and mass-radius relation and analyze the results in terms of Pearson correlation coefficients $r$. We find that the value of $K_{\rm sym}$ is strongly correlated with the radius of both a canonical and a massive star ($|r|>0.9$). If we impose that all known astronomical constraints on the neutron star radii must be satisfied, we deduce $-150 < K_{\rm sym}<0$. As a result, the symmetry energy as a function of the density is consistently found to have an inflection point at $ρろー_0<ρろー<2ρろー_0$. We take the opportunity to report that the neutron skin thickness of $^{208}$Pb shows no correlation at all with the neutron star radii ($|r|<0.1$), in contrast with studies which focus on the role of $L$ only. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: 24 pages including supplementary material, 2 tables and 3 figures

Constraining the density dependence of the symmetry energy with nuclear data and astronomical observations in the KIDS framework

Authors: Hana Gil, Young-Min Kim, Panagiota Papakonstantinou, Chang Ho Hyun

Abstract: The KIDS framework for the nuclear equation of state (EoS) and energy density functional (EDF) offers the possibility to explore symmetry-energy (SE) parameters such as J (value at saturation density), L (slope), Ksym (curvature) and so on independently of each other and of assumptions about the effective mass. Here we examine the performance of EoSs with different SE parameters in reproducing nuc… ▽ More The KIDS framework for the nuclear equation of state (EoS) and energy density functional (EDF) offers the possibility to explore symmetry-energy (SE) parameters such as J (value at saturation density), L (slope), Ksym (curvature) and so on independently of each other and of assumptions about the effective mass. Here we examine the performance of EoSs with different SE parameters in reproducing nuclear properties and astronomical observations in an effort to constrain especially L and Ksym or the droplet-model counterpart Ktau. Assuming a standard EoS for symmetric matter, we explore several points on the hyperplane of (J,L,Ksym or Ktau) values. For each point, the corresponding EDF parameters and a pairing parameter are obtained for applications in spherical even-even nuclei. This is the first application of KIDS EDFs with pairing correlations. The EoSs are tested successively on properties of closed-shell nuclei, along the Sn isotopic chain, and on astronomical observations, in a step-by-step process of elimination and correction. A small regime of best-performing parameters is determined. The results strongly suggest that Ksym is negative and no lower than -200MeV, that Ktau lies between roughly -400 and -300MeV and that L lies between 40 and 65MeV with L<55MeV more likely. Correlations between symmetry-energy parameters are critically discussed. Predictions for the position of the neutron drip line and the neutron skin thickness of selected nuclei are reported. They are only weakly affected by the choice of effective mass values. Parts of the drip line can be sensitive to the SE parameters. The results underscore the role of Ktau and of precise astronomical input. Better constraints are possible with precise fits to nuclear energies and, in the future, more-precise input from astronomy. △ Less

Submitted 29 April, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: 15 pages incl. 7 figures, 4 tables; this version of the preprint (v2) corresponds to the published version (some comments added, a few typos corrected)

Journal ref: Phys. Rev. C 103, 034330 (2021)

Daejeon16 interaction with contact-term corrections for heavy nuclear systems

Authors: Panagiota Papakonstantinou, James P. Vary, Youngman Kim

Abstract: The Daejeon16 two-nucleon interaction is employed in many-body approaches based on the mean-field approximation. The perturbative character of Daejeon16 is verified by comparing results for 16O from the Hartree-Fock (HF) approximation and from the no-core shell model and by examining the magnitude of perturbative corrections to the HF energy in light and heavy nuclei. In order to approximately des… ▽ More The Daejeon16 two-nucleon interaction is employed in many-body approaches based on the mean-field approximation. The perturbative character of Daejeon16 is verified by comparing results for 16O from the Hartree-Fock (HF) approximation and from the no-core shell model and by examining the magnitude of perturbative corrections to the HF energy in light and heavy nuclei. In order to approximately describe energies and radii across the nuclear chart, a phenomenological correction in the form of a two-plus-three-nucleon contact interaction is introduced. With fitted parameters we achieve a very good description of medium-mass nuclei in terms of energy and size and also in terms of the centroid energy of the giant monopole resonance and the dipole polarizability calculated within the random-phase approximation (RPA). Our results provide further justification for the use of Daejeon16 augmented with phenomenological corrections as an effective interaction of perturbative character in a variety of applications. △ Less

Submitted 9 July, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: v3: 29 pages incl. 5 figures, 4 tables; corresponds to published version

Journal ref: J. Phys. G: Nucl. Part. Phys. 48 085105 (2021)

Neutron drip line in the deformed relativistic Hartree-Bogoliubov theory in continuum: Oxygen to Calcium

Authors: Eun Jin In, Youngman Kim, Panagiota Papakonstantinou, Seung-Woo Hong

Abstract: The location of the neutron drip line, currently known for only the lightest elements, remains a fundamental question in nuclear physics. Its description is a challenge for microscopic nuclear energy density functionals, as it must take into account in a realistic way not only the nuclear potential, but also pairing correlations, deformation effects and coupling to the continuum. The recently deve… ▽ More The location of the neutron drip line, currently known for only the lightest elements, remains a fundamental question in nuclear physics. Its description is a challenge for microscopic nuclear energy density functionals, as it must take into account in a realistic way not only the nuclear potential, but also pairing correlations, deformation effects and coupling to the continuum. The recently developed deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) aims to provide a unified description of even-even nuclei throughout the nuclear chart. Here, the DRHBc with the successful density functional PC-PK1 is used to investigate whether and how deformation influences the prediction for the neutron drip-line location for even-even nuclei with 8<=Z<=20, where many isotopes are predicted deformed. The results are compared with those based on the spherical relativistic continuum Hartree-Bogoliubov (RCHB) theory and discussed in terms of shape evolution and the variational principle. It is found that the Ne and Ar drip-line nuclei are different after the deformation effect is included. The direction of the change is not necessarily towards an extended drip line, but rather depends on the evolution of the degree of deformation towards the drip line. Deformation effects as well as pairing and continuum effects treated in a consistent way can affect critically the theoretical description of the neutron drip-line location. △ Less

Submitted 2 July, 2020; originally announced July 2020.

Comments: 9 pages incl. 7 figures

Deformed relativistic Hartree-Bogoliubov theory in continuum with point coupling functional: examples of even-even Nd isotopes

Authors: Kaiyuan Zhang, Myung-Ki Cheoun, Yong-Beom Choi, Pooi Seong Chong, Jianmin Dong, Lisheng Geng, Eunja Ha, Xiaotao He, Chan Heo, Meng Chit Ho, Eun Jin In, Seonghyun Kim, Youngman Kim, Chang-Hwan Lee, Jenny Lee, Zhipan Li, Tianpeng Luo, Jie Meng, Myeong-Hwan Mun, Zhongming Niu, Cong Pan, Panagiota Papakonstantinou, Xinle Shang, Caiwan Shen, Guofang Shen , et al. (13 additional authors not shown)

Abstract: The aim of this work is to develop the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) theory based on the point-coupling density functionals and extend it to provide a unified description for all even-even nuclei in the nuclear chart by overcoming all possible challenges. The nuclear superfluidity is considered via Bogoliubov transformation. Densities and potentials are expan… ▽ More The aim of this work is to develop the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) theory based on the point-coupling density functionals and extend it to provide a unified description for all even-even nuclei in the nuclear chart by overcoming all possible challenges. The nuclear superfluidity is considered via Bogoliubov transformation. Densities and potentials are expanded in terms of Legendre polynomials to include the axial deformation degrees of freedom. Sophisticated relativistic Hartree-Bogoliubov equations in coordinate space are solved in the DiracWoods-Saxon basis to consider the continuum effects. Numerical checks are performed from light nuclei to heavy nuclei. The techniques to construct the DRHBc mass table for even-even nuclei are explored. The DRHBc theory is extended to study heavier nuclei beyond magnesium isotopes. Taking Nd isotopes as examples, the experimental binding energies, two-neutron separation energies, quadrupole deformations, and charge radii are reproduced rather well. The deformation and continuum play essential roles in the description of nuclear masses and prediction of drip-line nuclei. By examining the single-particle levels in the canonical basis and their contributions to the total density, the thickness of the neutron skin, the particles number in continuum, and the Coulomb barrier, the exotic structures including the neutron skin and the proton radioactivity are predicted. △ Less

Submitted 30 July, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: 45 pages, 17 figures, and 1 table, accepted version for Physical Review C

Journal ref: Phys. Rev. C 102, 024314 (2020)

Analysis of nuclear structure in a converging expansion scheme

Authors: Hana Gil, Young-Min Kim, Chang Ho Hyun, Panagiota Papakonstantinou, Yongseok Oh

Abstract: In the framework of the KIDS generalized energy density functional (EDF), the nuclear equation of state (EoS) is expressed as an expansion in powers of the Fermi momentum or the cubic root of the density ($ρろー^{1/3}$). Although an optimal number of converging terms was obtained in specific cases of fits to empirical data and pseudodata, the degree of convergence remains to be examined not only for h… ▽ More In the framework of the KIDS generalized energy density functional (EDF), the nuclear equation of state (EoS) is expressed as an expansion in powers of the Fermi momentum or the cubic root of the density ($ρろー^{1/3}$). Although an optimal number of converging terms was obtained in specific cases of fits to empirical data and pseudodata, the degree of convergence remains to be examined not only for homogeneous matter but also for finite nuclei. One goal of the present work is to validate the minimal and optimal number of EoS parameters required for the description of homogeneous nuclear matter over a wide range of densities relevant for astrophysical applications. The major goal is to examine the validity of the adopted expansion scheme for an accurate description of finite nuclei. To this end we vary the values of the high-order derivatives of the EoS, namely the skewness of the energy of symmetric nuclear matter and the kurtosis of the symmetry energy, at saturation and examine the relative importance of each term in $ρろー^{1/3}$ expansion for homogeneous matter. For given sets of EoS parameters determined in this way, we define equivalent Skyrme-type functionals and examine the convergence in the description of finite nuclei focusing on the masses and charge radii of closed-shell nuclei. The EoS of symmetric nuclear matter is found to be efficiently parameterized with only 3 parameters and the symmetry energy (or the energy of pure neutron matter) with 4 parameters when the EoS is expanded in the power series of the Fermi momentum. Higher-order EoS parameters do not produce any improvement, in practice, in the description of nuclear ground-state energies and charge radii, which means that they cannot be constrained by bulk properties of nuclei. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 15 pages, 10 figures

Journal ref: Phys. Rev. C 100, 014312 (2019)

Nuclear structure and the nucleon effective mass: explorations with the versatile KIDS functional

Authors: Panagiota Papakonstantinou, Hana Gil

Abstract: The connection from the structure and dynamics of atomic nuclei (finite nuclear system) to the nuclear equation of state (thermodynamic limit) is primarily made through nuclear energy-density functional (EDF) theory. Failure to describe both entities simultaneously within existing EDF frameworks means that we have either seriously misjudged the scope of EDF or not fully taken advantage of it. Ente… ▽ More The connection from the structure and dynamics of atomic nuclei (finite nuclear system) to the nuclear equation of state (thermodynamic limit) is primarily made through nuclear energy-density functional (EDF) theory. Failure to describe both entities simultaneously within existing EDF frameworks means that we have either seriously misjudged the scope of EDF or not fully taken advantage of it. Enter the versatile KIDS Ansatz, which is based on controlled, order-by-order extensions of the nuclear EDF with respect to the Fermi momentum and allows a direct mapping from a given, immutable equation of state to a convenient Skyrme pseudopotential for applications in finite nuclei. A recent proof-of-principle study of nuclear ground-states revealed the subversive role of the effective mass. Here we summarize the formalism and previous results and present further explorations related to giant resonances. As examples we consider the electric dipole polarizability of 68Ni and the giant monopole resonance (GMR) of heavy nuclei, particularly the "fluffiness" of 120Sn. We find that the choice of the effective mass parameters and that of the compression modulus affect the centroid energy of the GMR to comparable degrees. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Comments: 8 pages, incl. 2 figures; prepared for the Procs. 27th Annual Symposium of the Hellenic Nuclear Physics Society (HNPS2018), June 8-9 2018, Athens. Greece

Women Scientists Who Made Nuclear Astrophysics

Authors: C. V. Hampton, M. Lugaro, P. Papakonstantinou, P. G. Isar, B. Nordström, N. Özkan, M. Aliotta, A. Ćiprijanović, S. Curtis, M. Di Criscienzo, J. den Hartogh, A. S. Font, A. Kankainen, C. Kobayashi, C. Lederer-Woods, E. Niemczura, T. Rauscher, A. Spyrou, S. Van Eck, M. Yavahchova, W. Chantereau, S. E. de Mink, E. Kaiser, F. -K. Thielemann, C. Travaglio , et al. (2 additional authors not shown)

Abstract: Female role models reduce the impact on women of stereotype threat, i.e., of being at risk of conforming to a negative stereotype about one's social, gender, or racial group. This can lead women scientists to underperform or to leave their scientific career because of negative stereotypes such as, not being as talented or as interested in science as men. Sadly, history rarely provides role models… ▽ More Female role models reduce the impact on women of stereotype threat, i.e., of being at risk of conforming to a negative stereotype about one's social, gender, or racial group. This can lead women scientists to underperform or to leave their scientific career because of negative stereotypes such as, not being as talented or as interested in science as men. Sadly, history rarely provides role models for women scientists; instead, it often renders these women invisible. In response to this situation, we present a selection of twelve outstanding women who helped to develop nuclear astrophysics. △ Less

Submitted 25 August, 2018; originally announced September 2018.

Comments: 5 pages; to appear in Springer Proceedings in Physics (Proc. of Intl. Conf. "Nuclei in the Cosmos XV", LNGS Assergi, Italy, June 2018)

From homogeneous matter to finite nuclei: Role of the effective mass

Authors: Hana Gil, Panagiota Papakonstantinou, Chang Ho Hyun, Yongseok Oh

Abstract: Recent astronomical observations, nuclear-reaction experiments, and microscopic calculations have placed new constraints on the nuclear equation of state (EoS) and revealed that most nuclear structure models fail to satisfy those constraints upon extrapolation to infinite matter. A reverse procedure for imposing EoS constraints on nuclear structure has been elusive. Here we present for the first t… ▽ More Recent astronomical observations, nuclear-reaction experiments, and microscopic calculations have placed new constraints on the nuclear equation of state (EoS) and revealed that most nuclear structure models fail to satisfy those constraints upon extrapolation to infinite matter. A reverse procedure for imposing EoS constraints on nuclear structure has been elusive. Here we present for the first time a method to generate a microscopic energy density functional (EDF) for nuclei from a given immutable EoS. The method takes advantage of a natural Ansatz for homogeneous nuclear matter, the Kohn-Sham framework, and the Skyrme formalism. We apply it to the realistic nuclear EoS of Akmal-Pandharipande-Ravenhall and describe successfully closed-(sub)shell nuclei. In the process, we provide predictions for the neutron skin thickness of nuclei based directly on the given EoS. Crucially, bulk and static nuclear properties are found practically independent of the assumed effective mass value - a unique result in bridging EDF of finite and homogeneous systems in general. △ Less

Submitted 19 July, 2019; v1 submitted 29 May, 2018; originally announced May 2018.

Comments: v3 (some typos corrected, esp. Table I and Fig. 4): 8 pages, incl. 5 figures, 1 table

Journal ref: Phys. Rev. C 99, 064319 (2019)

Structures in the nuclear electric-dipole spectrum

Authors: P. Papakonstantinou

Abstract: I undertake an anatomy of the nuclear electric-dipole spectrum and discuss selected observations. First, using data and calculations on 208Pb as a representative case, I identify five conspicuous structures in the spectra of heavy nuclei: (1) the giant dipole resonance (GDR), (2) the dipole compression mode, (3) the isoscalar low-energy dipole mode (IS-LED),(4) the isoscalar mode in the region of… ▽ More I undertake an anatomy of the nuclear electric-dipole spectrum and discuss selected observations. First, using data and calculations on 208Pb as a representative case, I identify five conspicuous structures in the spectra of heavy nuclei: (1) the giant dipole resonance (GDR), (2) the dipole compression mode, (3) the isoscalar low-energy dipole mode (IS-LED),(4) the isoscalar mode in the region of the GDR, and (5) a concentration of "pygmy" dipole strength (PDS) between the latter two. Next, compiling the data from Texas A&M I show that the energy of (4) follows roughly a A^{-1/6} dependence on the mass number. Finally, I summarize recent theoretical studies on the PDS in exotic nuclei, which suggest a strong influence of shell effects and loose binding. △ Less

Submitted 28 November, 2017; originally announced November 2017.

Comments: 6 pages, 2 figures, Proc. VIII International Symposium on EXOtic Nuclei (EXON 2016), Kazan, Russia, September 5-10, 2016

Journal ref: In Proceedings of the International Symposium on Exotic Nuclei, Eds. Yu. E. Penionzhkevich and Yu. G. Sobolev, World Scientific 2017; pp. 278-283

Three-nucleon forces and superfluidity in neutron matter

Authors: P. Papakonstantinou, J. W. Clark

Abstract: The existence of superfluidity of the neutron component in the core of a neutron star, associated specifically with triplet $P-$wave pairing, is currently an open question that is central to interpretation of the observed cooling curves and other neutron-star observables. Ab initio theoretical calculations aimed at resolving this issue face unique challenges in the relevant high-density domain, wh… ▽ More The existence of superfluidity of the neutron component in the core of a neutron star, associated specifically with triplet $P-$wave pairing, is currently an open question that is central to interpretation of the observed cooling curves and other neutron-star observables. Ab initio theoretical calculations aimed at resolving this issue face unique challenges in the relevant high-density domain, which reaches beyond the saturation density of symmetrical nuclear matter. These issues include uncertainties in the three-nucleon (3N) interaction and in the effects of strong short-range correlations -- and more generally of in-medium modification of nucleonic self-energies and interactions. A survey of existing solutions to the gap equations in the triplet channel shows that the separate or combined impacts of 3N forces, coupled channels, and mass renormalization range from moderate to strong to devastating, thus motivating a detailed analysis of the competing effects. In the present work we track the effects of the 3N force and in-medium modifications in the representative case of the $^3P_2$ channel, based on the Argonne V18 two-nucleon (2N) interaction supplemented by 3N interactions of the Urbana IX family. Sensitivity of the results to the input interaction is clearly demonstrated, while consistency issues arise with respect to the simultaneous treatment of 3N forces and in-medium effects. We consider this pilot study as the first step towards a systematic and comprehensive exploration of coupled-channel $^3P F_2$ pairing using a broad range of 2N and 3N interactions from the current generation of refined semi-phenomenological models and models derived from chiral effective field theory. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: 17 pages, incl. 6 figures; intended for a special issue of the Journal of Low Temperature Physics on superfluidity and pairing phenomena

From chiral NN(N) interactions to giant and pygmy resonances via extended RPA

Authors: P. Papakonstantinou, R. Trippel, R. Roth

Abstract: The properties of giant and pygmy resonances are calculated starting from chiral two-and three-nucleon interactions. The aim is to assess the predictive power of modern Hamiltonians and especially the role of the three-nucleon force. Methods based on the random-phase approximation (RPA) provide an optimal description of the modes of interest with minimal computational requirements. Here we discuss… ▽ More The properties of giant and pygmy resonances are calculated starting from chiral two-and three-nucleon interactions. The aim is to assess the predictive power of modern Hamiltonians and especially the role of the three-nucleon force. Methods based on the random-phase approximation (RPA) provide an optimal description of the modes of interest with minimal computational requirements. Here we discuss the giant resonances (GRs) of 40,48Ca isotopes and their low-energy dipole response. A comparison with previous results obtained with a transformed Argonne V18 two-nucleon potential points to certain improvements. △ Less

Submitted 27 April, 2017; originally announced April 2017.

Comments: 4 pages, 1 figure; from a talk at the 2016 Zakopane Conference on Nuclear Physics, September 2016

Journal ref: Acta Phys. Pol. B 48 (2017) 537

Nuclear Energy Density Functional for KIDS

Authors: Hana Gil, Panagiota Papakonstantinou, Chang Ho Hyun, Tae-Sun Park, Yongseok Oh

Abstract: The density functional theory (DFT) is based on the existence and uniqueness of a universal functional $E[ρろー]$, which determines the dependence of the total energy on single-particle density distributions. However, DFT says nothing about the form of the functional. Our strategy is to first look at what we know, from independent considerations, about the analytical density dependence of the energy o… ▽ More The density functional theory (DFT) is based on the existence and uniqueness of a universal functional $E[ρろー]$, which determines the dependence of the total energy on single-particle density distributions. However, DFT says nothing about the form of the functional. Our strategy is to first look at what we know, from independent considerations, about the analytical density dependence of the energy of nuclear matter and then, for practical applications, to obtain an appropriate density-dependent effective interaction by reverse engineering. In a previous work on homogeneous matter, we identified the most essential terms to include in our "KIDS" functional, named after the early-stage participating institutes. We now present first results for finite nuclei, namely the energies and radii of $^{16,28}$O, $^{40,60}$Ca. △ Less

Submitted 15 March, 2017; v1 submitted 14 November, 2016; originally announced November 2016.

Comments: 4 pages, 1 figure, Contribution to the Proceedings of the 2016 Zakopane Conference on Nuclear Physics

Density dependence of the nuclear energy-density functional

Authors: Panagiota Papakonstantinou, Tae-Sun Park, Yeunhwan Lim, Chang Ho Hyun

Abstract: The explicit density (rho) dependence in the coupling coefficients of the non-relativistic nuclear energy-density functional (EDF) encodes effects of three-nucleon forces and dynamical correlations. The necessity for a coupling coefficient in the form of a small fractional power of rho is empirical and the power often chosen arbitrarily. Consequently, precision-oriented parameterisations risk over… ▽ More The explicit density (rho) dependence in the coupling coefficients of the non-relativistic nuclear energy-density functional (EDF) encodes effects of three-nucleon forces and dynamical correlations. The necessity for a coupling coefficient in the form of a small fractional power of rho is empirical and the power often chosen arbitrarily. Consequently, precision-oriented parameterisations risk overfitting and loss of predictive power. Observing that the Fermi momentum kF~rho^1/3 is a key variable in Fermi systems, we examine if a power hierarchy in kF can be inferred from the properties of homogeneous matter in a domain of densities which is relevant for nuclear structure and neutron stars. For later applications we want to determine an EDF that is of good quality but not overtrained. We fit polynomial and other functions of rho^1/3 to existing microscopic calculations of the energy of symmetric and pure neutron matter and analyze the fits. We select a form and parameter set which we found robust and examine the parameters' naturalness and the resulting extrapolations. A statistical analysis confirms that low-order terms like rho^1/3 and rho^2/3 are the most relevant ones. It also hints at a different power hierarchy for symmetric vs. pure neutron matter, supporting the need for more than one rho^a terms in non-relativistic EDFs. The EDF we propose accommodates adopted properties of nuclear matter near saturation. Importantly, its extrapolation to dilute or asymmetric matter reproduces a range of existing microscopic results, to which it has not been fitted. It also predicts neutron-star properties consistent with observations. The coefficients display naturalness. Once determined for homogeneous matter, EDFs of the present form can be mapped onto Skyrme-type ones for use in nuclei. The statistical analysis can be extended to higher orders and for different ab initio calculations. △ Less

Submitted 27 July, 2017; v1 submitted 14 June, 2016; originally announced June 2016.

Comments: 14 pages, 4 figures; v2: extended version

Journal ref: Phys. Rev. C 97, 014312 (2018)

Proton pygmy resonances: predictions for N=20 isotones

Authors: Y. Kim, P. Papakonstantinou

Abstract: We study theoretically the low-energy electric-dipole response of N=20 isotones. We present results from a quasiparticle random-phase approximation (QRPA) and a continuum random-phase approximation (CRPA), and we compare them with results for the mirror Z=20 nuclei. According to our analysis, enhanced E1 strength is expected energetically well below the giant dipole resonance in the proton-rich is… ▽ More We study theoretically the low-energy electric-dipole response of N=20 isotones. We present results from a quasiparticle random-phase approximation (QRPA) and a continuum random-phase approximation (CRPA), and we compare them with results for the mirror Z=20 nuclei. According to our analysis, enhanced E1 strength is expected energetically well below the giant dipole resonance in the proton-rich isotones. Large amounts of E1 strength in the asymmetric N=20 isotones are predicted, unlike their equally asymmetric Z=20 mirror nuclei, pointing unambiguously to the role of structural effects such as loose binding. A proton-skin oscillation could develop especially in 46Fe. The proper description of non localized threshold transitions and the nucleon effective mass in mean-field treatments may affect theoretical predictions. We call for systematic theoretical investigations to quantify the role bulk-matter properties, in anticipation of measurements of E1 transitions in proton-rich nuclei. △ Less

Submitted 2 March, 2016; v1 submitted 8 September, 2015; originally announced September 2015.

Comments: 10 pages, incl. 9 figures and 2 tables; v2: some rephrasing and clarifications, corrected Fig. 6

Journal ref: Eur.Phys.J.A52:176,2016

Isoscalar and neutron modes in the E1 spectra of Ni isotopes and the relevance of shell effects and the continuum

Authors: P. Papakonstantinou, H. Hergert, R. Roth

Abstract: We study theoretically the electric dipole transitions of even Ni isotopes at low energies, using the self-consistent quasi-particle random-phase approximation (RPA) with the D1S Gogny interaction and a continuum-RPA model with the SLy4 Skyrme force. We analyze isoscalar states, isovector states, and the dipole polarizability. We define a reference value for the polarizability, to remove a trivial… ▽ More We study theoretically the electric dipole transitions of even Ni isotopes at low energies, using the self-consistent quasi-particle random-phase approximation (RPA) with the D1S Gogny interaction and a continuum-RPA model with the SLy4 Skyrme force. We analyze isoscalar states, isovector states, and the dipole polarizability. We define a reference value for the polarizability, to remove a trivial dependence on the mass number. We compare our results with data and other calculations, with a focus on collective states, shell effects, and threshold transitions. Our results support the presence of a strong isoscalar transition, with little or moderate E1 strength, as a universal feature of ordinary nuclei. In moderately neutron-rich Ni isotopes, namely 68Ni and neighboring isotopes, this transition is found bimodal due to couplings with surface neutrons. An adequate treatment of the continuum states appears essential for describing suprathreshold E1 strength, especially beyond 68Ni. Very exotic isotopes (N>50) are found highly polarizable, with practically all their E1 strength in the continuum. The dipole polarizability and the neutron skin thickness are influenced by shell structure in different ways, so they can appear anticorrelated. A comparison with existing results for lighter (Ca) and heavier (Sn) nuclei suggests that the so-called pygmy dipole strength is influenced strongly by shell effects and that, partly for that reason, its isospin structure depends on the mass region. △ Less

Submitted 8 September, 2015; v1 submitted 16 June, 2015; originally announced June 2015.

Comments: 14 pages, incl. 10 figures (v2) and 1 table; (v2: clarifications and minor restructuring to improve readability; one extra figure)

Properties and significance of the surface dipole mode

Authors: P. Papakonstantinou

Abstract: A strong isoscalar dipole resonance is known to be excited in a variety of nuclei, including isospin symmetric ones, at approximately 6-7 MeV. A series of theoretical studies and accumulating experimental evidence support an interpretation of the above dipole resonance as an elementary surface vibration. Obviously, such a mode is potentially as interesting as any collective excitation for a variet… ▽ More A strong isoscalar dipole resonance is known to be excited in a variety of nuclei, including isospin symmetric ones, at approximately 6-7 MeV. A series of theoretical studies and accumulating experimental evidence support an interpretation of the above dipole resonance as an elementary surface vibration. Obviously, such a mode is potentially as interesting as any collective excitation for a variety of reasons. In addition, though, it is found to account for the observed isoscalar segment of pygmy dipole strength. As discussed here, this has important implications for pygmy-strength interpretations and searches for genuine neutron-skin oscillations. △ Less

Submitted 5 August, 2014; originally announced August 2014.

Comments: 4 pages, incl. 1 figure. Adapted from a presentation at ARIS2014 (June 1-6, Tokyo, Japan)

Second random-phase approximation, Thouless' theorem and the stability condition reexamined and clarified

Authors: P. Papakonstantinou

Abstract: It has been revealed through numerical calculations that the Second Random Phase Approximation (SRPA) with the Hartree-Fock solution as its reference state results in 1) spurious states at genuinely finite energy, contrary to common expectation, and 2) unstable solutions, which within the first-order Random Phase Approximation correspond to real low-energy collective vibrations. In the present wor… ▽ More It has been revealed through numerical calculations that the Second Random Phase Approximation (SRPA) with the Hartree-Fock solution as its reference state results in 1) spurious states at genuinely finite energy, contrary to common expectation, and 2) unstable solutions, which within the first-order Random Phase Approximation correspond to real low-energy collective vibrations. In the present work, these shortcomings of SRPA are shown to not contradict Thouless' theorem about the energy-weighted sum rule, and their origin is traced to the violation of the stability condition. A more general theorem is proven. Formal arguments are elucidated through numerical examples. Implications for the validity of SRPA are discussed. △ Less

Submitted 30 July, 2014; originally announced July 2014.

Comments: 11 pages, incl. 5 figures; to appear in Phys. Rev. C

Journal ref: Phys. Rev. C 90, 024305 (2014)

Isospin properties of electric dipole excitations in 48Ca

Authors: V. Derya, D. Savran, J. Endres, M. N. Harakeh, H. Hergert, J. H. Kelley, P. Papakonstantinou, N. Pietralla, V. Yu. Ponomarev, R. Roth, G. Rusev, A. P. Tonchev, W. Tornow, H. J. Wörtche, A. Zilges

Abstract: Two different experimental approaches were combined to study the electric dipole strength in the doubly-magic nucleus 48Ca below the neutron threshold. Real-photon scattering experiments using bremsstrahlung up to 9.9 MeV and nearly mono-energetic linearly polarized photons with energies between 6.6 and 9.51 MeV provided strength distribution and parities, and an (αあるふぁ,αあるふぁ'γがんま) experiment at E_αあるふぁ=136 MeV… ▽ More Two different experimental approaches were combined to study the electric dipole strength in the doubly-magic nucleus 48Ca below the neutron threshold. Real-photon scattering experiments using bremsstrahlung up to 9.9 MeV and nearly mono-energetic linearly polarized photons with energies between 6.6 and 9.51 MeV provided strength distribution and parities, and an (αあるふぁ,αあるふぁ'γがんま) experiment at E_αあるふぁ=136 MeV gave cross sections for an isoscalar probe. The unexpected difference observed in the dipole response is compared to calculations using the first-order random-phase approximation and points to an energy-dependent isospin character. A strong isoscalar state at 7.6 MeV was identified for the first time supporting a recent theoretical prediction. △ Less

Submitted 3 February, 2014; originally announced February 2014.

Comments: 6 pages, 5 figures, as accepted in Phys. Lett. B

Journal ref: Phys. Lett. B 730 (2014) 288

Low-energy electric dipole response of Sn isotopes

Authors: P. Papakonstantinou, H. Hergert, V. Yu. Ponomarev, R. Roth

Abstract: We study the low-energy dipole (LED) strength distribution along the Sn isotopic chain in both the isoscalar (IS) and the isovector (IV, or E1) electric channels, to provide testable predictions and guidance for new experiments with stable targets and radioactive beams. We use the self-consistent Quasi-particle Random-Phase Approximation (QRPA) with finite-range interactions and mainly the Gogny D… ▽ More We study the low-energy dipole (LED) strength distribution along the Sn isotopic chain in both the isoscalar (IS) and the isovector (IV, or E1) electric channels, to provide testable predictions and guidance for new experiments with stable targets and radioactive beams. We use the self-consistent Quasi-particle Random-Phase Approximation (QRPA) with finite-range interactions and mainly the Gogny D1S force. We analyze also the performance of a realistic two-body interaction supplemented by a phenomenological three-body contact term. We find that from N=50 and up to the N=82 shell closure (132Sn) the lowest-energy part of the IS-LED spectrum is dominated by a collective transition whose properties vary smoothly with neutron number and which cannot be interpreted as a neutron-skin oscillation. For the neutron-rich species this state contributes to the E1 strength below particle threshold, but much more E1 strength is carried by other, weak but numerous transitions around or above threshold. We find that strong structural changes in the spectrum take effect beyond N=82, namely increased LED strength and lower excitation energies. Our results with the Gogny interaction are compatible with existing data. On this basis we predict that a) the summed IS strength below particle threshold shall be of the same order of magnitude for N=50-82, b) the summed E1 strength up to approximately 12 MeV shall be similar for N=50-82 MeV, while c) the summed E1 strength below threshold shall be of the same order of magnitude for N ~ 64 - 82 and much weaker for the lighter, more-symmetric isotopes. We point out a general agreement of our results with other non-relativistic studies, the absence of a collective IS mode in some of those studies, and a possibly radical disagreement with relativistic models. △ Less

Submitted 10 March, 2014; v1 submitted 9 October, 2013; originally announced October 2013.

Comments: v2: 13 pages, incl. 10 figures; extended discussion, esp. of collectivity; revised Fig.9; to appear in Phys.Rev.C

Journal ref: Phys. Rev. C 89, 034306 (2014)

Densities and energies of nuclei in dilute matter

Authors: P. Papakonstantinou, J. Margueron, F. Gulminelli, Ad. R. Raduta

Abstract: We explore the ground-state properties of nuclear clusters embedded in a gas of nucleons with the help of Skyrme-Hartree-Fock microscopic calculations. Two alternative representations of clusters are introduced, namely coordinate-space and energy-space clusters. We parameterize their density profiles in spherical symmetry in terms of basic properties of the energy density functionals used and prop… ▽ More We explore the ground-state properties of nuclear clusters embedded in a gas of nucleons with the help of Skyrme-Hartree-Fock microscopic calculations. Two alternative representations of clusters are introduced, namely coordinate-space and energy-space clusters. We parameterize their density profiles in spherical symmetry in terms of basic properties of the energy density functionals used and propose an analytical, Woods-Saxon density profile whose parameters depend, not only on the composition of the cluster, but also of the nucleon gas. We study the clusters' energies with the help of the local-density approximation, validated through our microscopic results. We find that the volume energies of coordinate-space clusters are determined by the saturation properties of matter, while the surface energies are strongly affected by the presence of the gas. We conclude that both the density profiles and the cluster energies are strongly affected by the gas and discuss implications for the nuclear EoS and related perspectives. Our study provides a simple, but microscopically motivated modeling of the energetics of clusterized matter at subsaturation densities, for direct use in consequential applications of astrophysical interest. △ Less

Submitted 1 May, 2013; originally announced May 2013.

Comments: 20 pages, incl. 12 figures

Journal ref: Phys. Rev. C 88, 045805 (2013)

Giant Resonances based on Unitarily Transformed Two-Nucleon plus Phenomenological Three-Nucleon Interactions

Authors: A. Günther, P. Papakonstantinou, R. Roth

Abstract: We investigate giant resonances of spherical nuclei on the basis of the Argonne V18 potential after unitary transformation within the Similarity Renormalization Group or the Unitary Correlation Operator Method supplemented by a phenomenological three-body contact interaction. Such Hamiltonians can provide a good description of ground-state energies and radii within Hartree-Fock plus low-order many… ▽ More We investigate giant resonances of spherical nuclei on the basis of the Argonne V18 potential after unitary transformation within the Similarity Renormalization Group or the Unitary Correlation Operator Method supplemented by a phenomenological three-body contact interaction. Such Hamiltonians can provide a good description of ground-state energies and radii within Hartree-Fock plus low-order many-body perturbation theory. The standard Random Phase Approximation is applied here to calculate the isoscalar monopole, isovector dipole, and isoscalar quadrupole excitation modes of the 40Ca, 90Zr, and 208Pb nuclei. Thanks to the inclusion of the three-nucleon interaction and despite the minimal optimization effort, a reasonable agreement with experimental centroid energies of all three modes has been achieved. The role and scope of the Hartree-Fock reference state in RPA methods are discussed. △ Less

Submitted 3 September, 2014; v1 submitted 25 March, 2013; originally announced March 2013.

Comments: v2: 11 pages, incl. 3 figures; extended discussion and outlook; to appear in J.Phys.G

Some aspects of the phase diagram of nuclear matter relevant to compact stars

Authors: F. Gulminelli, Ad-R. Raduta, J. Margueron, P. Papakonstantinou, M. Oertel

Abstract: Dense matter as it can be found in core-collapse supernovae and neutron stars is expected to exhibit different phase transitions which impact the matter composition and the equation of state, with important consequences on the dynamics of core-collapse supernova explosion and on the structure of neutron stars. In this paper we will address the specific phenomenology of two of such transitions, nam… ▽ More Dense matter as it can be found in core-collapse supernovae and neutron stars is expected to exhibit different phase transitions which impact the matter composition and the equation of state, with important consequences on the dynamics of core-collapse supernova explosion and on the structure of neutron stars. In this paper we will address the specific phenomenology of two of such transitions, namely the crust-core solid-liquid transition at sub-saturation density, and the possible strange transition at super-saturation density in the presence of hyperonic degrees of freedom. Concerning the neutron star crust-core phase transition at zero and finite temperature, it will be shown that, as a consequence of the presence of long-range Coulomb interactions, a clusterized phase is expected which is not accessible in the grand-canonical ensemble. A specific quasi-particle model will be introduced and some quantitative results relevant for the supernova dynamics will be shown. The opening of hyperonic degrees of freedom at higher densities corresponding to the neutron stars core also modifies the equation of state. The general characteristics and order of phase transitions in this regime will be analyzed in the framework of a self-consistent mean-field approach. △ Less

Submitted 15 February, 2013; originally announced February 2013.

Comments: arXiv admin note: substantial text overlap with arXiv:1206.4924, arXiv:1301.6955

Journal ref: International Summer School for Advanced Studies 'Dynamics of open nuclear systems' (PREDEAL12), Predeal : Romania (2012)

Phase diagram of neutron-rich nuclear matter and its impact on astrophysics

Authors: F. Gulminelli, Ad. R. Raduta, J. Margueron, P. Papakonstantinou, M. Oertel

Abstract: Dense matter as it can be found in core-collapse supernovae and neutron stars is expected to exhibit different phase transitions which impact the matter composition and equation of state, with important consequences on the dynamics of core-collapse supernova explosion and on the structure of neutron stars. In this paper we will address the specific phenomenology of two of such transitions, namely… ▽ More Dense matter as it can be found in core-collapse supernovae and neutron stars is expected to exhibit different phase transitions which impact the matter composition and equation of state, with important consequences on the dynamics of core-collapse supernova explosion and on the structure of neutron stars. In this paper we will address the specific phenomenology of two of such transitions, namely the crust-core solid-liquid transition at sub-saturation density, and the possible strange transition at super-saturation density in the presence of hyperonic degrees of freedom. Concerning the neutron star crust-core phase transition at zero and finite temperature, it will be shown that, as a consequence of the presence of long-range Coulomb interactions, the equivalence of statistical ensembles is violated and a clusterized phase is expected which is not accessible in the grand-canonical ensemble. A specific quasi-particle model will be introduced to illustrate this anomalous thermodynamics and some quantitative results relevant for the supernova dynamics will be shown. The opening of hyperonic degrees of freedom at higher densities corresponding to the neutron stars core modifies the equation of state. The general characteristics and order of phase transitions in this regime will be analyzed in the framework of a self-consistent mean-field approach. △ Less

Submitted 3 September, 2012; originally announced September 2012.

Comments: Invited Talk given at the 11th International Conference on Nucleus-Nucleus Collisions (NN2012), San Antonio, Texas, USA, May 27-June 1, 2012. To appear in the NN2012 Proceedings in Journal of Physics: Conference Series (JPCS)

Low-energy dipole strength and the critical case of 48Ca

Authors: P. Papakonstantinou, H. Hergert, V. Yu. Ponomarev, R. Roth

Abstract: Recent theoretical work has not led to a consensus regarding the nature of the low-energy E1 strength in the 40,44,48Ca isotopes, for which high-resolution (gamma,gamma') data exist. Here we revisit this problem using the first-order quasi-particle random-phase approximation (QRPA) and different interactions. First we examine all even Ca isotopes with N=14-40. All isotopes are predicted to undergo… ▽ More Recent theoretical work has not led to a consensus regarding the nature of the low-energy E1 strength in the 40,44,48Ca isotopes, for which high-resolution (gamma,gamma') data exist. Here we revisit this problem using the first-order quasi-particle random-phase approximation (QRPA) and different interactions. First we examine all even Ca isotopes with N=14-40. All isotopes are predicted to undergo dipole transitions at low energy, of large and comparable isoscalar strength but of varying E1 strength. Provided a moderate and uniform energetic shift is introduced to the results, QRPA with the Gogny D1S interaction is able to account for the (gamma,gamma') data, because, up to N=28, it yields a rather pure isoscalar oscillation. A neutron-skin oscillation is anticipated for N larger or equal to 30. This contradicts existing predictions that 44,48Ca develop a neutron-skin mode. Which theoretical result is correct cannot be resolved conclusively using the available data. We propose that alpha-scattering, possibly followed by an electroexcitation experiment, could resolve the situation and thereby help to improve the different models aspiring to describe reliably the low-energy dipole strength of nuclei. △ Less

Submitted 5 October, 2011; originally announced October 2011.

Comments: 7 pages, incl. 3 figures; PLB submitted

Journal ref: Phys. Lett. B 709 (2012) 270

Quasiparticle Random Phase Approximation with Interactions from the Similarity Renormalization Group

Authors: H. Hergert, P. Papakonstantinou, R. Roth

Abstract: We have developed a fully consistent framework for calculations in the Quasiparticle Random Phase Approximation (QRPA) with $NN$ interactions from the Similarity Renormalization Group (SRG) and other unitary transformations of realistic interactions. The consistency of our calculations, which use the same Hamiltonian to determine the Hartree-Fock-Bogoliubov (HFB) ground states and the residual int… ▽ More We have developed a fully consistent framework for calculations in the Quasiparticle Random Phase Approximation (QRPA) with $NN$ interactions from the Similarity Renormalization Group (SRG) and other unitary transformations of realistic interactions. The consistency of our calculations, which use the same Hamiltonian to determine the Hartree-Fock-Bogoliubov (HFB) ground states and the residual interaction for QRPA, guarantees an excellent decoupling of spurious strength, without the need for empirical corrections. While work is under way to include SRG-evolved 3N interactions, we presently account for some 3N effects by means of a linearly density-dependent interaction, whose strength is adjusted to reproduce the charge radii of closed-shell nuclei across the whole nuclear chart. As a first application, we perform a survey of the monopole, dipole, and quadrupole response of the calcium isotopic chain and of the underlying single-particle spectra, focusing on how their properties depend on the SRG parameter $λらむだ$. Unrealistic spin-orbit splittings suggest that spin-orbit terms from the 3N interaction are called for. Nevertheless, our general findings are comparable to results from phenomenological QRPA calculations using Skyrme or Gogny energy density functionals. Potentially interesting phenomena related to low-lying strength warrant more systematic investigations in the future. △ Less

Submitted 19 June, 2011; v1 submitted 1 April, 2011; originally announced April 2011.

Comments: 18 pages, 17 figures, 3 tables (RevTeX 4.1), v2: fixed typos & figures, as published

Journal ref: Phys.Rev.C83:064317,2011

Isoscalar dipole coherence at low energies and forbidden E1 strength

Authors: P. Papakonstantinou, V. Yu. Ponomarev, R. Roth, J. Wambach

Abstract: In 16O and 40Ca an isoscalar, low-energy dipole transition (IS-LED) exhausting approximately 4% of the isoscalar dipole (ISD) energy-weighted sum rule is experimentally known, but conspicuously absent from recent theoretical investigations of ISD strength. The IS-LED mode coincides with the so-called isospin-forbidden E1 transition. We report that for N=Z nuclei up to 100Sn the fully self-consiste… ▽ More In 16O and 40Ca an isoscalar, low-energy dipole transition (IS-LED) exhausting approximately 4% of the isoscalar dipole (ISD) energy-weighted sum rule is experimentally known, but conspicuously absent from recent theoretical investigations of ISD strength. The IS-LED mode coincides with the so-called isospin-forbidden E1 transition. We report that for N=Z nuclei up to 100Sn the fully self-consistent Random-Phase-Approximation with finite-range forces, phenomenological and realistic, yields a collective IS-LED mode, typically overestimating its excitation energy, but correctly describing its IS strength and electroexcitation form factor. The presence of E1 strength is solely due to the Coulomb interaction between the protons and the resulting isospin-symmetry breaking. The smallness of its value is related to the form of the transition density, due to translational invariance. The calculated values of E1 and ISD strength carried by the IS-LED depend on the effective interaction used. Attention is drawn to the possibility that in N-not-equal-Z nuclei this distinct mode of IS surface vibration can develop as such or mix strongly with skin modes and thus influence the pygmy dipole strength as well as the ISD strength function. In general, theoretical models currently in use may be unfit to predict its precise position and strength, if at all its existence. △ Less

Submitted 4 November, 2010; originally announced November 2010.

Comments: 9 pages, 6 figures, EPJA submitted

Journal ref: Eur.Phys.J.A47:14,2011

Large-scale second RPA calculations with finite-range interactions

Authors: P. Papakonstantinou, R. Roth

Abstract: Second RPA (SRPA) calculations of nuclear response are performed and analyzed. Unlike in most other SRPA applications, the ground state, approximated by the Hartree-Fock (HF) ground state, and the residual couplings are described by the same Hamiltonian and no arbitrary truncations are imposed on the model space. Finite-range interactions are used and thus divergence problems are not present. We… ▽ More Second RPA (SRPA) calculations of nuclear response are performed and analyzed. Unlike in most other SRPA applications, the ground state, approximated by the Hartree-Fock (HF) ground state, and the residual couplings are described by the same Hamiltonian and no arbitrary truncations are imposed on the model space. Finite-range interactions are used and thus divergence problems are not present. We employ a realistic interaction, derived from the Argonne V18 potenial using the unitary correlation operator method (UCOM), as well as the simple Brink-Boeker interaction. Representative results are discussed, mainly on giant resonances and low-lying collective states. The focus of the present work is not on the comparison with data, but rather on technical and physical aspects of the method. We present how the large-scale eigenvalue problem that SRPA entails can be treated, and demonstrate how the method operates in producing self-energy corrections and fragmentation. The so-called diagonal approximation is conditionally validated. Stability problems are traced back to missing ground-state correlations. △ Less

Submitted 9 October, 2009; originally announced October 2009.

Comments: 13 pages, incl. 9 figures, 1 table

Journal ref: Phys.Rev.C81:024317,2010

Second RPA and correlated realistic interactions

Authors: P. Papakonstantinou, R. Roth

Abstract: We examine the response of closed-shell nuclei using a correlated interaction, derived with the Unitary Correlation Operator Method (UCOM) from the Argonne V18 potential, in second RPA (SRPA) calculations. The same correlated two-body interaction is used to derive the Hartree-Fock ground state and the SRPA equations. Our results show that the coupling of particle-hole states to higher-order conf… ▽ More We examine the response of closed-shell nuclei using a correlated interaction, derived with the Unitary Correlation Operator Method (UCOM) from the Argonne V18 potential, in second RPA (SRPA) calculations. The same correlated two-body interaction is used to derive the Hartree-Fock ground state and the SRPA equations. Our results show that the coupling of particle-hole states to higher-order configurations produces sizable effects compared with first-order RPA. A much improved description of the isovector dipole and isoscalar quadrupole resonances is obtained, thanks in part to the more fundamental treatment of the nucleon effective mass offered by SRPA. The present work suggests the prospect of describing giant resonance properties realistically and consistently within extended RPA theories. Self-consistency issues of the present SRPA method and residual three-body effects are pointed out. △ Less

Submitted 27 May, 2008; originally announced May 2008.

Comments: 6 pages, incl. 5 figures

Journal ref: Phys. Lett. B 671 (2009) 356

Giant Resonances using Correlated Realistic Interactions: The Case for Second RPA

Authors: P. Papakonstantinou, R. Roth

Abstract: Lately we have been tackling the problem of describing nuclear collective excitations starting from correlated realistic nucleon-nucleon (NN) interactions. The latter are constructed within the Unitary Correlation Operator Method (UCOM), starting from realistic NN potentials. It has been concluded that first-order RPA with a two-body UCOM interaction is not capable, in general, of reproducing qu… ▽ More Lately we have been tackling the problem of describing nuclear collective excitations starting from correlated realistic nucleon-nucleon (NN) interactions. The latter are constructed within the Unitary Correlation Operator Method (UCOM), starting from realistic NN potentials. It has been concluded that first-order RPA with a two-body UCOM interaction is not capable, in general, of reproducing quantitatively the properties of giant resonances (GRs), due to missing higher-order configurations and long-range correlations as well as neglected three-body terms in the Hamiltonian. Here we report results on GRs obtained by employing a UCOM interaction based on the Argonne V18 potential in Second RPA (SRPA) calculations. The same interaction is used to describe the Hartree-Fock (HF) ground state and the residual interactions. We find that the inclusion of second-order configurations -- which effectively dress the underlying HF single-particle states with self-energy insertions -- produces sizable corrections. The effect appears essential for a realistic description of GRs when using the UCOM. We argue that effects of higher than second order should be negligible. Therefore, the UCOM-SRPA emerges as a promising tool for consistent calculations of collective states in closed-shell nuclei. This is an interesting development, given that SRPA can accommodate more physics than RPA (e.g., fragmentation). Remaining discrepancies due to the missing three-body terms and self-consistency issues of the present SRPA model are pointed out. △ Less

Submitted 20 September, 2007; originally announced September 2007.

Comments: 6 pages, incl. 1 figure; Proc. 26th Int. Workshop on Nuclear Theory, June 2007, Rila mountains, Bulgaria

Reduction of the RPA eigenvalue problem and a generalized Cholesky decomposition for real-symmetric matrices

Authors: P. Papakonstantinou

Abstract: The particular symmetry of the random-phase-approximation (RPA) matrix has been utilized in the past to reduce the RPA eigenvalue problem into a symmetric-matrix problem of half the dimension. The condition of positive definiteness of at least one of the matrices A+-B has been imposed (where A and B are the submatrices of the RPA matrix) so that, e.g., its square root can be found by Cholesky de… ▽ More The particular symmetry of the random-phase-approximation (RPA) matrix has been utilized in the past to reduce the RPA eigenvalue problem into a symmetric-matrix problem of half the dimension. The condition of positive definiteness of at least one of the matrices A+-B has been imposed (where A and B are the submatrices of the RPA matrix) so that, e.g., its square root can be found by Cholesky decomposition. In this work, alternative methods are pointed out to reduce the RPA problem to a real (not symmetric, in general) problem of half the dimension, with the condition of positive definiteness relaxed. One of the methods relies on a generalized Cholesky decomposition, valid for non-singular real symmetric matrices. The algorithm is described and a corresponding routine in C is given. △ Less

Submitted 8 February, 2007; originally announced February 2007.

Comments: 5 pages, to appear in EPL, http://www.iop.org/EJ/journal/EPL

Journal ref: Europhys.Lett.78:12001,2007

Nuclear collective excitations using correlated realistic interactions: the role of explicit RPA correlations

Authors: P. Papakonstantinou, R. Roth, N. Paar

Abstract: We examine to which extent correlated realistic nucleon-nucleon interactions, derived within the Unitary Correlation Operator Method (UCOM), can describe nuclear collective motion in the framework of first-order random-phase approximation (RPA). To this end we employ the correlated Argonne V18 interaction in calculations within the so-called "Extended" RPA (ERPA) and investigate the response of… ▽ More We examine to which extent correlated realistic nucleon-nucleon interactions, derived within the Unitary Correlation Operator Method (UCOM), can describe nuclear collective motion in the framework of first-order random-phase approximation (RPA). To this end we employ the correlated Argonne V18 interaction in calculations within the so-called "Extended" RPA (ERPA) and investigate the response of closed-shell nuclei. The ERPA is a renormalized RPA version which considers explicitly the depletion of the Fermi sea due to long-range correlations and thus allows us to examine how these affect the excitation spectra. It is found that the effect on the properties of giant resonances is rather small. Compared to the standard RPA, where excitations are built on top of the uncorrelated Hartree-Fock (HF) ground state, their centroid energies decrease by up to 1 MeV, approximately, in the isovector channel. The isoscalar response is less affected in general. Thus, the disagreement between our previous UCOM-based RPA calculations and the experimental data are to be attributed to other effects, mainly to a residual three-body force and higher-order configurations. Ground-state properties obtained within the ERPA are compared with corresponding HF and perturbation-theory results and are discussed as well. The ERPA formalism is presented in detail. △ Less

Submitted 15 September, 2006; originally announced September 2006.

Comments: 13 pages, incl. 7 figures and 1 table

Journal ref: Phys.Rev.C75:014310,2007

Nuclear Structure in the UCOM Framework: From Realistic Interactions to Collective Excitations

Authors: R. Roth, H. Hergert, N. Paar, P. Papakonstantinou

Abstract: The Unitary Correlation Operator Method (UCOM) provides a means for nuclear structure calculations starting from realistic NN potentials. The dominant short-range central and tensor correlations are described explicitly by a unitary transformation. The application of UCOM in the context of the no-core shell model provides insight into the interplay between dominant short-range and residual long-… ▽ More The Unitary Correlation Operator Method (UCOM) provides a means for nuclear structure calculations starting from realistic NN potentials. The dominant short-range central and tensor correlations are described explicitly by a unitary transformation. The application of UCOM in the context of the no-core shell model provides insight into the interplay between dominant short-range and residual long-range correlations in the nuclear many-body problem. The use of the correlated interaction within Hartree-Fock, many-body perturbation theory, and Random Phase Approximation gives access to various nuclear structure observables throughout the nuclear chart. △ Less

Submitted 7 August, 2006; originally announced August 2006.

Comments: 9 pages, 3 figures, invited talk at the 2nd Int. Conf. on "Collective Motion in Nuclei under Extreme Conditions" (COMEX 2), Sankt Goar, Germany

Journal ref: Nucl.Phys.A788:12-19,2007

Correlation energies in the random phase approximation using realistic interactions

Authors: C. Barbieri, N. Paar, R. Roth, P. Papakonstantinou

Abstract: The self-consistent random phase approximation (RPA) based on a correlated realistic nucleon-nucleon interaction is used to evaluate correlation energies in closed-shell nuclei beyond the Hartree-Fock level. The relevance of contributions associated with charge exchange excitations as well as the necessity to correct for the double counting of the second order contribution to the RPA ring summat… ▽ More The self-consistent random phase approximation (RPA) based on a correlated realistic nucleon-nucleon interaction is used to evaluate correlation energies in closed-shell nuclei beyond the Hartree-Fock level. The relevance of contributions associated with charge exchange excitations as well as the necessity to correct for the double counting of the second order contribution to the RPA ring summation are emphasized. Once these effects are properly accounted for, the RPA ring summation provides an efficient tool to assess the impact of long-range correlations on binding energies throughout the whole nuclear chart, which is of particular importance when starting from realistic interactions. △ Less

Submitted 3 August, 2006; originally announced August 2006.

Comments: 4 pages, 3 figs., revtex