A novel nonparaxial technique based on Pade approximant for modeling ultra short optical pulses has been discussed and verified. The characterization of the technique shows excellent performance in terms of accuracy, efficiency, and stability where the traditional paraxial time‐domain beam propagation method failed to achieve. The new method is very useful and efficient for modeling long optical pulse interaction when compared to the classical finite‐difference time‐domain. © 2007 Wiley Periodicals, Inc. Microwave Opt Technol Lett 49: 2481–2486, 2007; Published online in Wiley InterScience (www.interscience.wiley.com) DOI 10.1002/mop.22796