Abstract
We develop and employ a Thomas-Fermi-Dirac-Weizsäcker (TFD-λW) statistical model to investigate diatomic molecules and infinite molecular chains in strong magnetic fields. Such material may form on the surface of a neutron star where field strengths can reach B ∼ 1012 G. The standard magnetic Thomas-Fermi-Dirac kinetic, potential, and exchange energy functionals are supplemented by a gradient correction to the kinetic energy. This correction leads to a model of electronic structure determined by a nonlinear system of two coupled partial-differential equations with an eigenvalue. The numerical method used for solving this system in two spatial dimensions is detailed. In contrast to simpler Thomas-Fermi and Thomas-Fermi-Dirac treatments, molecules and chains can be bound in the TFD-λW model due to the gradient term. We present numerical solutions for a wide range of magnetic field strengths and elements to demonstrate the robustness, as well as the limitations, of the statistical approach. Our calculations qualitatively reproduce many of the results of detailed quantum mechanical treatments. For example, the fractional binding energy is greatest for low atomic numbers and for strong magnetic fields.
Original language | English (US) |
---|---|
Pages (from-to) | 233-241 |
Number of pages | 9 |
Journal | Astrophysical Journal |
Volume | 382 |
Issue number | 1 |
DOIs | |
State | Published - Nov 20 1991 |
Externally published | Yes |
Keywords
- Molecular processes
- Stars: magnetic
- Stars: neutron
ASJC Scopus subject areas
- Astronomy and Astrophysics
- Space and Planetary Science