In energy decomposition analysis of Kohn-Sham density functional theory calculations, the so-called frozen (or pre-polarization) interactionenergy contains contributions from permanent electrostatics, dispersion, and Pauli repulsion. The standard classical approach to separate them suffers from several well-known limitations. We introduce an alternative scheme that employs valid antisymmetric electronic wavefunctions throughout and is based on the identification of individual fragment contributions to the initial supersystem wavefunction as determined by an energetic optimality criterion. The density deformations identified with individual fragments upon formation of the initial supersystem wavefunction are analyzed along with the distance dependence of the new and classical terms for test cases that include the neon dimer, ammonia borane, water-Na+, water-Cl−, and the naphthalene dimer.