The neutron stars in low-mass X-ray binaries (LMXBs) are usually thought to be weakly magnetized objects accreting matter from their low-mass companions in the form of a disk. Albeit weak as compared to those in young neutron-star systems, the neutron-star magnetospheres in LMXBs can play an important role in determining the correlations between spectral and temporal properties. Parallel tracks appearing in the plane of kilohertz (kHz) quasi-periodic oscillation (QPO) frequency versus X-ray flux can be used as a tool to study the magnetosphere-disk interaction in neutron-star LMXBs. For dynamically important weak fields, the formation of a non-Keplerian magnetic boundary layer at the innermost disk truncated near the surface of the neutron star is highly likely. Such a boundary region may harbor oscillatory modes of frequencies in the kHz range. We generate parallel tracks using the boundary region model of kHz QPOs. We also present the direct application of our model to the reproduction of the observed parallel tracks of individual sources such as 4U 1608--52, 4U 1636--53, and Aql X-1. We reveal how the radial width of the boundary layer must vary in the long-term flux evolution of each source to regenerate the parallel tracks. The run of the radial width looks similar for different sources and can be fitted by a generic model function describing the average steady behavior of the boundary region in the long term. The parallel tracks then correspond to the possible quasi-steady states the source can occupy around the average trend.

We study the dependence of kHz quasi-periodic oscillation (QPO) frequency on accretion-related parameters in the ensemble of neutron star low-mass X-ray binaries. Based on the mass accretion rate, $\dot{M}$, and the magnetic field strength, $B$, on the surface of the neutron star, we find a correlation between the lower kHz QPO frequency and $\dot{M}/B^{2}$. The correlation holds in the current ensemble of Z and atoll sources and therefore can explain the lack of correlation between the kHz QPO frequency and X-ray luminosity in the same ensemble. The average run of lower kHz QPO frequencies throughout the correlation can be described by a power-law fit to source data. The simple power-law, however, cannot describe the frequency distribution in an individual source. The model function fit to frequency data, on the other hand, can account for the observed distribution of lower kHz QPO frequencies in the case of individual sources as well as the ensemble of sources. The model function depends on the basic length scales such as the magnetospheric radius and the radial width of the boundary region, both of which are expected to vary with $\dot{M}$ to determine the QPO frequencies. In addition to modifying the length scales and hence the QPO frequencies, the variation in $\dot{M}$, being sufficiently large, may also lead to distinct accretion regimes, which would be characterized by Z and atoll phases.