Jeans equations
The Jeans equations are a set of partial differential equations that describe the motion of a collection of stars in a gravitational field. The Jeans equations relate the second-order velocity moments to the density and potential of a stellar system for systems without collision. They are analogous to the Euler equations for fluid flow and may be derived from the collisionless Boltzmann equation. The Jeans equations can come in a variety of different forms, depending on the structure of what is being modelled. Most utilization of these equations has been found in simulations with large number of gravitationally bound objects.
History
The Jeans equations were originally derived by James Clerk Maxwell. However, they were first applied to astronomy by James Jeans in 1915 while working on stellar hydrodynamics. Since then, multiple solutions to the equations have been calculated analytically and numerically. Some notable solutions include a spherically symmetric solution, derived by James Binney in 1983 and asymmetric solutions found in 1995 by Richard Arnold.[1][2]
Mathematics
Derivation from Boltzmann equation
The Collisionless Boltzmann equation, also called the Vlasov Equation is a special form of Liouville' equation and is given by:[3]
[math]\displaystyle{ {\partial f \over \partial t}+v {\partial f \over \partial r}-{\partial \Phi \over \partial r}{\partial f \over \partial v}=0 }[/math] Or in vector form: [math]\displaystyle{ {\partial f \over \partial t}+\vec{v}\cdot\vec{\nabla}f-\vec{\nabla}\Phi\cdot{\partial f \over \partial \vec{v}}=0 }[/math]
Combining the Vlasov equation with the Poisson equation for gravity:[math]\displaystyle{ \nabla^2 \phi = 4\pi G \rho. }[/math]gives the Jeans equations.
More explicitly, If n=n(x,t) is the density of stars in space, as a function of position x = (x1,x2,x3) and time t, v = (v1,v2,v3) is the velocity, and Φ = Φ(x,t) is the gravitational potential, the Jeans equations may be written as:[4][5]
[math]\displaystyle{ \frac{\partial n }{\partial t} + \sum_i \frac{\partial(n \langle{v_i}\rangle)}{\partial x_i}=0, }[/math]
- [math]\displaystyle{ \frac{\partial(n \langle{v_j}\rangle)}{\partial t} + n \frac{\partial \Phi}{\partial x_j} + \sum_i \frac{\partial(n \langle{v_i v_j}\rangle)}{\partial x_i}= 0 \qquad (j=1, 2, 3.) }[/math]
Here, the ⟨...⟩ notation means an average at a given point and time (x,t), so that, for example, [math]\displaystyle{ \langle{v_1}\rangle }[/math] is the average of component 1 of the velocity of the stars at a given point and time. The second set of equations may alternately be written as
- [math]\displaystyle{ n \frac{\partial \langle{v_j}\rangle}{\partial t} + \sum_i n \langle{v_i}\rangle \frac{\partial{\langle{v_j}\rangle}}{\partial x_i} = -n \frac{\partial \Phi}{\partial x_j} - \sum_i \frac{\partial (n \sigma_{ij}^2)}{\partial x_i} \qquad (j=1, 2, 3.) }[/math]
where the spatial part of the stress–energy tensor is defined as: [math]\displaystyle{ \sigma_{ij}^2=\langle{v_i v_j}\rangle-\langle{v_i}\rangle \langle{v_j}\rangle }[/math] and measures the velocity dispersion in components i and j at a given point.
Some given assumptions regarding these equations include:
- The flow in phase space must conserve mass
- The density around a given star remains the same, or is incompressible
Notice that the Jeans equations contain 9 unknowns (3 average velocities and 6 stress tensor terms), but only 3 equations. This means that Jeans equations are not closed. To solve different systems, various assumptions are made about the stress tensor.[6]
Spherical Jeans equations
One fundamental usage of Jean's equation is in spherical gravitational bodies. In spherical coordinates, the equations are:[6]
[math]\displaystyle{ {\partial \rho\langle v_r \rangle \over \partial t} +{\partial \rho\langle v_r^2 \rangle \over \partial r} + {\rho\over r}[2\langle v_r^2\rangle -\langle v_\theta^2\rangle -\langle v_\phi^2\rangle] + \rho {\partial \Phi \over \partial r} =0 }[/math]
[math]\displaystyle{ {\partial \rho\langle v_\theta \rangle \over \partial t} +{\partial \rho\langle v_rv_\theta \rangle \over \partial r} + {\rho\over r}[3\langle v_rv_\theta\rangle +(\langle v_\theta^2\rangle -\langle v_\phi^2\rangle)\cot(\theta)] =0 }[/math]
[math]\displaystyle{ {\partial \rho\langle v_\phi \rangle \over \partial t} +{\partial \rho\langle v_rv_\phi \rangle \over \partial r} + {\rho\over r}[3\langle v_rv_\phi\rangle +2\langle v_\theta v_\phi \rangle\cot(\theta)] =0 }[/math]
Using the stress tensor with the assumption that it is diagonal and [math]\displaystyle{ \sigma_\theta^2 =\sigma_\phi^2 }[/math], can reduce these equations to a single simplified equation:
[math]\displaystyle{ {\partial (\rho \sigma_r^2 ) \over \partial r} + {2\rho\over r}[\sigma_r^2 - \sigma_\theta^2] +\rho{\partial \Phi \over \partial r} =0 }[/math]
Again, there are two unknown functions ([math]\displaystyle{ \sigma_r^2 (r) }[/math] and [math]\displaystyle{ \sigma_\theta^2(r) }[/math]) that require assumptions for the equation to be solved.
Applications
Jeans equation have found great utility in N-body simulation gravitational research.[7] The scale of these simulations can range in size from just our solar system to the entire universe. Using measurements of stellar number density and various kinematic values, parameters within the Jeans equations can be estimated. This allows for various analyses to be made through the lens of Jeans equations. This is particularly useful when simulating dark matter halo distributions, due to its isothermal, non-interactive behavior. Searches for structure in galaxy formation, dark matter formation, and universe formation can have observations supplemented with simulations using Jeans equations.
Milky Way dark matter halo
File:Artist’s impression of the expected dark matter distribution around the Milky Way.ogv An example of such an analysis is given by the constraints that can be placed on the dark matter halo within the Milky Way. Using Sloan Digital Sky Survey measurements of our Galaxy, researchers were able to simulate the dark matter halo distribution using Jeans equations.[8] By comparing measured values with Jeans equation simulation results, they confirmed the need for extra dark matter and placed limits on its ellipsoid size. They estimated the ratio of minor to major axis of this halo to be 0.47 [math]\displaystyle{ \pm }[/math] 0.14. This method has been applied to many other galactic halos [9] and have produced similar results regarding dark matter halo topology.
Simulation limitations
The limiting factor of these simulations however, has been the data required to approximate stress tensor parameter values that dictate the Jeans equations behavior. Additionally, some constraints can be placed on Jeans equation simulations in order to produce reliable results [10][11] Some of these limitations include a wavelength resolution requirement, variable gravitational softening, and a minimum vertical structure particle resolution.
See also
References
- ↑ Arnold, Richard (1995). "Axisymmetric solutions of the Jeans equation" (in en). Monthly Notices of the Royal Astronomical Society 276: 293–300. doi:10.1093/mnras/276.1.293. ISSN 1365-2966. Bibcode: 1995MNRAS.276..293A. https://academic.oup.com/mnras/article/276/1/293/974731/Axisymmetric-solutions-of-the-Jeans-equation.
- ↑ Cappellari, Michele (2020-06-01). "Efficient solution of the anisotropic spherically aligned axisymmetric Jeans equations of stellar hydrodynamics for galactic dynamics" (in en). Monthly Notices of the Royal Astronomical Society 494 (4): 4819–4837. doi:10.1093/mnras/staa959. ISSN 0035-8711. https://academic.oup.com/mnras/article/494/4/4819/5818356.
- ↑ "The Collisionless Boltzmann Equation". https://www.cv.nrao.edu/~jhibbard/students/CPower/dynamics/cbe/cbe.html.
- ↑ Galactic Dynamics (1 ed.). Princeton University Press. 1988. pp. 195–197. ISBN 0-691-08445-9. Bibcode: 1988gady.book.....B.
- ↑ Merritt, David (2013). Dynamics and Evolution of Galactic Nuclei. Princeton, NJ: Princeton University Press. https://openlibrary.org/works/OL16802359W/Dynamics_and_Evolution_of_Galactic_Nuclei.
- ↑ 6.0 6.1 Van Den Bosch, Frank. "Yale Astronomy Lecture Notes". http://www.astro.yale.edu/vdbosch/lecture7.pdf.
- ↑ Haas, Brian L.; Hash, David B.; Bird, Graeme A.; Lumpkin, Forrest E.; Hassan, H. A. (1994). "Rates of thermal relaxation in direct simulation Monte Carlo methods" (in en). Physics of Fluids 6 (6): 2191–2201. doi:10.1063/1.868221. ISSN 1070-6631. Bibcode: 1994PhFl....6.2191H. http://aip.scitation.org/doi/10.1063/1.868221.
- ↑ Loebman, Sarah R.; Ivezić, Željko; Quinn, Thomas R.; Governato, Fabio; Brooks, Alyson M.; Christensen, Charlotte R.; Jurić, Mario (2012-09-26). "Constraints on the Shape of the Milky Way Dark Matter Halo from Jeans Equations Applied to Sloan Digital Sky Survey Data" (in en). The Astrophysical Journal 758 (1): L23. doi:10.1088/2041-8205/758/1/l23. ISSN 2041-8205. Bibcode: 2012ApJ...758L..23L. https://iopscience.iop.org/article/10.1088/2041-8205/758/1/L23.
- ↑ Adams, Joshua J.; Gebhardt, Karl; Blanc, Guillermo A.; Fabricius, Maximilian H.; Hill, Gary J.; Murphy, Jeremy D.; van den Bosch, Remco C. E.; van de Ven, Glenn (2012-01-03). "The Central Dark Matter Distribution of NGC 2976". The Astrophysical Journal 745 (1): 92. doi:10.1088/0004-637x/745/1/92. ISSN 0004-637X. Bibcode: 2012ApJ...745...92A. http://dx.doi.org/10.1088/0004-637x/745/1/92.
- ↑ Nelson, Andrew F. (2006). "Numerical requirements for simulations of self-gravitating and non-self-gravitating discs" (in en). Monthly Notices of the Royal Astronomical Society 373 (3): 1039–1073. doi:10.1111/j.1365-2966.2006.11119.x. ISSN 0035-8711. Bibcode: 2006MNRAS.373.1039N. https://academic.oup.com/mnras/article-lookup/doi/10.1111/j.1365-2966.2006.11119.x.
- ↑ Truelove, J. Kelly; Klein, Richard I.; McKee, Christopher F.; Holliman Ii, John H.; Howell, Louis H.; Greenough, Jeffrey A. (1997). "The Jeans Condition: A New Constraint on Spatial Resolution in Simulations of Isothermal Self-gravitational Hydrodynamics" (in en). The Astrophysical Journal 489 (2): L179–L183. doi:10.1086/310975. Bibcode: 1997ApJ...489L.179T.
 |  |
