skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Testing the density matrix expansion against ab initio calculations of trapped neutron drops

Abstract

Microscopic input to a universal nuclear energy density functional can be provided through the density matrix expansion (DME), which has recently been revived and improved. Several DME implementation strategies are tested for neutron drop systems in harmonic traps by comparing to Hartree-Fock (HF) and ab initio no-core full configuration (NCFC) calculations with a model interaction (Minnesota potential). The new DME with exact treatment of Hartree contributions is found to best reproduce HF results and supplementing the functional with fit Skyrme-like contact terms shows systematic improvement toward the full NCFC results.

Authors:
 [1];  [1];  [2];  [3];  [4];  [5];  [5]
  1. Michigan State University, East Lansing
  2. Ohio State University
  3. ORNL
  4. University of Tennessee, Knoxville (UTK) & Oak Ridge National Laboratory (ORNL)
  5. Iowa State University
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1037058
DOE Contract Number:
DE-AC05-00OR22725
Resource Type:
Journal Article
Resource Relation:
Journal Name: Physical Review, C (Nuclear Physics); Journal Volume: 84; Journal Issue: 4
Country of Publication:
United States
Language:
English
Subject:
73 NUCLEAR PHYSICS AND RADIATION PHYSICS; CONFIGURATION; DENSITY MATRIX; FUNCTIONALS; HARMONICS; IMPLEMENTATION; NEUTRONS; NUCLEAR ENERGY; TESTING

Citation Formats

Bogner, S. K., Hergert, H., Furnstahl, R. J., Kortelainen, Erno M, Stoitsov, M. V., Maris, Pieter, and Vary, J. P. Testing the density matrix expansion against ab initio calculations of trapped neutron drops. United States: N. p., 2011. Web.
Bogner, S. K., Hergert, H., Furnstahl, R. J., Kortelainen, Erno M, Stoitsov, M. V., Maris, Pieter, & Vary, J. P. Testing the density matrix expansion against ab initio calculations of trapped neutron drops. United States.
Bogner, S. K., Hergert, H., Furnstahl, R. J., Kortelainen, Erno M, Stoitsov, M. V., Maris, Pieter, and Vary, J. P. 2011. "Testing the density matrix expansion against ab initio calculations of trapped neutron drops". United States. doi:.
@article{osti_1037058,
title = {Testing the density matrix expansion against ab initio calculations of trapped neutron drops},
author = {Bogner, S. K. and Hergert, H. and Furnstahl, R. J. and Kortelainen, Erno M and Stoitsov, M. V. and Maris, Pieter and Vary, J. P.},
abstractNote = {Microscopic input to a universal nuclear energy density functional can be provided through the density matrix expansion (DME), which has recently been revived and improved. Several DME implementation strategies are tested for neutron drop systems in harmonic traps by comparing to Hartree-Fock (HF) and ab initio no-core full configuration (NCFC) calculations with a model interaction (Minnesota potential). The new DME with exact treatment of Hartree contributions is found to best reproduce HF results and supplementing the functional with fit Skyrme-like contact terms shows systematic improvement toward the full NCFC results.},
doi = {},
journal = {Physical Review, C (Nuclear Physics)},
number = 4,
volume = 84,
place = {United States},
year = 2011,
month = 1
}
  • Cited by 23
  • An alternative approach to open-shell molecular calculations using the variational two-electron reduced-density-matrix (2-RDM) theory [Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)] is presented. The energy and 2-RDM of the open-shell molecule (or radical) are computed from the limit of dissociating one or more hydrogen atoms from a molecule in a singlet state. Because the ground-state energy of an 'infinitely' separated hydrogen atom in a given finite basis is known, we can determine the energy of the radical by subtracting the energy of one or more hydrogen atoms from the energy of the total dissociated system. The 2-RDM is constrained tomore » have singlet symmetry in all calculations. Two sets of N-representability conditions are employed: (i) two-positivity conditions, and (ii) two-positivity conditions plus the T{sub 2} condition, which is a subset of the three-positivity conditions. Optimization of the energy with respect to the 2-RDM is performed with a first-order algorithm for solving the semidefinite program within the variational 2-RDM method. We present calculations of several radicals near equilibrium as well as the dissociation curves of the diatomic radicals CH and OH.« less
  • Spin Orbit Coupling (SOC) is introduced to molecular ab initio density matrix renormalization group (DMRG) calculations. In the presented scheme, one first approximates the electronic ground state and a number of excited states of the Born-Oppenheimer (BO) Hamiltonian with the aid of the DMRG algorithm. Owing to the spin-adaptation of the algorithm, the total spin S is a good quantum number for these states. After the non-relativistic DMRG calculation is finished, all magnetic sublevels of the calculated states are constructed explicitly, and the SOC operator is expanded in the resulting basis. To this end, spin orbit coupled energies and wavefunctionsmore » are obtained as eigenvalues and eigenfunctions of the full Hamiltonian matrix which is composed of the SOC operator matrix and the BO Hamiltonian matrix. This treatment corresponds to a quasi-degenerate perturbation theory approach and can be regarded as the molecular equivalent to atomic Russell-Saunders coupling. For the evaluation of SOC matrix elements, the full Breit-Pauli SOC Hamiltonian is approximated by the widely used spin-orbit mean field operator. This operator allows for an efficient use of the second quantized triplet replacement operators that are readily generated during the non-relativistic DMRG algorithm, together with the Wigner-Eckart theorem. With a set of spin-orbit coupled wavefunctions at hand, the molecular g-tensors are calculated following the scheme proposed by Gerloch and McMeeking. It interprets the effective molecular g-values as the slope of the energy difference between the lowest Kramers pair with respect to the strength of the applied magnetic field. Test calculations on a chemically relevant Mo complex demonstrate the capabilities of the presented method.« less
  • The phase diagram of the two-leg t-J{sub z} ladder is explored, using the density-matrix renormalization group method. Results are obtained for energy gaps, electron density profiles, and correlation functions for the half filled and quarter filled cases. The effective Lagrangian velocity parameter v{sub {rho}} is shown to vanish at half filling. The behavior of the one-hole gap in the Nagaoka limit is investigated, and found to disagree with theoretical predictions. A tentative phase diagram is presented, which is quite similar to the full t-J ladder, but scaled up by a factor of about 2 in coupling. Near half filling amore » Luther-Emery phase is found, which may be expected to show superconducting correlations, while near quarter filling the system appears to be in a Tomonaga-Luttinger phase.« less