Star Formation Newsletter#345 30-35
Yusuke Aso (KASI)

Quantifying Variability of YSOs in the Mid-IR Over Six Years with NEOWISE

iv:2107.10751v1 [astro-ph.SR] 22 Jul 2021

Wooseok Park,1 Jeong-Eun Lee,1 Carlos Contreras Peña,2, 3 Doug Johnstone,4, 5 Gregory Herczeg,6, 7
Sieun Lee,1 Seonjae Lee,8 Anupam Bhardwaj,9 and Gerald H. Moriarty-Schieven4
1 School of Space Research, Kyung Hee University, 1732, Deogyeong-daero, Giheung-gu, Yongin-si, Gyeonggi-do 17104, Korea:jeongeun.lee@khu.ac.kr
2 Centre for Astrophysics Research, University of Hertfordshire, Hatfield AL10 9AB, UKØ 変光ß突発降着、⾮⼀様円盤、磁気リコネクション、連星、星の⾼温低温スポット。
3 School of Physics, Astrophysics Group, University of Exeter, Stocker Road, Exeter EX4 4QL, UK
4 NRC Herzberg Astronomy and Astrophysics, 5071 West Saanich Road, Victoria, BC, V9E 2E7, CanadaØ 全光度の変化を⾒るならMIR。
5 Department of Physics and Astronomy, University of Victoria, 3800 Finnerty Road, Elliot Building, Victoria, BC, V8P 5C2, Canada
6 Kavli Institute for Astronomy and Astrophysics, Peking University, Yiheyuan 5, Haidian Qu, 100871 Beijing, China
Ø Spitzer+WISEで若いほど原始星アウトバーストが多い。
7 Department of Astronomy, Peking University, Yiheyuan 5, Haidian Qu, 100871 Beijing, China

8 Department 0-IIIをカバー。
Ø NEOWISEはClassof Physics and Astronomy, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Korea
9 Korea Astronomy and Space Science Institute (KASI), 776 Daedeokdae-ro, Yuseong-gu, Daejeon 34055, Korea
Ø 0.5-6.5年のMIR測光。1700/5400 (原始星0/Iの55%, 円盤IIの33%, PMS IIIの15%)で変光。

ABSTRACT

Variability
in young(YSOs)
be 22
causedby various
time-dependent phenomena
Ø WISE W1 3.4µm, W2
4.6 stellarµm, objects
W3 12µm,canW4µm,6"-12"、
associated with star formation, including accretion rates, geometric changes in the circumstellar disks,
半年毎に全天,14 hydromagnetic
epochs, 10-20
exposures/epoch~数⽇。
stochastic
interactions
between stellar surfaces and inner disk edges, reconnections
withinthestellarmagnetosphere,andhot/coldspots on stellar surfaces. We uncover and characterizeØ W2で5 epochs, 不確かさ<0.2 mag, 座標が既知のものを使⽤。
⇠1700 variables from a sample of ⇠5400 YSOs in nearby low-mass star-forming regions using mid-IR
light curves obtained from the 6.5-years NEOWISE All Sky Survey. The mid-IR variability traces awide range of dynamical, physical, and geometrical phenomenon. We classify six types of YSO mid-IR
Ø 光度曲線全体の偏差/観測誤差>3、ΔW2>3σで変光同定。
variability based on their light curves: secular variability (Linear, Curved, Periodic) and stochastic
Ø Lomb-Scargle
periodogram解析と線形フィット。
variability(Burst, Drop, Irregular ). YSOs in earlier evolutionary stages have higher fractions ofvariables and higher amplitudes for the variability, with the recurrence timescale of FUor-type outbursts
Ø 雑⾳による誤報率が線形フィットで低ければ線形。
(defined here as W1 or W2 > 1 mag followed by inspection of candidates) of ⇠1000 years in the early
embedded protostellar phase. Known eruptive young stars and subluminous objects show fractions ofvariables similar to the fraction (⇠ 55%) found in typical protostars, suggesting that these two distinct
types are not distinct in variability over the 6.5-year timescale. Along with brightness variability, we
also find a diverse range of secular color variations, which can be attributed to a competitive interplay
between the variable accretion luminosity of the central source and the variable extinction by material
associated with the accretion process.

iv:2107.10751v1 [astro-ph.SR] 22 Jul 2021

Quantifying Variability of YSOs in the Mid-IR Over Six Years with NEOWISE
Wooseok Park,1 Jeong-Eun Lee,1 Carlos Contreras Peña,2, 3 Doug Johnstone,4, 5 Gregory Herczeg,6, 7
Sieun Lee,1 Seonjae Lee,8 Anupam Bhardwaj,9 and Gerald H. Moriarty-Schieven4
1 School of Space Research, Kyung Hee University, 1732, Deogyeong-daero, Giheung-gu, Yongin-si, Gyeonggi-do 17104, Korea:jeongeun.lee@khu.ac.kr
2 Centre for Astrophysics Research, University of Hertfordshire, Hatfield AL10 9AB, UK
3 School of Physics, Astrophysics Group, University of Exeter, Stocker Road, Exeter EX4 4QL, UK

4 NRC Herzberg
Secular(29%)
Astronomy and Astrophysics, 5071 West Saanich Road, Victoria, BC, V9E有意な線形
2E7, CanadaØ Linear。負勾配が多いà減衰が⻑い。
6 Kavli Institute for Astronomy and Astrophysics, Peking University, Yiheyuan 5, Haidian Qu, 100871 Beijing, China
7 Department of Astronomy, Peking University, Yiheyuan 5, Haidian Qu, 100871 Beijing, China
Ø 周期3-13年Curved。
8of Physics and Astronomy, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, KoreaØDepartment
周期<3年Periodic。
9 Korea
Astronomy and Space Science Institute (KASI), 776 Daedeokdae-ro, Yuseong-gu, Daejeon 34055, Korea

5 Department of Physics and Astronomy, University of Victoria, 3800 Finnerty Road, Elliot Building, Victoria, BC, V8P 5C2, Canada

ABSTRACT
Variability in young stellar objects (YSOs) can be caused by various time-dependent phenomena
associated with star formation, including accretion rates, geometric changes in the circumstellar disks,stochastic hydromagnetic interactions between stellar surfaces and inner disk edges, reconnections
within the stellar magnetosphere, and hot/cold spots on stellar surfaces. We uncover and characterize⇠1700 variables from a sample of ⇠5400 YSOs in nearby low-mass star-forming regions using mid-IR
Stochastic
(71%)
light
curves obtained
from the 6.5-years NEOWISE All Sky Survey. The mid-IR variability traces a
wide
range
of
dynamical,
physical, and geometrical phenomenon. We classify six types of YSO mid-IRØ 数epochs（1-2年）のみの変光。
variability based on their light curves: secular variability (Linear, Curved, Periodic) and stochastic
Ø (med- min
orIrregularmax)>0.8(max- evolutionary stages have higher fractions ofvariability(Burst,Drop,). YSOs in earlier
variables
and Burst
higher amplitudes
for the variability, with the recurrence timescale of FUor-type outbursts
min)or Drop。
(defined here as W1 or W2 > 1 mag followed by inspection of candidates) of ⇠1000 years in the earlyØ その他Irregular。
embedded
protostellar phase. Known eruptive young stars and subluminous objects show fractions ofvariables similar to the fraction (⇠ 55%) found in typical protostars, suggesting that these two distinct
types are not distinct in variability over the 6.5-year timescale. Along with brightness variability, we
also find a diverse range of secular color variations, which can be attributed to a competitive interplay
between the variable accretion luminosity of the central source and the variable extinction by material
associated with the accretion process.

の続き

iv:2107.10751v1 [astro-ph.SR] 22 Jul 2021

Quantifying Variability of YSOs in the Mid-IR Over Six Years with NEOWISE
Wooseok Park,1 Jeong-Eun Lee,1 Carlos Contreras Peña,2, 3 Doug Johnstone,4, 5 Gregory Herczeg,6, 7
Sieun Lee,1 Seonjae Lee,8 Anupam Bhardwaj,9 and Gerald H. Moriarty-Schieven4

減光優勢
⾚化

の続き

1 School of Space Research, Kyung Hee University, 1732, Deogyeong-daero, Giheung-gu, Yongin-si, Gyeonggi-do 17104, Korea:jeongeun.lee@khu.ac.kr
2 Centre for Astrophysics Research, University of Hertfordshire, Hatfield AL10 9AB, UK

School of Physics, Astrophysics Group, University of Exeter, Stocker Road, Exeter EX4 4QL, UKØ 4 NRC
Dropは不均⼀物質、Burstは粘性磁気不安定。Linearは降着
Herzberg Astronomy and Astrophysics, 5071 West Saanich Road, Victoria, BC, V9E 2E7, Canada
5 Department of Physics
円盤上の降着率変化。Periodic,
Curvedは連星や内円盤
and Astronomy, University of Victoria, 3800 Finnerty
Road, Elliot Building, Victoria, BC, V8P 5C2, Canada （減光も増える）
6 Kavli Institute for Astronomy and Astrophysics, Peking University, Yiheyuan 5, Haidian Qu, 100871 Beijing, Chinaの⾼さ変化。
7 Department of Astronomy, Peking University, Yiheyuan 5, Haidian Qu, 100871 Beijing, China
8 Department of Physics and Astronomy, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Korea
Ø 変光率は0/I, II, IIIで55, 33, 15%、強度変化も0/Iが~3倍
9 Korea Astronomy and Space Science Institute (KASI), 776 Daedeokdae-ro, Yuseong-gu, Daejeon 34055, Korea
⾼いà進化。
戻る
ABSTRACTØ FUor/EXorの60%, VeLLO/LLSの55%à6.5年スケール
Variability
in young stellar objects (YSOs) can be caused by various time-dependent phenomena
点の⾊は時刻で同じ。
associated with star formation, including accretion rates, geometric changes in the circumstellar disks,
しかしFUor/EXorはほぼSecular、VeLLO/LLSは半分がstochastic
hydromagnetic interactions between stellar surfaces and inner disk edges, reconnections
within Irregularà変光では分けられない。
the stellar magnetosphere, and hot/cold spots on stellar surfaces. We uncover and characterize
内壁が降着
3⇠1700 variables from a sample of ⇠5400 YSOs in nearby low-mass star-forming regions using mid-IR
light curves obtained from the 6.5-years NEOWISE All Sky Survey. The mid-IR variability traces awide range of dynamical, physical, and geometrical phenomenon. We classify six types of YSO mid-IR
Ø FUor（ΔW1, W2>3 mag, linear）から間隔を推定。
variability based on their light curves: secular variability (Linear, Curved, Periodic) and stochastic
N: 観測原始星数=735-1059
variability(Burst, Drop, Irregular ). YSOs in earlier evolutionary stages have higher fractions ofvariables
and higher amplitudes
t: 観測時間=6.5
yr for the variability, with the recurrence timescale of FUor-type outbursts(defined here as W1 or W2 > 1 mag followed by inspection of candidates) of ⇠1000 years in the earlyk: 観測FUor数=2-9embedded
protostellar phase. Known eruptive young stars and subluminous objects show fractions ofvariables
similar to the fraction (⇠ 55%) found in typical protostars, suggesting that these two distinct
τ: バースト間隔
types are not distinct in variability over the 6.5-year timescale. Along with brightness variability, weτ=Nt/(k+1)=500-2300
yr
also
find a diverse range of secular color variations,which can be attributed to a competitive interplay
between the variable accretion luminosity of the central source and the variable extinction by material
associated with the accretion process.

戻る

内壁が⾼くなり⾚化

Susanne Pfalzner1, 2 and Amith Govind1, 2
1 Jülich Supercomputing Center, Forschungszentrum Jülich, 52428 Jülich, Germany
2 Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany

8.10296v1 [astro-ph.SR] 23 Aug 2021

4

Close stellar flybys common in low-mass clusters(Received August 24, 2021)

Submitted to ApJ

Ø ⼩質量星形成領域の円盤で渦状腕やtruncationがよく⾒つかるが、flybyは稀と思われていた。
Ø それを確かめるためにクラスターをシミュレーション。
ABSTRACT4‒3
Ø 環境的影響は>10/pcで顕著。半質量半径を1.3pcに固定すると>1000個の⼤質量クラスターで
Numerous
protoplanetary
discs show distinct spiralarms features. While possibly caused by arange of processes, detailed pattern∝Ranalysis1.6-2.1
pointsat close stellar flybys as cause for some of them.
のみ可能。しかし太陽系近傍ではMなので、⼩質量ほど⾼密度になる。
chmSurprisingly, these discs reside in young low-mass clusters, where close stellar flybys are expected to
be rare. This fact motivated us to take a fresh look at the frequency of close flybys in low-mass
Pfalzner & Govind
clusters. In the solar neighbourhood, low-mass clusters have smaller half-mass radii than their more
massive counterparts. We show that this observational fact results in the mean and central stellar
density of low-mass clusters being approximately the same as in high-mass clusters, which is rarelyreflected in theoretical studies. We perform N-body simulations of the stellar dynamics in young
clusters obeying the observed mass-radius relation. Taking the mean disc truncation radius as a proxy
for the degree of influence of the environment, we find that the influence of the environment on discs
is more or less the same in low- and high-mass clusters. Even the fraction of small discs (< 10 au) is
グレーより上は⻘がnearly identical. Our main conclusion is that the frequency of close flybys seems to have been severely
underestimated for low-mass clusters. A testable prediction of this hypothesis is that low-mass clusters
⾚より多い。
should contain 10%–15% of discs smaller than 30 au truncated by flybys. These truncated discs should
King W=9プロファ
be distinguishable from primordially small discs by their steep outer edge.

イル。

Keywords: protoplanetary disks — local interstellar matter

クラスター質量と半径の相関1. INTRODUCTION

The observed diversity in the structure of protoplanetary discs (for example, Andrews 2020; van der Marel et al.

8.10296v1 [astro-ph.SR] 23 Aug 2021

the complexity and numerical challenges the formation of massive clusters still faces. Thus the simulation sets only inwhen all stars are formed, and we take the mass-radius relations discussed in section 2 at face value. Naturally, any
Closestellarflybyscommonlow-massclusters
7 and gradually grows tohigh-massclustermust go througha phase where
it starts initially with a small number of stars
Close stellar
flybyscommonininlow-massclusters
its final size. However, such cluster formation simulations are still limited to clusters containing ⌧ 1000 stars, at the
1, ratio2
1, 2we revert
forthethe mass
resolution
required
here.Govind
Therefore,
to the close
often usedof following
each star is recorded, meaning forSusanneeach least
star,
and
periastron
distanceof eachflybymethod
is listed.
In the stellar dynamics
Pfalzner
and
Amith
in
the
fully
formedclusters
instead.
The
limitationsof
thisapproacharediscussedin
sectionthe second step, the truncation disc size for each of the close flyby events is calculated. Afterwards, we use the mean 5.
1 Jülich Supercomputing Center,
We performof the dynamics
of clusters
N stars using the simulation code Nbody6++. In the
ForschungszentrumJülich, 52428
Germany
value of the final truncation
disc size to measurethesimulations
cluster environment’se↵ect Jülich,on containing
formingplanetarysystems.
2 Max-Planck-Institut fürsimulations,
a
star
is
represented
by
a
test
particle
with
a
given
mass, position and velocity. The particles’ positions
Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
are chosen so that the resulting stellar number density distribution obeys a King profiles (King 1966). The reason isDiscthat King 3.2.
profilesseemtruncation
to present young clusters well (Hillenbrand & Hartmann 1998; Harfst et al. 2010; Nürnberger(Received
August24, 2021)&
Petr-Gotzens
2002);
however,there
no generalThis
observational
whether
W=7, W=8 or W=9 is the best
We start with the artificial situation of all disc having a size of 1000
auisinitially.
value isagreement
typical forthe top-end
choice. Here we mostly use a King profile with W0 > 9, but also discuss in section 5 the role of the chosen profile forof the size of observed discs. (Andrews et al. 2018b).
The reason
Submitted
to ApJis that we only intend to obtain a measure of the
our results.
influence of the environment rather than aTherealistic
representationof thedisc masses
sizes infollowing
a specificThus
we obtain
individual
test particles are
assignedthe cluster.
initial massfunction(IMF) by Kroupa (2001)
the size that a standard disc is truncated to due to thecloseflybysofaparticularclusterenvironment.
Thistruncation
8
ABSTRACT
0.3>
0.01  m/M < 0.08
disc size is obtained the following way:<m
1.3
Numerous
protoplanetary
discs
show
distinctspiralarms
features.
While
by a
⇠(m)
/
(2)
m
0.08possibly
them/M<caused0.5of a parameterIn this work, all encounters were assumed to be parabolic, prograde and coplanar,using
results
>
:
2.3
range
of processes,
detailed
patternbyanalysis
pointsat closestellarflybys
as 0.5cause
forof them.
m truncation
 m/M
study
of star-discencounters
performed
Breslau etal. (2014).
In thiscase,
thedisc, somesize rdisc
after theØ シミュレーションは全星ができた状態で始め、2 Myr後にガスが無くなってクラスターが膨張。Nbody6++。
Ø King W=9に従って星を配置。星質量はあるIMFを使⽤。0.08-150 Msun。primordial singleのみ。
Ø ビリアル平衡を仮定してマクスウェル速度分布。
Ø モデルP: γ=1.7, K: γ=2.0, V: Rc=1.3 pc。
Surprisingly,
these discs
young low-mass clusters, where close stellar flybys are expected to
flyby
can be approximatedas reside in withthe lower mass limit set to 0.08M (hydrogen burning limit) and an upper mass limit of 150M , which is theØ flyby発⽣時の質量⽐と近星点距離を記録。シミュレーション後にトランケーションを評価。
be rare. This fact motivated us acceptedtoa freshlook(Figerat theof closeflybys in low-mass( takeuppermass limit
2005;frequency
Weidner & Kroupa2006).
0.32
Ø 初めは円盤半径は全て1000
au。
0.28· rall· clusters
M21
ifsingle
rdiscsmaller< rprevious
Initially,clusterstars ,arestars,
meaning
primordial
binaries
excluded.
This approximation is chosenclusters. In the solar neighbourhood,
low-masshave
half-mass
radii
thanastheir
more
peri
rdisc =for simplicity(3) inclusion of primordialintheconsecutive
processing
of
the
e↵ect
of
flybys
on
disc
size,
where
the
M1:
host
,
iffactrdiscresults
rprevious
massive counterparts. We show thatrprevious
this observationalin the M21=M2/M1,
mean and central stellar
binaries would significantly increase the complexity. The absence of primordial binaries might lead to underestimating
density of low-mass clusters beingejections
approximately
theduesame
as inhigh-massclusters,
is 1975).
rarelyfrom the cluster
to strong
few-bodyinteractions(Heggie which
1975; Hills
For a more detailed discussion
Ø V以外は⼤⼩質量で同じ。
where rperi is the periastron distance between
the
two
stars
involved
in
a
close
flyby
in
au.
M
=
M
/M
is
the mass
21
2
1
of omitting
binariesin this typeof study, see of
Olczakal. (2006).dynamics
Besides, binaries
can accelerate a cluster’s expansion
reflected in theoretical studies. We
perform
N-bodysimulations
theet stellar
in young
ratio ofPは13%が1000
the two stars, with M1 auのまま。
being theprocess
massandof thedisc-hostingstar
M(Kroupa
theKaczmarek
perturberetstar.
ris
2 that of
reducetheTakingencounterrateandfaster2001;
Ø ⼩質量Vは95%,
clusters obeying the observed
mass-radiusrelation.
themeandisc truncation
radiusal.as2011).
a previous
proxy
the disc size before the encounter.
Potentially existing initial mass segregation in the clusters was neglected. The occurrence of mass segregation withinfor thedegreeof influence
we debate,
find
that
the
thetheenvironment
discs
Ø Pは10-30%が30
au以下に。
clusters
underwhether
it influenceisthe
primordial—foundindiscthe dense
core onphaseof theInclinedor retrograde
flybys ofarethelessenvironment,
efficientis still
to removematerial
fromdisc, ofandtruncation
radii
arecluster formation process
is
more
or
less
the
same
in
lowandhigh-massclusters.
Eventhefraction
of
small
discs
(<
10
au)
is
therefore larger (Clarke & Pringle 1993; Bhandare et al. 2016). However, the di↵erence to coplanar flybys is notØ 配置をPlummer（〜⼩W）にしても同じ。
nearly identical.
Our mainconclusion
is thatthe frequency
flybys radii
seemshaveto haveseverely
extremely
large. Applying
the result
for coplanar
flybysmeans that oftheclosetruncation
to be been
regardedas lowerØ スパイラルが⾒つかったTaurusで<30
auはunderestimatedfor low-mass
A testablepredictionof thisis thatlow-massclusters
limits.
Similar, hyperbolic
flybys clusters.
are less e↵ectivein disctruncation
thanhypothesis
parabolic ones.
However,more
eccentric
クラスター中⼼にいて、割合も合う。
should contain
10%–15%more
of discs
smallerthan 30mass
au truncatedby flybys.
Thesefocusestruncateddiscs should
encounters
are significantly
commonin higher
clusters. Since
our study
on low-massclusters,
approximating
the flybyfrom
pathprimordiallyas parabolic issmall
likely discs
justified.
Moreover,applyedge.
these simplifications in all simulations
be distinguishable
by theirsteepweouterØ 端が鋭利なら光蒸発とも区別可能。
identically, so the relative di↵erences between the diverse cluster environments should not be a↵ected. Here, only aØ 今まで知られていないのは⼩質量の10%だ
summaryof themethod has been
provided;
weinterstellar
refer the reader
to Breslau et al. (2014) for more information on the
Keywords:
protoplanetarydisks— localmatterdisc
size
determination
duringclose
stellarencounters.
と統計的に不⾜、0.1 Myrで中⼼から離れる
100 au以下になるエンカウント
ため。
3.3. 1.Observational
uncertainties
INTRODUCTION
Close stellar flybys common in low-mass clusters

13

Figure 7. Location of the members of the Taurus star forming region according to Luhman (2018). The systems showing
spiral arms attributed to stellar flybys (Cuello et al. 2019; Ménard et al. 2020) are indicated as blue circles. Note, that someare located close to each other, so that they are better identified by their labels. The red symbols indicate the location of thestars for which Andrews et al. (2018a) determined disc sizes. Discs with sizes smaller than 30 au are indicated as red crosses,
whereas discs with sizes exceeding 30 au are shown as circles. The disc sizes for RW Aur A & B were taken from Rodriguez
et al. (2018b)

As mentioned
in section
the exponent
in the derived discs
mass-radiusrelations Andrews
for young2020;
stellarvan
clustersin theet al.
The observed
diversityin the1,structureof protoplanetary(for example,
der Marelsolar neighbourhood varies between.
= 1.7 — 2.0 and that of theconstant spans the range of C= 359 — 1263.

LETTER TO THE EDITOR

An unbiased NOEMA 2.6 to 4 mm survey of the GG Tau ring: First
detection of CCS in a protoplanetary disk
N. T. Phuong1,2 , A. Dutrey3 , E. Chapillon3,4 , S. Guilloteau3 , J. Bary5 , T. L. Beck6 , A. Coutens7 ,O. Denis-Alpizar8 , E. Di Folco3 , P. N. Diep2 , L. Majumdar9 , J.-P. Melisse3,4 , C.-W. Lee1,10 , V. Pietu4 ,T. Stoecklin11 , and Y.-W. Tang121

Korea Astronomy and Space Science Institute, 776 Daedeokdae-ro, Yuseong-gu, Daejeon, Korea
e-mail: tpnguyen@kasi.re.kr
2
Department of Astrophysics, Vietnam National Space Center, Vietnam Academy of Science and Techonology, 18 Hoang Quoc
Viet, Cau Giay, Hanoi, Vietnam
N. T. Phuong et al.: An unbiased NOEMA 2.6 to 4 mm survey of the GG Tau ring
3
Laboratoire d’Astrophysique de Bordeaux, Université de Bordeaux, CNRS, B18N, Allée Geo↵roy, Saint-Hilaire 33615,
Pessac, France
4
IRAM, 300 Rue de la piscine, 38406 Saint Martin d’Hères Cedex, France
5
Department of Physics and Astronomy, Colgate University, 13 Oak Drive, Hamilton, New York 13346, USA
6Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, Maryland 21218, USA
7
Institut de Recherche en Astrophysique et Planétologie, Université de Toulouse, UPS-OMP, CNRS, CNES,
9 Av. du Colonel Roche, 31028 Toulouse Cedex 4, France
8
Instituto de Ciencias Químicas Aplicadas, Facultad de Ingeniería, Universidad Autónoma de Chile, Av. Pedro de Valdivia 425,
7500912 Providencia, Santiago, Chile
9
School of Earth and Planetary Sciences, National Institute of Science Education and Research, HBNI, Jatni 752050, Odisha, India
10
University of Science and Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon 34113, Republic of Korea
11
Institut des Sciences Moléculaires, UMR5255-CNRS, 351 Cours de la libration, 33405 Talence, France
12
Academia Sinica Institute of Astronomy and Astrophysics, PO Box 23-141, Taipei 106, Taiwan

Ø PPDの化学は分⼦雲からとその場。GG Tau Aa-Ab1/b2の円盤は⼤きく、低温で重いので分⼦化学向き。
Ø NOEMA, 2.6-4.2 mm, dv=0.2-0.3 km/s

Received 27 July 2021 / Accepted 27 August 2021
ABSTRACT
Context. Molecular line surveys are among the main tools to probe the structure and physical conditions in protoplanetary disks(PPDs), the birthplace of planets. The large radial and vertical temperature as well as density gradients in these PPDs lead to acomplex chemical composition, making chemistry an important step to understand the variety of planetary systems.
Aims. We aimed to study the chemical content of the protoplanetary disk surrounding GG Tau A, a well-known triple T Tauri system.
Methods. We used NOEMA with the new correlator PolyFix to observe rotational lines at ⇠2.6 to 4 mm from a few dozen molecules.
We analysed the data with a radiative transfer code to derive molecular densities and the abundance relative to 13 CO, which wecompare to those of the TMC1 cloud and LkCa 15 disk.
Results. We report the first detection of CCS in PPDs. We also marginally detect OCS and find 16 other molecules in the GG Tauri
outer disk. Ten of them had been found previously, while seven others (13 CN, N2 H+ , HNC, DNC, HC3 N, CCS, and C34 S) are newdetections in this disk.
Conclusions. The analysis confirms that sulphur chemistry is not yet properly understood. The D/H ratio, derived from DCO+ /HCO+ ,DCN/HCN, and DNC/HNC ratios, points towards a low temperature chemistry. The detection of the rare species CCS confirms that
GG Tau is a good laboratory to study the protoplanetary disk chemistry, thanks to its large disk size and mass.
Key words. astrochemistry – molecular data – protoplanetary disks – stars: individual: GG Tau A

1. Introduction

layer just below it, while molecules are frozen on the dust

LETTER TO THE EDITOR

An unbiased NOEMA 2.6 to 4 mm survey of the GG Tau ring: First
detection of CCS in a protoplanetary disk
N. T. Phuong1,2 , A. Dutrey3 , E. Chapillon3,4 , S. Guilloteau3 , J. Bary5 , T. L. Beck6 , A. Coutens7 ,O. Denis-Alpizar8 , E. Di Folco3 , P. N. Diep2 , L. Majumdar9 , J.-P. Melisse3,4 , C.-W. Lee1,10 , V. Pietu4 ,T. Stoecklin11 , and Y.-W. Tang121

Korea Astronomy and Space Science Institute, 776 Daedeokdae-ro, Yuseong-gu, Daejeon, Korea
e-mail: tpnguyen@kasi.re.kr
2
Department of Astrophysics, Vietnam National Space Center, Vietnam Academy of Science and Techonology, 18 Hoang Quoc
Viet, Cau Giay, Hanoi, Vietnam
3
Laboratoire d’Astrophysique de Bordeaux, Université de Bordeaux, CNRS, B18N, Allée Geo↵roy, Saint-Hilaire 33615,
Pessac, France
4
IRAM, 300 Rue de la piscine, 38406 Saint Martin d’Hères Cedex, France
5
Department of Physics and Astronomy, Colgate University, 13 Oak Drive, Hamilton, New York 13346, USA
6Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, Maryland 21218, USA
7
Institut de Recherche en Astrophysique et Planétologie, Université de Toulouse, UPS-OMP, CNRS, CNES,
9 Av. du Colonel Roche, 31028 Toulouse Cedex 4, France
8
Instituto de Ciencias Químicas Aplicadas, Facultad de Ingeniería, Universidad Autónoma de Chile, Av. Pedro de Valdivia 425,
7500912 Providencia, Santiago, Chile
9
School of Earth and Planetary Sciences, National Institute of Science Education and Research, HBNI, Jatni 752050, Odisha, India
10
University of Science and Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon 34113, Republic of Korea
11
Institut des Sciences Moléculaires, UMR5255-CNRS, 351 Cours de la libration, 33405 Talence, France
12
Academia Sinica Institute of Astronomy and Astrophysics, PO Box 23-141, Taipei 106, Taiwan
Received 27 July 2021 / Accepted 27 August 2021
ABSTRACT
Context. Molecular line surveys are among the main tools to probe the structure and physical conditions in protoplanetary disks(PPDs), the birthplace of planets. The large radial and vertical temperature as well as density gradients in these PPDs lead to acomplex chemical composition, making chemistry an important step to understand the variety of planetary systems.
Aims. We aimed to study the chemical content of the protoplanetary disk surrounding GG Tau A, a well-known triple T Tauri system.
Methods. We used NOEMA with the new correlator PolyFix to observe rotational lines at ⇠2.6 to 4 mm from a few dozen molecules.
We analysed the data with a radiative transfer code to derive molecular densities and the abundance relative to 13 CO, which wecompare to those of the TMC1 cloud and LkCa 15 disk.
Results. We report the first detection of CCS in PPDs. We also marginally detect OCS and find 16 other molecules in the GG Tauri
outer disk. Ten of them had been found previously, while seven others (13 CN, N2 H+ , HNC, DNC, HC3 N, CCS, and C34 S) are newdetections in this disk.
Conclusions. The analysis confirms that sulphur chemistry is not yet properly understood. The D/H ratio, derived from DCO+ /HCO+ ,DCN/HCN, and DNC/HNC ratios, points towards a low temperature chemistry. The detection of the rare species CCS confirms that
GG Tau is a good laboratory to study the protoplanetary disk chemistry, thanks to its large disk size and mass.
Key words. astrochemistry – molecular data – protoplanetary disks – stars: individual: GG Tau A

1. Introduction

layer just below it, while molecules are frozen on the dustの続き

arXiv:2108.12401v1 [astro-ph.SR] 27 Aug 2021

arXiv:2108.12401v1 [astro-ph.SR

explicar las curvas de luz usando extinción de la radiación estelar y emisión del polvo
dentro del agujero de un disco transicional. Para el objeto Mon-296, fuimos exitosos.
Sin embargo,para Mon-314 y Mon-433
nuestraevidencia
sugiereque contienen un
INTERPRETATION
OF OPTICALAND
IR LIGHT
CURVES
disco pre-transicional.
Para Mon-1308
un nuevo
ajuste delIN
espectro
usando el código
FOR TRANSITIONALDISKSCANDIDATESNGC 2264
USING THE
EXTINCTEDSTELLAR
RADIATIONAND
3D de transferencia
radiativa
Hyperion,nos permite
concluir que
esteTHEobjeto tiene un
EMISSION
OFdeOPTICALLY
THIN DUSTEsto
INSIDE
THEconHOLE
disco completoen lugar
un disco transicional.
coincide
trabajo previo
E. Nagel,1 F. Gutiérrez-Canales,1 S. Morales-Gutiérrez,1 and A.P. Sousa2
sobre Mon-1308 y también conque no podemos encontrar un ajuste de las curvas deDraft version: August 30, 2021
luz solo usando la contribución del polvo dentro del agujero de un disco transicional.
RESUMEN
En la región de formación estelar ABSTRACT
NGC 2264 hay objetos catalogados como anfitriones de discos transicionales de acuerdo al modelaje de su espectro. Cuatro miemIn the stellar
forming
regionthere
arey enobjects
catalogued
as hosting abros deeste conjunto
tienenNGC
curvas 2264
de luz enel óptico
el infrarrojo
provenientes
los telescopios
CoRoT
y Spitzer. modeling.
En este trabajo,
tratamosde simultáneamentetransitionalde disk
accordingto spectra
Fourmembersof this set have opexplicar
lascurvasdeluz
usandoextincióndelaradiaciónestelaryemisióndel polvo
tical and infrared light curves coming from the CoRoT and Spitzer
telescopes. In
dentro del agujero de un disco transicional. Para el objeto Mon-296, fuimos exitosos.
this work,Sinweembargo,
try to para
simultaneously
explaintheevidencia
light curves
using
the extinctionofMon-314 y Mon-433
nuestrasugiere que
contienenunpre-transicional.
Para Mon-1308
un dustnuevo ajuste
deltheespectrocódigothe stellar disco
radiationand emissionof theinsideholeusandoof aeltransitional
disk.
3D
de
transferencia
radiativa
Hyperion,nospermiteconcluirqueeste
objetotieneun
For the object Mon-296, we were successful to do this. However, for Mon-314, anddisco completo en lugar de un disco transicional. Esto coincide con trabajo previo
Mon-433 ourthathost aencontrar
pre-transitional
disk.
For
sobreevidence
Mon-1308 ysuggests
también conquetheyno podemos
un ajuste de lascurvasde Mon-1308
a new spectra
fitting
the 3Ddelradiative
transfer
code
Hyperion,allow us to conluz solo
usandousing
la contribuciónpolvo dentrodel agujero
de undisco transicional.

clude that this object host a full-diskABSTRACT
instead of a transitional disk. This is in accord
to previousInwork
on Mon-1308
wecatalogued
cannot asfind
a fitaof the lightthe stellar
forming regionand
NGCwith2264thetherefact
are that
objects
hostingtransitional
disk
accordingtospectramodeling.
Fourmembers
of
this
set
have
opcurves only using the contribution of the dust inside the hole of a transitional
disk.

tical and infrared light curves coming from the CoRoT and Spitzer telescopes. In
this work, we try to simultaneously explain the light curves using the extinction of
Key Words:
dust, extinction — protoplanetary disks — Stars: pre-main sequence
the stellar radiation and emission of the dust inside the hole of a transitional disk.
For the object Mon-296, we were successful to do this. However, for Mon-314, andMon-433 our evidence suggeststhat they host a pre-transitional disk. For Mon-1308
1. INTRODUCTIONa new spectra fitting using the 3D radiative transfer code Hyperion, allow us to conclude that this object host a full-disk instead of a transitional disk. This is in accord
Spectraland photometric variability of young stellar objects (YSOs) is the usual
to previous work on Mon-1308 and with the fact that we cannot find a fit of the lightoutcome ofcurves
multiwavelength
campaigns
2015, 2016;
only using the contributionof the(Stauffer
dust inside etthe al.
hole2014,
of a transitional
disk. Cody et al.

分光モデルから遷移円盤と思われる4個の⽳の光度曲線とダスト放射を説明。
Mon-296はOK。Mon-314とMon-433はpre遷移円盤ではない。Mon-1308は完全な円盤。
2014; Morales-Calderón
et al. 2011). For young stars in NGC 2264, Stauffer et al.
Key Words: dust, extinction — protoplanetary disks — Stars: pre-main sequence(2014) extract the accretion burst dominated light curves (lcs), Stauffer et al. (2015)
1. INTRODUCTION
1 Departamento de Astronomı́a, Universidad
de Guanajuato, México.
2 Institut de Planétologie et d’Astrophysique de Grenoble, Université Grenoble Alpes, France.

Spectral and photometric variability of young stellar objects (YSOs) is the usual

Universidad Nacional Autónoma de México, Instituto de Ciencias Nucleares, AP 70-543, CDMX 04510, México

An orbital release model for the Orion BN/KL fingers

2 Observatoire de Grenoble, Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France

3 Universidad Nacional Autónoma de México, Instituto de Astronomı́a, AP 70-264, CDMX 04510, México

N ED
ED
IT M
ED A
MN
A U
N S
U C
SC R
RI I
PT PT

Universidad Nacional Autónoma de México, Instituto de Ciencias Nucleares, AP 70-543, CDMX 04510, México

2 Observatoire de Grenoble, Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France

3 Universidad Nacional Autónoma de México, Instituto de Astronomı́a, AP 70-264, CDMX 04510, México

ABSTRACT

We present a simple model in which the bullets that produce the “Orion fingers”(ejected by the BN/KL object)
are interpreted as protoplanets or low mass protostars
ABSTRACTWe
present
a simplein which theexplosion.
bullets that produce
the remnant“Orion fingers”
in orbit around a high mass star that
hasmodela supernova
As theof(ejected by the BN/KL object) are interpreted as protoplanets or low mass protostars
the SN explosion has only ain small
fraction of the mass of the pre-supernova star, theorbit around a high mass star that has a supernova explosion. As the remnant oforbiting objects then movetheaway
in freehastrajectories,
preserving
orbital velocity
SN explosion
only a small fractionof the mass theirof the pre-supernova
star, theorbiting
objects
then
move
away
in
free
trajectories,
preserving
theirorbitalvelocityat the time of release. We atshowthat a system of objects arranged in approximately
the time of release. We show that a system of objects arranged in approximately
co-planar orbits results in trajectories
with inmorphologicaland kinematicalcharacterco-planar orbits results
trajectories with morphologicaland kinematicalcharacteristics
resembling
the
Orion
fingers.
We
show
that,
under
the
assumptionof
constant
istics resembling the Orionvelocity
fingers.
We show that, under the assumption of constantmotions, the positions of the observed heads of the fingers can be used to revelocity motions, the positionsof theobservedof the fromfingers
used resulting
to reconstruct
the propertiesof theheads
orbital structure
whichcanthey beoriginated,inacompact
disk
with
an
outer
radiusof⇠
2.4
AU.
construct the properties of the orbital structure from which they originated, resulting
Key words:
stars:ofwinds,
outflows– ISM: jets and outflows – ISM: individual objects:
in a compact disk with an outerradius⇠ 2.4
AU.
Orion BN/KL

Key words: stars: winds, outflows – ISM: jets and outflows – ISM: individual objects:
Orion BN/KL
1

INTRODUCTION“bullets” at the tips of the fingers might have had a substantial slowing down along their trajectories (Rivera-Ortiz
et al. 2019a, b).
It has been proposed that this sudden event could havebeen:

Downloaded from https://academic.oup.com/mnrasl/advance-article/doi/10.1093/mnrasl/slab072/6328490 by KASI user on 21 Oc

A. C. Raga1? , P. R. Rivera-Ortiz2 , J. Cantó3 , A. Rodrı́guez-González1
A.
Castellanos-Ramı́rez3
1

ownloaded from https://academic.oup.com/mnrasl/advance-article/doi/10.1093/mnrasl/slab072/6328490 b

3
A.
Castellanos-Ramı́rez
1Orion fingersを作る弾丸が超新星爆発した⼤質量星を周る⼩質量原始星だというモデル。
One of the more intriguing flows in the Galactic ISM is com-U

IT

posed of the “fingers” radiating away from the Orion BN/KL
object (a massive embedded stellar object surrounded by acompact nebula, see Becklin & Neugebauer 1967 and Klein1 INTRODUCTION“bullets” at the tips• ofthe fingers might have had a submann & Low 1967). These fingers were first reported bya result of the breakup of a multiple stellar system
slowingalong theirtrajectories(Rivera-Ortiz
Allen & Burton (1993) in IR H2 andstantial[Fe II] emission,and down(see Reipurth
et al. 2010;
Reipurth & Mikkola
2012, 2015;
One of the more intriguing flows in the Galactic
have lengthsISM
of up isto com⇠ 20 (⇠ 7 ⇥ 1017 cm at a distance ofMoeckel&Goddi2012).
The
dynamicsofa
3or 4-body
et al. 2019a, b).
400thepc). Orion
The basesof these fingers, as well as a multitude
system can explain the high velocities of the stars near the
posed of the “fingers” radiating away from
BN/KL
It hasbeen proposedthat(observed
this sudden
event could
haveof other filaments are detected in interferometric
CO obserBN/KL object
in radio continuumby Rodrı́guez
object (a massive embedded stellar object
surrounded
by aZapata et al. 2017; Bally et al.
vations(Zapata et al. 2009;
et
al.
2005
and
by
Gómezetal.
2005,2008),
andin
principle
been:2017).1967
The COemissionof the individual fingers shows a
could trigger the “interstellar bullets that create the Orioncompact nebula, see Becklin & Neugebauer
andKlein-