Stephen St.Vincent - Swarthmore College
Research Home | Journal | Articles | To-Do | IDL code
Back to main images page
Here's a direct comparison of thermal broadening (top) versus no broadening (bottom). The line center moved to negative with broadening, so I don't think that this is right yet.
27 July 2006
Here are two X-ray light curve, with 40 random divisions.
The top image is the light curve for timestep 6 of the θ1 Ori C simulation, which is before infall is supposed to have started occurring. The bottom image is the corresponding emission measure contour plot. I can't readily explain the small dips near φ=0.05 and φ=0.95.
The top image is the light curve for timestep 13 of the θ1 Ori C simulation, which is before infall is supposed to have started occurring. The bottom image is the corresponding emission measure contour plot. It appears that most of the Hα emission is coming from that red bar on the magnetic axis, which is causing the horizontal wings on the light curve.
26 July 2006
Above are the light curves for two randomly-generated simulations. A little different, no?
Here's the color contour of zone 52. The black contours mark areas that do not emit Hα. 46% of total emission is from the initially occulted portion of the star.
Here's a shot of the light curve when all of the slices are of simulation zone 65. There are 100 slices, 100 timesteps. Looks pretty much like the zone 75 image.
Here's the color contour of zone 52. The black contours mark areas that do not emit Hα. It looks like most of the emission is from areas on the under side of the star. In fact, 71% of the total emission come from the initially occulted area.
This is the same as above, but for zone 52. It looks completely backwards. I'm going to investigate the cause of this today.
Here's the corresponding logarithmic plot for yesterday's linear plot. Again, they look almost identical, so I'm happy that this is working.
25 July 2005
Here's a light curve after I fixed everything. That hump in the middle is back... I really want to believe that it's because the highest-EM material is becoming uncovered starting at φ≈0.3. My rotation still seems to be working, so I don't really know what else it could be.
20 July 2006
The newest version of the GUI.
19 July 2006
Here's the same H-alpha test as yesterday, but with a linear y-axis. It looks almost the same, but I checked the values and everything is actually correct.
Here's an image showing the results of thermal broadening as currently implemented. It sure smoothed out, but the question in my mind is, did it smooth it too much? In the movies of this (of which there are not currently any up on my website), you can see the line center move slightly back and forth across zero.
18 July 2006
This is a test of the Hα variability. Since the plot from yesterday looked so choppy, David decided that I should do a test with the same slice in all regions (azimuthally symmetric) to see how it looked. So I made a simulation with 100 divisions, all of which were zone 75 from the 2D simulation. It looks a little odd, with the maxima occuring at φ=0.2 and 0.8. Also, there's a local minimum at φ=0, which is a little unexpected. My guess is that this is because at φ=0, all of the occulted material is emitting Hα. But most of the density is where the gas is emitting mostly X-rays but still some Hα, so these small, high-density regions may be the dominant source of Hα and thus produce the large global minimum at φ=0.5.
17 July 2006
This is a plot for Hα. Not sure how I feel about this, or if the temperature threshold is even correct.
Here's the emission measure vs time histogram, now that the rotation has been fixed. It looks better, but there's still a little hump in the middle near phase=0.5. I'm not really sure what that's from. Also, the values given for inclination and obliquity are incorrect. They've been fixed in the code now.
13 July 2006
Here's a quick test of the thermal broadening. I've set vlos=2000 km/s, T=108 K, EM=1054 cm-3, and there are bins every 0.3 km/s. It looks good, but the broadening is only really appreciable since the velocity is so large and the grid is so fine.
Here's the same image, but with two points. The curves mesh exactly as expected.
The top image does not use thermal broadening, while the bottom one does. See the difference? I can't.
11 July 2006
Here's what I've decided to do with the logarithmic line profiles. The axis ranges will be fixed to these same values. The red curve represents the temperature cutoff, and will only be displayed if maxthresh is greater than zero.
07 July 2006
Here's my first emission measure vs time histogram for the full 3D simulation. There are 40 divisions, 100 timesteps, obliquity is 42 degrees and inclination is 45 degrees.