Abstract. Phase mixing of relativistically intense longitudinal wave packets in a cold homogeneous unmagnetized plasma has been studied analytically and numerically using Dawson Sheet Model. A general expression for phase mixing time (ω p t mix ) as a function of amplitude of the wave packet(δ) and width of the spectrum(∆k/k) has been derived. It is found that phase mixing time crucially depends on the relative magnitude of amplitude "δ" and the spectral width "∆k/k". For ∆k/k ≤ 2ω 2 p δ 2 /c 2 k 2 , ω p t mix scales with δ as ∼ 1/δ 5 , whereas for ∆k/k > 2ω 2 p δ 2 /c 2 k 2 , ω p t mix scales with δ as ∼ 1/δ 3 , where ω p is the non-relativistic plasma frequency and c is the speed of light in vacuum. We have also verified the above theoretical scalings using numerical simulations based on Dawson Sheet Model.