Relativistic Jets & Oscillations in Magnetized Accretion Flows

Despite decades of research, four key mysteries about magnetic reconnection remain unsolved: (1) What produces the dissipation that occurs in the diffusion region? (2) What determines the aspect ratio of the diffusion region (whether it is fat or thin), which controls the reconnection rate? (3) What is the spark that initiates reconnection when magnetic fields build up over time? (4) What mechanism converts magnetic energy to particle energy? The MMS mission aims to answer all four of these fundamental questions about the microphysics of magnetic reconnection.

Plasma is highly electrified gas composed of ions and electrons, whose motion is driven by magnetic and electric fields. Magnetic reconnection is a fundamental plasma process where magnetic topology rearranges, converting magnetic energy into kinetic energy, thermal energy, and particle acceleration. Opposing magnetic field lines connect at X-points, causing particles to change direction and heat up. This process occurs universally in space environments including the Sun, Jupiter, laboratory settings, the heliosphere edge, and black hole accretion discs.

Understanding magnetic reconnection requires a multiscale and multiphysics approach connecting microscopic to largest scales. Many scale lengths enter into play, and answering reconnection questions requires looking at different descriptions: MHD versus kinetic, relativistic effects, geometry considerations (strong guide field vs double periodic), 2D versus 3D, boundary conditions, and the mutual influence between reconnection and turbulence. Spontaneous and driven models must be integrated to understand how necessary conditions for reconnection are self-consistently created.

Magnetic reconnection is a change in the topology of the magnetic field. It requires non-ideal effects that make the plasma non-perfectly conducting and break the frozen-in-field condition. These non-ideal terms can include resistivity (from particle collisions), finite electron mass (causing delayed response to electric field fluctuations), or particle gyroradius effects. For reconnection to occur, one of these effects must become important at small scales locally in a region called the diffusion region, where the magnetic field direction reverses and a strong out-of-plane current exists.

Magnetic reconnection violates the frozen field theorem when anti-parallel magnetic flux tubes collide, annihilate, and merge, jetting out perpendicular to inflow. At Earth, southward IMF enables dayside reconnection. The process involves three regions: inflow (flux tubes convected in), diffusion (where plasma decouples from field), and outflow (exhaust jet). Due to mass differences, ions and electrons decouple differently: ions have larger gyro radii and cross the separatrix further out, creating larger ion diffusion regions (~ion inertial lengths) and smaller electron diffusion regions (~electron skin depths). This decoupling creates Hall magnetic fields and current loops essential for reconnection dynamics.

Quasi-periodic oscillations (QPOs) are flickering patterns in X-ray light from astronomical objects like neutron stars and black holes, occurring near the inner edges of accretion disks; QPOs appear as broader peaks in power spectra compared to sharp periodic pulsations, and their frequencies correlate with the mass of the central compact object, allowing astronomers to measure black hole masses by analyzing how the 'congestion zone' of hot gas near the event horizon oscillates at different rates depending on the black hole's mass.

QPOs originate from shock oscillations in the accretion flow: (1) Shocks form when gravitational and rotational forces balance, causing matter to virtually stop; (2) When cooling time scales become comparable to compressional heating time scales, oscillations develop even if Rankine-Hugoniot conditions predict steady solutions; (3) QPO frequencies match theoretical predictions when cooling-to-heating time scale ratio falls within 0.5-1.5; (4) No QPOs are observed in soft states when the corona nearly disappears; (5) This provides a direct link between temporal variability and underlying accretion physics, allowing constraints on cooling rates and flow dynamics.

Quasi-periodic oscillations (QPOs) appear as broad peaks in power spectra with quality factors Q = f₀/Δf of 1-10. They arise from damped driven oscillators where turbulent driving combines with damping to produce characteristic frequencies. In Cygnus X-1, QPO frequencies correlate smoothly with spectral properties, suggesting they originate from the inner accretion disk region (20-70 Schwarzschild radii). The disk may act as a noise source, while the corona/jet represents the response, creating a damped driven oscillator system. Multiple independent processes can coexist, producing complex power spectra requiring careful modeling.

Quasi-periodic oscillations (QPOs) are periodic variations in X-ray brightness observed around black holes. In this case, scientists discovered extremely precise QPOs with periods starting at 18 minutes and decreasing to just 7.1 minutes within two years. These short-period oscillations (less than 10 minutes) had never been observed around any supermassive black hole before, representing an entirely new astrophysical phenomenon that required new theoretical explanations.

Quasi-periodic oscillations (QPOs) are nearly periodic signals in X-ray binaries that reveal fundamental physics of compact objects. Using the damped harmonic oscillator model, pure sine waves produce sharp frequency peaks while damping broadens them proportionally to dissipation strength. QPOs are classified by frequency: high-frequency QPOs at hundreds to thousands of hertz arise from hot blobs in Keplerian orbits at the inner disk edge. Low-frequency QPOs at 0.1-10 Hz may result from Lense-Thirring precession of misaligned hot inner flows around spinning black holes. The Lense-Thirring mechanism involves precession of the entire hot inner flow region, where magnetic field breaks allow coherent precession around the black hole's spin axis. QPO evolution correlates with spectral state transitions: hard states show high variability without QPOs, intermediate states show emerging QPOs with frequency wandering, and soft states show decreased variability with no QPOs. This systematic evolution indicates QPO mechanisms require both corona and disk components simultaneously.

This comprehensive section covers the theoretical and numerical foundations of computational astrophysics. It begins with astrophysical plasmas, which constitute over 99.9% of universal matter, and explains how plasma conditions vary dramatically across objects (length scales from 10^10 to 10^24 cm, density from stellar interiors to extra-galactic medium). The major challenge is scale separation—no single numerical method covers from kinetic to fluid scales. Different theoretical descriptions apply at different scales: kinetic description (particle-in-cell) for small scales, hybrid kinetic models, two-fluid description, and MHD for largest scales. MHD is derived from two-fluid equations by assuming low frequency, large length scales, and ignoring electron effects. The section then transitions to numerical methods, explaining the finite volume method (transforming continuous equations into discrete relations using Gauss's theorem, guaranteeing conservation), and the Riemann solver concept (Godunov, 1861) for hyperbolic fluxes. In Euler equations, this generates left rarefaction, middle contact, and right shock waves. In MHD, the structure is more complex with fast waves (Alfvén waves), slow waves, and entropy waves.

Simulating jets requires modeling gravitational fields and hydrodynamics including Kelvin-Helmholtz and Rayleigh-Taylor instabilities. Codes like Enzo, Gizmo, Gadget, and Arepo handle these calculations. Resolution elements are typically 10 light-years across, much larger than our solar system. Adding more physics dramatically increases computational cost, requiring supercomputers. Simulations start simple and build complexity incrementally.

Computational astrophysics emerged as computers enabled simulations of processes too complex for manual calculation. Early applications included modeling star formation, evolution, and supernovae. The Atlas of Peculiar Galaxies cataloged unusual shapes that puzzled astronomers until simulations revealed they result from galaxy collisions. This computational approach resolved major mysteries by showing that peculiar galaxies are normal galaxies distorted by gravitational interactions, demonstrating how simulations extend human understanding beyond direct observation.

MHD instabilities arise from complex frequencies causing exponential growth or damping. Key examples include gravitational instability (self-gravity-driven collapse), thermal instability (fingering at interfaces), Rayleigh-Taylor instability (shear-driven mixing), Kelvin-Helmholtz instability (relative motion between fluids), and magneto-rotational instability (angular momentum transport in accretion discs). These instabilities drive structure formation in astrophysical systems. Three fundamental concepts govern numerical method quality: convergence (truncation error decreases with resolution), consistency (discrete equations converge to correct PDE solution), and stability (roundoff errors remain bounded). Roundoff error from finite precision arithmetic is generally uncontrollable, while truncation error from approximating differential equations can be minimized. Higher-order schemes reduce truncation error but may sacrifice shock-capturing properties. First-order schemes are unavoidable for discontinuities, so high-order methods must maintain good shock handling.

Computational astrophysics uses numerical simulations to study cosmic phenomena that cannot be directly observed due to their enormous scale and extremely long time scales, such as star formation and supernova explosions. The Center for Computational Astrophysics (CfCA) at the National Astronomical Observatory of Japan provides researchers with specialized computing resources including supercomputers, computing servers, GRAPE gravitational accelerators, and GPU servers to perform these simulations. These computational tools enable scientists to model complex astronomical processes like galaxy formation, stellar evolution, and supernova diversity by simulating gravitational interactions and other physical phenomena that would be impossible to study through direct observation alone.