There is no doubt that electric fields of a specific frequency and intensity could excite certain vibrational modes of a macromolecule, which alters its mode coupling and conformation. Motivated by recent experiments and theories, we study the mode coupling between the Fe-CO mode and CO-stretch mode and vibration energy transfer among the active site and proteins in carboxyhemoglobin (HbCO) under different electric fields using the quasistatic two-dimensional infrared spectra. This study uses iron-porphyrin-imidazole-CO and two distal histidines in HbCO as the subsystem. The potential energy and dipole moment surfaces of the subsystem are calculated using an all-electron ab initio (B3LYP-D3(BJ)) method with the basis set Lanl2dz for the Fe atom and 6-31G(d,p) for C, H, O, and N atoms. Although the subsystem is reduced dimensionally, the anharmonic frequency and anharmonicity of the CO-stretch mode show excellent agreement with experimental values. We use the revealing noncovalent interaction method to confirm the hydrogen bond between the H-epsilon atom of the His63 and the CO molecule. Our study confirms that the mode coupling between the Fe-CO mode and CO-stretch mode does not exist when the subsystem is free of electric field perturbation, which is coupled when the electric field is -0.5142 V/nm. In addition, with the increases of distance between the active site and the His92, there is no vibrational energy transfer between them when the electric field is 1.028 V/nm. We believe that our work could provide new ideas for increasing the dissociation efficiency of the Fe-CO bond and theoretical references for experimental research.