Magnetized neutron stars in X-ray binaries show a variety of spectacular timing and spectral properties. The strong NS magnetosphere dominates the accretion flow beyond the magnetospheric radius and channels matter to flow along the field lines on to the magnetic poles of the NS. X-ray emission from this region of the NS is modulated at the NS spin frequency and is observed as periodic X-ray pulsations.As part of my study of such accretion powered X-ray pulsars, I am interested in understanding the accretion flow-induced variability signatures like Quasi Periodic Oscillations (QPOs), Cyclotron Resonant Scattering Features (CRSF), flares and outbursts. Characterizing such properties in various samples of low mass and high mass X-ray binary systems will not only enable constraining the various theoretical emission models and flow instability models, but also allow for understanding the various possible evolutionary pathways for Pulsars in general.