We propose a general variational fermionic many-body wavefunction that generates an effective Hamiltonian in quadratic form which can then be exactly solved. The theory can be constructed within density functional theory framework and a self-consistent scheme is proposed for solving the exact density functional theory. We apply the theory to structurally disordered system and a symmetric and asymmetric Hubbard dimer and corresponding lattice models and the the single fermion excitation spectra show a persistent gap due to the fermionic entanglement induced pairing condensate. For disordered system, density of state at the edge of the gap diverges in the thermodynamic limit, suggesting a topologically ordered phase and a sharp resonance is predicted as the gap is not dependent on the temperature of the system. For the symmetric Hubbard model, the gap for both half filling and doped case suggests quantum phase transition between the AFM and SC is a continuous phase transition.