Symbolic computation of the Hartree-Fock energy from a chiral EFT three-nucleon interaction at N2LO

Abstract We present the first of a two-part Mathematica notebook collection that implements a symbolic approach for the application of the density matrix expansion (DME) to the Hartree–Fock (HF) energy from a chiral effective field theory (EFT) three-nucleon interaction at N2LO. The final output from the notebooks is a Skyrme-like energy density functional that provides a quasi-local approximation to the non-local HF energy. In this paper, we discuss the derivation of the HF energy and its simplification in terms of the scalar/vector–isoscalar/isovector parts of the one-body density matrix. Furthermore, a set of steps is described and illustrated on how to extend the approach to other three-nucleon interactions. Program summary Program title: SymbHFNNN Catalogue identifier: AEGC_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEGC_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 96 666 No. of bytes in distributed program, including test data, etc.: 378 083 Distribution format: tar.gz Programming language: Mathematica 7.1 Computer: Any computer running Mathematica 6.0 and later versions Operating system: Windows Xp, Linux/Unix RAM: 256 Mb Classification: 5, 17.16, 17.22 Nature of problem: The calculation of the HF energy from the chiral EFT three-nucleon interaction at N2LO involves tremendous spin–isospin algebra. The problem is compounded by the need to eventually obtain a quasi-local approximation to the HF energy, which requires the HF energy to be expressed in terms of scalar/vector–isoscalar/isovector parts of the one-body density matrix. The Mathematica notebooks discussed in this paper solve the latter issue. Solution method: The HF energy from the chiral EFT three-nucleon interaction at N2LO is cast into a form suitable for an automatic simplification of the spin–isospin traces. Several Mathematica functions and symbolic manipulation techniques are used to obtain the result in terms of the scalar/vector–isoscalar/isovector parts of the one-body density matrix. Running time: Several hours

[1]  H. Witała,et al.  Three-nucleon forces from chiral effective field theory , 2002, nucl-th/0208023.

[2]  E. Epelbaum,et al.  The Two-nucleon system at next-to-next-to-next-to-leading order , 2004, nucl-th/0405048.

[3]  Isovector splitting of nucleon effective masses, ab-initio benchmarks and extended stability criteria for Skyrme energy functionals , 2006, nucl-th/0607065.

[4]  Paul-Henri Heenen,et al.  Self-consistent mean-field models for nuclear structure , 2003 .

[5]  B. Gebremariam,et al.  Improved density matrix expansion for spin-unsaturated nuclei , 2009, 0910.4979.

[6]  S. Bogner,et al.  Model independent low momentum nucleon interaction from phase shift equivalence , 2003, nucl-th/0305035.

[7]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[8]  J. Dobaczewski,et al.  Local density approximation for proton-neutron pairing correlations: Formalism , 2004 .

[9]  Low-momentum interactions with smooth cutoffs , 2006, nucl-th/0609003.

[10]  H. Sagawa,et al.  Effective pairing interactions with isospin density dependence , 2007, 0712.3644.

[11]  T. Skyrme.  CVII. The nuclear surface , 1956 .

[12]  H. Sagawa,et al.  Extended Skyrme interaction: I. Spin fluctuations in dense matter , 2009, 0905.1931.

[13]  P. Ring,et al.  Relativistic nuclear energy density functionals: Adjusting parameters to binding energies , 2008, 0809.1375.

[14]  D. Brink,et al.  Time-dependent hartree-fock theory with Skyrme's interaction , 1975 .

[15]  D. R. Entem,et al.  Accurate charge dependent nucleon nucleon potential at fourth order of chiral perturbation theory , 2003 .

[16]  Stephen Wolfram,et al.  The Mathematica Book , 1996 .

[17]  T. Lesinski,et al.  The tensor part of the Skyrme energy density functional. I. Spherical nuclei , 2007, 0704.0731.

[18]  B. Gebremariam,et al.  Symbolic integration of a product of two spherical Bessel functions with an additional exponential and polynomial factor , 2010, Comput. Phys. Commun..

[19]  Markus Kortelainen,et al.  Local nuclear energy density functional at next-to-next-to-next-to-leading order , 2008 .

[20]  John W. Negele,et al.  Density-matrix expansion for an effective nuclear Hamiltonian , 1972 .

[21]  S. Goriely,et al.  First Gogny-Hartree-Fock-Bogoliubov nuclear mass model. , 2009, Physical review letters.