We present a newly developed Monte Carlo scheme to predict bulk surfactant concentrations and surface tensions at the air water interface for various surfactant interfacial coverages. Since the concentration regimes of these systems of interest are typically very dilute (<< 10(-5) mol. frac.), Monte Carlo simulations with the use of insertion/deletion moves can provide the ability to overcome finite system size limitations that often prohibit the use of modem molecular simulation techniques. In performing these simulations, we use the discrete fractional component Monte Carlo (DFCMC) method in the Gibbs ensemble framework, which allows us to separate the bulk and air water interface into two separate boxes and efficiently swap tetraethylene glycol surfactants C10E4 between boxes. Combining this move with preferential translations, volume biased insertions, and Wang Landau biasing vastly enhances sampling and helps overcome the classical "insertion problem", often encountered In non-lattice Monte Carlo simulations. We demonstrate that this methodology is both consistent with the original molecular thermodynamic theory (MTT) of Blankschtein and co-workers, as well as their recently modified theory (MD/MTT), which incorporates the results of surfactant infinite dilution transfer free energies and surface tension calculations obtained from molecular dynamics simulations.