Silicon photonic networks are revolutionizing computing systems by improving the energy efficiency, bandwidth, and latency of data movements. Optical modulators, such as microresonators and Mach-Zehnder Interferometers (MZIs), are the basic building blocks of silicon photonic networks. This work proposes a SPICE-compatible electro-optical co-simulation model, BOSIM, to systematically study optical modulators using PN, PIN, and metal-insulator-silicon (MIS) capacitor device technologies. BOSIM holistically models both transient and steady state properties, such as switching speed, power, transmission spectrum, area and carrier distribution. BOSIM is validated by the measured data from eight research groups and companies. Compared to microresonators, BOSIM shows MZIs are fast, with a high extinction ratio and large bandwidth but in the sacrifice of loss, energy, and area. Using a PIN diode over a PN diode can save area, but retain the loss and energy, while an MIS capacitor has the shock response of carrier distribution in a narrow range and is marginalized gradually. For instance, an MZI can achieve a 2.5X bit rate, 6.06X extinction ratio, 71.04X 3dB bandwidth, but costs at least 1.93X passing loss, 1.46X energy consumption and 16.67X area, compared with microresonator.