In this paper experimental equilibrium data of the system CO2 - bergamot essential oil are reported at 323 K and 343 K, in the pressure range 7.8-13.1 MPa. Furthermore equilibrium data for the subsystem CO2 - linalyl acetate at 323 K are reported in the pressure range 7.0-10.3 MPa. The equilibrium data are predicted by means of two thermodynamic models based on Peng - Robinson and PC-SAFT equations of state. In particular the bergamot essential oil is described as a mixture of four components: limonene, linalool, linalyl acetate and beta-caryophyllene. In the models regression parameters are calculated from data of binary subsystems, CO2 - limonene, CO2 - linalool and CO2-beta-caryophyllene, and CO2 - linalyl acetate. Therefore both of the models are predictive with regard to the bergamot oil and show quite a good agreement (especially the Peng Robinson one) with respect to the experimental data. (C) 2016 Elsevier B.V. All rights reserved.