The treatment of quantum thermodynamic systems beyond weak coupling and Markovian regimes is of increasing relevance, yet extremely challenging.  The evaluation of thermodynamic quantities in these regimes requires a nonperturbative knowledge of the bath dynamics, which in turn ytpically relies on heavy numerical simulations.  To tame these difficulties, considering thermal bosonic baths linearly coupled to the open system, one can consider expressions for heat, work, and average system-bath interaction energy that only involve the autocorrelation function of the bath and two-time expectation values of system operators.  In [F. Albarelli, B. Vacchini and A. Smirne, arXiv:2407.17886 (2024)] this perspective has been combined with the pseudomode approach, to numerically evaluate these relevant thermodynamic quantities.  It was shown in particular that this method allows for an efficient numerical evaluation of thermodynamic quantities in terms of one-time expectation values of the open system and the pseudomodes.