Accurate characterization of the environmental effects on a quantum system remains a fundamental challenge in the theory of open quantum systems. In this talk, I will introduce the purified pseudomode approach developed by us recently. This method allows for efficient modeling and numerically exact simulation of general linear-Gaussian baths. Extensions of this method to model bath input-output and nonlinear system-bath interactions will also be discussed.