We continue exploring our understanding of the collective properties of X-ray binaries in the stellar fields (i.e., outside globular clusters) of normal galaxies, introduced in Paper I of this series, where we considered high-mass X-ray binaries (HMXBs). In this paper (Paper II of the series) and the companion paper (Paper III of the series), we consider low-mass X-ray binaries (LMXBs), whose evolutionary scenario is very different from that of HMXBs. In this paper, we consider the evolution of primordial binaries upto the stage where the neutron star just formed in the supernova explosion of the primary is in a binary with its low-mass unevolved companion, and this binary has circularized tidally, producing what we call a pre-low-mass X-ray binary (pre-LMXB). We study the constraints on the formation of such pre-LMXBs in detail (since these are low probability events), and calculate their collective properties and formation rate. To this end, we first consider the changes in the binary parameters in the various steps involved, viz., the common-envelope (CE) phase, the supernova, and the tidal evolution. This naturally leads to a clarification of the constraints. We then describe our calculation of the evolution of the distributions of primordial binary parameters into those of the pre-LMXB parameters, following the standard evolutionary scenario for individual binaries. We display the latter as both bivariate and monovariate distributions, discuss their essential properties, and indicate the influence of some essential factors on these. Finally, we calculate the formation rate of pre-LMXBs. The results of this paper are used in the next one (Paper III) to compute the expected X-ray luminosity function (XLF) of LMXBs, which is compared with observation.

Continuing our exploration of the collective properties of low-mass X-ray binaries (LMXBs) in the stellar fields of normal galaxies, we compute in this paper the expected X-ray luminosity function (XLF) of LMXBs, starting from the results of the previous paper in the series (Paper II). We treat separately two classes of LMXB evolution, the first being close systems whose initial orbital periods are below the bifurcation period, wherein the companion is on the main sequence when Roche-lobe contact occurs, the subsequent evolution is driven by angular-momentum loss from the system, and the second being wider systems whose initial orbital periods are above the bifurcation period, wherein the companion is on the giant branch when Roche-lobe contact occurs, and the evolution of these systems is driven by the nuclear evolution of the companion. We obtain model luminosity profiles L(t) for individual LMXBs of both classes. We then compute the LMXB XLF by folding in the inputs for the pre-LMXB collective properties and formation rates with the above luminosity profiles. Because of the long timescale on which LMXBs evolve, one needs to keep track of the evolution of the star-formation rate (SFR) on the same timescale, and we use star-formation histories given by canonical models. We compare the observed LMXB XLF with our computed one, keeping in mind that we have included only neutron-star systems in this work, so that there would be an unaccounted-for population of black-hole binaries at the high-luminosity end. We show that a qualitative similarity already exists between the two, and we discuss the role of the giant fraction, i.e., that fraction of all LMXBs which harbors a low-mass giant companion, on the shape of the XLF at the high-luminosity end.

We introduce a method for obtaining the X-ray luminosity function (XLF) and the binary-period distribution of populations of high-mass X-ray binaries (HMXBs) in the stellar fields (i.e. outside globular clusters) of normal galaxies. We start from standard distributions of the parameters of those primordial binaries which are the progenitors of HMXBs, and follow the transformation of these ditributions with the aid of a Jacobian formalism as the former evolve into the latter through the processes of the first mass transfer and the supernova (SN) that follows. We discuss the distributions of the post-SN binaries and the HMXBs. We show that our calculated model XLF has a differential slope $\approx -1.6$ with a flattening at low luminosities, in excellent agreement with observations. The calculated binary-period distribution, which basically has a slightly sloping plateau-like character at intermediate periods, with a rise to this plateau at shorter periods and fall-off from it at longer periods, is in agreement with the observed distribution within observational uncertainties. We discuss the physical origin of these distributions. We demonstrate that, while the effects of both (a) the distribution of the properties of the massive companion in the HMXBs, and (b) the primordial orbital distribution and the SN dynamics are important, the former appear to be dominant in determining the XLF, and the latter in determining the HMXB binary-period distribution. We discuss the possible roles of stellar-mass black holes and ultra-luminous X-ray sources (ULX) in the observed "universal" XLF of HMXBs.

High mass X-ray binary luminosity function (XLF) is an important tool for studying binary evolution processes and also the mass loss and consequent evolution in massive stars. We calculate the XLF for neutron star binaries using the standard scenario for formation and evolution of these systems. A one to one relation between primordial binary parameters and the HMXB parameters is established. The probability density function is then transformed using the standard Jacobian formalism. It is shown that the model successfully explains some basic properties of the observed XLF.

Discovery of the 6.7-hour periodicity in the X-ray source 1E 161348-5055 in RCW 103 has led to investigations of the nature of this periodicity. We explore a model for 1E 161348-5055, wherein a fast-spinning neutron star with a magnetic field $\sim 10^{12}$ G in a young pre-Low-Mass X-ray Binary (pre-LMXB) with an eccentric orbit of period 6.7 hr operates in the "propeller" phase. The 6.7-hr light curve of 1E 161348-5055 can be quantitatively accounted by a model of orbitally-modulated mass transfer through a viscous accretion disk and subsequent propeller emission (both Illarionov-Sunyaev type and Romanova-Lovelace et al type), and spectral and other properties are also in agreement. Formation and evolution of model systems are shown to be in accordance both with standard theories.