This work presents
a new method how to obtain approximate analytical
solutions for the Euler equation for second-row homonuclear dimers.
In contrast to the well-known Kohn–Sham method where a system
of N nonlinear coupled differential equations must be solved iteratively,
orbital-free density functional theory allows to access the minimizing
electron density directly via the Euler equation. For simplified models,
here, an atom-centered monopole expansion with one free parameter,
solutions of the electron density can be obtained analytically by
solving the Euler equation at the bond critical point. The procedure
is exemplarily carried out for N2, C2, and B2, yielding bound molecules with an internuclear distance of
2.01, 2.43, and 3.07 bohr, respectively.