Molecular expressions for thermodynamic properties and derivatives of the Gibbs energy up to third order in the isobaric-isothermal (NpT) ensemble are systematically derived using the methodology developed by Lustig for the microcanonical and canonical ensembles [J. Chem. Phys. 100, 3048 (1994)10.1063/1.466446; Mol. Phys. 110, 3041 (2012)10.1080/00268976.2012.695032]. They are expressed by phase-space functions, which represent derivatives of the Gibbs energy with respect to temperature and pressure. Additionally, expressions for the phase-space functions for temperature-dependent potentials are provided, which, for example, are required when quantum corrections, e.g., Feynman-Hibbs corrections, are applied in classical simulations. The derived expressions are validated by Monte Carlo simulations for the simple Lennard-Jones model fluid at three selected state points. A unique result is that the phase-space functions contain only ensemble averages of combinations of powers of enthalpy and volume. Thus, the calculation of thermodynamic properties in the NpT ensemble does not require volume derivatives of the potential energy. This is particularly advantageous in Monte Carlo simulations when the interactions between molecules are described by empirical force fields or very accurate ab initio pair and nonadditive three-body potentials.