星形成ゼミ資料
Dec.8, 2017 古屋玲（徳島⼤学・教養教育院）

THE FLUXCOMPENSATOR: MAKING RADIATIVE TRANSFER MODELS OF HYDRODYNAMICALSIMULATIONS DIRECTLY COMPARABLE TO REAL OBSERVATIONS
Christine M. Koepferl1,2 and Thomas P. Robitaille1,3

1 Max Planck
Christine
M.for
Koepferland Thomas
P. RobitailleInstitute
Astronomy,Königstuhl
17, D-69117 Heidelberg, Germany
2 University Observatory Munich, Scheinerstr. 1, D-81679 Munich, Germany

radiative transfer codese. g. HYPERION

3 Freelance Consultant, Headingley Enterprise and Arts Centre, Bennett Road Headingley, Leeds LS6 3HN(Received 22nd June 2017; Accepted 14th August 2017)
The Astrophysical Journal

ABSTRACT
When modeling astronomical objects throughout the universe, it is important to correctly treat thelimitations of the data, for instance finite resolution and sensitivity. In order to simulate these e↵ects,
and to make radiative transfer models directly comparable to real observations, we have developed
an open-source Python package called the FluxCompensator that enables the post-processing of
the output of 3-d Monte-Carlo radiative transfer codes, such as Hyperion. With the FluxCompensator, realistic synthetic observations can be generated by modeling the e↵ects of convolution with
arbitrary point-spread-functions (PSFs), transmission curves, finite pixel resolution, noise and reddening. Pipelines can be applied to compute synthetic observations that simulate observatories, suchas the Spitzer Space Telescope or the Herschel Space Observatory. Additionally, this tool can read inexisting observations (e. g. FITS format) and use the same settings for the synthetic observations. In
this paper, we describe the package as well as present examples of such synthetic observations.
FLUXCOMPENSATOR

observing device (telescope, detector, etc.)

interstellar/intergalactic
mediumastronomical object

2

1. INTRODUCTION

Theoretical simulations (e. g. hydrodynamical simulans) are currently used to reproduce the physical proses in astronomical objects, such as young stellar obs (YSOs), star-forming regions and galaxies. Howr, we cannot use these simulations directly to pret or explain observational features of such objects.
tead, full radiative transfer calculations need to beformed to properly take into account e↵ects, such as
al temperature variations or changes in the opticalth. These radiative transfer calculations provide ideed observations, which are somewhat closer to real obvations than the hydrodynamical simulations. In thejects but rather provides a generic set of geometries. Hyperion as other dust radiative transfer
codes, treats all dust interactions, such as absorption, re-emission and scattering. Hyperion is an
LTE code, and, in addition, can treat (simple) nonLTE approximations, such as the emission of polycyclic aromatic hydrocarbons (PAHs). It uses the
Lucy method, the modified random walk, peelingo↵ and raytracing.
• MOCASSIN
Mocassin, a 3-d photoionization and dust nonLTE Monte-Carlo radiative transfer code (for de-シュミレーション出⼒

FluxCompensator
PACSThe70mic透過曲線を畳み込み

⾚化補正

5

今後の開発計画は,• PSFの波⻑依存性も実装したい。
• ユーザーが指定した装置でのノイズレベ
ル評価機能など(ソースとの分離⼯程も含めるの意味か？)• 検出器がサチッたバ
ンドにも対応させる。
• Gas-radiativetransfer postprocessing (何のことかよくわからなかった)• ⼲渉計の空間フィル
ター機能をCASAとの連携で実装
PACS 70mic PSFを畳み込み
PACS 70 mic のpixel sizeでなまして、ノイズ付加• astropyのaﬃliatedFig. 3.— Snapshots during the ongoing pipeline: a) Hyperion output image slice closest to = 70 µm, b) image after convolution withACS 70 µm transmission curve, c) PACS 70 µm image after reddening, d) PACS 70 µm image with corresponding pixel resolution, e) PACS
package群に加えた
0 µm image after PACS PSF convolution, f) PACS 70 µm with Gaussian noise.
い。
he synthetic image, where Figure 3(a) was weighted
tion of energy. See Table 2 in Appendix A for ↵ andwith a PACS 70 µm filter from the FluxCompensator
atabase. The database currently has 24 transmision curves from Two Micron All-Sky Survey (2MASS),

values of the filters in the database and their references.
To weight the synthetic data with the spectral transmission, the FluxCompensator calculates the e↵ec-

INSIGHTS FROM SYNTHETIC STAR-FORMING REGIONS:
I. RELIABLE MOCK OBSERVATIONS FROM SPH SIMULATIONS
Christine M. Koepferl1,2 , Thomas P. Robitaille1,3 , James E. Dale4 , and Francesco Biscani1
1 Max Planck Institute for Astronomy, Königstuhl 17, D-69117 Heidelberg, Germany

2 Scottish Universities Physics Alliance (SUPA), School of Physics and Astronomy, University of St Andrews

North Haugh, St Andrews, KY16 9SS, UK

3 Freelance Consultant, Headingley Enterprise and Arts Centre, Bennett Road Headingley, Leeds LS6 3HN
4 University Observatory Munich, Scheinerstr. 1, D-81679 Munich, Germany(Received 29 January 2016; Accepted 1 March 2016)
The Astrophysical Journal Supplement

ABSTRACT
Through synthetic observations of a hydrodynamical simulation of an evolving star-forming region,
we assess how the choice of observational techniques a↵ects the measurements of properties which
trace star formation. Testing and calibrating observational measurements requires synthetic observations which are as realistic as possible. In this part of the paper series (Paper I), we explore di↵erent
techniques for how to map the distributions of densities and temperatures from the particle-based
simulations onto a Voronoi mesh suitable for radiative transfer and consequently explore their accuracy. We further test di↵erent ways to set up the radiative transfer in order to produce realistic
synthetic observations. We give a detailed description of all methods and ultimately recommend techniques. We have found that the flux around 20 µm is strongly overestimated when blindly coupling
the dust radiative transfer temperature with the hydrodynamical gas temperature. We find thatwhen instead assuming a constant background dust temperature in addition to the radiative transfer
heating, the recovered flux is consistent with actual observations. We present around 5800 realistic
synthetic observations for Spitzer and Herschel bands, at di↵erent evolutionary time-steps, distancesand orientations. In the upcoming papers of this series (Paper II, Paper III and Paper IV), we will
test and calibrate measurements of the star-formation rate (SFR), gas mass and the star-formation
efficiency (SFE) using our realistic synthetic observations.
1. INTRODUCTION

the last decade, progress in the field of star forhas been fueled by an increasing number of large/optical surveys of our Galaxy (e. g. from space

温度

温度

密度

シュミレーション
エッジからの動径距離Fig. 2.— Temperature relation in density and radius of SPH particles and a 2-d projection showing the temperature increase at the edgesof the simulation at time-step 122 (6.109 Myr).

2.2.1. Method
c1 — Neutral
SPH Particles
which
lie withinfrom
ourSPHselectedbox of 30 pc width. 7 The
Synthetic
Star-forming
Regions in- I.Box
Reliable Mock
Observations
Simulations

While the D14 SPH simulations have high-mass stellar
feedback mechanisms implemented,
the simulations do
フィードバック
not take into account the dust heating from the stars.
In D14 SPH simulations, the external
heating throughを⼊れた時刻the di↵use Galactic radiation field is approximated by
the Larson (2005) EOS8 rather than just an isothermal
approximation. The EOS contains the following regimes(following Equation 1 and Figure 1 of Dale et al. 2012):• line cooling processesP ⇠ ⇢3/4(6)

8 ⇢ < 5.5 ⇥ 10 19 g cm 3

• dust cooling processesP ⇠ ⇢1.0(7)

red points are (partly or fully) ionized particles (c. f. with
Dale et al. 2012, 2007) within the box. The horizontal
フィードバックあり
line at 104 K shows all the fully ionized SPH particles.
The other red points are partly ionized and lie within theselected box. The black points are neutral SPH particles
in the selected box. The state of these neutral particles
is described by the EOS which is describing the dependency of density and temperature (diagonal line: Eq. 6,
Eq. 8, Eq. 9; horizontal: Eq. 7).
In this work, we will use a dust radiative transfer code,
therefore we will only use SPH particles in the box which
are completely neutral since dust gets destroyed during
the ionization process (Diaz-Miller et al. 1998). Hereafter, we call the clipping for size and only neutral SPH
particles, resulting in the black particles in Figure 2(a),method c1.

2.2.2. Method c2 — Neutral SPH Particles in Box &8 5.5 ⇥ 10
g cm < ⇢ < 5.5 ⇥ 10
g cm ,
Temperature Limit
Fig. 3.— Temperature histogram of neutral SPH particles for time-step 024 (3.576 Myr), 032 (4.768 Myr), 072 (5.364 Myr) and 122
19

3

15

3

neutral SPHの時間進化(右パネルへ）ごとのTの頻度分布(6.109 Myr). Feedback gets switched on after time-step 32. The sharp spikes In
at 7.5
K are 2(b),SPH particles
the that
isothermal
phase.temperature
Figurewe caninsee
the SPH

INSIGHTS FROM SYNTHETIC STAR-FORMING REGIONS:
II. VERIFYING DUST SURFACE DENSITY, DUST TEMPERATURE & GAS MASS MEASUREMENTS
WITH MODIFIED BLACKBODY FITTING
Christine M. Koepferl1,2 , Thomas P. Robitaille1,3 , and James E. Dale4
1 Max Planck Institute for Astronomy, Königstuhl 17, D-69117 Heidelberg, Germany

2 Scottish Universities Physics Alliance (SUPA), School of Physics and Astronomy, University of St Andrews

North Haugh, St Andrews, KY16 9SS, UK

3 Freelance Consultant, Headingley Enterprise and Arts Centre, Bennett Road Headingley, Leeds LS6 3HN
4 University Observatory Munich, Scheinerstr. 1, D-81679 Munich, Germany(Received 29 January 2016; Accepted 23 June 2016)
The Astrophysical Journal

ABSTRACT
We use a large data-set of realistic synthetic observations (produced in Paper I) to assess how
observational techniques a↵ect the measurement physical properties of star-forming regions. In thispart of the paper series (Paper II), we explore the reliability of the measured total gas mass, dust
surface density and dust temperature maps derived from modified blackbody fitting of synthetic
Herschel observations. We found from our pixel-by-pixel analysis of the measured dust surface densityand dust temperature a worrisome error spread especially close to star-formation sites and low-density
regions, where for those ”contaminated” pixels the surface densities can be under/over-estimate-0.3(e)-236.0.6(y)-09.96(

⻩⾊の星：降着中の原始星,
⻘で⽩抜きの星：電離ガスのバブル内の星;
⽩コントア：Chi^2>5 temperatures. U
Fig. 6.— Diagnostic
maps of measured and intrinsic
values for dust surface density, dust meandensities and dust
wo 対
rows for distance D1 at 3 kpc and lower two rows for distance D2 at 10 kpc. Yellow stars are the projected positions of young, acc
tars応while blue stars are stars inside the ionized bubble. 2 > 5-contours are plotted in white. Displayed is the result for versio
ersion R2 will be available as分⼦雲全体
online figure.
星形成サイト星形成サイトの周辺領域

f the transition from gas to stars and can be expressed

munity,hence for entire
galaxies,
assumecharacy the star-formation
efficiency(SFE) atanda certaintime
FROM
SYNTHETIC
STAR-FORMING REGIONS:
tic timescales of INSIGHTSt⇤ = 100 Myr(Kennicutt& Evans
III. CALIBRATIONOF MEASUREMENTTECHNIQUESOF STAR-FORMATION RATES; Murphyet al. 2011; Calzetti
et al. 2010; Rieke
et al.
M⇤ (t)
SFE(t) = have not yet. been
(1) suffi). These techniques
testedM⇤M.
(t) +Koepferl
Mgas (t) 1,2 , Thomas P. Robitaille1,3 , and James E. Dale4
Christine
ly for (smaller,
local) individual regions. Therefore, Synthetic Star-forming Regions - III. Calibration of Measurement Techniques of Star-formation Rates
5
1 Max
Planck Institute for Astronomy, Königstuhl 17, D-69117 Heidelberg, Germany
Another
importantphysicalproperty
defininga instar2 Scottish
nalyzethereliabilityofthese
techniques
this secUniversities
Physics
Alliance(SUPA),
Schoolof Physics and Astronomy, University of St Andrews
ormingregion
is
the
rate
at
which
it
forms
stars.
The
North Haugh, St Andrews, KY16 9SS, UK
銀河研究で⽤いられているSFR
3 FreelanceverageSFR between
times t0Headingley
and t is definedas: and Arts Centre, Bennett Road Headingley, Leeds LS6 3HNConsultant,Enterprise
Figure 1, we plot
the intrinsic
instantaneous
SFR di4 University
Observatory Munich, Scheinerstr. 1, D-81679 Munich, Germanyy from the simulation
solid
curve),
we callM⇤ (t) M⇤(red(t0 )(ReceivedM⇤ (t)⇤ (twhich
0)29MJanuary
2016;
Accepted 8 March 2016)
SFR(t) ==, (2)t) from
wasAstrophysicalcalculated Journal
t t0now on, sincet⇤it The
sim ( t⇤ =the star-formation timescale t⇤ equal
to over
the timewhere t⇤ = t t0 is the characteristic
timescaleABSTRACT
of には、特徴的な時間尺度があり、これは
thet. There
To compare
it to
the meawhich
the simulation
SFR is measured.
exists a large
varietyextensive
set ofEachrealistic
synthetic
observations (produced in Paper I), we assess in thisdf Throughtechniques,whichcharacteristic
timescales
tracers
to anmeasurethe follow
SFR.
techniqueprobes
通常、δt*=
100
Myrが仮定される。
of100
the Myr,paper
series(Paper
III) howthechoiceof observational techniques a↵ects the measurementtarformation
withinaconvolve
certaincharacteristic
timescale
we
(see
Press
et
al. 1992)⇤art=tf⇤ ,star-formation
tracerto tracer.
can divide
⼀⽅、銀河スケールのSFRの計算にはlocal
rates
in Onestar-forming
regions. We test the accuracy of commonly used(whicht⇤ =variest)fromwitha (SFRs)
normalized
top-hatfunctionsimhedi↵erent
techniques
intotwo
categories:2
echniques
andconstructnewmethods
to extract
Bronstein
et al.
2005) ofwidth
100 Myr. Thethere-SFR, so that these findings can be applied toな星形成が考慮されることはほとんどない。
measurethe SFRin real regions
throughout
thehereafter
Milky Way. We investigate di↵use infrared SFR tracers
ng1.なので、あいだをつなぐ研究が必要。
convolved
function(red dashedcurve),
Indirect Techniquesuch ˆasthoseusing 24starsµm,ionize
70 µmtheandtotal infraredemission, which have been previously calibrated forYounghigh-massmaterialarounddlobal
SFR
(
t
=
100
Myr),
is
aroundthreeorderssimwithscales.
⇤ X-ray andgalaxyWe setup a toy(UV)
model
of a galaxyof
and show that the infrared emission is consistentthem
ultra-violet
radiation.
nitude
lower
since
t/
t
=
149
000
yr/100Myr⇡⇤ tracewith This
the intrinsic
SFR
using
extra-galactic
calibrated
laws (although the consistency does not prove
process canbe usedtostar formationin
観測的銀河研究で、SFRを推定するとき、Fig. 1.— Averaged
measured SFRs
with the techniques
SFR24, SFR70
and SFRIR using the total flux at 24 µm, 70 µm and the total
【上の図】
.
heir
reliability).
For local scales, we show that these techniques
producecompletely
unreliable
results
infrared flux, respectively, in comparison with the intrinsic SFR.
or 代表的な⽅法として、24mic,
single star-forming regions, which
are governedbydi↵erentcharacteristic
timescales.
70mic,
andfromthe simulation
averagedovera similar
timescale We
Weshowderive the total luminosity in the infrared under the
⾚は「δt*=
100
Myrの瞬間」ごとのintrinsicなSFR(計算値)、
ˆ
assumption
ow calibrationofthese
techniques
canbeimprovedforsinglestar-forming
regionsbyadjustingthe that our 7 bands from ⌫0 (IRAC 8 µm) to ⌫1
2.1. Method SFR24 — 24 µm Tracer
SFRsim ( t⇤ = 100 Myr). However, they should match if
total
IR光度から求めらることが多い。例え500 µm) cover most of the infrared emission. We
thetechniqueworked ofaccurately.
We provide the meaharacteristic timescale and the scaling factor and givesuggestions
new calibrationsof the (SPIRE
di↵use星印は背景光なしの場合、
integrate(using Simpson’s rule, see Bronstein et al. 2005,
sured
values
in
Table
1
of
AppendixA.
eke
et al. (2009)
empirically calibrated a SFR-totar-formation
tracers. We show that star-forming regions that are dominated by high-mass Press
stellar
et al. 1992) to calculate the flux density F⌫corr ( ):ば、24mic光度とのスケール則は、
nosity relationof 24a µm
emission:
Method SFR70 — 70 µm TracerZ⌫eedback
experience
rapid
drop in infrared emission下部の⾚い点線が銀河全体での換算値：通常観測されるSFR
once2.2.high-massstellar feedback is turned on,
1

Based on Spitzer data, Calzetti et al. (2010) empirL(IR) = d⌫F⌫corr ( )
(8)
which
impliesdi↵erentcharacteristic
timescales.
Moreover,we
explorethemeasuredSFRs
calculated
icallycalibrated
the
SFR
using
the
total
luminosityin
SFR24 µm ( t⇤ = 100 Myr)
43
42
1
⌫0
70
µm findabove ⌫L
1.4 ⇥ 10 ergspoints := 2.03⇥ 10population.
⇥(5)⌫ (70 µm)
irectly from the observed
youngstellarWe
thatthe measured
sources follow
1
M yrpace of star formation more directly thanWith the total infrared luminosity L(IR), we use Eq. 7
SFR
100 Myr)he evolutionarydi↵usestar-formation
tracers.(6)70 µm ( t⇤ == 5.88 ⇥ 10 44 ⇥to derive the SFRIR ( t⇤ = 100 Myr) for our set of syn1
M yr
⌫L⌫ (24 µm)
thetic observations. Again, in Figure 1 (right panel),⌫L⌫ (70 µm)⇥,
we compare the measurements (black) with the values
1
⇥
1 .
ergs
s
thethesimulations(red). Again the measurements
ergs s
1. INTRODUCTION
The measured star-formation
propertiesfrom
areobserwith combined background (black dots) are higher due
Notethat here
we will use PACS
70 µm to estimate
theand tovationalbenchmarkbetweenreality
Examthe complex
background. The o↵set is comparable
42
1
ar-forming
region,is transformed
h they foundto begasvalid
within 2.3 ⇥into
10 stars.
ergs s SFR
 70 µm ( t⇤ = 100 Myr) with Eq. 6 and again test the totheory.
thetechnique
SFR70 because Herschel bands withples ofincludegas mass
Mgas andSFR70
for properties
luminosities below
1039 ergsthes 1 . total
We show
e24dominating
mechanisms
and the im-we willhigh background have been used (see Paper II). For theµm)  5.0driving⇥ 1043 ergs
s 1 . Nevertheless,thein Figure1 (black
midtheresults
totalagain
stellarmass
M⇤symbols
in a instar-forming
region measurements(asbackground-free(black stars), the simuhe
stellararewithstill a total
hisdi↵erent
techniquefor feedback
syntheticmechanisms
observations
dle panel). We can see that the di↵erence between thelationspreadisagainlarger,which
is due to the MIR
summarized
Dalebackgroundet al. 2014).
The property whichmeasurementwith bycombined(black dots)
40 are di↵erent
1
nless
thethanfield, there
approaches
amongparts of the total infrared flux (see SFR24 and SFR70).
10 ergs s
to test whether the SFR24
andthe background-freemeasurements(black gas
stars)massisconnects
stellar
mass
and
total
is
the efficiency
Overall,
the technique SFRIR produces values of thers to addresstheseforquestions:
In the observathan for the method SFR24. This is again due
nique
can be usedsingle star-forming
regions. larger
SFR
(
t
IR
⇤ = 100 Myr) which are up to three orders ofof thethehightransition
gas toplanestars
tobackground fromof the Galactic(as andex- can be expressed
proach,thecharacteristic
observable
properties
ˆ sim ( t⇤ =magnitude higher than the simulated rate SFR
ing Eq. 5, we derive the SFR24 µm ( t⇤ = 100 Myr)plainedin PaperII) and a result ofefficiency
the non-uniform
back- at a certain timebythestar-formation(SFE)
100 Myr). We provide the measured values in Table 1 ofrming
are used
to study andcompareground during background correction. Again, similarly
ur set regionsof synthetic
observations
in MIPS
24 µm forAppendix A.
t:to what we found for the SFR24 technique, the mearegions in an attempt to disentangle the domi-に相当
縦の点線はmassive starのfeedbackを⼊れた時刻【結果】「δt* ごとの瞬間SFR」値と「通常観測されるSFR」の３桁のギャップ：「They should match if the technique worked
accurately.」としか説明ない。
３つのSFR推定法に顕著な差なし。

Fig. 2.— Instantaneous rate SFRsim ( t⇤ = t) from the simulation (red dots), rate convolved with di↵erent characteristic timescales t⇤(yellow to red) and the scaled measured emission (black stars). Numbers on top represent the time-step IDs and vertical lines highlightthe time-step of the simulations when feedback was switched on.

3.2.

1∗
MikitoKOHNO, Kazufumi
TORII2 , Kengo TACHIHARA1 , Tomofumi
Publications of the Astronomical Societyof Japan,(2014), Vol.
00, No. 0
2,3
2,3
1

FOREST Unbiased Galactic plane Imaging
survey with the Nobeyama 45-m telescope(FUGIN) : Molecular clouds toward W33 ;possible evidence for a cloud-cloud collision
triggering O star formation

Mikito KOHNO1 ∗ , Kazufumi TORII2 , Kengo TACHIHARA1 , Tomofumi
U MEMOTO 2,3 , Tetsuhiro M INAMIDANI2,3 , Atsushi N ISHIMURA1 , Shinji F UJITA1 ,Mitsuhiro M ATSUO2 , Mitsuyoshi YAMAGISHI4 , Yuya T SUDA5 , Mika K URIKI6 ,Nario K UNO6 , Akio O HAMA1 , Yusuke H ATTORI1 , Hidetoshi S ANO1,7 , Hiroaki
YAMAMOTO1 and Yasuo F UKUI1,7
1

arXiv:1706.07964v4 [astro-ph.GA] 25

20

U MEMOTO , Tetsuhiro M INAMIDANI , Atsushi N ISHIMURA , Shinji F UJITA1 ,Mitsuhiro M ATSUO2 , Mitsuyoshi YAMAGISHI4 , Yuya T SUDA5 , Mika K URIKI6 ,Nario K UNO6 , Akio O HAMA1 , Yusuke H ATTORI1 , Hidetoshi S ANO1,7 , Hiroaki
YAMAMOTO1 and Yasuo F UKUI1,7
1

Department of Physics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601,
Japan
2
Nobeyama Radio Observatory, National Astronomical Observatory of Japan (NAOJ),
National Institutes of Natural Sciences (NINS), 462-2, Nobeyama, Minamimaki, Minamisaku,
smaller than that in the larger cloud. Different color components in the collision schematics correspond to the different color on the
Nagano384-1305, Japan
Fig. 11. Schematic of a collision between two dissimilar clouds and cartoons position-velocity diagrams, where the gas density in3
Department of Astronomical Science, School of Physical Science, SOKENDAI (TheP
Graduate University for Advanced Studies), 2-21-1, Osawa, Mitaka, Tokyo 181-8588, Japan
4
Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency,
Smaller cloud
Chuo-ku, Sagamihara 252-5210, Japan
5
Graduate School of Science and Engineering, Meisei University, 2-1-1 Hodokubo, Hino,
Tokyo 191-0042, Japanv
6
Department of Physics, Graduate School of Pure and Applied Sciences, University ofTsukuba, 1-1-1 Ten-nodai, tsukuba, Ibaraki 305-8577, Japan
7
Institute for Advanced Research, Nagoya University, Furo-cho, Chikusa-ku, Nagoya
Larger cloudV
464-8601, Japan

Department of Physics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601,
Japan
2
Nobeyama Radio Observatory, National Astronomical Observatory of Japan (NAOJ),
National Institutes of Natural Sciences (NINS), 462-2, Nobeyama, Minamimaki, Minamisaku,
Nagano 384-1305, Japan3
Department of Astronomical Science, School of Physical Science, SOKENDAI (The
Graduate University for Advanced Studies), 2-21-1, Osawa, Mitaka, Tokyo 181-8588, Japan
4
Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency,
Chuo-ku, Sagamihara 252-5210, Japan
5
Graduate School of Science and Engineering, Meisei University, 2-1-1 Hodokubo, Hino,
Tokyo 191-0042, Japan
6
Department of Physics, Graduate School of Pureand Applied Sciences, University ofFig. 5. R3−2/1−0 maps of the (a) 35 km s−1 , (b) 45 km s−1 , and (c) 58 km s−1 cloud are shown in the color image. Contours show the 12 COJ =3–2
Tsukuba, 1-1-1 Ten-nodai, tsukuba, Ibaraki ′′305-8577, Japan∗ E-mail: mikito@a.phys.nagoya-u.ac.jp
emission, which was spatially smoothed 50 . Plots are same as Figure 1. The clipping level is 8σ.
7Institute for Advanced Research, Nagoya University, Furo-cho, Chikusa-ku, Nagoya
464-8601, Japan
Received 2017 July 10; Accepted 2017 October 25

Received 2017 July 10; Accepted 2017 October 25

Abstract

Turbulent layer

Abstract
thgis-fo-enil∗ E-mail: mikito@a.phys.nagoya-u.ac.jp

We observed molecular clouds in the W33 Larger
high-masscloud star-forming region associated with com-

Fig. 2. (a) Longitude-velocity diagram of the 12 COJ = 1–0 emission with the NANTEN2 telescope. The arrows indicate the four velocity components ; black
pact and extended H II regions using the NANTEN2 telescope as well as the Nobeyama 45We observed molecular
clouds in the W33 high-mass star-forming region associated with comarrow : 10 km s−1 , blue arrows : 35 km s−1 , 45 km s−1 , and 58 km s−1 . The dashed box indicate the W33 region as the range of Figure 2(b). (b)pact and extended H II regions using
the
NANTEN2
telescope
as
well
as
the
Nobeyama
45m telescope in the J =1–0 transitions of 12 CO, 13 CO, and C18 O as a part of the FOREST
Longitude-velocity diagram of the 12 COJ =1–0
emissionwith 18
the Nobeyama 45-m telescope. The arrows indicate the three velocity components ; 35 km s−1 ,12
13
m telescope
in
the
J
=1–0
transitions
of
CO,
CO,
and
C
O
as
a
part
of
the
FOREST−1
Unbiased Galactic plane Imaging survey with the Nobeyama 45-m telescope(FUGIN)legacy45 km s−1
, and 58
km sImaging
.
Dense
layer
Unbiased
Galacticplane
survey with the Nobeyama 45-m telescope (FUGIN) legacy

survey. We detected three velocity components at 35 km s−1 , 45 km s−1 , and 58 km s−1 . The
35 km s−1 and 58 km s−1 clouds are likely to be physically associated with W33 because of
the enhanced 12 CO J = 3–2 to J =1–0 intensity ratio as R3−2/1−0 > 1.0 due to the ultraviolet
irradiation by OB stars, and morphological correspondence between the distributions of molecular gas and the infrared and radio continuum emissions excited by high-mass stars. The two
clouds show complementary distributions around W33. The velocity separation is too large
to be gravitationally bound, and yet not explained by expanding motion by stellar feedback.
c 2014. Astronomical Society of Japan.
⃝

survey. We detected three velocity components at 35 km s−1 , 45 km
s−1 , and 58 km s−1 . The
Cavityv−1−1
35 km s and 58 km s clouds are likely to be physically associated with W33 because of
the enhanced 12 CO J = 3–2 to J =1–0 intensity ratio as R3−2/1−0 > 1.0 due to the ultraviolet
Smaller cloud
irradiation by OB stars, and morphological correspondence between the distributions of molecular gas and the infrared and radio continuum emissions excited by high-mass stars. The twoof the Astronomical
Societyaroundof Japan,(2014),
00, No.separation
0
clouds Publications
show complementary
distributions
W33.
TheVol.
velocity
is too large
to be gravitationally bound, and yet not explained by expanding motion by stellar feedback.
c 2014. Astronomical Society of Japan.
⃝コントアは野辺⼭で取得C18O(1-0)

Fig. 16. Velocity channel map of the C18 O J = 1–0 emission with a velocity step of 4.0 km s−1 with Nobeyama. Contours show the VLA 90 cm radio
continuum image. Plots are same as Figure 1.