![Page 1: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/1.jpg)
1
Statistical determination of chromospheric density structure
using RHESSI flares
Pascal Saint-HilaireSpace Sciences Lab, UC Berkeley
RHESSI WorkshopGenova, 2009/09/02
![Page 2: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/2.jpg)
2
Sato, 2005: (similar to Matsushita et al., 1992)
(L:14–23,M1:23–33, M2:33–53, H:53–93 keV)
L-M1: 460 ± 40 kmL–M2: 880 ± 90 kmL-H: 930 ± 150 km
Aschwanden, 2002:
![Page 3: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/3.jpg)
3
DATA:
• Visibilities from 830 RHESSI flares (took all flares “imageable above 25 keV”, as per Jim McTiernan’s “official RHESSI flare list”)
• SC 3-9, energy bins from 4-100 keV• Found position at different energies by forward
fitting the visibilities: his_vis_fwdfit.pro• No big difference if taking only events >50 keV
or other SC combinations. Better with longer accumulations (1-3 mins)
![Page 4: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/4.jpg)
4
Average height differences:
Average altitude of emission at energy hν, compared to altitude of 35 keV emission.
METHOD:• Assume footpoints emission at different energies extend radially• Choose two energies• Plot altitude differences r vs. altitude r for all flares (projected)• Obtain slope s, h = s RS
• Repeat for any position-energy pairs of data
r
r [RS] Thermal contamination
![Page 5: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/5.jpg)
5
dz
d
Kdz
dNzn
)(
For injected power-law of electrons in steeply increasing density gradients: Most of the emission at energy is emitted where the electrons have traveled a column density N≈ 2/2K (Brown, 1972; Brown et al., 2002; Xu et al., 2008)”. There is a dependence on .
90 keV
40 keV
10 keV
![Page 6: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/6.jpg)
6
CICM (H)VAL-C
Aschwanden et al (2002)
Href (35 keV)≈ 1.8 Mm
above photosphere
Comparisons with a few models:
•Arbitrary shift in the X-direction. •Last points to the right : ~15-20 keV: contamination by thermal plasma•Points on the left: energies > 80 keV: fewer statistics•Assumed uniform target ionization! High-energy emissions should actually yield ~2.8 times higher densities…It’s a fitting game…
![Page 7: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/7.jpg)
7
Forward-fitting attempts:
1. Double exponentials (0% ionized at low altitude, 100% ionized at higher altitude)
2. Assuming CICM (Caltech Irreference Chromospheric Model, Ewel et al., 1993) is “the truth” (!), trying to determine at which height the ionization level (suddenly) changes…
![Page 8: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/8.jpg)
8
Fitting double exponentials:
CICM (H)VAL-C (H)
Fit
Gabriel Corona, models A & D
HIGHLOW HhHIGH
HhLOW enen
hn //
8.2)(
hREF = 0.69 +/- 0.43 @ 35 keVhacc = 8.6 +/- 4.3; nacc = (1.3 +/- 2.2) x1010 cm-3 (!!)
Scale heights: HHIGH= 1.5 +/- 1.8 Mm; HLOW= 0.11 +/- 0.06 Mm
Additional pt
![Page 9: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/9.jpg)
9
Fitting step ionization:
Fit
VAL-C (H)
CICM (H)
h* = 1.38 +/- 0.09 Mm (at E*≈45 keV) (height of ionization change)
hacc = 27.3 +/- 61.3 Mm (acceleration height)
![Page 10: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/10.jpg)
10
Conclusions:
• Used 800+ RHESSI flares to statistically study chromospheric densities
• Double-exponential chromospheric fittings:– h35 keV = 0.7 +/- 0.4 Mm; hacc= 8.6 +/- 4.3 Mm– HLOW= 0.11 +/- 0.06 Mm, HHIGH ≈ 1 Mm
• Step ionization fitting: – h*=1.38 +/- 0.09 Mm
• Difficult to trust results from models with more than a few parameters: error bars still rather large
![Page 11: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/11.jpg)
11
![Page 12: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/12.jpg)
12
![Page 13: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/13.jpg)
13
Type III radio burst density (with the Nancay Radioheliograph, 1998):
![Page 14: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/14.jpg)
14
Density structure using dm radio bursts (preliminary):
•Same method•12’000 radio bursts over 10 years•164-432 MHz•νp ≈ 9000 ne
1/2
Forward-fitting an exponential atmosphere: Scale height H = 95.3 +/- 1.2 Mm
![Page 15: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/15.jpg)
15
![Page 16: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/16.jpg)
16
![Page 17: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/17.jpg)
17
![Page 18: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/18.jpg)
18
![Page 19: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/19.jpg)
19
Statistical approach: data from 830 flares…(inspired from Matsushita et al., 1992)
cossin)( 11 rhX
cossin12 hXXX
cossin)( 22 rhX
SR
h
r
h
X
X
11
•(X1,Y1) position of source at hν1 (e.g. arcsecs from sun center)•(X2,Y2) position of source at hν2
•Assume both sources are radially distributed above solar photosphere
Slope yields height difference.Tilt of magnetic field (and other) effects should average out!!
![Page 20: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/20.jpg)
20
1. Step ionization height(assuming CICM nH model)
“Ionization changes from 100% to 0% at an altitude of ~1.34 Mm”
CICM nH
CICM ne
FAL-P
Gabriel
![Page 21: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/21.jpg)
21
1. Step ionization height(100 MC runs)
h* = 1.38 +/- 0.09 Mm (height of step ionization change, E*≈45 keV)hacc = 27.3 +/- 61.3 Mm, i.e. is unreliable (height of acceleration region)
![Page 22: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/22.jpg)
22
2. Double exponential atmosphere(low: neutral; high: 100% ionized)
hREF = 0.8 @ 35 keVhacc = 6.2 Mm ; nacc = 2.9 x 1010 cm-3
Scale heights: HHIGH= 1.0 Mm; HLOW=126 km (Kontar et al., 2008: 140 +/- 30)
![Page 23: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/23.jpg)
23
2. Double exponential atmosphere(low: neutral; high: 100% ionized)
hREF = 0.69 +/- 0.43 @ 35 keVhacc = 8.6 +/- 4.3; nacc = (1.3 +/- 2.2) x1010 cm-3 (!!)
Scale heights: HHIGH= 1.5 +/- 1.8 Mm; HLOW= 0.11 +/- 0.06 Mm…not really enough reliable low-E points to determine well high altitude exp characteristics
![Page 24: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/24.jpg)
24
Adding low-energy points leads to problems… (as expected).
(nice fit, though…)
![Page 25: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/25.jpg)
25
“For a steep (power-law) accelerated electron spectrum propagating through an increasingly denser plasma, there is a peak in emission at energy hν= located at column density N ≈ 2/K. N is the column density required to stop electrons of initial energy .”(K actually varies with ionization level and with : see e.g. Kontar et al., 2002, Xu et al., 2008)
*1 keV
3 keV
10 keV
20 keV
50 keV
![Page 26: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/26.jpg)
26
Aschwanden et al., 2002:
![Page 27: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/27.jpg)
27
![Page 28: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/28.jpg)
28
Average altitude of emission at energy hν, compared to altitude emission at 35 keV.
Contamination from thermal emission in loops.(probably up to 30 keV…)
dz
d
Kdz
dNzn
)(
“Direct derivation or inversion”(assumes uniform ionization, here)
![Page 29: 1 Statistical determination of chromospheric density structure using RHESSI flares Pascal Saint-Hilaire Space Sciences Lab, UC Berkeley RHESSI Workshop](https://reader034.vdocument.in/reader034/viewer/2022051619/56649d575503460f94a35634/html5/thumbnails/29.jpg)
29
Direct derivation (100% ionization):
FAL model P
CICM radio mm limb
Gabriel corona, models A & D
My new measurements!
•Arbitrary shift in the X-direction. •Last points to the right : ~15-20 keV: contamination by thermal plasma•Points on the left: energies > 100 keV: fewer statistics•Assumed uniform target ionization! High-energy emissions should actually yield ~2.8 times higher densities…It’s a fitting game…