1 / 27

NLTE polarized lines and 3D structure of magnetic fields

NLTE polarized lines and 3D structure of magnetic fields. P.Mein, N.Mein, M.Faurobert, V.Bommier, J-M.Malherbe, G.Aulanier. M agnetic fields cross canopy regions, not easily investigated by extrapolations, between photosphere and chromosphere.

talor
Télécharger la présentation

NLTE polarized lines and 3D structure of magnetic fields

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. NLTE polarized lines and 3D structure of magnetic fields P.Mein, N.Mein, M.Faurobert, V.Bommier, J-M.Malherbe, G.Aulanier Magnetic fields cross canopy regions, not easily investigated by extrapolations, between photosphere and chromosphere. Full knowledge of the 3D structure implies diagnostics extracted from strong NLTE lines. Fortunately, the domain of ‘’weak field ’’ approximation is more extended for such lines (smaller Lande factor, broad lines). The data analysed below are obtained with THEMIS / MSDP and MTR in 589.6 NaI (D1) 610.27 CaI 630.2 FeI (for comparison)

  2. 1) D1 line and facular magnetic flux tubes Problems of filling factor, vertical gradients, MHD models Simulation of line profiles MULTI code with field free assumption, 1D model Instrumental profile included • - Quiet Sun = VAL3C model • - Circular polarization: I -V profile • Solid line: flux tube, dashed: quiet • Bisector for = +/-8, 16, 24, 32 pm • Weak field assumption B//

  3. 1 - 2D model flux tube compensating horizontal components of Lorentz forces Magnetic field Departures from equilibrium Formation altitudes of B// for  = +/- 8, 16, 24, 32 pm Bz(0,z) ~ exp(-z/h)Bz(x,z) ~ cos2(x/4d(z)) d(z) by constant flux Bx(x,z) by zero divergence P(x,z) compensates Lorentz horiz. comp. Vertical accelerations exceed solar gravity at high levels

  4. Simulation Smoothing by seeing effects convol cos2(x/4s) s=400 km No smoothing by seeing effects wings wings core core B// from tube center at  = 8, 16, 24, 32 pm Points at half maximum values(crosses) are in the same order as tube widths at corresponding formation altitudes Seeing effects ~ filling factor effects hide vertical magnetic gradients attube center

  5. Filling factors and slope-ratios of profiles flux-tube/quiet-sun I - V Stokes Vobs = f Stokes Vtube Tube Quiet Sun Zeeman shift of I -V profile: Zobs (dI/dl)obs = f Ztube (dI/dl)tube If f << 1, from core to wings Zobs = f Ztube (dI/dl)tube / (dI/dl)QS dI/dl Tube QS Decrease of observed B// in the wings Different models for tube and quiet sun !

  6. 2 - Model flux tube closer to magneto-static equilibrium Magnetic field Departures from equilibrium Formation altitudes of B//  = +/- 8, 16, 24, 32 pm Bz(0,z)2/20 ~ Pquiet(z) Bz(x,z) ~cos2(x/4d(z)) d(z) by constant flux Bx(x,z) by zero divergence P(x,z) = Pquiet – Bz(x,z)2/20 Departures from equilibrium never exceed solar gravity

  7. Observation Simulation Average of 6 magnetic structures Faculae near disk center (N17, E18) With seeing effects s=500 km core core wings wings Sections for = 8, 16, 24, 32 pm Qualitative agreement only: - tube thinner in line wings - apparent B// smaller in line wings (seeing effects) But impossible to increase the magnetic field and/or the width of the tube without excessive departures from equilibrium

  8. 3 - Conglomerate of flux tubes Magnetic field Departures from equilibrium

  9. Observation Simulation Seeing effects s = 700 km Better qualitative agreement (tube width) But magnetic field still too low Coronal magnetic field outside the structure? MHD models, including temperature and velocity fluctuations…? P. Mein, N. Mein,M. Faurobert, G. Aulanier and J-M. Malherbe, A&A 463, 727 (2007)

  10. 2) Fast vector magnetic maps with THEMIS/MSDP UNNOFIT inversions NLTE line 610.27 CaI + 630.2 FeI • Examples of fast MSDP vector magnetic maps and comparison with MTR results - How to reconcile high speed and high spectral resolution by compromise withspatial resolution in MSDP data reduction - Capabilities expected from new THEMIS set-up (32) and EST project (40) • Departures between 610.3 CaI and 630.2 FeI maps Gradients along LOS? sensitivity of lines? filling factor effects?

  11. Example of MSDP image (Meudon Solar Tower 2007, courtesy G. Molodij): In each channel, x and  vary simultaneously along the horizontal direction

  12. Compromise spatial resol / spectral resol interpolation in x,plane A , D 80 mA B , C --> E 40 mA cubic interpol --> F,G 20 mA

  13. Profile deduced from 16 MSDP channels + interpolation x,plane 610.27 CaI

  14. THEMIS / MSDP 2006 610.3 CaI UNNOFIT inversion THEMIS / MTR 2006 630.2 FeI 160’’ 70’’ 70’’ Aug 18, NOAA 904 S13 , W35 120’’

  15. Scatter plots Ca (MSDP) / Fe (MTR) f B// f Bt

  16. THEMIS/MSDP 2007 UNNOFIT inversion 630.2 FeI 610.3 CaI 160’’ 120’’ 120’’ June 11, NOAA 10960 S05, W52

  17. 610.3 CaI I Q/I U/I V/I

  18. THEMISMSDP

  19. f Bx f By Similar Bx and By similar f angles

  20. THEMIS MSDP • Bt 6103 < Bt 6302 • - Gradients along line of sight ? • B t more sensitive than B// to line center, 6103 saturated NLTE line? • stray-light effects? • instrumental profile not included? • filling factor effects? • further simulations needed ….. • comparisons with MTR data (not yet reduced) • Possible improvements: • - Include instrumental profile • - set-up 32 channels (2 cameras = effective increase of potential well) • better size of 6302 filter !

  21. Scanning speed for targets 100’’x160’’ 9 mn

  22. 3) Problems and plans: Gradients of B along LOS from NaD1, 610.3 Ca, 630.2 Fe, … Weak field approximation Disk center, no rotation of B along LOS: Stokes U = 0 1- Simple case: LTE, Milne Eddington, B vector and f independent of z V(l) ~ f Bl dI/dl Q(l) ~ f Bt2 d2I/dl2 For weak fields, line profile inversions provide only 2 quantities, f Bl and f Bt2

  23. Bt* f Bt* f1/2 Below a given level of c2 the range of possible solutions is larger for (Btransverse * f) than for (Btransverse * f ½) ?

  24. 2–NLTE, B function of z, f = 1, given solar model(parts of spots?) > Computation of response functions by MULTI code V(l) = S Bl(z) R(l,z) dz Q(l) = S Bt2(z) R’(l,z) dz ? > Formation altitudes: barycenters of response functions

  25. > Vertical gradients: Circular polarization: Bl (0), dBl/dz Bl = a + bz V(l)=SR(l,z) (a+bz) dz Instead of using individual points of the bisector, integrations along line profile to optmize signal/noise ratio. Choose functions with different weights at line-center and wings: Examples: w1(l) = +1 and -1 around line center, 0 elsewhere w2(l) = +1 and -1 in line wings, 0 elsewhere S1= SV(l) w1(l) dl S2= SV(l) w2(l) dl a, b Linear polarization: Bt ??? Bt2= a’+ b’z Q(l)= SR’(l,z) (a’+b’z) dz Bt2, dBt2/dz ??? > Application: comparisons between gradients from full profile of 1 line and 2 different lines

  26. 3 – NLTE, B and f functions of z, given solar model(flux tubes?) Example: flux tubes, NaD1 line (section 1) Bl = a + bz f = a+ bz V(l) = S R(l,z) (a+bz) (a+bz) dz Both unknown quantities b and b are present in the coefficient of z Impossible to determine separatelyf and Bl In particular, when flux tubes are not spatially resolved, the assumption of constant flux f Bl implies that gradients of f and Bl are compensated (b/a ~ - b/a) MHD model necessary in case of weak fields (see section 1). 4 – Possible extension of UNNOFIT to NLTE lines close to LTE ? Example: Analysis of depatures between UNNOFIT results and parameters used for synthetic profiles in case of 610.3 Ca …

More Related