Vibrational phase relaxation near gas-liquid and liquid-solid phase coexistence has been studied by molecular dynamics simulations of N-N stretch in N 2 . Experimentally observed pronounced insensitivity of phase relaxation from the triple point to beyond the boiling point is found to originate from a competition between density relaxation and resonant-energy transfer terms. The sharp rise in relaxation rate near the critical point (CP) can be attributed at least partly to the sharp rise in vibrationrotation coupling contribution. Substantial subquadratic quantum number dependence of overtone dephasing rate is found near the CP and in supercritical fluids. [S0031-9007(99)09318-7] PACS numbers: 61.25.Em, 64.70.Fx Vibrational phase relaxation (VPR) is a powerful tool to probe the interactions of a chemical bond with its surrounding medium [1]. Traditionally, VPR has been studied by the Raman scattering techniques, the most popular being the isotropic Raman line-shape studies [1]. The rate of VPR is proportional to the width of the usually obtained Lorentzian line shape. More recently, nonlinear optical techniques, such as fifth and seventh order responses from the vibrational mode, have been used to study this problem in greater depth [2(a)]. Many interesting new results have been obtained, such as the subquadratic quantum number dependence of overtone dephasing [2].The quantum number dependence of vibrational phase relaxation is usually explained by using the well-known Kubo-Oxtoby theory where the overtone vibrational line shape [I͑v͒ n0 ] for the ͑0, n͒ pair of vibrational levels is related to the normal-coordinate (Q) time correlation function for the nth overtone as [1]The normal-coordinate time correlation for the nth overtone ͓͗Q͑t͒Q͑0͒͘ n0 ] is in turn related to the instantaneous fluctuation in energy between the vibrational quantum levels n and 0, or frequency modulation ͓Dv n0 ͑t͔͒, bywhere v n0 is the vibrational frequency for the 0 to n transition. A cumulant expansion of Eq. (2) [3] leads to the expression of the normal-coordinate time correlation function in terms of the frequency-modulation time correlation function ͗Dv n0 ͑t͒Dv n0 ͑0͒͘ [1(a)]. The fluctuation in frequency gap [Dv n0 ͑t͒] derives contributions from several sources, including the atom-atom (AA) interactions, the resonant energy transfer (RT) between molecules and the vibration-rotation (VR) coupling [1]. While the atomatom interaction terms have been evaluated using analytical methods based on a hydrodynamic description [1(a)] and also at the molecular level using the mode coupling theory [4], the RT and VR contributions are notoriously difficult to evaluate analytically. The simulation work of Oxtoby et al.[5] amply demonstrated the complexity of the correlations and cross correlations of these two terms; even then, that study was restricted to only one thermodynamic state point (near the boiling point) of liquid N 2 . Perhaps as a consequence of the above, many interesting results have remained ill understood and many important que...