星形成ゼミ：2019年8⽉2⽇

古屋 玲 (徳島⼤学・教養教育院)

abundance ratios, and find that the abundance ratio is as high as 40 near the emission front, and decreases to the−3
typical valueALMA
in the solar
system of 5.5 in ofa smallangularStructures
scale of 4” = 2.6×10= 560
au. We alsofind that the in the ρ
Observations
Layereddue
to pcCO
Selective
Dissociation
12
13
I( CO(2–1))/I(
CO(2–1))intensity ratio is veryPDR
high (>21) in the flat-spectrum young stellar object, GY-51, locatedOphiuchi
A Plane-parallelin the PDR. The enhancement
indicates that
the CO isotopologues
1 of the ratios
2
3 the UV radiation
3 significantly aﬀects
3
4
5
M.
Yamagishi
,
C.
Hara
,
R.
Kawabe,F.
Nakamura,
T.
Kamazaki,T.
Takekoshi,
Y.
Shimajiri, H.
via selective dissociationintheoverallρOphiuchi
A
PDR,
and
that
the
ρ
Ophiuchi
A
PDR
has
a
plane-parallel
67
8
Nomura,S.
TakakuwaandJ.
DiFrancesco
structure.
1
Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, Chuo-ku, Sagamihara, Kanagawa
Accepted by ApJ
252-5210, Japan; 2 Graduate School of Science, The University of Tokyo, Bunkyo-ku, 113-0033, Tokyo, Japan; 3 Nahttps://arxiv.org/pdf/1903.02192
tional Astronomical Observatory of Japan, National Institutes of Natural Sciences, 2-21-1 Osawa, Mitaka, Tokyo 1814
58588,
Japan;
InstituteofAstronomy,
The
UniversityofTokyo,
2-21-1
Osawa,
Mitaka,Tokyo
181-0015,
Japan;
Lab主な結果
rho
Ophiuchi
Aの光解離領域(解離⾯をちょうど横から⾒ている)で観測されたCO同位体の⽐は,Syntheticoratoire
molecularline observationsof thefirsthydrostatic
core
from chemical
calcuAIM,
CEA/DSM-CNRS-Universit
Paris
Diderot,
IRFU/Service
dAstrophysique,
CEA
Saclay, F-91191 Gif12 CO(2–1) (red),
13
18
Fig.
3.—
Averaged
spectraofCO(2–1)(green),andCO(2–1)(blue)
emissions
6
lations sur-Yvette,selective
dissociationでうまく説明できる.
France;Department of Earth and Planetary Sciences, Tokyo
Institute of Technology, 2-12-1 Ookayama,
−1 ).
7
at
the
layeredstructure(thewhite-box
region
in
Fig.
2,
V
=2.2
km
s
1
1
lsr and Astronomy,
Meguro,1 ,Tokyo,152-8551,
Japan;
Department
Graduate
School of Science and EngineerAlison K. Young
Matthew
R. Bate, Tim
J. Harriesof1 Physics
and David
M. Acreman
8
ing,ofKagoshima
1-21-35
Korimoto,
Kagoshima,
Kagoshima890-0065,
Japan;
NRC
Herzberg Astronomy
11213
18 O(2–1) (blue) emissions
Department
Physics andUniversity,
Astronomy,
UniversityofExeter,
StockerRoad,
Exeter
EX4
4QL,
UK
Fig.
3.—
Averaged
spectraofCO(2–1)(red),
CO(2–1)(green),andC
観測and Astrophysics Research Centre, 5071 W Saanich Road, Victoria BC V9E 2E7, British Columbia,
Canada−1
at the layered structure (the white-box region in Fig. 2, Vlsr =2.2 km s ).
E-mail contact:ayoung at astro.ex.ac.uk
ALMA
12mアレイと森⽥アレイ,
J=2-1 of 12CO, 13CO and C18O, ⾓度分解能~1”E-mail contact: yamagish at ir.isas.jaxa.jp
The first stable object todevelop in thelow–mass star 18formation process has long been predicted to be the first
12
13
CO(J=2–1),
C Oyet(J=2–1),
and 1.3 mm
continuum
maps ofMore
the ρspecificOphiuchi A photohydrostatic We
coreanalyze(FHSC). CO(J=2–1),
Despite much eﬀort,
it has still
to be definitively
observed
in nature.
dissociation
region(PDR)obtainedwith ALMA.
Layeredthe
structuresof thethreefaint,
CO butisotopologues
with an angularobservationalsignatures are
requiredto−3enable observers
to distinguish
FHSC fromyoung,
more evolved
separationof 10”synthetic= 6.6×10spectral
pc =line
1400
au are clearly
detectedaroundtheBe+star,S1(i.e.,calculated
each frontfromof emission shifts
protostars. Here
we present
observations
for
CO,
SO,
CS
and
HCO
that
were
12
13
from the near to far side inorder ofCO, 13CO, and andC18 O).
We Carlo
estimate
the spatial
variations
X(CO)/X(C18 O)
radiation (magneto)hydrodynamical
models,chemicalmodellingMonteradiative
transfer.
HCO+of(1− 0)abundance
ratios, and find that the abundance ratio is as high as 40 near the emission front, and decreases to theand SO (87 −76 ) spectra of the FHSC show variations for observations at a low inclination which may allow a candidate−3
value infrom
the solarsystemof 5.5in a small
angularscaleFHSC
of 4” outflow= 2.6×10
pc = to560beau.
We also find that the
FHSC to betypicaldistinguisheda moreevolvedobject.
We findthat theis unlikely
detectable
12
13
CO(2–1))/I(
CO(2–1))intensitysources
ratio iswith
veryslowhighoutflows(>21) inthat
the flat-spectrum
stellarobject, GY-51, locatedwith ALMA,I(which
would discountthe observed
are currently young
identified
as candidate
the PDR.
enhancementof theratiosobservations
indicates thattheobserved
UV radiation
significantly
aﬀectsthe CO isotopologues
FHSCs. Weincompare
theTheresultsof simulated
ALMA
with
candidate
FHSCs andrecommend
viaandselective
the overall
ρ Ophiuchi
A PDR,
Oph A SM1N
N6-mmdissociation
as the mostinpromising
candidates
to follow
up. and that the ρ Ophiuchi A PDR has a plane-parallel
structure.
Accepted by MNRAS
Accepted by ApJ
https://arxiv.org/pdf/1906.01324
https://arxiv.org/pdf/1903.02192

72

Sequential Star Formation in the filamentary structures of the Planck Galactic cold
Synthetic molecular line observations of the first hydrostatic core from chemical calcuclump G181.84+0.31
lations
Lixia Yuan1,2,3 , Ming Zhu1,3 , Tie Liu4,5 , Jinghua Yuan1 , Yuefang Wu6 , Kee-Tae Kim4 , Ke Wang7,8 ,18 O)10,11
1
1,2,3 K. Young1 , Matthew
9 R.
Fig. 4.—X(13CO)/X(Cat V
km1 s−1 . Contours are the same
AlisonBate, Tim
J. abundanceHarries1 ratio
and maps
David
M.
Acremanlsr =2.2–2.6
Chenlin Zhou
, Ken’ichi Tatematsu
and
Nario
Kuno
in Fig.
1. ratio
The dashedtheStocker
direction
to calculatethetheabundance
1 13 CO)/X(C18 O) abundanceFig. 4.— X(
mapsline
at indicatesVoflsrExeter,=2.2–2.6
km Road,s−1used. Contours
same ratio profiles
Department of Physicsasand
Astronomy,
University
Exeter
EX4are4QL,
UKin Fig. 5.
as in Fig. 1.
Thecontact:dashedayoung
line indicatesthe direction used to calculate the abundance ratio profiles
E-mailat astro.ex.ac.uk

– 13 –

Fig. 5.— (top) Spatial variations of the X(13 CO)/X(C18 O) abundance ratio at Vlsr
(color-corded lines; left axis) and those of the Spitzer 4.5 µm (squares; right acontinuum intensities (diamonds; right axis) as a function of the distance from
calculated from averaged profiles in a rectangular region with a width of 1.′′ 5 aline in Fig. 4. The Spitzer 4.5 µm and 1.3 mm continuum data are normaliz
mum values of 23.9 MJy sr−1 and 1.0 mJy beam−1 , respectively. The dotted lX(13 CO)/X(C18 O) abundance ratio of 5.5 in the solar system (Anders & GreveFig. 2.— 1.3 mm continuum map and color-composite channel maps of 12 CO(2–1) (red), 13 CO(2–1) (green), and C18 O(2–1) (blue) emissions at Vlsr =1.0–3.0 km s−1 for the yellow dashed box in
Figure 1. Contours are the same as in Fig. 1. The color range in the 1.3 mm continuum map is
12 CO(2–1)
左上以外のパネルは12CO,
0–3 mJybeam−1 , while thechannel
peak intensity
is 37.7
beam−1 associatedwith
GY-51. The colorFig.
2.— 1.3 mm continuum map andcolor-composite
maps
of mJy(red), 13CO(2–range is common among the channel maps−1
and is linear from 0 % to 80% of the peak intensity in1)(green),and C18 O(2–1) (blue) emissionsatV=1.0–3.0kms
for the yellow dashed box in
13CO,
C18Oの三⾊図
lsr
this area (12 CO: 4.98 Jy beam−1 ; 13 CO: 2.43 Jy beam−1 ; C18 O: 1.32 Jy beam−1 ).

Figure 1. Contours are the same as in Fig. 1. The color range in the 1.3 mm continuum map iscascade.
compression
the andspectrumof density
fluctuations
2005; Kowaldeforms
et al. 2007)
accounts
for the formationof to as the shock propagation distance reduced by the compression
thefilament,
width ofcan be as
estimatedratio. In
Foraddition,
an obliquethedense
widthfilaments
can be estimatedbedifferent
from
that
of
turbulentvelocities(Beresnyak
et
al.
dense
filaments
in
MCs(Robertson&
Goldreich
2018;
Mocz
• シュミレーション結果の解析がなく,
オーダー評価のみ.
もっとも肝⼼な結論部は,
議論 2005;
the shockpropagation distance reduced by the compression& Burkhart
Kowal2018)
et al.4 . 2007) and accounts for the formation of (Eq.as(16)),
L the width can be estimated asratio. For an oblique filament,
Thefilaments
density contrast
the shockdense
in MCsacross(Robertson& isGoldreich 2018; Mocz
p
⇠
,
(18)
o
(Eq.
(16)),
4
2MA1& Burkhart 2018) . ⇢2 pL
⇡ the
2Mshock(16)
A1
The density contrast ⇢across
is
p
which
dependsonthemagnetization
andthe ,filament orienta⇠(18)o1
2MA1 by assuming aption. For the densest perpendicular filaments,⇢2
when the magnetic pressure⇡dominates
2MA1in the downstream re-(16) constant compression ratio determined by Ms1 (⇠ Ms ), we
which depends on the magnetization and the filament orienta⇢1gion, andcanapproximately have the filament width as (Eq. (17))⇢2tion. For the densest perpendicular filaments, by assuming a2⇡
Ms1(17) rewhen the magnetic pressure
dominatesin the downstream
constant compression ratioLdetermined by Ms1 (⇠ Ms ), we
⇢
1
• 結局,
,
(19)
gion,
and この研究だけからは明瞭な結論を引き出すのは困難と述べている.
? ⇠the filament2
canapproximately
have
width
as
(Eq.
(17))
M
s
⇢2dominates
when
thethermalpressure
inthedownstream
re2✓…fully understandthe
turbulence
dynamics
in multi-phase
ISM,alonga modelfor(upper
two-phaseturbulence⇡Mof(17)viewedofdensityfluctuations
at different
times,
z-axispanel)MHD
and x-axis(lower should
panel). be applied…
s1
gion, where (a)⇢1 Sliceand ⇢plots
aredensities
theupstream
and
downL
2
which depends on the compressibility of medium. Provided
⇢1
(19)
? ⇠
streammedia,
M
is
the
shock
AlfvénMach
numbercorre✓…densityA1
gradient is not always a good indicator L(で,⇠ 速度勾配は？)…
10 pc and Ms ⇠ 10 (Schneider
et2 ,al. 2013; Henshaw
M
when
the thermalpressurecomponent
dominatesofinupstream
the downstream
sponding
to the transverse
magnetic re- et al. 2016), we find that ? is of the ordersof 0.1 pc, consistent
field,whereand M
is ⇢theshocksonic ofMach
number (seegion,
⇢1s1anddensities
the upstream
and Apdown- withwhichobservations
by,one.g.,
et al.of(2011);
Koch
2 are
depends
theArzoumaniancompressibility
medium.
Providedpendixmedia,
A for thedetailedThe
above
compresstream
MA1
is the derivation).
shock AlfvénMachnumbercorre- & Rosolowsky
L ⇠ 10 pc(2015).
and Ms ⇠ 10 (Schneider et al. 2013; Henshaw
sion
ratiosareconsistentwithearlierstudiesonisothermal
sponding to the transverse component of upstream magnetic
et al. 2016), we find that ? is of the order of 0.1 pc, consistentshocks(see,e.g.,Draine2011).
In
highlysupersonicMHD
DISCUSSION
field, and Ms1 is the shock sonic Mach number (see Apwith observations5. by,e.g., Arzoumanian et al. (2011); Kochturbulence,theshock-compressed
densityvariesindifferent
pendix A for the detailed derivation). The above compres& Rosolowsky(2015).
Motivatedby ample
observations of interstellar density filturbulence
regimes.
sion ratios are consistent with earlier studies on isothermal aments, there have been numerous theoretical arguments and(1) M
MsDraine
1. 2011).
In a weakly
magnetized
medium,A e.g.,
shocks(see,
In highly
supersonic
MHD numerical experiments on theirformation mechanisms (e.g.,
5. DISCUSSION
the shock compression in MHD turbulence is similar to that
turbulence,the shock-compressed
densitydensity
varies enhancein different Lazarian 1993; Hennebelle 2013; Smith et al. 2014; Ntorin hydrodynamic
turbulence. The resulting
Motivatedbyample
observationsofinterstellar
density filmousi
et
al.
2014;
Inoue
&
Inutsuka
2016;
Smith
et
al.
2016;
turbulence
regimes.
ment solely depends on Ms1 (⇠ Ms ) (Eq. (17), see Padoan &aments,thereBanda-Barragánhave been numerous
theoreticalargumentsandFederrath
2016;
et
al.
2018).
In
this
work,
(1)
M
M
1.
In
a
weaklymagnetized
medium,As
Nordlund 2011; Federrath & Banerjee 2015; Federrath 2016;
numerical
experimentsdynamics
on theirofformationon the fundamental
Alfvénic mechanisms
turbulence, (e.g.,theshock compression
MHDWeturbulenceis similar
to that based
Robertson& Goldreich in2018).
note that the
weak magHennebelle
2013;
Smith et al.
2014; Ntorwe Lazarian
propose a 1993;
formationmechanismof low-density
parallelinnetic
hydrodynamic
turbulence.
The
resulting density
enhancefield can be amplified
by thedownstream
turbulence(Xumousiinettheal.ISM.
2014; Inoue & Inutsuka 2016; Smith et al. 2016;
ment
solely depends
on Ms1 (⇠ Ms ) (Eq. (17), see Padoan & filaments& Lazarian
2016, 2017a).
Federrath
2016;inBanda-Barragán
et numericalal. 2018).studies,
In this work,
Cold
HI filaments
the WNM. Based on(2) Ms2011;
MFederrath
MBanerjee
1 >2015;MA . Federrath
In a strongly
Nordlund2016;
A > 1 or &s
Saury et al.on(2014)
argued that thedynamics
dynamics of the
WNM andthe fundamental
Alfvénic
turbulence,magnetized&medium,
Eq.2018).
(16) applies
to thethat
shock
Robertson
Goldreich
We notethe compresweak mag- the based
CNM
are tightly
coupled, andthe velocity
be- parallel
we
proposea formation
mechanismof dispersion
low-density
netic field can be amplified by the downstream turbulence (Xu tween
CNM clumps
close to the velocity dispersion within
4 In 3D supersonic turbulence, due to the interaction of density inhomofilamentsin the is
ISM.
&geneities
Lazarian
2016,
2017a).
the extended
WNM.
It implies
thatthe Based
denser on
CNM
embed- studies,with multiple shocks and turbulent mixing (Pan & Scannapieco
ColdHI
filaments
in
the
WNM.
numerical(2) M
Mor Mofs 2D sheet-like
1 > Ms
A > 1instead
A . In a2010),
filamentary
structures
structures
arestrongly
more
ded in the etdiffuseWNMargued
is dynamically
mixed by the
turbual. (2014)
that the dynamics
of the
WNM andlikely to be generated.
magnetized
medium, Eq. (16) applies to the shock compres- lentSaury
velocities
in
the
WNM,
and
thus
the
cold
HI
filaments
the CNM are tightly coupled, and the velocity dispersion between CNM clumps is close to the velocity dispersion within
4 In 3D supersonic turbulence, due to the interaction of density inhomothe extended WNM. It implies that the denser CNM embedgeneities with multiple shocks and turbulent mixing (Pan & Scannapieco2010), filamentary structures instead of 2D sheet-like structures are moreded in the diffuse WNM is dynamically mixed by the turbulikely to be generated.
velocities
inalong
the WNM,
andpanel)thusandthex-axis
cold (right
HI filaments(b) Slice plots of density fluctuations superposed with magnetic field lines at t = lent
6 teddy,L, viewed
z-axis (leftpanel).
F IG . 2.— The structures of density fluctuations and magnetic fields.

These results may be due to the fact that the global collapse of the extended filamentary structure S242 is driven by
Sequential
Star
Formationin theof the Planck
Galactic
an edge eﬀect,
for which thefilament
edges collapse
first filamentary
and then furtherstructures
trigger star formation
activities inward.
We coldclump
G181.84+0.31
identify threesubstructuresin G181.84 and estimate their critical masses per unit length, which are ∼ 101±15 M⊙
1
4−1 1,2,3 , Ming Zhu−1
Yuan
Tie Liu4,5 , Jinghua
Yuan
, all
Yuefang
Wu6the, Kee-Tae, Ke
Wang7,8 ,pc−1 , 56±8LixiaM⊙ pc
and 28±4 M⊙ pc 1,3,, respectively.
These values
arelower thanobservedKim
values(∼ 2009
10,11Zhou, Ken’ichi
Tatematsu
and Nario Kuno
M⊙ pc−1 ), Chenlin
suggestingthat1,2,3
these
sub-structures
are gravitationallyunstable.

74

Accepted by MNRAS
https://arxiv.org/pdf/1905.12258

46

One
Solution to the Mass Budget Problem for Planet Formation: Optically Thick Disks
観測と使⽤データwith Dust Scattering

SCUBA-2 850 micron, Herschell, WISE, Spitzer（以上は連続波）, 野辺⼭45m鏡でN2H+ and HCO+輝線

Zhaohuan Zhu1 , Shangjia Zhang1 , Yan-Fei Jiang2 , Akimasa Kataoka3 , Tilman Birnstiel4 , Cornelis P.
Dullemond5, Sean M. Andrews6 , Jane Huang6 , Laura M. Pérez7 , John M. Carpenter8 , Xue-Ning Bai9 ,David J. Wilner6 and Luca Ricci10
1

Department of Physics and Astronomy, University of Nevada, Las Vegas, 4505 S. Maryland Pkwy, Las Vegas,
観測事実
NV
89154, USA; 2 Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA, USA; 34
National
Astronomical
ObservatoryofJapan,
Osawa
2-21-1,
Mitaka,Tokyo
181-8588,
Japan;
University Observatory,
1. フィラメントに沿って分⼦雲コアが形成されている⼀⽅, フィラメントが交差する中⼼領域ではすでFaculty of Physics, Ludwig-Maximilians-Universität München, Scheinerstr. 1, 81679 Munich, Germany; 5 Zentrum fürに若い星が存在する.
Astronomie,
Heidelberg University, Albert Ueberle Str. 2, 69120 Heidelberg, Germany; 6 Harvard-Smithsonian Center
7
for Astrophysics,
60 Garden Street, Cambridge, MA 02138,
Departamento
de=Astronomı́a,
Universidad de
2. フィラメントと中⼼領域で⼤局的な速度勾配(~1USA;kms/s
per pc, error
0.05 km/s)が存在する.
Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile; 8 Joint ALMA Observatory, Avenida Alonso deCórdova 3107, Vitacura, Santiago, Chile; 9 Institute for Advanced Study and Tsinghua Center for Astrophysics,
Tsinghua University, Beijing 100084, China; 10 Department of Physics and Astronomy, California State University
Northridge,
18111 Nordhoﬀ Street, Northridge, CA 91130, USA
著者の主張
E-mail contact: zhaohuan.zhu at unlv.edu

1. surveys
連鎖的星形成ALMA
have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass
2. globalundercollapsein exoplanets,
the assumption that the mm disk continuum emission is optically thin. This optically thinassumption seems to be supported by recent DSHARP observations where the measured optical depths of spatiallyresolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission
from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated.
An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined

igure 12. The skeletons of the filament network extracted with
he FilFinder, color-coded with the centroid velocities of the
CO+ (1-0) spectral line. Magenta ellipses represent the SCUBAcompact sources. Green contours show the distribution of theCO+ integrated intensity, levels are from 10% to 90% by theterval of 15% of the peak value (13 K km s 1 ).

radients along Fb and Fc in the northern end of EFS S242

4.3 Gravitational fragmentation of sub-filam

As shown in Figure 12, G181 contains sub-filame
9 dense cores embedding in them. To investigate
mentation of sub-filaments and hence the form
dense cores, we compare the sub-structures of G18
model of a collapsing isothermal cylinder, support
ther thermal or non-thermal motions. André et a
has suggested that thermally supercritical filamen
the mass per unit length is greater than the criti(Mline,crit ; i.e., Mline > Mline,crit ), are associatedstellar clumps/cores and star formation activities.
there are few Herschel pre-stellar clumps/cores andded protostars associated with thermally subcrit
ments (Mline < Mline,crit ).
For thermal support, the critical mass per un(Mline,crit = 2c2s /G; where cs is the isothermal soun
G is the gravitationally constant) is needed for ato be gravitationally unstable to radial contrac
fragmentation along its length (Inutsuka & Miyam
The thermal critical line mass Mline,crit is ⇠ 16 M
(Tgas /10 K) for gas filaments (André et al. 2014).
mentary structures of G181 have an observed mass
length of ⇠ 200 M pc 1 (Dewangan et al. 2017).
served value largely exceeds the critical value of ⇠ 1
pc 1 (at T = 10-30 K).
Taking into account additional supports from t
pressure, the total 1D velocity dispersion is calcula
the following equation (Fuller & Myers 1992)
s

v2

✓

◆

Synthetic molecular line observations of the first hydrostatic coreaがFC形成時, bとdは断熱収縮期, gは２回⽬のcollpase後で原始星期なので,gがbやdと区別がつくか？

(=進化段階を判定する指標となりうるか？)右と左のパネルで差が⾒えるか？ (＝どっちの⽅向から⾒ると確実に判定できるか？)
上下のパネルで優劣が顕著か？ (=どっちの分⼦を選ぶ？)

レポーター感想：n(H2)を知るために必要な塵熱輻射諸元の不定性も考えると, アバンダンス⼀桁の差を検出してもそれが有意であることを⽰すのは難しそう. ⼀⽅, 分⼦雲コア中⼼部に初めて光学的に厚い核ができたときの変化は劇的.

log [ n (HCO+) / n (H2) ]分⼦流の軸⽅向から覗き込んだとき (横軸の単位はAU)
CO

CN

-10.0
-12.0
-14.0
-16.0
-18.0
-20.0
1

1

+

-8.0

10

10

100

100

1000
NH3

a
b
d
g

-7.0
-9.0
-11.0
-13.0
-15.0
-17.0a-19.0b
-21.0d-23.0g
-25.0
1000
1

CS

log [ n (CS) / n (H2) ]

a
b
d
g

log [ n (CN) / n (H2) ]

-7.0
CS
-4.0
-9.0
-11.0MNRAS 000, 1–14 (2019) -5.0
-6.0
-13.0
-15.0
-7.0
-17.0
-8.0
-19.0
-9.0
-21.0
-10.0
-23.0
-11.0
-25.0
-12.0
1
10
100
1000
-13.0

log [ n (CO) / n (H2) ]log [ n (CS) / n (H2) ]Figure 5. Average abundances of selected species calculated for the RHD mode
hetic molecular line observationsof thefirstcore direction
15 (parallel to the rotation axis).
Figure 6. Sameas forFig. hydrostatic
5 but in the vertical

a
b
d
g

10

100

1000

log [ n (HCO+) / n (H2) ]log [ n (H2CN) / n (H2) ]log [ n (NH3) / n (H2) ]

HCO
-8.0
-8.0
-10.0
H2CN
HCO+
-8.0
-9.0
-11.0
-10.0
-9.0
-10.0
-12.0
-11.0
-10.0
-13.0
-12.0
-14.0
-12.0
-11.0
a
a
-15.0
-14.0
-13.0
-12.0
b
-16.0 ba
-14.0
-13.0d
-16.0-17.0 db
-15.0
-14.0g-18.0 g
d
-16.0
-18.0
-19.0
-15.0
1
10
100
1000
1
10
100 g
1000
-20.0
-7.0
-16.0
-8.0
1 a
10
100
1000
1
10
100
1000
OCS
SO
-9.0
-10.0
b
-7.0
+
-11.0
N
H
-12.0
2
d
Figure 5. Average abundances of selected
species calculated
for the RHD model perpendicular to the OCS
rotation axis.

a
b
d
g

)]

-9.0

(H2) ]2) ]n (H2) ]分⼦流やenvelopeを横からみたとき (横軸の単位はAU)

-9.0

a

bはFC形成後2000年, dはFC形成後4000年(2nd collapse直前)

Figure 17. Evolution of the SO (87 76 ) line profile in the MHD
model for RHD model (top) and MHD model (bottom). Eachspectrum is averaged over pixels within a central 0.35 arcsec
diameter
circular beam. Solid lines: i = 0 , dashed lines:
A.(53
K.AU)
Young
et al.
i = 90 .

Figure 19. HCO+ (1 0) Integrated intensity
Middle: RHD model; Bottom: MHD model

MHDのみ
分⼦流軸⽅向

RHDのみ

MHDのみ

17. Evolution of the SO (87

Figure 23. Synthetic ALMA spectra of the MHD model in the
Figure 23.
Synthetic
spectrathe integration,
MHD model iface–on direction.
Top:SO (87ALMA76 ), 8hours of
total
face–on direction.
Top: SO(87 76 ), 8spectra
hours total
out14 configuration.
It is possible
to distinguishearly integr
and
out14
configuration.
is outflow
possible develops.
to distinguishspectraearllate in FHSC
stage
at i = 0 asItthe
Bottom:MHD
late+in(1-0),
FHSC8stage
attotali = 0 integration,
as the outflow
develops.
model HCO
hours
0.21
⇥ 0.25 Bottom:
arcsec+HCOgrows(1-0),
8 hours
total integration,
0.21 ⇥ and0.25 abeam. Asmodel
the FHSC
emissionincreasesat the line centre
beam.
the+ FHSC
grows
emissionincreases
atdiscernible
the line
centr
Figure19. AsHCO(1the0)blue
Integrated
intensity
snapshot(d).
becomesbrighter
than
peak
and
this
isforstill+
Figure20.theHCO(10) andspectrafor
RHDdiscermodbecomesbrighter
than
blue
peak
this
is
still
Middle:RHD
model;
Bottom:MHD
model76 ) line profile in the MHD in the synthetic interferometric spectra.
model(bottom) within
a 0.67 arcsec aperturein the synthetic
interferometric
spectra.