Optimization of pulsed gas tungsten arc welding (pulsed GTAW) process parameters was carried out to obtain optimum weld bead geometry with full penetration in welding of stainless steel (304L) sheets of 3 mm thickness. Autogenuous welding with square butt joint was employed. Design of experiments based on central composite rotatable design was employed for the development of a mathematical model correlating the important controllable pulsed GTAW process parameters like pulse current (Ip), pulse current duration (Tp), and welding speed (S) with weld bead parameters such as penetration, bead width (W), aspect ratio (AR), and weld bead area of the weld. The developed models were checked for adequacy based on ANOVA analysis and accuracy of prediction by conducting a confirmation test. Weld bead parameters predicted by the models were found to confirm observed values with high accuracy. Using these models, the main and interaction effects of pulsed GTAW process parameters on weld bead parameters were studied and discussed. Optimization of pulsed GTAW process parameters was carried out to obtain optimum bead geometry using the developed models. A quasi-Newton numerical optimization technique was used to solve the optimization problem and the results of the optimization are presented.