“…It has been shown that the constrained free energy
as a functional of given density fields n ( r , t ) = ( n 1 ( r , t ), n 2 ( r , t )) T can be expressed by considering fluctuating potential fields ϕ ( r , t ) = ( ϕ 1 ( r , t ), ϕ 2 ( r , t )) T , which are conjugate to n ( r , t ), in addition to an adjusted potential field ϕ dft l ( r , t ) similar to that of the equilibrium DFT. 37,38 Extending the previous result 24,28 to the expression for two-component systems (see Appendix B for details), we have
with the following constraint imposed by the canonical ensemble:
where the total number of either anions or cations is equally N . The free-energy functional F [ n , ϕ ] in the exponent of eqn (A3) is defined using the grand potential of the primitive model with an imaginary external field i ϕ ( r ) applied, and can be divided into two parts (see Appendix B for details): F [ n , ϕ ] = F [ n ,0] + Δ F [ n , ϕ ].The free-energy functional F [ n ,0] in the absence of fluctuating potential reduces to the intrinsic Helmholtz free energy, a key thermodynamic quantity in the equilibrium DFT.…”