Phosphorus-related species are not known to be as omnipresent in space as hydrogen, carbon, nitrogen, oxygen, and sulfur-bearing species. Astronomers spotted very few P-bearing molecules in the interstellar medium and circumstellar envelopes. Limited discovery of the P-bearing species imposes severe constraints in modeling the P-chemistry. In this paper, we carry out extensive chemical models to follow the fate of P-bearing species in diffuse clouds, photon-dominated or photodissociation regions (PDRs), and hot cores/corinos. We notice a curious correlation between the abundances of PO and PN and atomic nitrogen. Since N atoms are more abundant in diffuse clouds and PDRs than in the hot core/corino region, PO/PN reflects <1 in diffuse clouds, MUCH LESS-THAN1 in PDRs, and >1 in the late warm-up evolutionary stage of the hot core/corino regions. During the end of the post-warm-up stage, we obtain PO/PN > 1 for hot core and <1 for its low-mass analog. We employ a radiative transfer model to investigate the transitions of some of the P-bearing species in diffuse cloud and hot core regions and estimate the line profiles. Our study estimates the required integration time to observe these transitions with ground-based and space-based telescopes. We also carry out quantum chemical computation of the infrared features of PH3, along with various impurities. We notice that SO2 overlaps with the PH3 bending-scissoring modes around similar to 1000-1100 cm(-1). We also find that the presence of CO2 can strongly influence the intensity of the stretching modes around similar to 2400 cm(-1) of PH3.