User:Tohline/Appendix/Ramblings/Additional Analytically Specified Eigenvectors for ZeroZero Bipolytropes
From VisTrailsWiki
Contents 
Searching for Additional Eigenvectors of ZeroZero Bipolytropes
This chapter is an extension of two accompanying discussions: The original discovery and detailed derivation; and the more readable, summary.
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
In our accompanying summary, we have demonstrated how analytically specified eigenvectors can be constructed for the mode labeled, . This was done by specifying , then solving a quartic equation for . Shortly after completing this summary chapter, we noticed that an alternate approach may be to specify , then solve for ; and this path may be simpler because it may only involve solution of a quadratic equation. (Actually, we later have realized that the relevant equation is cubic, rather than quadratic. This is nevertheless simpler than the quartic equation.) If this proves to be the case, then it may also be possible to analytically construct eigenvectors of additional modes. Let's see.
CAUTION! All of the "core" eigenfunction expressions are suspect because (I think!) they have been developed herein under the assumption that , whereas the relation should read,

Seek Alternate Solution to Mode21
Setup21
According to STEP 4 in our accompanying summary discussion, we need to solve the following "derivative matching" expression:






where, recognizing that,


















Here, we assume that is specified and seek the corresponding value of . Given that the LHS of this matching relation is known once has been specified, in order to simplify notation we will also define,



Then the matching relation becomes,




































Solve Cubic Equation
Here we draw from a separate discussion of solutions to a cubic equation.
Using in place of , this "derivative matching" relation can be written in the form of a standard cubic equation. Specifically,
where,
As is well known and documented — see, for example Wolfram MathWorld or Wikipedia's discussion of the topic — the roots of any cubic equation can be determined analytically. In order to evaluate the root(s) of our particular cubic equation, we have drawn from the utilitarian online summary provided by Eric Schechter at Vanderbilt University. For a cubic equation of the general form,
a real root is given by the expression,
where, and (There is also a pair of imaginary roots, but they are irrelevant in the context of our overarching astrophysical discussion.) 
Upon evaluation, we have found that the expression inside of the square root is negative over the region of parameter space that is of most physical interest. Hence, we need to call upon a separate discussion in which the cube root of complex numbers was discussed.
We'll shift to Wolfram's notation; specifically,
Then, after defining, and Wolfram states that the three roots of the cubic equation are (the first one being identical to the "real" root identified above),
Now, whenever is intrinsically negative, we need to treat both and as complex numbers. If we define,
then we can write,
As is explained in this online resource, both and must formally submit to three separate roots tagged by the integer index, . Working only with the root for both, we find that the above expressions for the three roots of our cubic equation become,

We have deduced empirically that is the root that is physically relevant in our case. That is to say, for a given ,
In turn, for a given value of , the corresponding value of is obtained via the relation, .
The righthand panel of Figure 1 presents a plot of this function; actually what has been plotted is the inverted relation, . The open, red circular markers trace the portion of the function that provides physically viable solutions, in the sense that the corresponding value of lies between the values, negative one and three; the filled, light=blue circular markers identify roots of the cubic equation that are not physically viable.
In the lefthand panel of Figure 1, we redisplay a plot that has been discussed in an accompanying chapter. It contains a plot (blue markers) of the same same function but, this time, as determined from the root of a quartic equation. In order to illustrate more clearly that the two curves are the same, we have plotted the quartic solution (small, purple circular markers) on top of the cubic solution in the righthand panel.
Illustration21
Analytically Definable Eigenvectors in Bipolytropes Quantum Numbers: 


Analyzable Model Sequence  One Example Eigenfunction 
TopLeft Panel: Plotted points show how the location of the core/envelope interface, , varies with — where is the adiabatic exponent of the envelope — in equilibrium models that are amenable to analytic modal analysis for quantum numbers, . Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the core [see bottomleft panel] falls inside (alternatively, outside) the physically viable range, namely, . The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.
BottomLeft Panel: Plotted points show how varies with over the physically viable parameter range, . Both axes have been flipped so that incompressible models lie on the left/bottom while isothermal models lie on the right/top. The core is more compressible than the envelope in models that lie above and to the left of the blackdashed, diagonal line. The yellow circular marker identifies the same example model as it does in the topleft panel.
TopRight Panel: Displays — as a function of the fractional radius, — the analytically determined eigenfunction for the mode in the model identified by the yellow circular marker in both lefthand panels, for which,
and, correspondingly (via our separate solution of the governing quartic equation), . (This same plot appears in Figure 1 of our accompanying summary.) Specifically, over the radial interval, , the green markers trace the core's contribution to the combined eigenfunction, namely,



and, over the radial interval, , the purple markers trace the envelope's contribution to the combined eigenfunction, namely,



where the coefficients, , are as defined above in terms of the parameter, . The corresponding eigenfrequency is, from the perspective of the core,



and, from the perspective of the envelope,



where the relevant density ratio is, .
BottomRight Panel: The green and purple markers present the same eigenfunctionamplitude information, , as in the TopRight panel, but on a logarithmic scale. Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,



where, for plotting purposes, we have used, , and have set to a value that ensures that at the outer edge. In this type of logamplitude plot, the eigenfunction's various nodes — that is, radial locations where passes through zero — are highlighted; here, specifically, there is one node inside the core and one node resides in the envelope. Using the vertical coordinate to represent, instead, the configuration's massdensity normalized to its central value, , the solid black line segments trace the unperturbed density distribution throughout this specific bipolytrope. Throughout the core, ; then, at the location of the interface , the density abruptly drops to its envelope value .
NOTE: As may be ascertained from our general discussion of the structural properties of bipolytropes, equilibrium "zerozero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, ; and for any chosen value of , the envelope/core density ratio can, in principle, be set to any value, . We have not, however, been able to analytically solve the relevant pair of linearadiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is limited by the constraint,



The blackdotted curve in the BottomRight Panel displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this blackdotted curve for any equilibrium model — such as the example illustrated here — that has an analytically prescribable radialoscillation eigenvector.
Try Mode22
Setup22
In this case we need to replace the core eigenfunction segment that was specified in STEP 3 in our separate discussion of mode with the following:



So, following the procedure outlined in STEP 4 of our separate discussion, the requirement that the first derivatives match at the interface becomes,



where,






So, as described above, in order to determine for any specified value of — or, equivalently, — we simply need to replace with in the definition of the coefficient of each term in the cubic expression, where,



Finally, as explained in our summary discussion, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,



Illustration22
Analytically Definable Eigenvectors in Bipolytropes Quantum Numbers: 


Analyzable Model Sequence  One Example Eigenfunction 
TopLeft Panel: Plotted points show how the location of the core/envelope interface, , varies with — where is the adiabatic exponent of the envelope — in equilibrium models that are amenable to analytic modal analysis for quantum numbers, . Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the core [see bottomleft panel] falls inside (alternatively, outside) the physically viable range, namely, . The red marker that is farthest to the left identifies the model whose analytically determined eigenfunction is displayed, as an example, in the righthand panels.
BottomLeft Panel: Plotted points (only 2, here!) show how varies with over the physically viable parameter range, . Both axes have been flipped so that incompressible models lie on the left/bottom while isothermal models lie on the right/top. The core is more compressible than the envelope in models that lie above and to the left of the blackdashed, diagonal line.
TopRight Panel: Displays — as a function of the fractional radius, — the analytically determined eigenfunction for the mode in the model identified by the red circular marker that is farthest to the left in both lefthand panels, for which,
and, correspondingly, . Specifically, over the radial interval, , the green markers trace the core's contribution to the combined eigenfunction, namely,



and, over the radial interval, , the purple markers trace the envelope's contribution to the combined eigenfunction, namely,



where the coefficients, , are as defined above in terms of the parameter, . The corresponding eigenfrequency is, from the perspective of the core,



and, from the perspective of the envelope,



where the relevant density ratio is, .
BottomRight Panel: The green and purple markers present the same eigenfunctionamplitude information, , as in the TopRight panel, but on a logarithmic scale. Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,



where, for plotting purposes, we have used, , and have set to a value that ensures that at the center of the configuration. In this type of logamplitude plot, the eigenfunction's various nodes — that is, radial locations where passes through zero — are highlighted; here, specifically, there are two nodes inside the core and none in the envelope, although one of the nodes in the core lies just inside the core/envelope interface. Using the vertical coordinate to represent, instead, the configuration's massdensity normalized to its central value, , the solid black line segments trace the unperturbed density distribution throughout this specific bipolytrope. Throughout the core, ; then, at the location of the interface , the density abruptly drops to its envelope value .
NOTE: As may be ascertained from our general discussion of the structural properties of bipolytropes, equilibrium "zerozero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, ; and for any chosen value of , the envelope/core density ratio can, in principle, be set to any value, . We have not, however, been able to analytically solve the relevant pair of linearadiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is limited by the constraint,



The blackdotted curve in the BottomRight Panel displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this blackdotted curve for any equilibrium model — such as the example illustrated here — that has an analytically prescribable radialoscillation eigenvector.
Try Mode31
Setup31
In this case we need to replace the envelope eigenfunction segment that was specified in STEP 2 in our separate discussion of mode with the following: For the case of , this means that, throughout the envelope, the eigenfunction is,



where, the values of the newly introduced coefficients,



























Then, after defining,



the matching relation becomes,






























Root of Quartic Equation
To solve this equation analytically, we follow the Summary of Ferrari's method that is presented in Wikipedia's discussion of the Quartic Function to identify the roots of an arbitrary quartic equation.
First, we adopt the shorthand notation:
where, in our particular case,

Now, define,
Then the four roots of the quartic equation are,
It is this fourth root that interests us, here. 
We have determined empirically that, in our specific case, the quantity, is negative over the range of physically interesting values of . Hence, the quantity, , is necessarily complex. Let's work carefully through a determination of and, by consequence, , in this situation.
where,
We therefore can state that, in the complex plane, the three roots of this expression are,
where,
Focusing on the simplest root, we have,
Because this expression does not contain an imaginary component, we understand that is real. 
Finally, as explained in our summary discussion, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,






where, the last expression follows from plugging in the desired mode's quantum numbers, , and, again, using the relation,






Illustration31
Analytically Definable Eigenvectors in Bipolytropes Quantum Numbers: 


Analyzable Model Sequence  One Example Eigenfunction 
TopLeft Panel: Plotted points show how the location of the core/envelope interface, , varies with — where is the adiabatic exponent of the envelope — in equilibrium models that are amenable to analytic modal analysis for quantum numbers, . Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the core [see bottomleft panel] falls inside (alternatively, outside) the physically viable range, namely, . The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.
BottomLeft Panel: Plotted points show how varies with over the physically viable parameter range, . Both axes have been flipped so that incompressible models lie on the left/bottom while isothermal models lie on the right/top. The core is more compressible than the envelope in models that lie above and to the left of the blackdashed, diagonal line. The yellow circular marker identifies the same example model as it does in the topleft panel.
TopRight Panel: Displays — as a function of the fractional radius, — the analytically determined eigenfunction for the mode in the model identified by the yellow circular marker in both lefthand panels, for which,
and, correspondingly, . Specifically, over the radial interval, , the green markers trace the core's contribution to the combined eigenfunction, namely,



and, over the radial interval, , the purple markers trace the envelope's contribution to the combined eigenfunction, namely,



where the coefficients, , are as defined above in terms of the parameter, . The corresponding eigenfrequency is, from the perspective of the core,



and, from the perspective of the envelope,



where the relevant density ratio is, .
BottomRight Panel: The green and purple markers present the same eigenfunctionamplitude information, , as in the TopRight panel, but on a logarithmic scale. Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,



where, for plotting purposes, we have used, , and have set to a value that ensures that at the outer edge. In this type of logamplitude plot, the eigenfunction's various nodes — that is, radial locations where passes through zero — are highlighted; here, specifically, there is one node inside the core and two nodes reside in the envelope. Using the vertical coordinate to represent, instead, the configuration's massdensity normalized to its central value, , the solid black line segments trace the unperturbed density distribution throughout this specific bipolytrope. Throughout the core, ; then, at the location of the interface , the density abruptly drops to its envelope value .
NOTE: As may be ascertained from our general discussion of the structural properties of bipolytropes, equilibrium "zerozero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, ; and for any chosen value of , the envelope/core density ratio can, in principle, be set to any value, . We have not, however, been able to analytically solve the relevant pair of linearadiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is limited by the constraint,



The blackdotted curve in the BottomRight Panel displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this blackdotted curve for any equilibrium model — such as the example illustrated here — that has an analytically prescribable radialoscillation eigenvector.
Check Surface Boundary Condition
Let's determine analytically the logarithmic derivative of the envelope segment of the eigenfunction, at the surface of the configuration. Given that,



where,









we have,






This means that,


















which precisely matches the desired surface boundary condition for that has been detailed in an accompanying summary discussion.
Material that appears after this point in our presentation is under development and therefore
may contain incorrect mathematical equations and/or physical misinterpretations.
 Go Home 
Try Mode32
Setup32
Here, we set the quantum numbers to . In this case, as when addessing Mode22, above, the core eigenfunction segment should be,



and, as when addressing Mode31, above, throughout the envelope, the eigenfunction should be,



where,









In establishing the interface matching conditions, we should be able to follow the same steps as in our search for "31" modes, but replace with,



And, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,






Illustration32
Analytically Definable Eigenvectors in Bipolytropes Quantum Numbers: 


Analyzable Model Sequence  One Example Eigenfunction 
TopLeft Panel: Plotted points show how the location of the core/envelope interface, , varies with — where is the adiabatic exponent of the envelope — in equilibrium models that are amenable to analytic modal analysis for quantum numbers, . Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the core [see bottomleft panel] falls inside (alternatively, outside) the physically viable range, namely, . The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.
BottomLeft Panel: Plotted points show how varies with over the physically viable parameter range, . Both axes have been flipped so that incompressible models lie on the left/bottom while isothermal models lie on the right/top. The core is more compressible than the envelope in models that lie above and to the left of the blackdashed, diagonal line. The yellow circular marker identifies the same example model as it does in the topleft panel.
TopRight Panel: Displays — as a function of the fractional radius, — the analytically determined eigenfunction for the mode in the model identified by the yellow circular marker in both lefthand panels, for which,
and, correspondingly, . Specifically, over the radial interval, , the green markers trace the core's contribution to the combined eigenfunction as defined by the expression for , given above; and, over the radial interval, , the purple markers trace the envelope's contribution to the combined eigenfunction as defined by the expression for , given above.
The corresponding eigenfrequency is, from the perspective of the core,



and, from the perspective of the envelope,



where the relevant density ratio is, .
BottomRight Panel: The green and purple markers present the same eigenfunctionamplitude information, , as in the TopRight panel, but on a logarithmic scale. Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,



where, for plotting purposes, we have used, , and have set to a value that ensures that at the enter of the configuration. In this type of logamplitude plot, the eigenfunction's various nodes — that is, radial locations where passes through zero — are highlighted; here, specifically, there is one node inside the core and one node resides in the envelope. Using the vertical coordinate to represent, instead, the configuration's massdensity normalized to its central value, , the solid black line segments trace the unperturbed density distribution throughout this specific bipolytrope. Throughout the core, ; then, at the location of the interface , the density abruptly drops to its envelope value .
NOTE: As may be ascertained from our general discussion of the structural properties of bipolytropes, equilibrium "zerozero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, ; and for any chosen value of , the envelope/core density ratio can, in principle, be set to any value, . We have not, however, been able to analytically solve the relevant pair of linearadiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is limited by the constraint,



The blackdotted curve in the BottomRight Panel displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this blackdotted curve for any equilibrium model — such as the example illustrated here — that has an analytically prescribable radialoscillation eigenvector.
Try Mode33
Setup33
Here, we set the quantum numbers to . In this case, drawing from our summary of Sterne's presentation, the core eigenfunction segment should be,



and, as when addressing Mode31, above, throughout the envelope, the eigenfunction should be,



where,









In establishing the interface matching conditions, we should be able to follow the same steps as in our search for "31" modes, but replace with,



And, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,






Figure 4: Results for mode  
[[850pxResults for (ell,j) = (3,3)]] 

Illustration33
Analytically Definable Eigenvectors in Bipolytropes Quantum Numbers: 


Analyzable Model Sequence  One Example Eigenfunction 
TopLeft Panel: Plotted points show how the location of the core/envelope interface, , varies with — where is the adiabatic exponent of the envelope — in equilibrium models that are amenable to analytic modal analysis for quantum numbers, . Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the core [see bottomleft panel] falls inside (alternatively, outside) the physically viable range, namely, . The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.
BottomLeft Panel: Plotted points show how varies with over the physically viable parameter range, . Both axes have been flipped so that incompressible models lie on the left/bottom while isothermal models lie on the right/top. The core is more compressible than the envelope in models that lie above and to the left of the blackdashed, diagonal line. The yellow circular marker identifies the same example model as it does in the topleft panel.
TopRight Panel: Displays — as a function of the fractional radius, — the analytically determined eigenfunction for the mode in the model identified by the yellow circular marker in both lefthand panels, for which,
and, correspondingly, . Specifically, over the radial interval, , the green markers trace the core's contribution to the combined eigenfunction as defined by the expression for , given above; and, over the radial interval, , the purple markers trace the envelope's contribution to the combined eigenfunction as defined by the expression for , given above.
The corresponding eigenfrequency is, from the perspective of the core,



and, from the perspective of the envelope,



where the relevant density ratio is, .
BottomRight Panel: The green and purple markers present the same eigenfunctionamplitude information, , as in the TopRight panel, but on a logarithmic scale. Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,



where, for plotting purposes, we have used, , and have set to a value that ensures that at the outer edge. In this type of logamplitude plot, the eigenfunction's various nodes — that is, radial locations where passes through zero — are highlighted; here, specifically, there are three nodes inside the core and one node resides in the envelope. Using the vertical coordinate to represent, instead, the configuration's massdensity normalized to its central value, , the solid black line segments trace the unperturbed density distribution throughout this specific bipolytrope. Throughout the core, ; then, at the location of the interface , the density abruptly drops to its envelope value .
NOTE: As may be ascertained from our general discussion of the structural properties of bipolytropes, equilibrium "zerozero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, ; and for any chosen value of , the envelope/core density ratio can, in principle, be set to any value, . We have not, however, been able to analytically solve the relevant pair of linearadiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is limited by the constraint,



The blackdotted curve in the BottomRight Panel displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this blackdotted curve for any equilibrium model — such as the example illustrated here — that has an analytically prescribable radialoscillation eigenvector.
Summary Plots
Related Discussions
© 2014  2021 by Joel E. Tohline 