We investigate Rabi oscillation of an atom ensemble in Gaussian spatial distribution. By using the ultrafast laser interaction with the cold atomic rubidium vapor spatially confined in a magnetooptical trap, the oscillatory behavior of the atom excitation is probed as a function of the laser pulse power. Theoretical model calculation predicts that the oscillation peaks of the ensemble-atom Rabi flopping fall on the simple Rabi oscillation curve of a single atom and the experimental result shows good agreement with the prediction. We also test the the three-pulse composite interaction Rx(π/2)Ry(π)Rx(π/2) to develop a robust method to achieve a higher fidelity population inversion of the atom ensemble.PACS numbers: 32.80. Qk, 32.80.Wr, 42.65.Re Rabi oscillation is a fundamental concept in physics with a significant pedigree first discovered in the context of nuclear magnetic resonance (NMR) [1][2][3] and later extended to atomic physics and quantum optics [4,5]. In the presence of an oscillatory driving field E(t) = A(t) cos(ωt), a two-state quantum system undergoes a cyclic change of Bloch vector ρ manifested by the precessionabout an effective torque Ω = (−µA(t)/2 , 0, δ), where µ is the transition dipole moment between the two energy states, A(t) is the field envelope, and δ is the frequency detuning under the slowly-varying envelope approximation [4]. This generic feature of Rabi oscillation is universally found in a vast variety of material systems ranging from simple atoms and molecules [6][7][8][9][10] When a two-state atom interacts with a resonant (δ = 0) laser pulse, the dynamics of the excited state probability, which we may refer to as single-atom Rabi oscillation (SARO), is represented bywhere Θ o is the pulse area defined by Θ o = µA(t)dt/ . Since the pulse area is subject to both the pulse duration and the electric-field envelope, Rabi oscillations of an ultra-short time scale can be implemented by ultrafast optical interaction at a strong-enough laser intensity regime. However, the spatial extent of the laser beam over the laser-atom interaction region inevitably causes spatial average effect that often leads to vanishing of the oscillatory behavior. To overcome this problem, homogenizing the spatial profile of laser beams [22,23] and * Electronic address: jwahn@kaist.ac.kr adapting chirped laser interaction [24] have been considered. This paper aims quantitative analysis of spatially averaged Rabi oscillation. For this, we use the atom ensemble localized in a magneto-optical trap (MOT) [25] interacted with ultrafast laser pulses. As a theoretical model to investigate the spatially inhomogeneous interaction, we consider a Gaussian laser beam propagating along z direction. The pulse area in Eq. (2) is then represented in the polar coordinate system aswhere r = x 2 + y 2 , w(z) is the beam waist at z, w o = w(0) is the minimal beam waist, Θ o is the maximal pulse area, and Θ z = w o Θ o /w(z). When we assume the atom density profile in the MOT is also a Gaussian, i.e., ρ(r, z) = ρ o e −(r 2 +z 2 ...