User:Tohline/MathProjects/EigenvalueProblemN1
From VisTrailsWiki
Contents 
Find Analytic Solutions to an Eigenvalue Problem
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
Note from J. E. Tohline to Students with Good Mathematical Skills: This is one of a set of welldefined research problems that are being posed, in the context of this online H_Book, as challenges to young, applied mathematicians. The astronomy community's understanding of the Structure, Stability, and Dynamics of stars and galaxies would be strengthened if we had, in hand, a closedform analytic solution to the problem being posed here. A solution can be obtained numerically with relative ease, but here the challenge is to find a closedform analytic solution. As is true with most meaningful scientific research projects, it is not at all clear whether this problem has such a solution. In my judgment, however, it seems plausible that a closedform solution can be discovered and such a solution would be of sufficient interest to the astronomical community that it would likely be publishable in a professional astronomy or physics journal. At the very least, this project offers an opportunity for a graduate student, an undergraduate, or even a talented highschool student (perhaps in connection with a mathematics science fair project?) to hone her/his research skills in applied mathematics. Also, I would be thrilled to include a solution to this problem — along with full credit to the solution's author — as a chapter in this online H_Book. Having retired from LSU, I am not in a position to financially support or formally advise students who are in pursuit of a highereducation degree. I would nevertheless be interested in sharing my expertise — and, perhaps, developing a collaborative relationship — with any individual who is interested in pursuing an answer to the mathematical research problem that is being posed here.
The Challenge
Formally, this is an eigenvalue problem.

Context
The challenge posed above is one of a set of closely related eigenvalue problems that arise in the context of the study of the pulsating stars and the governing 2^{nd}order ODE is often referred to as the Linear Adiabatic Wave Equation (LAWE). In the most general context, the LAWE takes the form,



where, and are both functions of that have different prescriptions for each specified astrophysics problem — see the table of examples presented below — and primes denote differentiation with respect to . The symmetries associated with this broad set of eigenvalue problems can perhaps be better appreciated by rearranging terms in the LAWE to obtain,



or, equivalently,



or, equivalently,



Properties of Analytically Defined Astrophysical Structures  

Model  
Uniformdensity  
Linear  
Parabolic  
Polytrope 
Drawing the expressions for , , and from the last row of this table and plugging them into this generic form of the LAWE leads to the specific statement of the astrophysically motivated eigenfunction problem presented above — inside the blueframed box. As is discussed in the subsection that follows, an analogous eigenvalue problem whose analytic solution is known comes from plugging expressions presented in the first row of this table into the generic form of the LAWE.
Analogous Problem with Known Analytic Solutions
Here is an analogous problem whose analytic solution is known. Anyone interested in tackling the challenge, provided above, should study — and even extend — the known set of solutions of this analogous problem. This exercise should provide at least partial preparation for addressing the above challenge.
Statement of the Problem
As above, the task here is to find one or more analytic expression(s) for the (eigen)function, — and, simultaneously, the unknown value of the (eigen)frequency, — that satisfies the following 2^{nd}order, ordinary differential equation:
Try a Polynomial Expression for the Eigenfunction
Let's guess that the proper eigenfunction is a polynomial expression in . Specifically, let's try a solution of the form,
truncated at progressively higher and higherorder terms.
Lowestorder mode (Mode 0)
Try,
in which case,



and,



So, the governing 2^{nd}order ODE reduces to,



which will be satisfied as long as, We conclude, therefore, that the eigenvector defining the lowestorder (the simplest) solution to the governing ODE has an eigenfunction given by,



with a corresponding eigenfrequency whose value is,



Second Guess
Try,
in which case,



and,



Plugging this trial eigenfunction into the governing 2^{nd}order ODE gives,






But this expression can be satisfied for all values of only if , in which case the trial eigenfunction reduces to the earlier solution, . We conclude, therefore, that our "second guess" does not generate a new solution to this eigenfunction problem.
Third Guess
Try,
in which case,



and,



Plugging this trial eigenfunction into the governing 2^{nd}order ODE therefore gives,






This relation will be satisfied for all values of if both expressions inside the square brackets are simultaneously zero, that is, if,



and, simultaneously,



 Mode 1:
 , in which case,
where,
Therefore,
Hence,
and, to within an arbitrary normalization factor,
Astrophysical Context
A star that is composed entirely of a gas for which the pressure varies as the square of the gas density will have an equilibrium structure that is defined by, what astronomers refer to as, an polytrope. Inside such an equilibrium structure, the density of the gas will vary with radial position according to the expression,
where, is the density at the center of the star, and,



where, is the radius of the equilibrium star. Notice that, according to this expression, the density will drop to zero when , in which case, . If a star of this type is nudged out of equilibrium — for example, squeezed slightly — in such a way that it maintains its spherical symmetry, the star will begin to undergo periodic, radial oscillations about its original equilibrium radius. The 2^{nd}order ODE whose solution is being sought in the above challenge is the equation that describes the behavior of these oscillations. In particular, the function,
describes the relative amplitude of the oscillation as a function of position, , within the star, and gives the frequency of the oscillation.
Related Discussions
© 2014  2019 by Joel E. Tohline 