In present work, we will propose the spectroscopic model for the analysis of Raman gain spectrum based on oscillator theory and its application to the design of optical fiber Raman amplifiers with the multiwavelength pumping scheme. The oscillatory lineshape functions can simply describe the wavelength dependence of Raman gain in optical fibers and are very useful for the estimation of the gain bandwidth, Raman lasing, noise performance, and amplification processes in Raman amplifiers. The concept of an actual band based on the lineshape function is useful for the design of fiber Raman amplifiers with multiple wavelengths pumping. Proposed modeling allows us to analyze fiber Raman amplifier with combined multiwavelength pumping source for the extension of amplification banrwidth to L-band, which has the broad bandwidth over 80 nm and low gain ripple less than 0.5 dB.Fiber Raman amplifiers (FRA) are being deployed in almost every new long-haul fiberoptic transmission systems [1]. Despite the simple architecture of a Raman amplifier, many factors in the design of amplifiers such as pump-to-pump power transfer, signal-to-signal power transfer, pump depletion (saturation), double Rayleigh scattering (multipath interference), and amplifier spontaneous noise should be considered [2]. The simplistic design of the multiple pumping Raman amplifier results to a very poor amplifier when the evenly spaced eight pumps with equal pump power of 120 mW gives the FRA with the gain variation more than 10 dB [2]. Numerical optimization methods for the optimal determination of the power and wavelength of each pump source based on the coupled equations such as simulated annealing algorithm [3], single-layer feedforward neural networks [4], and so on [5] have been investigated. The theoretical analysis of the Raman gain coefficient with one pumping wavelength in optical fiber and comprehensive analysis of its scaling with wavelength, modal overlap, and material composition are also reported [6].In present work, we will propose the spectroscopic model for the analysis of Raman gain spectrum based on oscillator theory and its application to the design of wideband optical FRA.Pumping saturation is not negligible in distributed optical fiber Raman amplifier with multiple wavelengths pumping. The pumps depletion caused by interaction between the pumps and the amplified signals is described by a well-known set of coupled equations [2][3][4][5]8]. When the frequency dependence of the Raman gain gR(v) is modeled as a frequency independent constant or simple single peak, for example as a triangle, the numerical methods may give the unsatisfactory results because the solutions of coupled equation are very sensitive to absolute values ofequation coefficients.According to the semi-classical derivation [6], spontaneous Raman scattering (SRS) can be described in terms of classical optical fields and the molecular system can be treated as a forced harmonic oscillator. This model reveals the connection between the differential polarizability (aa, /aq...