A Fortran 90 program to solve the Hartree-Fock equations for interacting spin-1/2 Fermions confined in Harmonic potentials

H. K. Pal; A. Shukla

doi:10.1016/j.cpc.2008.04.004

paper

A Fortran 90 program to solve the Hartree-Fock equations for interacting spin-1/2 Fermions confined in Harmonic potentials

arXiv:0807.3444 · doi:10.1016/j.cpc.2008.04.004

Abstract

A set of weakly interacting spin-1/2 Fermions, confined by a harmonic oscillator potential, and interacting with each other via a contact potential, is a model system which closely represents the physics of a dilute gas of two-component Fermionic atoms confined in a magneto-optic trap. In the present work, our aim is to present a Fortran 90 computer program which, using a basis set expansion technique, solves the Hartree-Fock (HF) equations for spin-1/2 $F er mi o n sco n f in e d b y a t h r ee - d im e n s i o na l ha r m o ni cosc i l l a t or p o t e n t ia l, an d in t er a c t in g w i t h e a c h o t h er v ia p ai r - w i se d e l t a - f u n c t i o n p o t e n t ia l s . A dd i t i o na l l y, t h e p r o g r am c ana l so a cco u n t f or t h ose anha r m o ni c p o t e n t ia l s w hi c h c anb ee x p r esse d a s a p o l y n o mia l in t h e p os i t i o n o p er a t or s$ x,y $, an d$ z$. Both the restricted-HF (RHF), and the unrestricted-HF (UHF) equations can be solved for a given number of Fermions, with either repulsive or attractive interactions among them. The option of UHF solutions for such systems also allows us to study possible magnetic properties of the physics of two-component confined atomic Fermi gases, with imbalanced populations. Using our code we also demonstrate that such a system exhibits shell structure, and follows Hund's rule.

19 pages and 2 figures

A Fortran 90 program to solve the Hartree-Fock equations for interacting spin-1/2 Fermions confined in Harmonic potentials

Abstract

References in corpus (13)