Properties of constraint-based single-point approximate kinetic energy functionals

We present an analysis and extension of our constraint-based approach to orbital-free (OF) kinetic-energy (KE) density functionals intended for the calculation of quantum-mechanical forces in multiscale molecular-dynamics simulations. Suitability for realistic system simulations requires that the OF-KE functional yield accurate forces on the nuclei yet be computationally simple. We therefore require that the functionals be based on density-functional theory constraints, be local, be dependent at most upon a small number of parameters fitted to a training set of limited size, and be applicable beyond the scope of the training set. Our previous modified-conjoint generalized-gradient-type functionals were constrained to producing a positive-definite Pauli potential. Though distinctly better than several published generalized-gradient-approximation-type functionals in that they gave semiquantitative agreement with Born-Oppenheimer forces from full Kohn-Sham results, those modified-conjoint functionals suffer from unphysical singularities at the nuclei. Here we show how to remove such singularities by introducing higher-order density derivatives and analyze the consequences. We give a simple illustration of such a functional and a few tests of it.