text: The output beam is collected by a detector and its intensity  is calculated using
  where  is the visible frequency,  is the IR frequency and  is the SFG frequency. The constant of proportionality varies across literature, many of them including the product of the square of the output frequency,  and the squared secant of the reflection angle, .  Other factors include index of refractions for the three beams.
  The second order susceptibility has two contributions
  where  is the non-resonating contribution and  is the resonating contribution.  The non-resonating contribution is assumed to be from electronic responses. Although this contribution has often been considered to be constant over the spectrum, because it is generated simultaneously with the resonant response, the two responses must compete for intensity. This competition shapes the nonresonant contribution in the presence of resonant features by resonant attenuation. Because it is not currently known how to adequately correct for nonresonant interferences, it is very important to experimentally isolate the resonant contributions from any nonresonant interference, often done using the technique of nonresonant suppression.
  The resonating contribution is from the vibrational modes and shows changes in resonance.  It can be expressed as a sum of a series of Lorentz oscillators
  where  is the strength or amplitude,  is the resonant frequency,  is the damping or linewidth coefficient (FWHM), and each  indexes the normal (resonant vibrational) mode. The amplitude is a product of , the induced dipole moment, and , the polarizability.  Together, this indicates that the transition must be both IR and Raman active.
  The above equations can be combined to form
  which is used to model the SFG output over a range of wavenumbers.  When the SFG system scans over a vibrational mode of the surface molecule, the output intensity is resonantly enhanced.  In a graphical analysis of the output intensity versus wavenumber, this is represented by Lorentzian peaks.  Depending on the system, inhomogeneous broadening and interference between peaks may occur.  The Lorentz profile can be convoluted with a Gaussian intensity distribution to better fit the intensity distribution.