cps web site - 破壊を考慮した 惑星形成mosir/pub/2009/2009-10-19/01...2009/10/19 · f om...
TRANSCRIPT
破壊を考慮した惑星形成
イエナ大学小林 浩
09/10/13 14:09Google マップ - 地図検索
ページ 1/1http://maps.google.co.jp/
©2009 Google - 画像 ©2009 TerraMetrics, 地図データ ©2009 PPWK, Tele Atlas -
google map
内容
•惑星形成における破壊
• cratering v.s. catastrophic•破壊を考慮した惑星形成
•まとめ
•議論:ガス惑星形成
惑星形成標準モデル
105-6yr
105-7yr
106-8yr
109yr
破壊の重要性• 惑星形成後期になると、微惑星の衝突速度が原始惑星の重力により大きくなり、微惑星は衝突•破壊を起こす。
• 破片は次々に衝突を繰り返し、小さくなる。• 小さい天体はガス抵抗により消失する。 • その結果、衝突カスケードは天体の面密度を減少させる。
破壊の種類crateringcatastrophic
QD*
破片のモデル
mL m1 e
mass (m)
diffe
rent
ial n
umbe
r
m
ejecta
remnant0.01 1 100
1cratering catastrophic
mem1
mLm1
4
0.5
Eimp / m1QD*
面密度減少
• cratering と catastrophicのどちらが重要?
•一回の衝突破片の分布の影響は?
解析的に面密度減少時間を導出し、モデル依存性を明らかにする!
基礎方程式
F (m) = −ΩK
∞
mdm1
∞
0dm2m1f(m,m1,m2)Pcol ns(m1)ns(m2)
+ΩK
m
0dm1
∞
0dm2m1[1− f(m,m1,m2)]Pcol ns(m1)ns(m2),
∂mns(m)∂t
+∂F (m)
∂m= 0
質量保存則
質量フラックス (Tanaka et al. 1996)
m1 f (m,m1,m2): m1とm2の衝突の結果、m1から生成される質量mより小さい天体の総質量
減少時間
Σ =Σ0
1 + t/τdep
τdep =m
13max
(2− α)2Σ0h0ΩK
v(mmax)2
2QD(mmax)
−α+1
×− ln +
12− b
s1(α) + s2(α)
−1
.
衝突時間 (e.g., Wyatt et al. 2007)
α = (11+3p)/(6+3p) (e.g., O’Brien and Greenberg 2003)
仮定: v2/QD* >> 1 かつ v2/QD* ∝ m-p
1.65 1.7 1.75 1.8 1.852
610
20
60100
200
s1 v.s. s2
s2
s1
α破片の分布などは、τdep を変えない。
α =11 + 3p
6 + 3p
v2
QD
∝ m−p
(e.g. O’Brien and Greenberg 2003; Kobayashi and Tanaka 2009)
mass depletion
Σ =Σ0
1 + t/τdep
τdep =m
13max
(2− α)2Σ0h0ΩK
v(mmax)2
2QD(mmax)
−α+1
×− ln +
12− b
s1(α) + s2(α)
−1
.
cratering v.s. catastrophic
1.65 1.7 1.75 1.8 1.8512
61020
60100200
α
cratering
catastrophic
s 24-50 倍crateringの方が効く
α =11 + 3p
6 + 3p
v2
QD
∝ m−p
(e.g. O’Brien and Greenberg 2003; Kobayashi and Tanaka 2009)
統計的シミュレーション
質量
個数
シミュレーションとの比較破壊しきい値
v = 3 km/s
初期総質量: MEarth
(Benz and Asphaug 1999)
104 106 10810 4
10 3
10 2
10 1
100
104 106 10810 4
10 3
10 2
10 1
100
with crateringno cratering
Time [yr]
面密度/初
期面密度
QD = 3.5× 107
r1
1 cm
−0.38
+0.3
ρ
1 g cm−3
r1
1 cm
1.36[erg/g]
惑星形成での破壊
統計的シミュレーション
•質量座標を細かいビンに切り、ビンの中の天体の数の衝突による進化を計算する。•ビンの天体の速度(離心率、傾斜角)も天体の質量分布の進化に応じて変化すすので、同時に計算する。•動径方向にもメッシュを切り、天体のガス抵抗の移動による変化も計算する。
惑星形成中の質量分布
4 2 0 2 40
10
20
30
log10( r / 1 km )
log 1
0 N>
r
0 year
5.2AU
4 2 0 2 40
10
20
30
log10( r / 1 km )
log 1
0 N>
r
3.2 x 104 year
5.2AU
惑星形成中の質量分布
4 2 0 2 40
10
20
30
log10( r / 1 km )
log 1
0 N>
r
2.6 x 105 year
5.2AU
惑星形成中の質量分布
4 2 0 2 40
10
20
30
log10( r / 1 km )
log 1
0 N>
r
5.0 x 105 year
5.2AU
惑星形成中の質量分布
4 2 0 2 40
10
20
30
log10( r / 1 km )
log 1
0 N>
r
1.2 x 106 year
5.2AU
惑星形成中の質量分布
破壊を考慮した惑星形成
原始惑星
微惑星破片
消失τdep
τgrowKokubo and Ida (2002)
Kobayashi and Tanaka (2009)
最終原始惑星質量
1024 1026 1028105
106
107
M [g]
times
cale
[yr]
grow
dep
τgrow >> τdep
no more growth
原始惑星質量[g]
破壊込みの惑星成長原始惑星(M)の成長
面密度(Σ)の減少
dM
dt=
M
τgrow(M,Σ)
dΣdt
= − Στdep(M,Σ)
Kokubo and Ida (2002)
Kobayashi and Tanaka (2009)
惑星の最終質量
1 2 3 4 5 6 7 8910 20 30
0.001
0.01
0.1
1
10
isolation mass
distance [AU]
final
em
bryo
mas
s [M
]
ガス中ガス無し
林モデル
計算条件
ガスの寿命1千万年
1億年後
シミュレーション
結論•衝突カスケードによる面密度減少時間を導出。
• Crateringはcatastrophicより4-50倍重要である。
•天体消失時間を用いて原始惑星質量を見積もった。
•最終質量は火星質量程度。
ガス惑星形成
ガス惑星とデブリ円盤
•惑星形成時の破片がデブリ円盤を作る(e.g., Kenyon and Bromley 2008)
•ガス惑星とデブリ円盤が観測されている(Kalas et al. 2008)
Optical Images of an ExosolarPlanet 25 Light-Years from EarthPaul Kalas,1* James R. Graham,1 Eugene Chiang,1,2 Michael P. Fitzgerald,3 Mark Clampin,4Edwin S. Kite,2 Karl Stapelfeldt,5 Christian Marois,6 John Krist5
Fomalhaut, a bright star 7.7 parsecs (25 light-years) from Earth, harbors a belt of cold dust with astructure consistent with gravitational sculpting by an orbiting planet. Here, we present opticalobservations of an exoplanet candidate, Fomalhaut b. Fomalhaut b lies about 119 astronomical units (AU)from the star and 18 AU of the dust belt, matching predictions of its location. Hubble Space Telescopeobservations separated by 1.73 years reveal counterclockwise orbital motion. Dynamical models of theinteraction between the planet and the belt indicate that the planet’s mass is at most three times that ofJupiter; a higher mass would lead to gravitational disruption of the belt, matching predictions of its location.The flux detected at 0.8 mm is also consistent with that of a planet with mass no greater than a few timesthat of Jupiter. The brightness at 0.6 mm and the lack of detection at longer wavelengths suggest thatthe detected flux may include starlight reflected off a circumplanetary disk, with dimension comparableto the orbits of the Galilean satellites. We also observe variability of unknown origin at 0.6 mm.
About 15% of nearby stars are surroundedby smaller bodies that produce copiousamounts of fine dust via collisional ero-
sion (1). These “dusty debris disks” are analogsto our Kuiper Belt and can be imaged directlythrough the starlight they reflect or thermalemission from their dust grains. Debris disksmay be gravitationally sculpted by more massiveobjects; their structure gives indirect evidencefor the existence of accompanying planets [e.g.,(2, 3)]. Fomalhaut, an A3V star 7.69 pc fromthe Sun (4), is an excellent example: A planetcan explain the observed 15 AU offset betweenthe star and the geometric center of the belt,as well as the sharp truncation of the belt’s inneredge (3, 5–7). With an estimated age of 100 to300 million years (My) (8), any planet aroundFomalhaut would still be radiating its forma-tion heat and would be amenable to direct de-tection. The main observational challenge is thatFomalhaut is one of the brightest stars in thesky (apparent visual magnitude mV = 1.2 mag);to detect a planet around it requires the use ofspecialized techniques such as coronagraphy toartificially eclipse the star and suppress scatteredand diffracted light.
Detection of Fomalhaut b. Coronagraphicobservations with the Hubble Space Telescope(HST) in 2004 produced the first optical im-age of Fomalhaut’s dust belt and detected sev-
eral faint sources near Fomalhaut (6). Fomalhaut’sproper motion across the sky is 0.425 arc secper year in the southeast direction, which meansthat objects that are in the background will ap-pear to move northwest relative to the star. To
find common proper motion candidate sources,we observed Fomalhaut with the Keck II 10-mtelescope in 2005 and with HST in 2006 (9). InMay 2008, a comprehensive data analysis re-vealed that Fomalhaut b is physically associatedwith the star and displays orbital motion. Follow-up observations were then conducted at GeminiObservatory at 3.8 mm (9).
Fomalhaut b was confirmed as a real astro-physical object in six independent HST obser-vations at two optical wavelengths (0.6 mm and0.8 mm; Fig. 1 and table S1). It is comovingwith Fomalhaut, except for a 0.184 T 0.022 arcsec (1.41 T 0.17 AU) offset between 2004 and2006 (DT = 1.73 years) corresponding to 0.82 T0.10 AU year−1 projected motion relative toFomalhaut (9). If Fomalhaut b has an orbitthat is coplanar and nested within the dustbelt, then its semimajor axis is a ≈ 115 AU,close to that predicted by Quillen (7). An objectwith a = 115 AU in near-circular Keplerianmotion around a star with mass 2.0 times thatof the Sun has an orbital period of 872 yearsand a circular speed of 3.9 km s–1. The sixKeplerian orbital elements are unconstrainedby measurements at only two epochs; however,by comparing the deprojected space velocity(5.5–0.7
+1.1 km s–1) with the circular speed, we
RESEARCHARTICLES
1Astronomy Department, University of California, Berkeley, CA94720, USA. 2Department of Earth and Planetary Science,University of California, Berkeley, CA 94720, USA. 3Institute ofGeophysics and Planetary Science, Lawrence Livermore NationalLaboratory, Livermore, CA 94551, USA. 4Exoplanets and StellarAstrophysics Laboratory, Goddard Space Flight Center, Greenbelt,MD 20771, USA. 5MS 183-900, Jet Propulsion Laboratory,California Institute of Technology, Pasadena, CA 91109, USA.6Herzberg Institute for Astrophysics, Victoria, British ColumbiaV9E 2E7, Canada.
*To whom correspondence should be addressed. E-mail:[email protected]
Fig. 1. HST coronagraphic image of Fomalhaut at 0.6 mm, showing the location of Fomalhaut b (whitesquare) 12.7 arc sec radius from the star and just within the inner boundary of the dust belt. All the otherapparent objects in the field are either background stars and galaxies or false positives. The fainter lower halfof the dust belt lies behind the sky plane. To obtain an orientation with north up and east left, this figureshould be rotated 66.0° counterclockwise. The yellow circlemarks the location of the star behind the occultingspot. The yellow ellipse has a semimajor axis of 30 AU at Fomalhaut (3.9 arc sec) that corresponds to the orbit ofNeptune in our solar system. The inset is a composite image showing the location of Fomalhaut b in 2004 and2006 relative to Fomalhaut. Bounding Fomalhaut b are two elliptical annuli that are identical to those shownfor Fomalhaut’s dust belt (6), except that here the inner and outer annuli have semimajor axes of 114.2 and115.9 AU, respectively. The motion of Fomalhaut b therefore appears to be nested within the dust belt.
www.sciencemag.org SCIENCE VOL 322 28 NOVEMBER 2008 1345
CORRECTED 2 JANUARY 2009; SEE LAST PAGE
on O
cto
ber
14, 2009
ww
w.s
cie
ncem
ag.o
rgD
ow
nlo
aded fro
m
破壊(デブリ円盤)を考慮したガス惑星形成が重要
Kalas et al. Science 322,1345 (2008)
ガス惑星形成の臨界コア質量
で静水圧平衡の大気を持てなくなりMcore ≥ 10M⊕ガス惑星形成 (e.g., Ikoma et al. 2001)
0.1 1 10
0.001
0.01
0.1
1
10
100
原始惑星質量
円盤質量 [0.008 Msun]惑星質量 [M
Earth] 見積もりと同程度の
天体ができる。
臨界惑星核質量は超えられない。
臨界質量
at 5 AU
ガス無し
ガスあり
5-30AU
大気の効果
log rp[cm]2 4 6 8 10
Rc: collisional radius Rcore: target radiusrp: projectile radius
Inaba and Ikoma (2003)
100 102 104 106100
102
104
106
108
小さい天体の重要性
天体半径[cm]
落下時間 [年
]
smin = 1cm
at 5 AU
Kenyon and Bromley (2009)
0.1 1 10
0.001
0.01
0.1
1
10
100
ガス惑星のための惑星核形成
円盤質量 [0.008 Msun]
惑星質量[M
Earth] 臨界質量
at 5 AU
今後の研究
•大気の効果でガス惑星はできる(Inaba et al. 2003)。
•どのような条件で大気の効果は有効か。
•大気の効果を含めた、原始惑星の質量はどのように表現されるか。
ご清聴ありがとうございました。
Thank you for your attention.
Danke für Ihre Aufmerksamkeit.