400 likes | 499 Views
Modeling Disks of sgB[e] Stars. Jon E. Bjorkman Ritter Observatory. Dusty Hot Star Winds. Hot stars with dust : B[e] WR Novae and Supernovae Wind must cool below condensation temperature Dust forms at large distances
E N D
Modeling Disks of sgB[e] Stars Jon E. Bjorkman Ritter Observatory
Dusty Hot Star Winds • Hot stars with dust: • B[e] • WR • Novae and Supernovae • Wind must cool below condensation temperature • Dust forms at large distances • Problem: density must be large enough that reaction rates are faster than flow times Zickgraf, et al. 1986
Eta Carinae Morse & Davidson 1996
General Wind Kinematics • Radial Momentum Equation • Radial Motion
General Wind Kinematics • Azimuthal Motion
Bi-Stability & Disks • Ionization shift at low latitudes • Higher mass loss • Lower terminal speed Lamers & Pauldrach 1991
Rotationally Induced Bi-Stability Terminal Speed Mass Loss Need additional factor Pelupessy et al. 2000
Rotating Stellar Winds Low Density,High V∞ High Density, Low V∞ Bjorkman & Cassinelli, 1993
Ionization Structure Krauss & Lamers 2003
WCD Inhibition Radial Force Only Non-radial Force Effects Owocki, Cranmer, & Gayley 1996
WCDs and Be Stars • Non-radial line forces (prevent disk formation) • Outflow speed too large (~400 km/s) • Density too small (to explain IR excess) • Disk “leaks” • Material falls back onto star • Material flows outward through disk • Must put material into orbit • Must remove radiative acceleration
Magnetic Channeling Cassinelli et al. 2002 Owocki & ud-Doula 2003
Asymmetric Mass Ejections • Kroll’s gravity filter: • Point “explosion” • Material thrown backward falls onto star • Material thrown forward goes into orbit Stellar Bright Spot Model Spot + Line-Force Cutoff Owocki 2003
Keplerian (Orbiting) Disks • Fluid Equations • Vertical scale height (Keplerian orbit) (Hydrostatic) (Scale height)
Viscosity in Keplerian Disks • Viscosity • Diffusion Timescale (eddy viscosity) Lynden-Bell & Pringle 1974
Viscous Decretion Disk • Lee, Saio, Osaki 1991
Disk Temperature Flat Reprocessing Disk Flared Reprocessing Disk
Disk Winds Model for HD 87643 Oudmaijer et al. 1998
Power Law Approximations • Keplerian Decretion Disk • Flaring
Monte Carlo Radiation Transfer • Divide stellar luminosity into equal energy packets • Pick random starting location and direction • Transport packet to random interaction location • Randomly scatter or absorb photon packet • When photon escapes, place in observation bin (frequency and direction) REPEAT 106-109 times
MC Radiative Equilibrium • Sum energy absorbed by each cell • Radiative equilibrium gives temperature • When photon is absorbed, reemit at new frequency, depending on T
T Tauri Disk Temperature Whitney, Indebetouw, Bjorkman, & Wood 2004
T Tauri Disk Temperature Water Ice Snow Line Methane Ice
Effect of Disk on Temperature • Inner edge of disk • heats up to optically thin radiative equilibrium temperature • At large radii • outer disk is shielded by inner disk • temperatures lowered at disk mid-plane • Does not solve dust formation problem; requires • high density at condensation radius • additional opacity interior to condensation radius
Dust Formation • Reaction network • Timescales • Condensation Condition Porter 2003 (Gail & Sedlmayr 1988)
SED Viscous Decretion Viscous Decretion Bi-stability Bi-stability Porter 2003
NLTE Monte Carlo RT • Gas opacity depends on: • temperature • degree of ionization • level populations • During Monte Carlo simulation: • sample radiative rates • Radiative Equilibrium • Whenever photon is absorbed, re-emit it • After Monte Carlo simulation: • solve rate equations • update level populations and gas temperature • update disk density (solve hydrostatic equilibrium) determined by radiation field
sgB[e] Density (pure H model) Bi-stability Viscous Decretion
Gas (Electron) Temperature Bi-stability Viscous Decretion
Dust Temperature Bi-stability Viscous Decretion
Mid-Plane Temp Bi-stability Viscous Decretion Rdust = 1300 R* Rdust = 400 R*
Density Bi-stability Viscous Decretion
sgB[e] Model SED Bi-stability Viscous Decretion l (mm) l (mm)
IR Spectroscopy Roche, Aitken, & Smith 1993
Dust Properties Large Dust Grains Wood, Wolff, Bjorkman, & Whitney 2001
YSO (GM Aur) SED • Inner Disk Hole = 4 AU Rice et al. 2003
Line-Blanketed Disk Opacity Bjorkman, Bjorkman, & Wood 2000
Conclusions • Bi-Stability: • Pros: • Provides better shielding for dust formation • Cons: • Requires small condensation radius • Viscous Decretion • Pros: • Slow outflow enables much larger condensation radius • Disk wind may produce low velocity outflow • Cons: • Dust optical depth is much too small • Generally, • need to increase disk outflow rate (without increasing free-free excess) • Or provide more shielding to decrease condensation radius