Some time back, a finite-temperature Schrodinger equation was obtained to describe the pairing of an electron and a proton in a medium of such particles at finite temperature  . This work was a follow-up of an earlier paper where the Coulomb potential was temperature-generalized for the first time  . In these papers the approach followed was somewhat intuitive, leaving open the question of a derivation from first principles. Hence, beginning with the Bethe-Salpeter equation in vacuum, such a derivation was given in  . In an approximation scheme, which should work at low particle densities and high temperatures, this equation gives a temperature- dependent bound state spectrum for the electron-proton pair. In the temperature range around 106 K, the deepest of the states in the spectrum have binding energies in the keV range and can withstand the background thermal agitation of the medium. The transitions from the short-lived excited states to the deepest ones in this spectrum lead to spectral lines in the soft X-ray region. An application of the approach to the flaring regions of the Sun therefore leads to the prediction of such lines in the flare spectra. In the present note, we report on this matter: the calculated lines at a certain temperature from three Balmer-like series are essentially all seen in the flare data. We also point out that, since many of these lines are identifiable with the lines which in the conventional approach are presumed to originate from the low first ionization potential (FIP) elements, our approach seems to provide an explanation of the FIP effect which has been extensively discussed in the literature, e.g., in  and  . The details are given in the following.
2. The Finite-Temperature Schrodinger Equation (FTSE) and Its Solutions
The equation for the bound state at temperature T, or the FTSE, is given by
where W is the binding energy,
with ma and mb as the electron and the proton mass, respectively. The function Q(W, p) is given by
and k is the Boltzmann constant. Using the approximation
we can write
Equation (7) can be solved in this approximation to give 
where l is the angular momentum quantum number, m the azimuthal quantum number and are given in terms of Gegenbauer polynomials:
the coefficients ak in Equation (8) obey a difference equation given by
where we have put
The parameters ci are given by
Equation (10) is a second order difference equation and will in general have two solutions; the dominant and the dominated. The eigenvalues of the equation are those values of (or W) for which the dominant solution vanishes (for a review see   ), and they can be computed easily by the Hill determinant method; see  -  and  . The generic computer program for such computations is already given in  , and can be adapted to the present situation without any difficulty. On substituting the numerical values for ma (electron mass), mb (proton mass) and α (the fine structure constant), we can use the program to solve for W at any given value of T. What emerges is a discrete spectrum for W: for each l, we get an increasing set of W values which can be labeled by a serial number n (which plays the role of the principal quantum number in the case of the hydrogen atom). The results are illustrated in Table 1. The typical dimensions of the Hill determinant for stability of levels up to the third decimal place (in eV) are k = 400 for n = 1, l = 0; k = 16,000 for n = 10, l = 0, k = 80,000 for n = 20, l = 0; and so on.
It is now important to make sure that the solutions obtained are indeed consistent with the approximation
given in Equation (5). This is done by verifying that the expectation values of over the desired solution corresponding to any Wn,l are quite small, say < 0.01. If not so, the solution must be discarded. It may be noted that since ma is much smaller than mb, it is sufficient to check that is quite
small. To do this, we have to first calculate.. The details are given in Appendix A. We quote here only the result:
Note that M in the summations S1 and S2 denotes the highest k value for which ak ≠ 0. is the magnitude of the binding energy at serial number n and angular momentum l, corresponding to which is being calculated. Note further that is already calculated from Equation (10) by the Hill determinant method, which also determines the value of M. This method, however, does not require an explicit knowledge of the coefficients ak. The latter are therefore yet to be determined. This can be done conveniently by using the backward Miller algorithm   which is reviewed in Appendix A, where we also record the steps that need to be taken to calculate ak and whence can be calculated. We then find that the smallness criterion is satisfied
Table 1. Energy eigenvalues of Equation (10) at temperature T = 4.26 × 106 K for 1 ≤ n ≤ 30 and l = 0, 1, and 2 obtained via the Hill determinant method. Note that energy level (n = 1, l = 0) will not be physically realizable at this temperature (see the text for details). Typical dimensions of the Hill determinant for stability up to the third place of decimal are: k = 400 (n = 1, l = 0); k = 16,000 for (n = 10, l = 0).
for any given Wn,l only if the temperature T is above a certain value Tn, and furthermore, T1 > T2 > T3…etc. The relevant results are displayed in Table 2. We see from this table that the solution corresponding to W1,0 is acceptable only for T > T1 ≈ 107 K, that corresponding to W20 or W21 only for, and so on. Clearly, all the solutions corresponding to Wn,l for n > 2 are now acceptable if T > T2. We must also keep in mind that a bound state with binding energy W will survive in a medium at temperature T only if is appreciably greater than (3kT/2), the latter being the average kinetic energy of thermal motion per particle in the medium. Thus, for instance, at a temperature T = 4.26 × 106 K, the deepest possible states that can form in our spectrum are at the level n = 2, i.e., W2,0 and W2,1. We note that for these levels. So the states at these levels will not only form, but also survive in the hot medium. Of course the states with n > 2 will increasingly have and will therefore behave as very short-lived excited states. The transitions from these n > 2 excited states to the n = 2 states should lead, in view of selection rule, to the following three Balmer- like series of spectral lines:
The wavelength λ of the spectral line corresponding to any transition involving the energy difference in eV is given by
Substituting for any allowed from the above series, and using numbers from Table 1, we can check that the resulting wavelengths fall in the soft X-ray region.
3. Application to Solar Flares
Let us now apply our approach to Solar flares. The latter are appropriate for such an application as they occur in a medium with low particle densities (≈1012/cm3) at temperatures around a few million degrees Kelvin and are known to emit soft X-ray lines. Furthermore, since the flaring phenomenon is a prolonged affair, we may expect that between the initial (growing) phase and the final (decaying) phase, there should be a period over which the flare burns at a reasonably constant temperature (with fluctuations of, say, not more than). If we knew this temperature, we could calculate the spectral lines of the kind given by Equation (13) and then look for them in the soft X-ray line data from flares. Of course this temperature is not quite known to us. So what we could do is to calculate our spectral lines at different temperatures between 106 K and 107 K, to see if there is any temperature in this range at which our calculated lines are reproduced in the flare data. As will be seen in the following, we find that there indeed is such a temperature given by T = 4.26 × 106 K. The binding energy spectrum at this temperature is given in Table 1. From here, the wavelengths of spectral lines corresponding to series A, B, and C, given by Equation (13), can easily be calculated. They are given in Table 3. The notation () means that the transition is from the level to. These transitions are referred to in columns 2 and 6 and the corresponding calculated wavelengths are given in columns 3 and 7.
Table 2. Temperatures for decreasing values of for several energy levels calculated via Miller’s algorithm.
Table 3. Calculated wavelength λ(Ȧ) at temperature T = 4.26 × 106 K for the emission lines comprising the generalized Balmer series A, B, and C which characterize our pair spectrum. Except for a very small difference in wavelength for the first set of lines at (3.0 → 2.1) and (3.2 → 2.1), the series A and C are degenerate. Therefore for the rest of the transitions in these series, only those from the former are given in the table. The observed wavelengths (good to within) are from the data on Solar flares by Phillips et al.  ; those marked with (+) are from the data of McKenzie et al.  .
We now turn to the relevant experimental information. The data between 5 and 20 Ȧ were obtained by Phillips et al.  from a flare on August 25, 1982 and from another flare on November 5, 1980. The work of Acton et al.  covers the range between 11 and 94 Ȧ. These observations were made on an M-class flare on July 13, 1982. They covered a much larger range of wavelengths than covered by Phillips et al., but missed out bunches of lines, especially between 12.82 and 13.45 Ȧ. McKenzie et al.  also provide useful data between 8 and 22 Ȧ. In the range common to all these three groups, the data of  are the most exhaustive. As to the accuracy of the wavelengths quoted by these authors, Acton et al. estimate it to be not better than. The other authors do not give direct estimates of their own, but their discussions suggest that it would not be safe to violate the limits set by Acton et al. Keeping this in mind, we now compare the emission lines from these data with those calculated from our pair spectrum at T = 4.26 × 106 K. Any line from the data which is within of a calculated line is identified with the latter. We may mention here that a change of up to does not make any change in our calculated wavelengths up to the third decimal place.
The experimental results in Table 3 (columns 4 and 8) are taken from Phillips et al.  , except in the range not covered by them, i.e., for λ(3.1 → 2.0), λ(3.0 → 2.1) and λ(3.2 → 2.1), which are from McKenzie et al. (the corresponding numbers from Acton et al. are 21.80 Ȧ and 21.60 Ȧ, not different from those of McKenzie et al. within the accuracy limit of). An exceptional situation exists in the case of λ(4.1 → 2.0): the line close to its calculated value, 16.947 Ȧ, is not seen by Phillips et al., though it is within their range. The experimental value 16.956 Ȧ is from McKenzie et al.
A look at Table 3 will show that except for the single line at λ(14.1 → 2.0) = 13.576 Ȧ [the observed wavelength closest to which is λ = 13.551 Ȧ, so that ∆λ = 0.025 which is just outside the accuracy limit of], series B is observed in its entirety. The single missing line could also account for the line at λ(12.0 → 2.1) = 13.589 Ȧ in Series A. For the latter series, however, two additional lines λ(5.0 → 2.1) and λ(6.0 → 2.1) are not observed. Other than this, all the lines in this series are also observed. Series C is degenerate with Series A except for the first line λ(3.2 → 2.1) = 21.613 Ȧ, the counterpart of the line λ(3.0 → 2.1) = 21.608 Ȧ. These two lines are nearly degenerate and correspond to the observed line at 21.602 Ȧ  or 21.60 Ȧ  .
4. Possible Connection with the FIP Effect
It may be noted that 29 distinct lines from the Solar flare data have been used up in the above as evidence for our generalized Balmer series. The data of course contain hundreds of lines over a wide range in the X-ray region. Our object is not to suggest that they all originate from our pairing mechanism. In fact, following the pioneering works of Grotrian and Edlen, we take it for granted that these lines follow from the so-called forbidden transitions in highly ionized high-Z atoms, the presence of which in stellar plasmas is a natural consequence of their having been formed in the interior of stars. The sources of a great many of these lines have accordingly been identified and the effort in that direction continues. Nevertheless, we may mention that, of the 29 lines from the data which we matched with our calculated lines, eight remain unidentified (marked [?] in Table 3) and two have questionable identification (marked [?] after the possible identification).
The conventional approach is thus not all-encompassing. In fact, the lack of proper identification in this approach of as many as eight (possibly 10) lines out of 29 does suggest that these lines may have a different origin. But what evidence could there be for such an additional origin for the other 19 lines which have been identified in the conventional approach? Interestingly enough, some evidence is indeed there. As already noted, it comes from what is called the FIP effect   . The effect is widely observed, but we shall confine ourselves here to the part which relates to the emission lines from the Solar flares. According to conventional wisdom underlying the classical stellar atmospheres theory, the relative abundances of various elements in a star are not expected to show any variation in its upper layers (unlike in the interior where thermonuclear-process gradients exist). The relative abundance pattern for elements from the Solar photosphere should thus not be different from that obtained from the Solar flares or from the Solar coronal active regions. It turns out that this is not so. Elements with low FIP appear to be relatively more abundant in regions with T ≥ 106 K than in the photosphere (T ≈ 6000 K). How does one infer this? One compares the relative intensity pattern of emission lines for the same elements from a) the high temperature sources like the flares and the coronal active regions and b) the photosphere. One finds that the intensities of lines corresponding to the low FIP elements from the former are anomalously enhanced. This may be interpreted to mean that these low FIP elements are relatively more abundant in the high temperature sources, which is in disagreement with the classical stellar atmospheres theory. Let us now look at the situation from the point of view of this study. We observe that of the 19 lines predicted by us, which are also identifiable as lines from known elements, almost all can be traced to a single low FIP element Fe (see Table 3) in its ionized state given by Fe XIX. In the light of the FIP effect, these are thus the lines which show anomalous enhancement. Clearly, such enhancement can now be attributed to an additional mechanism for their origin viz. the pairing mechanism considered herein. Note that such enhancement will not be possible for any lines originating from the photosphere which exists at a temperature of about 6000 K at which the pairing does not take place. This opens up the possibility of leaving undisturbed the generally accepted (and empirically corroborated) picture of the uniformity of relative abundances of elements in the upper reaches of the Sun.
It may also be pointed out that the detailed theoretical calculations  of the intensity ratios of lines, specifically from Fe XIX, do not seem to match with the observed data from the flares. This mismatch between theory and observation may be another pointer to the existence of an additional mechanism for the origin of the above lines.
In the context of the present study, we note that the generally accepted and empirically corroborated picture of the uniformity of elements in the upper reaches of the Sun, e.g., the photosphere and the Solar flares and coronal regions, would imply that the relative intensities of spectral lines from various elements would not show any significant variation from one region to another. It turns out that this is not so. The reason is: while 19 of the 29 lines in the data analyzed here and attributed to our pairing mechanism at T = 4.26 × 106 K can also be identified with those from a single low FIP elements Fe XIX, their intensities in the flares region are found to be anomalously enhanced as compared with the intensities of lines from the photosphere. It then makes sense to conclude that this enhancement is due to the existence of another mechanism operative in the region of flares, but not in the photosphere. As has been argued above, our pairing mechanism takes place at temperatures exceeding about a million K and not at temperatures around 5000 K that characterize the photosphere. We further note that 8 (possibly 10) of the 29 lines in the data analyzed above are not identifiable (or have questionable identification) on the basis of transitions from Fe XIX or any other elements. This lends support to the view that they may well be due to the additional pairing mechanism presented here.
The author is grateful to Dr. Bernhard Haisch for pointing out the relevance of the FIP effect to the present work. He is also grateful to Dr. J.J. Drake and Dr. Anand Bhatia for providing some valuable information. Finally, it is a matter of great pleasure to thank Dr. G.P. Malik without whose encouragement the present investigation would not have been completed.
We evaluate here the expectation value, given by
where, in view of Equations (8) in the text,
and is given by Equation (9).
we obtain, after going through some elementary algebra,
With the help of the standard result
and the recurrence relation
we then get
Proceeding in a similar manner, we also obtain
We thus obtain 
Note that M is the summations in S1 and S2 denotes the highest k value for which ak ≠ 0, and is the magnitude of the binding energy corresponding to principal quantum number n and angular momentum l for which is being calculated. Note further that is already calculated vide Equation (10) by the Hill
determinant method, which also determines the value of M. The method does not require an explicit knowledge of coefficients ak, which can be determined by using the backward Miller algorithm  . To this end, we define
whence Equation (10) can be written as
We thus have
These equations have to be supplemented by
Equation (A14a) then determines, and given bl (which can be taken as unity because the ratio S2/S1 above will be independent of this choice), are determined via Equation (A14b). Equation (11) then enables the coefficients to be calculated.
Thus knowing Wn,l corresponding to any value of T as calculated from Equation (10), and the corresponding values of l and M, we can calculate the coefficients ak and, finally, using Equation (A10), the value of.
 Malik, G.P., Pande, L.K. and Varma, V.S. (1991) Astrophysical Journal, 379, 788-795.
 Laming, J.M., Drake, J.J. and Widing, K.G. (1995) Astrophysical Journal, 443, 416-422.
 Hautot, A. (1986) Physical Review D, 33, 437-443.
 Znojil, M. (1986) Physical Review D, 34, 1224-1228.
 Biswas, S.N., Datta, K., Saxena, R.P., Srivastava, P.K. and Varma, V.S. (1971) Physical Review D, 4, 3617-3620.
 Biswas, S.N., Datta, K., Saxena, R.P., Srivastava, P.K. and Varma, V.S. (1973) Journal of Mathematical Physics, 14, 1190-1195.
 Gautschi, W. (1967) SIAM Review, 9, 24-82.
 Malik, G.P., Malik, U. and Varma, V.S. (1993) Astrophysics and Space Science, 199, 299-315.
 Phillips, K.J.H., Fawcett, B.C., Kent, B.J., Gabriel, A.H., Leibbacher, J.W., et al. (1982) Astrophysical Journal, 256, 774-787.
 Acton, L.W., Bruner, M.E., Brown, W.A., Fawcett, B.C., Schweizer, W. and Speer, R.J. (1985) Astrophysical Journal, 291, 865-878.
 McKenzie, D.L., Landecker, P.B., Broussard, R.M., Rugge, H.R., Young, R.M., Feldman, U., et al. (1980) Astrophysical Journal, 241, 409-416.
 Bhatia, A.K., Fawcett, B.C., Lemen, J.R., Mason, H.E. and Phillips, K.J.N. (1989) Monthly Notices of Royal Astronomical Society, 240, 421-444.
 Laming, J.M. (2015) Living Rev. Solar Physics, 12, 2.