In-depth analysis of bar formation mechanisms of disk galaxies in halos of different concentrations

Using N-body simulations, this study reveals that bar formation mechanisms in disk galaxies vary significantly with halo concentration: highly concentrated halos favor swing-amplified multi-arm modes followed by particle trapping with minimal angular momentum transfer, lowly concentrated halos are dominated by linearly unstable modes triggering strong corotation resonance and rapid slowdown, while intermediate concentrations exhibit a combination of all these processes.

The Gist

Imagine a galaxy as a giant, spinning pizza dough made of stars. Sometimes, instead of being a perfect circle, this dough stretches out into a long, thick bar shape in the middle. Astronomers have known for a long time that these "bars" exist, but they've been arguing about how they get there.

This paper acts like a cosmic detective story. The author, T. Worrakitpoonpon, ran computer simulations to see how these bars form in galaxies living inside different types of "dark matter halos" (invisible, heavy clouds of gravity that hold the galaxy together). The key variable was how concentrated this invisible cloud was: was the heavy stuff packed tightly in the center, or was it spread out loosely?

Here is the breakdown of the three different stories the paper found, explained with everyday analogies:

1. The "Crowded Room" Scenario (High Concentration)

The Setup: Imagine the galaxy is in a very dense, crowded room where the gravity is tightly packed in the center.
What Happens:

• The Chaos First: When the simulation starts, the stars don't immediately form a bar. Instead, they get excited and form multiple spiral arms (like a 3-armed or 4-armed pinwheel). This is caused by a mechanism called "swing amplification," which is like a crowd doing "the wave" in a stadium—it amplifies small ripples into big spirals.
• The Bar is Stuck: These multi-arm spirals are so dominant that they actually block the bar from forming early on. It's like trying to build a straight wall in a room where everyone is dancing in a circle; the wall can't get started.
• The Slow Build: Eventually, the spiral arms die down. Only then does the bar start to grow. But it doesn't grow because of a magical resonance; it grows mechanically. Imagine a few stars getting "trapped" in the gravity of a tiny seed of a bar and slowly dragging other stars into line with them.
• The Result: A short, slow-growing bar. Because it didn't rely on a "resonance" (a perfect timing match between the bar and the stars), it doesn't slow down much, and it doesn't transfer much energy to the surrounding halo. It's a quiet, mechanical process.

2. The "Open Field" Scenario (Low Concentration)

The Setup: Now, imagine the galaxy is in a vast, open field where the gravity is spread out loosely.
What Happens:

• The Instant Bar: In this environment, the "swing amplification" (the wave) doesn't happen. Instead, the galaxy immediately picks the fastest-growing instability: a two-armed bar. It's like a straight line forming instantly because there's no crowd to distract it.
• The Resonance Engine: Because this bar forms so fast and is so strong right from the start, it immediately locks into a "resonance" with the stars. Think of this like a swing set: if you push the swing at exactly the right moment (resonance), it goes higher and higher.
• The Result: A strong, fast-growing bar that quickly grabs a lot of energy from the galaxy. This causes the bar to slow down significantly over time as it dumps that energy into the surrounding halo. It's a high-energy, explosive start.

3. The "Middle Ground" Scenario (Intermediate Concentration)

The Setup: This is the galaxy in a medium-density environment.
What Happens:

• The Mix: This scenario is a hybrid. It starts with some spiral arms (like the crowded room), but because the gravity isn't too tight, the two-armed bar mode can also grow early.
• The Teamwork: The spiral arms and the bar work together. The spiral arms help "seed" the bar, and then the bar takes over, locking into resonance.
• The Result: A bar that has features of both the other two types. It's stronger than the "Crowded Room" bar but not quite as explosive as the "Open Field" bar.

How to Tell Them Apart (The Detective Work)

The paper explains that if you look at a real galaxy, you can tell which story happened by looking at two things:

1. The Shape of the Surroundings:
   • Crowded Room: The bar is short, and the stars around it look very round and circular (like a smooth dough).
   • Open Field: The bar is long, and the stars around it look stretched out and oval-shaped (like the dough was pulled).
2. The "Speed" of the Stars:
   • By analyzing how the stars move (kinematics), astronomers can see if the bar is "resonating" (slowing down fast) or just "trapping" stars (slowing down slowly).

The "Shearing vs. Rigid" Ratio

Finally, the author introduces a simple math tool (called $\Xi$) to predict which scenario will happen.

• Think of the galaxy as a spinning record.
• Shearing: The tendency for the outer parts to spin slower than the inner parts, which tries to stretch things out (like pulling taffy).
• Rigid Motion: The tendency for the whole thing to spin together like a solid wheel.
• The Verdict: If the "stretching" force (shearing) is stronger than the "spinning together" force, you get the Crowded Room scenario (spirals first, then a slow bar). If the "spinning together" force wins, you get the Open Field scenario (instant bar).

Summary

The paper concludes that there isn't just one way a galaxy bar forms. It depends entirely on how heavy and concentrated the invisible gravity cloud is around the galaxy.

• Tight Gravity = Spirals first, then a slow, mechanical bar.
• Loose Gravity = Instant, fast, resonant bar.
• Medium Gravity = A mix of both.

This helps astronomers understand why some galaxies look different from others and gives them new tools to figure out the history of the galaxies they observe in the night sky.

Technical Summary

Technical Summary: In-depth analysis of bar formation mechanisms of disk galaxies in halos of different concentrations

Problem Statement
While the prevalence of barred structures in disk galaxies is well-documented, the detailed physical processes governing their formation across different parameter ranges remain incompletely understood. Previous research has largely bifurcated bar formation into two paradigms: the internal (secular) scheme driven by global unstable non-axisymmetric modes, and the external scheme driven by tidal interactions. Within the secular framework, the interplay between three key theoretical mechanisms—linear instability, swing amplification, and orbital (particle) trapping—has not been unified. Specifically, it is unclear where and when each mechanism dominates, and how the central mass concentration of the host dark matter halo influences the transition between these processes.

Methodology
The authors employ N-body simulations using the GADGET-2 code to investigate bar formation in disk-halo systems with varying halo concentrations.

• System Setup: The models consist of a live Hernquist halo and a three-dimensional exponential disk. The disk mass is fixed at $9.3 \times 10^9 M_\odot$, while the halo mass is $25 M_d$. The central mass concentration is varied by changing the halo scale radius ($r_h$) from 35 kpc to 75 kpc.
• Initial Conditions: The initial radial Toomre $Q$ parameter is set to a minimum ($Q_{min}$) of 1.1 to ensure bar instability occurs within the simulated timescale. The disk and halo are populated with $2 \times 10^6$ and $3 \times 10^6$ particles, respectively.
• Analysis Techniques: The study utilizes morphological analysis (isodensity contours, radial Fourier amplitudes $\tilde{A}_m(R)$), kinematical analysis (Fourier spectrograms $\omega_F$ of $m=2$ and $m=3$ modes), and dynamical analysis (probability distribution functions of trapped particle angular frequencies $P(\omega_i)$, angular momentum transfer $\Delta L_z$, and pattern speed evolution $\omega_b$). A new dimensionless parameter, the shearing-to-rigid ratio ($\Xi$), is introduced to quantify the competition between shearing forces and rigid mode stability.

Key Results
The simulations reveal three distinct bar formation scenarios dependent on the halo concentration (represented by $r_h$):

1. High Concentration (Low $r_h$, e.g., R35): Dynamics-Based Bar Formation

   • Mechanism: In highly concentrated halos, strong shearing promotes the rapid growth of multi-arm modes ($m=3, 4$) via swing amplification in the early stages ($<3$ Gyr). These modes dominate and suppress the growth of the bi-symmetric ($m=2$) linear bar modes.
   • Evolution: Once the multi-arm modes decay, the bar grows mechanically through particle trapping (orbital resonance with the bar potential) rather than through resonance-driven linear instability.
   • Characteristics: The corotation resonance is absent throughout the evolution. Consequently, there is minimal disk-halo angular momentum transfer, and the bar pattern speed decreases only modestly. The resulting bar is relatively short, embedded in a dense, circular background, with vanishing bi-symmetric Fourier amplitude beyond the bar end.

2. Low Concentration (High $r_h$, e.g., R75): Linear-Mode-Based Bar Formation

   • Mechanism: In dilute halos, shearing is mild. The fastest-growing modes are the linearly unstable, bi-symmetric ($m=2$) modes of a single uniform frequency. Swing amplification plays no significant role.
   • Evolution: These linear modes trigger the corotation resonance immediately. The resonance is maintained throughout the evolution, driving rapid bar growth and significant angular momentum transfer to the halo.
   • Characteristics: The bar pattern speed slows down rapidly. The bar is surrounded by a non-vanishing bi-symmetric Fourier amplitude extending beyond the bar, creating an elliptical environment.

3. Intermediate Concentration (Medium $r_h$, e.g., R50): Mixed Mechanism

   • Mechanism: This regime exhibits a combination of the previous two. Swing-amplified two-armed modes arise first, followed by the triggering of resonance-supported linear bar modes. Particle trapping also contributes but is less dominant than in the high-concentration case.
   • Characteristics: The resulting bar shows mixed morphological and kinematical features, being stronger than the high-concentration case but sharing some structural similarities with the low-concentration case.

Key Contributions and Significance

• Unified Framework: The paper provides a unified picture of bar formation, demonstrating that the dominant mechanism is not static but depends critically on the halo's central mass concentration. It clarifies that swing amplification and linear instability are distinct entities that can act counter-actively (high concentration), constructively (intermediate), or independently (low concentration).
• Role of Resonance: The study highlights that a bar can form and grow without the support of corotation resonance (the "dynamics-based" scenario), challenging the assumption that resonance is a requisite for all bar evolution. This explains the existence of short, fast-rotating bars (nuclear bars) observed in the universe.
• Diagnostic Tools: The authors propose specific observational diagnostics to distinguish between these formation histories:
  • The ratio of ellipticities of outer isodensity contours ($E = e_{0.2}/e_{0.4}$).
  • The radial profile of the bar amplitude $\tilde{A}_2(R)$ (specifically the presence or absence of bi-symmetric modes beyond the bar).
  • Fourier spectrograms to identify the presence or absence of resonance and the nature of the dominant modes.
• Shearing-to-Rigid Ratio ($\Xi$): The introduction of the parameter $\Xi$ offers a quantitative way to describe the transition between formation regimes, linking the Coriolis force (shearing) against centrifugal forces (rigid modes).

Limitations and Modesty
The authors acknowledge that their analysis is restricted to the $Q_{min} \approx 1$ regime and bulge-less models. They note that for higher $Q$ values or systems with massive bulges, pressure forces and cavity narrowing might alter these tracks. Furthermore, while the parameter $\Xi$ successfully describes the transition in their models, pinpointing a precise critical value for $\Xi$ that delineates the regimes is currently impractical due to sensitivity to the choice of pattern speed. The study focuses on the intrinsic dynamics of isolated disks, leaving the complex interplay with external tidal forces for future unified investigations.