The study of a two-level molecular system with permanent dipole moments (PDM) immersed in a thermal bath and interacting with an electromagnetic field depends on the inclusion of the rotating wave approximation (RWA). A methodology based on cumulant expansions to obtain the average in the Fourier components associated with the coherences and populations, calculated by the optical stochastic Bloch equations, is employed. With these equations, we have modeled the optical responses for the four-wave mixing signal and we have demonstrated that inside the RWA the PDM cannot be measured.