“…However, tracking back to a BEC described by the GP equation in the meanfield approximation, one can realize that the original Ermakov equation obtained in non-interacting case needs to be modified in the TF limit or in the case of a time-dependent interaction [21,26] . To remedy it, the variational approximation [30] (which is equivalent to moment method [31]), can be complemented by the concept of STA, for studying the dynamics of BECs [32][33][34], valid for the range from zero to small atomic interaction, with the implication on the quantum speed limits and quantum thermodynamics [35,36]. As a matter of fact, the accuracy of the variational approximation depends on the presumed ansatz in terms of nonlinearity [33,34].…”