A definition of nonequilibrium free energy $mathcal{F}_{textsc{s}}$ is proposed for dynamical Gaussian quantum open systems strongly coupled to a heat bath and a formal derivation is provided by way of the generating functional in terms of the coarse-grained effective action and the influence action. For Gaussian open quantum systems exemplified by the quantum Brownian motion model studied here, a time-varying effective temperature can be introduced in a natural way, and with it, the nonequilibrium free energy $mathcal{F}_{textsc{s}}$, von Neumann entropy $mathcal{S}_{vN}$ and internal energy $mathcal{U}_{textsc{s}}$ of the reduced system ($S$) can be defined accordingly. In contrast to the nonequilibrium free energy found in the literature which references the bath temperature, the nonequilibrium thermodynamic functions we find here obey the familiar relation $mathcal{F}_{textsc{s}}(t)=mathcal{U}_{textsc{s}}(t)- T_{textsc{eff}} (t),mathcal{S}_{vN}(t)$ {it at any and all moments of time} in the systems fully nonequilibrium evolution history. After the system equilibrates they coincide, in the weak coupling limit, with their counterparts in conventional equilibrium thermodynamics. Since the effective temperature captures both the state of the system and its interaction with the bath, upon the systems equilibration, it approaches a value slightly higher than the initial bath temperature. Notably, it remains nonzero for a zero-temperature bath, signaling the existence of system-bath entanglement. Reasonably, at high bath temperatures and under ultra-weak couplings, it becomes indistinguishable from the bath temperature. The nonequilibrium thermodynamic functions and relations discovered here for dynamical Gaussian quantum systems should open up useful pathways toward establishing meaningful theories of nonequilibrium quantum thermodynamics.