The Effect of Gravitational Stratification on Kink Oscillations in Curved Coronal Loops

Titolo: L'effetto della stratificazione gravitazionale sulle oscillazioni di kink in loop coronali curvi

1. Problema e Contesto

La sismologia coronale utilizza le onde MHD (magnetoidrodinamiche), in particolare le onde di kink, per inferire indirettamente parametri fisici del Sole, come il campo magnetico coronale. Tuttavia, i modelli teorici attuali si basano spesso su cilindri magnetici canonici che trascurano due fattori cruciali presenti nella realtà: la curvatura dei loop coronali e la stratificazione gravitazionale (variazione della densità con l'altezza).
Sebbene studi precedenti abbiano esaminato separatamente l'effetto della curvatura o della stratificazione, non è ancora chiaro come la combinazione di questi due fattori influenzi le frequenze delle oscillazioni di kink, specialmente distinguendo tra le due polarizzazioni osservate: verticale e orizzontale. La mancanza di un modello che integri simultaneamente gravità e curvatura limita l'accuratezza delle inversioni sismologiche.

2. Metodologia

Gli autori hanno condotto simulazioni numeriche tridimensionali (3D) all'interno del quadro della MHD ideale utilizzando il codice PLUTO.

• Configurazione del Modello: È stato creato un loop coronale curvo (geometria toroidale approssimata) immerso in un campo magnetico potenziale. Il sistema include una stratificazione gravitazionale realistica, dove la gravità diminuisce con l'altezza secondo la legge di Newton.
• Condizioni Iniziali: È stato costruito un loop con densità aumentata rispetto al mezzo circostante (rapporto di densità 2:1 a $z=0$) e temperatura interna di 1 MK. Il sistema è stato lasciato rilassare per 2000 secondi fino a raggiungere un equilibrio magnetostatico stabile.
• Perturbazioni: Sono state eccitate diverse modalità di oscillazione tramite perturbazioni di velocità iniziali:
  • Modo fondamentale: Perturbazioni sia verticali ($\hat{e}_v$) che orizzontali ($\hat{e}_h$).
  • Primi armonici (overtones): Perturbazioni con profilo $\sin(\theta)\cos(\theta)$ per eccitare i primi overtones.
• Analisi: Le frequenze di oscillazione sono state determinate tracciando lo spostamento Lagrangiano di particelle fluide (tracer passivi) posizionate all'apice del loop e a un quarto della sua lunghezza. I risultati numerici sono stati confrontati con l'approssimazione WKB (Wentzel-Kramers-Brillouin) e con la frequenza di Alfvén locale.

3. Risultati Chiave

Lo studio ha rivelato differenze sostanziali nel comportamento delle oscillazioni a seconda della polarizzazione e dell'armonica:

• Modo Fondamentale (Polarizzazione Verticale):

  • La frequenza di oscillazione misurata si discosta dall'approssimazione WKB di circa il 18%.
  • Tuttavia, la frequenza misurata corrisponde quasi perfettamente alla frequenza di Alfvén locale vicino all'apice del loop.
  • La deviazione dal WKB è attribuita alla non uniformità assiale indotta dalla gravità, che crea gradienti più ripidi nella frequenza di Alfvén e rende l'approssimazione WKB meno affidabile per questa polarizzazione.

• Modo Fondamentale (Polarizzazione Orizzontale):

  • La frequenza mostra una deviazione molto minore dall'approssimazione WKB, pari a circa il 7%.
  • La frequenza misurata corrisponde alla frequenza di Alfvén locale situata circa a un quarto della lunghezza del loop.
  • In questo caso, l'oscillazione avviene in una direzione dove la frequenza di Alfvén è più uniforme, permettendo al WKB di funzionare bene.

• Primi Armonici (Overtones):

  • Per entrambi i tipi di polarizzazione (verticale e orizzontale), le frequenze dei primi overtones sono ben descritte dall'approssimazione WKB (deviazione massima ~4%).
  • Le frequenze di entrambi gli overtones corrispondono alla frequenza di Alfvén locale vicino a un quarto del loop.

• Rapporto dei Periodi ($P_1/P_2$):

  • Il rapporto tra il periodo fondamentale e il primo overtone ($P_{fun}/P_{1st}$) è stato analizzato per diagnosticare la scala di altezza della densità.
  • Per la polarizzazione verticale, il rapporto supera 2 (valore misurato: 2.35), un fenomeno attribuito all'espansione della sezione trasversale del loop percepita solo da questa polarizzazione a causa della curvatura e della gravità.
  • Per la polarizzazione orizzontale, il rapporto è 1.94.
  • L'analisi mostra che l'eterogeneità del campo magnetico deve essere inclusa nei modelli di inversione per stimare correttamente la scala di altezza della densità.

4. Contributi e Significato Scientifico

Questo lavoro rappresenta un passo avanti significativo per la sismologia coronale spaziale:

1. Validazione di Nuovi Indicatori Sismologici: Lo studio dimostra che la frequenza di Alfvén locale in punti specifici del loop (apice per la polarizzazione verticale, un quarto per quella orizzontale) è un predittore più accurato della frequenza di oscillazione rispetto all'approssimazione WKB globale, specialmente in presenza di forte stratificazione gravitazionale.
2. Diagnostica del Campo Magnetico Spazialmente Dipendente: I risultati suggeriscono che, osservando simultaneamente le oscillazioni di kink con diverse polarizzazioni, è possibile sondare la forza del campo magnetico in diverse posizioni lungo il loop (ad esempio, all'apice e a un quarto della lunghezza), purché la densità locale sia nota. Questo permette una mappatura più dettagliata del campo magnetico coronale.
3. Raffinamento dei Modelli di Inversione: Lo studio evidenzia i limiti dell'uso esclusivo della teoria WKB per le polarizzazioni verticali in loop stratificati e sottolinea la necessità di includere l'eterogeneità del campo magnetico quando si utilizzano i rapporti dei periodi per diagnosticare la stratificazione della densità.
4. Implicazioni Osservative: Con l'avvento di osservazioni ad alta risoluzione e multi-angolo (es. da SDO, Solar Orbiter), la capacità di distinguere le polarizzazioni e i periodi associati diventa cruciale per applicare questi nuovi metodi di inversione alla realtà osservativa.

In conclusione, il lavoro fornisce un quadro teorico e numerico robusto per interpretare le oscillazioni di kink in loop reali, migliorando la nostra capacità di diagnosticare le condizioni fisiche della corona solare.