Editing Apps/HayashiNaritaMiyama82 (section)

=Rotationally Flattened Isothermal Structures=
{| class="HNM82" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black"
|- 
! style="height: 150px; width: 150px; background-color:#ffeeee;" |[[H_BookTiledMenu#Spheroidal_.26_Spheroidal-Like|<b>Hayashi, Narita<br />&amp; Miyama's<br />Analytic Sol'n<br />(1982)</b>]]
|}
[https://ui.adsabs.harvard.edu/abs/1982PThPh..68.1949H/abstract Hayashi, Narita &amp; Miyama] (1982; hereafter HNM82) discovered an analytic solution to the equations that govern the structure of rotationally flattened, self-gravitating isothermal gas clouds.  Their solution describes a family of centrally condensed models whose degree of flattening ranges from a spherical structure to an infinitesimally thin disk.  For several reasons, I consider this to be one of the most remarkable discoveries &#8212; and, hence, one of the most significant papers &#8212; related to the structure of self-gravitating systems that was published in the decade of the '80s.  First, as has been [[Apps/MaclaurinSpheroids#Maclaurin_Spheroids_(axisymmetric_structure)|remarked earlier]], there is no particular reason why one should guess ahead of time that the equilibrium properties of ''any'' rotating, self-gravitating configuration should be describable in terms of analytic functions.  When dealing with compressible equations of state, such analytic solutions are rare even in the context of spherically symmetric structures, so it is impressive that HNM82 found a solution for rotationally flattened, isothermal configurations.  Second, about six months earlier the same year, [https://ui.adsabs.harvard.edu/abs/1982ApJ...259..535T/abstract Alar Toomre (1982)] published an independent discovery and strikingly independent derivation of exactly the same family of rotationally flattened, isothermal models.<sup>1</sup>  Third, these two independent derivations were motivated by a desire to better understand two quite different astrophysical environments:  The research of HNM82 was focused on star-forming gas clouds while Toomre's research was focused on the structure of elliptical galaxies and the dark-matter halos around spiral galaxies.  

Despite my assertion that HNM82 is one of the most significant papers to be published over the past few decades, citation indexes reveal that it is not widely referenced.  In large part I attribute this  to the fact that HNM82 was published in a Japanese journal ([https://academic.oup.com/ptp ''Progress of Theoretical Physics'']) that has only fairly recently made its archival articles available to the open-access, [https://ui.adsabs.harvard.edu Astrophysics Data System].

==Governing Relations==

As has been derived [[AxisymmetricConfigurations/SolutionStrategies|elsewhere]], for axisymmetric configurations that obey a barotropic equation of state, hydrostatic balance is governed by the following ''algebraic'' expression:
<div align="center">
<math>~H + \Phi_\mathrm{eff} = C_\mathrm{B} ,</math>
</div>
where <math>~C_\mathrm{B}</math> is the Bernoulli constant,
<div align="center">
<math>~\Phi_\mathrm{eff} \equiv \Phi + \Psi ,</math>
</div>
and <math>\Psi</math> is the relevant centrifugal potential.  For self-gravitating configurations, this algebraic expression must be satisfied in concert with a self-consistent solution of the Poisson equation, which for axisymmetric configurations may be written as,
<div align="center">
<math>
~\frac{1}{\varpi} \frac{\partial }{\partial\varpi} \biggl[ \varpi \frac{\partial \Phi}{\partial\varpi} \biggr] + \frac{\partial^2 \Phi}{\partial z^2} = 4\pi G \rho . 
</math><br />
</div>

==Supplemental Relations==
Following HNM82, we supplement the above-specified set of governing equations with an isothermal equation of state,
<div align="center">
<math>
~P = c_s^2 \rho ,
</math>
</div>
where <math>~c_s</math> is the sound speed.  Hence, to within an additive constant, the relevant expression for the enthalpy is,
<div align="center">
<math>
~H = \int \frac{dP}{\rho} = c_s^2 \ln\biggl(\frac{\rho}{\rho_0}\biggr) .
</math>
</div>

Also following HNM82, we impose a steady-state velocity flow-field that is described by a fluid with uniform rotational velocity, <math>v_\varphi</math>.  Drawing from our table of example [[AxisymmetricConfigurations/SolutionStrategies#SRPtable|''Simple rotation profiles'']], the centrifugal potential that describes this chosen flow-field is given by the expression,
<div align="center">
<math> 
~\Psi(\varpi) = - v_\varphi^2 \ln\biggl(\frac{\varpi}{\varpi_0}\biggr) .
</math> 
</div>

==Summary==
In summary, for the particular structures of interest to HNM82, the algebraic relation governing hydrostatic balance is,
<div align="center">
<math>
\Phi(\varpi,z) = C_\mathrm{B} - c_s^2 \ln\biggl[\frac{\rho(\varpi,z)}{\rho_0}\biggr] + v_\varphi^2 \ln\biggl(\frac{\varpi}{\varpi_0}\biggr) .
</math>
</div>
This must be solved in conjunction with the Poisson equation, which specifies the second key relationship between {{ Math/VAR_NewtonianPotential01 }} and {{ Math/VAR_Density01 }}, namely,
<div align="center">
<math>
\frac{1}{\varpi} \frac{\partial }{\partial\varpi} \biggl[ \varpi \frac{\partial \Phi(\varpi,z)}{\partial\varpi} \biggr] + \frac{\partial^2 \Phi(\varpi,z)}{\partial z^2} = 4\pi G \rho(\varpi,z) . 
</math><br />
</div>

==Solution==
===Rationalizing the Approach===
HNM82 provide a very clear and detailed description of the approach that they took to solving the above-identified set of simplified governing relations.  (Note that Hayashi is credited with deriving the analytic solution.)  In very general terms, one can understand the thought process that must have been going on in Hayashi's mind:
# For an isothermal gas cloud, the enthalpy is necessarily a logarithmic function of the density.
# For spherically symmetric, isothermal configurations, a solution to the governing relations exists in which {{ Math/VAR_Density01 }} can be expressed as a [[SSC/Structure/PowerLawDensity#Isothermal_Equation_of_State|power-law function]] of the radius, specifically, <math>~\rho \propto r^{-2}</math>; hence, the enthalpy displays a logarithmic dependence on the distance.
# Although one could attempt to derive equilibrium structures having a wide range of [[AxisymmetricConfigurations/SolutionStrategies#SRPtable|''Simple rotation profiles'']], it would seem wisest to select a centrifugal potential function that at least has the same ''form'' as the enthalpy; hence, the choice was made to impose <math>~v_\varphi = \mathrm{constant}</math> so that, like the enthalpy, the centrifugal potential would exhibit a logarithmic dependence on the distance.
# Three birds, so to speak, can be killed with one stone by ''guessing'' a 2D equilibrium density profile of the form <math>~\rho(\varpi,z) = g(\varpi,z)/\varpi^{2}</math>:  The <math>~\varpi^{-2}</math> dependence can be combined strategically with <math>~\varpi^{-2}</math> dependence of the centrifugal potential; although it depends on <math>\varpi</math> instead of <math>r</math>, the density profile will at least ''resemble'' the spherical solution; and &#8212; certainly the most critical realization &#8212; the Poisson equation, which for this problem is a 2D elliptic PDE, can be rewritten as a 1D ODE and solved analytically!

===HNM82 Derivation===

HNM82 ''guess'' a density distribution of the form (see their Eq.2.3),
<div align="center">
<math>
~\frac{\rho(\varpi,z)}{\rho_0} = g(\varpi,z) \biggl(\frac{\varpi}{\varpi_0}\biggr)^{-2} .
</math>
</div>
The algebraic expression defining hydrostatic balance then becomes,
<div align="center">
<math>
~\Phi(\varpi,z) = C_\mathrm{B} - c_s^2 \ln g(\varpi,z) + (2c_s^2 + v_\varphi^2) \ln\biggl(\frac{\varpi}{\varpi_0}\biggr) ;
</math>
</div>
and, after multiplying both sides by <math>~\varpi^2</math>, the Poisson equation becomes,
<div align="center">
<math>
~\varpi \frac{\partial }{\partial\varpi} \biggl[ \varpi \frac{\partial \Phi(\varpi,z)}{\partial\varpi} \biggr] + \varpi \frac{\partial}{\partial z} \biggl[ \varpi \frac{\partial \Phi(\varpi,z)}{\partial z} \biggr] = (4\pi G\rho_0 \varpi_0^2) g(\varpi,z) . 
</math><br />
</div>
Plugging the expression for <math>~\Phi</math> into the Poisson equation gives (see Eq. 2.4 of HNM82),
<div align="center">
<math>
~\varpi \frac{\partial }{\partial\varpi} \biggl[ \varpi \frac{\partial \ln g(\varpi,z)}{\partial\varpi} \biggr] + \varpi \frac{\partial}{\partial z} \biggl[ \varpi \frac{\partial \ln g(\varpi,z)}{\partial z} \biggr] = - \biggl[ \frac{4\pi G\rho_0 \varpi_0^2}{c_s^2} \biggr] g(\varpi,z) . 
</math><br />
</div>
HNM82 realized that since this equation "is invariant to the scale changes of both <math>~\varpi</math> and <math>~z</math>, it has a conformal solution such that <math>~g</math> is a function of <math>~z/\varpi</math> alone."  In particular, as HNM82 pointed out, by making the substitution,
<div align="center">
<math>
~\zeta \equiv \sinh^{-1}\biggl(\frac{z}{\varpi}\biggr)= \ln \biggl[ \frac{r+z}{\varpi} \biggr] ,
</math>
</div>
the above, 2D elliptic PDE (Poisson equation) can be written as the following, 1D second-order ODE:
 <div align="center">
<math>
~\frac{d^2\ln g(\zeta)}{d\zeta^2} = - \biggl[ \frac{4\pi G\rho_0 \varpi_0^2}{c_s^2} \biggr] g(\zeta) .
</math>
</div>
As presented by HNM82, the solution to this Poisson equation that meets the most physically reasonable boundary conditions at <math>~\zeta = 0</math> (''i.e.,'' <math>~g</math> is finite and <math>~dg/d\zeta = 0</math>) is,
<div align="center">
<math>
~g(\zeta) = \biggl[ \frac{c_s^2}{2\pi G\rho_0 \varpi_0^2} \biggr] \frac{\gamma^2}{\cosh^2(\gamma\zeta)} ,
</math>
</div>
where <math>~\gamma = [1 + v_\varphi^2/(2c_s^2)]</math> is a parameter that identifies an individual equilibrium structure from the ''family'' of allowed solutions.

==Structural Properties==

By analogy with the solution that was derived for a [[SSC/Structure/PowerLawDensity#Isothermal_Equation_of_State|spherically symmetric isothermal structure with a power-law density distribution]], we can associate the scale length <math>~\varpi_0</math> with the characteristic density <math>~\rho_0</math> at that location through the relation,
<div align="center">
<math>
~\rho_0 = \frac{c_s^2}{2\pi G \varpi_0^2 } ,
</math>
</div>
in which case we can write,
<div align="center">
<math>
~g(\zeta) = \frac{\gamma^2}{\cosh^2(\gamma\zeta)} .
</math>
</div>
With this definition in hand, the equilibrium models discovered by HNM82 exhibit the following properties.

* <span id="Density"><font color="red">Density Distribution</font></span>: 
:Confirming the expression presented as Eq. (3.1) in HNM82, the 2D density distribution for models with different values of the dimensionless parameter <math>~\gamma</math> is, 
<table align="center" border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
~\rho(\varpi,z) 
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
~\rho_0 g(\varpi,z) \biggl(\frac{\varpi}{\varpi_0}\biggr)^{-2} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
~\biggl[ \frac{c_s^2}{2\pi G \varpi_0^2} \biggr]\biggl(\frac{\varpi_0}{\varpi}\biggr)^{2} \frac{\gamma^2}{\cosh^2(\gamma\zeta)}  .
</math>
  </td>
</tr>
</table>

* <span id="Potential"><font color="red">Gravitational Potential</font></span>: 
:As presented in Eq. (3.4) of HNM82, to within an additive constant the 2D potential distribution for models with different values of the dimensionless parameter <math>~\gamma</math> is, 
<table align="center" border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
~\frac{\Phi(\varpi,z)}{2c_s^2} 
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
~\frac{C_\mathrm{B}}{2c_s^2} - \frac{1}{2} \ln g(\varpi,z) + \biggl(1 + \frac{v_\varphi^2}{2c_s^2} \biggr) \ln\biggl(\frac{\varpi}{\varpi_0}\biggr)
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
~\frac{C_\mathrm{B}}{2c_s^2} - \frac{1}{2} \ln \biggl[ \frac{\gamma^2}{\cosh^2(\gamma\zeta)} \biggr] + \gamma \ln\biggl(\frac{\varpi}{\varpi_0}\biggr) 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
~\frac{C_\mathrm{B}}{2c_s^2} + \ln \biggl[ \frac{1}{\gamma} \biggr] + \ln \biggl[ \cosh (\gamma\zeta) \biggr] + \ln\biggl(\frac{\varpi}{\varpi_0}\biggr)^\gamma 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
~C' + \ln \biggl[ \varpi^\gamma \cosh (\gamma\zeta) \biggr] .
</math>
  </td>
</tr>
</table>

:HNM82 point out that this last expression also can be written in terms of <math>~r</math> and <math>~z</math> as follows:

<table align="center" border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
~\frac{\Phi(\varpi,z)}{2c_s^2} 
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
~C' + \ln \biggl[ \frac{1}{2}(r+z)^\gamma + \frac{1}{2}(r-z)^\gamma \biggr] .
</math>
  </td>
</tr>
</table>


<table border="1" align="center" cellpadding="5" width="420px">
<tr>
  <th align="center">3D Renderings of Isothermal Disks from [https://ui.adsabs.harvard.edu/abs/1982PThPh..68.1949H/abstract Hayashi, Narita &amp; Miyama (1982)]</th>
</tr>
<tr><td align="center">
[[File:HNM82ThreeDisksAB.png|400px|Three HNM82 disks]]
</td></tr>
<tr>
  <td align="left">(Top) Edge-on view of three isothermal disks; (Bottom) Cut-away view of highest-density region of the same three disks.  Values of flattening parameter, <math>~\gamma</math>, shown in bottom-right corner of each frame.  
</table>