Friday, May 31, 2019

Selected Questions - Answers From All Experts Astronomy Forum (Spiral Galaxy Shape)

Question:  I understand something called 'density waves' is responsible for shaping spiral galaxies like the Milky Way. Can you explain in detail how this happens?  - Margie, Mission Viejo, CA

Answer:  First, it's well to bear in mind that the particular density model (for a given galaxy) will vary depending on the conditions for that galaxy. In principle then, there can be differing density wave models, and these pertain to a number of variables, factors such as: how tightly wound the spiral is (there are different grades which are assigned, e.g. Sa, Sb, Sc etc.), the degree of axial symmetry of the galaxy, and the modeling assumptions - in particular the potential-gravitational fields (V(r)) imposed on the system which determine the locations of orbital resonance in conjunction with the equations used. 

The gravitational potential energy is defined according to:

V(r) = - GMm/r

where G is the Newtonian gravitational constant, m is the unit or elemental mass within the galaxy, M is associated with the central mass concentration, and r is the distance from m to M. The negative sign, as usual, indicates a bound system.
Interestingly, the use of density wave model development is largely contingent on the Boltzmann equation which is also used in space plasma physics. (The equation is used to describe the evolution of the velocity distribution of particles in space and time, cf. 

d f(x,v,t)/dt+ v*grad(x) f(x,v,t) +a*grad(v) f(x,v,t)=Col

Now, technically, the Boltzmann eqn. is applied to FLUIDS and for that purpose, the galaxies to which density wave approaches are applied are modeled firstly in the fluid format. (It is easier when dealing with an agglomeration of some 100 or 200 billion separate stars and associated orbits to think of them as comprising a "fluid" as opposed to say, 100 billion separate bodies to be treated in a 100 billion -body problem of celestial mechanics!) 

The referencing of stars, their locations and movements meanwhile embodies particulate approaches that are more kinematical in nature (but often less amenable to consistency with the density wave approach). Orbital assumptions, declarations are not simple by any means, and merely because a source says or asserts that "The stars in the inner part of a galaxy move faster than the density wave/s and the stars in the outer parts move slower than the wave/s." should not be taken too literally without posing a lot of further questions. (One could argue here that "apples" and "oranges" are being compared because the two entities, stars and density waves arise from differing backgrounds - kinematic-particle based and fluid mechanical, wave based.)

For example, what class is the spiral? How tightly wound? One must recognize too that an orbit in a spiral galaxy that appears closed (e.g. elliptic) in one reference frame may not be so in another. As an example, assume the (polar) coordinates for a galactic rotating frame are given as (r, φ) with:

dφ/dt = dΘ/dt - r

where  Ω r  is the angular velocity of the rotating frame. 

Then orbits are described by a Hamiltonian, H (the addition of the potential and kinetic energies of the system)  The point is that H can change depending on the coordinates, and what is presented for the previous frame as H = E - J Ω  may well be different for another frame.

Second, we see from this that the question as to why the spiral pattern is not affected by stars much further out cannot really be properly answered unless a full vetting of the assumed density waves for the particular galaxy is presented. In this sense, one recognizes that a full analysis of density waves for a galaxy - call it "Barred G1"- is needed before one can say stars in a given G1 region (e.g. inner or outer) "move faster or more slowly" than the waves at that place. We need to know then: the physical conditions for the establishment of the density waves at location r1 in G1 and r10 in G1 where the r's denote radial distances from the center with r10 = 10 (r1).

Another problem (or complexity) in dealing with density wave models is the fact they are mainly based on the mode chosen for particular dynamical wave equations that can be applied to the fluid framework. (In generic dynamical terms, a "mode" is a standing wave that can be supported by a disk of given dimensions, mass.) More broadly, most astronomers who work in this specialist area use the term interchangeabley with Fourier m-component. (And it should be understood here that one of the main tools is Fourier analysis of the waves, but alas Fourier analysis is only taught usually to those who take advanced Calculus or analysis courses).

As an example, a particular Fourier coefficient, call it a_n, applicable to a wave - may be defined: 

a_n = 1/π    ò  f(x) cos mx dx

 Integrated from (-π) to π.

Where m is the Fourier m-component. 

What types of modes can one have in these models? One is the "global" or m= 1 mode. Then there are the unstable (m= 2) modes.

Whether one mode or another appears (or is used in a spiral galaxy modelling) is critical since it may well determine at what stage a barred spiral develops, if at all. Alas, another complexity enters here since mode analysis is not simply a stand alone but also incorporates a subtle aspect called "marginal stability analysis" wherein one will solve for a quantity Q and if it is very close to 1 one has the case of marginal stability and tightly wound modes or in the case of spirals, around the Sa class. The trouble is that when one seriously incorporates any heating of the disk for whatever reason (say a massive central black hole sucking up matter and generating much radiation) then the desired values of Q are soon out of range, making it impossible for a given spiral structure to sustain itself.

Lastly, whenever one considers density waves in galaxies, it's important to bear in mind there remain enormous stumbling blocks even when applied to the simplest models of galactic disks (e.g. "zero thickness" disks). One of these arises from potential (V(r) - seen earlier) theory. Thus, the perturbed gravitational field at one location depends on the *density perturbation* at every other location. How will you know, ab initio, that the density perturbation at location r,φ, z say, does not accelerate the associated wave (in the fluid rest frame) to a higher velocity than any stars at the same or near location? You don't unless you investigate! What does it mean to "investigate"? It means a full bore mathematical modelling procedure to locate where all the (so-called "Lindblad") orbital resonances are, which ones can speed up the star, and also where the Landau damping regions are (which can impose a retardation of the waves).

In short, what I've shown is that density wave analysis as applied to galaxies is a field almost to itself in terms of being amenable to general understanding. And complexity is often compounded because the mere posing of a question to do with a particular spiral galaxy's form almost always introduces a number of tacit assumptions that may not be applicable at all. This, of course, is the difficulty when dealing in generalities, as opposed to specific cases, examples.


No comments: