As fourth-generation refrigerants with great development prospects, hydrofluoroolefins (HFOs) can be mixed with other refrigerants, such as carbon dioxide (CO2), to form refrigerant mixtures with low global warming potential (GWP) and zero ozone depleting potential (ODP) while retaining the advantages of each component. Refrigerants can work together to achieve complementary benefits. Combinations of CO2 and HFOs can strengthen the thermodynamic properties of CO2 while inhibiting the flammability of HFOs. At present, relatively few studies have been conducted on the CO2 + 3,3,3-trifluoropropene (R1243zf) mixture. Besides experimental approaches, molecular simulation has grown in importance as a way to determine thermodynamic and transport properties in recent years. In this study, the Gibbs Ensemble Monte Carlo (GEMC) method was used to simulate the vapor-liquid equilibrium (VLE) properties of the CO2 + R1243zf binary mixture in the temperature range of 273.15 to 313.15 K, in which the three-site rigid TraPPE force field and a fully flexible transferable all-atom force field were selected to describe CO2 and R1243zf, respectively. By comparing the GEMC simulation results with the experimental data, it was found that the average deviation of pressure is 2.33%, the average deviation of liquid phase molar fraction Δx is 0.0099, and the mean deviation of gaseous phase molar fraction Δy is 0.0204. The simulation results accord well with the experimental data and the fitting data of the correlation model in the literature, indicating that the molecular models used for CO2 and R1243zf can reliably predict the VLE properties. Finally, the critical parameters of the mixture at a temperature of 313.15 K were predicted, and the radial distribution functions (RDFs) of pure R1243zf and the mixture at 273.15 K and 3.5 MPa were calculated and analyzed by molecular dynamics (MD) simulation.