While the relatively complex chemical composition of Prussian blue analogues (PBAs) allows a high degree of flexibility in tuning the catalytic performances toward CO2-based copolymerization, little is known regarding the influence of the crystal phase. Herein, we investigate the impact of the crystal phase tuning of the representative Zn−Co-based PBA on CO2−epoxide copolymerization. The utilization of oxygen-containing tert‑butanol molecules (tBuOH) replaces water molecules from the Zn active center within the PBA-based nanosheets, inducing the crystal transformation from conventionally cubic to rhombohedral. The rhombohedral phase can be formed as clusters dispersed within the amorphous matrix, exhibiting the most active and selective CO2−epoxide copolymerization performance with a turnover number of 31,250 g/g Zn compared to that of long-range ordered cubic and rhombohedral counterparts. Such an improvement is ascribed to the abundant exposed Zn-O active centers enabled by the replacement of tBuOH molecules and rhombohedral coordination. Density functional theory calculation indicates different average adsorption energies for cyclohexene oxide (CHO) and CO2 (-5.45 vs. -0.02 eV) on the (001) surface, implying the initialization of copolymerization from the adsorption of CHO onto the Zn-O active centers. This work provides a feasible approach to enable on-demand manipulation of catalytic performances of emerging PBAs.