The enzymatic reaction−distillation coupling process is a promising technology to realize the chiral resolution of (R/S)-1-phenylethanol with ethyl butyrate and ethanol as acyl donor and byproduct, respectively. In order to design such a coupling process, the vapor−liquid equilibrium (VLE) data are necessary. However, these data were incomplete in literatures. In this work, the isobaric VLE experiments were conducted at 101.3 kPa for the binary systems of 1-phenylethanol + ethyl butyrate, 1-phenylethanol + ethanol, 1-phenylethyl butyrate + ethyl butyrate, 1phenylethyl butyrate + ethanol, and 1-phenylethyl butyrate +1-pheny ethanol. Additionally, the saturated vapor pressure of 1-phenylethyl butyrate was measured in the temperature range 420−500 K and correlated with the Antoine equation. Both the Non-Random Two Liquid (NRTL) and Wilson models were used for the correlation of experimental VLE data. Both models showed satisfactory accuracy, although the Wilson model performed slightly better than the NRTL model. This work would provide important thermodynamic information for the design of enzymatic reaction−distillation coupling process to realize the chiral resolution of (R/S)-1-phenylethanol.