Magnetorotational instability

Last updated

The magnetorotational instability (MRI) is a fluid instability that causes an accretion disk orbiting a massive central object to become turbulent. It arises when the angular velocity of a conducting fluid in a magnetic field decreases as the distance from the rotation center increases. It is also known as the Velikhov–Chandrasekhar instability or Balbus–Hawley instability in the literature, not to be confused with the electrothermal Velikhov instability. The MRI is of particular relevance in astrophysics where it is an important part of the dynamics in accretion disks.

Contents

Gases or liquids containing mobile electrical charges are subject to the influence of a magnetic field. In addition to hydrodynamical forces such as pressure and gravity, an element of magnetized fluid also feels the Lorentz force where is the current density and is the magnetic field vector. If the fluid is in a state of differential rotation about a fixed origin, this Lorentz force can be surprisingly disruptive, even if the magnetic field is very weak. In particular, if the angular velocity of rotation decreases with radial distance the motion is unstable: a fluid element undergoing a small displacement from circular motion experiences a destabilizing force that increases at a rate which is itself proportional to the displacement. This process is known as the Magnetorotational Instability, or "MRI".

In astrophysical settings, differentially rotating systems are very common and magnetic fields are ubiquitous. In particular, thin disks of gas are often found around forming stars or in binary star systems, where they are known as accretion disks. Accretion disks are also commonly present in the centre of galaxies, and in some cases can be extremely luminous: quasars, for example, are thought to originate from a gaseous disk surrounding a very massive black hole. Our modern understanding of the MRI arose from attempts to understand the behavior of accretion disks in the presence of magnetic fields; it is now understood that the MRI is likely to occur in a very wide variety of different systems.

Discovery

The MRI was first noticed in a non-astrophysical context by Evgeny Velikhov in 1959 when considering the stability of Couette flow of an ideal hydromagnetic fluid. [1] His result was later generalized by Subrahmanyan Chandrasekhar in 1960. [2] [3] This mechanism was proposed by D. J. Acheson and Raymond Hide (1973) to perhaps play a role in the context of the Earth's geodynamo problem. [4] Although there was some follow-up work in later decades (Fricke, 1969; Acheson and Hide 1972; Acheson and Gibbons 1978), the generality and power of the instability were not fully appreciated until 1991, when Steven A. Balbus and John F. Hawley gave a relatively simple elucidation and physical explanation of this important process. [5]

Physical process

A simple model of MRI MRI spring model.png
A simple model of MRI

In a magnetized, perfectly conducting fluid, the magnetic forces behave in some very important respects as though the elements of fluid were connected with elastic bands: trying to displace such an element perpendicular to a magnetic line of force causes an attractive force proportional to the displacement, like a spring under tension. Normally, such a force is restoring, a strongly stabilizing influence that would allow a type of magnetic wave to propagate. If the fluid medium is not stationary but rotating, however, attractive forces can actually be destabilizing. The MRI is a consequence of this surprising behavior.

Consider, for example, two masses, mi ("inner") and mo ("outer") connected by a spring under tension, both masses in orbit around a central body, Mc. In such a system, the angular velocity of circular orbits near the center is greater than the angular velocity of orbits farther from the center, but the angular momentum of the inner orbits is smaller than that of the outer orbits. If mi is allowed to orbit a little bit closer to the center than mo, it will have a slightly higher angular velocity. The connecting spring will pull back on mi, and drag mo forward. This means that mi experiences a retarding torque, loses angular momentum, and must fall inward to an orbit of smaller radius, corresponding to a smaller angular momentum. mo, on the other hand, experiences a positive torque, acquires more angular momentum, and moves outward to a higher orbit. The spring stretches yet more, the torques become yet larger, and the motion is unstable! Because magnetic forces act like a spring under tension connecting fluid elements, the behavior of a magnetized fluid is almost exactly analogous to this simple mechanical system. This is the essence of the MRI .

A more detailed explanation

To see this unstable behavior more quantitatively, consider the equations of motion for a fluid element mass in circular motion with an angular velocity In general will be a function of the distance from the rotation axis and we assume that the orbital radius is The centripetal acceleration required to keep the mass in orbit is ; the minus sign indicates a direction toward the center. If this force is gravity from a point mass at the center, then the centripetal acceleration is simply where is the gravitational constant and is the central mass. Let us now consider small departures from the circular motion of the orbiting mass element caused by some perturbing force. We transform variables into a rotating frame moving with the orbiting mass element at angular velocity with origin located at the unperturbed, orbiting location of the mass element. As usual when working in a rotating frame, we need to add to the equations of motion a Coriolis force plus a centrifugal force The velocity is the velocity as measured in the rotating frame. Furthermore, we restrict our attention to a small neighborhood near say with much smaller than Then the sum of the centrifugal and centripetal forces is

(1)

to linear order in With our axis pointing radial outward from the unperturbed location of the fluid element and our axis pointing in the direction of increasing azimuthal angle (the direction of the unperturbed orbit), the and equations of motion for a small departure from a circular orbit are:

(2)

(3)

where and are the forces per unit mass in the and directions, and a dot indicates a time derivative (i.e., is the velocity, is the acceleration, etc.). Provided that and are either 0 or linear in x and y, this is a system of coupled second-order linear differential equations that can be solved analytically. In the absence of external forces, and , the equations of motion have solutions with the time dependence where the angular frequency satisfies the equation

(4)

where is known as the epicyclic frequency. In our solar system, for example, deviations from a sun-centered circular orbit that are familiar ellipses when viewed by an external viewer at rest, appear instead as small radial and azimuthal oscillations of the orbiting element when viewed by an observer moving with the undisturbed circular motion. These oscillations trace out a small retrograde ellipse (i.e. rotating in the opposite sense of the large circular orbit), centered on the undisturbed orbital location of the mass element.

The epicyclic frequency may equivalently be written which shows that it is proportional to the radial derivative of the angular momentum per unit mass, or specific angular momentum. The specific angular momentum must increase outward if stable epicyclic oscillations are to exist, otherwise displacements would grow exponentially, corresponding to instability. This is a very general result known as the Rayleigh criterion (Chandrasekhar 1961) for stability. For orbits around a point mass, the specific angular momentum is proportional to so the Rayleigh criterion is well satisfied.

Consider next the solutions to the equations of motion if the mass element is subjected to an external restoring force, where is an arbitrary constant (the "spring constant"). If we now seek solutions for the modal displacements in and with time dependence we find a much more complex equation for

(5)

Even though the spring exerts an attractive force, it may destabilize. For example, if the spring constant is sufficiently weak, the dominant balance will be between the final two terms on the left side of the equation. Then, a decreasing outward angular velocity profile will produce negative values for and both positive and negative imaginary values for The negative imaginary root results not in oscillations, but in exponential growth of very small displacements. A weak spring therefore causes the type of instability described qualitatively at the end of the previous section. A strong spring on the other hand, will produce oscillations, as one intuitively expects.

The spring-like nature of magnetic fields

The conditions inside a perfectly conducting fluid in motion is often a good approximation to astrophysical gases. In the presence of a magnetic field a moving conductor responds by trying to eliminate the Lorentz force on the free charges. The magnetic force acts in such a way as to locally rearrange these charges to produce an internal electric field of In this way, the direct Lorentz force on the charges vanishes. (Alternatively, the electric field in the local rest frame of the moving charges vanishes.) This induced electric field can now itself induce further changes in the magnetic field according to Faraday's law,

(6)

Another way to write this equation is that if in time the fluid makes a displacement then the magnetic field changes by

(7)

The equation of a magnetic field in a perfect conductor in motion has a special property: the combination of Faraday induction and zero Lorentz force makes the field lines behave as though they were painted, or "frozen," into the fluid. In particular, if is initially nearly constant and is a divergence-free displacement, then our equation reduces to

(8)

because of the vector calculus identity Out of these 4 terms, is one of Maxwell's equations. By the divergence-free assumption, . because B is assumed to be nearly constant. Equation 8 shows that changes only when there is a shearing displacement along the field line. To understand the MRI, it is sufficient to consider the case in which is uniform in vertical direction, and varies as Then

(9)

where it is understood that the real part of this equation expresses its physical content. (If is proportional to for example, then is proportional to )

A magnetic field exerts a force per unit volume on an electrically neutral, conducting fluid equal to Ampere's circuital law gives because Maxwell's correction is neglected in the MHD approximation. The force per unit volume becomes

(10)

where we have used the same vector calculus identity. This equation is fully general, and makes no assumptions about the strength or direction of the magnetic field. The first term on the right is analogous to a pressure gradient. In our problem it may be neglected because it exerts no force in the plane of the disk, perpendicular to The second term acts like a magnetic tension force, analogous to a taut string. For a small disturbance it exerts an acceleration given by force divided by mass, or equivalently, force per unit volume divided by mass per unit volume:

(11)

Thus, a magnetic tension force gives rise to a return force which is directly proportional to the displacement. This means that the oscillation frequency for small displacements in the plane of rotation of a disk with a uniform magnetic field in the vertical direction satisfies an equation ("dispersion relation") exactly analogous to equation 5 , with the "spring constant"

(12)

As before, if there is an exponentially growing root of this equation for wavenumbers satisfying This corresponds to the MRI. Notice that the magnetic field appears in equation 12 only as the product Thus, even if is very small, for very large wavenumbers this magnetic tension can be important. This is why the MRI is so sensitive to even very weak magnetic fields: their effect is amplified by multiplication by Moreover, it can be shown that MRI is present regardless of the magnetic field geometry, as long as the field is not too strong.

In astrophysics, one is generally interested in the case for which the disk is supported by rotation against the gravitational attraction of a central mass. A balance between the Newtonian gravitational force and the radial centripetal force immediately gives

(13)

where is the Newtonian gravitational constant, is the central mass, and is radial location in the disk. Since this so-called Keplerian disk is unstable to the MRI . Without a weak magnetic field, the flow would be stable.

For a Keplerian disk, the maximum growth rate is which occurs at a wavenumber satisfying is very rapid, corresponding to an amplification factor of more than 100 per rotation period. The nonlinear development of the MRI into fully developed turbulence may be followed via large scale numerical computation.

Applications and laboratory experiments

Interest in the MRI is based on the fact that it appears to give an explanation for the origin of turbulent flow in astrophysical accretion disks (Balbus and Hawley, 1991). A promising model for the compact, intense X-ray sources discovered in the 1960s was that of a neutron star or black hole drawing in ("accreting") gas from its surroundings (Prendergast and Burbidge, 1968). Such gas always accretes with a finite amount of angular momentum relative to the central object, and so it must first form a rotating disk — it cannot accrete directly onto the object without first losing its angular momentum. But how an element of gaseous fluid managed to lose its angular momentum and spiral onto the central object was not at all obvious.

One explanation involved shear-driven turbulence (Shakura and Sunyaev, 1973). There would be significant shear in an accretion disk (gas closer to the centre rotates more rapidly than outer disk regions), and shear layers often break down into turbulent flow. The presence of shear-generated turbulence, in turn, produces the powerful torques needed to transport angular momentum from one (inner) fluid element to another (farther out).

The breakdown of shear layers into turbulence is routinely observed in flows with velocity gradients, but without systematic rotation. This is an important point, because rotation produces strongly stabilizing Coriolis forces, and this is precisely what occurs in accretion disks . As can be seen in equation 5 , the K = 0 limit produces Coriolis-stabilized oscillations, not exponential growth. These oscillations are present under much more general conditions as well: a recent laboratory experiment (Ji et al., 2006) has shown stability of the flow profile expected in accretion disks under conditions in which otherwise troublesome dissipation effects are (by a standard measure known as the Reynolds number) well below one part in a million. All of this changes, however, are when even a very weak magnetic field is present. The MRI produces torques that are not stabilized by Coriolis forces. Large scale numerical simulations of the MRI indicate that the rotational disk flow breaks down into turbulence (Hawley et al., 1995), with strongly enhanced angular momentum transport properties. This is just what is required for the accretion disk model to work. The formation of stars (Stone et al., 2000), the production of X-rays in neutron star and black hole systems (Blaes, 2004), and the creation of active galactic nuclei (Krolik, 1999) and gamma ray bursts (Wheeler, 2004) are all thought to involve the development of the MRI at some level.

Thus far, we have focused rather exclusively on the dynamical breakdown of laminar flow into turbulence triggered by a weak magnetic field, but it is also the case that the resulting highly agitated flow can act back on this same magnetic field. Embedded magnetic field lines are stretched by the turbulent flow, and it is possible that systematic field amplification could result. The process by which fluid motions are converted to magnetic field energy is known as a dynamo (Moffatt, 1978); the two best studied examples are the Earth's liquid outer core and the layers close to the surface of the Sun. Dynamo activity in these regions is thought to be responsible for maintaining the terrestrial and solar magnetic fields. In both of these cases thermal convection is likely to be the primary energy source, though in the case of the Sun differential rotation may also play an important role. Whether the MRI is an efficient dynamo process in accretion disks is currently an area of active research (Fromang and Papaloizou, 2007).

There may also be applications of the MRI outside of the classical accretion disk venue. Internal rotation in stars (Ogilvie, 2007), and even planetary dynamos (Petitdemange et al., 2008) may, under some circumstances, be vulnerable to the MRI in combination with convective instabilities. These studies are also ongoing.

Finally, the MRI can, in principle, be studied in the laboratory (Ji et al., 2001), though these experiments are very difficult to implement. A typical set-up involves either concentric spherical shells or coaxial cylindrical shells. Between (and confined by) the shells, there is a conducting liquid metal such as sodium or gallium. The inner and outer shells are set in rotation at different rates, and viscous torques compel the trapped liquid metal to differentially rotate. The experiment then investigates whether the differential rotation profile is stable or not in the presence of an applied magnetic field.

A claimed detection of the MRI in a spherical shell experiment (Sisan et al., 2004), in which the underlying state was itself turbulent, awaits confirmation at the time of this writing (2009). A magnetic instability that bears some similarity to the MRI can be excited if both vertical and azimuthal magnetic fields are present in the undisturbed state (Hollerbach and Rüdiger, 2005). This is sometimes referred to as the helical-MRI, (Liu et al., 2006) though its precise relation to the MRI described above has yet to be fully elucidated. Because it is less sensitive to stabilizing ohmic resistance than is the classical MRI, this helical magnetic instability is easier to excite in the laboratory, and there are indications that it may have been found (Stefani et al., 2006). The detection of the classical MRI in a hydrodynamically quiescent background state has yet to be achieved in the laboratory, however.

The spring-mass analogue of the standard MRI has been demonstrated in rotating Taylor–Couette / Keplerian-like flow (Hung et al. 2019).

Related Research Articles

<span class="mw-page-title-main">Angular momentum</span> Conserved physical quantity; rotational analogue of linear momentum

In physics, angular momentum is the rotational analog of linear momentum. It is an important physical quantity because it is a conserved quantity – the total angular momentum of a closed system remains constant. Angular momentum has both a direction and a magnitude, and both are conserved. Bicycles and motorcycles, flying discs, rifled bullets, and gyroscopes owe their useful properties to conservation of angular momentum. Conservation of angular momentum is also why hurricanes form spirals and neutron stars have high rotational rates. In general, conservation limits the possible motion of a system, but it does not uniquely determine it.

Continuum mechanics is a branch of mechanics that deals with the deformation of and transmission of forces through materials modeled as a continuous medium rather than as discrete particles. The French mathematician Augustin-Louis Cauchy was the first to formulate such models in the 19th century.

<span class="mw-page-title-main">Torque</span> Turning force around an axis

In physics and mechanics, torque is the rotational analogue of linear force. It is also referred to as the moment of force. The symbol for torque is typically , the lowercase Greek letter tau. When being referred to as moment of force, it is commonly denoted by M. Just as a linear force is a push or a pull applied to a body, a torque can be thought of as a twist applied to an object with respect to a chosen point; for example, driving a screw uses torque, which is applied by the screwdriver rotating around its axis. A force of three newtons applied two metres from the fulcrum, for example, exerts the same torque as a force of one newton applied six metres from the fulcrum.

<span class="mw-page-title-main">Equations of motion</span> Equations that describe the behavior of a physical system

In physics, equations of motion are equations that describe the behavior of a physical system in terms of its motion as a function of time. More specifically, the equations of motion describe the behavior of a physical system as a set of mathematical functions in terms of dynamic variables. These variables are usually spatial coordinates and time, but may include momentum components. The most general choice are generalized coordinates which can be any convenient variables characteristic of the physical system. The functions are defined in a Euclidean space in classical mechanics, but are replaced by curved spaces in relativity. If the dynamics of a system is known, the equations are the solutions for the differential equations describing the motion of the dynamics.

<span class="mw-page-title-main">Angular velocity</span> Direction and rate of rotation

In physics, angular velocity, also known as angular frequency vector, is a pseudovector representation of how the angular position or orientation of an object changes with time, i.e. how quickly an object rotates around an axis of rotation and how fast the axis itself changes direction.

In physics, angular acceleration is the time rate of change of angular velocity. Following the two types of angular velocity, spin angular velocity and orbital angular velocity, the respective types of angular acceleration are: spin angular acceleration, involving a rigid body about an axis of rotation intersecting the body's centroid; and orbital angular acceleration, involving a point particle and an external axis.

The vorticity equation of fluid dynamics describes the evolution of the vorticity ω of a particle of a fluid as it moves with its flow; that is, the local rotation of the fluid. The governing equation is:

<span class="mw-page-title-main">Dynamo theory</span> Mechanism by which a celestial body generates a magnetic field

In physics, the dynamo theory proposes a mechanism by which a celestial body such as Earth or a star generates a magnetic field. The dynamo theory describes the process through which a rotating, convecting, and electrically conducting fluid can maintain a magnetic field over astronomical time scales. A dynamo is thought to be the source of the Earth's magnetic field and the magnetic fields of Mercury and the Jovian planets.

<span class="mw-page-title-main">Rigid body dynamics</span> Study of the effects of forces on undeformable bodies

In the physical science of dynamics, rigid-body dynamics studies the movement of systems of interconnected bodies under the action of external forces. The assumption that the bodies are rigid simplifies analysis, by reducing the parameters that describe the configuration of the system to the translation and rotation of reference frames attached to each body. This excludes bodies that display fluid, highly elastic, and plastic behavior.

In classical mechanics, the Laplace–Runge–Lenz (LRL) vector is a vector used chiefly to describe the shape and orientation of the orbit of one astronomical body around another, such as a binary star or a planet revolving around a star. For two bodies interacting by Newtonian gravity, the LRL vector is a constant of motion, meaning that it is the same no matter where it is calculated on the orbit; equivalently, the LRL vector is said to be conserved. More generally, the LRL vector is conserved in all problems in which two bodies interact by a central force that varies as the inverse square of the distance between them; such problems are called Kepler problems.

In classical mechanics, Euler's rotation equations are a vectorial quasilinear first-order ordinary differential equation describing the rotation of a rigid body, using a rotating reference frame with angular velocity ω whose axes are fixed to the body. Their general vector form is

In quantum physics, the spin–orbit interaction is a relativistic interaction of a particle's spin with its motion inside a potential. A key example of this phenomenon is the spin–orbit interaction leading to shifts in an electron's atomic energy levels, due to electromagnetic interaction between the electron's magnetic dipole, its orbital motion, and the electrostatic field of the positively charged nucleus. This phenomenon is detectable as a splitting of spectral lines, which can be thought of as a Zeeman effect product of two relativistic effects: the apparent magnetic field seen from the electron perspective and the magnetic moment of the electron associated with its intrinsic spin. A similar effect, due to the relationship between angular momentum and the strong nuclear force, occurs for protons and neutrons moving inside the nucleus, leading to a shift in their energy levels in the nucleus shell model. In the field of spintronics, spin–orbit effects for electrons in semiconductors and other materials are explored for technological applications. The spin–orbit interaction is at the origin of magnetocrystalline anisotropy and the spin Hall effect.

<span class="mw-page-title-main">Rotation around a fixed axis</span> Type of motion

Rotation around a fixed axis or axial rotation is a special case of rotational motion around an axis of rotation fixed, stationary, or static in three-dimensional space. This type of motion excludes the possibility of the instantaneous axis of rotation changing its orientation and cannot describe such phenomena as wobbling or precession. According to Euler's rotation theorem, simultaneous rotation along a number of stationary axes at the same time is impossible; if two rotations are forced at the same time, a new axis of rotation will result.

The theoretical and experimental justification for the Schrödinger equation motivates the discovery of the Schrödinger equation, the equation that describes the dynamics of nonrelativistic particles. The motivation uses photons, which are relativistic particles with dynamics described by Maxwell's equations, as an analogue for all types of particles.

The angular velocity tensor is a skew-symmetric matrix defined by:

In general relativity, Lense–Thirring precession or the Lense–Thirring effect is a relativistic correction to the precession of a gyroscope near a large rotating mass such as the Earth. It is a gravitomagnetic frame-dragging effect. It is a prediction of general relativity consisting of secular precessions of the longitude of the ascending node and the argument of pericenter of a test particle freely orbiting a central spinning mass endowed with angular momentum .

<span class="mw-page-title-main">Symmetry in quantum mechanics</span> Properties underlying modern physics

Symmetries in quantum mechanics describe features of spacetime and particles which are unchanged under some transformation, in the context of quantum mechanics, relativistic quantum mechanics and quantum field theory, and with applications in the mathematical formulation of the standard model and condensed matter physics. In general, symmetry in physics, invariance, and conservation laws, are fundamentally important constraints for formulating physical theories and models. In practice, they are powerful methods for solving problems and predicting what can happen. While conservation laws do not always give the answer to the problem directly, they form the correct constraints and the first steps to solving a multitude of problems. In application, understanding symmetries can also provide insights on the eigenstates that can be expected. For example, the existence of degenerate states can be inferred by the presence of non commuting symmetry operators or that the non degenerate states are also eigenvectors of symmetry operators.

In quantum mechanics, magnetic resonance is a resonant effect that can appear when a magnetic dipole is exposed to a static magnetic field and perturbed with another, oscillating electromagnetic field. Due to the static field, the dipole can assume a number of discrete energy eigenstates, depending on the value of its angular momentum (azimuthal) quantum number. The oscillating field can then make the dipole transit between its energy states with a certain probability and at a certain rate. The overall transition probability will depend on the field's frequency and the rate will depend on its amplitude. When the frequency of that field leads to the maximum possible transition probability between two states, a magnetic resonance has been achieved. In that case, the energy of the photons composing the oscillating field matches the energy difference between said states. If the dipole is tickled with a field oscillating far from resonance, it is unlikely to transition. That is analogous to other resonant effects, such as with the forced harmonic oscillator. The periodic transition between the different states is called Rabi cycle and the rate at which that happens is called Rabi frequency. The Rabi frequency should not be confused with the field's own frequency. Since many atomic nuclei species can behave as a magnetic dipole, this resonance technique is the basis of nuclear magnetic resonance, including nuclear magnetic resonance imaging and nuclear magnetic resonance spectroscopy.

<span class="mw-page-title-main">Accretion disk</span> Structure formed by diffuse material in orbital motion around a massive central body

An accretion disk is a structure formed by diffuse material in orbital motion around a massive central body. The central body is most frequently a star. Friction, uneven irradiance, magnetohydrodynamic effects, and other forces induce instabilities causing orbiting material in the disk to spiral inward toward the central body. Gravitational and frictional forces compress and raise the temperature of the material, causing the emission of electromagnetic radiation. The frequency range of that radiation depends on the central object's mass. Accretion disks of young stars and protostars radiate in the infrared; those around neutron stars and black holes in the X-ray part of the spectrum. The study of oscillation modes in accretion disks is referred to as diskoseismology.

References

  1. Velikhov, E. P. (1959), "Stability of an Ideally Conducting Liquid Flowing Between Cylinders Rotating in a Magnetic Field", J. Exptl. Theoret. Phys., vol. 36, no. 5, pp. 1398–1404, Bibcode:1959JETP....9..995V
  2. Chandrasekhar, S. (1960), "The stability of non-dissipative Couette flow in hydromagnetics", Proc. Natl. Acad. Sci., vol. 46, no. 2, pp. 253–257, Bibcode:1960PNAS...46..253C, doi: 10.1073/pnas.46.2.253 , PMC   222823 , PMID   16590616
  3. "Astrophys" . Retrieved 6 August 2023.
  4. Acheson, D. J.; Hide, R. (1973), "Hydromagnetics of Rotating Fluids", Reports on Progress in Physics, vol. 36, no. 2, pp. 159–221, Bibcode:1973RPPh...36..159A, doi:10.1088/0034-4885/36/2/002, S2CID   250881777
  5. Balbus, Steven A.; Hawley, John F. (1991), "A powerful local shear instability in weakly magnetized disks. I – Linear analysis. II – Nonlinear evolution", Astrophysical Journal, vol. 376, pp. 214–233, Bibcode:1991ApJ...376..214B, doi:10.1086/170270

Further reading